压剪载荷作用下TB6钛合金的动态力学性能

邹学韬 张晓晴 姚小虎

引用本文:
Citation:

压剪载荷作用下TB6钛合金的动态力学性能

    作者简介: 邹学韬(1993-),男,硕士研究生,主要从事冲击动力学研究. E-mail: zouxt-scut@outlook.com;
    通讯作者: 姚小虎, yaoxh@scut.edu.cn
  • 中图分类号: O347.3

Dynamic Behavior of TB6 Titanium Alloy under Shear-Compression Loading

    Corresponding author: YAO Xiaohu, yaoxh@scut.edu.cn
  • CLC number: O347.3

  • 摘要: 钛合金以其轻质高强的优异力学性能被广泛应用于航空航天领域。使用Instron万能材料试验机和分离式霍普金森压杆,对TB6钛合金进行准静态和动态力学性能实验,得到了压缩、拉伸和压剪载荷作用下TB6钛合金的准静态和动态应力-应变曲线,构建了单轴压缩和纯剪切两种应力状态下的Johnson-Cook动态本构模型。结果表明,TB6钛合金的屈服应力表现出明显的拉压不对称性、应变率强化和热软化效应。使用拉压不对称因子,修正了von Mises屈服准则,修正的屈服准则可很好地预测TB6钛合金的准静态和动态屈服行为。
  • 图 1  压剪试件和帽型试件的应力状态

    Figure 1.  Stress of shear-compression and cap-type specimens

    图 2  试件应变云图

    Figure 2.  Nephograms of strain distribution in specimens

    图 3  准静态正应力-应变曲线

    Figure 3.  Quasi-static normal stress-strain curve

    图 4  准静态切应力-应变曲线

    Figure 4.  Quasi-static shear stress-strain curve

    图 5  SHPB实验典型波形

    Figure 5.  Typical waves of SHPB experiments

    图 6  SHPB动态压缩应力-应变曲线

    Figure 6.  Stress-strain curves of TB6 titanium alloy under SHPB dynamic compression

    图 7  绝热温升曲线

    Figure 7.  Time history of temperature rise during adiabatic process

    图 8  不同应变率下J-C模型预测与实验结果的比较

    Figure 8.  Stress-strain cures at different strain rates from experiments and J-C model fitting

    图 9  不同应变率下的初始屈服面

    Figure 9.  Initial yield surface at different strain rates

    图 10  应变率3500 s-1下的后继屈服面

    Figure 10.  Succeeding yield surface at strain rate of 3500 s-1

    表 1  J-C模型参数拟合结果

    Table 1.  Fitting results of J-C model parameters

    ConditionA/MPaB/MPanm
    Uniaxial compression939326.60.240.661 18
    Pure shear495 92.50.210.655 48
    下载: 导出CSV
  • [1] 欧阳德来. TB6和TA15钛合金β锻组织演变及动态再结晶行为研究 [D]. 南京: 南京航空航天大学, 2011: 1-2
    [2] 吴琳, 王克鲁, 鲁世强. 基于逐步回归法的TB6钛合金本构关系研究 [J]. 热加工工艺, 2010, 39(8): 29–35 doi: 10.3969/j.issn.1001-3814.2010.08.010
    WU L, WANG K L, LU S Q. Study on constitutive relationship of TB6 alloy based on stepwise regression method [J]. Hot Working Technology, 2010, 39(8): 29–35 doi: 10.3969/j.issn.1001-3814.2010.08.010
    [3] 雷力明, 黄旭, 黄利军. 铸态TB6钛合金热变形行为及本构关系 [J]. 中国有色金属学报, 2010, 20(Suppl 1): 377–380
    LEI L M, HUANG X, HUANG L J. Hot deformation behavior and constitutive relationship of as-cast TB6 alloy [J]. The Chinese Journal of Nonferrous Metals, 2010, 20(Suppl 1): 377–380
    [4] 段园培, 黄仲佳, 余小鲁. 基于摩擦修正的TB6合金流变应力行为研究及本构模型建立 [J]. 稀有金属, 2014, 38(2): 202–209
    DUAN Y P, HUANG Z J, YU X L. Flow stress behavior and constitutive model of as-cast TB6 titanium alloy based on friction correction [J]. Chinese Journal of Rare Metals, 2014, 38(2): 202–209
    [5] WU Y, LIU J, WANG H, et al. Effect of stress ratio on very high cycle fatigue properties of Ti-10V-2Fe-3Al alloy with duplex microstructure [J]. Journal of Materials Science & Technology, 2018, 34(7): 1189–1195.
    [6] LI J, LI F G, MA X K, et al. Micromechanical study of the forged Ti-1023 titanium alloy by micro-indentation [C]//Key Engineering Materials. Trans Tech Publications, 2018, 765: 160-165.
    [7] RITTEL D, LEE S, RAVICHANDRAN G. A shear-compression specimen for large strain testing [J]. Experimental Mechanics, 2002, 42(1): 58–64. doi: 10.1007/BF02411052
    [8] 周刚毅, 董新龙, 付应乾. 不同加载状态下TA2钛合金绝热剪切破坏响应特性 [J]. 力学学报, 2016, 48(6): 1353–1361
    ZHOU G Y, DONG X L, FU Y Q. An experimental study on adiabatic shear behavior of TA2 titanium alloy subject to different loading condition [J]. Chinese Journal of Theoretical & Applied Mechanics, 2016, 48(6): 1353–1361
    [9] JIN T, ZHOU Z, SHU X, et al. Effects of strain rate on PMMA failure behavior [J]. Applied Physics A, 2016, 122(1): 7. doi: 10.1007/s00339-015-9526-0
    [10] 宋立, 胡时胜. SHPB数据处理中的二波法与三波法 [J]. 爆炸与冲击, 2005, 25(4): 368–373 doi: 10.3321/j.issn:1001-1455.2005.04.014
    SONG L, HU S S. Two-wave and three-wave method in SHPB data processing [J]. Explosion and Shock Waves, 2005, 25(4): 368–373 doi: 10.3321/j.issn:1001-1455.2005.04.014
    [11] JOHNSON G R, COOK W H. A constitutive model and data for metals subjected to large strains, high strain rates and high temperatures [C]//Proceedings of the 7th International Symposium on Ballistic, 1983: 541-547.
    [12] 徐天平, 王礼立, 卢维娴. 高应变率下钛合金Ti-6AI-4V的热-粘塑性特性和绝热剪切变形 [J]. 爆炸与冲击, 1987, 7(1): 1–8
    XU T P, WANG L L, LU W X. The thermo-visco plasticity and adiabatic shear deformation for a titanium alloy Ti-6Al-4V under high strain rates [J]. Explosion and Shock Waves, 1987, 7(1): 1–8
    [13] CAZACU O, PLUNKETT B, BARLAT F. Orthotropic yield criterion for hexagonal closed packed metals [J]. International Journal of Plasticity, 2006, 22(7): 1171–1194. doi: 10.1016/j.ijplas.2005.06.001
    [14] KHAN A S, YU S, LIU H. Deformation induced anisotropic responses of Ti-6Al-4V alloy Part II: a strain rate and temperature dependent anisotropic yield criterion [J]. International Journal of Plasticity, 2012, 38(4): 14–26.
  • [1] 尚兵王彤彤 . 一种立式分离式霍普金森压杆实验装置研制. 高压物理学报, 2018, 32(4): 043201-1-043201-5. doi: 10.11858/gywlxb.20170658
    [2] 王立闻庞宝君陈勇张凯 . 高温处理后活性粉末混凝土动力学行为及本构模型研究. 高压物理学报, 2012, 26(4): 361-368. doi: 10.11858/gywlxb.2012.04.001
    [3] 贾宇刘彦梁晓璐郑腾 . 高强度合金钢30CrMnMoRE/30CrMnSi的动态力学性能. 高压物理学报, 2018, 32(4): 044102-1-044102-7. doi: 10.11858/gywlxb.20170659
    [4] 李俊玲卢芳云傅华赵玉刚谭多望 . 某PBX炸药的动态力学性能研究. 高压物理学报, 2011, 25(2): 159-164 . doi: 10.11858/gywlxb.2011.02.012
    [5] 姚小虎任会兰林荣张晓晴 . 聚合物泡沫材料动态力学性能及其能量吸收研究. 高压物理学报, 2012, 26(5): 531-536. doi: 10.11858/gywlxb.2012.05.008
    [6] 王政倪玉山曹菊珍张文 . 冲击载荷下混凝土本构模型构建研究. 高压物理学报, 2006, 20(4): 337-344 . doi: 10.11858/gywlxb.2006.04.001
    [7] 许金余刘健范飞林任韦波席阳阳杨坤 . 高温SHPB冲击实验技术及其应用. 高压物理学报, 2013, 27(1): 57-62. doi: 10.11858/gywlxb.2013.01.008
    [8] 陈青山苗应刚郭亚洲李玉龙 . 比较93钨合金材料的3种本构模型. 高压物理学报, 2017, 31(6): 753-760. doi: 10.11858/gywlxb.2017.06.010
    [9] 李定远朱志武卢也森 . 冲击加载下42CrMo钢的动态力学性能及其本构关系. 高压物理学报, 2017, 31(6): 761-768. doi: 10.11858/gywlxb.2017.06.011
    [10] 李磊张先锋吴雪高飞刘闯 . 不同硬度30CrMnSiNi2A钢的动态本构与损伤参数. 高压物理学报, 2017, 31(3): 239-248. doi: 10.11858/gywlxb.2017.03.005
    [11] 刘海峰韩莉 . 二维骨料随机分布混凝土的动态力学性能数值模拟. 高压物理学报, 2016, 30(3): 191-199. doi: 10.11858/gywlxb.2016.03.003
    [12] 肖大武李英雷胡时胜 . 高温高应变率下纯锆的本构模型研究. 高压物理学报, 2009, 23(1): 46-50 . doi: 10.11858/gywlxb.2009.01.008
    [13] 李雪艳李志斌张舵 . 考虑温度效应的泡沫铝准静态压缩本构模型. 高压物理学报, 2018, 32(4): 044103-1-044103-6. doi: 10.11858/gywlxb.20170642
    [14] 陈丁丁卢芳云林玉亮蒋邦海 . 某含铝PBX压缩性能的应变率与温度效应. 高压物理学报, 2013, 27(3): 361-366. doi: 10.11858/gywlxb.2013.03.007
    [15] 傅华彭金华李俊玲唐孝容谭多望 . 一种适合于Hopkinson杆的实验测试方法探索. 高压物理学报, 2014, 28(4): 423-428. doi: 10.11858/gywlxb.2014.04.006
    [16] 雷忠琦姚小虎龙舒畅常建虎王海波 . 聚苯乙烯泡沫材料的动态压缩特性. 高压物理学报, 2019, 33(2): 024202-1-024202-10. doi: 10.11858/gywlxb.20180655
    [17] 刘海峰吴萍 . 具有不同界面相厚度的混凝土动态特性. 高压物理学报, 2017, 31(3): 249-260. doi: 10.11858/gywlxb.2017.03.006
    [18] 林玉亮卢芳云卢力 . 石英压电晶体在霍普金森压杆实验中的应用. 高压物理学报, 2005, 19(4): 299-304 . doi: 10.11858/gywlxb.2005.04.003
    [19] 赵习金卢芳云王悟李英华林玉亮 . 入射波整形技术的实验和理论研究. 高压物理学报, 2004, 18(3): 231-236 . doi: 10.11858/gywlxb.2004.03.007
    [20] 林玉亮卢芳云王晓燕赵习金 . 低密度聚氨酯泡沫压缩行为实验研究. 高压物理学报, 2006, 20(1): 88-92 . doi: 10.11858/gywlxb.2006.01.017
  • 加载中
