柔性多体撞击系统的精细H∞控制研究

本文在柔性多体撞击系统动力学方程的基础上 ,将撞击看成是一种扰动 ,设计了一个稳定的控制器以使撞击扰动对控制输出的影响达到最小的目的 ,在分析过程中采用 H∞ 控制原理对该问题进行了研究。在计算中 ,将精细积分法应用于状态动力学方程和 Riccati方程的求解中 ,使计算精度得到大幅度提高。最后通过数值例题说明了本文方法的有效性     

基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。     

细长三角翼滚转/侧滑耦合运动效应分析

通过耦合求解NS方程组和刚体动力学方程组,数值研究了80°后掠三角翼自由滚转/自由侧滑的耦合运动特性。通过与单自由度翼摇滚特性比较,发现滚转/侧滑耦合条件下三角翼的摇滚现象更为剧烈:在相同的滚转角下,耦合运动的滚转角速度更大,具有相同角速度时,耦合运动所达到的滚转角更大。研究表明:滚转/侧滑耦合条件下,涡流的非对称迟滞振荡仍是维持三角翼周期性等幅摇滚振荡的流动机理;摇滚过程中法向力的侧向分力是产生侧滑运动的主要原因;左滚右侧滑、右滚左侧滑是三角翼滚转/侧滑双自由度耦合运动的作动机制;侧滑运动的引入,导致迎角和侧滑角随滚转角迟滞变化,强化了三角翼摇滚过程中的迟滞效应,加剧流动结构的非对称特性。     

三维可压缩 Navier-Stokes 方程的间断Galerkin 有限元方法研究

拓展了二维间断 Galerkin(DG)有限元方法研究,将该数值方法用于三维可压缩欧拉方程和 Navier-Stokes方程的求解。基于六面体网格单元,采用插值方法将物面的四边形面网格单元构造为弯曲面网格单元,更好地表述了真实物面特征;物面边界相邻体网格单元相应构造为高阶体网格单元,其余体网格单元采用八节点六面体单元,以较小的计算代价使网格满足 DG 方法计算需求。通过对三维带 bump 管道内流、圆球绕流以及旋转流线体绕流进行的数值求解,验证了边界弯曲方法的可行性及 DG 方法的高精度特性。此外,由于采用了隐式计算方法,仅需较少的时间步就能迭代收敛。     

直升机缆位系统建模与仿真

针对直升机绞车收放时缆位保持问题,建立了负载出水回收进舱阶段缆位的动力学模型。由于浮力与洋流作用突然消失,负载出水后发生摆动,将负载的摆动作为带有理想约束的质点系动力学系统,采用拉格朗日方程建立动力学模型。考虑直升机悬停提升的情况进行仿真,仿真结果表明建立的数学模型能够较好的模拟负载摆动时的缆位变化。     

燃烧室污染物生成大涡模拟模型的研究

采用亚网格动力学模型和联合概率密度函数(PDF)亚网格模型对污染物NO的生成进行大涡模拟.计算中采用的数学模型有:k方程亚网格尺度模型来估算亚网格湍流粘性;亚网格EBU燃烧模型估算化学反应速率;在交错网格体系下, 采用SIMPLE算法和混合差分格式来求解差分方程.数值分析两种不同亚网格污染物生成模型和两种进口气流温度对NO排放的影响, 计算结果与实验数据比较表明:两种亚网格污染物生成模型都能较好的预估燃烧室污染物NO的生成, 但亚网格动力学模型的计算工作量要比联合概率密度函数亚网格模型少得多, 更适合工程应用.      

一类高阶有理差分动力系统

讨论了高阶有理差分方程xn 1=(αxn-k 1)/(-1 xnxn-1…xn-k 1),for n=0,1,…,其中α<0,k∈{1,2,…},x-k 1,…,x0都是实数,并且使得x-k 1…x0≠1。给出了这类方程解的明确表达式,并对其解的稳定性和收敛性也做了进一步讨论。     

新型航空润滑油油膜拖动力计算研究

在自制的试验装置上测量了新型航空润滑油油膜拖动力;给出了高剪切率弹流润滑下的Ery-ing特征应力、极限剪切应力、剪切弹性模量的新型航空油流变参数拟合计算式;确定了该油本构方程;建立了该油油膜拖动力计算式, 并进行了试验验证.试验结果表明:该油在试验条件下表现为粘弹性;拖动力的理论计算值与试验测量值相当吻合.给出的拖动力计算方法可用于指导润滑分析领域的工程计算.      

内埋式弹舱舱门气动载荷计算分析研究

采用中心有限体积法进行空间离散和四步Runge-Kutta显式推进方法进行时间推进数值求解了非定常欧拉方程;采用Delaunay方法生成内埋式弹舱舱盖打开时的非结构网格.数值模拟了内埋式轰炸机弹舱单开舱门在分别打开8°和45°时的复杂流场.通过对上述复杂流场的分析,发现作用于舱门的气动载荷具有使其趋向于关闭的特性,且有随打开角度的增大而减小的趋势.     

陶瓷纤维梯度增强活塞的梯度方程研究

 功能梯度材料零件具有单质材料零件无法比拟的理化性能优势,然而由于材料分布复杂以及对功能梯度材料本身性能研究不充分,使性能分析存在很多困难。论文应用复合材料热性能理论,采用有限元分析软件ADINA,分析陶瓷纤维梯度增强活塞（材料梯度方程的参数不同）的温度分布和应力分布,结果表明陶瓷纤维梯度层可以明显改变活塞温度分布,缓和由于热膨胀系数不匹配,在陶瓷纤维增强层与活塞本体交界处产生的应力。根据计算结果拟合出温度峰值、整体应力峰值和层间应力峰值与方程系数之间的曲线,并加以验证。