基于非确定性的翼型鲁棒反设计方法

实际的翼型是在一个非确定性的状态范围内工作的,为此应用鲁棒控制理论建立了翼型的鲁棒气动反设计方法及数值求解过程,基本的优化器为基于伴随方程的位流优化方法。对守恒型位流方程应用高效AF 2格式求解,对非守恒型的伴随方程采用了旋转差分和高效AF 3格式,提高了梯度计算的精度和效率。优化算例表明:本文方法所得到的翼型在一个较宽的工作状态范围内都具有良好的性能。     

中心或反中心对称线性方程组的缩减算法

推广了文［１，２］的概念，系统地叙述了中心对称和反中心对称矩阵的性质，导出了一类求解系数矩阵为中心对称和反中心对称线性方程组及超定方程组的缩减算法。这类算法可节省７０％－７５％的计算量，大大提高了计算效率     

求解二维双曲型守恒律的一类MmB差分格式

本文利用ＭｍＢ差分格式的思想，将文（１）中所得差分格式改进为ＭｍＢ差分格式，对步长比为ＣＦＬ限制虽然较文（１）严格，但有效地防止了非物理振荡，而且在数值计算过程中，ＣＦＬ可适当地放大，数值实验结果表明，改进的ＭｍＢ差分格式处理间断解的能力令人满意。     

评述B—P统一弹—粘塑性本构模型

综述了Ｂｏｄｎｅｒ－Ｐａｒｔｏｍ统一弹－粘塑性本构模型的基本理论，包括流动法则、运动方程和内变量演变方程三个部分。该模型是基于位错动力学的思想提出来的，认为材料变形的应变率由弹性和非弹性两部分组成，它们始终耦合在一起，采用内变量来表示材料抵抗塑性变形的非弹性状态，内变量包括各向同性硬化和方向性硬化两种。该模型不需要屈服准则和加卸载准则，具有较强的模拟能力，能模拟循环硬化与循环软化、非比例加载引起的     

双马赫反射的数值模拟与光学干涉定量测量的比较(M_s=4.62)

本文用欧拉方程计算了楔形角为40°,M_s=4.62双马赫反射的流场,计算采用了矢通量分裂和 TVD 两种格式,并考虑了γ=1.4和高温平衡气体两种情况。流场的密度场测量采用双曝光全息、M-Z、差分三种干涉技术。壁面密度分布的计算与测量结果比较表明:三种干涉方法测量结果差别不大,都处在γ=1.4和平衡气体的计算结果范围内。     

一类奇异H∞问题的降阶控制器的解释

一类奇异H_∞控制问题可以利用系统矩阵的结构特征通过线性矩阵不等式或2代数黎卡提方程方法直接显式构造得到降阶控制器.进而利用系统不稳定的不变零点可以构造出更低阶次的控制器.基于2代数黎卡提方程方法,特别是系统Hamilton矩阵的结构特点和系统可观测性理论,论述了不变零点在构造系统降阶控制器中所起的作用,即通过选择控制器的自由参量,使某些不稳定的不变零点关于虚轴对称的点构成系统控制器不可观的稳定极点,从而实现了控制器进一步降阶.     

翼型周线布线性涡的有限基本解法

在翼型周界线上布线性变强度分布涡,在一定边界条件下建立求解涡强度的线性方程组。经过试探,确定了比较合适的控制点位置。为满足库塔条件而令后缘处涡强度为零,但导致了方程组的矛盾性。采用最小二乘法求解矛盾方程组,可获得误差最小的唯一解。本方法简便,易于应用。计算结果与试验值或精确解相比,有良好的准确性,可用于一般工程计算。     

任意形状壳体非线性非K－L普适理论

从理性力学的观点导出任意形状壳体非线性非Ｋ－Ｌ理论的普适方程，由此即得Ｋ－Ｌ理论的普适方程，以及各种近似方程，并可将现有文献上的方程作为其近似特例而包含在其中．     

用泰勒级数求解非线性代数和微分方程组

提出一种求解非线性代数方程和非线性常微分方程的新方法。所论方程组被转换成增量形式，而解被表示成泰勒级数。在代数方程（多项式形式）的情况下，此算法归结为求解一系列线性方程组，而且在每一增量步中系数矩阵是不变的；在常微分方程组的情况下，此算法归结为一组递归算式，并不要求解方程组。此方法固有的自动走步功能，可保证得到收敛的解，从而大大节省计算时间。     

高超声速验证飞行器助推分离段流场数值研究

为了模拟有相对运动的多体非定常绕流，发展了基于弹簧近似的非结构动网格方法及耦合动网格的Euler解算器，作为方法验证，模拟了俯爷振动的NACA0012翼型绕流，计算结果与实验及文献结果非常接近。采用二维近似外形，对高超声速验证习行器助推分离段流场进行了数值研究，得到了不同时刻由于多体相对运动形成的干扰流场结构以及分离过程的气动力参数。