一种跨声速定常流场求解加速方法

跨声速定常流场的隐式求解相当于使用牛顿迭代法求解一个非线性方程组。为满足牛顿迭代收敛性的要求,通常需要对所求解问题进行全局化处理。在同伦延拓的框架内,提出了一种基于拉普拉斯算子的方程延拓方法,提高了定常流场隐式求解收敛速度。针对定常流场通常初始化为均匀来流的特点,一方面利用拉普拉斯算子的椭圆性加快边界条件信息向流场内部的传播,另一方面利用拉普拉斯算子的线性和正定性改善延拓问题的正则性,综合两者增加拟牛顿算法的稳定性,提高可用CFL数,最终达到提高流场求解效率的目的。由于流场问题的复杂性和非线性,难以通过理论分析得出先验的最优非线性求解策略。因此,通过无黏NACA0012翼型、湍流RAE2822翼型和三维ONERA M6机翼等算例的数值实验,研究了拉普拉斯项参数对收敛效率的影响,给出了效率较优的参数组合,验证了本文方法在跨声速情况下相对于经典伪时间推进法可以节约20%以上的CPU计算时间。      

有热辐射边界条件的热传导问题

本文研究在有热辐射的非线性边界条件下,在圆柱体内的温度分布问题。本问题可以化为非线性积分方程问题,并最终地利用非线性代数方程组求解,从而得到了温度场的级数形式的解析解。     

嵌入式大气数据传感系统测压孔布局研究CSCD

嵌人式大气数据传感系统是一种依靠嵌入在飞行器前端的压力传感器阵列来测量飞行器表面的压力分布,并由此压力分布通过求解非线性方程组间接获得飞行参数的先进飞行数据传感系统。针对嵌入式大气数据传感系统测压孔布局进行研究,通过三点法对各种布局的数据解算收敛性进行了分析,并在白噪声干扰下对各种布局下的解算精度进行了验证。     

Euler方程的分裂型通量分裂双时间步隐式方法

传统隐式方法有格式复杂、计算量大等缺点,在Euler方程的差分离散过程中,利用算子分裂思想,结合通量分裂法、双时间步法等隐式离散方法,构造了一种更简单的分裂型隐式计算方法.通过对典型空气动力学问题的计算,检验了该方法的有效性和可靠性,并对其性能做了具体讨论.该方法具有稳定性好、时间步长约束小等隐式格式的普遍优点,同时具有格式简单、程序易实现等优点;避免了传统隐式方法单步推进时的方程组常规求解及矩阵求逆过程,计算量小;比LU-SGS方法收敛速度快.      

一种功能梯度Timoshenko梁的损伤识别方法

为了获得功能梯度材料的高精度损伤识别方法,本文基于动力学方法,通过对状态空间变量进行变量替换,求得了沿轴向指数分布的功能梯度Timoshenko梁的传递矩阵,通过分析裂纹对结构局部柔度的影响,采用扭转弹簧模拟裂纹对结构局部柔度的贡献,建立了功能梯度Timoshenko梁的表面裂纹传递矩阵,并且推导了复杂边界条件下多跨梁的理论模型。通过将非线性方程组转化为单一目标函数优化问题,并将增广拉格朗日算法与差分进化算法相结合对结构进行损伤识别。计算实例表明,本文提出的方法具有精度高、收敛快等特点,且适用于复杂边界条件下多损伤模型的损伤识别。      

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

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

基于润滑理论的二维积冰数值模拟

将基于润滑理论的水膜-冰层模型推广到在任意二维截面上建立的正交曲线坐标系,得到描述二维型线表面水膜流动和积冰的偏微分控制方程.采用一种隐式-显式有限差分格式建立控制方程组的代数形式,给出求解非稳态控制方程组的数值方法.为验证模型和数值解法的有效性,在典型的航空积冰和传输线积冰环境下对翼型和传输线表面的积冰算例进行数值模拟.将水膜-冰层模型的冰形计算结果与传统Messinger模型的模拟结果以及冰风洞试验结果进行对比,低温明冰条件下采用当前方法计算的翼型表面冰形接近于Messinger模型的模拟冰形曲线;相对高温条件下的计算结果比传统Messinger模型更为精确.低速积冰环境下Messinger模型难以模拟的传输线表面积冰同样可以采用水膜-冰层模型进行有效 预测.     

刚性方程的二阶Gear方法的注记

构造了一个求解刚性方程的二阶差分格式,它是隐式的,并采用简单迭代求解,对于稳定系统,迭代收敛性条件与步长无关,而对不稳定系统,迭代收敛性条件对步长却有所限制,利用这一限制性条件去克服解刚性方程的二阶Gear方法时出现的“危险性”。     

球坐标系六片网格下三维定态行星际太阳风模拟

采用二阶MacCormack差分格式, 利用稳态的磁流体(MHD)方程组在球坐标系六片网格下模拟研究了行星际太阳风. 六片网格系统能有效避免极区奇性和网格收敛性. 迭代按径向方向推进求解, 很大程度上减少了计算量, 节约了计算时间. 内边界条件根据太阳与行星际观测确定, 比较测试了5种内边界条件, 模拟给出了1922卡林顿周的背景太阳风结构. 几种内边界条件所得模拟结果与行星际观测基本吻合. 太阳风速度采用McGregor 等的经验公式给出, 磁场由水平电流片(HCCS)模型得到, 密度和温度分别根据动量守恒和气压守恒得到, 研究表明采用这样的边界条件模拟结果最佳.      

