塑性有限元几种方法的收敛过程.docx

弹塑性有限元法绘图并说明：1）全量割线刚度法的收敛过程；2）加快收敛的全量割线刚度法的收敛过程；3）全量初应力法的收敛过程；4）（增量）切线刚度法的收敛过程；6）增量初应力法的收敛过程。（1）假设材料初始形态为全弹性，施加全部荷载，由割线刚度E’(1)可求得应力A(1)、应变；通过应变求得曲线对应的相当应力值B(1)，再由曲线求得下一次迭代的割线刚度E’(2)；施加全部荷载，根据割线刚度E’(2)求得应力A(2)、应变；由应变求得曲线对应的相当应力值B(2)，再由曲线求得下一次迭代的割线刚度E’(3)；依次循环直到A(n)与B(n)满足收敛的条件。（2）假设材料初始形态为全弹性，施加全部荷载，由割线刚度E’(1)求得应力A(1)、应变，由应变求得曲线对应的相当应力值B(1)，根据本次迭代求得的和的大小来确定的正负，计算出下次迭代使用的。第一次迭代后得到,则根据求得的进行第二次迭代，B(2)由曲线上对应的点得到，通过比较、的大小来确定的正负，若,则若,则按之前的方程进行计算；直到、满足收敛的条件循环停止。（3）一次性加上全部的外载荷，形成节点载荷列矩阵，按全弹性进行迭代，计算位移和应变，获得应力，根据应变状态判断单元是否进入塑性状态；由求得曲线对应的相当应力值，处于弹性状态单元；对塑性单元计算初始应力及节点初载荷列矩阵，并组装成总的节点塑性初载荷列矩阵；将所得节点塑性初载荷列矩阵带入全结构节点平衡方程进行第二次迭代，过程与第一次相同，不同之处在于全结构节点平衡方程右侧增加了节点塑性初载荷列矩阵；随着迭代次数增加，塑性元提供的初应力也不断加大，位移及应变也不断加大，直到两次迭代应变差达到一定精度。（4）切线刚度法将外载荷分为若干增量段，然后将载荷逐段施加。由刚度E’(1)（完全弹性）求得在载荷作用下的位移；第二加载段开始，刚度为曲线所对应点的切线斜率E’(2)，作用下产生位移，总位移为，根据曲线求出所对应点的切线斜率供下次循环使用；对于加载过程中出现的过渡元，通过加权系数公式给出加权平均刚度矩阵。（5）（5）增量的初应力法增量的初应力法将外载荷分为若干增量段，逐段施加载荷增量。增量的初应力法在每一段内按照全量初应力法的步骤进行计算，在各增量段相邻两次迭代应变增量差达到一定精度时，进入下一个增量段的迭代过程；第一次刚度矩阵为完全弹性刚度矩阵，第二增量段刚度矩阵由曲线所对应点的切线斜率决定，以此类推。