后掠机翼欧拉方程数值模拟

本文对后掠机翼用保角转绘加剪切变换生成C-H型网格,并在这类网格拓扑上,用Jameson的三维欧拉方程有限体积法,四步Runge-Kutta时间推进格式,研制出可供分析三维后掠机翼亚、跨、超声速绕流的计算程序。本方法的特色是改进了机翼表面网格点的分布,使机翼后缘与网格线一致。在数值计算中,通过附面层位移厚度计算,引入了附面层粘性修正,改善了计算结果。     

磁粉探伤用环形试件的几何尺寸——磁粉探伤原理之四十三

根据磁偶极子链近似理论对Betz环的分析结果得悉:D、d、a、h中只有3个是独立的。由上式的数值计算证明:英国的BS 6072-81标准试件设计合理,而日本的JISG 0565-82中B型对比试件的内径过小,无法达到使用目的。本文给出了改进后的合理内径尺寸。     

辐射及化学非平衡流耦合场计算方法

把一种新的三维热辐射计算方法应用到高超声速非平衡流场的数值模拟过程中:在非结构网格上对高超声速化学非平衡流场的数值模拟过程中,耦合了辐射输运方程来求解辐射热流。控制方程为含有化学反应源项和辐射源项的三维N-S方程,数值离散格式采用了Jam eson有限体积法,辐射输运方程采用有限体积法求解。化学反应模型为7组元7反应模型,并采用了化学反应特征时间步长。本文分别对RAM C-Ⅱ和MU SES-C模型进行了计算,计算结果与参考文献的结果吻合。     

一类增量式支持向量机的分析

针对训练数据发生增量改变时,标准一类支持向量机的批处理算法需要重新进行训练,不适合在线增量环境学习的问题,提出一种详细的增量式标准一类分类向量机算法,并通过理论分析对该算法的可行性和有限收敛性进行了证明,确保该算法的每步调整都是可靠的,并确保该算法通过有限步调整最终收敛到问题的最优解。在标准数据集上的实验结果验证了理论分析的正确性。     

解大规模无约束优化的新有限储存对称秩1校正算法

提出了一个新的有限储存对称秩1校正算法,该算法把自适应调比对称秩1校正与有限储存、非单调线搜索技术相结合。算法中的下降方向由有限储存逆自适应调比对称秩1校正产生从而简化计算。将算法与著名的LBFGS进行了数值比较,结果表明,新算法能较好地解一类大规模无约束优化问题。     

用遗传算法精确计算圆度误差

提出了一种应用遗传算法计算满足最小区域法的圆度误差的新思路,并对传统的遗传算法提出了一些改进。采用实数值编码,其计算结果精确度非常高,理论上可以获得全局最优解;保留上一代种群中适应度最好的个体到下一代,可以确保解的收敛性;对基于实数值编码的繁殖算子、交叉算子、变异算子给出了具体的操作方法。仿真结果表明,用改进的遗传算法求解圆度误差,简单明了,收敛速度快,在计算机上容易实现。     

尖前缘任意翼型的超音速和高超音速俯仰安定性导数

本文推广了文[1]的欧拉方程摄动法,讨论了尖前缘任意翼型俯仰安定性导数的计算方法,并应用该方法计算了双圆弧翼型的俯仰安定性导数,其数值结果与文[2]的超音速线化位流方程有限差分法的数值结果比较,更接近于实验值。在激波附体的条件下,本文的方法可用于计算大迎角尖前缘任意翼型的超音速和高超音速俯仰安定性导数。     

爆炸波与超声速运动物体相互作用的数值研究

针对在弹道靶及爆炸波装置上完成的实验,用数值方法来研究在比热为常数和层流情况下爆炸波与超声速飞行物体的相互作用.为了准确模拟爆炸波,首先通过有限体积法求解轴对称N-S方程计算了爆炸波流场,分析比较了数值干涉图和实验干涉图以及数值计算和实验得到的三个截面的密度分布曲线.其次以获得的轴对称流场为爆炸波初场,采用重叠网格技术,通过有限体积法求解三维层流N-S方程数值计算了在比热为常值和层流情况下爆炸波与超声速飞行物体的相互作用,得到了数值干涉图并与实验得到的干涉图进行了比较,分析了爆炸波对飞行物体表面压力的影响.     

钢球碰撞柔性杆动力分析的有效数值方法(英文)

柔性物体的碰撞动力学是工程实际中的一个重要问题。获得问题的动力响应是一个有挑战性的任务,通常通过数值计算来完成。首次采用离散奇异卷积法分析碰撞动力学问题,并通过钢球碰撞柔性杆的算例来研究几种数值方法的有效性。数值方法包括常规的有限元法、离散奇异卷积法、谱有限元法和提出的改进谱有限元法,将各种方法获得的算例的数值结果与理论解比较来研究方法的有效性。研究结果表明,采用恰当的常规有限元长度后,提出的改进谱有限元法给出的接触应力的结果最精确。另外还发现,离散奇异卷积法可以用于求解碰撞问题。     

对特征向量导数之动柔度法补充

一种计算很多特性向量导数的动柔度法在结构处于自由状态（具有刚体运动）下，因其工作方程之系数阵为满阵，使得计算效率不如结构处于约束状态那样好。本文移频动柔度法可以消除原动柔度法的这一缺陷，因此它的工作方程之系数陈总是具有刚度阵那样的带状特点。另外，由于对动柔度采用的移频步骤与特征方程求解时的移频技术是一致的。故而对方法的程序化极为有利。     

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

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

