不可压扰动方程高精度对称紧致差分数值解法及应用

从粘性不可压扰动方程一阶改型形式出发,对其实现了高精度对称紧致差分离散,就导出的扰动线性特征值问题给出了一个高铲双重迭代局部解法,以相同精度将特征值和特征函数同时得到。通过不可平面Ｐｏｉｓｅｕｉｌｌｅ流时间稳定性算例详细对比显示了算法良好的谱分辨能力和较弱的网格依赖性,并结合复矩阵广义特征隐式移ＱＺ算法获取了一个 扰动特征值谱计算结果。     

圆管内表面周向三维裂纹应力强度因子的新解法

对圆管内表面周向三维裂纹应力强度因子的求解,发展了传统的能量差率方法,同时研究了新兴的双重边界元理论,并运用这两种新解法,分别获得了关于裂纹张开位移的伯努利方程的半数值半工程封闭解和基于非连续裂尖奇异单元离散的应力强度因子数值解。并讨论了应力强度因子与裂纹形状a/c及裂纹深度a的关系,两者比较分析发现,其计算结果规律是一致的,且能量差率方法更保守,也更安全。     

某型航空发动机后支承动刚度的有限元计算

支承动刚度是影响发动机振动特性的主要因素之一,但目前航空发动机支承结构的刚度系数一般只能靠经验给出一个范围,这严重影响到转子动力特性计算结果的准确性。根据某型航空发动机的实际结构,建立发动机后支撑结构的动力学模型,应用有限元软件MSC.Nastran进行频率响应计算,再由所得位移曲线求出相应的动刚度曲线。本文所采用的计算方法和得到的结果对发动机转子动力特性的研究具有一定的参考价值。     

随机集的概率假设密度粒子滤波

多目标跟踪问题中,当目标数已知时,可以用概率数据互联(PDA)或联合概率数据互联(JPDA)算法。而当目标数未知或随时间变化时,需要对不同目标数的跟踪进行比较。可以把目标集看作随机集进行讨论,目标数N是随机变量。随机集的跟踪通过有限集统计(FISST)理论来完成。文中讨论了用粒子滤波实现跟踪随机集的方法。实验表明,在杂波环境下,粒子滤波可以稳健跟踪目标状态和目标数。     

基于SFEM的SRM纤维缠绕结构的可靠性分析

采用纽曼级数展开的蒙特卡罗随机有限元方法(SFEM),分析了固体火箭发动机(SRM)复合材料缠绕的燃烧室在内压正态变化下的应力响应,然后利用应力强度干涉原理分析计算其可靠性.     

激光冲击残余应力场的有限元分析

建立模拟激光冲击残余应力场的非线性弹塑性有限元模型,预测AISI 304奥氏体不锈钢厚、薄靶材经激光冲击后表面以及内部的残余应力场分布.模拟过程中,考虑材料高应变率(106s-1)下的力学性能和激光诱导冲击波的精确加载,分析激光功率密度、激光光斑尺寸、激光脉宽、激光冲击次数、激光单面与双面冲击等激光冲击参数以及初始残余应力等因素对不锈钢靶材残余应力分布的影响.数值模拟结果与X射线衍射法测量值吻合较好.在有初始残余拉应力(250MPa)存在的情况下,激光冲击仍能使304奥氏体不锈钢靶材形成残余压应力层,说明激光冲击工艺可提高奥氏体不锈钢焊接构件的抗应力腐蚀开裂能力.     

直升机桨叶刚柔耦合特性及计算方法分析

为了研究桨叶刚柔耦合特性,采用中等变形梁模型处理桨叶弹性变形,分别采用小转角和有限转角两种方法处理桨叶绕铰的刚性运动。以桨叶扬起下坠碰撞问题为研究对象。计算分析表明:(1)弹性桨叶存在较明显的刚性运动和弹性变形耦合;(2)挥舞角较小时,两种处理方法计算结果一致,均与试验数据吻合较好;(3)挥舞角较大时,刚柔耦合增强,两种处理方法的计算值差别较为明显;(4)桨叶弯曲刚度对桨叶刚柔耦合特性有较明显影响,刚度越小,耦合越强;(5)桨叶如果比较柔软,小转角处理方法计算值相对有限转角处理方法,相位滞后、幅值增加均较为明显,计算桨叶动响应时需引起注意。      

基于Euler／N—S方程的跨音速非线性静气动弹性问题研究

在C-H网格的基础上,采用Jameson的中心差分有限体积法求解Euler/N-S方程,采用结构影响系数法计算结构的弹性变形,用三角元面积加权法和常体积转换法(CVT)实现流固耦合,计算了弹性后掠机翼在跨音速状态的静气动弹性问题     

平移速度的变化对搅拌焊接材料流动的影响

基于有限元方法对不同搅拌头平移速度下搅拌摩擦焊接过程中焊接构件Al6061-T6材料的三维流动问题进行分析,研究了搅拌头平移速度的变化对等效塑性应变的分布和材料流动的影响。结果表明,搅拌头前方的材料向上涌起,被旋推到搅拌头后方并开始向下运动,这一过程是促使搅拌摩擦焊接顺利完成的主要原因。材料的等效塑性应变等值线与不同焊接区域的边界具有较好的对应关系。随着搅拌头平移速度的增加,垂直于焊缝方向上的截面上,等效塑性应变减小,但是材料的流动速度会随之增加。较低的搅拌头平移速度有利于焊缝中心线两侧材料的融合。     

求解Euler方程的无网格与有限体积混合算法研究

研究了无网格方法与有限体积法相结合的、求解Euler方程的混合算法.与单一的无网格方法相比较,由于在大部分计算区域采用了有限体积法,使得发展的算法在运算效率上能与有限体积法相当;同时在物体附近嵌入了无网格区,使得在处理几何外形上更具灵活性.有关算法涉及的无网格与网格区域搭接处理,通过在交界面处引入辅助卫星点和网格单元构成边界信息,实现了区域间的流动信息传递.Euler方程的空间导数分别在两区域用有限体积法和无网格法离散近似,时间方向都采用四步显式Runge-Kutta格式推进求解,数值模拟了喷管内流和绕翼型外流,并分别与整体有限体积法和整体无网格方法进行了比较.算例展示出,提出的混合算法能有效捕捉激波间断,且两区域等值线过渡光滑,算法效率如预期与有限体积法相当,表明该方法是可行的.