在原始程序中修改流动方程，加入背应力项，引入运动硬化项，从而可以描述多晶金属循环加载中的包辛格效应，背应力的演化遵循AF模型，并使用原始的PAN模型描述滑移系统的硬化行为，为了表征多晶的疲劳寿命，引入两类疲劳指示因子分别为：累计塑性滑移，累计能量耗散