ABAQUSVUMAT初學(xué)者用戶子程序小例子

1、Oeser ipon:GehavisUse*上一次發(fā)過ABAQUS初學(xué)者用戶子程序小例子，給學(xué)習(xí) UMAT的初學(xué)者帶來了一定的幫助?，F(xiàn)在用到VUMAT,發(fā)現(xiàn)網(wǎng)上 這種小例子很少，關(guān)于VUMAT的資料也不多。摸索了一天，做個 VUMAT的小例子供大家分享。實(shí)例：簡單的平面平板拉伸，材料本構(gòu)模型采用隨動強(qiáng)化模型,E=210000MPa,尸0.3,樂=200MPn, E'=10000MPa。左端約束,右端施加位移載荷V = 0.2mm。Author: xueweek Y1材料本構(gòu)模型：在 Property 中定義材料時，在 General 卜選中 User Material.輸入 21000

2、0、0.3、200、 10000,以上兩個數(shù)值代表E、v,、人、E'。在用戶子程序中代表看PROPS、PROPS(2)、 PROPS(3)、PROPS(4)。然后在Gniwl下選中Depvar,由于該例子中使用了五個狀態(tài)變量 (背應(yīng)力張彊和累積塑性應(yīng)變變量)，因此在第一項(xiàng)中輸入大于5的數(shù)值即可。另外還需要 輸入密度。fienel MecMesl bermai QttoDcMUs lAawnalUw«type MxMriol日B Uu vnc»«nMitric etWMxhaniolConstants1 210000290c,Dmu：Edit SilWp M

3、ncrialMst-al “hzo”Variable rtimber comrollin clement dete<<fi Sbg"/t*pi c4o 刃.Njm»cr c< wiv6x dc«xdcrtOxriptio 代MateUI Behftticrs*xrd Mzh.rM：dl Rcrmd QtherUm tootfitu c deptoitnt do3CMNumber H <<Hd “HoBc：Dataaoi建模人家都會，故省略2 ABAQUS中STEP的設(shè)置由于VUMAT需要用到Explicit求解，因此需要在step步驟

4、中設(shè)置explicit選項(xiàng)，如下圖, 其設(shè)置可以用默認(rèn)設(shè)置。 Create StepName: Step-1Procedure type: GeneralInsert new step after#Direct cyclic Dynamic, ImplicitDynamic, ExplicitDynamic, Temp-disp, ExplicitGeostaticHeat transferMass diffusionCancelSoilsContinue.3 ABAQUS調(diào)用VUMAT用戶子程序同UMAT用戶子程序的調(diào)用方法。在Job Manager中點(diǎn)擊Edit選項(xiàng)，在General選項(xiàng)的

5、最后一項(xiàng)中選擇自己建立好的用戶 子程序文件。（注：用戶子程序文件可以使用文本編輯器進(jìn)行編輯，當(dāng)然也可以用Fortran 編譯器，如果對自己的用戶子程序文件的語法不放心，町以先用Fortan編譯器進(jìn)行編譯， 不過編譯前要先建立Project,關(guān)于Fonran編譯，這里不再介紹）。完成后，點(diǎn)擊submit即可進(jìn)行分析。3(Avg: 75%)4 結(jié)果以下兩張圖分布是用戶子程序和ABAQUS自帶的材料模型（Standard求 解器）得到的應(yīng)力云圖，可以看出兩種圖形基本相同。+2.4970*02 *2.4519*02+2.3619*02 4-2.31bO*U2 *2 2700*02 +2.2250*02

6、 +2.1&00*02 4-2.1350*02 +2.0590*02 +2.0449*02 +1.9990*02 *1 9540*02ODB:ld» .ka<w9Txplkii6.1O 1 Mon 0ct31 21u)9:J0 GMTM)8:(M) 2011V 3tep 9tep-l八 Ineraneni 3859 SSp T“ I KOPrsnary Vrn S, LlirsDeformed V ar U E efznnaLion Scale Faclcu +2 l2!e0l#3/ Mises (Avg: 75%)*2.4459-4-02 +2.4030+02 42.

