子程序subroutinevumat.docVIP

subroutine vumat(1 k4 A$ K- v. c3 F! YC Read only -9 n) Z5 Z3 D3 v* C |0 F F? ???1??nblock, ndir, nshr, nstatev, nfieldv, nprops, lanneal,. n, O. _1 b* c# K \? ???2??stepTime, totalTime, dt, cmname, coordMp, charLength, p3 k+ ]: C( b? ???3??props, density, strainInc, relSpinInc,8 w8 |/ W ?0 }4 K4 B/ x- T# ]? ???4??tempOld, stretchOld, defgradOld, fieldOld,1 m f7 r- z f; O f7 p K9 i? ???3??stressOld, stateOld, enerInternOld, enerInelasOld,! a. s+ r7 ]% Z s; F. r e? ???6??tempNew, stretchNew, defgradNew, fieldNew,; J# H( } }- s! yC Write only -0 p; ^0 N1 } ^ Z4 j( V? ???5??stressNew, stateNew, enerInternNew, enerInelasNew )+ _6 e1 m% }- K/ B$ f? ? ? ? include vaba_param.inc% R# b# z% \+ e; _8 m8 C S$ z! E3 Q6 T+ _ Z- PC! ~* Q: i- t; ~! V6 b LC? ? MODEL FOR THE SHEAR DEFORMATION AND CRACK OF POLYMERS! R$ M/ o! C- Z+ A5 v \C? ? SHEAR MODEL: SPRING+BURGERS ELEMENT+PLASTIC ELEMENT(D-P EYRING)$ `% l# O g q8 m7 i3 _C? ? CRACK MODEL: S11S11_CR WHEN SKK0 A9 p! B5 {/ U+ N! IC? ? : W0 T/ J3 F. G3 R9 ]6 hC? ? WRITTEN BY??ZHANG CHUNYU(CHUNYU@NUS.EDU.SG)) b6 w3 e. }( yC? ? , N3 I) [2 H8 Z) y4 l* u2 ZC All arrays dimensioned by (*) are not used in this algorithm+ \ S s3 R??p, y6 G+ T? ?? ? - t% X( r$ m; e- T, P! h? ?? ???dimension props(nprops), density(nblock),$ p5 @! t8 i2 m9 e C? ???1??coordMp(nblock,*),) F/ P/ @! Y) R9 O. c? ???2??charLength(*), strainInc(nblock,ndir+nshr),0 Y d- l9 n: D. P% S$ ^ r? ???3??relSpinInc(*), tempOld(*),8 ~9 h# H/ w6 Q. G: A? ???4??stretchOld(*), defgradOld(nblock,ndir+nshr), T( W _$ O5 P) g2 U ~7 X? ???5??fieldOld(*), stressOld(nblock,ndir+nshr),$ F, e. n/ z6 [? ???6??stateOld(nblock,nstatev), enerInternOld(nblock),* d6 H! A, Q# P? ???7??enerInelasOld(nblock), tempNew(*),, m6 P8 K3 K/ f? ???8??stretchNew(*), defgradNew(nblock,ndir+nshr), fieldNew(*),$ K, H! N3 \; \2 L# `? ???9??stressNew(nblock,ndir+nshr), stateNew(nblock,nstatev),7 N2 b??Q f- H% \* c4 H8 a? ???1??enerInternNew(nblock), enerInelasNew(nblock) l8 }$ V!