擴展黃永剛原始晶體塑性程序加入AF背應力模擬金屬疲勞問題

參考文獻：《Low-cycle fatigue life prediction of a polycrystalline nickel-base superalloy using crystal plasticity modelling approach》

在原始程序中修改流動方程，加入背應力項，引入運動硬化項，從而可以描述多晶金屬循環加載中的包辛格效應

背應力的演化遵循AF模型

并使用原始的PAN模型描述滑移系統的硬化行為

為了表征多晶的疲勞壽命，引入兩類疲勞指示因子分別為

一：累計塑性滑移

二：累計能量耗散