三维瞬态导热温度场边界元计算及一阶数学奇异的处理

推导了三维瞬态导热温度场边界元解法的积分方程和离散代数方程，并编制了相应的计算机程序，采用坐标变换的解析消去法，处理边界元的一阶数学奇异点，从而提高了算法的数值精度。用本程序计算了正方体和圆柱环瞬态导热温度场，数值解与解析解的吻合程度令人满意。     

基于特征的锻件分类及特征有限元法

综合分析锻件后引入特征技术并归纳锻件的形状及变形特点 ,提出了锻件特征概念并给出锻件形状特征的一般定义和构成 ;采用面向对象技术表示特征 ,提出了锻件特征类的概念模式和具体的参数化数据模型 ;针对有限元分析中复杂多变的约束条件和载荷状况两个关键性因素的表达 ,提出特征类的关联思想。在基于特征的锻件分类基础上 ,把锻件变形归纳为“挤”或“压”两种特征元 ,提出了采用基于特征的有限元法 ,并综合自适应FEM、并行FEM、智能型FEM共同求解大型复杂锻件的变形模拟问题。     

直升机坠毁响应的有限元模拟

 本文概述用非线性有限元技术对直升机结构的坠毁过程进行数值模拟。所分析的结构包括飞行员座椅、座椅减震器、地板盒段和滑撬。问题被当作结构在规定初速度或脉冲载荷作用下的非线性（大变形与弹塑性）动反应。用特殊设计的单元模拟减震器及座椅滑轨,用间隙单元模拟结构与地面的撞击。     

翼型风洞侧壁干扰的数值模拟研究

运用Navier-Stokes数值模拟对翼型模拟试验对风洞侧壁干扰进行模拟，将简单代数湍流模型扩展用于机翼／风洞侧壁结合区流动，分析了风洞实验侧壁干扰问题的形成机理和影响翼型实验侧壁干扰的各种因素，如翼型展长、风洞侧壁边界层厚度及侧壁边界层抽吸等，对实验结果的影响，得出了一些有用的结论。计算格式空间采用中心有限体积离散，时间采用多步Runge-Kutta时间步长格式进行积分。计算结果证明了该方法的可行性和优越性。     

二维对流—扩散方程的控制体有限元法数值模拟

推导了二维对流—扩散方程的控制体有限元公式，其主要特点是所得离散化方程适合于逐点迭代法及共轭梯度法进行求解，使得该方法更适合于求解复杂的几何域及节省更多的计算机内存。通过几个算例表明，该方法对于特别大的Ｐｅｃｌｅｔ数其精度也是很高的     

复合材料叠层结构剖面翘曲修正理论

 首先假设复合材料叠层结构的面内位移沿厚度方向为直线分布,然后根据剪切变形进行连续修正,获得逐步升高的各阶翘曲位移分布函数。文中以叠层梁为例,应用最小势能原理导出其静力方程和力的边界条件。在长度方向则采用升阶谱位移函数,进行了数值计算。所得叠层梁的位移和应力分布的数值结果与精确解相比,表明本文方法有很好的收敛性。     

缘板及叶冠间的接触状况及其对叶片位移和应力的影响

1.计算方法及其实验验证 本文用分析法见文献[1],这里以WS9高压二级涡轮转子叶片(简称WS9叶片)为例来说明方法的特点。该转子轮盘上均匀安装110片叶片,组成具有循环对称接触结构的叶片环,有接触表面的循环对称结构。据此,可将图1所示的一片叶片取作基本段分析,但同时必须在该叶片一般有限元刚度方程中引入缘板及叶冠的循环对称接触表面(叶冠的A1侧面与A2侧面,以及缘板的B1侧面与B2侧面)处的循环对称接触     

复合材料旋翼桨叶剖面扭转刚度计算

 1、各向异性弹性柱体扭转的理论 复合材料旋翼桨叶的扭转问题简化为各向异性的弹性柱体的扭转问题,忽略桨叶的根部影响,是一个弹性柱体的自由扭转问题。 假设各向异性的弹性柱体,其正截面的形状不受限制（图1）可以是单连通的亦可以是复连通的（图2）。柱体内部不受外力（包括重力）作用。柱体周围侧面也是悬空的,即     

基于变域变分有限元的翼型反设计

反问题作为一种自由(未知)边界问题近年来得到了广泛的研究。变域变分通过把自由边界结合在变分泛函中，使其与求解流场的控制方程结合起来，从而自由边界求解和流场分析可以完全耦合进行。本文采用变域变分和可变形有限元来设计可压缩和不可压缩翼型，计算显示即使在初始翼型的一侧为三角形时仍能得到满意的结果。文中还采用了两种计算网格(H型和O型)来设计翼型，比较了它们的精度，计算表明本方法完全能用于翼型的反设计。     

关于厚薄板弯曲分析的内位移杂交元法

 本文通过引入单元内位移和优化处理的单元应力场,以计入剪切变形的板弯曲能量泛函为基础,导出适用于厚、薄板分析的四边形杂交应力元。内力场和内位移的合理选择避免了单元零能模式的出现、避免了解的剪切自锁现象。该优化单元显示了较好的通用性和收敛性。对关于剪切自锁的RCI判据提出了反例、说明了它的局限性。