有限元计算模型的试验修正

本文阐述了有限元计算模型试验修正的一种方法,推导了求特征值和特征向量对结构参数导数的二阶矩阵摄动近似解,并在直升机桨叶的有限元计算模型的试验修正中得到成功应用。     

高速碰撞中Lagrange有限元方法及其应用

根据连续介质力学守恒方程和有限元离散方法,针对三维高速碰撞现象的特性,应用瞬时最小闰能原理,完整地推导了高速碰撞有限元基本列式和离散方程,并给出了高速碰撞问题的中心差分求解算法,讨论了非线性Ｌａｇｒａｎｇｅ动力有限元程序数值计算中应重视的问题,本文给出的典型算例说明本方法对高速碰撞问题分析的有效性。     

基于Euler方程解摄动的超声速旋转导数计算方法

基于笛卡尔网格和采用Euler方程基本流场解,通过Riemann不变量摄动,对旋转扰动引起的物面法向速度构建摄动关系,发展了一种超声速旋转导数快速计算方法。文中给出了具体的推导过程,可以看出由此得到的超声速旋转导数,仅与定常欧拉解基本流相关。因此,一旦CFD计算得到定常欧拉解,就可一次性地快速计算出感兴趣的超声速旋转导数。为了验证该计算方法,本文选用国际动导数标准模型Basic Finner算例进行了计算,计算结果与试验数据、文献数据具有很好的一致性,展示出本文方法计算耗时少、精度高、工程实用性强等优点,且该算法基于笛卡尔网格,特别适合复杂气动外形飞行器超声速旋转导数的快速计算。     

用配点法求结构瞬态响应

本文运用配点法,对用有限元建立的动力方程求结构的瞬态响应。文中提供两种计算方法;(一)逐步配点法;(二)最小二乘配点法。为了与各种逐步积分法作比较,文末有一算例,表明按这两种方法计算,其误差是较小的。     

基于热流体动力润滑分析的滑动轴承平衡位置和动力特性系数研究(英文)

建立了考虑润滑油粘热效应的滑动轴承热流体动力润滑分析模型,采用有限元法分析滑动轴承油膜压力和温度分布。对有限元方程关于扰动项求导直接得到油膜压力对于各扰动项的偏导数,并且通过油膜压力对于扰动项的偏导数构建油膜力Jacobi矩阵。采用Newton-Raphson方法确定给定静载下的轴承平衡位置,同时不需要额外的计算即可通过油膜力Jacobi矩阵获得滑动轴承动力特性系数。通过算例研究了轴承几何参数对滑动轴承动力系数的影响。分析结果表明,油膜温度对轴承刚度系数的影响大于阻尼系数;在轴承偏心率较小时,轴承承载力和刚度系数随转速的增加而增大;而在轴承偏心率较大时,轴承承载力和刚度系数随转速的增加而减小。     

求解织物热湿耦合方程的控制体-时域递归展开算法

为求解织物热湿传递耦合方程,给出了一种基于控制体积法和时域递归展开的求解织物热湿传递耦合方程的算法。首先,在时域对方程变量和参数进行级数展开,然后使用控制体积法对方程空间域进行离散,从而化连续性的非线性微分方程组为一系列的递归形式的线性代数方程组。文中给出了问题的求解步骤和算例。结论表明该方法的预测结果具有不依赖时间步长的特点。     

欧拉-伯努利梁平面超大挠性变形问题变分法求解

对于欧拉-伯努利悬臂梁平面超大挠性变形问题,由于其复杂的非线性几何方程,以位移为基本变量进行求解时,通常只能采用如多重打靶、微分求积等数值方法求得梁上离散点的位移值。本文研究了欧拉-伯努利悬臂梁平面超大挠性变形问题变分法求解理论。通过假设多项式形式的梁的曲率试函数以及常数中心线应变,基于欧拉-伯努利悬臂梁的基本假设,推导出了相互耦合的位移函数的精确表达式,并基于变分法理论和三角函数级数展开,推导出欧拉-伯努利梁的非线性控制方程组。利用迭代法对非线性控制方程组中的未知参数进行求解,最终得到欧拉-伯努利悬臂梁的位移函数的解析表达式。利用有限元计算结果对提出的变分法求解理论进行验证,并分别计算了欧拉-伯努利悬臂梁在自由端集中力及位移约束情况下的大变形。算例表明,基于本文的变分法求解理论,利用6个未知参数,即能够精确预测欧拉-伯努利悬臂梁在自由端集中力及位移约束下的超大挠性变形,该研究成果为欧拉-伯努利悬臂梁的超大变形问题提供了新的求解方法。     