图(10)表(1)
计量
  • 文章访问数:  239
  • 阅读全文浏览量:  243
  • PDF下载量:  7
出版历程
  • 收稿日期:  2019-01-17
  • 录用日期:  2019-01-28
  • 网络出版日期:  2019-03-21
  • 刊出日期:  2019-04-01

压剪载荷作用下TB6钛合金的动态力学性能

    作者简介:邹学韬(1993-),男,硕士研究生,主要从事冲击动力学研究. E-mail: zouxt-scut@outlook.com
    通讯作者: 姚小虎, yaoxh@scut.edu.cn
  • 华南理工大学土木与交通学院,广东 广州 510641

摘要: 钛合金以其轻质高强的优异力学性能被广泛应用于航空航天领域。使用Instron万能材料试验机和分离式霍普金森压杆,对TB6钛合金进行准静态和动态力学性能实验,得到了压缩、拉伸和压剪载荷作用下TB6钛合金的准静态和动态应力-应变曲线,构建了单轴压缩和纯剪切两种应力状态下的Johnson-Cook动态本构模型。结果表明,TB6钛合金的屈服应力表现出明显的拉压不对称性、应变率强化和热软化效应。使用拉压不对称因子,修正了von Mises屈服准则,修正的屈服准则可很好地预测TB6钛合金的准静态和动态屈服行为。

