LS-DYNA材料模型及参数

1、材料模型及参数目录1基本的状态方程21.1EOS_JWL22.2 EOS_GRUNEISEN22.3 EOS_LINEAR_POLYNOMIAL32.材料模型32.1 MAT_HIGH_EXPLOSIVE_BURN3RDX5HMX5TNT51.2 MAT_NULL 5空气6水61.3 MAT_JOHNSON_COOK7紫铜8钢81.4 MAT_PLASTIC_KINEMA TIC9钢10高导无氧铜 10土壤101.5 MAT_STEINBERG 10高导无氧铜 121.6 MAT_JOHNSON_HOLMQUIST_CERAMICS 12B4C陶瓷141.7 MAT_JOHNSON_HOLMQ

2、UIST_CONCRETE14混凝土 143其它材料参数15LY12CZ铝合金15 一（回到目录）许鸿杰一整理收集15主要材料模型及参数1基本的状态方程1.1EOS_JWLViiriuhkEOSIDARR.3R2OMEGED'VOT>peIFFFFFFFTheWL equiition of $2出 defines the presfiure asand is usually u&ed for detonation products of high explosives.2.2EOS_GRUNEISENR/marko:The Grunei&en «quaf

3、jon. of state with cubic shock velocky-psutjcle velocity defines pressure for compressed maierjals asi-(e -iftrctiid for expanded materiials 抽where C js the intercept of thecurve: S 八 §2sl and Sj are the coefficieiiB of the slope of the vs»vp curve;芾卜 is Ehc Grunciscn garni ms; u i原 the fi

4、rst order valumc carrecliiQn lo 布 and 白 ± E - L2.3EOS_LINEAR_POLYNOMIAL(对EOS_GRUNEISEN进行线性化)L The linear pylynomia cqualiori. of sUlLc is lincfir in inlcmH energy. The pressure 冠 given by;尸=£, + 0* + G 犷 + J+ 仁+ C# + G * -where tenins 匚工口工 iund are to zero if <0, andl is the ratio of cu

5、rrentP口PiJdensity to initial density.2.材料模型2.1MATHIGHEXPLOSIVEBURNt iird KurmiilCilid JVariableMIDRODPCJBETAKGSIGYTypeFFFFI-FEJurn thictiflns, /, which mukiply (he equiUinns nf stares fnr high esplnsives, cnntrol rhe release ofchemjcal energy Rir simulating delonutions. Al ；uiy time, the pressure in

6、 a high explosive elrnitjit is given by;产"#kJ匕昌where /二，后 the pfeshurc Fruit】 the cquuliun u( lule (either types 2 ur 3j. V A the rekitjve volujue and E is (he huernal energy density per unir initial vnlume.In the inilializat ion ph use, j ifh tinp licTie t匕 computed for t LM?h e lenitnl by dix

7、idinp the distimce from the detenatien point to the center of the element hv the detonation velocity f). f multiple detonarmn points Lire defined, the closest deronLit ion poinr detemiinei The bum trjcdcn 7 is taken ns rhe niiuhnumwlTereMr)也.where 【:丫 依 die Ch upi n lui -J on g net re I at i v e vul

8、mtit iuid t is currenl tiJiir. If * exceeds 1, il is reset tt I. This CLikLilLUiun of rhe bum fniction usimlly requirts seveml time step% for 7 lo r? uch unity, thereby spreading the bum front over wieral elemeiirji. After reaching unit/ / is held constant. Th is burr traction cakualien is based an

9、work by Wilkins 1964 and is also discussed by Giroux 11973.If the beta burn upHcm g used* SETA=L0, uny volutiwiric cunipresion will caiid dctvtiatiun andiand £ im nm enmputed.f prngramniecl bum is used, BHTA=2.0, the explosive model will behave as an elastic perfetlly p lust it mate ri j I 什 th

10、e bulk E 口dulus* shciir Eudulus. .uid yield blrcy Lire <J±ined. Ttiere'ore» with this opden the explosive marerial can compress xvirheut causing detonation.As an option jhe hifh explosive marerjal c；in behave as an eLjstlc perteei ly-pasdc solid prior to detonation. Jn this eti!>

11、e we update the stereos tensor to an elcisdc trial stress= < +小再+ J： 口 *+之吃加where <r is rhe shear mcclulusT end 3 is rhe devinrnnc srrnin rate. The vnn Miyield enndirion is given by:I三where rhe seccnd stress invariant. Z. defined in terms ol the deviataric stress ccirponents js4二*的and the yiel