用加权余量法建立三维涡流场的有限元方程

本文应用加权余量法推导了三维涡流场的有限元离散化方程,使得三维涡流场的有限元方程表达方便,概念清楚,计算较为简捷。     

纤维增强复合材料层合结构的非线性热稳定性分析

考虑到纤维增强复合材料的弹性常数和热应变随温度变化的非线性关系，推导出复合材料层合结构在热载荷作用下大挠度的有限元列式和几何非线性屈曲分析的理论表达式，利用局部线性化，得到了既考虑几何非线性又考虑热非线性的屈曲方程，本文提出了一种实现双重非线性问题的计算方法，使该问题的工程计算成为可能。     

求解织物热湿耦合方程的控制体一时域递归展开算法

为求解织物热湿传递耦合方程,给出了一种基于控制体积法和时域递归展开的求解织物热湿传递耦合方 程的算法.首先,在时域对方程变量和参数进行级数展开,然后使用控制体积法对方程空间域进行离散,从而化 连续性的非线性微分方程组为一系列的递归形式的线性代数方程组.文中给出了问题的求解步骤和算例.结论 表明该方法的预测结果具有不依赖时间步长的特点.     

具有随机刚度薄板小挠度问题的数值分析

从薄板小挠度问题的控制方程出发，将结构的随机参数统一到刚度随机场中，采用局部平均法对刚度随机场进行离散，并将其融于二阶差分方法的列式中，给出了求解刚度随机时薄板上各点位移的数字特征的随机有限元法。该方法不但可以解决复杂边界下的薄板小挠度问题，而且进一步解决了结构的可靠度问题。文中最后对算例进行了计算，用解析解验证了本文方法的正确性。     

轮辐框架式结构机匣的有限元强度计算方法

本文阐述在航空发动机上使用的轮辐框架式承力机匣,在受各种复杂载荷作用下的有限元强度计算方法。讨论了机匣力学模型的简化和建立有限元模型的三种方法,并对某发动机的扩散机匣做了强度计算,计算结果与试验结果基本相符。     

圆柱壳热应力的理论与实验研究

本文根据薄壳的经典理论,导出了在受有轴对称的多项式温度分布条件下,圆柱壳的位移、内力和热应力的理论公式。并以三次多项式温度分布为例,计算了三种边界条件(自由端、固支端、简支端)下,理论的和有限元的热应力曲线,对固支端进行了热应力实验,绘制了实验曲线。最后,对三种情况作了分析对比,结果表明,理论与有限元计算曲线相当吻合,计算与实验曲线分布趋势一致。     

基于随机有限元法空间梁板结构系统可靠性分析

飞行器工作环境复杂，影响其结构安全的不确定因素很多，并且结构失效会造成巨大的经济损失，甚至危及生命，所以对其进行结构可靠性分析是非常必要的。为了使可靠性指标计算简化，现有文献主要采用杆板结构来模拟飞行器的空间结构，并且主要基于有限元来求可靠性指标。本文采用空间梁元与加筋板格元（SP）来模拟其三维空间结构，模型更符合实际模型。同时考虑材料的强度、梁元件截面积、板元件厚度和外载荷均为随机变量，采用随机有限元法分析了随机变量对结构节点位移的影响，并应用可靠性基本理论对该结构系统进行结构可靠性分析。最后编制了相应的软件，计算了两个算例，结果证明该方法对飞行器结构可靠性分析具有重要的作用。     

基于SFEM的结构元件疲劳可靠性分析

