组合杂交三角形单元的加权能量正交关系

组合变分原理可以增强杂交元方法解的稳定性。建立热传导方程基于区域分解的组合杂交有限元方法,给出单元上温度梯度插值为线性、但温度插值为协调线性插值与非协调二次插值之和的组合杂交三角形单元,并通过数值实验验证理论结果的正确性。结果表明:分片线性温度梯度插值的散度(热源)与非协调温度插值是加权能量正交的;组合杂交三角形元刚度矩阵等同于协调的三角形线性元刚度矩阵,即非协调部分无温度增强特性。     

瞬态热应力问题的组合杂交有限元方法

求解瞬态热应力问题的杂交元方法面临Ladyshenskaya-Babuska-Brezzi(LBB)条件的检查,其稳定性难以保证,从而提出瞬态热应力问题的组合杂交有限元方法。建立瞬态热应力问题的基于区域分解的组合杂交变分原理及其有限元离散;通过数值算例验证求解瞬态热应力问题的组合杂交元的数值性能。结果表明:在大宽厚比网格下,八节点六面体组合杂交元(CHH(0-1))可以获得与二十节点六面体元(B20)精度相当的位移计算结果和更好的应力计算结果。     

基于非协调模式的几何非线性广义杂交退化壳单元

依据一个在局部坐标系下放松单元间连续性要求的广义变分原理,建立了基于非协调模式的几何非线性广义杂交退化壳单元。数值结果表明,构造的四结点退化壳元拥有很好的性能,可避免各种自锁问题。     

复合材料层合板的杂交有限元方法

结合复合材料修正后的H-R混合变分原理,直接借助应力-应变关系,推导了新的应力模式,建立了复合材料层合板的杂交等参有限元列式.利用Mathematica语言编程进行数值实例分析,其计算结果与相关文献的精确解以及Abaqus软件建模分析结果对比,实例证明该方法所得到的各个静力学量更接近精确解,并且可用较少的网格划分得到较精确的解.     

一种新的含转角自由度的四节点杂交元

基于组合杂交变分原理,使用含转角自由度的位移插值和CH(0-1)元应力模式,得到CHD(0-1)单元.它既满足杂交方法所要求的inf-sup条件,又有较高的位移、应力精度;和普通的八节点二次单元相比计算花费少;对单元几何形状不敏感;能很好地模拟弯曲变形.数值试验结果支持了本文结论.     

关于计入剪切效应的板弯曲杂交模式的若干问题

本文从Reissner变分原理出发,构造一计入剪切效应的9自由度板弯曲杂交元。通过离散型Kirchhoff约束的实施,把厚板公式直接推广于薄板计算而不出现解的“闭锁”现象。文中提出的能量调整原理有效地解决了厚板元素中的弯、剪匹配问题,从而实现了厚薄板分析的统一性。文章最后给出计算实例。     

运动激波自适应网格算法中权函数问题的研究

变分原理的自适应网格技术被应用到运动激波问题的求解上,在解的大梯度区自动加密网格,从而非常成功地算出了激波。通过分析发现权函数选择和取值对自适应网格技术至关重要     

基于重叠网格技术和多重网格算法的悬停旋翼粘性绕流数值模拟

发展了一种基于重叠网格的多重网格算法,并应用于悬停旋翼粘性绕流的Navier-Stokes方程数值模拟.多重网格算法在重叠网格上实施主要面临如下两个问题:粗网格间几何关系的建立,即洞单元、边界单元的确定;洞边界存在的情况下如何进行粗细网格间流场信息的传递.本文针对以上问题给出了相应的解决方案.采用两层网格的V循环,应用有限体积法对一跨声速悬停旋翼粘性流场进行了数值计算.计算结果表明该多重网格方法可以显著加快悬停旋翼粘性流场数值解收敛速度,加速比约为3.     

应用复合型网格求解叶栅气动设计杂交型问题

