基于非结构网格流场计算的网格重排序

提出基于非结构网格流场计算时的网格重排序,优化了网格信息在内存中的存贮位置,使得计算过程中调用数据的处理速度加快,有效提高了计算效率。以M6机翼及DLR-F4翼身组合体的跨音速无粘流场为算例,分别通过显式四步龙格-库塔推进格式和隐式LU-SGS推进格式来验证网格重排序后对计算效率的提高作用。原始网格均用Delaunay方法生长,并分别使用了两种重排序方法进行比较。计算结果表明:重排序后的计算效率可以提高到未排序时的3到5倍,很具实用价值。     

叶轮机气弹耦合方程求解的时间推进方法评估

为了提高叶轮机颤振预测时的计算效率,同时兼顾求解的稳定性,有必要选取最适合的时间推进方法。基于降阶结构动力学方程的双向流固耦合方法常用于压气机气弹性能分析,针对气动弹性控制方程中结构动力学方程考查了多种时间推进方法对计算结果的影响。通过单自由度弹簧系统的数值求解探讨了不同时间推进方法的特点;选取实际压气机算例NASA Rotor 67进行颤振流固耦合分析,进行大量数值试验找到了各时间推进方法需要的最大时间步长,以此为依据对比出不同时间推进方法的计算效率。结果表明：采用单自由度弹簧系统和压气机颤振耦合求解得出的结论基本一致,说明时间推进方法的性能更多地与方法本身的数学性质有关,但是流固耦合求解时流场和结构的信息交互会带来一定影响。以实际颤振算例叶片瞬态响应振幅的对数衰减率作为衡量精度的标准,经典龙格-库塔方法是颤振耦合计算达到相同精度耗时最少的方法,其次是4阶隐式Adams方法和Newmark方法。     

一种圆轨道航天器外热流通用计算方法

文章提出了一种圆轨道航天器外热流通用计算方法,采用解析法计算太阳直射热流,采用数值积分法求解行星红外和反照,无需求解行星中心到太阳矢量与航天器表面法线在局部球坐标系的经度差,能有效简化计算过程。在此基础上选取一种特殊圆轨道,对常用姿态(+Z对地,-Z对日,偏航)六面体卫星的瞬态外热流进行分析,并将偏航姿态瞬态外热流分析结果与热分析软件计算结果进行了对比,对比结果充分说明了此方法的准确性,它能快速准确地计算圆轨道航天器常用姿态外热流,可为各领域航天器热控设计和热分析提供重要参考。     

复合固体推进剂松弛模量与蠕变柔量转换关系

以粘弹理论为基础，采用Volterra算子，建立了复合固体推进剂松弛模量和蠕变柔量之间的相互转换关系。结果表明：由松弛模量计算蠕变柔量的理论值与实验值吻合；该转换表达式避免了数值积分法求解蠕变柔量值产生舍人误差和采用拉氏变换的繁琐计算。研究方法可以用于固体推进剂其它力学参数的转换计算，如体积松弛模量和体积蠕变柔量之间的转换计算。     

弹射救生数值仿真及不利姿态下救生性能分析

 国内弹射座椅救生性能预测分析主要依赖地面弹射试验。由于国内相关试验设备的缺乏,不利姿态下弹射座椅的救生性能很难通过试验获得,因此国内不利姿态下弹射座椅救生性能的分析研究比较薄弱。采用数值仿真的方法对弹射座椅的救生性能进行预测和分析,针对弹射座椅弹射过程的4个阶段分别建立了数学模型,其中自由飞阶段采用四元数法代替欧拉速率方程,以解决传统的六自由度方程出现的奇异性问题。采用四阶龙格-库塔法对数学模型进行求解,计算结果与地面弹射试验结果吻合较好。在此基础上,对某新型弹射座椅的弹射姿态轨迹进行计算,分析研究了弹射座椅在不利姿态下的救生性能,得到不利姿态下弹射座椅救生性能不佳的主要原因是座椅高速稳定性差以及飞机向下的牵连速度和飞机滚转角的影响,并根据分析结果提出了改善弹射座椅不利姿态下救生性能的基本方法。     

旋翼桨叶绕流Euler方程数值模拟

本研究应用格心格式有限体积法求解欧拉方程,模拟绕旋翼浆叶的流动。旋翼尾流的作用通过自由尾流分析来求解局部的诱导下洗速度以修正翼型攻角。选择O—H型式生成弦向和展向网格。通过算例计算为进一步寻求有效的Euler或NS方程数值解积累了经验。     

直升机尾传动系统扭转振动建模与特性

将直升机传动系统简化为轴段和当量圆盘的串联系统,建立了直升机尾传动系统扭转振动的等效多自由度动力学模型.模型中考虑了啮合齿轮对的综合啮合误差激励和尾减齿轮的啮合刚度.针对系统的扭转动力学方程,求得了系统的扭转振动响应,分析了直升机尾传动系统在轴的不同扭转刚度和齿轮的不同啮合刚度下的扭转振动的特性,结果表明:与尾斜轴相联的当量圆盘的扭转角位移始终比与水平轴相联的当量圆盘的扭转角位移的数值大,即与尾斜轴相联的尾减输出齿轮振动大于输入齿轮;当轴的扭转刚度变化时,水平轴相联的当量圆盘与尾斜轴相联的当量圆盘的扭转角位移变化的趋势相反;啮合刚度对系统扭转角位移的影响比较大,在建模时应当给予重视.      

非线性结构动响应的Taylor级数解法

Taylor级数方法是结构动力分析中一种新的时间积分方法,它在求解线性问题方面的理论和应用已经比较完善和成熟.将Taylor方法进一步用于非线性结构动响应的求解,对于非线性项可以表示为多元多项式的结构动响应问题,建立了Taylor级数方法的理论,给出了递归求解通式.通过对典型方程的求解,阐述了Taylor级数方法的应用.算例表明,Taylor级数方法解决非线性结构动响应问题是行之有效的.     

基于降维算法和Edgeworth级数的结构可靠性分析

针对工程实际中存在功能函数为隐式或高维非线性的复杂结构,本文提出了一种基于降维算法和Edgeworth级数的可靠性分析方法。利用降维算法将n维函数展开为n个一维函数,经变量转换后变量都相互独立且服从均值为0、方差为0.5的正态分布,再结合Gauss-Hermite积分方法计算出一维函数的原点矩,从而得到结构功能函数的中心矩,将所得的矩信息应用到Edgeworth级数展开式中,给出功能函数的累积分布函数表达式,计算得到结构的失效概率。该方法避免了功能函数对变量梯度的要求,仅需少量的确定性重分析计算。数值算例结果表明了本方法的有效性和正确性。      

高维数值积分的蒙特卡罗方法

给出了对于任意概率密度函数产生随机数的一种方法,同时对随机数进行均匀性及独立性检验,将产生的随机数用于计算高维数值积分的蒙特卡罗平均值方法,得到了一种计算高维数值积分的改进平均值方法,并进行复化.最后,给出了几个数值算例以验证方法的有效性.