高温环境下二维正交编织C∕SiC复合材料壁板振动模态演化.pdf 12页

航空学报ActaAeronauticaetAstronauticaSinicaJul252017VoI．38No7SSN1000—6893ON11-1929／Vhttp：／／hkxb．buaa．edu．cahkxb@buaaeduca高温环境下二维正交编织振动模态演化贺尔铭*，陈兵，曹存显西北工业大学航空学院，西安710072C／SiC复合材料壁板摘要：为了研究高温环境下二维正交编织C／SiC复合材料壁板的固有振动特性随温度的变化规律，进一步揭示模态跃迁和丢失现象。首先建立了二维正交编织复合材料的细观非均匀有限元模型，基于细观模型采用体积平均法计算得到了均匀化后材料的模量、热膨胀系数和热传导系数等宏观性能参数。在此基础上，分别研究了均匀温度和线性非均匀温度载荷下，四边简支复合材料板在屈曲前后固有频率及模态振型随温度的变化规律。研究结果表明：在均匀温度场下，如果仅考虑热载荷对材料模量的影响时，随着温度的升高，各阶固有频率逐渐降低，但降低的幅度不大；仅考虑热应力对结构刚度影响时，在临界屈曲温度后固有频率反而上升，并且随着温度的升高出现了频率交错，其对应振型发生了跃迁现象；当同时考虑热物性和热应力影响时，各阶频率值变化趋势与仅考虑热应力时类似，只是对应温度点的频率值有所下降。而在线性非均匀温度场下，通过对比不同工况下的模态振型发现，随着温度的升高某些模态振型还会出现丢失，并且这种丢失是从低阶到高阶依次发生的。关键词：二维C／SiC；复合材料壁板；热模态；固有频率；模态跃迁和丢失现象中图分类号：V214．8文献标识码：A文章编号：1000—6893(2017)07—220553—12高超声速飞行器在某些飞行条件下会造成严重的气动加热问题，气动加热通过附面层传到结构内部导致结构温度急剧升高，使飞行器表面壁板结构面临着较为严峻的高温环境。而二维正交编织C／SiC复合材料具有高比强度、高比刚度、耐高温等优异性能，同时具有平行碳布面内力学性能相同的优点，特别适合于制备典型的热结构件，因而得到了广泛关注。高温对结构的影响主要表现在2个方面：一方面会改变材料的机械性能使材料的弹性模量降低；另一方面温度梯度导致受约束结构内部产生热应力，这些影响都会导致结构刚度发生变化，从而直接影响到结构的固有频率，因此在进行高速飞行器结构热振动分析时必须考虑温度效应的影响。近年来，国内外许多学者开始关注壁板结构的热振动特性研究。早在1989年Spain等[1]就关注到温度升高会使结构固有频率降低的现象；随后Shukla等[2]研究了温度对复合材料曲板前10阶固有频率的影响规律，结果发现：随着温度升高，前5阶频率缓慢上升而后5阶频率却缓慢下降，并且在某一温度时各阶频率大小顺序发生了改变，即出现“频率跃迁”现象；夏巍和杨智春[3]研究了复合材料平板的热模态特性，发现屈曲前结构各阶频率都会出现下降，而且对于两对边简支的正方形复合材料板频率下降得更严重，但是其研究并未考虑到结构非线性引起的刚度非线性收稿日期：2016—06—22；退修日期：2016—10—19；录用日期：2017—03—07：网络出版时间：2017-04—0113：32网络出版地址：WWW．cnki．net／kcms／detail／111929V．201704011332．006．html基金项目：国家自然科学基金(51675426)；航天科技创新基金(2014KC010043)*通讯作者E-mail：heerming@nwpu．edu．ca嘲用格武：贺尔铭，陈兵．营存显．高温环境下二维正交编织c／Sic复合材料壁板振动模态演化!如，航空学报．2017，38(7)：220553jHEEM．CHENB．CAOcX．Vibrationmodeevolutionof2Dwovenc／s|ccompositepanels{nhotenvironment!JIActaAeronu．aticaetAstronauticaSinica，2017，38t7)：220553．220553．1 航空学报项，因此只适合结构发生热屈曲前的固有振动特性分析；刘芹等¨j对线性温度场作用下的层合板进行了非线性热振动分析，得到了基频随温度的变化关系。对于其他高阶固有频率随温度的变化并未进行深入研究；Singh和Panda一刀研究了考虑高阶几何非线性项时，双曲复合材料板的自由振动问题，杨莉等¨o分析了不同因素对圆柱曲板振型的影响，研究发现温度场对功能梯度材料的频率影响不大。在实验研究方面，许多学者对不同材料板的固有频率随温度变化做了大量研究，其实验结果都表明，随着温度的升高结构的前几阶固有频率出现了下降u1]。史晓鸣和杨炳渊一”-通过对瞬态加热环境下变厚度板的热模态分析，发现温度使结构的固有频率出现下降，但持续加热时固有频率又会出现回升；黄世勇和王智勇。1L研究了根部固支的变厚度铝合金板的热模态，发现随着温度升高结构的固有频率出现了下降，但其并未对振型演化规律进行研究。对于模态振型的研究，Cheng等。121通过实验分别研究了温度载荷对平板和加筋板模态振型的影响，发现随着温度的升高无论是平板还是加筋板模态振型都出现了交错现象。上述研究虽然给出了结构固有频率与温度变化的关系，但是并未全面考虑不同参数对复合材料板固有振动特性的影响规律，而且对于频率演化的同时振型出现跃迁、丢失等这些现象并未从理论上进行研究，只有少数文献指出在实验中发现了这些现象，对于出现振型跃迁、丢失的机理也没有进行更多的分析。目前对于复合材料热振动特性的研究大多数都是基于层合板或常见的金属板，对于二维正交编织C／SiC复合材料板的很少。本文对高温环境下二维正交编织C／SiC复合材料板的热振动特性进行了深入分析，通过细观建模得到了材料宏观等效参数并与实验值进行了对比发现误差较小，在此基础上，分别研究了均匀温度场和非均匀温度场下．四边简支复合材料壁板结构的同有频率及模态振型的演化规律。通过对温度载荷作用下复合材料板结构进行热模态分析，得到了结构振型频率随温度的变化曲线，并利用有限元仿真再现了二维正交编织材料板热振动实验中出现的模态振型跃迁及丢失现象，还运用摄动法对模态振型跃迁的机理进行了解释。1复合材料板的细观建模及热屈曲分析二维正交编织复合材料内部结构复杂，其模量等性能参数随温度变化的试验数据缺乏。因此通过精确建立编织材料的细观模型，运用有限元方法模拟预测材料的宏观等效性能参数就十分必要了。本文从细观模型出发，运用均匀化假设预测了等效后的材料性能参数。在此基础上，完成了壁板结构的热屈曲分析及热模态分析。本文所选用的复合材料细观模型是参考已有文献u副中对二维正交编织c／sic复合材料的建模方法，通过在扫描电镜图片中取点、拟合纤维曲线从而建立如图1所示的精确代表体积单元(RVE)模型，运用结构有限元软件PATRAN中的布尔运算得到基体模型。对于真实材料中的空隙模拟，这里是通过在有限元建模时分别在纤维束、基体以及纤维束与基体之间预制裂纹来实现的，最终得到如图2所示的有限元细观计算模型(其中纤维体积含量为40％左右，基体为45％左右，空隙为15％左右)。在对材料宏观参数进行预测时，目前使用最多的方法为通过均匀化假设将其等效为一个宏观220553．21RVE几何模型R、7Egeometricmodel图2有限元细观计算模型Fig．2FiniteelementmicroscopiccalculationI"110( 航空学报均匀材料，并施加周期性边界条件从而得到宏观材料参数，随后将预测值与实验值进行对比，验证所建立细观模型的有效性口。15]。本文在进行材料参数预测时也采用类似的方法，不同的是本文所建的有限元模型是基于真实电镜扫描照片，对纤维束的走势也不是简单地用折线代替。在进行材料模量预测时，选取组分材料的弹性模量随温度T的变化如图3和图4所示[161(其中E和E。分别为某温度下和常温下材料的弹性模量)，采用前面所述的细观力学方法，计算得到了不同单胞形式的等效模量参数(如表1所示，其中E，。为沿着纤维方向的模量值，E3。为与纤维垂直方向的模量值，其他参数方向的定义与此相同，v为泊松比，G为剪切模量)，对比发现预制孔隙后预测值与实验数据吻合良好。同时改变组分模量参数，得到了不同温度下二维正交编织C／SiC复合材料的参数，与实验值(如表2所示)相比误差较小。使用表1中的材料属性建立了四边简支的复合材料壁板结构，其几何尺寸为：0．61TI×0．61TI×0．003m，并采用Shell单元建立如图5所示的有限元模型。通过有限元计算得到薄板在均匀温度载荷下的临界屈曲温度为34℃。四边简支板的临界屈曲温度也可由式(1)计算。【o裔刀℃图4SiC基体弹性模量随温度的变化曲线Fig．4CurveofelasticmodulusofSiCmatrixvstemperature表1复合材料参数的预测值与实验值Table1Predictiveandexperimentalvaluesofcompositematerialparameters表2不同温度下复合材料参数的实验值Table2ExperimentalvalueofcompositematerialparametersatdifferenttemperaturesEll／GPaE33／GPaG23／GPav1212080．488．165．843．60．250．35220553．32198∞蚰“驼n牡7O0O72∞m∞卯 航空学报式中：a、b和h分别为薄板的长、宽和厚度；a为热膨胀系数；埘、挖与模态振型分布有关；D为与弹性模量、剪切模量、泊松比和板厚度有关的弯曲刚度。由式(1)计算得到临界屈曲温度为34．31℃，与有限元计算结果十分接近。在后续热模态计算时由于温度范围变化大，涉及到屈曲前后的整个过程，这样屈曲后的热模态分析需重新建模。这里提取特征值屈曲分析得到的1阶屈曲模态，将其以l％的比例更新到结构中，通过计算得到该结构的屈曲平衡路径如图6所示。图5壁板有限元模型Fig．5FiniteelementmodelofpanelBucklingcoefficient图6简支板屈曲平衡路径由图6可见，屈曲前结构的法向位移较小，屈曲后法向位移突然增加，但是增加的幅度越来越缓．在临界屈曲载荷后由于结构刚度有所增大，所以法向位移增加变缓。通过对复合材料壁板的屈曲分析，可以得到不同温度载荷下的热应力和应变场，在后面的热振动模态分析时可以直接将得到的应力应变场作为载荷引入进来。2壁板结构热模态分析温度载荷对结构模态的影响主要表现在2个方面：①温度的改变会使材料的弹性模量发生变化，随着温度的增加材料的弹性模量会下降口7。1叫；②温度梯度的出现会使结构内部产生变形从而产生热应力。材料特性的下降和热应力都会导致结构刚度出现变化，从而直接影响到结构的固有频率以及振型分布，因此，在后面的热模态分析时本文主要也是从这两方面出发的。Malekzadeh等一!‘。利用Hamilton原理推导出了系统在热环境下的运动方程，本文基于能量原理和变分方法的Hamilton原理，建立复合材料壁板在热环境下的振动方程为[：6：][汁憾kbZ1卜陈》r足N。一kjTL0(2)式中：k，为剪切刚度；正。为剪切矫正因子；肌n和k。分别为平板结构的质量矩阵和刚度矩阵；k。为温度引起的刚度变化；‰，为耦合刚度项；w。，为结构的位移场；w。为平板的位移场；k№为结构的面内位移在弯曲自由度上的附加刚度；k。，为温度应变在板弯曲自由度上的附加刚度；P。。r和P。r分别为温度载荷引起的附加力和附加力矩。对于正交各向异性的复合材料平板结构，刚度矩阵中的耦合刚度项k。，一0，并且有Pm一0，则式(2)可以变为13Llcklinge[1uilibⅢ1nⅥ矾h。fsimplyⅧpp。m—Jkmw,,一pm△TIm。lm^品一+(七，+忌2)w6—0220553．4(3)一1llJ忱‰r1．、lll●●●。、、，●●●J]lllJ0O 航空学报式中：七，一k。+忌；k，为不计及温度影响的复合材料壁板结构的刚度矩阵；k。一kNm+k。r为由于温度载荷引起的刚度矩阵的修正项。在实际的求解中，首先求出式(3)中的第1个静力学方程，得到结构的位移场W。，并将结构变形施加到结构上，即为动力学分析的刚度矩阵引入了修正项足。，在此基础上就可以对某一特定温度场下复合材料壁板的振动特性进行分析，这样就在动力学分析中考虑了热效应的影响。2．1均匀温度载荷下壁板结构热模态分析在均匀温度场下进行模态分析时，本文从3个方面综合考虑了热物性和热应力对结构热模态的影响：①仅考虑热载荷引起的材料弹性模量变化时，结构固有频率以及模态振型的变化规律；②仅考虑由热载荷引起的热应力时，结构固有频率以及模态振型的变化规律；③综合考虑两种因素时，结构固有频率以及模态振型的变化规律。在进行热模态分析时，选取参考温度为T一20℃，选取的结构为进行屈曲分析后更新过的有限元模型，采用有限元分析软件ANSYS，通过定义材料属性表来实现材料模量等随温度的变化，进行热分析得到温度分布、热应力等数值，随后运用该软件的模态分析模块进行热模态分析。只考虑温度载荷对材料参数的影响时，计算得到前10阶频率．厂随温度T的变化如图7所墨示，其中部分阶次曲线完全重合。可以看出，各阶频率随着温度的升高而降低，但没有出现频率交错现象，而且频率降低幅度也比较小，当温度升高到1000℃时1阶频率也只下降了9．3Hz。可以发现，在1000℃附近频率随着温度的增加反而增加，从表2中数据可看出这是由于该温度范围内材料模量随温度升高而增加导致的。仅考虑温度载荷对模量影响时得到的频率随温度变化不是很明显，为了区别不同因素对结构固有频率变化规律的影响，仅考虑热应力影响时选取弹性模量为一常数。计算时仅考虑热载荷引起的热应力对结构模态变化规律的影响。图8为仅考虑热应力时前10阶频率随温度的变化曲线。从图中可以看出在临界屈曲温度前，各阶固有频率均随着温度的升高而降低，且幅度比较大；达到临界屈曲温度时，复合材料板趋于失稳，第1阶频率接近于零；跨过临界屈曲温度后，前6阶频率开始升高且在临界屈曲温度附近增幅较大，同时随着温度的升高出现了频率交错现象，且某些相同阶次的频率值不再相同。图9～图11给出了不同温度下前6阶的模态振型，对比图9和图10、图10和图11中的振型演化顺序发现，临界屈曲温度后结构模态发生了跃迁现象(如图12所示)。在247℃时，第1阶与第2阶振型次序发生了交换，从图8中也可以看出在403℃附近也出现了一次跃迁。当温度达到图7仅考虑热物性时频率随温度的变化曲线Fig．7Frequency—temperaturecurveswhenonlyconsideringeffectoftemperatureonelasticmodulus220553．5 航空学报1000℃时，第1阶振型和第2阶振型又出现了交错，这与图8是吻合的。这些跃迁现象表明：在温度载荷作用下，复合材料板的模态振型会出现交错，也就是说温度载荷会使高阶模态振型提前出现。这里还给出了1100℃时前6阶模态振型(如图13所示)，从图中可以看出此时振型变化已经十分复杂。刀℃图8仅考虑热应力时频率随温度的变化曲线Fig．8Frequencytemperaturecurveswhenonlyconsideringthermalstressrr1●■L—r1●●Ljrr]r11●LJr1●LJr1●LJL—r1●●，L—_一r—1●●●L。jr1●■◆LJ(a)，=421I-Iz(b)f2=1336Hz(c)正21364Hz(d)f422255Hz(e)f5=2821l-lz(D‘22867Hz图929℃时前6阶模态振型图(仅考虑热应力)Fig．9Firstsixvibrationmodeshapesat29℃(whenonlyconsideringthermalstress)r1●■L—r]L—rr1r11●L—(a)f,=1595Hz(b)正=l597Hz(c)Z=1788Hz(d)‘22226Hz(e)‘22595Hz(f)‘22788Hz图10247℃时前6阶模态振型图(仅考虑热应力)Fig．10Firstsixvibrationmodeshapesat247℃(whenonlyconsideringthermalstress)r11●■L。Jr1●■一Jr1r1●LJr1I-．：一卜一r●11(a)f；=1895Hz(b)f,=2133Hz(c)六22134Hz(d)‘22392Hz(e)六22853Hz(f)‘。3323Hz图11403℃时前6阶模态振型图(仅考虑热应力)Fig．11Firstsixvibrationmodeshapesat403℃(whenonlyconsideringthermalstress)220553—6 航空学报”。●●-1t。⋯，一i。——爿、|k!—～⋯一———一t±———————i。：I图12模态跃迁现象Fig．12Modejumpingphenomenonr●]L一一0"一(d)Z2354．5I-lz(e)fs=383．8Hz(f)f6=,427．1Hz图131100oC时前6阶模态振型图(仅考虑热应力)Fig．13Firstsixvibrationmodeshapesat100。C(whenonlyconsideringthermalstress)从式(1)中可以看出，四边简支壁板结构的临界屈曲温度不仅与材料的特性、尺寸等参数有关，还与模态振型分布有关。笔者在Nastran非线性模块下进行了热模态分析，得到了不同阶次频率对应的临界屈曲温度，但是由于该情况在工程实际中并未得到广泛认可，因此此处不再进行深入研究，仅给出部分阶次频率随温度的变化曲线，如图14所示。综合考虑温度对材料参数和热应力的影响时，计算得到频率随温度的变化曲线如图15所示，与仅考虑热应力时的结果(图8所示)对比可以发现，当温度较低时模量的变化对频率的影响比较小，但是在温度较高时由于模量的降低导致频率出现了明显的变化，但也只是在各个温度点Temperature．／℃图14考虑后屈曲时各阶频率随温度变化曲线Fig．14Frequencytemperaturecurveswhenconsideringpostbuckling220553—7 图15同时考虑热应力和热物性时频率随温度的变化曲线Fig．15Frequency—temperaturecurveswhenconsideringthermalstressandeffectoftemperatureonelasticmodulus处频率值出现了下降，模量变化对结构整体模态振型变化趋势的影响不大。2．2非均匀温度场下壁板结构热模态分析在对非均匀温度场下的壁板结构进行热模态分析时选取的参考温度为T一20℃，首先在板的两端施加表3所示的不同温度边界条件，进行热分析得到温度分布，随后将得到的温度场以载荷的形式施加在结构上，再对整体结构进行静力学分析得到应力场，再将得到的热应力场作为初始边界条件施加在结构上进行模态分析。FirstfoorvibrationmodeshapesinCase3FirstfourvibrationmodeshapesinCase5FirstfourvibrationmodeshapesinCase7计算得到不同工况下的模态振型如图16所示，表3非均匀的线性温度载荷工况Table3Non—uniformlineartemperatureloads图16不同温度下的模态振型丢失Fig．16Vibrationmodeshapelosingunderdifferentthermalconditions220553—8 航空学报对比图中振型发现：随着工一Z边温度的增加，模态振型出现了丢失现象，工况3时1阶模态振型不再是1个半波的振型，变为2个半波的振型，工况5时1阶振型已经变为4个半波的振型。随着边界温度的增加，振型丢失也越来越严重。但是值得注意的是，这种丢失不是杂乱的，而是按照从低阶到高阶的顺序依次丢失的。在线性非均匀温度场下，随着温度的升高振型丢失越来越严重。结构处于较高的基础温度环境时，结构的模态振型出现了更复杂的变化。这里给出了2种不同的工况(如表4所示)，计算得到了这2种工况下的模态振型如图17和图18所示，由于此时模态振型变化复杂，这里仅给出了前4阶振型。对比图17和图18可以发现，当整体环境温度较高时，也会出现振型丢失现象。表4不同温度载荷工况Table4Differenttemperatureloadconditionsrl1巨；r卜1卜◆_1■1瞄l●卜●J。：l_LJ。、j1oj‘?i黎嗲夺I1K’1．1阎雾，◇_卜◆一ij0图17工况A时前4阶模态振型Fig．17FirstfourvibrationmodeshapesinCaseA图18工况B时前4阶模态振型Fig．18FirstfourvibrationmodeshapesinCaseB3模态跃迁和丢失机理分析结构模态的跃迁及丢失现象以往主要发生在一些模态密集结构，像大型空间网壳、网架、索穹顶结构、叶片轮盘等具有周期对称特点的结构。这些结构由于其基频较低且模态比较密集，很小的参数扰动就有可能导致某些模态阶次的跃迁。在研究此类结构模态跃迁的机理时，应用较多的还是摄动法，将外界因素的影响作为一个小参数摄动引入振动特征方程，通过分析摄动前后特征值的变化来解释模态跃迁的机理。下面就基于摄动法分析，对随着环境温度的升高上述复合材料简支板的模态跃迁机理作出解释。在不考虑阻尼时，系统振动的特征方程为K。，Hj“一怂“M。H∥(4)式中：各量均为未受到热载荷之前的原始结构物理量；K0和M0分别为结构的刚度矩阵和质量矩阵；硝’为第i阶特征值，即固有频率，H∥为该特征值对应的特征向量，即固有模态振型。当结构受到高温载荷后，结构的质量、刚度矩阵都会发生变化，在这里引入一个微小量e，用它220553—9 航空学报来表示温度载荷对质量、刚度矩阵的改变，这样变化后的质量矩阵和刚度矩阵可以表示为K—Ko+￡K1(5)114一Mo+sMl(6)式中：sK，和eM。分别为刚度矩阵和质量矩阵的变化量，下标1表示1阶摄动。当质量矩阵和刚度矩阵由于温度变化而发生改变时，对应的特征值和特征向量也发生了变化。将特征值和特征向量按e展开成幂级数，有H“’一“j。+￡H{。+O(￡2)(7)|：l“’=|；L5。+￡A{。4-O(￡2)(8)其中：“；。和AP分别为特征向量和特征值的1阶摄动量，对于2阶及更高阶摄动量这里当作小量不予以考虑。将得到的式(5)～式(8)代入式(4)中可以得到：．KoH；。+K1H5。一A：。YloH{。4-A：。MlH5。4-A；。M。“；。(9)根据展开定理，将“P在原系统的特征向量H：站上展开则可以得到：比P一掣Hjl’+c跏52’4-⋯+c～∥一∑c‰≯s=l(10)将原系统的特征矢量归一化为(Hj。)1Ko(Hj。)1一d。A：订(11)(H55’)1Mo(H5。)1一艿。(12)将式(10)代入式(9)可以得到特征值的1阶摄动表达式为Af。=(“j。)1K。Hj”一A：。(u(o。)1M，“5。(13)a1’一一百1(“护)1M1uo(。(5一i)(14)厶d。=一可r上而(H≯)1K，H铲一^0A5。(Hj。)1M1H：。(s≠i)(15)从式(13)可以看出，相同的参数摄动对不同阶次特征值的摄动结果却是不相同的。当相邻的2个特征值满足硝’≤Aj汁”，并且1阶摄动满足e(AP一世_1’)>A护硝_1’时，就可以得到A“’≥A““’。而在振动分析时，对于模态的阶次是严格按照特征值(固有频率)大小由小到大排列的，所以摄动后频率的变化导致低阶模态频率高于高阶模态频率，即发生了模态跃迁现象。4结论本文研究了高温环境下二维正交编织c／sic复合材料简支板的固有振动特性，通过分析不同因素对结构固有频率的影响，得到了不同温度下结构固有频率随温度的变化曲线，同时还对结构模态振型的演化规律进行了深入研究。利用有限元法仿真重现了在高温环境下复合材料壁板结构振动实验中出现的模态振型跃迁及丢失现象，并运用摄动理论对这些现象的发生机理进行了解释。得到以下结论：1)均匀温度场下，仅考虑材料参数随温度变化时，二维正交编织C／SiC复合材料板各阶固有频率有所下降，但是下降的幅度比较小。2)假设材料参数为一定值，仅考虑热应力对结构刚度影响时，二维正交编织c／sic复合材料板在临界屈曲温度(34℃)下前lo阶固有频率随着温度升高而降低；而在临界屈曲温度后，随着温度升高结构固有频率又开始增加。在大约247℃和403℃附近发生了2次模态振型跃迁现象。从图8还可以看出在850℃附近，第1、2阶振型又出现了交错现象，即高阶模态振型在较低阶频率处出现。3)与仅考虑热应力时的结果相比，综合考虑热物性和热应力时结构各阶频率值在每个温度点有所下降而且随着温度的升高下降幅度更加明显，但是总体变化趋势和仅考虑热应力时是一致的。4)在线性非均匀温度场下，当基础温度较低时，随着温度升高复合材料壁板结构的低阶模态振型出现了丢失，高阶模态振型提前出现，而且这种模态丢失现象是按照由低阶振型到高阶振型的顺序依次发生的。当基础温度较高时，前面的很多阶振型都已丢失。与较低基础温度下的振型相比，此时的模态振型也更为复杂。参考文献r1]SPAINCV，SOISTMANNDL，LINVILLETW．Integrationofthermaleffectsintofiniteelementaero—。thermo—_elasticanalysiswithillustrativeresults[c]∥sixthNationalAeroSpacePlaneTechnologySymposium，1989：PaperNO．9．r2]SHUKI，AA，GORDONRW，HOLLKAMPJJ．Numer一220553．10 航空学报icalinvestigationofthesnap—throughresponseofacurved，clamped—clampedplatewiththermalandrandomloading[C]／／AIAA／ASME／ASCE／AHS／ASCStructures，StructuralDynamics，andMaterialsConference，AIAA／ASME／ASCE／AHS／ASCAdaptiveStructuresConference[11]10th．Reston：AIAA，2008．[3]夏巍，杨智春．热环境下复合材料壁板的振动特性分析[J]．应用力学学报，2005，22(3)：359—363．XIAW，YANGZC．Vibrationanalysistocompositepan一[4][5]E6][7][8][9][10]elsinthermalenvironment[J]．ChineseJournalofAppliedMechanics，2005，22(3)：359363(inChinese)．刘芹，任建亭，姜节胜，等．复合材料层合板非线性热振动分析[J]．动力学与控制学报，2005，3(1)：78—83．LUIQ，RENJT，JIANGJS，eta1．Nonlinearthermalvibrationcharacteristicanalysiscompositelaminatedplates[J]．JournalofDynamics＆Control，2005，3(1)：78—83(inChinese)．SINGHVK，PANDASK．Nonlinearfreevibrationanal—ysisofsingle／doublycurvedcompositeshallowshellpanels[J]．Thin—WailedStructures，2014，85：341—349．杨莉，郝育新，杨绍武，等．热环境下FGM圆柱曲板的模态分析[J]．材料热处理学报，2014，35(8)：198—205．YANGI。，HAOYX，YANGSW，eta1．Modalanalysisoffunctionallygradedcylindricalpanelinthermalfield[J]．TransactionsofMaterials8LHeatTreatment，2014．35(8)：198-205(inChinese)．GENGQ，LIH，LIY．Dynamicandacousticresponseofaclampedrectangularplateinthermalenvironments：Ex—perimentandnumericalsimulation[J]．JournaloftheAcousticalSocietyofAmerica，2014，135：2674—2682．吴大方，王岳武，蒲颖，等．高超声速飞行器复合材料翼面结构1100℃高温环境下的热模态试验[J]．复合材料学报，2015，32(2)：323—331．WUDF，WANGYW，PUY，eta1．Thermalmodaltestofcompositewingstructureinhightemperatureenviron—mentsupto1i00℃forhypersonicflightvehicles[J]．ActaMateriaeCompositaeSinica，2015，32(2)：323—331(inChinese)．程昊，吴振强，李海波，等．高温环境下典型结构热模态分析[c]∥北京力学会第20届学术年会．北京：力学学会，2014．CHENGH，WUZQ，LIHB，eta1．Modalanalysisoftypicalstructureinhotenvironment[C]ffThe20thAnnualMeetingofBeijingSocietyofTheoreticalandAppliedMe—chanics．Beijing：BeijingSocietyofTheoreticalandAp—pliedMechanics，2014(inChinese)．史晓鸣，杨炳渊．瞬态加热环境下变厚度板温度场及热模态分析[J]．计算机辅助工程，2006，15(B09)：15-18．[12][13][14][15][16][17][18][19]220553．11SHIXM，YANGBY．Temperaturefieldandmodeanal—ysisofflatplatewiththermalenvironmentoftransientheating[J]．ComputerAidedEngineering，2006，15(B09)：15—18(inChinese)．黄世勇，王智勇．热环境下的结构模态分析[J]．导弹与航天运载技术，2009，29(5)：50—52．HUANGSY，WANGZY．Thestructuremodalanalysiswiththermalenvironment[J]．Missiles＆SpaceVehicles，2009，29(5)：5052(inChinese)．CHENGH，LIHB，ZHANGW，eta1．Effectsofradia—tionbeatingonmodalcharacte“sticsofpanelstructures[J]．JournalofSpacecraft＆Rockets，2015，52(4)：18．吕毅．平纹编织C／SiC复合材料力学性能的数值模拟[D]．西安：西北工业大学，2009：4051．LVY．Thenumericalsimulationofthemechanicalprop—ertiesofplainweavec／siccomposites[D]．Xi’an：NorthwesternPolytechnicalUniversity，2009：40—51(inChinese)．张超，许希武，郭树祥．二维二轴1X1编织复合材料细观结构模型及力学性能有限元分析[J]．复合材料学报，2011，28(6)：215-222．ZHANGC，XUXW，GUOSX．Microstructuremodelandfiniteelementanalysisofmechanicalpropertiesof2D1×1biaxialbraidedcomposites[J]．ActaMateriaeCom—positaeSinica，2011，28(6)：215-222(inChinese)．严雪，许希武，张超．二维三轴编织复合材料的弹性性能分析[J]．固体力学报，2013，34(2)：140—151．YANX，XUXW，ZHANGC．Analysisofelasticprop—ertiesof2Dtriaxialbraidedcomposites[J]．ChineseJour—nalofSolidMechanics，2013，34(2)：140—151(inChi—nese)．张歌，张程煜，乔生儒，等．2D-C／SiC高温弹性模量的预测[J]．机械强度，2013，35(1)：94—99．ZHANGG，ZHANGCY，QIA0SR，eta1．Predictionofelasticmodulusatelevatedtemperaturefor2D—C／SiC[J]．JournalofMechanicalStrength，2013，35(1)：94—99(inChinese)．WANGHL，ZHANGCY，LIUYS，eta1．Temperaturedependencyofinterlaminarshearstrengthof2D-C／SiCcomposite[J]．Materials&Design，2012，36：172176．马彦，陈朝辉．1800℃热处理对PIP法C／SiC复合材料结构和性能的影响[J]．材料工程，2015，43(4)：98—101．MAY，CHENZH．Effectof1800℃annealingonmi—crostructuresandpropertiesofC／SiCcompositesfabrica—tedbyprecursorinfiltrationandpyrolysis[J]．JournalofMaterialsEngineering，2015，43(4)：98—101(inChi—nese)．SUNZ，CHENX，SHAOH，eta1．Numericalmodeling 航空学报ofoxidized2DC／SiCcompositesinairenvironmentsbelow900℃：Mierostructureandelasticproperties[J]．AppliedCompositeMaterials，2016，23(4)：卜21．[20]MALEKZADEHP，SHAHPARISA，ZIAEEHR．Three——dimensionalfreevibrationofthickfunctionallygra—-dedannularplatesinthermalenvironment[J"]．JournalofSoundandVibration，2010，329(4)：425-442．(责任编辑：徐晓)Vibrationmodeevolutionof2DwovenC／SiCcompositepanelsinhotenvironmentHEErming*，CHENBing，CAOCunxianSchoolofAeronautics，NorthwesternPolytechnicalUniversity，Xi’an710072，ChinaAbstract：TostudytheIawofthevariationofnaturalvibrationfrequenc。lesof2DwovenC／SiCcompositepanelsinhotenvi-ronment，andfurtherrevealthemodejumpingandlosingphenomenainthethermalmodetest,themechanicsmodules，ther。malexpansioncoefficientandthermalconductivityof2DwovenC／SiCcompositepanelsafterhomogenizationarecalculatedbyvolume-averagedmethodbasedonitsmicroscopicmechanicsmodel．Thevibrationmodeevolutionsofthesimply‘suppor。tedrectangularpanelunderunifOrmtemperatureloadsandlinearnon-uniformtemperatureloadsarestudiedbasedonthemacroscopicmaterialpropertiesobtainedabove，andtherelationshipsbetweenhightemperatureandthenaturalfrequency，modeshapesareobtainedfromthemodeanalysis．TheresearchresultsshowthatundertheuniformtemperatureJcads，whenweconsideronlytheeffectoftemperatureontheelasticmodulus，thenaturalfrequencydecreaseslightlyasthetern—peratureincreases．Whenthethermalstresscausedbyhightemperatureisconsideredseparately，itisfoundthatthenaturalfrequencyincreasesafterthecriticalbuckledtemperature，andasthetemperatureincreases，themodeshapebeingstag-geredwhichresultsinsomemodeshapesappearearly，inotherwords，modejumpingoccursinthosemoments．Ifbothofthetwofactorsareconsidered，thenaturalfrequencywillbesmallerateachtemperaturecomparedwithonlyconsideringthethermalstresscausedbyhightemperature．However，thetrendsaresimilar．Underthenon-uniformhightemperatureloads，itshowsthatsomemodalshapesmaybelostbycomparingthemodeshapesunderdifferentoperatingconditions，andthelossisnotmessy，butgoesfrOmlow—ordertohigh-orderfrequency．Keywords：2DC／SiC；compositepanel；thermalmode；naturalfrequency；modejumpingandlosingphenomennaReceived：2016—06—22；Revised：2016-10—19；Accepted：2017-03—07；Publishedonline：2017—04—0113：32URL：WWW．cnki．net／kcms／detail／111929V201704011332．006．htmlFoundationitems：NationalNaturalScienceFoundationofChina(51675426)；AerospaceScienceandTechnologyInnovationFund(2014KC010043)*Correspondingauthor．E-mail：heerming@nwpueduca220553．12