高温高压环境下金属材料的本构关系研究

摘要

高温高压条件下固体材料的本构关系是研究化学炸药爆轰、高能离子束或激光束以及高速撞击作用下,固体中状态量变化规律和运动特性必不可少的材料特性方程。为了精准描述和预测上述的状态量变化和运动学特征,本文根据陈大年等人[高压物理学报,19(2005)193]和彭建祥[博士论文,中国工程物理研究院研究生部,2006]的见解,着重分析了目前受到广泛应用的由Steinberg等人提出的[J.Appl.Phys.,51(1980)1498:J.Appl.Phys.,65(1989)1528]金属本构模型中存在的固有问题,并提出了一个新的金属本构函数表示式。获得的主要结果有:
　　 (1)在全面分析了Burakovsky等人[Phys.Rev.B,71(2005)184118]、Nadal等人[J.Appl.Phys.,93(2003)2472]和俞宇颖等人[Chin.Phys.ett.,17(2008)64]针对Steinberg等人的金属材料本构模型提出的若干质疑及修改意见的基础上,提出了一个新的适用于密堆集晶体金属的本构模型,具体的函数表示式是()式中,G和Y为剪切模量和屈服强度;β和n为硬化系数,εp为塑性应变;T0为参考温度,通常取零温或室温。式中的G(P,T0),((e)G/(e)T)P,β,n和εP等值均可通过实验和理论计算方法求得。这个模型具有如下优点:
　　 --与Steinberg等人、Nadal等人和俞宇颖等人的意见相比较,本模型中G随P的变化没有采用在一阶导数项截断的近似表示,而是采用了G(P)的全函数表示形式;
　　 --与Steinberg等人、Nadal等人、Burakovsky等人的意见相比较,本模型中G随T的变化虽然仍采用了一阶导数项的近似表达,但抛弃了把((e)G/(e)T)P用((e)G/(e)T)0代替的做法,而是在向士凯等人[Phy.Rev.B70(2004)174102]第一原理计算结果启示下采用了把((e)G/(e)T)P,表示为以P为参数的函数式;
　　 (2)以LY12铝合金为例,本文给出的本构函数式是()()
　　 这组本构函数式的有效性得到了以下检验结果的支持:
　　 --上列G(P,T)本构函数式是基于Sin'ko等人[J.Phys.:Condens.Matter,4(2002)6989]和向士凯等人的第一原理计算结果导出的,因而具有坚实物理基础。此外,它的有效性还得到了Hügoniot压缩线上声速测量结果和等熵压缩线上剪切模量计算结果(通过有限形变理论的Birch-Murnaghan方程)的支持;
　　 --上列Y(P,T)本构模型的构建尽管在其物理构思上沿袭了Steinberg等人的做法,但由于本文G(P,T)本构函数式的不同,故在Y(P,T)本构函数式的具体形式上是与Steinberg等人的有所不同的。它的有效性得到了本文冲击加载实验结果(采用了Asay和Chhabildas基于各向同性硬化模型的数据处理方法)的支持。特别值得提出的是:实验结果再一次支持了G(P,T)/Y(P,T)≈常数近似关系在高压区的有效性(对于LY12铝合金,这个高压区的阀值是Pc≥30GPa);但若关注到Pc的低压区,则G(P,T)/Y*(P,T)≈常数近似关系(式中Y*=Y/[1+β(εp+εp1)]n)也是成立的。其实,这个结果也正是Steinberg等人构建他们本构模型的初衷。
　　 --为了综合检验本文G(P,T)和Y(P,T)本构函数式在描述冲击加载下LY12铝合金状态量变化及其运动特性的有效性,文中利用冲击/卸载、冲击/冲击/卸载、冲击/再冲击以及冲击/卸载/冲击等四类复合加载过程的样品/窗口界面速度剖面实验结果进行了检验,发现模型预测值与实测值的符合程度是令人满意的,从而全面认定了本文本构模型的有效性。
　　 (3)在发展实验技术方面,本文采用Asay等人[SAND-85-2009]建议的组合靶实验方案,解决了前人未能给出～20GPa以上压力区冲击/再冲击实验数据的难题(本文中给出了高达～70GPa预冲击压力点的实验数据),从而为满足今后精准测定更高压力区屈服强度的要求打通了一条新的技术途径。
　　 最后应该说明,本文提出的高温高压区金属本构模型与前人提出的相比,尽管它具有更为坚实的物理基础和更少的近似假设,但仅能看着是初步的结果,因为它的普适性还需要通过更多类别金属的检验。即使对LY12铝合金来说,也还需要更多的实验数据和理论计算结果的检验,因为,只有这样才能对本构函数式中的有关材料参数值进行精确化修正。为此目的,发展或优化现有的实验和理论计算方法也是值得重视的。在发展实验技术方面,看来进一步发展高温高压条件下的激光-超声(Laser-ultrasonic)和静态压缩下的屈服强度测量技术是必要的。

目录

文摘

英文文摘

第1章 绪论

1.1 研究背景

1.2 本构关系研究历史的回顾—类本构关系的提出

1.3 两类本构模型压力适用范围的研究

1.4 本文研究内容及章节安排

第2章 高温高压下金属剪切模量本构关系及铝的G(P,T)函数表达式

2.1 对SCG本构模型中G(P,T)函数式有效性的几个有代表性的质疑及改进意见

2.1.1 俞宇颖、谭华等人的意见

2.1.2 Nadal和Poac等人的意见

2.1.3 Burakovsky和Preston的意见

2.1.4 本节小结

2.2 本文对构制剪切模量本构函数关系式的考虑

2.2.1 问题的提出

2.2.2 一种新的剪切模量本构函数式

2.3 铝的剪切模量本构关系式

2.3.1 G(P,T0)函数式的确定

2.3.2((e)G/(e)T)p函数式的确定

2.4 GA1(P,T)本构关系式有效性的检验

2.4.1 用冲击压缩线上剪切模量随压力变化实测数据的检验

2.4.2 用等熵压缩线上剪切模量随压力变化的有限形变理论计算数据的检验

2.5 本章小结

第3章 高温高压下金属屈服强度本构关系及铝的Y(P,T)函数表达式

3.1 一维应变条件下固体材料形变过程及其弹塑性本构响应

3.2 高压区固体材料屈服强度测量的基本原理

3.2.1 从实测U1,(t)计算样品中声速C(t)

3.2.2 从U1(t)和C(t)数据计算D点之后样品卸载波中状态参量σ,ρ,和U

3.2.3 Asay和Chhabidas(AC)技术中的基本假设

3.3 本文对Al在高压区屈服强度的实验测量

3.3.1 背景情况说明

3.2.2 本文对冲击/再冲击实验的具体安排

3.3.3 冲击/再冲击实验测量结果

3.3.4 冲击/再冲击加载条件下LY12铝中的声速及状态参量σ,ρ,U,θ的计算

3.3.5 屈服强度的计算

3.4 Y/G=const关系的验证

3.5 小结

第4章 铝的G(P,T)和Y(P,T)函数式有效性的综合检验

4.1 引言

4.2 对本文提出的G(P,T)和Y(P,T)函数式预测能力的综合考核

4.2.1 与LY12铝合金冲击/再冲击实验结果的比较

4.2.2 与LY12铝合金冲击/再冲击/卸载实验结果的比较

4.2.3 与LY12铝合金加载/卸载实验结果的比较

4.2.4 与LY12铝合金加载/完全卸载/再加载实验结果的比较

4.3 本章小结

第5章 全文总结

参考文献

致谢

附录 攻读博士期间发表的论文