拉格朗日求解器的计算过程大致如下:
- 在初始计算中,通过对节点速度积分得到新节点位置;
- 通过新节点的坐标位置,计算新单元密度和单元应变率;
- 通过应变,新的应力可以计算得到;
- 对于各向同性材料,单元变形能分成两个独立的部分;
- 体积改变是由球应力张量引起;
- 形状改变是由于偏应力张量引起;
- 球应力由状态方程(EOS) 求解;
- 偏应力是由胡克定律(弹性)和塑性屈服准则(本构关系)求解 - 即 AUTODYN 的强度模型;
- 再加上失效模型就构成了一个完整的材料模型。
- 一旦计算单元的内部应变率被决定,每一个节点的受力能被计算;在这个阶段,任何从边界条件或者与其它物体作用的外力作为最后的节点力;
- 知道节点力,可以计算节点加速度;积分节点加速度可以得到节点新的速度;
- 完整的一个计算周期,重复这样的周期达到要求的求解时间或循环次数。
计算流程图如下: