《abaqus材料子程序》由会员分享，可在线阅读，更多相关《abaqus材料子程序（8页珍藏版）》请在金锄头文库上搜索。
1、 各向同性材料损伤本构模型SUROUTINE UMAT(SRES,SAE,DDSDE,SE,PD,S, R,DSDDT,DRPLDE，DRPLT， + STN，DSRN,TIME，DIME，TEMP,DTP,PREDE,DPRE,CMNME, + D,HR，NTEN,NSTATV,RS,NPROS,CORD,DROT，PNEDT， + CELET,DGRD0,DGRD1,NL，NPT，LAYE,KPT,KSTP,KNC） INCUDE ABA_ARAM.INC CHRACTR0 CNAME DIMNIN SRS(NTEN),STATEV（NSTATV), + DDSDDE(NTENS,NTE）
2、,DDDT(NTNS）, + RPLDE(NTE）,STRN(NTENS)，DSTAN（NES), +IME(),PREDEF(1),DPED(1),POS(NPROPS), +CORDS(）,DR(3,3)，DGRD0（,3）,DGRD(3,) DIMENSON TRANT(6)，TSTRANT()，P() DIMNSN LD_STRES(6） IMENSION DOLDSTRESS(6),DSTRES(6） DIMENSN C(，），D(,),STRESS(6）,BSRESS(),ROO(3), + FMND(6），DDMD()，DDM(6,6),ATEP(6）,ATEMP2(） ARAM
3、ER (ZE=0.,ONE1.D0,TWO=2D0,OUR=.D,HALF 0.5D0)C tartC IF（NPROS.LT.)THEC WRIE(,*)* RROR: MT REQUIRES *PROP=2 SOPC NDF E =POPS(1) V1 =POPS（2） G12 =PROPS(1)/TWO/(ONE+POS（2）)C Critcal vaues f stesse =PROP（) CPRPS(4) XS=OPS(5) GXPROPS(6) !rure energyn matrix TA=0.01C Curent rain DOI =1, NTENS STRANT() =STR
4、(I) + DSR(I) END DO Stiffness I =1,6 O= 1, 6 （I,J)ERO D DO END O ATEMP = (1+V12）（-WO*V12） (1,1） = E1*(1V2)/ATEM C(2,2) 11*(1V12)ATMP （3,3)=E11*(1-V2)ATEP C(1，2) E11*V12/TM C(1,3) = 11*12/AEP C（2,3)=E11*V12/AMP C(4，） 12 C(5,5)= G12 C(6,6) = G12 DO I =2, 6 D J ,1 C(I,） = (J,) ED DO ED DOC Critialvlues
5、 srains XET=XT/(C(1,)-2*V1C(1,2) ECXC(1，)2*V*(1,)） XES=X/C(4,) MOD SAEV(1) CStrn initiatincition A11 = TNT()*T+S(2)*TOSRNT(3)*TWO A12 1 / XT / XC A1 TAT(1)+STRANT(2)STAN(3) A2 =(XEC - XT）/ XEC /* 21 A31 STR()*TWO+STRA（）*TWO+SRAN（6）*O A2 = 1 /XE*TO = 12+ +A2 B11= SAN(2）*TOC 12= B11/E/XC B21 = TRANT(2
6、） B22 = (XEC ET) XC XET * B1 B3 = SRANT(5）*TWO 32 1 XES*TW 1= B1 + B22 + B3 11 = STANT(3)*TWOC C12 = C11/ XE XEC C1 SRNT(3) C22 = （XEC - XET)/ ET C1C 1=STRANT(6)*WO C32 = C / XES*TWOC C1 1 C22 C32 STTEV（2）=AC STATE（3)=B1C TTEV(4）=C FM ZEO F (A.G.EO) THN MN=SQT(A1）C F (1.TONE) HEN FMN =M+SQT(B) IF(C1
7、.GT.ONE) TH C NFMN+SQRT(C) C END F NDIF EDI STTV（5)FMN C wre(*,*) FM DM =EO DMDFM= ZE DO ，6 DFMDE(I） ZER DMDE() =ER EN DO F (FN .T. ONE) TENC CCULATE DM， DDMDFMNC WRE（,*）FM T1 = (C(1,1）-2V2*C(1,2) XE*2*CELENT /GX T = (ONE FMN) T DM= - XP(2)/FN RITE(6,*）T ,T1, T2, , M,C wrte(,*)DC CACLATE TEDERIVATIV
8、EO AMAE VARBLEWITHET TOFAILUEC ITERIO DMDFN = (N / MN T1) * (ONE - D）C CALCULAEDFMNE ( .GT. MLD)TEN DE(1) = AFMN*(WO*STRAT(）EC-T)XT/XEC DDE(2)= HALFMN*（WSTNT()+XEC-XET)/XETEC FDE(）= HALF/FMN*(TO*STRANT（)+X-XE）/XET/EC DFMN(4) = ON/FMN*TW*STRAN(4)/XE*W DFMND(5） =ONE/FMN*TWO*RANT(5）/XS*TW FMND() ONE/FM
9、TWO*TRAN(6)XES*TWO DI =1, DME(I)=DFMNDE(I） * DDMDFN N DO ED I EDIF D =MAX (M, DMD)C write(6,*) DM SAVE THEOLD RESS TO O_TR DO I =1, NTENS OLD_STRES(I） = STRES（I) D C Effctiv stiffnss DO I =1, D J = 1,6 D（I,)C(I,J) E DO ND DO F（DM.NE.ZER） TEN C(1，1） = (NE DM)*C(1,） CD(1,) = (OE- M）（1，2) CD（2,) = (1,2) C(2,2) = (ON - DM)C(2,） C(1，3) =（OE D)*C(1，3） C(3,1) = CD(1,) D（2，3）= (O D)C(，3) C(3，2) = C(2，3) D(,） = (NE - D)*C（4，） C(5,5)=
