v1.1版本功能更新如下:
- 支持荷载类型:位移加载、表面压强加载(集中力我舍弃了,建议使用位移加载)
- 单元类型:C3D4、C3D10、C3D8、C3D20,后续的更多三维单元类型将逐步添加
- 内核求解:C3D20 单元增加了 14 点积分求解技术,保持精度的前提下,减少原有的高斯积分数量
- 稀疏矩阵存储格式:使用三元组(Triplet)形式,即行索引数组、列索引数组和非零值数组,分别对应矩阵中非零元素的位置和它们的值,大大提升了整体刚度矩阵的组装效率
- 线性方程组求解技术:当系统自由度个数小于 1W 时,采用直接矩阵求解技术;当系统自由度个数大于 1W 时,采用 PCG 迭代求解技术,预条件使用Incomplete Cholesky,可以大大提升刚度方程的求解效率
- 边界条件:使用罚函数技术计算出节点位移后,一键解出节点的支反力
- 前处理:支持一键提取 Abaqus-inp 文件的节点、单元、节点集、单元集、surface 等信息,供 MFEA 前处理使用
- 配置文件:使用 json 文件作为程序的配置文件,参数化控制材料信息、边界条件、后处理等
- 后处理:在 json 文件中添加后处理的控制选项,可在 Matlab 环境下进行云图绘制,可选择:场变量、单元边界颜色、面颜色、colormap、discretize 离散程度等参数
- 计算结果文件(.dat),记录各节点的坐标、位移、支反力的值
- vtk 文件,可导入 Paraview 中进行可视化显示,或者使用 Python 的 Pyvista 库进行可视化也都是可以的
- log 文件,记录计算各个过程的具体时间明细以及一些模型信息。
具体各个模块的功能,将在明天做详细介绍,今天先测试了几个案例,整套程序已上传至知识星球中,后台回复:星球
,即可加入。