abaqus1用戶材料子程序

19-19-#EPLAS(K1)=EPLAS(K1)+THREE*FLOW(K1)*DEQPL/TWOEELAS(K1)=EELAS(K1)-THREE*FLOW(K1)*DEQPL/TWO150CONTINUEDO160K1=NDI+1,NTENSSTRESS(K1)=FLOW(K1)*SYIELDEPLAS(K1)=EPLAS(K1)+THREE*FLOW(K1)*DEQPLEELAS(K1)=EELAS(K1)-THREE*FLOW(K1)*DEQPL160CONTINUEEQPLAS=EQPLAS+DEQPLSPD=DEQPL*(SYIEL0+SYIELD)/TWORPL=PROPS(3)*SPD/DTIMECCJACOBIANCEFFG=EG*SYIELD/SMISESEFFG2=TWO*EFFGEFFG3=THREE*EFFG2/TWOEFFLAM=(EBULK3-EFFG2)/THREEDO220K1=1,NDIDO210K2=1,NDIDDSDDE(K2,K1)=EFFLAM210CONTINUEDDSDDE(K1,K1)=EFFG2+EFFLAM220CONTINUEDO230K1=NDI+1,NTENSDDSDDE(K1,K1)=EFFG230CONTINUEDO250K1=1,NTENSDO240K2=1,NTENSDDSDDE(K2,K1)=DDSDDE(K2,K1)+FLOW(K2)*FLOW(K1)1*(EFFHRD-EFFG3)240CONTINUE250CONTINUEENDIFENDIFCCSTORESTRAINSINSTATEVARIABLEARRAYCDO310K1=1,NTENSSTATEV(K1)=EELAS(K1)STATEV(K1+NTENS)=EPLAS(K1)310CONTINUESTATEV(1+2*NTENS)=EQPLASCRETURNENDCCCSUBROUTINEUSERHARD(SYIELD,HARD,EQPLAS,TABLE)CINCLUDE'ABA_PARAM.INC'CDIMENSIONTABLE(3)CCGETPARAMETERS,SETHARDENINGTOZEROCA=TABLE(1)B=TABLE(2)EN=TABLE(3)HARD=0.0CCCALSULATECURRENTYIELDSTRESSANDHARDENINGRATECIF(EQPLAS.EQ.0.0)THENSYIELD=AELSEHARD=EN*B*EQPLAS**(EN-1)SYIELD=A+B*EQPLAS**ENENDIFRETURNEND18?5?2UMATHT(包含材料的熱行為)SUBROUTINEUMATHT(U,DUDT,DUDG,FLUX,DFDT,DFDG,STATEV,TEMP,$DTEMP,DTEMDX,TIME,DTIME,PREDEF,DPRED,CMNAME,NTGRD,NSTATV,$PROPS,NPROPS,COORDS,PNEWDT,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)CINCLUDE'ABA_PARAM.INC'CHARACTER*80CMNAMECDIMENSIONDUDG(NTGRD),FLUX(NTGRD),DFDT(NTGRD),$DFDG(NTGRD,NTGRD),STATEV(NSTATV),DTEMDX(NTGRD),TIME(2),$PREDEF(1),DPRED(1),PROPS(NPROPS),COORDS(3)CCCOND=PROPS(1)SPECHT=PROPS(2)CCINPUTSPECIFICHEATDUDT=SPECHTDU=DUDT*DTEMPU=U+DUCCINPUTFLUX=-[K]*{DTEMDX}DOI=1,NTGRDFLUX(I)=-COND*DTEMDX(I)ENDDOCCINPUTISOTROPICCONDUCTIVITYCDOI