本发明专利技术涉及一种气动弹性分析方法，提供了一种高超声速飞行器热环境下气动弹性力学特性分析方法。使用活塞理论对高超声速飞行器全机有限元模型进行频域非定常气动力计算，在此基础上考虑了在高超声速热环境中，模型受到的气动热效应，忽略气动热输入对气动力输入和弹性力输入的弱耦合效应，仅考虑气动热载荷作用下的结构温度的稳态特性，在进行定常气动力求解之后采用参考焓法求得飞行器表面的热流密度，进而计算出飞行器表面的稳态温度分布，并求得此时结构的实际等效刚度矩阵，并采用工程方法求解出临界颤振速度。本发明专利技术解决了气动热环境下超声速飞行器气动弹性分析问题，通过对颤振速度的分析从而提高高超声速飞行器气动弹性性能。
【技术实现步骤摘要】
一种高超声速飞行器热环境下气动弹性力学特性分析方法
本专利技术涉及一种航空航天领域的气动弹性分析方法，具体地说，是高超声速飞行器考虑热效应情况下的气动弹性分析方法。
技术介绍
在气动弹性学科的研究初期，通常采用的方法是以结构和流体动力学模型为基础进行线性化假设，颤振边界问题归根到底是非自伴系统模型的特征值求解。非自伴系统模型基于颤振运动微分方程而构建，对于位于颤振临界点附近的简谐振动特性分析也能获得可靠的结果。该方法的不足之处是，需要对模型进行线性化假设，系统中若带有非线性的气动弹性问题，就显无能为力了。伴随着飞行器理论的发展，结构上采用高升阻比轮廓、机身采用智能材料、运用智能变形技术，气流三维效应、结构间的相互干扰等非线性气动力因素，使得材质、负荷、装配缝隙和外形柔性形变等结构非线性效应变得更加突出。高超声速气动热弹性问题是典型的多场耦合问题，通过对此一次求解而得出结果在工程实践中可望而不可及。在气动热环境下，原本由气动力、惯性力和弹性力三者相互作用的耦合关系，又会加入气动热耦合，使这四者相互影响形成不同程度的耦合作用关系。高超声速气动热弹性力学问题求解的关键之一就是获得适合于气动弹性分析的非定常气动力。非线性的引入，使得非定常气动力求解方法的使用受到极大限制，比如：基于势流理论和模态叠加法只能局限于时域求解微分方程，获得系统响应曲线，从而对系统的稳定性进行判别。时域分析的好处是可以让模型脱离结构/气动非线性的限制，进行多场耦合状态下的时域仿真，并且能够结合工作模态分析中的最小二乘法曲线拟合、小波变换法、分段滑动、自回归滑动平均法等具体手段，从而确定系统频率、阻尼系数等气动弹性的重要临界参数。但是计算过程复杂、工作量大，尤其在时域方法中做CFD/CSD的直接耦合，对计算机的配置提出了很高的要求，限制了其在工程实践中的使用范围。从国内外研究资料现状可以发现，对大高超声速飞行器进行热气动弹性分析方面的研究工作还有待发掘，利用现有商用有限元软件，结合工程计算方法，对具有真实翼型的高超声速飞行器全机模型进行热气动弹性分析对于高超声速飞行器的发展有着很强的现实意义，怎样使用热颤振分析理念，对高超声速飞行器进行热气动弹性计算分析已成为了一个新的命题。
技术实现思路
本专利技术针对上述现有技术状况而设计提供了一种高超声速飞行器热环境下气动弹性力学特性分析方法，其目的是解决在气动热环境下，对高超声速飞行器进行气动弹性分析，判断高超声速飞行器气动弹性性能满足技术要求与否，进而提高高超声速飞行器气动弹性性能。本专利技术所采用的技术方案包括以下步骤：(1)建立高超声速飞行器全机几何模型和有限元模型，包括以下内容：(a)在三维造型软件中建立高超声速飞行器全机几何模型，几何模型采用有利于机身和超燃冲压发动机实现一体化布局设计的二维升力体外形，同时机身和全动平尾采用了双楔形薄翼型，垂尾采用的是梯形翼型，其中最易发生颤振的平尾部分与机身由转轴相连；(b)在将几何模型读入MSC.Patran之后，建立有限元模型并对模型进行有限元网格划分，根据全机的组成部件分别处理，机身蒙皮、垂尾表面蒙皮和平尾表面蒙皮等结构部分选择壳单元建模，并采用MSC.Patran中的四边形单元进行自动划分，机身内部和平尾内部采用梁结构，并赋予梁单元属性，飞机头部、垂尾前缘和平尾前缘部分采用实体单元建模，使用MSC.Patran中MatchParasoildFaces/NeighborSolidList命令的智能控制功能可以很好的保障各部分之间单元的连续性并使用MSC.Patran中的四面体单元进行自动划分；(c)在材料的选用上，机身主要选用了高强度、高密度、耐高温的金属材料结构，外层用氧化铝隔热层覆盖以保护其内部结构，其中在机头和平尾前部都采用了C-C复合材料以耐受高温，机头C-C材料后方是Densalloy180钨合金，同时对刚度和密度进行了适当折减，在逐一对各个部件划分有限单元并赋予材料属性的基础上，对各部件加入集中质量，使各个部件的质量特性更加合理，系统及燃油质量等，用集中质量卡(CONM2)施加于相应质心上，并用MPC元约束在主盒段上，全动平尾部分采用的是Haynes230镍合金蒙皮和梁结构骨架，用BEAM单元建立刚轴大梁来模拟水平安定面和升降舵的刚度特性，CONM2单元模拟结构质量，水平安定面与升降舵之间采用约束X、Y、Z三个方向自由度的MPC来模拟铰链连接，采用BAR、ROD单元来模拟作动器和摇臂；(2)针对高超声速飞行器全机有限元模型，求解其在初始状态下振动模态：根据模型的对称性，取半机身，根部固支，通过MSC.Nastran求解序列SOL103对模型在初始状态下进行振动模态分析，再用MSC.Patran进行后处理分析得到的前n阶模态；(3)计算模型表面受到的热流密度场分布：机翼前缘部分受到的气动加热现象最为明显，高超声速飞行器薄型机翼机身的前缘部分和翼面部分可以忽略厚度因素看做平板来进行计算：参考粘性系数μ*通过萨特兰表达式求解得出：参考密度ρ*通过状态方程表达式求解得出：同时可以得到参考雷诺数，如下：斯坦顿数通过雷诺比拟式求解得出：最终高超声速飞行器表面热流通过以下公式求得：Qaero＝St*ρ*V∞cρ(Tr-Tω)沿着机身x方向将机身分成若干等份，假设每份所处的热流密度一定，则热流场为若干个不相互连续的离散温度值构成，在初始环境温度一定时，计算得到的热流场分布；(4)将热流密度场加载到有限元模型上，计算其有限元模型的稳态温度场分布：考虑流经飞行器表面的气流，飞行器表面任一点的τω等价与相应的可压流的τω，不可压缩气流的温度升高至一个给定的参考温度值T*：T*＝0.5Tω+0.22Tr+0.28T∞其中：T*表示飞行器表面温度，表示边界层外缘温度，由气动力计算部分参数可以得到，Tr为恢复温度，根据算出的热流密度场，在MSC.Patran中将其加载在有限元模型表面，同时加上适当的辐射值加以平衡，通过MSC.Nastran求解序列SOL153对其进行温度分布计算，得到其有限元模型的稳态温度场分布；(5)根据稳态温度场分布情况得到稳态热应力变形后的结构参数：令结构初始单元线性刚度矩阵为K0，为热应力作用生成的多余刚度为Kσ，则结构受热效应作用后的总的应力刚度矩阵等效为K＝K0+Kσ则在热效应作用下的结构振动表达式为式中，M表示质量矩阵，表示振型，ω为振频，此时的动力学方程即描述了为结构在热载荷作用下结构振动特性，即简化为求上式的广义特征值问题，根据得到的温度分布情况，进一步的通过MSC.Nastran求解序列SOL153对对结构进行静态分析，得到稳态热应力变形后的结构参数，包括刚度矩阵，此时的刚度矩阵即为考虑热效应下的等效刚度矩阵；(6)建立气动网格模型：应用MSCFlightloads中的气动弹性模块，将模型划分成数目适合的气动分区，每个气动分区又会分成数目适合的气动片，获得气动网格参数；(7)根据结构的刚度矩阵等结构参数和划分的气动网格参数，使用活塞理论进行线性的频域非定常气动力计算：初始状态下的气体压力大小为p∞，密度为ρ∞，声速大小为a∞，假设在等熵条件下其中，γ为气体比热比。当地流音速大小为两边同时积分可以得到无穷远处的p＝p∞，且此时本文档来自技高网...
【技术保护点】
一种高超声速飞行器热环境下气动弹性力学特性分析方法，其特征在于，包括以下步骤：(1)建立高超声速飞行器全机几何模型和有限元模型，包括以下内容：(a)在三维造型软件中建立高超声速飞行器全机几何模型，几何模型采用有利于机身和超燃冲压发动机实现一体化布局设计的二维升力体外形，同时机身和全动平尾采用了双楔形薄翼型，垂尾采用的是梯形翼型，其中最易发生颤振的平尾部分与机身由转轴相连；(b)在将几何模型读入MSC.Patran之后，建立有限元模型并对模型进行有限元网格划分，根据全机的组成部件分别处理，机身蒙皮、垂尾表面蒙皮和平尾表面蒙皮等结构部分选择壳单元建模，并采用MSC.Patran中的四边形单元进行自动划分，机身内部和平尾内部采用梁结构，并赋予梁单元属性，飞机头部、垂尾前缘和平尾前缘部分采用实体单元建模，使用MSC.Patran中Match Parasoild Faces/Neighbor Solid List命令的智能控制功能可以很好的保障各部分之间单元的连续性并使用MSC.Patran中的四面体单元进行自动划分；(c)在材料的选用上，机身主要选用了高强度、高密度、耐高温的金属材料结构，外层用氧化铝隔热层覆盖以保护其内部结构，其中在机头和平尾前部都采用了C‑C复合材料以耐受高温，机头C‑C材料后方是Densalloy180钨合金，同时对刚度和密度进行了适当折减，在逐一对各个部件划分有限单元并赋予材料属性的基础上，对各部件加入集中质量，使各个部件的质量特性更加合理，系统及燃油质量等，用集中质量卡(CONM2)施加于相应质心上，并用MPC元约束在主盒段上，全动平尾部分采用的是Haynes230镍合金蒙皮和梁结构骨架，用BEAM单元建立刚轴大梁来模拟水平安定面和升降舵的刚度特性，CONM2单元模拟结构质量，水平安定面与升降舵之间采用约束X、Y、Z三个方向自由度的MPC来模拟铰链连接，采用BAR、ROD单元来模拟作动器和摇臂；(2)针对高超声速飞行器全机有限元模型，求解其在初始状态下振动模态：根据模型的对称性，取半机身，根部固支，通过MSC.Nastran求解序列SOL103对模型在初始状态下进行振动模态分析，再用MSC.Patran进行后处理分析得到的前n阶模态；(3)计算模型表面受到的热流密度场分布：机翼前缘部分受到的气动加热现象最为明显，高超声速飞行器薄型机翼机身的前缘部分和翼面部分可以忽略厚度因素看做平板来进行计算：参考粘性系数μ*通过萨特兰表达式求解得出：μ*=(T*288.15)1.5398.55T*+110.4×1.7894×10-5]]>参考密度ρ*通过状态方程表达式求解得出：ρ*=p∞RT*]]>同时可以得到参考雷诺数，如下：Rex*=ρ*V∞xμ*]]>斯坦顿数通过雷诺比拟式求解得出：St*=cf*21(Pr)2/3]]>最终高超声速飞行器表面热流通过以下公式求得：Qaero＝St*ρ*V∞cρ(Tr‑Tω)沿着机身x方向将机身分成若干等份，假设每份所处的热流密度一定，则热流场为若干个不相互连续的离散温度值构成，在初始环境温度一定时，计算得到的热流场分布；(4)将热流密度场加载到有限元模型上，计算其有限元模型的稳态温度场分布：考虑流经飞行器表面的气流，飞行器表面任一点的τω等价与相应的可压流的τω，不可压缩气流的温度升高至一个给定的参考温度值T*：T*＝0.5Tω+0.22Tr+0.28T∞其中：T*表示飞行器表面温度，表示边界层外缘温度，由气动力计算部分参数可以得到，Tr为恢复温度，根据算出的热流密度场，在MSC.Patran中将其加载在有限元模型表面，同时加上适当的辐射值加以平衡，通过MSC.Nastran求解序列SOL153对其进行温度分布计算，得到其有限元模型的稳态温度场分布；(5)根据稳态温度场分布情况得到稳态热应力变形后的结构参数：令结构初始单元线性刚度矩阵为K0，为热应力作用生成的多余刚度为Kσ，则结构受热效应作用后的总的应力刚度矩阵等效为K＝K0+Kσ则在热效应作用下的结构振动表达式为式中，M表示质量矩阵，表示振型，ω为振频，此时的动力学方程即描述了为结构在热载荷作用下结构振动特性，即简化为求上式的广义特征值问题，根据得到的温度分布情况，进一步的通过MSC.Nastran求解序列SOL153对对结构进行静态分析，得到稳态热应力变形后的结构参数，包括刚度矩阵，此时的刚度矩阵即为考虑热效应下的等效刚度矩阵；(6)建立气动网格模型：应用MSC Flightloads中的气动弹性模块，将模型划分成数目适合的气动分区，每个气动分区又会分成数目适合的气动片，获得气动网格参数；(7)根据结构的刚度矩阵等结构参数和划分的气动网格参数，使用活塞理论...
【技术特征摘要】
1.一种高超声速飞行器热环境下气动弹性力学特性分析方法，其特征在于，包括以下步骤：(1)建立高超声速飞行器全机几何模型和有限元模型，包括以下内容：(a)在三维造型软件中建立高超声速飞行器全机几何模型，几何模型采用有利于机身和超燃冲压发动机实现一体化布局设计的二维升力体外形，同时机身和全动平尾采用了双楔形薄翼型，垂尾采用的是梯形翼型，其中最易发生颤振的平尾部分与机身由转轴相连；(b)在将几何模型读入MSC.Patran之后，建立有限元模型并对模型进行有限元网格划分，根据全机的组成部件分别处理，机身蒙皮、垂尾表面蒙皮和平尾表面蒙皮结构部分选择壳单元建模，并采用MSC.Patran中的四边形单元进行自动划分，机身内部和平尾内部采用梁结构，并赋予梁单元属性，飞机头部、垂尾前缘和平尾前缘部分采用实体单元建模，使用MSC.Patran中MatchParasoildFaces/NeighborSolidList命令的智能控制功能可以很好的保障各部分之间单元的连续性并使用MSC.Patran中的四面体单元进行自动划分；(c)在材料的选用上，机身主要选用了高强度、高密度、耐高温的金属材料结构，外层用氧化铝隔热层覆盖以保护其内部结构，其中在机头和平尾前部都采用了C-C复合材料以耐受高温，机头C-C材料后方是Densalloy180钨合金，同时对刚度和密度进行了适当折减，在逐一对各个部件划分有限单元并赋予材料属性的基础上，对各部件加入集中质量，使各个部件的质量特性更加合理，系统及燃油质量，用集中质量卡CONM2施加于相应质心上，并用MPC元约束在主盒段上，全动平尾部分采用的是Haynes230镍合金蒙皮和梁结构骨架，用BEAM单元建立刚轴大梁来模拟水平安定面和升降舵的刚度特性，CONM2单元模拟结构质量，水平安定面与升降舵之间采用约束X、Y、Z三个方向自由度的MPC来模拟铰链连接，采用BAR、ROD单元来模拟作动器和摇臂；(2)针对高超声速飞行器全机有限元模型，求解其在初始状态下振动模态：根据模型的对称性，取半机身，根部固支，通过MSC.Nastran求解序列SOL103对模型在初始状态下进行振动模态分析，再用MSC.Patran进行后处理分析得到的前n阶模态；(3)计算模型表面受到的热流密度场分布：机翼前缘部分受到的气动加热现象最为明显，高超声速飞行器薄型机翼机身的前缘部分和翼面部分可以忽略厚度因素看做平板来进行计算：参考粘性系数μ*通过萨特兰表达式求解得出：参考密度ρ*通过状态方程表达式求解得出：同时可以得到参考雷诺数，如下：斯坦顿数通过雷诺比拟式求解得出：最终高超声速飞行器表面热流通过以下公式求得：Qaero＝St*ρ*V∞cρ(Tr-Tω)沿着机身x方向将机身分成若干等份，假设每份所处的热流密度一定，则热流场为若干个不相互连续的离散温度值构成，在初始环境温度一定时，计算得到热流场分布；(4)将热流密度场加载到有限元模型上，计算其有限元模型的稳态温度场分布：考虑流经飞行器表面的气流，飞行器表面任一点的τω等价与相应的可压流的τω，不可压缩气流的温度升高至一个给定的参考温度值T*：T*＝0.5Tω+0.22Tr+0.28T∞其中：T*表示飞行器表面温度，表示边界层外缘温度，由气动力计算部分参数可以得到，Tr为恢复温度，根据算出的热流密度场，在MSC.Patran中将其加载在有限元模型表面，同时加上适当的辐射值加以平衡，通过MSC.Nastran求解序列SOL153对其进行温度分布计算，得到其有限元模型的稳态温度场分布；(5)根据稳态温度场分布情况得到稳态热应力变形后的结构参数：令结构初始单元线性刚度矩阵为K0，热应力作用生成的多余刚度为Kσ，则结构受热效应作用后的总的应力刚度矩阵等效为K＝K0+Kσ则在热效应作用下的结构振动表达式为式中，M表示质量矩阵，表示振型，ω为振频，此时的动力学方程即描述了结构在热载荷作用下结构振动特性，即简化为求上式的广义特征值问题，根据得到的温度分布情况，进一步的通过MSC.Nastran求解序列SOL153对结构进行静态分析，得到稳态热应力变形后的结构参数，包括刚度矩阵，此时的刚度矩阵即为考虑热效应下的等效刚度矩阵；(6)建立气动网格模型：应用MSCFlightloads中的气动弹性模块，将模型划分成数目适合的气动分区，每个气动分区又会分成数目适合的气动片，获得气动网格参数；(7)根据结构的刚度矩阵结构参数和划分的气动网格参数，使用活塞理论进行线性的频域非定常气动力计算：初始状态下的气体压力大小为p∞，密度为ρ∞，声速大小为a∞，假设在等熵条件下其中，γ为气体比热比，当地流音速大小为两边同时积分可以得到无穷远处的p＝p∞，且此时v∞＝0，得到得扰动压力为
【专利技术属性】
技术研发人员：马金玉，余胜东，
申请(专利权)人：温州职业技术学院，
类型：发明
国别省市：浙江;33
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。