贾 林，李从富，邹学韬，姚小虎

（1. 中国航发上海商用航空发动机制造有限责任公司，上海 201108；2. 华南理工大学土木与交通学院，广东 广州 510641）

TC4 钛合金是一种 α+β型中强度钛合金，具有较高的强度和较优异的塑性，在航空、航天、船舶以及兵器领域应用广泛[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 钛合金平板鸟撞数值模型，并将模拟结果与鸟撞实验进行对比验证。

1 钛合金平板鸟撞实验

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

图 1 实验装置示意图Fig. 1 Schematic of experimental apparatus

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

图 2 高速摄影机布局Fig. 2 High speed camera arrangement

图 3 靶板尺寸Fig. 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 的鸟撞实验结果Fig. 4 Results of bird strike experiments at the speed of 149, 167 and 180 m/s

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

2 钛合金平板鸟撞数值计算

2.1 计算模型

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

图 5 TC4 钛合金平板破坏照片Fig. 5 Failure of TC4 titanium alloy plate

图 6 鸟弹的几何尺寸Fig. 6 Geometry of bird ball

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

表 1 鸟体材料参数Table 1 Material parameters of bird body

表 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。

在TC4 钛合金平板上选取6 个具有代表意义的观测点，如图8 所示，其中观测点S1、S2和S3沿轴向分布，S4、S5和S6沿对角线方向分布。

图 7 鸟撞数值计算模型Fig. 7 Numerical model of bird strike

2.2 计算结果

图 8 观测点布局Fig. 8 Distribution of observation points

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

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

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

图 9 149 m/s 鸟撞等效应力云图Fig. 9 Equivalent stress nephograms of bird impacting with the velocity of 149 m/s

图 10 TC4 钛合金平板鸟撞变形结果Fig. 10 Deformation of TC4 titanium alloy plate under the bird impact

表 3 TC4 钛合金平板鸟撞最大位移Table 3 Maximum displacement of titanium alloy plate impacted by a bird

图 11 180 m/s 工况下计算得到的TC4 钛合金平板破坏过程Fig. 11 Failure process of TC4 titanium alloy plate calculated at 180 m/s

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

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

图 12 位移-时间曲线的计算和实验结果对比Fig. 12 Comparison of calculated displacement-time curves with experimental results

图 13 应变-时间曲线的计算与实验结果对比Fig. 13 Comparison of calculated strain-time curves with experimental results

3 结 论

通过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 动态本构和损伤模型中，在本构中同时考虑拉伸和剪切两种应力状态。该本构能够准确地模拟鸟撞平板问题。