航空发动机空心风扇叶片扩散连接焊缝建模与优化分析技术

航空发动机空心风扇叶片扩散连接焊缝建模与优化分析技术The Modeling and Optimization for Diffusion Bonding Seam of Hollow Fan Blade in AeroEngine柴象海1,2，侯亮1,2（1.中航商用航空发动机有限责任公司设计研发中心，上海市200241；2.上海商用飞机发动机工程技术研究中心，上海市200241）摘要: 扩散连接钛合金空心瓦伦结构已经被成功地应用于航空发动机零部件，如风扇叶片、导流叶片等结构。

研究发现，瓦伦结构的几何特征对零件的抗冲击强度，例如抗鸟撞性能，有很大的影响，通过优化设计空心瓦伦结构来提升零件的抗冲击性能一直是工业界追求的目标。

本文针对典型的用于发动机空心风扇叶片的三层板空心瓦伦结构，采用HyperWorks工具包，通过基于试验的仿真优化研究了几个重要的瓦伦结构特征参数对抗冲击性能的影响。

首先通过有限元数值模拟和逆向分析的手段估算出钛合金空心瓦伦焊缝的失效强度。

然后通过数值模拟的结果得到了相邻瓦伦夹角、焊缝长度与扩散连接焊焊缝在冲击载荷下失效强度之间的关系。

本项研究可以为钛合金空心瓦伦结构零部件设计和抗冲击强度校核提供参考。

关键词:扩散连接；碰撞试验；HyperWorks建模优化；鸟撞；瓦伦结构Abstract: Diffusion-bonded titanium hollow components with a Warren girder internal structure, such as fan blade, OGV, etc., has been successfully used on aircraft engines. It was found that the geometric features of such Warren girder structure have significant effect on the impact strength of the component. Such design has provided the possibility for the designer to improve bird-strike resistance of the component by optimizing its internal Warren girder structure based HyperWorks. At first, a set of ballistic impact tests with titanium hollow panels with various internal geometric parameters were conducted, and the failure stresses of the diffusion-bonding area of the hollow structure were estimated through an inverse method based on test-analysis correlation. Then, a qualitative relation between the failure stress of the diffusion bonds and two important geometric parameters, skew angle of the girder and length of the diffusion bond, was established through numerical simulations. This study provided useful reference for the optimal design of components with Warren girder hollow structures.Key words:diffusion bonding; impact test; HyperWorks simulation optimizasition; bird-strike; warren structure1 概述伴随着涡轮喷气式发动机的研究与发展，扩散连接钛合金空心瓦伦结构零部件在国内外各大航空发动机中得到了越来越广泛的应用。

20世纪80年代，罗·罗公司采用钛合金三层板结构的超塑成形/扩散连接(SPF/DB)组合制造工艺的宽弦空心风扇叶片研制成功，并于1992年成功完成了在航空发动机上的整机试验。

与带凸肩的窄弦实心风扇叶片相比，宽弦空心风扇叶片的成功应用，使得风扇转子叶片数量减少了大约1/3，转子重量也减轻了10%～30%，蜂窝夹芯板的空腔结构改善了叶片的抗振性能，提高了叶片的气动效率，增加了叶片的喘振裕度。

此外，伴随着空心瓦伦结构零部件的深入应用，研究者发现在减轻发动机重量的同时，还可以通过空心结构的优化设计改善零部件强度性能。

对于空心风扇叶片而言，通过优化设计空心瓦伦结构来提升其抗鸟撞性能一直是工业界追求的目标[1]。

焊缝失效是空心瓦伦结构在冲击载荷下最重要失效模式之一。

理想情况下，扩散焊的焊缝强度与母体材料是非常接近的，但在实际工程应用中有诸多因素对焊缝强度产生影响[2]。

首先是实际工艺过程中难以保证扩散连接部位100%的焊合率，而针对某一特定的工艺流程，尽管这些空心瓦伦结构的焊合率实际已经在产品的检测许可范围内，不同空心瓦伦结构也会导致实际焊合率的不同。

另一方面，空心瓦伦结构的几何特征本身也会影响到焊缝边缘的应力集中情况，进而影响焊缝的失效模式[3-6]。

因此，空心瓦伦结构的抗冲击失效强度是工艺过程与结构几何特征综合作用的结果。

对于设计人员来说，找到一个相对简单有效的空心瓦伦结构参数与抗冲击强度的关联关系对于提升设计水平无疑是十分重要的。

本文首先针对发动机零部件所常用的钛合金三层空心瓦伦结构，采用HyperMesh建模工具，建立含有扩散连接焊缝的空心平板数值仿真模型。

然后，以焊缝强度作为设计变量，通过尺寸优化，得到钛合金空心瓦伦焊缝失效强度。

最后，对焊缝强度进行分析，得到在冲击载荷下相邻瓦伦夹角、焊缝长度以与扩散连接焊焊缝失效强度之间的关系。

