共查询到19条相似文献,搜索用时 176 毫秒
1.
2.
3.
基于CFD的带控制舵导弹的动导数计算 总被引:4,自引:4,他引:0
介绍了基于计算流体动力学(CFD)的阻尼导数预测方法以及其在复杂带控制舵导弹外形的俯仰/滚转阻尼导数预测中的应用.方法采用求解三维非定常ALE(arbitrary Lagrangian-Eulerian)形式RANS(Reynolds-averaged Navier-Stokes)方程加湍流模型的方法得到模型强迫运动的气动力时间历程,再通过绕模型非定常运动过程得到的气动力迟滞曲线计算得到俯仰及滚转阻尼导数.为了验证采用数值方法的正确性,针对基本带翼导弹外形,采用强迫模型按正弦形式振动,计算得到了超声速范围内马赫数从1.5到2.5的俯仰及滚转稳定性参数,通过与风洞试验及弹道数据结果的比较发现,两者吻合较好. 相似文献
4.
本文从非定常Euler方程出发,应用谐振摄动法计算了球锥体超声速、高超声速下的俯仰阻尼导数和。Euler方程的求解利用SCM方法和MacCormack的两步算法。文中用任意攻角下不同俯仰轴的气动导数转换公式来简化计算,用较粗网格可以得到较精确的结果,程序编制相对简单。一次计算可以得到定常绕流解及一系列钝头率、不同俯仰轴下的多种动导数结果,在ELXSI6400机上约需20~40分钟。所得结果与试验结果符合程度令人满意。 相似文献
5.
本文提出了LTRAN2的一种改进方案BTRAN2用来分析翼型的各种频率的非定常运动。用E-O调转换格式和ADI方法求解了完全的二维非定常跨音速小扰动位流方程,并用单调转换的AF2格式计算了定常跨音速小扰动方程,以此作为非定常计算的初场。本文给出了带后缘正弦振荡襟翼的NACA 64 A006翼型的绕流和做正弦俯仰振荡的NACA 64 A010翼型的绕流计算结果,它们与Euler方程解或实验数据很吻合。 相似文献
6.
本文采用Euler方程和TVD格式数值求解了小幅振动旋成体的非定常跨声速绕流问题,,假设在绕小幅振旋成体的非定常气流中,物理量的变化为平均定常扰动物理量的高阶小量,然后对Euler方程进行摄动展开,导出前两阶(零阶及一阶)摄动量的控制方程和相庆的边界条件,并进行数值求解。文中仅给出了振动旋成体的非定常压力分布。 相似文献
7.
8.
9.
飞行器大迎角下俯仰静、动导数的数值计算 总被引:8,自引:0,他引:8
采用一种适用于非定常流动计算的逐次超松弛(SOR)时间推进方法耦合求解非定常Navier-Stokes方程和强迫运动方程,数值模拟了NACA0012翼型、弹道外形(HBS)和有翼导弹等飞行器在不同初始迎角下做强迫俯仰振荡的黏性动态流场,给出动态气动力及力矩系数的时间历程;采用积分法获得飞行器的静导数和动导数,将计算结果与实验结果及半经验理论预测方法进行比较并分析飞行器在大迎角下动态特性的非线性特征;计算结果与实验值相吻合,表明本文采用的非定常数值模拟方法在大迎角下能够较准确地模拟飞行器的动态特性,具有工程应用价值。 相似文献
10.
考虑地面效应的翼型动态特性数值模拟 总被引:1,自引:0,他引:1
为了研究考虑地面干扰的飞行器非定常气动特性,基于滑移网格技术,通过求解二维非定常Euler方程,对不同近地高度下的NACA0012翼型的动态振动特性进行了数值模拟。分析了存在地面干扰时的升力系数以及俯仰力矩系数迟滞环变化,并结合强迫振动法,计算了翼型各个近地高度下的纵向组合动导数。结果表明,地面效应不仅对定常流场产生影响,更显著地影响了非定常气动力及力矩,近地高度越小,这种效应越明显,使得升力系数的迟滞环面积变小,幅值减小,而俯仰力矩系数的迟滞环变化不规律,力矩系数的动导数随着近地高度减小而增大(绝对值减小),纵向的动态稳定性产生损失。地面效应干扰在飞行器非定常气动研究中应该引起重视。 相似文献
11.
二维非定常欧拉方程的有限元解法 总被引:1,自引:0,他引:1
本文通过求解时间相关的欧拉方程,对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法,并用VanLeer逆风通量分裂格式来计算控制体通量。隐式格式的采用,大大提高了计算效率;同时,通过插值使得空间离散达到了二阶精度。另外,本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动,从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了NACA0012翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好, 相似文献
12.
13.
一、基本理论 为了能够准确地计算绕弹性体,特别是跨音速范围的非定常气动力,有必要建立一种快速生成的随时间变化的贴体计算网格,即通过变换ζ=ζ(x,y,t),η=η(x,y,t),ζ=ζ(x,y,t),τ=t (1)将边界形状随时间变化的物理流场(x,y,z)变换到便于计算的矩形计算空间(ζ,η,ζ)内。这样,与计算空间内计算节点相对应的流场节点的位置是随时间改变的。因此在进行流场计算时,不仅要确定流场内每一时刻计算节点的位置,同时还应确定其随时间的变化速率(x?,y?,z?)。在此主要以二维问题进行讨论。 相似文献
14.
为了从时间离散精度的角度评估傅里叶时间谱方法的实际计算效率,对用该方法求解二维跨音速流动问题的解进行了精度研究。通过求解Euler方程模拟俯仰振动NACA0012翼型的周期性非定常流动。结果表明:傅里叶时间谱方法可以使用很大的时间间隔较为准确地模拟有激波存在的周期性跨音速流动;对于此类流动问题,相比二阶向后差分公式法,傅里叶时间谱方法在预测翼型表面压力系数积分量特别是升力系数的时间变化规律方面具有明显的计算效率优势;用于模拟翼型表面压力系数自身的时间变化规律时,傅里叶时间谱方法具有不低于二阶向后差分公式法的计算效率。 相似文献
15.
基于非结构重叠网格的二维外挂物投放模拟 总被引:3,自引:0,他引:3
本文发展了基于非结构重叠网格的动态网格技术,并提出了新的相邻单元搜索算法。非结构重叠网格集非结构网格和重叠网格的优点于一身,适合处理复杂外形和多体相对运动问题。新的相邻单元搜索算法在碰到物面时能继续搜索下去,搜索效率大大提高。非定常Euler方程的求解,采用有限体积法,时间上采用双时间步长推进。为验证算法的有效性,先对绕振荡NACA0012翼形的流动进行了模拟;然后通过耦合求解刚体三自由度运动方程,模拟了二维外挂物投放过程,结果表明非结构重叠网格在处理多体相对问题时具有优越性。 相似文献
16.
复合振动翼型跨声速非定常流计算 总被引:2,自引:1,他引:2
在文献[1]、[2]的基础上,本文把绕翼弦中点作旋转振动的翼型跨声速流计算推广到翼型作水平、垂直方向的平移振动以及平移与旋转振动的叠加。本方法主要特点是网格与翼型一起运动,对无穷远及翼面边界条件作了精确的处理以及计算网格数很少。本文结果对进一步研究旋翼、叶轮转子等的非定常气动力将有重要的意义。 相似文献
17.
用N-S方程计算翼型非定常粘性大攻角绕流 总被引:1,自引:0,他引:1
本文给出了用NS方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二级时均可压缩完全NS方程;湍流模型采用双层代数涡粘性模型 ̄[1]。使用近似因式分解ADI差分格式离散求解,网格是用保角变换方法生成的相对翼型固定的C型网格。本文给出两类典型非定常绕流数值模拟结果:翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较,表明了本方法准确、高效的特点。 相似文献
18.
介绍了一种从频域出发模拟二维周期性非定常翼型绕流的数值方法.考虑到很多工程中遇到的非定常流动问题都是周期性的,因此非定常流场中的守恒变量可以展开为傅立叶级数,其系数随空间位置的不同而改变.将这样的谐波假设引入欧拉方程,经过处理的方程不显式地含有时间导数项,实际求解的方程将具有定常方程的形式.在使用多重网格技术后,计算可以很快达到收敛.对NACA64A010翼型典型的周期性非定常绕流进行的数值模拟表明,这种方法在能够给出较为准确的计算结果的同时,与传统的双时间方法相比,可以大幅提高计算效率. 相似文献
19.
本文采用MacCormack两步显式格式,用有限体积法求解了二元跨声速欧拉流。推导了物面边界条件,采用了特征远场边界条件及远场环量修正。利用保角变换方法生成O型贴体网格,并修改得到了一种在激波处局部加密的半自适应贴体网格。采用多重网格及焓阻尼加速收敛技术计算了NACA0012翼型的跨声速气动特性,得到了十分满意的收敛过程和计算结果。 相似文献