缓变过渡段几何外形对实验段内气动力特性的影响

本文提出了一种计算高速风洞三维过渡段的气动特性的方法。本方法是以一维管流解为基础,并通过坐标变换,将变系数的线化速度势方程化为波动方程(对超音速流)或拉普拉斯方程(对亚音速流),然后用有限基本解法求解。本文着重讨论了超音速流求解的过程。在超音速流时,对矩形实验段区域,则采用分离变数法求解。作为例子,计算了由圆变方的超音速流过渡段,得出了三维变化对实验段气流参数的影响。     

应用间断伽辽金方法求解二维欧拉方程

在二维非结构网格上,构造了规范正交基作为测试函数,应用不同精度的间断伽辽金方法(DGM)数值求解二维Euler方程,其中数值通量项采用Lax-Friedrichs格式计算.通过对绕NACA0012翼型的跨音速流场的数值模拟,分析了该方法在捕捉间断解方面的特性;同时借鉴有限体积法(FVM)中的一些技术,对求解方案进行了优化.把间断伽辽金方法和有限体积法计算的结果进行比较,认为前者能够更好地处理间断解问题.     

在运动突风作用下亚声速机翼非定常气动力数值计算

本文应用运动涡有限基本解方法求得亚声速机翼在运动突风作用下非定常气动力的数值解。对于不同平面形状的机翼进行了计算。计算结果表明,该方法解类似的问题是有效的。     

三维非结构网格的欧拉方程聚合多重网格法

介绍一种基于三维非结构网格的聚合多重网格技术,用于求解Euler方程定常解,加速收敛效果十分明显。求解Euler方程的基本算法采用有限体积中心差分格点格式。聚合多重网格技术按一定的规则将细网格的若干控制体积聚合在一起,组成一个较粗的网格。数值实验表明本文多重网格法十分有效。     

某些偏微分方程的插值函数微分解法

对某些偏微分方程进行数值解,在规则区域或通过坐标变换能变成规则的区域内,广泛使用有限差分法;在具有不规则边界的区域内,有限元法得到了迅速发展。这里以二维问题为例,用从插值函数导出函数近似微分的方法对方程进行数值解。这实际上是有限差分法的一种推广,又能象有限元法那样适应求解区域和边界的不规则性。     

分析有限大内部开裂板的复变-变分方法及其应用

 本文利用复交函数方法导出了含内部裂纹有限大板的应力及位移的全场解。用变分原理处理边界条件,求解了含中心裂纹、孔边裂纹、偏心裂纹有限大板的应力强度因子,并对用铆接加强环加固以后的孔边开裂板进行了分析和计算。算例表明,本方法收敛快,省机时。     

用有限基本解计算喷气襟翼气动特性

本文介绍使用有限基本解方法对喷气襟翼作气动力分析计算的方法。在DJS-6计算机上对矩形机翼在不同攻角,不同喷气角和不同喷气动量系数情况下作了一系列计算,计算的结果表明,用有限基本解方法来计算喷气襟翼的气动特性是可取的。对二元矩形喷气襟翼的计算结果与国外的某些计算和实验结果作了比较,对三元矩形喷气襟翼的计算是本文对这一问题的一个新的尝试。     

同伦方法中辅助映射的影响

简述了同伦方法中解曲线跟踪算法，它是求解非线性问题的一种非常有效的方法，并说明辅助映射的构造形式将影响非线性问题的求解。根据同伦方法的基本思想和求解过程，提出了采用“相似”准则设计辅助映射的方法，即辅助映射要与原映射同类型，使待求函数较为容易地同化变形到辅助映射，同时求到问题的解。该方法已应用于一些实际课题中。效果很好。     

飞机稳态尾旋平衡点的综合求解方法研究

为求解飞机稳态尾旋平衡点,推导了飞机稳态尾旋的动力学方程,并针对该方程给出了三种求解方法:作图法、数值求解法以及在前两种方法的基础上所提出的一种综合求解方法。综合求解方法首先使用作图法得出方程的所有近似解,并用稳态尾旋判据判断解是否为稳定解,然后将稳定近似解作为数值求解法的初值进行精确求解。运用综合求解方法对算例飞机进行了稳态尾旋平衡点的计算,算法收敛迅速,结果精确,表明该方法具有很好的工程应用价值。     

跟踪问题最优控制律精细积分

构成有限时间最优跟踪系统的控制律需要求解 Riccati微分方程及外部控制输入向量满足的微分方程 ,前者是非线性矩阵微分方程 ,后者是变系数线性微分方程。在结构力学与最优控制的模拟理论基础上所发展的精细积分方法借鉴了计算结构力学中的算法 ,可以精确有效地求解这些微分方程。这种方法的特点之一在于步长幅度变化较大时 ,Riccati微分方程的数值解仍可以保持很高的精度 ,并且变系数线性微分方程的求解亦可纳入其体系而不必用通常的差分方法。本文介绍了用精细积分方法求解这些方程的过程 ,并给出了数值算例。     

注塑模熔体非牛顿流动的数值模拟