7、3600402 k 十2.31%十02 2.275e-t02 乜233M02 - 4-2.1909*02 k 7 148902*2.10594-02 4-2.0639+02 42.02004 02 tl.978B-tO2ODB: .Job.2>A Ab2<iiKnndir4(fjn.l Men Orf 3】21:48:49 GMT*fl« 00 2011X SltTi StefuiIncrfTiMit I:Step丁m I 000Pnmarz V u： S, Mi?"Defotriei V v: U D GfotnaUn S c £e F actor:

8、 *2.1259-HJ15 VUMAT子程序?qū)τ诔鯇W(xué)者來說，需要注意的是，F(xiàn)ORTRAN對干程序語言格式上的要求。例如，對于 FORTRAN語言，前六個字符必須空出來，等等。檢査語法最好的方法就是在FORTRAN編 譯器上進(jìn)行編譯。SUBROUTINE VUMAT(C Read only -1 nblock, ndir. nslu; listatev, iifieldv, npiops, lanneal,2 stepTiine, totalTnue, dt, ciniiame, coordMp. chaiLength.3 props, density, straiiilnc, relSpmli

9、ic,4 tempOld, stretchOld, defgradOld. fieldOld,3 stiessOld, stateOld. eneilnteinOld, eneilnelasOld,6 tempNew. suetcliNew, defgiadNew, fieldNew, C Write only 5 stressNew, stateNew. eneiliiternNew, enerliielasNew)cmclude ,vaba_paiam.inc,CC J2 Mises Plasticity with kuieinatic liaidenmg for planeC stmin

10、 case.C Elastic predictoi; radial coirector algorithmCC The state variables are stored as:C STATE(*J) back stiess component 11C STATE(*,2) = back suess component 22C STATE(*3) = back suess component 33C STATE(*,4) = back suess component 12C STATE(*,5) = equivalent plastic strainccC All arrays dimens

11、ioned by (*) are not used in this algoritlun duneiision props(nprops). density(nblock),1 cooidMp(nblock)、2 charLengtli(*), stiauilnc(nblock衛(wèi)diLiislii).3 relSpinInc(*), tempOld(*),4 stietchOld(*), defgradOld(*),5 fieldOld(*), suessOld(nblock,ndu+nslu),6 state01d(nblock,nstatev). enerlntemOld(nblock),

