高温环境下发动机叶片-硬涂层阻尼结构动力学建模与分析.pdf

第42 卷第1期2024年2 月中国民航大学学报JOURNAL OF CIVIL AVIATION UNIVERSITY OF CHINAVol.42 No.1February2024高温环境下发动机叶片-硬涂层阻尼结构动力学建模与分析李双宝,李芳b（中国民航大学a.科技创新研究院；b.中欧航空工程师学院，天津30 0 30 0）摘要：硬涂层作为一种阻尼结构可以有效地抑制发动机叶片的有害振动，其动力学研究对于设计和优化复合结构有着重要的意义。本文探讨了高温环境下发动机叶片-硬涂层阻尼结构的动力学模型，对其进行动力学分析，研究固有频率、振动响应幅值等动力学特性。采用数值模拟分析硬涂层、温度对复合结构动力学特性影响。利用有限元仿真对结果进行了验证，得到硬涂层对于发动机叶片有减振作用，在高温环境下依然有效。关键词：高温环境；发动机叶片；硬涂层；动力学模型中图分类号：V231;TB535.1Dynamic modeling and analysis of engine blade-hard coating damping structure(a.Science and Technology Innouation Division;b.Sino-European Institute of Aviation Engineering,CAUC,Tianjin 300300,China)Abstract:Hard coating is a kind of damping structure can effectively restrain harmful vibration of engine blade,and its dy-namic research is of paramount significance for the design and optimization of composite structure.The dynamicmodel of engine blade-hard coating damping structure under high temperature environment is discussed,and con-ducts dynamic analysis on it to study the dynamic characteristics such as inherent frequency,vibration responseamplitude,etc.The numerical simulation is adopted to analyze the impact of hard coating and temperature on dy-namic characteristics of composite structure.The finite element simulation is adopted for verification of results toconclude that the hard coating can reduce vibration of engine blade,which is still effective under high temperatureenvironment.Key words:high temperature environment;engine blade;hard coating;dynamic model文献标志码：Aunder high temperature environmentLI Shuangbao,LI Fangh文章编号：16 7 4-559 0(2 0 2 4)0 1-0 0 30-10有关研究表明叶片作为航空发动机的核心零部件,其由振动造成的疲劳失效占失效总故障的6 8%。近几年研究发现，硬涂层阻尼结构可以应用于叶片减振，表现为其可以改变结构固有频率从而使其远离共振带并且减小共振应力和响应幅值2-4。硬涂层是指金属基、陶瓷基或金属陶瓷混合的复合涂层,其有较高硬度的同时还有耐高温的特性5,所以将硬涂层应用于结构减振是可行的。最初的硬涂层研究聚焦在飞机蒙皮、叶片等薄壁结构6-8 上。张波成等9 研究了硬涂层阻尼结构用于飞行器蒙皮的减振设计及优化方法，研究指出通过对硬涂层合理设计可以获得最优的减振性能。Luo等0 研收稿日期：2 0 2 2-11-15；修回日期：2 0 2 3-0 3-15作者简介：李双宝（19 7 8 一），男，河北隆化人，教授，博士，研究方向为非线性系统动力学究了硬涂层复合板的优化方法，获得了阻尼材料的最优布局。这些研究主要证实了硬涂层对于薄壁结构的减振作用并设计了优化方法。之后，也有一些学者进行了硬涂层对圆环薄壁结构影响的研究。Song等分析了不同环形区域上硬涂层的振动和阻尼效应，得到涂敷前后固有频率变化趋势相同，但硬涂层覆盖位置对固有频率有不同影响的结论。Sun等12 研究了硬涂层对悬臂圆柱壳的影响，分析表明，加人硬涂层结构呈现固有频率降低、模态阻尼比升高以及共振响应降低的效果。以上研究表明硬涂层对多种薄壁结构均具有阻尼效应。随着发动机叶片所处环境愈发恶劣，研究高温环基金项目：国家自然科学基金项目（12 17 2 37 6）(9)第42 卷第1期境下叶片振动是必不可少的。胡君逸等13研究了热环境下四边固支板的固有特性和激励响应。Zhang等4研究了双曲面壳在热环境中的振动特性，发现双曲面壳的固有频率随温度的增加而降低,随厚度的增加而增加。但目前对于高温环境下的研究主要集中在壁板,对于梁在高温环境下以及加人硬涂层之后的研究很少。研究发现通过合理制备硬涂层可以实现热稳定性15-17 ,所以研究高温环境下硬涂层对于叶片的减振效果是非常必要的。本文主要通过建立高温环境下发动机叶片-硬涂层阻尼结构动力学模型，分析温度、硬涂层材料以及尺寸参数对其动力学的影响，为涂层阻尼结构的合理设计以及叶片振动抑制提供理论支撑。1动力学模型建立以涡轮叶片为例进行高温环境下发动机叶片-硬涂层阻尼结构动力学建模与分析。在叶片基底表面制备一层与基底界面结合性良好的硬涂层材料，得到复合结构,如图1所示,其中:h,和h分别表示基底和硬涂层厚度；硬涂层基底界面与x轴距离用d表示；上表面红色箭头表示温度为T，的温度环境。叶片基底图1力学模型Fig.1Mechanical model使用复模量表示叶片基底、硬涂层和复合结构的弹性模量E,=E,(1+in.)E=E.(1+in)E=E(1+in)式中：i表示虚数单位；E。E和E分别为叶片基底、硬涂层和复合结构的复模量;E.E。E、n n n 为对应的弹性模量以及损耗因子。通过分析复合结构力矩、弯矩等平衡条件,列出平衡方程。通过平衡方程18 推导得到复合结构弹性模量以及损耗因子表达式为E=1+2(2h+3ht2h3,eh/.(1+ehi)(1+hi)3ehl_3+6hi+4h?+2eh+eht(n-n)+ns李双宝，李芳：高温环境下发动机叶片-硬涂层阻尼结构动力学建模与分析9.zdz=-EI 0dx式中:b为复合结构的宽度;h=hs+he为复合结构总bh3高度;I=表示截面二次矩。忽略剪切变形，将=129()代人式(5),得到弯矩M 与度的关系M=-EI 29(x,L)x2式中：9(x,t)为挠度;EI称为梁的抗弯刚度。在高温环境下，根据热应力理论19,热环境中梁在没有约束时的伸长量为l =T I式中:为线性热膨胀系数;T=T-T。表示温度变化,T为变化后温度,T。为变化前温度;l为复合结构的长度。坐标x处热应变为8xT=硬涂层EESh式中：Qxr为热应力;Nr为热轴力S=hb表示复合结构横截面面积。l表示热轴力造成的一个微端dx的hX-31-式中:e=E./E;hi=he/hso位移u相对x轴坐标的变化率等于应变8 x，产生沿x轴方向的应力为0,=Ee,=Eoudx在截面内积分计算合力矩，得到弯矩2M=b伸长量，即Al=/dAlr=1.8xdx=JoEhb由于该微端在x方向存在约束条件,则自由伸长量和由热轴力N造成的伸长量总和为0,代人式（7）、式(9)得到ATI=Nrl(1)Ehb则热轴力表达式2 0 为Nr=EhbT根据Euler-Bernoulli梁理论，以及忽略截面绕中性轴转动惯量和剪切变形，得到温度变化情况下梁在外部激励F作用下的微分方程为2EI9(a11+ps0(a.t)dx20 x2(2)N.0g=F2式中：p为密度;t为时间。简谐激振力会引起叶片强迫(3)振动,采用外部激励F表达式为F=Fosin(wt）。式(12)(4)(5)(6)(7)(8)NTdx=NilEhb(10)(11)at2(12)-32-与不考虑温度条件下微分方程区别主要体现在含有热轴力的参数项。约束边界条件为29(x,t)239(x,t)=0 x=lx=l考虑NiCrAIY目前应用广泛并且可以延缓结构寿命2 1,所以选用 NiCrAIY 来进行复合结构减振分析研究，叶片基底则选用Q235,其材料参数的取值如表1 所示2 。表1材料参数的取值Tab.1Parameter value of thematerial材料弹性模量/GPaQ235200.0NiCrAIY56.92动力学分析采用以上材料,化简式(12)可得9()+(100(g-N,9()=atpsFps根据振型叠加法，n阶阵型叠加分离变量后横向位移变量拥有表达式9(x,t)=20(x),()ni=1式中:和均为分离变量;表示只与x有关的函数，q表示只与t有关的函数。将式(16)代人式(15)可得:(x)g(t)+1-(Elg;(x)g,(t)-Nr9(x)g,(t)=F(17)pS式中（x）、（x）分别表示:的2 阶导和4阶导。4对式(17)进行化简，根据边界约束条件引人悬臂梁振动特征方程cos 入;cosh 入;+1=0式中：入，为悬臂梁特征值。从而得到悬臂梁振型p:(x)=cosh(4x)-cos(4cosh 入;+cos 入;sinh 入;+sin 入;根据正交性得到中国民航大学学报9;(x=0(13)x=0(14)密度/(kgm)泊松比7850.00.32.840.70.30 x40 x2x）(sinh(ix)-sin(ix)2024年2 月jix)dx=j=i0jix)dx=j=i式中,k;=六，和表示阶次。计算过将式(19)代入并使用三角函数与双曲函数泰勒展开进行化简得到(21)216j=i39824式中,令K=723化简得到Eh(1+jn)12i=1EhATZq(t)=1T.(x)e(a)da+phZi.(t).0:(x)g)(x)dx=F,(t)=1式中:F(t)06(15)Fioexp(iwt),其中,Fo为外部激励幅值,Fo为i阶激励幅值。当i=时,式(2 2)可以写成(16)Eh(1 in)k lq()-EhaATg(0)K+phl(t)=12F(t)化简得到phlq:(t)+F(t)根据方程特性可以得到固有频率表达式为4Eh34k,l-EhT KW;12Eh34-EhATK1121Iph化简得到方程(18)ji(t)+w0ig(t)+jm假设式(2 6)的解为(19)q:(t)=A;exp(i0)expi(wt-.)式中A为响应幅值。代人式(2 4)得到(20).0;(x)p(x)dx,通0ji8式(17)两边同乘以 s;()并9q:(t);(x)p;(x)dx-Fosin(wt)g;(x)dx=(23)Eh(1+in)k,1-EhaATK a(t)=12112Eh3ph:(t)=F(t)(22)(24)1phl4(25)(26)(27)第42 卷第1期李双宝，李芳：高温环境下发动机叶片-硬涂层阻尼结构动力学建模与分析-33-4Eh3-wA;+w,A;+jn12ph从而可以得到幅值和相位角的表达式分别为FoA;=十4Eh312ph;=tan*12W;之后可以根据所得公式进行动力学分析。3数值模拟3.1加入硬涂层固有频率作为结构最重要的特性,对于研究结构在受到外界激励产生运动时的动力学特性有着重要意义,另外,若振动能量降低,最明显的变化是响应幅值的降低。所以在研究系统减振时,研究幅频响应曲线的变化情况也是必不可少的。根据上述对方程化简得到的结果表达式,用Matlab对所得结果进行编程分析,只有叶片基底he=0时,不考虑温度影响(T=0),得到叶片基底幅频响应曲线如图2 所示。0-2-4-6/V-8-10-12-140500 1 000 1 500 2 000 2 500 3 000 3 500图2 叶片基底幅频响应曲线Fig.2Amplitude frequency response of the blade base从图2 可知,叶片基底的1阶固有频率为16 5.19 Hz，2阶固有频率为10 35.4Hz。加人0.1mm硬涂层后，得到复合结构振动响应如图3所示。从图3中可得：加入硬涂层后各阶固有频率均有所增加，在固有频率没有明显变化的同时整体振动响应幅值呈现减小的趋势。这主要是因为硬涂层内部存在细小不规则颗粒,在振动过程中颗粒之间相互摩擦消耗能量,使得动能降低,所以振动响应幅值降低。复合结构与叶片基底振动趋势没有太大改变，主要体现在固有频率增加以及振动响应幅值降低。其中1阶固有频率变化范围最小，2 阶变化幅度逐渐增大。具体变化如图4所示。0A(28)：(29)42Eh12phf/Hz一叶片基底-2-加人硬涂层-4-6-8-10-12(30)-1405001 0001 50020002.50030003500f/Hz图3口叶片基底与复合结构幅频响应曲线Fig.3 Amplitude frequency response of blade base andcomposite structure-1.5一叶片基底-2.0-加人硬涂层-2.5-3.0-3.5-4.0-4.5-5.0-5.5-6.0-6.5150155160165170175180f/Hz(a)1阶-5.0-5.5-6.0-6.5-7.07.5-8.0-8.5-9.0-9.5100010201 0401 06010801100f/Hz(b)2阶图41阶和2 阶响应曲线Fig.4Response curves of first order and second order通过对比图4的(a)和(b)可以发现,固有频率增大和响应幅值减小的幅度都随阶数增大而增大。具体表现为1阶固有频率从16 5.19 Hz变化到16 7.48 Hz，变化范围在2 Hz左右;2 阶固有频率从10 35.4Hz变化到10 49.7 Hz,变化范围在14Hz左右。3.2改变材料参数在3.1节中,通过对比加人硬涂层前后结构参数、各阶固有频率和振动响应幅值变化规律，得到加入硬涂层后具有增大固有频率以及减小响应幅值的效果，主要是由于硬涂层减振作用。本节具体讨论硬涂层弹性模量、损耗因子以及硬涂层厚度等材料参数对于结果的影响。1一叶片基底-加人硬涂层34结构参数主要分析复模量变化情况,其中实部表示弹性模量，虚部表示损耗模量，损耗因子为损耗模量与弹性模量之比。首先得到复合结构弹性模量随硬涂层弹性模量变化情况如图5所示。1.871.861.85(eduo1)/a1.841.831.821.811.805.56.06.57.07.58.08.59.09.510.0Ec/(10Pa)图5E随Ec变化趋势Fig.5Changes of EwithEc之后得到复合结构弹性模量随硬涂层厚度变化情况如图6 所示。1.91.81.7(ed mot)/a1.61.51.41.31.21.112345678910hc/(10-m)图6E随hc变化趋势Fig.6Changes of Ewith hc由图5和图6 可知，复合结构弹性模量随硬涂层弹性模量增加而增加，随硬涂层厚度增加而减小，但整体复合结构弹性模量小于叶片基底。图7 描述了复合结构损耗因子随硬涂层损耗因子的变化情况。4.03.83.63.43.2(-01)3.02.82.62.42.22.00.005Fig.7图8 描述了硬涂层厚度对复合结构损耗因子的影响。中国民航大学学报5.55.04.5(01)/l4.03.53.02.52.012345678910hc/(10-m)图8 n随hc变化趋势Fig.8Changes of n with hc由图7 和图8 可以得到复合结构损耗因子随硬涂层损耗因子和厚度增加而增加。图9 描述了不同硬涂层损耗因子对复合结构幅频响应曲线的影响。通过图9 可以发现，硬涂层损耗因子对复合结构固有频率没有影响，但振动响应幅值随硬涂层损耗因子增大而逐渐减小。根据式(2 9)可知，硬涂层损耗因子增大使得响应幅值减小与理论结果一致。-1.5-2.02.5-3.0uu/-3.5-4.0-4.5-5.0-5.5-6.0160图9 复合结构在不同nc下的响应Fig.9 The response of the composite structure under different nc图10 描绘了硬涂层弹性模量对结果的影响。虽然加入硬涂层后复合结构弹性模量降低，但质量等也发生改变,最后得到固有频率升高的结果。与损耗因子不同的是，改变弹性模量不仅响应幅值减小，固有-1.5-2.0-2.5-3.0uu/-3.5-4.0-4.50.0150.025nc图7 m随nc变化趋势Changes of n with nc2024年2 月0-0.00990.01990.0399160170f/HzEc/GPa056.960.070.00.0350.045175-5.0-5.5-6.0160图10 复合结构在不同Ec下的响应Fig.10 The response of the composite structure under different Ec165f/Hz170175第42 卷第1期频率也有所增加。可以得出复合结构固有频率只与硬涂层弹性模量有关，而响应幅值不只与弹性模量有关还受到损耗因子影响。图11展示了硬涂层厚度对于幅频响应曲线的影响，可以看出，复合结构固有频率随硬涂层厚度增加而增大，响应幅值随硬涂层厚度增加而减小。-1.5-2.0-2.5-3.0-3.5wu-4.0-4.5-5.0-5.5-6.0-6.5160图11复合结构在不同hc下的响应Fig.l1 The response of the composite structure under different hc由于硬涂层厚度足够小，所以对于复合结构质量等影响可以忽略不计,增加硬涂层厚度可以近似看作弹性模量与损耗因子的改变。3.3加入温度影响由于涡轮叶片所处高温环境主要是在10 0 0 左右,因此主要研究叶片在0 10 0 0 温度变化范围内的变化趋势，以及温度对振动的具体影响以及机理。加人温度影响因素,研究温度对复合结构振动响应结果的影响。不同温度下叶片基底振动响应的结果如图12 所示。0-12上-3-4二-5-6-7-8-950图12 叶片基底在不同温度T下的响应Fig.12The response of the blade base under different T由图12 可以看出，随着温度的升高，固有频率逐渐减小,振动响应幅值增大。很明显,结构的动态响应随着温度的增加整体地向低频范围移动,并且在不同的温度环境中响应的整体特征保持不变。这种变化趋势主要由热环境中自然振动的变化决定。固有频率随着温度的升高而降低,为研究其具体规律,需要判断其在0 10 0 0 范围内具体变化情况。对李双宝，李芳：高温环境下发动机叶片-硬涂层阻尼结构动力学建模与分析hc/mm00.10.2164168J/Hz100150f/Hz-35-固有频率和温度表达式进行分析得到图13的曲线。180160140120100806040200图13不同温度T下1阶固有频率Fig.13First order natural frequency under different T172176T/一0-200500-950-1200200250500180从图13可看出，在温度0 2 0 0 0 范围内1阶固有频率一直呈现减小趋势，而研究环境在10 0 0 左右，所以可以不考虑其他情况。接下来分析硬涂层以及温度对于叶片基底振动的联合影响。与常温下进行对比,发现不论是叶片基底还是复合结构，固有频率均随温度升高呈减小的趋势。将叶片基底、复合结构以及温度影响下复合结构共同进行对比，如图14所示。-1.5-2.0-2.5-3.0-3.5二-4.0A-4.5-5.0-5.5-6.0160图14叶片基底与复合结构在不同温度T下的响应Fig.14The response of the blade base and composite structure underdifferent T从图14中可以发现，加入硬涂层后固有频率增加，但温度升高固有频率降低，在某一特定温度下复合结构固有频率会低于叶片基底固有频率。4ANSYS有限元分析首先建立叶片基底和复合结构有限元模型，并且划分合理的网格之后进行模态分析，从而得到各阶模态频率及模态振型，之后进行瞬态温度分析耦合力学，模态分析得到其在高温环境下动力学响应。1000T/一叶片基底T=20叶片基底T=50叶片基底复合结构-T=20复合结构T=50复合结构上165f/Hz15001702 000175-36-4.1叶片基底参考发动机叶片及材料振动疲劳试验方法2 31准则，建立有限元模型结构，尺寸为长10 0 mm，宽2 0 mm，叶片基底厚2 mm,硬涂层厚0.1mm。其中叶片基底采用Q235，硬涂层阻尼材料采用NiCrAIY，设置材料参数。划分网格如图15所示，得到节点数为37 56 8。设置左端固定、右端自由的约束边界条件，之后进行模态分析可以得到各阶固有频率和振型，图16振型幅值357.82Max318.06278.30238.54198.79159.0379.515119.2739.757oMin振型幅值450.49Max400.43350.38300.33250.27200.22150.16100.1150.054oMin振型幅值359.95Max319.95279.96239.97199.97159.98119.9879.98939.997oMin由图16 可以看出，结构第1、2、4、5阶表现为弯曲振动，第3、6 阶表现为扭转振动，弯曲振动最大变形位于叶片最右端，最小变形位于叶片的最左端，扭转振动的最大变形位于叶片的两侧，最小变形位于叶片中间。振型幅值354.62Max315:22275.82236.42197.01157.61118.2178.805oMin39.403振型幅值351.89Max312.79273.69234.59195.50156.40117.3078.19839.099LoMin振型幅值356.88Max317:22277.57237.92198.26158.61118.9679.306oMin39.653中国民航大学学报展示了叶片基底各阶振型图。振型幅值357.97Max318.20278.42238.65198.87159.100.0020.0010.0030.00(a)1阶振型0.0020.0030.0010.00(c)3阶振型0.0020.0010.0030.00(e)5阶振型0.00015.0007.50022.500(a)1 阶振型30.000(mm)0.00015.0007.50022.500(c)3阶振型0.00015.0007.50022.500(e）5阶振型Fig.17 Vibration mode of composite structure at each order2024年2 月015.00030.000(mm)7.50022.500图15划分网格图Fig.15Gridding diagram40.00(mm)119.3279.54939.775oMin振型幅值354.03Max314.69275.35236.02196.68157.35118.0140.00(mm)78.67339.336oMin振型幅值469.61Max417.44365.26313.08260.90208.72156.5440.00(mm)104.3652.179oMin图16 叶片基底各阶振型图Fig.16Vibration mode of blade base at each order4.2加入硬涂层在叶片基底的基础上加人0.1mm硬涂层,硬涂层损耗因子通过设置阻尼来体现，选择Material DependentDamping中的DampingRatio来进行设置。加人硬涂层之后,固有频率与振动幅值皆有所变化,结果如图17 所示。振型幅值354.83Max315.40275.98236.55197.13157.70118.2830.000(mm)78.850Min39.425振型幅值447.14Max397.46347.78298.41198.73149.0599.36549.683oMin39.757振型幅值465.72Max413.97362.23310.48258.73206.9930.000(mm)155.24103.4951.747oMin图17 复合结构各阶振型图Z0.0020.0010.0030.00(b)2阶振型0.0020.0010.0030.00(d)4阶振型0.0020.0010.0030.00(f)6阶振型0.00015.0007.50022.500(b)2阶振型0.00015.0007.50022.500(d)4阶振型0.00015.0007.50022.500(f)6 阶振型40.00(mm)40.00(mm)40.00(mm)30.000(mm)30.000(mm)30.000(mm)第42 卷第1期从图16 和图17 可以看出，加入硬涂层后3阶和4阶振型发生改变。对比叶片基底与加人硬涂层后复合结构固有频率、响应幅值变化分别如表2 和表3所示。表2 叶片基底与复合结构各阶固有频率对比Tab.2 Comparison of natural frequency of blade base andcomposite structure at each order固有频率/Hz阶次叶片基底1阶164.972阶1 031.023阶1 586.834阶1 586.945阶2.883.406阶4 849.66表3叶片基底与复合结构各阶响应幅值对比Tab.3Comparison of response amplitude of blade base andcomposite structure at each order阶次响应幅值/mm叶片基底复合结构1阶357.822阶357.973阶450.494阶354.035阶359.956阶469.61分析图16 和图17 及表2 和表3结果，可发现结果与Matlab一致,且加人硬涂层后3阶、4阶固有频率以及振型发生对调。对其进行谐响应分析，如图18 所示。9.248 810-23.053 710-21.008 210-2uu3.328 910-31.099 110-33.628 910-1.198 210-43.956 010-51.306 1x10-040 80120160200 240280f/Hz(a)叶片基底8.050210-33.611 410-31.620 110-37.268 010-4=3.2605x104-1.462 710-46.561 810-52.943710-51.320 6x10-5L040 80120160200 240 280J/Hz(b)复合结构图18 叶片基底与复合结构各阶定向变形对比Fig.18 Directional deformation of blade base and compositestructure at eachorder根据图18 得到加人硬涂层后振动响应幅值减小，更好地验证了上述结果。李双宝，李芳：高温环境下发动机叶片-硬涂层阻尼结构动力学建模与分析复合结构167.24.1 045.151 584.201 607.572.922.504912.40354.62354.83351.89447.14356.88465.72-37-4.3加入温度影响分别对叶片基底和复合结构进行瞬态温度分析再耦合力学进行分析。对材料导热系数、热膨胀系数和比热容等进行设置，对结构施加如图19 所示温度变化过程。1000875750625500375250125250图19 温度变化过程Fig.19Change process of temperature瞬态温度分析后得到每个时刻的温度分布，图2 0为5s时温度分布情况，导入温度、耦合力学与模态分析,分析结果。温度/512.50Max509.92506.75503.87500.99498:13493:24492.36489.49486.6iMin温度/512-50Max509:12505:25502.37499.00495:62492.25488.87485.50482.12MinFig.20 Temperature distribution at 5 s根据图2 0 得到叶片基底与复合结构5s时温度分布情况近似均在50 0 附近，耦合力学与模态分析得到模态分析结果如表4和表5所示。表45s时固有频率对比Tab.4Comparison of natural frequency at 5 s阶次1阶2阶3阶4阶5阶6阶根据表4和表5得到高温环境下硬涂层对于复合结构依然具有增大固有频率,降低响应幅值的作用。160t/sM0.00(a)叶片基底0.0040.00(mm)20.00(b)复合结构图2 0 5s时温度分布固有频率/Hz叶片基底复合结构165.14167.511 031.301.046.501587.101584.801 590.101 615.202.883.702.928.204.857.904934.5017030.00(mm)15.0040-38-Tab.5Comparison of response amplitude at 5 s阶次1阶2阶3阶4阶5阶6阶之后进行谐响应分析，如图2 1所示。7.684.703.463 101.560 700.703320.316960.142.840.064.370.029 010.013 07040 80120160200240280f/Hz(a)叶片基底3.951 101.937 700.950280.466 040.228550.112.090.054970.026960.0132204080120160200240280f/Hz(b)复合结构图2 15s时1阶总变形Fig.21 First order total deformation at 5 s从图2 1可以看出，固有频率增大，响应幅值降低，更好地验证了上述结论。当温度再升高，10 s时温度分布情况如图2 2 所示。温度/964.9Max959.62962.26956.98954.34M951.70946.42949.06941.13Min943.77温度/964.9Max961.79958.69955.58952.47946.26949.37943.15940.05936.94MinFig.22TTemperature distribution at 10 s中国民航大学学报表55s时响应幅值对比从图2 2 可以看出，10 s时温度升高到高于9 0 0，叶片基底与复合结构温度分布情况无明显区别，研究响应幅值/mm900下硬涂层作用是否依然有效。耦合力学、模态分叶片基底复合结构357.86354.71357.94354.94354.35351.36451.17447.22359.96357.16469.86465.88Empotedfos-Tpeati0.0015.00(a)叶片基底0.0040.00(mm)20.00(b)复合结构图2 2 10 s时温度分布2024年2 月析得到高温下叶片基底和复合结构固有频率以及响应幅值对比如表6 和表7 所示。表6 10 s时固有频率对比Tab.6Comparison of natural frequency at 10 s阶次1阶2阶3阶4阶5阶6阶Tab.7Comparison of response amplitude at 10 s阶次1阶2阶3阶4阶5阶6阶对比表6 和表7 发现，与表4和表5得到结果一样，硬涂层依然具有增大固有频率，降低响应幅值的作用。对比表4和表6 以及表5和表7 发现,随着温度的升高，叶片基底和复合结构均表现出固有频率降低、响应幅值增大的效果，与硬涂层单独作用效果正好相反。之后进行谐响应分析，如图2 3所示。5.997702.788 201.296 200.602.570.280120.130 220.06 0540.02.8140.01 3080408012016020024028030.00(mm)f/HzZ(a)叶片基底5.011 202.385901.136 00tu0.540 860.25752MA0.122610.05 8380.02.7790.01 32304080120160200 240280f/Hz(b)复合结构图2 310 s时1阶总变形Fig.23First order total deformation at 10 s固有频率/Hz叶片基底复合结构165.00167.331 031.101 046.101 587.001 584.601587.301 611.702.883.402.927.404.850.904.925.80表7 10 s时响应幅值对比响应幅值/mm叶片基底复合结构357.91354.75357.91354.91354.20351.29451.46447.58359.96357.16470.09466.11第42 卷第1期从图2 3看出，复合结构与叶片基底相比，固有频率增大且响应幅值减小，对比图2 1和图2 3发现，温度升高，固有频率降低，响应幅值增大，更好验证了上述结论。5结语本文主要研究高温下发动机叶片-硬涂层阻尼结构动力学建模与分析。研究结果表明，硬涂层复模量对于复合结构固有频率和振动响应幅值有影响：其中，复合结构弹性模量随硬涂层弹性模量增加而增加，又因质量等改变从而增大固有频率；而硬涂层损耗因子通过增大复合结构损耗因子减小振动响应幅值。加人温度影响后发现，随着温度的升高，叶片基底和复合结构均表现出固有频率降低、响应幅值增大的效果，与硬涂层单独作用效果正好相反。若对其合理利用，两者共同作用下可以使得复合结构避开共振区且降低响应幅值，达到优化效果。参考文献：1】陈辉.硬涂层阻尼叶片振动疲劳特性研究D大连：大连理工大学,2 0 19.2 PATSIAS S,SAXTON C,SHIPTON M.Hard damping coatings:an ex-perimental procedure for extraction of damping characteristics and modu-lus of elasticityJ.Materials Science and Engineering:A,Structural Ma-terials,2004,370(1/2):412-416.3 TORVIK P J,HENDERSON J P.Influence of glass content on dampingproperties of plasma-sprayed mixtures of zirconia and glassJ.Journalof Materials Engineering and Performance,2012,21(7):1405-1415.4 SUN J,KA