资源描述:
《子程序subroutine vumat》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、subroutinevumat(1k4A$K-v.c3F!YCReadonly-9n)Z5Z3D3v*C'
2、0F"F 1 nblock,ndir,nshr,nstatev,nfieldv,nprops,lanneal,.n,O._1b*c#K& 2 stepTime,totalTime,dt,cmname,coordMp,charLength,'p3k+]:C(b 3 props,density,strainInc,relSpinInc,8w8
3、/W&?0}4K4B/x-T#] 4 tempOld,stretchOld,defgradOld,field
4、Old,1m'f7r-z"f;O'f7p&K9i 3 stressOld,stateOld,enerInternOld,enerInelasOld,!a.s+r7]%Z&s;F.r"e 6 tempNew,stretchNew,defgradNew,fieldNew,;J#H(}&}-s!yCWriteonly-0p;^0N1}&^"Z4j(V 5 stressNew,stateNew,enerInternNew,enerInelasNew)+_6e1m%}-K/B$f include'vaba_param.inc'%R#b#z%+e;_8m8C"S$z
5、!E3Q6T+_'Z-PC!~*Q:i-t;~!V6b'LC MODELFORTHESHEARDEFORMATIONANDCRACKOFPOLYMERS!R$M/o!C-Z+A5v'C SHEARMODEL:SPRING+BURGERSELEMENT+PLASTICELEMENT(D-P&EYRING)$`%l#O'g'q8m7i3_C CRACKMODEL:S11>S11_CRWHENSKK>0'A9p!B5{/U+N!IC :W0T/J3F.G3R9]6hC WRITTENBY ZHANGCHUNYU(CHUNYU@NUS.EDU.SG))b6w3e.}(yC ,N
6、3I)[2H8Z)y4l*u2ZCAllarraysdimensionedby(*)arenotusedinthisalgorithm+'S"s3R p,y6G+T -t%X(r$m;e-T,P!h dimensionprops(nprops),density(nblock),$p5@!t8i2m9e'C 1 coordMp(nblock,*),)F/P/@!Y)R9O.c 2 charLength(*),strainInc(nblock,ndir+nshr),0Y"d-l9n:D.P%S$^'r 3 relSpinInc(*),tempOl
7、d(*),8~9h#H/w6Q.G:A 4 stretchOld(*),defgradOld(nblock,ndir+nshr),"T(W"_$O5P)g2U"~7X 5 fieldOld(*),stressOld(nblock,ndir+nshr),$F,e.n/z6[ 6 stateOld(nblock,nstatev),enerInternOld(nblock),*d6H!A,Q#P 7 enerInelasOld(nblock),tempNew(*),,m6P8K3K/f 8 stretchNew(*),defgradNew(nblock
8、,ndir+nshr),fieldNew(*),$K,H!N3;2L#` 9 stressNew(nblock,ndir+nshr),stateNew(nblock,nstatev),7N2b Q"f-H%*c4H8a 1 enerInternNew(nblock),enerInelasNew(nblock)'l8}$V!?0TC straincomponentsstoredasstatevariables!P&I)N2b4C,`'h!c&H#F dimensioneelas(ndir+nshr),eplas(ndir+nshr),deplas(nd
9、ir+nshr))c;z9y3e"N/K*^1j5h dimensionevisco(ndir+nshr),devisco(ndir+nshr),veint(ndir+nshr)$y,c"J/Q2z dimensioneffstrn(ndir+nshr),ps(ndir),an(ndir,ndir),s(ndir+nshr)8r2x!A5F2} dimensiondr(3,3),f_old(6),f_new(6)8r&Q+z#