双曲型守恒律方程的一个大时间步长二阶TVD差分格式

本文根据Harten,A.的大时间步长差分分格式构造思想,为双曲型守恒律方程弱解计算构造了一个2K 3点大时间步长二阶显式差分格式——LTS-LF格式,得到了其在CFL限制K下为总变差不增差分格式(TVD格式)。文章按照Roe的方法推广格式到方程组情形,并就Burger’s方程和Euler方程组黎曼问题进行数值试验,结果令人满意。     

液体非线性晃动试验与等效力学模型仿真

航天器的在轨机动可能导致液体燃料出现包括非线性晃动在内的复杂动力学行为。采用运动脉动球模型（MPBM）对液体非线性旋转晃动问题进行了研究,并通过开展地面试验验证了MPBM仿真结果的准确性。首先,在MPBM的建模过程中以晃动回复力的形式考虑重力的影响,具体是通过虚功原理将重力作用引入到脉动球的能量变化过程。然后,搭建地面晃动试验平台,实现了对球型贮箱晃动现象的准确观测以及对晃动力和晃动力矩的实时测量。最后,以试验实测的加速度激励和晃动力为参照,采用MPBM进行了仿真与对比,仿真结果较好地匹配了试验结果。文中对液体晃动的合理抽象进一步地提升了等效力学模型对复杂晃动行为的仿真能力,为实现航天器更为高效、准确的在轨控制提供基础。     

翻管数值参数分析法研究

采用塑性流动理论，薄壳无力矩理论及幂函数强化模型，对求解轴对称塑性平面应力问题的数值参数法作了进一步的研究，推导出成形力学分析的一阶微分方程组，在此基础上，将该微分方程相应用于薄壁圆管轴压外翻成形的稳态流动分析，计算时采用了四阶龙格库塔法，并获得了应力，应数值解，本文通过大量的计算，分析了边界条件对计算结果的影响，并将理论计算结果与实验作了比较。研究结果表明，轴对称塑性平面应力问题可以转换为一阶常     

关节铰柔性多体系统动力学研究

针对关节铰柔性多体系统建立了Ｒ／Ｗ体系的动力学方程，并进一步推广到非树及约束型多体系统，本文采用适合描述弹性变形和进行弹性体离散的刚性连体浮动坐标系，用相对描述方法来分解部件的运动，用部件模态方法对弹性体进行离散，从达朗伯原理出发建立了Ｒ／Ｗ体系的系统动力学方程。方程对于运动的分解，部件的弹性离散以及对系统的结构描述都比较直观和简单，系统的广义坐标数目也小。方程推广到非树系统和约束系统后，采用奇异     

任意回转面叶栅C网格的数值生成

首先尝试采用任意曲面椭圆方程，数值生成任意回转面上叶栅Ｃ网格，经过对边界条件的仔细处理，生成的网格可保证在周期外边界和重合内边界上网格线导数连续，在叶片表面网格线相互垂直，提出了结合初始网格线分布，控制环绕叶片网格线分布的方法。     

二维突扩管流的有限分析解

运用有限分析法求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，获得了二维突扩管流的数值解。对于雷诺数为２５，５０，１００，１５０，２００的来流，预报的流动图谱具有相似性，回流区域及再附距离的计算也获得了比较满意的结果     

用 Euler 方程计算某些特殊的二维分离流动

采用非定常Ｅｕｌｅｒ方程、总体时间步长、隐式时间推进，计算了几种特殊类型钝体的二维有分离绕流。一类是光滑的圆柱体，在跨音速时计算出因激波诱导分离而周期性地脱落在尾迹上的非对称涡排。另一类是带尖角的物体，如三角形柱体，计算出从尖角处脱落在尾迹上的分离旋涡，形成卡门涡街。对这两类物体，计算的平均阻力系数，分离涡脱落的Ｓｔｒｏｕｈａｌ数以及流态都与实验接近。结果表明，对于这类绕流问题，用非定常Ｅｕｌｅｒ方程计算，具有一定工程实用意义。     

基于分布式粒子群算法的翼型优化设计

采用求解N-S方程作为优化算法中的CFD分析方法,基于标准粒子群优化算法(PSO),将其与遗传算法中的选择机制相结合,形成了一种改进的基于自然选择的粒子群算法(SELPSO),以提高算法的求解精度和改善算法的全局收敛性。为改善串行粒子群算法效率低,耗机时等缺点,文中将分布式计算引入到优化设计过程中,实现了基于分布式粒子群算法的翼型设计优化系统,设计实践表明,文中发展的优化算法对优化设计系统质量和效率都有着大幅度的提高,在工程中具有很好的实用价值。     

泛函级数在导弹主翼振动分析上的应用

为了研究突风载荷对导弹主翼的垂直振动及俯仰倾角的影响,提出了对非线性动态系统做建模研究的Voherra泛函方法的任意高阶核估计方法.该方法在核理论基础上,构造线性空间,将求解Volterra泛函各阶核的问题转化为求输出观测向量在希尔伯特空间子空间上的分量,利用线性空间中向量内积的求解而间接辨识出复杂的非线性动态系统.相对于其他在时域或频域估计Volterra核的理论,该方法数学基础清晰、计算量不随辨识精度增高而增加、能够对任意高阶核进行估计,可对强非线性动态系统进行辨识,回避了计算量很大的计算流体力学(CFD).虽然对模型的喂给数据量有一定要求,但与随机和伪随机信号的随机响应法、Lee-Schetzen法和其他的统计量辨识方法比较,不需要对喂给数据量和概率分布有太高要求.通过该方法辨识出来的各阶核,可以立即得到系统的一阶幅相频响应和各阶广义频率响应.在NACA64-008A翼型加载线性扫频调幅拉伸力,将系统响应与辨识模型响应及Ijee-Schetzen法在不同喂给数据量下辨识结果做了对比和说明.     

透明塑料板件注射-压缩成型模拟研究

在非线性数值算法基础上,采用水平集方程、动量方程与能量方程进行多物理场全耦合计算的方法,对非对称几何板件的注射-压缩成型过程进行了模拟分析.结果表明,注射-压缩工艺较传统注射工艺能够显著降低制成件收缩率;成型过程中熔体的纤维流线长度与温度变化会对制件收缩率产生影响.较短的纤维流线,较均匀的温度分布能够获得更小、更均匀的...