umat二次开发超弹性本构

1、APPENDIXNeo-Hookea n Hyperelatic Material User Subrout ine This program is based on the derivation of hyperelastic material constitutive model in Secti on . A stress and strain relati on shipwas derived from the n eo-Hookea nhyperelastic material con stitutive model that is n ormally represe nted as

2、 the strain en ergy with s t r a inin v a r i a nts.s u b r o u t i neC Read only (unm odifiable)variables -1 n block, n dir, nshr,nstatev,nf ieldv, n props, lann e a I ,2 st ep Ti m e , t o t a I T i m e , d t, cmn ame,coordMp,charLe ngth,3 props, den s i t y , s t r a inlnc, r e I S p inlnc,4 t e

3、m p O I d , stretch Old , defgradOld, f i e I d O I d ,defg r a d New,5 stressOld,stateOld, enerln t e r nO Id, enerln e I a s O I d ,t e mp N e w , stretchNew, 6C Wr iteonly (modif i a b I e) v a r i i a b I e s7 stressNew, st a t e N e w , enerln t e r nN e w , enerln e I a s N e w )CincI u deCdim

4、e nsion props (n props), den si t y (n block), coo rd Mp (n block,*),1 c h a r L e ngth(n block) , s t r a i nlnc(n block ,n d i r + nshr),2 r e I S p inlnc(n block ,nshr),tempOld( nblock),4 defg radOld (n block ,n d i r + nshr + n s h r),5 fieldOld( nblock ,nf ieldv), stressOld( nblock ,n dir + nsh

5、r), 6 s t a t e O I d ( nblock ,n s t a t e v ) , enerln t e r nOld(n block),7 enerln e I a s O I d ( nblock), tempNew( nblock)IdNNew (nc k ,n d i rr)stressNew( nblc)c k , nd i r-+ nshr),st a t e N e w ( nblock,n st at e v )e n e rl ntn b I lockif(cmn(1 :.eq.VUMATO')call VUMATO (n block, ndin sh

6、r, n st at ev,nf ieldv, n props, lannealst ep Ti m e ,t o t aITi me,c o o r d M p , charLe ngthprops,d en c , r e I S p i n I n c76)2 st ep Ti m e , t o t a I T i m e ,d t, cmn a m e,coordMp, charLe ngth,stretch Old , defgradOld, f i e I d O I d ,s t a t e O I d , enerln t e r nO Id, enerln e I a s

7、O I d , stretchNew, def grad New, f i e I d N e w ,stressNew, s t a t e N e w , enerln t e r nN e w , enerln e I a s N e w ) else i f (cmn a m e ( 1 :nf ieldv, n props, lann eal,call VU M AT 1 ( nblock, n d i r, n shr, n st at e v,3 props, den s i t y , s t r a inlnc, r e I S p inlnc4 tempOld, stret

8、ch Old , defgradOld, fieldOld,5 stressOld, s t a t e O I d , enerln t e r nO Id, enerln e I a s O I d ,6 tempNew, stretchNew, defgrad New, f i e I d N e w ,7 stressNew, s t a t e N e w , enerln t e r nN e w , enerln e I a s N e w )endiendCsub rout inevumatOCReadonly* n block, n dir, nshr, nstatev,nf

9、 ieldv, n props, lann e a I ,* st ep Ti m e , t o t a I T i m e , d t, cmn ame,coordMp,charLe ngth,* props, den s i t y , s t r a inlnc, r e I S p inlnc,* t e m p O I d , stretch Old , defgradOld, f i e I d O I d ,* stressOld,stateOld, enerln t e r nO Id, enerln e I a s O I d ,* tempNew, stretchNew,

10、 defgrad New, f i e I d N e w ,CWriteonly* stressNew, st a t e N e w , enerln t e r nN e w , enerln e I a s N e w )CineludeCdime nsion coord Mp ( nblock,*), ch ar Le ngth (n b I o ck) , props (n props),1 den s i t y (n block), s t r a i nlnc(n block ,n d i r + nshr),2 r e I S p inlnc(n block ,nshr),

11、tempOld( nblock),3s t r e t c h O I d (n block ,n d i r + n s h r),4 defgradOld (n block ,n dir + nshr + n shr),5 f i eld Old(n block ,nf ieldv), stress Old(nblock,n d i r + nshr),6 stateOld( nblock,n statev), enerln t e r nOld(n block),7 enerln e I a s O I d ( nblock), tempNew( nblock),8 stretchNew

12、( nblock ,n d i r + nshr),9 defgradNew( nblock ,ndir + nshr + nshr),1 field New (n block ,nf ieldv),2 stressNew( nblock, ndir + nshr),stateNew( nblock, nstatev),3 enerln t e r nN e w ( nblock), enerln e I a s N e w ( nblock)Cdime nsion d e v i a ( nblock ,n d i r + n s h r), 1 BBar (nblock,4), stret

13、chNewBar (n b I o c k , 4 ) , in t v ( 2 )Cc h a ra c te r *80 c m n a m ep a ram e t e r (zero =, one =,t w o =jth r e e =o u r=,h a Ireal C 1 0 , D 1 2k,twomu,amu,alamda,hydro,v onM ises, maxShear,1 m i d S t r a in, maxPri ncipalStra in1 1 8Cin t v (1)=ndirin t v (2)=nshrCi f (n d i r .ne. 3 .or.

14、 nshr .ne. 1 ) t h e ncall xp lb_abq err(1, 'Sub routi ne VUMAT is impleme nte d ' Q. zero) the ndokblockt r a ce 1 s t r a inIn c(k , 1 )+ s t r a inIn c(k,2 ) + s t r a inIn c(k,3)tressNew(k,1)stressOld(k,1)rainIn c ( k , 1 )racerainIn c ( k , 2 )race+ twomu*stra inIn c ( k , 3 )+alamda*tr

15、ace1 stressNew(k,4)=stressOld(k,4)* +twomu*strainlnc(k,4 )C w r i t e ( 6 , * ) totalTime,k,defgradNew(k, 1),stretchNew(k,1), C 1stressNew(k,2),stressNew(k,3)5stressNew(k,4) enddod okblockC JACOBI AN OF STRETCH TENSOR (U is sy mm etric and in local axis)Cdet = stretchNew(k,3)*119 1 (stretchNew(k, 1

16、) * stretch New ( k, 2)-stretchNew(k, 4)*two)scale = det* * (-ONE/THREE)stretchNewBar(k, 1)=stretchNew(k, 1)*scalestretchNewBar(k,stretch NewBar(k3retchNewBar(k,2 ) = stretchNew(k,3 ) = stretchNew(k,2 ) * scale3 ) * scaleCALCULATELEFTY-GREEN TENSOR (BBBar(k,1) = stretchNewBar(k, 1 )*two+stretchNewBa

17、r(k, 4) * * t w oBBar(k,2)=stretchNewBar(k, 2)*two+stretchNewBar(k, 4)3)*B B a r(k,4)= stretchBar (k, 1 ) *stretchNewBar(k,4) +stretchNewBretchNewBar(kTRCALBBaLATERESStensr(k,1 ) +Ba2) +BBar(k,3)wo*(det- one)EG* (BBar(k,2)-stressNew(k,1 )=EG*(BBar(k,1 )-TRBBar/Three) + PR stressNew(k,2)TRBBar/Three)

18、 + P R stressNew(k,3) = EG* (BBar(k,3)-TRBBar/Three) + P R stressNew(k,4) = EG*4)C Update the specific in ter nalhalfstressPower1( stressOld(k, 1 )+stressNew(k, 1)s t r a i nln c(k,1)(stressOld(k,2)+stressNew(k,2)* s t r a i nln c ( k , 2 )(stressOld(k, 3) + stressNew(k,3)s t r a inIn c ( k , 3 )(st

19、ressOld(k,4)+stressNew(k,4)* s t r a inIn(k,4enerln t e r nN e w ( k )=enerln t e r nO ld(k)+ stressPower/ den s i t y ( k )C S t r a ins un d e r c o r o t a t i onal coord in a t e sstateNew(k,1)=stateOld(k, 1 )tr a inlnc(stateNew(k,2)=stateOld(k,2)tr a inlnc(2)stateNew(k,3)=stateOld(k,3)tr a inln

20、c(3)4)nine(k,4)Calcu latehydro = (stressNew(k, 1) + stressNew(k,2) + 12 0stressNew(k, 3)/3.-hydro=stressNew(k, k1 )dok1 = n dir+1, ndir + nshrstressNew(kvon Mises + d e v i a ( k , kn d i r +1 n d i r + n s h rC w r i t e ( 6 , *) totalTime,defgradNew(k, 4),stretchNew(k,4)C 1 ,defgradNew (k, 3), def

21、gradNew(k,4) ,defgradNew(k, 5)C,det,TRBBarC 1, stressNew( k, 1 ),stressNew(k,2),stressNew(k,3), st ressNew(k,4)CCFailureCriteriaCm i d S t r a in 二 stateNew(k,1) + stateNew(k,2) maxShear = sqrt(stateNew(k,1 ) m i d S t r a in )*2.2.)stateNew(k,4)* *(m i d S t r a inmaxPr in cipalStra inm i d S t r a

22、 i n + m axShearm a x P r i ncipalStrai nmaxShearm i d S t r a invonsubo u t i n eCReadv u m a t 1o n 1 v.n block n dir,n sh r, n s t a t e v ,n f i e I d v,n p r o p s, 1 anneal* st ep Ti m e ,totalT i m e, d t, c m n am e,c h a r L e ngthprop s ,dtempOl d,s t r e t c h O 1 1 d , df ieldOld2stateOl

23、d, enerln t e r nO Id, enerln elasOld, stretch New, defgrad New, fieldNew,* stressNew, st a t e N e w , enerln t e r nN e w , enerln e I a s N e w ) ineludedime nsion coord Mp ( nblock,*), ch ar Le ngth (n block), props (n props),1 den s i t y (n block), s t r a i nlnc(n block ,n d i r + nshr),2 r e

24、 I S p inlnc(n block ,nshr),tempOld( nblock),3 s t r e t c h O I d (n block ,n d i r + n s h r),4 defg radOld (n block ,n d i r + nshr + n s h r),5 f i eld Old(n block ,nf ieldv), stress Old(nblock,n d i r + nshr),6 s t a t e O I d ( nblock ,n s t a t e v ) , enerln t e r nOld(n block),7 enerln e I

25、a s O I d ( nblock), tempNew( nblock),8 stretchNew( nblock ,n d i r + n s h r),9 defgradNew( nblock ,ndir + nshr + nshr),1 f i e I d N e w ( nblock ,nf ieldv),2 stressNew( nblock, ndir + nshr),stateNew( nblock, nstatev),3 enerln t e r nN e w ( nblock), enerln e I a s N e w ( nblock)Cdime nsion d e v

26、 i a ( nblock ,n d i r + n s h r),1 BBar (nblock,4), stretchNewBar (nblock,4),intv(2)Ccmn a m e,three =character* 80param e t e r (zero =, one =, t w o =*four=,half real C 1 0 , D 1 ,ak,twomu5amu,alamda, hydro,vo nM ises(n dir3 .or. n s h rim plem ente d ' Q. zero) the ncall xp lb_abq err(1, 

27、9;Sub routi ne VUMAT ist r a ce 1s t r a inIn c(k , 1 )+ s t r a inInc(k,2)st r a inIn c(k,3)Old(k,1 )u * s t r a inIn c ( k , 1 )+ alamda*trace1rainracet rainraces tress4)tressNew(k,1)t worns t r a in In c ( k ,C w r i t e ( 6 , *) totalTime,k,defgradNew(k, 1 )3stretchNew(k,1 )3C 1 stressNew(k,2),s

28、tressNew(k53),stressNew(k,4)C JACOBI AN OF STRETCH TENSOR (U is sy mm etric and in localaxis)det=stretchNew(k,3 )*1 (stretchNew(k, 1 ) * stretch New ( k, 2)-stretchNew(k, 4)*two) sea=s tst r etst r eta r ( k,2 )=s t2 )*st r et3 )=s t3 )*st r eta r ( k,4 )=s t4 )*CC CALCULATE LEFT CAUCHY-GREEN TENSOR

29、 ( BiB Ba r ( k , 1 ) = stretchNewBar(k,t w o + st r et c h N ew Ba r ( k,C CALCU LATEB Ba r ( k , 2) = stretchNewBar(k,STRESS ten s o r2 ) * * t w o + stretch New Bar ( k ,4) * * t w o4) +stretchNewBar(k,2)*stretchNewBar(k, 4),3 )=3 ) *BB a r ( k , 4 )= =s t r e t chNewBar(k,TRBBar = BBar(k,1 ) + BBar(k,2