资源描述：
《umat程序代码及注释.docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、SUBROUTINEUMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,1RPL,DDSDDT,DRPLDE,DRPLDT,2STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,3NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,4CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)CINCLUDE'ABA_PARAM.INC'CCHARACTER*8
2、0CMNAMEDIMENSIONSTRESS(NTENS),STATEV(NSTATV),1DDSDDE(NTENS,NTENS),2DDSDDT(NTENS),DRPLDE(NTENS),3STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),4PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3)CDIMENSIONSTRANT(6),CFULL(6,6),CDFULL(6,6),DCDDFV(6,6),1DCDDMV(6,
3、6),DCDDDV(6,6),DFFDE(6),DFMDE(6),DFDDE(6),2TEMP1(6),TEMP2(6),TEMP3(6),TEMP4(6),TEMP5(6),TEMP6(6),3TEMP7(6,6)PARAMETER(ZERO=0.D0,ONE=1.D0,TWO=2.D0,THREE=3.D0)C输入材料弹性参数e1=props(1)e2=props(2)e3=props(3)u12=props(4)u13=props(5)u23=props(6)g12=props(7)g13=props(8)g23=props(9)C输
4、入材料的强度参数，用于失效判断。xt=props(10)xc=props(11)yt=props(12)yc=props(13)s12=props(14)s13=props(15)s23=props(16)zt=props(17)zc=props(18)eta=props(19)粘性系数gc1=props(20)复合材料的1方向的临界能量释放率gc2=props(21)gc3=props(22)!CFULL无损伤下的刚度阵u21=(u12/e1)*e2u31=(u13/e1)*e3u32=(u23/e2)*e3fact=1-u12*u21-
5、u23*u32-u13*u31-&2*u21*u32*u13factor=fact/(e1*e2*e3)CFULL(1,1)=(1-u23*u32)/(e2*e3*factor)CFULL(2,2)=(1-u13*u31)/(e1*e3*factor)CFULL(3,3)=(1-u12*u21)/(e1*e2*factor)CFULL(1,2)=(u12+u32*u13)/(e1*e3*factor)CFULL(2,3)=(u32+u12*u31)/(e1*e3*factor)CFULL(1,3)=(u31+u21*u32)/(e3*e2*
6、factor)CFULL(3,1)=CFULL(1,3)CFULL(3,2)=CFULL(2,3)CFULL(2,1)=CFULL(1,2)CFULL(4,4)=g12CFULL(5,5)=g13CFULL(6,6)=g23CCRECOVERELASTICSTRAINS更新应变CDO10K1=1,NTENSSTRANT(K1)=STATEV(K1)+DSTRAN(K1)10CONTINUECCSTORESTRAINSINSTATEVARIABLEARRAY存储应变值，用于下一步的应变更新CDO20K1=1,NTENSSTATEV(K1)=S
7、TRANT(K1)20CONTINUE!damagestrain计算各种失效情况下发生失效时的应变极限值epsilonxt=xt/CFULL(1,1)epsilonxc=xc/CFULL(1,1)epsilonyt=yt/CFULL(2,2)epsilonyc=yc/CFULL(2,2)epsilonzt=zt/CFULL(3,3)epsilonzc=zc/CFULL(3,3)epsilons12=s12/CFULL(4,4)epsilons13=s13/CFULL(5,5)epsilons23=s23/CFULL(6,6)!fiberte
8、nsedamage失效准则if(STRANT(1).GE.0.0)theneft=(STRANT(1)/epsilonxt)**2+(STRANT(4)/epsilons12)*