GPS导航定位中周跳探测和修复方法

提出了一种基于双频观测数据的单历元周跳探测和修复方法.该方法分三步实现:(1)将上一历元的模糊度代入当前历元的观测方程,消去模糊度参数;(2)利用拟准检定法检测当前历元是否有异常观测,若有,则判断出发生周跳的卫星;(3)利用LAMBDA方法计算出周跳,与前一历元的模糊度相加,即为当前历元的模糊度.该方法能准确判定出异常观测,即发生周跳的卫星.文中算例与其他方法作了比较,结果验证该方案是可行性的.     

一种快速计算电磁各向异性介质目标RCS的新方法

从电磁理论中的体等效原理出发,推导了任意电磁各向异性介质目标电磁散射的体积分方程。在此基础上,由SWG基函数和伽略金法建立其矩阵方程。在迭代求解过程中,应用了快速计算阻抗元素的等效偶极子法和加速矩阵向量积的快速偶极子法两种快速算法,有效地降低了传统矩量法计算电磁各向异性介质目标RCS的时间。同时,在快速偶极子法中,远场组间互阻抗元素无需显式计算,并且形式简单的聚集、转移和发散函数可现用现算,从而大量节省了计算机内存。数值仿真结果表明:本文使用的方法在分析电磁各向异性介质目标的电磁散射特性中是非常有效的,不仅有高的计算效率和低的内存需求,而且有良好的数值精度。     

基于变马赫数的五孔探针三维插值方法

探讨了五孔探针气动数据的插值方法,为提高插值精度,开发了基于传统线性插值法的三维线性插值法。该方法把同一探针在不同马赫数下的校准图形成三维数据库,将实验数据通过三维图进行插值。并使用改进前后的两种插值方法分别对校准风洞测得的数据进行整理,对比结果证明:在实验工况连续变化的情况下,三维线性插值法在插值精度上要优于传统线性插值法。特别是在来流马赫数变化较大时,该方法可以改善单一校准文件处理造成的数据精度问题,可用于自动化流场采集系统,为流场高速高精度采集奠定基础。     

应力极值分布、强度威布尔分布的可靠性

本文运用应力—强度干涉理论,推导了应力为Ⅰ型极小值分布,强度为威布尔分布的可靠度计算公式,并对冗长的计算公式进行简化,在简化公式的基础上,运用一定的数学技巧,改变积分公式中的积分变量和上下限。将被积函数化成在某一区域内的可积函数。采用de Boor编制的一种严谨的自适应Romberg外推格式的FORTRAN程序进行数值积分。对应予不同的组合参数,给出应力服从Ⅰ型极小值分布,强度服从威布尔分布的可靠度数值。本文最后讨论了服从这两种分布的组合参数的变化对可算度数值变化的影响。     

标准孔板流量计内部流场的CFD数值模拟

通过在标准孔板流量计中引入CFD数值模拟,为流出系数的获取提供了新途径.对不可压缩流体在不同流量、不同直径比、不同孔板轴向厚度和不同流动介质下的内部流场进行了数值模拟计算,并将计算出的流出系数与根据ISO公式计算出的流出系数进行了分析对比.结果表明,随着结构参数和工作条件的改变,流出系数都会随之发生变化,但ISO公式对孔板厚度的变化不太敏感.CFD数值模拟可以作为标准孔板流量计的辅助设计与标定手段,以进一步提高孔板流量计计量的准确性.     

基于MPP编程环境的一种并行子空间迭代法

子空间迭代法是科学与工程计算中求解广义特征值问题的有效方法 ,针对向量机和共享内存的多处理机 ,前人已成功地作了并行处理。文中给出了适合 MPP大规模并行计算机的并行子空间迭代法。该算法将广义特征值问题转换为一般特征值问题 ,其计算工作量主要体现在矩阵乘法 ,通过对该方法作并行处理 ,使矩阵求逆及一部分乘法运算转换为各结点机上三角形方程组的并行求解。在大规模并行计算机 PA R95上结合 J8- II机翼的动力特性问题对该算法作了数值试验 ,结果说明所给算法是非常有效的     

微分求积单元法分析变曲率井内受径向约束管柱的非线性屈曲

基于Love的空间弯扭杆平衡方程,通过引入井壁约束条件导出了平面变曲率井内受径向约束管柱的平衡方程。采用微分求积单元和增量迭代法求解了曲率线性变化的井内管柱的非线性屈曲问题,通过与有限元计算结果的对比验证了所构建方法和编写程序的正确性,而且还表明微分求积单元法有方法简单、易于实施,计算量少、精度较高等优点。计算结果表明,等曲率井内管柱屈曲的临界载荷明显大于曲率线性变化井内管柱屈曲的临界载荷,另外,变曲率井眼的曲率变化对管柱弯矩、井壁约束力有显著的影响。     

基于旋转坐标系的高阶间断有限元方法非定常湍流数值模拟

基于高阶间断有限元方法(Discontinuous Galerkin method,DGM),对旋转非惯性系下耦合了修正的一方程S-A模型的RANS方程进行了离散求解。为了在稀疏网格上获得更贴近真实的物面形状,使用了多层高阶弯曲网格方法对物面进行拟合。非定常时间推进采用了隐式双时间步方法,每个时间步产生的线性系统采用预处理的方法,即广义最小残差方法(Generalized minimal residual method,GMRES)来求解。计算了旋转圆柱绕流以及经典翼型振荡算例的升力和力矩迟滞曲线,与实验结果以及前人的计算结果对比验证了本文方法的正确性和有效性。