研究注塑模流动的数值模拟。用有限元-有限差分方法模拟非牛顿流体的压力场和温度场,用流体网格法（FAN）模拟熔体锋面的推进过程,这种方法具有通用性好和容易实现自动推进的优点。在486微机上对几个算例进行了模拟计算。比较了牛顿和非牛顿流体的流动特点和锋面推进自动处理的效果,得到了满意的结果。     

单列叶栅紊流流场有限分析数值研究

本文将有限分析方法用于曲线座标系上紊流 N-S方程数值计算,研究了高雷诺数时叶栅粘性紊流流场。有限分析方法在网格单元上对非线性偏微分程进行线化处理,在解析边界条件下求出线化偏微分方程的解析解,以此解构造离散方程。有限分析方法能够根据对流的方向和大小自动改变格式系数,并具有数值扩散小、精度高和稳定性好等优点。本文以 k-ε紊流模型模化紊流,以壁面函数方法处理近壁区流动参数。      

N-S方程有限分析数值研究

 本文将有限分析方法用于曲线座标系上紊流N-S方程的数值计算。分别计算了单列和串列叶栅内部流场,计算中采用了k-ε紊流模型和壁面函数。计算结果与试验结果相比较,吻合程度令人满意。有限分析方法在网格单元上对N-S方程进行线化处理,以解析边界条件作为约束,得出解析解,在解析解基础上构造离散代数方程。有限分析方法的最大特点是可以适应对流速度大小和方向,自动调整格式系数,因而具有数值扩散小和稳定性高等优点。     

飞机层合风挡鸟撞击有限元数值模拟

建立了鸟撞击层合风挡三维有限元分析模型,对层合材料采用具有失效模型的率型本构模型,用非线性接触-碰撞有限元程序模拟了飞机层合风挡的鸟撞击动响应过程,给出了层合风挡鸟撞的计算实例,讨论了层合风挡结构被鸟冲击变形破坏的机理,数值模拟结果为层合风挡抗鸟撞击叠层设计提供了理论依据。     

挤压油膜阻尼器-滑动轴承-转子系统非线性动力特性的数值分析及实验研究

 采用非线性模型,研究挤压油膜阻尼器 -滑动轴承 -转子系统的运动稳定性和分岔特性。利用有限元离散变分不等原理和基于线性规划余原理的迭代方法,求解自由边值条件下的油膜力。采用模态综合技术中的缩减系统线性自由度的方法,并提出基于 Newmark法的预估 -校正 -局部迭代的方法,求解系统响应。采用 Poincaré映射法,将周期解的确定转化为 Poincaré截面上不动点的确定,然后分析不动点的稳定性。最后通过实验转子,从实验与数值计算进行对比分析,验证了挤压油膜阻尼器的减振特性及其提高系统运动稳定性的特性,同时也说明了上述方法可以有效地分析这类系统的局部非线性动力特性。     

非对称铺层复合材料板室温形状的非线性有限元分析

由于一般铺层复合材料板壳存在拉伸-弯曲、拉伸-扭转等耦合效应,在温度变化时将出现垂直于中面方向的翘曲变形,给工程使用和理论分析带来了极大的不便。但另一方面,利用非对称铺层复合材料板壳可以有效地减轻结构的重量,并可能提供利用平面模具制造曲面层合壳的新工艺。     

大攻角非线性问题的局部线化方法

 本文引用局部线化的概念,提出一种数值计算方法,用以研究小展弦比、大攻角带前缘或侧缘分离流机翼的双重非线性问题,即基本方程的非线性（包含了速势方程的所有二阶非线性项）和边界条件的非线性（离体涡位置不能预先给定）。该方法适用于亚临界下的所有Mach数,计算简单,工作量小。通过算例计算表明,该方法计算结果与实验值是一致的。     

跨声速翼型绕流计算边界层反方法的研究

本文用有限差分方法,通过有粘/无粘迭代计算了二元翼型的跨声速绕流问题。在边界层粘性区域内考虑了层流、转捩及湍流流动,当边界层内出现分离时,使用边界层反方法,采用代数湍流模型。算例表明,对激波/边界层弱干扰和强干扰情况,该方法的结果与风洞实验结果吻合良好,对于求解边界层小分离流场是一种好的近似方法。     

高阶精度间断Galerkin有限元方法研究

在二维非结构网格上,对高阶精度间断Galerkin有限元方法求解跨音速欧拉方程进行研究。运用间断有限元理论,采用施密特正交化多项式基函数对流场解进行近似描述,使用HLLC近似黎曼解方法计算网格单元边界处的数值通量,积分项通过高斯积分求解,时间推进采用经典四步Runge-Kutta方法,并引入斜率限制器,抑制数值振荡。对NACA0012翼型跨音速无粘流动进行数值模拟,数值结果表明：间断Galerkin有限元方法具有较高的精度,较小的数值耗散和较强的激波捕捉能力。     

计算机翼跨声速压力分布的格林函数方法

本文给出一种用格林定理将跨声速小扰动速势方程化成积分方程式,然后进行数值求解的方法,文中对场元强度进行松弛迭代,一般只需20次左右即达收敛。算例表明,该方法结果与差分法结果一致。