复合材料层板在弹性支持下的非线性稳定性

利用能量变分原理和非线性几何方程，建立具有弹性约束的复合材料层板，在面内载荷作用下的非线性稳定性控制方程组，通过运用广义傅里叶级数法对控制方程组进行求解，得到相应的截荷－中心挠度曲线。计算结果表明：非对称层合板在面内载荷作用下，一般产生的是向变形，但在一定的铺层顺序，边界条件、受载方式等情况下，也可能发生分叉型屈曲。     

月球探测器软着陆机构展开过程的运动学分析

软着陆机构是月球探测器实现软着陆的关键部件。针对一种四腿式软着陆机构方案,论述了该软着陆机构展开的工作原理,分析了机构的奇异性,推导了机构不发生奇异的条件,建立了展开机构的运动学方程组,采用Newton-Raphson法求解了该非线性运动学方程组,获得了主支撑杆运动参数的变化规律。     

吴方法在平面并联机构位置正解中的应用

采用吴方法对平面并联机构位置正解问题进行了研究.吴方法是一种求解非线性方程组的数学机械化方法,采用这种方法任何非线性方程组都可以在有限步内得到解决.在给出了吴方法基本原理的基础上,对本问题进行了求解,并将原始方程组转化成为一个三角化的方程组.其中单变量方程的次数为6次,说明平面并联机构可以有6个不同的位姿.最后用数值实例进行了验证,给出了计算结果.吴方法在这一问题中的应用,为求解其它机构学难题提供了新途径.     

三维重力波非线性传播数值模型的一种求解格式

采用交替方向隐格式（ADI）求解三维可压大气中重力波非线性传播的数值模型，经过几个波周期后，重力波波包括估然沿着线性理论射线方 稳定地传播，波包保持着很好的形状，内点和边界保持稳定；水平风与垂直运动的振幅随着高度的增加而增加。这说明ADI格式是模拟三维可压大气中重力波传播的一种可靠的格式，模拟结果表明，在非线性情况下，重力波仍可保持性传播的某些特征，非线性效应有现在对波包括传播速度的修正。     

轨道转移计算中的逐次消元迭代法

本文针对导弹飞行轨道转移的一个算例,叙述应用消元法思想构造数值求解非线性方程组的迭代算法的一般原理,具体算法步骤和一些计算结果。给出算法成功的理论依据,最后给出实现算法时的一些变形。     

曲线纤维壁板屈曲/后屈曲建模与快速分析方法

曲线（VAT）纤维复合材料壁板相比广泛应用的直线纤维形式具有更优的面内稳定性,作为机翼壁板,在同等质量时具有更高的抗屈曲潜力。为深入研究纤维路径对于曲线纤维壁板稳定性的影响规律,从各向同性薄板的理论出发,推导曲线纤维壁板在面内载荷下的稳定性分析方法;通过Airy应力函数和拉格朗日乘子描述边界条件,建立曲线纤维壁板适用于任意位移及载荷边界条件的单一变分方程,避免了非线性平衡方程和非线性相容方程间由于反复迭代对求解速度的制约。基于冯卡门大变形方程发展了曲线纤维壁板后屈曲状态下的非线性稳定性问题求解模型,并采用瑞利-里兹法建立了屈曲/后屈曲一体化半解析快速求解框架,该框架的求解精度与商用软件MSC.Nastran一致,但求解时间远低于商业软件;利用此优势,可以快速分析给定任意位移边界条件下的曲线纤维壁板屈曲响应特性,并得到纤维路径的影响规律。     

特定消谐逆变器开关角度的复形调优解法

针对特定消谐式逆变器开关角度难于求解的问题,提出了一种基于约束条件下n维极值复形调优法的开关角度求解方法,分析了目标函数和约束条件的设计方法.该算法对开关角度初值的选取只需考虑约束条件,大大增加了开关角度初值的取值范围;求解过程中不存在解发散或奇异无解问题,与迭代法解非线性超越方程组相比可靠性更高;通过约束条件的限制,使求解过程自动消除了可能出现的窄脉冲;通过目标函数和加权因子可灵活的进行开关角度优化.算例和实验结果验证了方法的正确性和有效性.      

超声速进气道可压及不可压流动数值模拟

采用Chio-Merkle预处理矩阵对可压NS(Navier-Stokes)方程组进行时间预处理,分析了预处理方法的物理和数学背景.用有限体积方法,结合LU-SGS(Lower-Upper-Symmetric-Gauss-Seidel)隐式时间积分和AUSM+(P)(Precondition Advection Upstream Splitting Method)格式、中心差分2种空间离散格式,求解该预处理NS方程组.通过圆弧凸包流动和二维方腔顶盖驱动流动的数值试验表明,该方法克服了传统时间迭代方法模拟低速流动时的刚性问题,加速了收敛过程,可以同时有效地模拟可压及不可压流动.用此方法与块结构化网格技术相结合,进行了在不同飞行马赫数、攻角和出口反压条件下,混压式轴对称超声速进气道三维流场数值模拟.计算结果表明:该方法能准确地捕获复杂激波系,清楚地揭示进气道三维流动现象;进气道性能随工作条件的变化规律,与理论分析一致.      

关于求解多元非线性方程组的(n 1)点顺序正割法的局部收敛性

本文提出了求解多元非线性方程组的(n 1)点顺序正割法的局部收敛条件,并且给出了收敛性的一种证明方法。