!!!Duncan-ChangModel
!!!Symbols:c-粘滞力,Fai-内摩擦角,
Sf-破坏强度(pl-p3)f,Ds-应力水平,
Pa-大气压,P3-围压
*CREATE,Duncan-Chang
!建立邓肯张模型的宏文件
*afun,deg
!将单位转化为角度
Pa=1e5
!给定大气压
*set,P1,-ArrS3(i)!
*set,Pl,-ArrS2(i)
*set,P3,-ArrSl(i)
*if,P3,LT,0.1*Pa,then
P3=0.1*Pa
*endif
Sf=2*(c*cos(Fai)+((P2+P3)/2)*sin(Fai))/(1-sin(Fai))
!Mohr-Coulomb破坏强度
Ds=(P1-((P2+P3)/2))/Sf
*if,Ds,GT,0.95,then
Ds=0.95
*endif
*if,MaxPs(i),LT,P1-(P2+P3)/2,then
Ei=k*Pa*(((P2+P3)/2)/Pa)**n
Et=Ei*(1-Rf*Ds)**2
MaxPs(i)=P1-(P2+P3)/2
*elseif,MaxPs(i),GE,P1-(P2+P3)/2
Et=Kur*Pa*(((P2+P3)/2)/Pa)**n
*endif
mp,ex,i,Et
!修改单元i的Et
mp,nuxy,i,Mu
*END