建立在流线坐标系上的叶栅气动设计反命题与杂交命题是在映象平面(由流线坐标与周向坐标构成H上求解的.以往在映像平面只采用H型网格,本文首次将H加C型复合网格用于这类反命题计算,使叶栅前缘附近网格构造趋于合理,从而明显提高了计算精度.本文采用的解法是基于叶轮机气动变分原理的有限元方法.     

径流式叶轮机S_2流面半反命题、A型杂交命题变分原理族及对混流式的推广及其应用

(一)文[4、5]建立了轴流式叶轮机S_2流面半反命题、A型杂交命题的变分原理,为应用有限元方法和其它变分直接解法提供了较完备的理论基础。但是由于相应尤拉方程组解的不定性,这两族变分原理并不适用于径流式叶轮机。     

一种基于有限元法和弹性接触理论的齿轮啮合刚度改进算法

提出了一种将有限元法和弹性接触理论相结合的齿轮啮合刚度计算方法.该方法利用子结构法提取齿面原始柔度矩阵并分离出接触点弯曲变形,根据线弹接触变形解析公式计算接触变形,通过求解非线性变形协调方程得到齿轮时变啮合刚度和齿面载荷分布.以一对齿轮副为例,计算的啮合刚度与航空标准计算结果相差在6%以内.该方法发挥了有限元法在预测物体整体变形方面的优势,同时结合弹性接触理论能够准确计算局部接触变形的优点,与常规有限元法相比,能够有效地提升计算效率.由于接触变形问题的非线性,啮合刚度随总啮合力增加呈现非线性增大的趋势.      

结构分析程序中几个问题的处理方法

JY851程序是一个复合材料、金属混合结构应力分析及强度校核程序。从1986年投入使用以来,已推广到几个厂、所使用,计算了几个型号的机翼、垂尾、方向舵等复合材料结构。现就编写元素刚度矩阵中所遇到的几个问题的处理方法作一简单的介绍。 1.复合材料板弹性模量阵的计算 复合材料板弹性模量阵A、B、D的计算中,关键是求偏轴模量阵[ō],层压板中第     

机翼和翼身融合体跨音速绕流有限元算法的研究

从全速位方程出发研究了机翼和翼身融合体跨音速绕流的有限元算法。给出适用于有限元法和复杂外形的块结构网格生成技术和使单元系数矩阵积分数目减少2/3的方法;并发展了一种内存量少的总体系数矩阵存贮方式。为了提高有限元离散方程组解的精度和加速收敛,提出α—GMRES算法代替传统全速位方程有限元解法中所习用收敛慢的线松弛迭代解法。通过机翼和翼身融合体实例考核,效果是好的,为飞行器先进气动布局的数值计算提供了一种新的有效算法。     

圆柱筒轴向载荷下单元尺寸效应数值模拟

利用显式非线性有限元方法对薄壁圆柱筒轴向载荷下的动力学响应进行了研究,模拟了圆柱筒在轴向冲击作用下的动力屈曲变形及其发展过程。分析比较了不同单元尺度下冲击模型的数值计算解与理论解,揭示了在进行有限元模拟时存在的网格尺寸效应,为研究有限元网格尺寸与有限元求解精度的内在联系提供参考;为在保证数值解满足工程实际精度要求的前提下,确定合理的单元尺度,提高有限元分析效率进行了有益的探索。     

一种确定模型误差位置的方法

本文首先由解析有限元模型的模态矩阵,得到未测量自由度与已测量自由度模态之间的变换关系,利用该关系对解析模型进行缩聚。然后针对缩聚后的特征方程,利用解析刚度矩阵和质量矩阵以及实测的模态参数确定模型误差位置。模拟结果表明,该方法可行。     

二维振荡机翼跨声速非定常绕流的变域变分有限元解

本文首先对文献（１－３）提出的二维振荡机翼含激波跨声速非定常绕流的变分原理进行了改造，使之能适用于时间推进法；构造了三维时空的可变节点有限元来捕获自由尾涡面，而跨声流中的激波采用人工密度法捕获。在远场边界上使用简化的无反射条件和新型非定常Ｋｕｔｔａ条件。用本文方法求解了作俯爷振荡的ＮＡＣＡ６４Ａ０１０翼型的跨声速绕流，计算结果令人满意。本文方法可以推广到三维机翼及二维和三维叶栅的同类气动力问题上去     

高维局部非线性转子-轴承动力系统的稳定性和分岔

采用模态缩减方法对具有非解析轴承的高维局部非线性转子系统进行自由度降阶。提出基于变分不等方程的有限元互补解的的数值方法,在几乎不增加计算量的情况下,使得实际轴承油膜力的Jacobian矩阵可与油膜力自身计算同时完成,并取得协调一致的精度。结合打靶法和Floquet理论对实际转子-轴承系统的非线性不平衡响应及其分岔行为进行计算分析,数值结果表明,所提出的方法不仅极大地降低了计算量,而且具有足够高的精度。     

有限元的一种并行算法及实现

当前,由于有限元素法大都使用串行算法,处理时间较长。本文根据有限元素法的特点和采用的并行计算机结构,提出了一种并行有限元算法,能大大提高其处理速度,加快处理周期。且算法简单,实现方便。     

非对称铺层复合材料板室温形状的非线性有限元分析

由于一般铺层复合材料板壳存在拉伸-弯曲、拉伸-扭转等耦合效应,在温度变化时将出现垂直于中面方向的翘曲变形,给工程使用和理论分析带来了极大的不便。但另一方面,利用非对称铺层复合材料板壳可以有效地减轻结构的重量,并可能提供利用平面模具制造曲面层合壳的新工艺。