民用飞机水上迫降分析模型和数值仿真

模拟计算了飞机水上迫降的冲击载荷.利用显式积分法求解离散的扭格朗日有限元方程,并采用欧拉有限体积法求解欧拉控制方程,通过一般耦舍(General coupling)算法实现流固耦合计算,运用Dytran进行数值仿真计算.飞机结构入水时压力在初期就达到峰值,然后衰减,在峰值过后会出现小幅渡动.相同条件下弹性体材料模型的压力峰值要比刚体小.4种工况下飞机着水时的压力峰值中以5°攻角最小,随着仰角增大,压力峰值以及峰值出现的区域随之改变.     

二维非结构网格的生成及其Euler方程解

 采用推进阵面法生成三角形非结构网格,用多重网格技术求解泊桑方程以取得结构背景网格.提出了一种高效的边界剖分法,并用它来生成初始阵面.另外,还提出了一种新颖的数据结构——数组链接表,用它来存贮网格数据.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼型和多段翼型的算例,并与实验结果进行了比较和分析.     

PWM变频器供电的感应电机传热特性

为了探究电动机在脉宽调(Pulse width modulation,PWM)变频器供电下的传热特性,近一步揭示变频参数（调制比）对电动机温升分 布的影响特征,以一台采用PWM变频器供电下55 kW驱动用感应电机为例,基于流体力学及传热学基本原理,结合电机通风结构特征,建立外部包裹有空气域的三维流热耦合求解域模型,并采用有限体积元法对电机内的温度场进行了数值研究。此外,针对PWM不同调制比控制条件下电动机全域内的传热特性进行了对比分析,结果表明：两个不同调制比控制条件下电机求解域 内各主要部件温升分布趋势大致相同;径向上,转子部分温升较高,在气隙位置温升出现阶跃式变化,定子部分温升较低;轴向上,电机各部分近风端温升较低,远风端温升较高;周向上,定子轭部温升呈波浪式变化。     

可控流场尺度预混湍流燃烧器及其火焰结构分析

为了研究单一湍流场参数对预混湍流火焰结构的影响,以及拓宽湍流场的强度和尺度范围,发展了一套可变结构的预混湍流燃烧器。采用恒温型热线风速仪标定流场,得到了一系列湍流参数。流场标定结果表明:该燃烧器能显著拓宽湍流强度和尺度范围,并能利用不同几何结构产生多种可控流场,实现研究单一湍流参数对湍流燃烧速度和火焰结构影响的目的。选用有代表性的15种湍流孔板组合结构,利用OH-PLIF燃烧激光诊断技术,开展了湍流燃烧实验,结果表明:湍流强度的增大（1 < u'/S_L,0 < 10）使得湍流火焰分区扩展到了薄层反应区,火焰面破碎程度明显增强,孤岛结构明显增多。高宏观雷诺数下,积分尺度的增长对湍流燃烧速度起抑制作用,可能存在临界宏观雷诺数Re_c,能够表现流体惯性力占主导地位的程度,决定积分尺度对湍流燃烧速度的影响效果。积分尺度能量大,扰动能力强,故积分尺度越大,火焰体积越大;但过高的湍流强度会使火焰面褶皱更加剧烈,小尺度叠加在大尺度上的程度增强,最终也使火焰体积显著增大,掩盖了积分尺度对火焰体积的影响,说明积分尺度（表征大尺度）不如湍流强度（表征叠加小尺度的程度）对火焰放热率影响大。     

二维非结构网格的一个TVD型有限体积方法

考虑双曲型守恒律方程,对二维非结构三角形网格给出一种 TVD型有限体积方法,主要思想是在一阶单调格式的基础上,在每一个单元上对变量作单调线性重构函数,时间离散采用二阶 Runge-Kutta方法。通过计算分析了该方法的精度,对平面激波反射和空穴流动的计算结果表明该方法是成功的。     

结构/直角切割网格下时域有限体积法在计算电磁中的应用研究

 采用有限体积法直接求解时域麦克斯韦方程组，研究了二维圆柱、翼型和三维球体、立方体、双球体等典型完全导电目标的电磁散射问题。空间离散使用近似黎曼解构建网格通量，并针对二维结构网格、三维自适应直角切割网格分别给出重构方法，时间方向采用二步龙格库塔法。计算与理论、实验或文献结果符合良好。     

2.5维机织复合材料疲劳寿命预测方法

针对疲劳载荷作用下的2.5维机织复合材料,建立了疲劳寿命预测方法.该方法主要包括单胞模型、疲劳失效判定准则和材料性能退化方法3部分.选取单胞模型为研究对象,利用三维有限元技术进行应力分析;引入改进的三维Hashin疲劳失效准则和Mises准则作为纤维束和树脂基体的疲劳失效判据;采用刚度性能突降准则描述疲劳失效后的材料性能,采用考虑纤维体积分数影响的剩余刚度和剩余强度退化模型描述失效前材料的性能.通过疲劳寿命预测值与试验值的对比,验证了疲劳寿命预测方法的有效性.研究表明:经向拉-拉疲劳寿命随经纱纤维体积分数增大而增加,纬向拉-拉疲劳寿命受纬纱纤维体积分数影响较小.      

Determination of thermal expansion coefficients for unidirectional fiber-reinforced composites

In the present work, the coefficients of thermal expansion（CTEs） of unidirectional（UD）fiber-reinforced composites are studied. First, an attempt is made to propose a model to predict both longitudinal and transverse CTEs of UD composites by means of thermo-elastic mechanics analysis. The proposed model is supposed to be a concentric cylinder with a transversely isotropic fiber embedded in an isotropic matrix, and it is subjected to a uniform temperature change. Then a concise and explicit formula is offered for each CTE. Finally, some finite element（FE） models are created by a finite element program MSC. Patran according to different material systems and fiber volume fractions. In addition, the available experimental data and results of other analytical solutions of CTEs are presented. Comparisons are made among the results of the cylinder model,the finite element method（FEM）, experiments, and other solutions, which show that the predicted CTEs by the new model are in good agreement with the experimental data. In particular, transverse CTEs generally offer better agreements than those predicted by most of other solutions.     

有限体积法的分析与改进

本文对以Jameson为代表提出的常用的一种有限体积法进行了较深入的分析，指出了它的适用范围，网格对计算精度的影响并提出了改进措施。算例表明本文提出的分析与改进措施是正确、有效的。     

含运动边界二维多体干扰流场的数值模拟

本文采用二阶精度Godonov扩展方法求解非定常：Euler方程，以基于格心的自适应非结构网格的显式有限体积法为基础，空间和时间都是二阶精度，使用HLLC近似黎曼解的方法计算网格单元边界处的守恒量通量，对作俯仰振动的NACA0012翼型绕流及激波管中运动、静止圆柱的流场进行了数值模拟，并对数值模拟的结果进行了分析，为进一步对含多体分离流场进行数值模拟打下了基础。