Abaqus子程序UMAT和VUMAT的例子

[e-mail]*******************************************************************************UMATFORABAQUS/STANDARDINCORPORATINGELASTICBEHAVIOURFORPLANE****STRAINANDAXI-SYMMETRICELEMENTS.*************************************************************************************************************************************************************USERSUBROUTINESUBROUTINEUMAT(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*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=0.D0,ONE=1.D0,TWO=2.D0,THREE=3.D0,+SIX=6.D0,NINE=9.D0,TOLER=0.D-6)CDIMENSIONDSTRESS(3)CC--------------------------------------------------------------------CSPECIFYMATERIALPROPERTIESCE=PROPS(1)XNUE=PROPS(2)CCCSETUPELASTICITYMATRIXCE11=E/(1.0-XNUE*XNUE)E12=E*XNUE/(1.0-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)=0.0DDSDDE(3,1)=0.0DDSDDE(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'vaba_param.inc'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=0.d0,one=1.d0,two=2.d0,1third=1.d0/3.d0,half=0.5d0,op5=1.5d0)ce=props(1)xnu=props(2)twomu=e/(one+xnu)e11=e/(1.0-xnu*xnu)e22=e*xnu/(1.0-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.0-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)stressNew(k,3)=1.0write(*,*)stressNew(k,3)=,stressNew(k,3)strainInc(k,3)=-xnu/(1.0-xnu)1*(strainInc(k,1)+strainInc(k,2))enddoendifCreturnend