共查询到17条相似文献,搜索用时 703 毫秒
1.
跨音速颤振计算方法研究 总被引:7,自引:3,他引:4
采用快捷的动态弹性变形技术生成三维块结构动态贴体网格。结构运动方程采用二级精度的龙格—库塔时间推进。气动力求解则用非定常欧拉方程的双时间有限体积推进,外时间为物理时间,与结构运动方程同步;对每个真实物理时间步采用五步Runge Kutta时间推进进行虚拟定常迭代。针对可压流颤振计算中存在的质量不相似问题,采用变质量、变刚度的方法计算出质量匹配点处的颤振速压值,并根据颤振速压随质量或刚度倍数的变化趋势,得出可用的飞机跨音速颤振速压。对一飞机简化外形的跨音速颤振特性计算得到了合理的计算结果。 相似文献
2.
3.
两种跨声速气动弹性问题分析研究 总被引:1,自引:0,他引:1
在非结构运动网格基础上,采用中心有限体积法进行空间离散和双时间方法进行时间推进求解非定常欧拉方程.通过与气动力方程的联立求解,在时域内用四步龙格-库塔方法求解结构运动方程.分析和研究了二维嗡鸣和三维机翼颤振这两种跨声速非线性气动弹性问题.二维嗡鸣问题的研究考虑了翼面-舵面系统的缝隙间网格运动、缝隙对嗡鸣的影响和扰流片对嗡鸣的抑制.耦合多自由度Lagrange结构运动方程数值模拟了三维机翼的颤振问题.通过跨声速标模算例AGARD445.6机翼的颤振计算,计算的颤振临界速度与实验值有5%左右的误差,验证本方法的正确性.由于本方法是在对外形具有良好普适性的非结构动网格基础上完成的,具有良好的工程实践价值. 相似文献
4.
机翼跨音速颤振的频域计算方法是以给定机翼模态分布下机翼上各点的模态值作为运动幅值,以三维非定常Euler方程为控制方程,采用有限体积法和双时间推进,求解三维机翼简谐运动下的非定常气动力。所求得的气动力作为已知值运用于颤振方程,利用v-g法进行求解。对得到的一系列的阻尼、速度和频率进行了线性插值,从而得到颤振速度和颤振频率。 相似文献
5.
为了提高叶轮机颤振预测时的计算效率,同时兼顾求解的稳定性,有必要选取最适合的时间推进方法。基于降阶结构
动力学方程的双向流固耦合方法常用于压气机气弹性能分析,针对气动弹性控制方程中结构动力学方程考查了多种时间推进方
法对计算结果的影响。通过单自由度弹簧系统的数值求解探讨了不同时间推进方法的特点;选取实际压气机算例NASA Rotor 67
进行颤振流固耦合分析,进行大量数值试验找到了各时间推进方法需要的最大时间步长,以此为依据对比出不同时间推进方法的
计算效率。结果表明:采用单自由度弹簧系统和压气机颤振耦合求解得出的结论基本一致,说明时间推进方法的性能更多地与方
法本身的数学性质有关,但是流固耦合求解时流场和结构的信息交互会带来一定影响。以实际颤振算例叶片瞬态响应振幅的对
数衰减率作为衡量精度的标准,经典龙格-库塔方法是颤振耦合计算达到相同精度耗时最少的方法,其次是4阶隐式Adams方法
和Newmark方法。 相似文献
6.
跨声速风扇的流固耦合的颤振分析(英文) 总被引:5,自引:0,他引:5
发展了时间推进的叶片颤振的数值方法,采用了时间推进求解流体和固体相互作用的过程。气动模型是基于求解三维雷诺平均的N-S方程,采用了多块结构化网格的有限体积格式,对流和耗散通量的计算使用了二阶迎风格式和中心格式。在变形的动网格上流体运动的守恒型方程的求解采用双时间步,隐式格式和多重网格方法。叶片振动采用了振型叠加的线性气动弹性模型。该方法在气动弹性标准算例4进行了验证,并用于求解跨音速风扇的颤振问题。 相似文献
7.
主要的研究工作是建立了以三维非定常N—S方程为主管方程的计算非定常气动力的高效计算方法,采用颤振方程和N—S方程耦合求解,时间推进的方法计算结构响应特性。发展了一种对既带有结构非线性又含有气动非线性的多重非线性的复杂颤振系统响应特性的分析方法。对全动操纵面转轴带有硬弹簧和间隙型非线性刚度的颤振系统进行了分析,计算结果表明这类系统具有复杂的极限环特性。 相似文献
8.
提出了一种高效的颤振预测分析方法,将之应用于跨声速颤振边界分析及结构参数影响研究中。本方法基于频域颤振分析 V-g 方法,为气弹系统提供一定的人工阻尼使之保持简谐运动状态,从而将结构动力学方程转换到频域内。然后通过高效的谐波平衡法得到一系列简谐运动频率下的气动力描述函数矩阵,结合频域结构方程,将气弹系统的稳定性问题转化为广义特征值求解。结果表明:本方法计算得到的颤振边界与高精度的时间推进方法非常吻合,分析效率有了明显的提升,而且当结构参数发生变化后,只需进行若干次广义特征值求解即可得到新的颤振边界,无需像时间推进方法一样开展大量的气动/结构耦合数值模拟。 相似文献
9.
基于预估校正和嵌套网格的虚拟飞行数值模拟 总被引:1,自引:0,他引:1
针对导弹虚拟飞行数值模拟问题,发展了空气动力学/飞行力学数值计算方法和软件。控制方程为非定常雷诺时均Navier-Stoker(RANS)方程和刚体六自由度运动方程;流场求解器为有限体积法结构网格求解器,时间推进采用双时间步法,湍流模型为Spalart-Allmaras一方程模型;采用Adams预估校正法实现飞行力学方程与流场控制方程的耦合计算;使用嵌套网格方法模拟多体运动。首先模拟了美国国家航空航天局(NASA)窄条翼导弹模型纵向虚拟飞行,研究耦合方式和时间步长的影响。仿真结果表明,双时间步三阶Adams耦合方法,同等精度下可以显著增大时间步长,缩短仿真时间。最后,采用该方法模拟了导弹自由摇滚特性和纵向虚拟飞行,模拟结果与试验值吻合较好。 相似文献
10.
通过交替求解流体运动方程和结构动力方程来计算时域内二维机翼的结构动力响应,并研究了结构参数对系统稳定性的影响,使用控制中的相平面图对计算结果进行分析,来判断系统的结构稳定性。计算中流体运动方程采用二维非定常Euler方程,并用Jameson提出的有限体积方案,四步Runge-Kutta时间推进求解;结构模型采用二自由度PAP模型。经分析,计算结果合理。 相似文献
11.
12.
13.
欧拉方程计算非定常流动的一种数值方法 总被引:4,自引:0,他引:4
提出了利用欧拉方程计算非定常流动的一种数值方法。在利用中心差分的有限体积法和双时间推进法求解非定常欧拉方程过程中,通过在固定物面边界上满足完全的非定常边界条件,计算出非定常流场。算例结果表明,本文方法具有高效、简便的特点,而且能描述激波的非定常运动和变化情况。 相似文献
14.
气动结构耦合是研究静气动弹性问题的关键所在。采用三维Euler方程为控制方程,以有限元方法计算结构弹性变形,用细平板样条法(TPS)实现气动结构数据间的转换,进行机翼静气动弹性数值的模拟。算例结果表明,本方法是一种具有很高计算精度和工程实用性的气动弹性模拟方法。 相似文献
15.
16.