English Abstract

  • TB6钛合金是我国在美国Ti-1023钛合金(Ti-10V-2Fe-3Al)基础上自主开发的一种高强韧近β型合金,既有诸多亚稳β钛合金的优点,又有(α+β)两相钛合金的固有特性,具有比强度高、断裂韧性好、淬透截面大、各向异性小、锻造温度低、抗应力腐蚀能力强等优点[1],广泛应用于航空航天飞行器结构中。航空航天飞行器工作时处于高速状态,受力复杂,且面临鸟撞等高速冲击的威胁,因此有必要研究其在压剪等复杂应力状态下的动态力学行为。

    近年来,众多学者对TB6钛合金的本构关系开展了研究,加深了对其力学行为的理解。吴琳等[2]基于参数逐步回归法建立了TB6钛合金在不同温度范围的Arrhenius幂函数型本构关系。雷力明等[3]通过研究铸态TB6钛合金的高温热变形行为,建立了考虑应变的流变应力影响的Arrhenius双曲正弦型本构关系。段园培等[4]基于摩擦修正了TB6钛合金流变应力曲线,建立了Arrhenius双曲正弦型热流变本构关系。Wu等[5]通过研究TB6钛合金在超高周疲劳状态下的S-N曲线和微观裂纹,揭示了应变比和应力水平对S-N曲线和裂纹萌生机制的影响。Li等[6]通过显微压痕试验,研究了TB6钛合金的弹性模量、显微硬度等力学参数,并采用幂函数型应力-应变方程建立了微压痕塑性本构关系。当前关于TB6钛合金的研究主要集中在高温环境下低应变率压缩加载的热流变应力-应变关系和疲劳性能,缺少高应变率以及压剪等复杂应力状态下的本构关系和屈服准则研究。

    本工作通过准静态和霍普金森压杆(Split Hopkinson Pressure Bar, SHPB)实验,研究TB6钛合金在压剪复杂应力状态下的动态力学性能,通过分析其力学性能的应变率相关性,建立高应变率范围的Johnson-Cook本构模型,最后通过引入拉压不对称因子修正von Mises屈服准则,得到一种适用于描述TB6钛合金高应变率屈服行为的修正屈服准则。

    • 为探讨不同应力状态下TB6钛合金的准静态和动态力学性能,分别对圆柱压缩试样、压剪试样(SCS)[7](开槽角度15°、30°和45°)和帽型纯剪切试样开展准静态和SHPB动态实验。试样材料为冷轧退火TB6钛合金板材,使用线切割和CNC雕刻机加工成圆柱压缩试样、压剪试样和帽型纯剪切试样,如图1所示。在本实验中,压剪试件的尺寸为:截面直径d为4 mm,H为1.25 mm,斜槽厚度T为2 mm;帽型纯剪切试样的尺寸为:w1=4.00 mm,w2=4.10 mm,h=2.00 mm,试样的厚度δ=5.00 mm。

      图  1  压剪试件和帽型试件的应力状态

      Figure 1.  Stress of shear-compression and cap-type specimens

    • Rittel等[7]对圆柱试样加工不同角度的斜槽,实现金属材料大应变范围的压剪应力加载。Jin等[9]分析了压剪试样的应力状态,如图1(a)所示。将内力沿斜槽分解成正应力和剪应力,得到正应力和剪应力的表达式

      $\sigma = {F_{{\rm{UTM}}}}{\cos ^2}\alpha /({{dT}}) $

      $\tau = {F_{{\rm{UTM}}}}\cos \alpha \sin \alpha/({{dT}}) $

      式中:FUTM为压缩载荷,στ分别为正应力和剪应力,d为截面直径,T为斜槽厚度,α为斜槽角度。与应力相同,将位移沿斜槽分解,得到正应变和剪应变的表达式

      ${\varepsilon _{\rm{n}}} = \Delta D {\cos ^2}\alpha /{H}$

      $\gamma = \Delta D \cos \alpha \sin \alpha /{H}$

      式中:$\varepsilon_{\rm n}$$ \gamma$分别为正应变和剪应变,ΔD为试样压缩位移。

      帽型纯剪切试样如图1(b)所示。受帽型纯剪切试样在加载过程中的几何特征和结构变形的影响,剪切应力区的应力状态并不是理想的纯剪切状态,剪切区应力状态如图1(b)右图所示。周刚毅等[8]研究发现,实验设计时帽型试样满足平面应变条件,当受迫剪切区受剪方向接近最大剪切方向时,可以认为该帽型试样在加载过程中近似处于平面纯剪切状态。

      剪切区剪应力和剪应变为

      $\tau = {F_{{\rm{UTM}}}}/{A} = {F_{{\rm{UTM}}}}/({{2\delta h}})$

      $\gamma = \Delta \theta \approx \tan \theta = \Delta D/\Delta w$

      式中:A为剪切区的剪切面积,δ为帽型试样厚度,h为设计剪切区高度,Δw为剪切变形区宽度。

    • SHPB通过调整入射波形实现材料的恒应变率动态加载,广泛应用于金属、陶瓷、PMMA等材料的高应变率动态力学性能研究。SHPB主要由空气枪、子弹、入射杆、透射杆、试样、光电测速器、超动态应变仪、数字化存储器、吸能器等部分组成。

      对于金属等塑性明显的韧性材料,根据一维应力波假定,可利用二波法[10]得到试样两端的动态压缩载荷FUTM(t)和压缩位移ΔD(t)

      ${F_{{\rm{UTM}}}}(t) = {E_{{\rm{bar}}}}{A_{{\rm{bar}}}}{\varepsilon _{\rm{t}}}(t)$

      $\Delta D(t) = - 2{C_0}\int_0^{{t_0}} {{\varepsilon _{\rm{r}}}} (t){\rm{d}}t$

      式中:$\varepsilon_{\rm r}$$ \varepsilon_{\rm t}$分别为入射杆和透射杆上测得的反射应变和透射应变,EbarAbar分别为杆的弹性模量和截面面积,C0为杆中应力波波速。

      在压剪和纯剪切等复杂应力加载过程中,应变并非单向应变,因此用等效应变量化不同加载状态的等效应变率。假设材料不可压缩,得到Mise等效应变和等效应变率为

      ${\varepsilon _{{\rm{eq}}}} = \sqrt {\frac{2}{3}{\varepsilon _{ij}}{\varepsilon _{ij}}} = \sqrt {{\varepsilon _{\rm n}^2} + \frac{1}{3}{\gamma ^2}} $

      ${\dot \varepsilon _{{\rm{eq}}}} = \frac{{3{\varepsilon _{\rm{n}}}{{\dot \varepsilon }_{\rm{n}}} + \gamma \dot \gamma }}{{\sqrt {9\varepsilon _{\rm{n}}^2 + 3{\gamma ^2}} }}$

      式中:$\varepsilon_{\rm n}$为正应变,$ \gamma $为剪应变,${{\dot \varepsilon }_{\rm{n}}}$为正应变率,$\dot \gamma$为剪应变率。

    • 所有准静态压缩实验都在Instron万能材料试验机上进行,应变率为0.001 s–1。为了获得压剪试样压剪应力区和帽型纯剪切试样剪切区的应变,采用工业相机拍摄准静态压缩实验过程中试样的变形过程,利用GOM Correlate软件获得试样的应变场分布。

      图2为利用二维图像相关法(DIC-2D)得到的压剪试样压剪应力区和帽型纯剪切试样剪应力区的准静态应变场。可见,压剪试样斜槽端部出现应变集中,而压剪应力区、纯剪切应力区的应变近似均匀,可以对试样进行SHPB动态压缩实验。选取相同平面上的3个点,计算平均应变。对应变-时间曲线进行线性插值,并与载荷-时间曲线对应,得到单轴压缩、压剪、纯剪切和拉伸状态的应力-应变关系。

      图  2  试件应变云图

      Figure 2.  Nephograms of strain distribution in specimens

      准静态正应力-应变曲线如图3所示:准静态压缩下,TB6钛合金的单轴压缩屈服应力为939 MPa,进入塑性阶段后应力强化,强度达到1152 MPa;单轴拉伸屈服应力为773 MPa,塑性阶段强度为893 MPa;准静态单轴压缩和拉伸明显出现不对称性;随着斜槽角度的增大,正应力减小。

      图  3  准静态正应力-应变曲线

      Figure 3.  Quasi-static normal stress-strain curve

      准静态切应力-应变曲线如图4所示:纯剪切加载时剪应力最大,屈服应力为495 MPa;塑性阶段同样出现应力强化现象,强度达到567 MPa;斜槽角度越大,加载过程中剪切作用越强,剪应力越大。

      图  4  准静态切应力-应变曲线

      Figure 4.  Quasi-static shear stress-strain curve

    • SHPB动态压缩实验中,子弹、入射杆和透射杆的长度分别为200、1200和1200 mm,杆和子弹的直径均为14 mm。采用TB6圆形薄片作为入射波整形片,实现试样加载过程中的变形均匀和应力平衡。图5显示实验得到的典型波形。可见,单轴压缩反射波出现材料热软化效应所特有的波形,而压剪加载的反射波未出现明显的热软化现象。

      图  5  SHPB实验典型波形

      Figure 5.  Typical waves of SHPB experiments

      图6为动态单轴压缩、压剪和纯剪切3种应力状态下SHPB动态压缩实验得到的应力-应变曲线,其中图6(a)为单轴压缩应力状态下的正应力-应变曲线,图6(b)为纯剪切应力状态下的剪应力-应变曲线,图6(c)图6(e)图6(g)为不同斜槽角度压剪应力状态下的正应力-应变曲线,图6(d)图6(f)图6(h)为不同斜槽角度剪切试样的剪应力-应变曲线。3种应力状态的正应力具有明显的应变率相关性,屈服强度随应变率的升高而增加。由图6(a)可知,高应变率单轴压缩加载下,材料进入塑性阶段后,应力首先开始强化,随着应变增加,绝热温升引起的热软化效应使材料强度下降,并保持应力平台。对比不同斜槽角度的正应力-应变曲线和剪应力-应变曲线可知,随着斜槽角度的增大,正应力分量的屈服强度减小,而剪应力分量的屈服强度增大,说明通过改变斜槽角度可以改变试样的应力状态。

      图  6  SHPB动态压缩应力-应变曲线

      Figure 6.  Stress-strain curves of TB6 titanium alloy under SHPB dynamic compression

    • 工程上常用热-黏塑性本构模型描述金属材料的动态力学响应,常见的热-黏塑性本构模型有Zerrilli-Armstrong模型、Johnson-Cook(J-C)模型[11]和Steinberg-Guinan模型。考虑到J-C模型能够较好地描述材料的应变硬化、应变率强化和热软化效应,本研究采用J-C模型对TB6钛合金单轴压缩和纯剪切动态应力-应变曲线进行拟合,其表达式为

      $\sigma = (A + B{\varepsilon ^n})(1 + C\ln {\dot \varepsilon ^*})(1 - {T^{*m}})$

      式中:σ为塑性流变应力;A为参考应变率下的屈服应力;Bn为拟合参数;$\varepsilon$为塑性应变;C为应变率敏感系数;${\dot \varepsilon ^{\rm{*}}} = \dot \varepsilon /{\dot \varepsilon _0}$$\dot \varepsilon$${\dot \varepsilon _0}$分别为等效应变率和参考应变率;T*=(TT0)/(TmT0),T0TTm分别为参考温度、变形温度和材料的熔点,T0Tm分别取293 K和1800 K。

      单轴压缩J-C动态本构拟合时,以准静态应变率$\dot \varepsilon$=0.001 s-1为参考应变率,A=939 MPa,拟合得到B=326.6 MPa,n=0.24,C=–0.016 33+8.56×10–6$\dot \varepsilon$。采用同样方法拟合纯剪切J-C本构,$\dot \varepsilon$=0.001 s–1A=495 MPa,拟合得到B=92.5 MPa,n=0.21,C=0.0227–1.375×10–5$\dot \varepsilon$+2.995×10–9${\dot \varepsilon ^2}$

      金属材料在高应变率压缩变形过程中伴随着显著的热-力耦合效应,导致应力的热软化,在较小的时间和空间尺度上该过程为绝热过程。准静态压缩曲线可看作等温曲线,高应变率下的动态压缩曲线可看作绝热曲线,高应变率下的应变硬化率实际上是应变硬化率效应和绝热温升引起的热软化效应的综合反映[12]。理论上认为,金属材料的绝热温升由塑性变形过程中的塑性功转化而成,温升ΔT

      $\Delta T = \beta {W_{\rm{p}}}/({{\rho {c_p}}})$

      ${W_{\rm{p}}} = \int_0^{{\varepsilon _{\rm{p}}}} {{\sigma _{\rm{p}}}{\rm{d}}{\varepsilon _{\rm{p}}}} $

      式中:β为热转化效率,一般取0.9;Wp为塑性功;$ \rho$为材料密度,取4.47 g/cm3cp为定压比热容,取0.576 J/(g∙K);σpεp分别为塑性应力和塑性应变。

      图7所示,动态加载过程中绝热温升与塑性应变近似呈线性关系,且与应变率相关。以往人们在研究材料的温升软化时只考虑温升与应变率的关系,使得拟合J-C模型参数m时无法准确地表示屈服应力。绝热温升曲线的拟合公式为

      图  7  绝热温升曲线

      Figure 7.  Time history of temperature rise during adiabatic process

      ${\text{单轴压缩}}\quad\quad\quad\quad\quad\quad\Delta {T_{{\rm{uc}}}} = f(\dot \varepsilon )\varepsilon = (809.95 - 0.49\dot \varepsilon + 1.823 \times {10^{ - 4}}{\dot \varepsilon ^2} - 2.004 \times {10^{ - 8}}{\dot \varepsilon ^3})\varepsilon \quad\quad\quad\quad\quad\quad $

      ${\text{纯剪切}}\quad\quad\quad\quad\quad\quad\;\Delta {T_{{\rm{ps}}}} = f(\dot \varepsilon )\varepsilon = (112.31 + 0.0077\dot \varepsilon - 6.79 \times {10^{ - 7}}{\dot \varepsilon ^2} + 1.326 \times {10^{ - 9}}{\dot \varepsilon ^3})\varepsilon \quad\quad\quad\quad\quad\quad$

      高应变率加载过程中变形温度T等于绝热温升ΔT,将(14)式和(15)式代入(11)式,拟合参数m,取平均值。J-C模型拟合参数列于表1

      ConditionA/MPaB/MPanm
      Uniaxial compression939326.60.240.661 18
      Pure shear495 92.50.210.655 48

      表 1  J-C模型参数拟合结果

      Table 1.  Fitting results of J-C model parameters

      图8为不同应变率加载下单轴压缩和纯剪切J-C模型拟合结果与SHPB动态压缩实验曲线的对比。拟合曲线很好地表现了TB6的初始屈服应力和应力的热软化现象,且相对偏差在5%以内。可以认为,该本构模型有效地预测了TB6钛合金在室温下的单轴压缩和纯剪切高应变率加载。拟合得到的完整J-C本构关系为

      图  8  不同应变率下J-C模型预测与实验结果的比较

      Figure 8.  Stress-strain cures at different strain rates from experiments and J-C model fitting

      ${\sigma _{\rm{s}}} = (939 + 326.6{\varepsilon ^{0.24}})[1 + ( - 0.016\;33 + 8.56 \times {10^{ - 6}}\dot \varepsilon )\ln {\dot \varepsilon ^*}]\left[ {1 - {{\left( {\frac{{\Delta {T_{{\rm{uc}}}}}}{{1780}}} \right)}^{0.661\;18}}} \right]$

      ${\tau _{\rm{s}}} = (495 + 92.5{\varepsilon ^{0.21}})[1 + (0.0227 - 1.375 \times {10^{ - 5}}\dot \varepsilon + 2.995 \times {10^{ - 9}}{\dot \varepsilon ^2})\ln {\dot \varepsilon ^*}]\left[ {1 - {{\left( {\frac{{\Delta {T_{{\rm{ps}}}}}}{{1780}}} \right)}^{0.655\;48}}} \right]$

    • Tresca最大剪应力准则和von Mises畸变能准则是目前应用最广泛的各向同性屈服准则。然而,大量实验结果表明,钛合金表现出明显的拉压不对称性,传统的Tresca和von Mises屈服准则无法准确地描述钛合金在复杂应力状态下的屈服行为。材料在不同应力状态下的力学响应与偏应力第三不变量J3有关。Cazacu等[13]基于J3偏应力不变量提出了拉压不对称因子,由此描述材料力学的不对称行为。Khan等[14]针对Ti-6Al-4V钛合金提出一种新型的拉压不对称因子G(σ)。本研究基于拉压不对称因子G(σ)和von Mises屈服准则,构建一种新的修正von Mises屈服准则。屈服函数为

      $\phi = f(\sigma )G(\sigma ) = 1$

      $f(\sigma ) = 3{J_2}/{\sigma _{\rm{s}}^2}$

      $G(\sigma ) = \exp \left[ { - c(\xi + 1)} \right] = 1$

      式中:$f(\sigma ) = 1$为von Mises屈服函数;σs为单轴压缩屈服应力;$\xi = {\rm{cos}}\left( {3\theta } \right) =$$ \displaystyle\frac{{27}}{2} \cdot \displaystyle\frac{{{J_3}}}{{{(3{J_2})^{3/2}}}}$为Lode参数,其中θ为Lode角;c为拉压不对称系数。

      纯剪切加载时,$\xi =0$J2=${\tau _{\rm{s}}^2}$,于是(18)式简化为

      ${\rm{exp(}} - c{\rm{)}} = {\sigma _{\rm{s}}^2}/(3{\tau _{\rm{s}}^2})$

      式中:$\tau_{\rm s}$为纯剪切屈服剪应力。对(21)式等号两边取自然对数,得到拉压不对称系数c

      $c = - 2\ln \left(\frac{{{\sigma _{\rm{s}}}}}{{\sqrt 3 {\tau _{\rm{s}}}}}\right)$

      由此得到屈服函数

      $\phi = \frac{{3{J_2}}}{{{\sigma _{\rm{s}}^2}}}{\rm{exp}}\left[ {2\ln \left(\frac{{{\sigma _{\rm{s}}}}}{{\sqrt 3 {\tau _{\rm{s}}}}}\right)(\xi + 1)} \right] = 1$

      联合(18)式、(19)式和(23)式,可以得到完整的等向强化拉压不对称屈服准则。

      图9中的初始屈服面所示,准静态加载下屈服面能够准确地预测TB6钛合金的单轴压缩、纯剪切和单轴拉伸3种应力状态的屈服应力,并能很好地预测压剪应力状态的屈服应力;在高应变率范围内,屈服面仍能很好地预测材料的屈服应力。由图10可见,高应变率加载下TB6钛合金首先经历应力强化阶段,随后塑性流变应力未出现强化,很好地预测了材料单轴压缩、压剪和纯剪切3种应力状态下塑性阶段的应力流动。

      图  9  不同应变率下的初始屈服面

      Figure 9.  Initial yield surface at different strain rates

      图  10  应变率3500 s-1下的后继屈服面

      Figure 10.  Succeeding yield surface at strain rate of 3500 s-1

    • 利用Instron万能材料试验机和SHPB,对TB6钛合金进行了单轴、压剪和纯剪切3种应力状态的准静态和动态加载实验,得到了不同应力状态的准静态和动态应力-应变曲线,通过分析实验结果,得到如下结论:(1)TB6钛合金具有良好的塑性,准静态压缩下表现出明显的应力强化效应,高应变率单轴压缩下温升热软化效应显著,塑性应力在强化后迅速软化,而压剪加载和纯剪切加载下,材料没有表现出明显的热软化现象;(2)针对单轴加载和纯剪切加载,通过改进的J-C本构模型拟合,得到了与实验结果拟合程度较好的J-C参数,模型预测结果与实验结果的相对误差小于5%;(3)根据Khan等提出的拉压不对称因子G(σ)和von Mises屈服准则,提出一种新的修正von Mises屈服准则,该屈服准则能够很好地描述TB6钛合金的动态压剪力学行为和拉压不对称性。

参考文献 (14)

目录

    /

    返回文章
    返回