2 扩散连接钛合金空心瓦伦结构建模空心瓦伦结构的抗冲击强度很大程度上取决于焊缝的强度。

而后者与空心瓦伦结构的几何特征紧密相关。

为了分析随着几何特征变化而产生的焊缝强度的变化规律，需要考虑焊缝长度、沿弦长方向的几何形状、以及超塑成形对焊缝受力特性的影响等因素，本文所定义的钛合金空心瓦伦结构的几何特征参数如图1所示。

图1 空心风瓦伦结构几何特征参数(1) 空心瓦伦结构的面板厚度T1；(2) 成形后的瓦伦芯板厚度T3；(3) 相邻瓦伦的夹角α；(4) 芯板与面板的单个扩散连接部位长度D。

本文根据超塑成形/扩散连接工艺所能加工的钛合金空心瓦伦结构参数范围，选取不同空心瓦伦结构参数，设计了等长度直线焊缝钛合金三层空心平板冲击试验件。

试验件采用超塑成形/扩散连接成形工艺加工。

空心平板试验件如图2所示，由上面板、芯板和下面板组成。

图2 钛合金三层空心平板试验件本文三层瓦伦结构空心平板采用8节点实体模型。

单元应变采用泰勒级数展开算法求解，并且通过特定的处理，使应变分量避免剪切自锁。

仿真模型考虑变形过程中的材料硬化现象，同时引入结构变形时材料高应变率效应，并采用RBE2单元仿真焊缝连接。

通过该方法建立的空心平板有限元模型，能够考虑工件结构对焊缝强度的影响，对不同结构下焊缝失效临界载荷进行定量设置。

通过试验对网格密度、材料失效参数和焊缝失效强度进行标定，能够有效提高仿真分析精度和仿真结果的可靠性。

图3为等长度焊缝的钛合金三层瓦伦结构空心平板有限元计算模型，面板两层单元，芯板一层单元，模型包含8976 单元，9635 节点。

图3 空心平板有限元模型本文采用显式分析求解器计算空心平板模型瞬态冲击问题，可以处理两个部件之间发生接触破坏并扩展失效的情况。

3扩散连接焊缝强度优化分析3.1 扩散连接钛合金空心瓦伦结构试验件高速冲击试验由于超塑成形/扩散连接瓦伦结构焊缝强度影响因素的复杂性，本文通过对空心平板在不同载荷工况下进行冲击试验，对焊缝在高速冲击载荷下脱焊失效的机理进行了研究。

试验设备采用150mm 口径高速气炮。

平板冲击试验针对三种不同的空心结构进行。

试验中，空心平板以不同的初始速度垂直撞击铝合金靶板。

空心平板在冲击载荷作用下，呈现出沿焊缝撕裂，面板向后卷曲的典型损伤情况，如图4所示。

图4 钛合金三层空心平板试验件冲击后的损伤情况3.2 空心瓦伦结构焊缝失效参数优化分析空心瓦伦结构焊缝失效应力的确定是本文工作的重点。

由于结构的特殊性，瓦伦结构的焊缝强度难以通过简单的试验直接测量。

本文首先对不同试验工况下空心平板高速冲击试验建立数值仿真模型，以焊缝强度作为设计变量，通过尺寸优化，拟合出钛合金空心瓦伦焊缝的失效参数，即试验件对应的扩散连接焊缝拉伸失效应力和剪切失效应力。

3.3.1 优化分析计算根据工程经验，选择芯板角度、芯板/面板厚度、非空心段焊缝长度作为优化变量。

1）基于空心平板变形的优化目标函数选择焊缝拉伸失效应力和剪切失效应力作为寻优的输入变量，空心平板试验件的变形作为寻优目标函数的构成要素，通过迭代计算使得模拟分析结果与空心平板试验变形形态相一致，逆向得到焊缝缝拉伸失效应力和剪切失效应力数值。

寻优问题由以下公式描述：m k y y y x F k obj obj k k ,,2,1,)(2K =⎟⎟⎠⎞⎜⎜⎝⎛−=∑ω (1)n k x x x f y n k L L 2,1),,(21==(2)n i c x g i i K ,2,1,)(=≤(3)p i s x s m u i i m l i K ,2,1,)()(=≤≤ (4)公式（1）表示目标函数，由空心平板不同位置变形量累加得到。

公式（2）描述了数值模拟输入变量与空心平板变形量之间的关系，通过手动分析调试，可以固化材料模型、网格密度、接触刚度等输入变量对模拟结果的影响，而仅选择空腔处焊缝拉伸失效应力1σ和剪切失效应力1τ，前后缘处焊缝拉伸失效应力2σ和剪切失效应力2τ作为输入变量。

公式（3）表示约束函数，公式（4）给出了输入变量取值范围。

其中，obj y 是由试验得到的空心平板变形量，k y 是数值模拟分析结果，k ω是权重。

)(m l i s and)(m u i s 是输入变量基于工程经验得到的上下限。

当焊缝拉伸和剪切失效参数变化时，k y 即空心平板试验件变形量随之改变。