擴展黃永剛原始晶體塑性程序加入AF背應力模擬金屬疲勞問題參考文獻:《Low-cycle fatigue life prediction of a polycrystalline nickel-base superalloy using crystal plasticity modelling approach》在原始程序中修改流動方程,加入背應力項,引入運動硬化項,從而可以描述多晶金屬循環加載中的包辛格效應
對于本例中所采用的現象學晶體塑性本構關系,共有114個狀態變量,分別如下:SDV1:HomogenizationCountSDV2:GrainCountSDV3:CrystalliteCountSDV4~7:quaternion(取向四元數)SDV8~11:deviation from initial orientation as axis (1-3) and angle in