4.4　本构模型计算的收敛算法

模型的开发、验证和应用离不开算法的保证。为了把模型编入有限元计算程序,必须在输入的每一个增量环节的计算中不断更新应力、应变和硬化变量等[206](Crisfield,1991)。除非增量微乎其微地小,否则必然会在计算过程中造成一定的无法容忍的误差。而实际上为了追求计算效率,输入的增量又不可能过于小,这样,只通过向前试算式的欧拉算法格式便无法保证误差的消除。数值积分时能否保证应力、应变和硬化变量等的不断更新,对于计算精度和准确性影响巨大。在模型的数值分析应用中,饱和土模型的数值积分无论是显式算法还是隐式算法都已比较成熟。