Abaqus子程序UMAT和VUMAT的例子Abaqus子程序UMAT和VUMAT的例子Abaqus子程序UMAT和VUMAT的例子Abaqus子程序UMAT和VUMAT的例子編制僅供參考審核批準生效日期地址：電話：傳真:郵編:From*******************************************************************************UMATFORABAQUS/STANDARDINCORPORATINGELASTICBEHAVIOURFORPLANE****STRAINANDAXI-SYMMETRICELEMENTS.*************************************************************************************************************************************************************USERSUBROUTINESUBROUTINEUMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,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''CCHARACTER*80CMNAMECCDIMENSIONSTRESS(NTENS),STATEV(NSTATV),1DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),2STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),3PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3)CCPARAMETER(M=3,N=3,ID=3,ZERO=,ONE=,TWO=,THREE=,+SIX=,NINE=,TOLER=CDIMENSIONDSTRESS(3)CC--------------------------------------------------------------------CSPECIFYMATERIALPROPERTIESCE=PROPS(1)XNUE=PROPS(2)CCCSETUPELASTICITYMATRIXCE11=E/*XNUE)E12=E*XNUE/*XNUE)EG2=E/(ONE+XNUE)EG=EG2/TWOCCDETERMINESTRESSINCREMENTCDSTRESS(1)=E11*DSTRAN(1)+E12*DSTRAN(2)DSTRESS(2)=E11*DSTRAN(2)+E12*DSTRAN(1)DSTRESS(3)=EG*DSTRAN(3)CCUPDATESTRESSCDOK=1,NTENSSTRESS(K)=STRESS(K)+DSTRESS(K)ENDDOCCDETERMINEJACOBIANCDDSDDE(1,1)=E11DDSDDE(2,2)=E11DDSDDE(1,2)=E12DDSDDE(2,1)=E12DDSDDE(1,3)=DDSDDE(3,1)=DDSDDE(3,3)=EGCCRETURNEND**

*********************************************************************VUMATFORABAQUS/ExplicitINCORPORATINGELASTICBEHAVIOUR****FORSHELLELEMENTSunderPLANESTRESS**********************************************************************USERSUBROUTINEsubroutinevumat(CReadonly-1nblock,ndir,nshr,nstatev,nfieldv,nprops,lanneal,2stepTime,totalTime,dt,cmname,coordMp,charLength,3props,density,strainInc,relSpinInc,4tempOld,stretchOld,defgradOld,fieldOld,3stressOld,stateOld,enerInternOld,enerInelasOld,6tempNew,stretchNew,defgradNew,fieldNew,CWriteonly-5stressNew,stateNew,enerInternNew,enerInelasNew)Cinclude''dimensionprops(nprops),density(nblock),1coordMp(nblock,*),2charLength(*),strainInc(nblock,ndir+nshr),3relSpinInc(*),tempOld(*),4stretchOld(*),defgradOld(*),5fieldOld(*),stressOld(nblock,ndir+nshr),6stateOld(nblock,nstatev),enerInternOld(nblock),7enerInelasOld(nblock),tempNew(*),8stretchNew(*),defgradNew(*),fieldNew(*),9stressNew(nblock,ndir+nshr),stateNew(nblock,nstatev),1enerInternNew(nblock),enerInelasNew(nblock)Ccharacter*80cmnameCCparameter(zero=,one=,two=,1third=/,half=,op5=ce=props(1)xnu=props(2)twomu=e/(one+xnu)e11=e/*xnu)e22=e*xnu/*xnu)Cif(stepTime.eq.zero)thendok=1,nblockstressNew(k,1)=stressOld(k,1)1+e11*strainInc(k,1)+e22*strainInc(k,2)stressNew(k,2)=stressOld(k,2)2+e22*strainInc(k,1)+e11*strainInc(k,2)stressNew(k,4)=stressOld(k,4)+twomu*strainInc(k,4)strainInc(k,3)=-xnu/1*(strainInc(k,1)+strainInc(k,2))enddoelsedok=1,nblockstressNew(k,1)=stressOld(k,1)1+e11*strainInc(k,1)+e22*strainInc(k,2)stressNew(k,2)=stressOld(k,2)2+e22*strainInc(k,1)+e11*strainInc(k,2)stressNew(k,4)=stressOld(k,4)+twomu*strainInc(k,4)