12、d stress is cf. If yielding has occurred, i.e0, the deviator ic trial stress is scaled to obtain rhe final dev j at curie stress ar time in-1:If <0uhenBeforedettwkiiion pressure in given h(y rheexprebsionwhere K is ilie bulk, nuxlulus. Once the explosive mjierhil delonaten：；md rhe material behave

13、s like a g口也RDX密度：1.69E+3 kg/m3; D: 8310m/s ; Pcj : 30.45 GpaA:850 Gpa： B: 18 Gpa：R” 4.6:R2 1.3: w°.38: a:10MJ/kgFor (g-cm-us)*MATHIGHEXPLOSIVEBURN1.698.3100.3015*EOSJWL8.500.184.61.30.3810 e-021.00HMX密度：1.891 E+3 kg/m3, D： 9910m/s. Pcj： 42Gpa,A: 778.3 Gpa： B: 7.1 Gpa：Ri ： 4.1;R2: 1.00; w。： 30;

14、E。： 1。5 MJ/kgFor (g-cm-us):*MA T_HIGH_EXPLOSIVE_BURN11.899.9100.420*EOS JWL17.7830.0714.21.00.30 10.5 e-021.00TNT密度：1.63 E+3 kg/m3; D： 6930 m/s； Pcj:27 Gpa：A: 371.2 Gpa； B: 3.21 Gpa；R1 ： 4.15;R2: 0.95; w0: 30 E0:4.29 MJ/kgFor (g-cm-us)*MA T_HIGH_1*EOS JWL_EXPLOSIVE_1.63_BURN6.9300.27013.7130.07434.1

15、50.950.307.0 e-021.001.2MAT_NULLGiid IVariabkMIDRDPCMUTERGCCERODPRType1FFFFFFFDefauhsnonenonen.Q0.0口加0.CIo.no-jnThe null materiuil iinist be u久d with an equLiliun or-state Pressure cutofF is negative in tension. A (deviator ic) viscous ntrei?s of the formi 岳 compuied for non zero p u here 匕(he devit

16、iiciric striiLn ratt 口 is the dyiumiic vhcosity with unit of Pascals.econd-2. The null rnateria h迎覆 no shear stiffhess and huurgl口写 control rniList be used with g;netit care. In some upplicuitiurii£"hu dl 匕 fiiullt hourgjiiss coerUcient might le；id tu igniriCciiit -energy I mses- In genera

17、l for fluidfX the hcurglass coefficient QM should be small (in the nmge LOE-4 tn 1.0E- 6 in the SI unit sysletn for the ytiintkird deRiult 1HQ chuicr).3. Th-e Null Hiuieria tia$ mu yield sltdiiglh. mid behiivesi in a flu id-like munri栏 r,4bThe pressure HiitpiT. PC5, ran st be tie fined to n low fur

18、u miilcTiiL】 to b*nidiiierkLillyo cavitute. Inother words, when a material undergnes dilararinn above certain magnitude, it should no longer be ubk to rtfsis.t this dilalatLon. Sinc dil.LiiLiE.ion stress or pressure is nefLilive, ing PC limit to a verjr small negative number would allow for the mate

19、r in I to cavitate once the pressure in the mater till cnes below this neEatire value.空气*MAT_NULLRO=1.25 kg/m3, PC= -1.0pa (<0)MU=1.7456E-5 （动力粘性系数）*EOS_LINEAR_POLYNOMIAL0,0.4 , 0.4 , 01 , 0 Gpa, 0 Gpa,0 Gpa,253312.5.1.0For (g-cm-us)*MA T_NULL3 0.125e-02-1.0E-12*EOS GRUNEISEN1.749E-70000000000000

20、30.34440000000/* EOS_LINEAR_POLYNOMIAL0000000001.400000300000.40.402.5000E-61*MAT_NULLVo：iFor (g-cm-us):*MA T_NULL10.998*EOS GRUNEISEN-1.0E-110.8684E-5000000000000011.65001.92-0.096000000.35000001.3MAT_JOHNSON_COOKCanl 11345678VariableMIDROGEPRDTFVPTyptr.FFFFFFDtluullnonenon匚WETnuntnuncG.OU.OCura) 2

21、VririhlcAliNcTRHFSOTypeFFFFFFFFDernuiinone.QiO.Q0.0nonenoneTICfliCn*TficCard JVariableTypeCPPCSPALLITDID2DSD4FFFFFFFFDtfaiJldune2.00.00.00.00.00.0Johnson and Cook express the flow stress astr =+一whereA,日* 匚，n, and ni = input CDtiitanls芭产 effective plastic strainp*= effective plasiic strain rate for

22、Fo = I .v'1 £融T = hnniclngous temperature =Constants fora variety of materials are provided in | Jolmson iuid Cock L 9831. A fully viscuplastic finrrnulatino is nptinnnl fVPj which incorporates The ratt equatinns wirhin the yield uirfhce. An jdd lL iuniil cost is incurTed but the imprui fin

23、cnt in rrulti clui he dTLunut ic.Duet。nonlinearity in the dependence of flow stress nn plastic strain, an nccumre value nf The flew stress iieqij irrs 让eruiun for Lhe increment in plustic si ruin. However, by using T uylor series expaiisiun with line jrizutjon ubuut the current tindeh can solve for

24、仃¥ with iutTicient ticcur比y to avoid iteratioD.The strain 跳 fract Lire is given by紫铜EX=1.19*MA T_JOHNSON_COOK18.960000.460.900E-03 2.920E-030.310 0.250E-011.09 0.1356E+04210 0.100E-050.383E-05 -9.00E+003.000.003.000.000.000.000.00*EOS GRUNEISEN10.3941.4890.000.002.020.470.001.00钢EX=2.0*MA T_JOH

25、NSON_COOK29.960000.460.900E-03 2.920E-030.310 0.25E-011.09 0.136E+04210 0.100E-050.383E-05 -9.00E+003.000.003.000.000.000.000.00*EOS GRUNEISEN10.3941.4890.000.002.020.470.001.001.4 MAT_PLASTIC_KINEMATICVariuhhMIDK0EPREICYETANUETATy产IE-l1FFFDef圳 liAneiti>ne的£0.00.0GirJ 2Variuhle£KCSRPFSV

26、PTypeFFFFDefaultnol ufciInol useda used0.0XEnwk 写:Strain rale 诂 lie counted for using the Cowper and Symonds model which scales the yield Sitress with the factorwhere e is the 立ruin rate. A fily viscoplastic fomiuhtioni 区 Dpriomd which incorporates the Cowper and Symonds formulaticin within the yiel

27、d surface. An iddiriiomil cast k incurred but the improvement is results can be drainatjc- To ignore s-train rate eftects &et both SRC and SRP to zero.Kinematic, isotropic or a comb in inion of kmematDc and isotropic huirdening may be specified by varying , between 0 and L For equal to O and T r

28、espectively kinematic and isotropic hunleniig mre obtirneLl 泗 sliown hi Figure 20,2, For iuirnpic liartlenirig,邛 =I. Material Mud el I2B "MATJSOTROPlC_ELASTI.C_PLASTIC, requires fe&s EtorHge and is more efficient. Whenever possible M ater ial 12 k recoinmen.ded for sol id e lenientsT but fb

29、r shel l elements it is 扪力 accurate and. thus material 12 is not rec am in end in this case,Flpjiirr 20-2. Ellie-pliistic behiivior 5 ith kineniLilic and inotropic hardening where 巾 w】ti I Lire unde formed and deformed lengths of uniLixil tension specimen. k the ilope of the biliiiear stress wiriiin

30、 curve.*MA T_PLASTIC_KINEMATIC47.832.070.300 0.400E-025.00E-021.000.000.000.00高导无氧铜*MAT_PLASTIC_KINEMATIC1 8.931.170.350000 0.400E-02 0.100E-021.000.000.000.00土壤*MAT_PLASTIC_KINEMA TIC1 7.80 2.1 0.300 0.023 0. 0240.1.000 0 01.5MAT_STEINBERG*MAT_STElNBEttGThis is Materia Type IL This material is avai

31、l lab le for modeling muteriai，deform ing at very high strain加士,(>IO5) Lind can be used with solid elemenu. The yield strength is 曲 Fuinction of temperature und pressure. An equation of Mutu determines Ihe pressure.Canl IVariableMIDROWSI0BETANGAMA叩MTyiK1FFFFFFFCunJ 2VuiablcBDFHFATMOGAMOSAT刑FFFFFF

32、pFViriuNcPCWALLRFFLAGMMWMMXECOECI丁伊FFFFFFFFCard!VKri：iN$BC224就5EC6印T掰REOTyi*FFFFFFFFMIDXfaneruJ ld«ntJtkaim. A umqlie number hgio be.ROMLks.s density.CiDB.ask sheur rmxlulus.SCO* 亮 e def in In* equ at lons.H ETA仅 3蚪 JeflH ISl g «L|UU|UI.N% hr dcftfiing equ 皿iovh.GaAIAyi, iilil iiiil plusl

33、ic s Lruiin1 十的 deriri in£ «qu ill icicli.5KjM务m fee d-efining equntians.FCivme w ide Tujjll？-l ,e*30)EPALLS|xiJI type:bq. 0.0； dctakMcwjr.SQ-1XI! P.Pm，EQ 工让 I r n-nai 舌 fuij eknw-nt spalls and tens»nr p < 0( u never -nlkiwKli.白3 3.D: p < -fc：ui ri-ietneni 即Nk uM<, Js euver u

34、lhwed.RPRr. I fit, # OXX A k dcH: detloed.FLMjSrl1,0 fcr Jl cmHkmils fwthe co-kl ccrupress ic>n energy fiL Drfuuk hMMNlivn or Hmiii- Option al p cr Tl| rnioinium vulud.klkiXLncn Hrnu- UptliMkDl 肛4 mOKriTlUm Vxilue.RCh”ECOIJ corhprliin mtH忖险的八面间画#j,BbP 际e 加fin i n串 网u iifiams.BPb； 5ec drFin in g e

35、qumi 口口 加Hh石虺 dcrWi j.|£, equ ufiicwi L.FC 盆把 defiHliag 叫 uui Im £.AAtom it weight (if 0*j0P R nwM be defined hTMOTj11sx4*诵初 芸 fqUuHi,抗M,GAMO%, see defining equal ion 筝.SA遍 sw defining eqjwifliws.吐Prai w -fff (de fwu 11®-1 ,c+30)Users who have an interest in this mcslel are encouraged

36、 to study the paper by Steinberg and Guinan which provides, the theoretical basis., Another useful reference 悟 the KOVEC user's m tuiuiiLn lerms of the foregoing: input pi,ir;iiniieter5. we define the she lit modulusi G. before the miiteriaJ melts 皿：where p is the pressure, v is the relative vol

37、ume.民 Is rhe cold compression energy:900 * cxp(rtv)X = - VhMd Eg is the undting energy'Emis)=耳i 幻 +3F4rTll1ixFwhich is in terms of tht melting emptraliirt Tmix lLind the niellin lempcrutnrc til p - pn,Tmn .In the uhove equation KJ is defined byAwhere R us the gus const Lint iind A is the Eiionii

38、c weight If Rr is not defined* LS YN A c ompu te s it with R in. the cm-gnim -jjiicrosecond system of units.The yield strengrh Q,除 given by:(7 = %： + tip V' -A -auollr ' " -* 'I 3* J高导无氧铜*MAT_STEINBERG28.930.4770.120E-02 36.00.4500.000.640E-022.832.830.377E-03 0.100E-02 63.50.179E+0

39、4 2.021.50-9.003.000.000.000.000.00*EOS_GRUNEISEN20.3941.490.000.000.000.000.000.000.000.002.020.000.000.000.000.4700.001.001.6MAT_JOHNSON_HOLMQUIST_CERAMICSThis is Material Type I 10. This Johnson-Holmquist Plusricity I Jamap? Medel is iietnl for modeling, centmits. glas and other brittle imaterial

40、s. A more dehiikd descripth>n can be found jn a puper hy Juhnson Lind Hultiiquiiil 11993).CiEid IVaiiibkMIDROGABCMNT*1FFFFFFFVuniiHcEPSITSFMAXHEI.PHELBETAT心FFFFFFCiinl 3V nn ableDI132KIK2K3F5FFFFFFMIDMaflerljL kdeniillcuiLun. A unique number lin m be clhuiem.EPSIRefiereTice Urain山七RODen的1Mi瞩 imun

41、i ten tile srnen£lh BGSli 阳 mEulus加MAKhfuxijiium nBimaliedl frisclunal Klrmlh. (if Eq.nildifr；iultE to 1 e2Oi.AIviiuKl nnmiBliT'ed Woi骐h ExunmetfrHELHugEiiUI 1dgLie LiiniE.kFractured rwimiuJi 才ed sljenglh pnnmidt-erFUEL片0 mpi而上ru ul1he Hugorn 虞 «las it llmh.CSlreiiiglh puriiLKter | for

42、 sLruin rul<e depcndnKif)BETAFnK-tKMi of clastic energy loss con varied to hydruW 印仃即；MIrhiLMUrwJ uifvngih pu/u£nier (prehaure eKpubetic)DIRwnmMcr for pl口hIk 5H-ni nto frhiTt.Nlimacl s-trench pariimctei (pres-sure Mponcni.IBP4inictcr for plnsl k; sir41i.11 m HtjcEtcKIFirst pre玉*uru 匚口匕fkifri

43、l (tKjuivuiknt ta the bulk moduluEl>.K2Snrurid pnfssure CDeiririmil.K3Hl.lc coriuj&K |k.lbulk mcdulUhFSRiilurt triteriiLFS <0 0 hiH if p" + tB <0 Hensile fnihjgh FS < no failure idefiujH).FS> 口 fail if Itie 师min > FS.R enuirlcs:The equivalent stress lor aceraniic-type m ut

44、eri al is given by<T = <7； -"行： 一#；wheretr； =+，/(l + fin e)represents the jnt3cL undamaged behaviour,D=represents the uecu mu luted damage bused upim the increase in plustiv strain per computation al cycle uoid the plastic stmin to fracture于.二洌产4户)M小面represents rhe damaged behaviouR n each ca&e, the i h，indicates ai normalized quondEyT the stresseji being norma