将结构元件的弹性模量、泊松比、几何尺寸定义为随机场,并将随机场离散成随机变量,将抗拉强度定义为随机变量,采用Monte Carlo随机有限元法(MCSFEM),给出了元件危险点应力的随机分布。由危险点应力分布、材料的p-S-N曲线和Miner理论给出了概率名义应力法,计算常幅载荷下元件的寿命分布,并利用剩余寿命模型计算元件的可靠度。给出了两个算例,结果表明:计算结果与试验结果吻合良好。     

基于有限元法的齿轮强度接触分析

基于有限元法的基本算理,利用有限元软件,从静态和动态两个方面对齿轮进行接触分析,计算出齿根的最大应力,通过与经验公式计算结果对比,有如下结论：1）基于有限元法的齿轮接触分析是切实可行的计算方法,静态和动态模拟数据均比较合理;2）动态接触计算模型计及了啮合过程中的冲击效应和摩擦的影响,相对于静态算法而言,能更为真实的模拟齿轮啮合的真实状况;3）静态接触模型能准确快速地模拟稳态运行状态下的齿轮强度。齿轮的接触强度的合理分析,对改善齿轮传动性能和提高齿轮设计制造水平有一定的意义。     

接触问题数值分析方法的研究现状与发展

接触问题涉及到应力集中、边界非线性、甚至材料或几何非线性，问题很复杂但又十分重要。随着电子计算机的发展，有限元与边界元理论与数值方法也取得了很大进展。目前，这两种方法已成为工程接触问题十分有效的研究和数值分析手段。本文综述了接触问题中的有限元法与边界元法及其工程应用方面所取得的进展，着重围绕小变形弹性接触、小变形弹塑性接触以及大变形接触等方面的内容。此外，文中还对今后研究的发展方向作了展望。二十多年来，在小变形弹性接触问题中的有限元法和边界元法发展的比较成熟，并在工程中获得广泛应用；小变形弹塑性接触问题中的有限元法发展的也比较成熟。然而，它们在如下方面还必须开展进一步的研究，例如，不同材料非线性模式的接触、不同蠕变模式的接触、大变形弹塑性接触、基于实验力学的反迭代分析、以及复杂的工程构件的接触分析等。     

Stokes问题的集中质量非协调有限元法

针对其平面有界凸区域上的一类非定常不可压Stokes方程，提出了一种新的有限元方法，借助于所谓的速度一压力混合型公式，讨论了非定常不可压Stokes方程的质量集中非协调有限元逼近格式(全离散情形)。首先，给出了非定常不可压Stokes方程的质量集中非协调Galerkin有限元逼近的全离散格式，其次，对所讨论问题的解与其所给出的离散问题的解之间的误差进行了分析研究；最后，利用Stokes投影算子的性质和离散的LBB条件，得到了非定常不可压Stokes方程关于速度L2模和能量模及压力L2模方面的最优阶误差估计。     

A FINITE ELEMENT ADAPTIVE METHOD FOR SOLVING GENERALIZED STOKES PROSLEM

ＡＦＩＮＩＴＥＥＬＥＭＥＮＴＡＤＡＰＴＩＶＥＭＥＴＨＯＤＦＯＲＳＯＬＶＩＮＧＧＥＮＥＲＡＬＩＺＥＤＳＴＯＫＥＳＰＲＯＳＬＥＭＤ．Ａ．ＦｏｋｉｎＹａｎｇＺｕｏｓｈｅｎｇ（ＤｅｐａｒｔｍｅｎｔｏｆＡｅｒｏｄｙｎａｍｉｃｓ，ＮＵＡＡ２９ＹｕｄａｏＳｔｒｅｅ...     

有限元频响分析中网格划分技术研究

阐述了利用有限元进行频响分析过程中存在的一个突出问题；简单介绍了波动理论的一些基本公式，并提出有限元频响分析时划分网格的参考标准；首先通过对简单结构(包括梁、矩形板)的有限元计算以及与解析解的对比初步验证了该方法的有效性，然后给出了两个工程算例的应用情况。