鸟撞冲击下TC4钛合金平板的变形和破坏

贾林 李从富 邹学韬 姚小虎

引用本文:
Citation:

鸟撞冲击下TC4钛合金平板的变形和破坏

    作者简介: 贾 林(1988-),男,硕士,工程师,主要从事发动机强度试验技术研究. E-mail:jialin20060598@126.com;
    通讯作者: 姚小虎, yaoxh@scut.edu.cn
  • 中图分类号: O347.3

Deformation and Destruction of TC4 Titanium Alloy Plate under the Bird Impact

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

  • 摘要: 为了分析航空用TC4钛合金板材在鸟撞冲击下的动态响应,通过3D-DIC动态变形场测试技术,研究了鸟撞冲击过程中TC4钛合金平板变形场,并基于显式有限元分析软件ABAQUS,建立了鸟撞数值计算模型。将考虑拉压不对称性的修正von Mises屈服准则引入TC4钛合金材料的Johnson-Cook动态本构模型和动态损伤模型中,采用光滑粒子流体动力学方法(SPH)建立了鸟体模型。数值计算和鸟撞实验的对比结果表明,计算应变与实验结果吻合很好,验证了TC4钛合金鸟撞冲击数值分析模型的合理性和可靠性。
  • 图 1  实验装置示意图

    Figure 1.  Schematic of experimental apparatus

    图 2  高速摄影机布局

    Figure 2.  High speed camera arrangement

    图 3  靶板尺寸

    Figure 3.  Size of target board

    图 4  149、167和180 m/s的鸟撞实验结果

    Figure 4.  Results of bird strike experiments at the speed of 149, 167 and 180 m/s

    图 5  TC4钛合金平板破坏照片

    Figure 5.  Failure of TC4 titanium alloy plate

    图 6  鸟弹的几何尺寸

    Figure 6.  Geometry of bird ball

    图 7  鸟撞数值计算模型

    Figure 7.  Numerical model of bird strike

    图 8  观测点布局

    Figure 8.  Distribution of observation points

    图 9  149 m/s鸟撞等效应力云图

    Figure 9.  Equivalent stress nephograms of bird impacting with the velocity of 149 m/s

    图 10  TC4钛合金平板鸟撞变形结果

    Figure 10.  Deformation of TC4 titanium alloy plate under the bird impact

    图 11  180 m/s工况下计算得到的TC4钛合金平板破坏过程

    Figure 11.  Failure process of TC4 titanium alloy plate calculated at 180 m/s

    图 12  位移-时间曲线的计算和实验结果对比

    Figure 12.  Comparison of calculated displacement-time curves with experimental results

    图 13  应变-时间曲线的计算与实验结果对比

    Figure 13.  Comparison of calculated strain-time curves with experimental results

    表 1  鸟体材料参数

    Table 1.  Material parameters of bird body

    Density/(kg·m−3Elastic modulus/GPaPoisson’s ratioYield stress/MPaFailure strainTangent modulus/MPa
    9280.0680.490.691.255
    下载: 导出CSV

    表 2  TC4钛合金材料参数[10-11]

    Table 2.  Parameters of TC4 titanium alloy[10-11]

    ρ/(g·cm−3)μTm/KA/MPaB/MPanCm
    4.4300.331 878106010900.8840.01171.1
    E/GPa${\dot \varepsilon _{{0}}}$/s−1D1D2D3D4D5
    1354 × 10−4−0.0900.2700.4800.0143.870
    下载: 导出CSV

    表 3  TC4钛合金平板鸟撞最大位移

    Table 3.  Maximum displacement of titanium alloy plate impacted by a bird

    Velocity/(m·s−1)MethodMaximum displacement/mm
    S1S2S3S4S5S6
    149Sim.655136645523
    Exp.686038635325
    167Sim.756048766527
    Exp.805947766324
    180Sim.1191016311810354
    Exp.11790531039550
    下载: 导出CSV
  • [1] 黄旭, 朱知寿, 王红红. 先进航空钛合金材料与应用 [M]. 北京: 国防工业出版社, 2012.
    HUANG X, ZHU Z S, WANG H H. Advanced aeronautical titanium alloys and applications [M]. Beijing: National Defense Industry Press, 2012.
    [2] GUIDA M, MARULO F, MEO M, et al. SPH-Lagrangian study of bird impact on leading edge wing [J]. Composite Structures, 2011, 93(3): 1060–1071. doi: 10.1016/j.compstruct.2010.10.001
    [3] DE VUYST T, VIGNJEVIC R, CAMPBELL J C. Coupling between meshless and finite element methods [J]. International Journal of Impact Engineering, 2005, 31(8): 1054–1064. doi: 10.1016/j.ijimpeng.2004.04.017
    [4] 张鼎逆, 上官倩芡, 刘富. 基于SPH方法的LY12-CZ铝合金平板鸟撞模型 [J]. 江苏大学学报: 自然科学版, 2016, 37(4): 418–422. doi: 10.3969/j.issn.1671-7775.2016.04.008
    ZHANG D N, SHANGGUAN Q Q, LIU F. Model of bird impact on LY12-CZ aluminum alloy plate based on SPH [J]. Journal of Jiangsu University (Natural Science Edition), 2016, 37(4): 418–422. doi: 10.3969/j.issn.1671-7775.2016.04.008
    [5] 刘军, 李玉龙, 刘元镛. 基于SPH方法的叶片鸟撞数值模拟研究 [J]. 振动与冲击, 2008, 27(9): 90–93. doi: 10.3969/j.issn.1000-3835.2008.09.022
    LIU J, LI Y L, LIU Y Y. Numerical simulation study of bird-impact on a blade using SPH method [J]. Journal of Vibration and Shock, 2008, 27(9): 90–93. doi: 10.3969/j.issn.1000-3835.2008.09.022
    [6] 刘富, 张嘉振, 童明波, 等. 2024-T3铝合金动力学实验及其平板鸟撞动态响应分析 [J]. 振动与冲击, 2014, 33(4): 113–118. doi: 10.3969/j.issn.1000-3835.2014.04.021
    LIU F, ZHANG J Z, TONG M B, et al. Dynamic tests and bird impact dynamic response analysis for a 2024-T3 aluminum alloy plate [J]. Journal of Vibration and Shock, 2014, 33(4): 113–118. doi: 10.3969/j.issn.1000-3835.2014.04.021
    [7] LIU J, LI Y L, GAO X S. Bird strike on a flat plate: experiments and numerical simulations [J]. International Journal of Impact Engineering, 2014, 70: 21–37. doi: 10.1016/j.ijimpeng.2014.03.006
    [8] 姚小虎, 韩强, 赵隆茂, 等. 飞机圆弧风挡的抗鸟撞击问题研究 [J]. 华南理工大学学报(自然科学版), 2007, 35(2): 6–12. doi: 10.3321/j.issn:1000-565X.2007.02.002
    YAO X H, HAN Q, ZHAO L M, et al. Investigation into arc windshield of aircraft subjected to bird impact [J]. Journal of South China University of Technology (Natural Science Edition), 2007, 35(2): 6–12. doi: 10.3321/j.issn:1000-565X.2007.02.002
    [9] 邹学韬, 张晓晴, 姚小虎. 压剪载荷作用下TB6钛合金的动态力学性能 [J]. 高压物理学报, 2019, 33(2): 024206. doi: 10.11858/gywlxb.20190713
    ZOU X T, ZHANG X Q, YAO X H. Dynamic behavior of TB6 titanium alloy under shear-compression loading [J]. Chinese Journal of High Pressure Physics, 2019, 33(2): 024206. doi: 10.11858/gywlxb.20190713
    [10] 惠旭龙, 牟让科, 白春玉, 等. TC4钛合金动态力学性能及本构模型研究 [J]. 振动与冲击, 2016, 35(22): 161–168. doi: 10.13465/j.cnki.jvs.2016.22.024
    HUI X L, MU R K, BAI C Y, et al. Dynamic mechanical property and constitutive model for TC4 titanium alloy [J]. Journal of Vibration and Shock, 2016, 35(22): 161–168. doi: 10.13465/j.cnki.jvs.2016.22.024
    [11] KAY G. Failure modeling of titanium 6Al-4V and aluminum 2024-T3 with the Johnson-Cook material model: DOT/FAA/AR-03/57 [R]. Washington: U.S. Department of Transportation, Federal Aviation Administration, 2003: 1–11.
  • [1] 胡八一董庆东韩长生张海平 . TC4钛合金及40Cr钢破片中绝热剪切带的TEM分析. 高压物理学报, 1996, 10(1): 37-43 . doi: 10.11858/gywlxb.1996.01.006
    [2] 龚自正傅世勤经福谦 . 一维应变压缩下TC4钛合金绝热剪切研究. 高压物理学报, 1992, 6(3): 180-191 . doi: 10.11858/gywlxb.1992.03.004
    [3] 张志春强洪夫高巍然 . 光滑粒子流体动力学有限元法接触算法研究. 高压物理学报, 2011, 25(2): 97-103 . doi: 10.11858/gywlxb.2011.02.001
    [4] 孙欢腾李名锐周刚马坤舒孝鸿 . 杆式钨合金弹超高速撞击薄靶的能量损耗. 高压物理学报, 2019, 33(6): 064106-1-064106-9. doi: 10.11858/gywlxb.20190732
    [5] 强洪夫张国星王广黄拳章 . SPH方法在宽速域岩石侵彻问题中的应用. 高压物理学报, 2019, 33(5): 055105-1-055105-9. doi: 10.11858/gywlxb.20180621
    [6] 李晓杰莫非闫鸿浩张程娇 . 爆炸焊接斜碰撞过程的数值模拟. 高压物理学报, 2011, 25(2): 173-176 . doi: 10.11858/gywlxb.2011.02.014
    [7] 缪广红王章文李亮江向阳刘文震程扬帆汪泉余勇马宏昊沈兆武 . 爆炸复合边界效应的数值模拟. 高压物理学报, 2017, 31(1): 93-96. doi: 10.11858/gywlxb.2017.01.014
    [8] 强洪夫孙新亚王广陈福振 . 混凝土薄板侵彻贯穿问题的SPH数值模拟. 高压物理学报, 2019, 33(2): 024101-1-024101-10. doi: 10.11858/gywlxb.20180634
    [9] 庄仕明丰树平王春彦孙承纬 . 高应变率下TC4及TC9钛合金的动态断裂. 高压物理学报, 1995, 9(2): 96-106 . doi: 10.11858/gywlxb.1995.02.003
    [10] 王春奎刘小苹郑融 . LY-12铝高温凝聚态动力学性质研究高温冲塞剪切破坏强度的测量. 高压物理学报, 1992, 6(2): 108-115 . doi: 10.11858/gywlxb.1992.02.004
    [11] 陈军田占东曾代朋张震宇段卓平 . PBX 9501炸药冲击起爆初期反应流体动力学特性研究. 高压物理学报, 2013, 27(3): 331-336. doi: 10.11858/gywlxb.2013.03.003
    [12] 马天宝费广磊张文耀 . 三维多物质弹塑性流体动力学Euler方法的并行算法研究及程序测试. 高压物理学报, 2011, 25(6): 508-513. doi: 10.11858/gywlxb.2011.06.005
    [13] 陈浩陈雄周长省薛海峰 . 超音速等离子体炬的磁流体动力学数值研究. 高压物理学报, 2015, 29(3): 185-190. doi: 10.11858/gywlxb.2015.03.004
    [14] 王春奎 . LF-6及LY-12铝高温凝聚态动力学性质研究──高温剪切模量的测量. 高压物理学报, 1994, 8(3): 213-219 . doi: 10.11858/gywlxb.1994.03.009
    [15] 张树华 . TC4、16Mn合金及Al2O3陶瓷的高温弹性模量. 高压物理学报, 1995, 9(2): 133-138 . doi: 10.11858/gywlxb.1995.02.008
    [16] 张树华 . 钛合金(TC-6)二阶和三阶弹性常数. 高压物理学报, 1989, 3(2): 172-175 . doi: 10.11858/gywlxb.1989.02.012
    [17] 邹学韬张晓晴姚小虎 . 压剪载荷作用下TB6钛合金的动态力学性能. 高压物理学报, 2019, 33(2): 024206-1-024206-9. doi: 10.11858/gywlxb.20190713
    [18] 李健郭晓璇马胜国李志强辛浩 . AlCrFeCuNi高熵合金力学性能的分子动力学模拟. 高压物理学报, 2020, 34(1): 011301-1-011301-9. doi: 10.11858/gywlxb.20190762
    [19] 孙立志李治源吕庆敖 . 四层剪切运动间隔靶板对射流干扰的数值模拟. 高压物理学报, 2013, 27(3): 423-430. doi: 10.11858/gywlxb.2013.03.017
    [20] 韩德才周铁农潘永阁 . 库特(Couette)剪切流准确解的一般形式. 高压物理学报, 2000, 14(3): 229-234 . doi: 10.11858/gywlxb.2000.03.013
  • 加载中
图(13)表(3)
计量
  • 文章访问数:  529
  • 阅读全文浏览量:  362
  • PDF下载量:  5
出版历程
  • 收稿日期:  2020-02-24
  • 录用日期:  2020-03-17
  • 网络出版日期:  2020-07-09
  • 刊出日期:  2020-08-05

鸟撞冲击下TC4钛合金平板的变形和破坏

    作者简介:贾 林(1988-),男,硕士,工程师,主要从事发动机强度试验技术研究. E-mail:jialin20060598@126.com
    通讯作者: 姚小虎, yaoxh@scut.edu.cn
  • 1. 中国航发上海商用航空发动机制造有限责任公司,上海 201108
  • 2. 华南理工大学土木与交通学院,广东 广州 510641

摘要: 为了分析航空用TC4钛合金板材在鸟撞冲击下的动态响应,通过3D-DIC动态变形场测试技术,研究了鸟撞冲击过程中TC4钛合金平板变形场,并基于显式有限元分析软件ABAQUS,建立了鸟撞数值计算模型。将考虑拉压不对称性的修正von Mises屈服准则引入TC4钛合金材料的Johnson-Cook动态本构模型和动态损伤模型中,采用光滑粒子流体动力学方法(SPH)建立了鸟体模型。数值计算和鸟撞实验的对比结果表明,计算应变与实验结果吻合很好,验证了TC4钛合金鸟撞冲击数值分析模型的合理性和可靠性。

English Abstract

  • TC4钛合金是一种$\alpha + \beta $型中强度钛合金,具有较高的强度和较优异的塑性,在航空、航天、船舶以及兵器领域应用广泛[1]。鸟撞问题是飞机结构在起降过程中面临的主要威胁之一。鸟撞发生时,鸟体在毫秒级时间内瞬间冲击机体结构。高速冲击产生的巨大能量将导致机身结构严重损伤,从而引发伤亡事故。因此,国际适航标准要求所有向前部件在使用前必须分析其抗鸟撞性能。

    大量的研究表明,鸟体在高速冲击作用下表现出明显的流体流动飞溅特性。近年来,诸多学者针对飞机结构的抗鸟撞性能开展了大量研究工作。普遍认为,高速鸟撞冲击问题是一个应变率相关的流固耦合问题。目前结构抗鸟撞性能的数值分析方法主要有3种:拉格朗日有限元法(Lagrangian finite element)、任意拉格朗日-欧拉法(Arbitrary Lagrange-Euler,ALE)和光滑粒子流体动力学方法(Smooth particle hydro-dynamic,SPH)。采用拉格朗日有限元法分析高速鸟撞冲击时,因鸟体结构变形大,致使单元发生畸变,故拉格朗日有限元法只适合模拟低速鸟撞。SPH方法是一种基于拉格朗日技术的自适应无网格粒子法,将其与有限元方法进行耦合,可在流固耦合问题求解中展现显著的优势[2-3]。由于SPH粒子在空间相互独立,因此SPH法比拉格朗日有限元法和ALE法更适于解决高速鸟撞冲击问题[4]。例如:刘军等[5]通过对比鸟撞平板叶片实验结果和SPH法及拉格朗日有限元法数值分析结果,发现SPH方法与实验结果更接近;刘富等[6]采用SPH方法进行了2024-T3铝合金平板抗高速鸟撞冲击性能研究,得到了与实验结果相近的模拟结果;Liu等[7]通过不同速度的平板鸟撞冲击实验和数值分析,研究了适用于不同鸟撞速度的鸟体模型;姚小虎等[8]通过鸟撞圆弧风挡实验和数值计算,分析了风挡玻璃在鸟撞冲击过程中的损伤破坏。

    本研究采用三维图像相关法(3 dimensional digital correlate,3D-DIC),分析TC4钛合金平板高速鸟撞过程中的变形场,基于SPH方法和TC4钛合金的Johnson-Cook动态损伤模型,建立TC4钛合金平板鸟撞数值模型,并将模拟结果与鸟撞实验进行对比验证。

    • 鸟撞实验装置由鸟弹发射系统、TC4钛合金靶板、速度测试系统、照明系统和高速摄像系统组成。实验装置如图1所示。本实验使用的鸟弹为长L = 228 mm、直径D = 114 mm的明胶弹。鸟弹由空气炮发射,利用激光测速仪记录发射速度,激光测速仪的系统误差小于0.5%。为了解析TC4钛合金靶板背面的三维变形场,在靶板背面设置两台I-SPEED 716型高速摄影机,拍摄帧率设置为104帧每秒。位于靶板正面的两台SA-X型高速摄影机记录鸟弹飞行轨迹和撞击靶板时的响应,保证鸟弹垂直撞击TC4钛合金靶板。实验开始之前,进行调焦、视场校准和同步设置。将4台高速摄影机的触发开关通过BNC线引至操作间,其中用于动态3D-DIC测量的两台相机使用转接头连接,以实现同步触发。高速摄影机布局如图2所示。

      图  1  实验装置示意图

      Figure 1.  Schematic of experimental apparatus

      图  2  高速摄影机布局

      Figure 2.  High speed camera arrangement

      试件材料为TC4钛合金平板,尺寸为600 mm × 600 mm × 1.6 mm。通过均匀分布的16颗M10螺栓及4.0 mm厚的夹具,将试件固定在试验工装上,夹具尺寸与螺栓分布如图3所示。

      图  3  靶板尺寸

      Figure 3.  Size of target board

      鸟撞实验共设3个发射速度,分别为149、167和180 m/s。每组进行4次重复实验。图4显示了3种速度工况下鸟撞实验结果。图4中第1行的3幅图为平板正面高速摄影图像,可以看出:鸟弹包裹在弹托中由炮管发射,在空气阻力和实验舱入射口的作用下,鸟弹和弹托在撞击TC4钛合金平板前完全分离。弹托保证了鸟弹在发射过程中的整体形状和结构不受炮管内高压气体的破坏,弹托与鸟弹的完全分离消除了弹托对TC4平板鸟撞响应的影响。图4中第2行和第3行图像分别显示了TC4钛合金平板的正面和背面鸟撞冲击结果。发射速度为149 m/s的4次实验中,TC4钛合金平板均未发生破坏;发射速度为167 m/s的4次实验中,2次发生破坏,2次未发生破坏;而发射速度为180 m/s的4次实验中,平板均发生破坏。

      图  4  149、167和180 m/s的鸟撞实验结果

      Figure 4.  Results of bird strike experiments at the speed of 149, 167 and 180 m/s

      图5为TC4钛合金平板破坏照片。鸟体撞击平板后产生的拉伸波向外传播,在螺栓处产生剪切作用,平板发生了剪切破坏。

      图  5  TC4钛合金平板破坏照片

      Figure 5.  Failure of TC4 titanium alloy plate

    • 数值计算采用的鸟体几何模型与实验相同,为两端半球状、中间圆柱体的胶囊状柱体,长径比L/D = 2,如图6所示。鸟体模型的质量为1.8 kg。采用SPH单元模拟高速鸟撞冲击过程中的鸟体流体状飞溅,鸟体材料参数列于表1

      Density/(kg·m−3Elastic modulus/GPaPoisson’s ratioYield stress/MPaFailure strainTangent modulus/MPa
      9280.0680.490.691.255

      表 1  鸟体材料参数

      Table 1.  Material parameters of bird body

      图  6  鸟弹的几何尺寸

      Figure 6.  Geometry of bird ball

      高速鸟撞实验过程中,靶板夹具和支撑架的刚度足够大,夹具和支撑架只发生线弹性变形,因此采用钢材的线弹性本构模型描述。TC4钛合金平板在高速鸟撞冲击载荷作用下发生了大变形和损伤破坏。实验发现,TC4钛合金平板的主要破坏形式是剪切破坏,因此在数值仿真计算中需要考虑剪应力的影响。大量实验表明,钛合金材料具有拉压不对称性,需要对von Mises屈服准则进行修正。本研究将Johnson-Cook动态本构模型和Johnson-Cook损伤失效模型引入邹学韬等[9]提出的von Mises修正本构框架中。该本构可以表征TC4钛合金在强冲击载荷作用下的塑性流动应力和损伤破坏行为。Johnson-Cook动态本构模型的表达式为

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

      式中:${\sigma _{\rm{s}}}$为塑性流动应力;A为参考应变率下的屈服应力;Bn为应变强化系数;$\varepsilon $为等效塑性应变;C为应变率敏感系数;${\dot \varepsilon ^*}{\rm{ = }}\dot \varepsilon /{\dot \varepsilon _0}$为无量纲应变率,${\dot \varepsilon _0}$为参考应变率;${T^*}{\rm{ = (}}T - {T_{\rm{r}}})/({T_{\rm{m}}} - {T_{\rm{r}}})$为无量纲温度,T为温度,${T_{\rm m}}$为材料熔化温度,${T_{\rm{r}}}$为参考温度;m为温度软化系数。

      考虑到TC4钛合金材料的拉压不对称性,引入拉压不对称因子$G(\sigma )$修正von Mises屈服面,得到屈服函数为

      $\phi {\rm{ = }}f(\sigma )G(\sigma ) = 1$

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

      $G(\sigma ) = \exp [ - c(\xi + 1)] = 1$

      式中:$f(\sigma ) = 1$ 为von Mises屈服面函数;$\xi =\cos (3\theta ) = \dfrac{{27}}{2} \cdot \dfrac{{{J_3}}}{{{{(3{J_2})}^{3/2}}}}$为Lode参数,其中$\theta $为Lode角,J2为偏应力第二不变量,J3为偏应力第三不变量;$c$为不对称系数,可以通过两种简单应力状态求解。本研究中,不对称系数$c$

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

      式中:${\tau _{\rm{s}}}$为纯剪切加载的屈服应力。于是,根据相关联的流动法则,即屈服面函数等于塑性势函数$g{\rm{(}}\sigma {\rm{)}}$,得到该本构的增量表达式

      $\left\{ \begin{aligned} & {\rm{d}}{\sigma _{ij}} = {{{D}}_{\rm{e}}}{\rm{(d}}{\varepsilon _{ij}} - {\rm{d}}\varepsilon _{ij}^{\rm{p}}) \\ & {\rm{d}}\varepsilon _{ij}^{\rm{p}} = {\rm{d}}\lambda \dfrac{{\partial g}}{{\partial {\sigma _{ij}}}} \\ & {\rm{d}}\lambda =\dfrac{{g\left( {\sigma _{ij}^{{\rm{trail}}}} \right) + \dfrac{{\partial g}}{{\partial \dot \varepsilon }}{\rm{d}}\dot \varepsilon }}{{\dfrac{{\partial g}}{{\partial {\sigma _{ij}}}}{Y_{ij}} - \dfrac{{\partial g}}{{\partial \sigma _{{\rm{eq}}}^{\rm{p}}}}H}} \\ & {Y_{ij}} = {{{D}}_{\rm{e}}}\dfrac{3}{{{\sigma _{\rm{s}}^2}}}{\rm{exp}}\left[ {2\ln \alpha (\xi + 1)} \right]\left( {\dfrac{{\partial {J_2}}}{{\partial {\sigma _{ij}}}} + 2{J_2}\ln \alpha \dfrac{{\partial \xi }}{{\partial {\sigma _{ij}}}}} \right) \\ & H=\sqrt {\dfrac{2}{3}\dfrac{{\partial g}}{{\partial {\sigma _{ij}}}}\dfrac{{\partial g}}{{\partial {\sigma _{ij}}}}} \end{aligned} \right.$

      式中:De为弹性矩阵,${\rm{d}}\lambda $为塑性流动因子,$\varepsilon _{ij}^{\rm{p}}$为塑性应变,$\sigma _{ij}^{{\rm{trail}}}$为试探应力,$\sigma _{{\rm{eq}}}^{\rm{p}}$为等效塑性应力。

      Johnson-Cook损伤失效模型为

      ${\varepsilon _{\rm{f}}}{\rm{ = [}}{D_1}{\rm{ + }}{D_2}\exp ({D_3}{\sigma ^*})](1 + {D_4}\ln {\dot \varepsilon ^*})(1 + {D_5}{T^*})$

      式中:${\varepsilon _{\rm{f}}}$为失效应变;D1D5为材料参数;${\sigma ^*}{\rm{ = }}p/{\sigma _{\rm{eq}}}{\rm{ }}$,其中p为静水压力,${\sigma _{\rm{eq}}}$为等效应力。

      数值计算所使用的本构模型参数列于表2[10-11],其中E为弹性模量,ρ为密度,μ为泊松比。

      ρ/(g·cm−3)μTm/KA/MPaB/MPanCm
      4.4300.331 878106010900.8840.01171.1
      E/GPa${\dot \varepsilon _{{0}}}$/s−1D1D2D3D4D5
      1354 × 10−4−0.0900.2700.4800.0143.870

      表 2  TC4钛合金材料参数[10-11]

      Table 2.  Parameters of TC4 titanium alloy[10-11]

      TC4钛合金平板高速鸟撞的数值计算有限元模型如图7所示。TC4靶板、夹具和M10螺栓均采用C3D8R六面体八节点减缩积分单元模拟。通过建立一般接触,计算鸟体撞击TC4钛合金靶板以及螺栓和靶板之间的接触。夹具通过16颗M10螺栓固定在支架上,在数值计算中对螺栓进行固支约束。鸟体速度分别设置为149、167和180 m/s。

      图  7  鸟撞数值计算模型

      Figure 7.  Numerical model of bird strike

      在TC4钛合金平板上选取6个具有代表意义的观测点,如图8所示,其中观测点S1S2S3沿轴向分布,S4S5S6沿对角线方向分布。

      图  8  观测点布局

      Figure 8.  Distribution of observation points

    • 图9为鸟撞速度为149 m/s时TC4钛合金平板的等效应力云图。鸟体撞击平板后,鸟体前端受到冲击压缩后解体并呈流体状飞溅,鸟体后端仍保持固体状态。平板受鸟体冲击后产生拉伸波,并向平板四周传播。1.66 ms时鸟体完全解体,鸟体撞击的冲击能量完全耗散,此时TC4钛合金平板的应力、应变和位移达到最大值,随后开始一定程度回弹。

      图  9  149 m/s鸟撞等效应力云图

      Figure 9.  Equivalent stress nephograms of bird impacting with the velocity of 149 m/s

      图10为鸟撞击后TC4钛合金平板变形的数值计算结果和3D-DIC实验结果对比。图10中149 m/s和167 m/s工况下的最大位移(Smax)图像选自未破坏实验,180 m/s工况下的最大位移图像选自平板破坏飞出前(2.00 ms前)。鸟撞过程中,TC4钛合金平板的变形大,对角线方向隆起,隆起处亮度明显增大,使得平板部分区域被遮挡,同时也遮挡了高速摄影机,因此出现部分区域未追踪到变形场的问题。由图10可知,计算得到的最大位移场与实验结果吻合较好。3种工况下数值仿真和实验得到的观测点最大位移如表3所示。

      Velocity/(m·s−1)MethodMaximum displacement/mm
      S1S2S3S4S5S6
      149Sim.655136645523
      Exp.686038635325
      167Sim.756048766527
      Exp.805947766324
      180Sim.1191016311810354
      Exp.11790531039550

      表 3  TC4钛合金平板鸟撞最大位移

      Table 3.  Maximum displacement of titanium alloy plate impacted by a bird

      图  10  TC4钛合金平板鸟撞变形结果

      Figure 10.  Deformation of TC4 titanium alloy plate under the bird impact

      图11为计算得到的180 m/s工况下TC4钛合金平板破坏过程中的等效塑性应变云图。从图11可以看出:0.40 ms时,位于轴线上的4颗螺栓附近开始出现裂纹;0.68 ms时,平板与夹具接触处进入塑性阶段;1.04 ms时,平板对角线和夹具接触处开始起裂,并沿着夹具边缘和对角线方向扩展;2.00 ms时,最先起裂的4颗螺栓孔处裂纹贯穿。对比可见,计算得到的螺栓孔处的损伤和破坏形式与实验结果基本一致。

      图  11  180 m/s工况下计算得到的TC4钛合金平板破坏过程

      Figure 11.  Failure process of TC4 titanium alloy plate calculated at 180 m/s

      计算与实验得到的位移-时间曲线对比如图12所示。从图12中可以看出,计算得到的位移变化趋势及大小与实验结果基本吻合,表明本研究使用的Johnson-Cook动态本构和损伤失效模型对于模拟TC4钛合金高速鸟撞冲击问题是比较准确的。

      图  12  位移-时间曲线的计算和实验结果对比

      Figure 12.  Comparison of calculated displacement-time curves with experimental results

      图13对比了6个观测点的应变时程曲线。数值计算得到的6个观测点应变与实验数据的整体吻合度较高。从图13中可以看出,最靠近鸟撞点的观测点S1S4的等效应变在0~0.2 ms内增大,0.2~1.0 ms内保持平稳,1.0~1.4 ms再次增大,1.4 ms后再次保持不变,呈现双台阶模式。其余观测点均未表现出此双台阶模式。观测点S1S2的等效应变出现双台阶的原因在于这两个点位于鸟弹半径范围之内。鸟弹撞击TC4钛合金平板瞬间,应变瞬间增大;0.2~1.0 ms内应力波向边界传递并在边界处反向,此时S1S2区域内材料包裹着鸟弹运动,因此应变出现平台段;1.0 ms时,边界反射的应力波再次到达S1S2区域,使得应变再次增大。

      图  13  应变-时间曲线的计算与实验结果对比

      Figure 13.  Comparison of calculated strain-time curves with experimental results

    • 通过3D-DIC实验和数值计算方法,研究了1.8 kg鸟体高速撞击1.6 mm厚TC4钛合金平板的动态响应和损伤破坏,得到了较精确、有效的有限元模型,并得到如下结论。

      (1)1.6 mm厚的TC4钛合金在1.8 kg鸟体高速撞击下的临界破坏速度为167 m/s。撞击过程中平板内部未破坏,而螺栓和夹具处发生剪切破坏。

      (2)3D-DIC测试技术能够比较准确地测定鸟撞冲击过程中TC4钛合金平板的变形场。高速冲击过程中平板的变形较大,易出现光线遮挡和反光,需要设置补充高速摄影机。

      (3)实验表明,鸟撞冲击后TC4钛合金平板破坏主要为螺栓等边界处的剪切破坏。将修正的von Mises屈服准则引入Johnson-Cook动态本构和损伤模型中,在本构中同时考虑拉伸和剪切两种应力状态。该本构能够准确地模拟鸟撞平板问题。

参考文献 (11)

目录

    /

    返回文章
    返回