Steinberg-Guinan本构模型在平面冲击实验应用.pdfVIP

一一一．墨生!!!垡：鱼!!翌塑查塑垡型查垩堂壁圭塞堕主盟堕旦 胡晓棉潘吴朱建士 北京应用物理与计算数学研究所 1引言 材料的弹塑性本构关系是联系着应力与应变、应变率、温度以及变形历史等之间的关系。为了准确 地描述材料的动态力学行为，必须全面考虑应变率、加载历史、高温、高压、大变形以及内部损伤的作 用。适合于数值计算程序的本构模型必须能准确描述材料的力学本构行为及相应的物理过程。对此在长 期的研究中，多种本构用来描述金属材料的动态响应。常用的模型一般分为三类：经验模型、物理模型 (-Lund)模型。 研究表明，无论那一种本构模型都有其自身的特点和优势。我们研究的目的决定了对本构关系的研 究的特殊需要：既要达到一定的计算精度，又要能获取特殊材料的相关参数，同时还要易于在总体程序 中实现。有些模型虽然对某些材料的计算可以精度相当高，但模型复杂、参数物理意义不明确，参数确 定也无一定理论和实验基础。在无法用大量实验数据支持的条件下，利用材料的结构和物理理论．来获 各种物理问题及软件中的广泛应用，我们重点研究如何确定相关参数。 2 steinberg书uinan模型及其参数 2．1 steiⅡber窑．GuiⅡan模型 stejnbe唱．Guin∞提出了在材料熔化前，剪切模量和流变应力的表达式分别为： p=po【l+bpVV3一h盯一300)】 (1c) Y=[Y0(1+Di9)“】【l+bpV7一h(T一300)】 (1d) i9为等效塑性应变，V为相对体积，它与比容v的关系，为v=p。v，v为比容。 且：Y0(1+p云-)“≤Y～：TTm时，Y=肛=o，T_=T-oViexp(2丫o(1一V))。 这里．Yo．Y赫，p，n。№，b，h，1‰和po为给定介质常数。 d￡：一de；-d￡：是由运动方 等效塑性应变i’=Jd￡9，其中，d￡-=胄de：de：，(i，j=x‘e)：并且·d￡：2 程给出的总应变增量：d￡：是由Hook’s定律获得的应变增量。 2．2 steinberg．Guinan模型参数的物理意义 steinbe曙．Guinan模型参数中。剪切模量对压力和温度的导数的确定最为关键。 M．w second-0rder Guinan和D．J．steinbe曜通过对高阶柔度矩阵和刚性矩阵的讨论，由二阶弹性常数(the eIastic constants．SOEC)，及三阶弹性常数(theThird-orderelastic 素的各向同性多晶材料的剪切模量对压力和温度的导数。即steinbe嚷一GuinⅫ本构模型中与温度效应和 压力效应相关的参数。有效的多晶材料弹塑性常数必须由单晶数据采用适当的平均方法得到。 对于非立方晶体，要求知道所有的三阶弹性常数。这是很难满足的条件，好在我们要解决的材料大 第八届全国冲击动力学学讨论会术 多是立方晶体。但上述计算方法的实现还是需要一个过程．目前这一计算我们还实现不了。尽管如此 T=3000K)。以满足目前的总体设计需要。 Gm ¨u(吼+{“。)(I+6T：) ¨ ¨ 2“舻6 2(声+；“)拶 6 ”E一 舻；鲁暴专忙；南!掣一÷飞=2(晚+2¨。)瓯’一志…r=鲁2(凰+2p．．)‘ 岛 竺竺堡兰，B％矾一， …， 4(I+吉’ 其中丸是冲击速度～粒子速度关系的斜率。剪切模量与体模量、杨氏模量的关系为：A：』堕，E为 ‘ 9占一E Y0ng’s模量。 (觌=忐障)，