12、7 enerliielasOld(nblock), tempNew(*),8 stietcliNew(*), defgradNew(*), fieldNew(*),9 stiessNew(nblockjidirnslu), stateNew(nblocknstatcv),1 enerliiteinNew(nblock), enerlnelasNew(nblock)ccliaiactei*80 cmiiainecpaiameter( zero = 0., one = 1., two = 2,tliree = 3,1 tlurd = one/tliree, half = .5, twoThirds

13、 = two/three,2 tlueeHalfs = 1.5 )ce =props(l)xnu = props(2)yield = piops(3)haid props(4) Ctwomu =e/( one + xnu )tluemu = tlueeHalfs * twomusixniu thiee * twomualamda = twomu * ( e - twonw ) /( sixnm - two * e )term = one / (twomu * ( one + liard/thiemu ) coni = sqrt( twoTlmds )cdo 100 1 = 1 ,nblockC

14、C Tnal stiesstrace = strauilnc(ij) + strauihic(i.2) + suamlnc(i,3)sigl = stressOld(iJ) + alamda* trace + tromu *strauilnc(i, 1) sig2 = stiessOld(i,2) + alainda*trace + twromu*strainlnc(i,2) sig3 stiessOld(i3) + alainda*trace + tromu*stiamlnc(i,3) sig4 = stiessOld(i,4)+ turomu su auihic(i,4)CC Trial

15、stiess measured fiom die back stress si = sigl - stateOld(i,l)s2 = sig2 - stateOld(i,2)s3 = sig3 stateOld(i,3) s4 = sig4 stateOld(i,4)CC Deviatonc pait of tiial stress measured fiom the back stress smean = third * (si + s2 + s3 ) dsl = si - smeands2 = s2 - smeands3 s3 smeancC Magiumde of the deviato

16、nc trial stress differencedsmag = sqrt( dsl*2 + ds2*2 + ds3*2 + 2.*s4*2 )cC Check for yield by deternunmg the factor for plasticity.C zero for elastic, one for yield radius = coni * yield facyld zero if( dsniag radius .ge. zero ) facyld = one cC Add a protective addition factor to prevent a divide b

17、y zeroC when dsniag is zero. If dsniag is zero, we will not have exceeded C the yield stress and facvld will be zero.dsmag = dsmag + ( one - facyld )cC Calculated increment m gaimiia (this explicitly includes theC time step)diff = dsniag - radius dgaimiia = facyld term * diff cC Update equivalent pl

18、astic strain deqps = coni * dgaimna stateNew(i,5) stateOld(i,5) + deqps cC Divide dganuna by dsniag so that the deviatonc stresses areC explicitly convened to tensors of unit magnitude m theC fbllowmg calculations dgaiiuna = dgamma dsniagCC Update back stressfactor hard * dgamma * nvoTluids stateNew

19、(i,l) = stateOld(iJ) + factor * dsl stateNew(i,2) = stateOld(i,2)十 factor * ds2 stateNew(i,3) = stateOld(i,3) + factor * ds3stateNew(i,4) = stateOld(i,4) + factor * s4cC Update the stiessfactor = twomu * dganmia stiessNew(i,l) = sigl - factor * dsl stiessNew(iJ) = sig2 - factor * ds2 stiessNew(i3) =

20、 sig3 - factor * ds3 stressNew(h4) = sig4 factor * s4CC Update the specific intenial energy stiessPower = half * (1( stressOld(iJ)+stressNew(i,l) )*stranilnc(i, 1)1+( stressOld(i,2)+stressNew(i,2) )*strauihic(h2) 1+( stressOld(i3)+stressNevv(i,3) )*stiaHiIiic(i3)1+ hvo*( stiessOld(i.4)+stressNew(i,4

21、) )*stiauilnc(i.4)cCeneilnteniNew(i) = eneiliiternOld(i)C 1+ suessPower / density(i)CC Update the dissipated inelastic specific energy plasticWorklnc = dgainina * half * (1( stressOld(i, 1 )+stiessNew(i, 1) )*dsl1+(stressOld(i.2)+stressNew(i,2) )*ds21+( stressOld(i3)+stressNew(i3) )*ds31+ two*( sues

22、sOld(i.4)+stressNew(i,4) )*s4 )CeneilnelasNew(i) eneiliielasOld(i)C 1+ plastic Woikliic / density(i)100 contmuecxetuinend6 INP輸入文件*Heading* Job name: Job-1 Model name: Model-1* Generated by： Abaqus/CAE 6.10-1*Preprmt, echo=NQ model=NO, histoiy=NO, contact=NO* PARTS*Pait, nanie=Pan-l*Node1,10.,52,8.4

23、8214245,53,6.96428585,54,5.44642878,55,3.92857146,562.41071439,5乙0.892857134,5.&-0.625,5.9.214285707,5.10.-3.66071439.5.lh-5.17857122,5.12,-6.69642878,5.13,-8.21428585,5.14,-9.73214245,5.15,-11.25,5.16.-12.7678576.5.17,-14.2857141,5.18.-15.8035717,5.19.-17.3214283,5.20.8392849.5.21.-203571434,5.

24、22.-21.875,5.23,23.3928566,5.24.-24.9107151,5.25.-26.4285717,5.26.-27.9464283,5.27,-29.4642849.5.2&-30.9821434,5.29.-32.5,5.30.10.,3.53L8.48214245,3.532,6.96428585,3.533.5.44642878,3.534.3.92857146,3.535,2.41071439,3.536,0.892857134,3.537,0.625,3.53&-2.14285707,3.539.-3.66071439.3.540.-5.178

25、57122,3.541.-6.69642878,3.542.-8.21428585.3.543,9.73214245,3.544.-11.25,3.545.-12.7678576.3.546.-14.2857141,3.547,-15.8035717,3.54&-17.3214283,3.549.-18.8392849,3.550.-203571434,3.551.-21.875,3.552,23.3928566,3.553,-24.9107151,3.554.-26.4285717,3.555,-27.9464283,3.556.-29.4642849.3.557,30.982143

26、4、3.55&-32.5,3.559.10.,2.60.8.48214245,2.61,6.96428585.2.62,5.44642878,2.63.3.92857146,2.764.2.41071439,2.65.0.892857134,2.66.0625,2.67,214285707,2.6&-3.66071439,2.69.-5.17857122,2.70.6.69642878.2.7L-8.21428585,2.72,-9.73214245,2.73.-11.25,2.74.-12.7678576.2.75,-14.2857141,2.76.-15.8035717,2

27、.77.-17.3214283,2.7&-18.8392849,2.79.-20.3571434,2.80.-21.875,2.81.-23.3928566,2.82.-24.9107151.2.83,26.4285717,2.84.-27.9464283,2.85.-29.4642849.2.86.30.9821434,2.8人-32.5,2.8&10,0.589.&48214245,0.590.696428585,0.59h5.44642878,0.592,3.92857146,0.593.2.41071439,0.594.0.892857134,0.595,0.6

28、25,0.596.-2.14285707,0.597.-3.66071439.0.598.-5.17857122,0.599.-6.69642878,0.5100.8.21428585、0.510L-9.73214245,0.5102,-11.25,0.5103.12.7678576,0.5104.-14.2857141,0.5105.-15.8035717,0.5106.-17.3214283,0.5107,-18.8392849,0.510&-20.3571434,0.5109.21.875,0.5110,-23.3928566,0.5Hl,-24.9107151,0.5112,2

29、6.4285717,0.516-27.9464283,0.5114.29.4642849,0.5115,30.98214310.5116,-32.5,0.5117,10.,-1.11&8.48214245,L119,6.96428585,-1.120.5.44642878,-1.121,3.92857146,-1.122,2.41071439，-1.123,0.892857134,-1.124,-0.625,-1.125.-2.14285707,-1.126.3.66071439,-1.127,-5.17857122,-1.12&66964287&-1.129.-8.2

30、1428585,-1.130.-9.73214245,-1.13L-11.25,-1.132,-12.7678576,-1.133,-14.2857141,-1.134.-15.8035717,亠135,-173214283,-1.136.-18.8392849,-1.137,20.3571434,亠138.-21.875,-1.139,23.3928566,-1.140.24.9107151,亠14k-26.4285717,-1.142.-27.9464283、143,-29.4642849,-1.144.-30.9821434,-1.145.-32.5,-1.146.10.,-2.5147

31、,8.48214245 2514&6.96428585,-2.5149.5.44642878,-2.5150.3.92857146,25151,2.41071439,-2.5152,0.892857134,-2.5153.-0.625,-2.5154.-2.14285707,25155,-3.66071439.-2.5156.-5.17857122,25157,6.69642878,-2.5158.-8.21428585,-2.5159.9.73214245、25160.-11.25,2516k-12.7678576,-2.5162.-14.2857141,-2.5163、-15.80

32、35717,-2.5164.-17.3214283,-2.5165.-18.8392849,-2.5166.20.3571434,-2.5167,-21.875,2516&-23.3928566,-2.5169.-24.9107151,-2.5170.-26.4285717,-2.517L-27.9464283,-2.5172,-29.4642849,-2.5173,30.9821434,-2.5174.-32.5,25Element, tvpe-CPE4Rh h 2,31、302,2,3,32,3198,9,10.It12,13,14,15,16,17,1&19,20,21,

33、22,23,24,25,26,27,2&29,30.31，32,33,34,35,36,37,3&39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,5&59.3.4.33,324.5,34.335,6.35,346.7,36、35乙&37,36&9.3&37210.39、3810.11,40,3911,241,4012.13.42、4113,14.43,4214,15,44、4315.16,45、4416.17.46、4517,1&47,4618.19,4&4719.

34、20.424820.21.50、4921.22.51.5022.23,52,5123.24.53、5224.25,54、5325,26.55,5426.27、56、5527,2&57,562&29,5&5730.31、60、5931.32.61.6032,33,62,6133.3第63、6234.35,616335,36,65,6436.37.66,6537,3&67,663&39.68.6739.40,69.6840,41,70、6941.42,7L7042.43.72,7143,44,73,7244.45,74,7345.46,75,7446.47,

35、76.7547,4&77,764&49,7&7749.50.79、7850.5h80、7951.52.81、8052,53,82,8153,54、83,8254.55.84,8355,56,85,8456.57,86.8557,5&87,8659.60,89,8860.61.90,8961.62.91,901160.62,63.92, 9161,63.64.93, 9262,64.65,94、9363、65.66,95, 9464.66.67.96. 9565,67,6&97, 9666.6&69,9& 9767,69.70.99、986

36、8,70.71,100. 9969.71,72,10L 10070、72,73.102, 10171,73,74,103, 10272,74.75、104. 10373,75,76.105. 10474,76.77,106. 10575,77,7&107. 10676,7&79.108. 10777,79.80,109. 1087&80,81.110, 10979,81,82,111,11080.82.83,112, 11181,83.84,113, 11282,84.85.114.11383,85.86,115. 11484.86.87,116, 11585,88.8

37、9,11& 11786.89.90,119, 11887,90.91.120. 11988,91,92.121, 1208292,93、122,12190.93.94,123.12291,94.95,124 J 2392,95,96,125. 12493,96,97.126. 12594,97,9&127,12695,9&99,12&96.99. J100,129. 12897,100. 101, 130, 1299& 10U 102, 131, 13099,102,103, 132, 131100.103, 104, 133, 132101.104,

38、105, 134, 133102.105, 106, 135, 134103.106.107, 136, 135104.107, 10& 137, 136105.108, 109. 13& 137106.109, 110. 139. 138107.110, 111, 140, 139108,111,112,141,140109.112, 113,142, 141110.113, 114, 143, 142111.114.115, 144, 143112.115, 116, 145, 144113,117, 118,147, 146114, 118, 119, 148, 1471

39、15.119, 120, 149, 148116.120, 121,150, 149117.121.122.151.150118.122.123.152.151119.123.124.153.152120.124, 125, 154, 153121.125, 126, 155, 154122,126,127,156, 155123,127, 128,157, 156124,12 & 129. 158, 157125.129, 130,159,158126.130, 131, 160, 159127.131, 132,161,160128.132, 133,162, 161129.133

40、, 134, 163, 162130.134,135,164,163131.135, 136, 165, 164132.136, 137,166, 165133.137, 138, 167, 166134.138.139, 168, 167135.139, 140, 169, 168136.140, 141,170,169137,141,142,171, 17013& 142, 143,172, 171139,143, 144,173,172140,144,145, 174, 173*Nset, nset=_PickedSet 10, mtemal, generate1, 174,1*Elset. elset=_PickedSet 10. internal, generate1, 140、1* Section: Section-1Solid Section, elset=_PickedSet 10. controls=EC-l, ma