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

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

广义协调条件的薄板弯曲三角形元和矩形元

其中兀_p和π_(mp)是最小势能原理和修正势能原理的泛函;w是单元内挠度;(?)、_(?)_n、(?)_s是单元间的位移;Q_n、M_n、M_(ns)是单元周线(?)A_e上的边界力。如果单元内部位移场在单元周线上的位移[w、((?)w)/((?)n)、((?)w)/((?)s)]和单元边界位移[(?)、(?)_n、(?)_s]~T满足广义协调条件     

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

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

基于非协调单元的液体火箭发动机推力室热结构分析

用巴兹公式和发动机的热力计算数据,得到了推力室内的燃气壁面对流传热系数和燃气温度的分布,并用有限单元法计算了推力室的温度场.在温度场计算基础上,采用Wilson非协调单元,对推力室的应力场进行了计算,其中位移函数采用二阶形式,将计算精度提高到二阶.结果表明:相比常规单元,Wilson非协调单元在计算结构比较复杂,温度梯度比较大的结构时,精度比较高,结果合理.      

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

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

旋转变厚度圆盘的样条积分方程法

一、前言 对于旋转变厚度圆盘的应力分析,一般可用有限元法或差分法,但相对于边界元来说精度较差,若要提高精度就必须将单元分细,随之而来的是方程阶数的增加和计算时间的增长。为了克服上述缺点,本文提出了一种新的样条积分方程法,直接利用弹性力学二维问题的环基本解建立积分方程并用样条函数进行求解。这样只需划分少量单元,就可获得精度很高的数值解,可大大节约计算机内存和计算时间。     

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

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

用广义乘子法求解航天器最优平面再入轨迹

从纵向平面无量纲航天器质点运动方程出发，引入了能量参数对运动方程进行推导，使得运动方程和优化问题易于处理。通过将优化变量转化为分段连续的线性函数和Runge—Kutta数值计算方法，将轨迹优化问题转化为非线性规划问题，应用广义乘子法对其进行数值分析，给出了不同航程要求下．总加热量最小，满足终端要求以及过载、动压、热流约束的航天器最优轨迹，并分析了其特点。     

孔边应力单元与层合板孔边应力精化分析方法研究

构造了一个满足孔边应力边界条件、层间应力连续条件以及层合板表面应力边界条件的孔边应力单元。并在位移型有限元计算所得位移基础上,提出了一个精确分析孔边应力的有限元方法。首先,运用莫雷阿(G.Morera)应力函数理沦,构造了柱坐标中满足平衡条件的三维应力场的一般表达式,并使之满足前述应力条件。然后,采用上述表达式并利用修正余能原理与拉格朗日乘子法,建立了由已知节点位移分析层合板孔边应力的基本方程。算例表明,所得计算结果与相应解析解或实验结果相符。     

旋转对称子结构罚单元分析

叙述了罚单元法和罚单元刚度矩阵的推导过程。列举了应用罚单元法处理旋转对称子结构边界位移协调和不同类型单元连接点位移协调的实例 ,验证了该方法的合理性和有效性     

极坐标系下谱元方法求解不可压缩Navier-Stokes方程

本文利用有限元的思想并结合谱方法的精度提出求解偏微分方程的谱元方法,并将谱元方法应用到极坐标系中;详细推导了在极坐标系下的谱元方法的具体计算公式,求解了极坐标系下的简单椭圆型二阶偏微分方程;并结合时间分裂方法应用于同心旋转圆筒间流体的流动问题,具体求解了原始变量速度和压力的不可压缩NavierStokes方程,均取得了满意的结果.     

一种同位网格上的迎风有限元方法

通过Ｇｒｅｅｎ定理将对流项变量从微分算子中分离出来,从插值函数入手引入迎风格式,是对强对流问题有限元计算中对流项变量的一种新的处理方法。按这种方法采用局部斜迎风格式及速－压同位网格公式,构成了一种对高Ｒｅｙｎｏｌｄｓ数流体流动数值模拟比较有效的有限元方法。     

Residual stress determination in friction stir butt welded joints using a digital image correlation-aided slitting technique

This paper presents an incremental cutting method for evaluating the longitudinal residual stresses in a butt welded thin plate via combining the traditional residual stress measurement methods and the advanced optical technique.The proposed approach,which can be called digital image correlation (DIC)-aided slitting technique,introduces a successive extension slot to a speci men and employs the DIC technique to measure the released displacement profiles of the cutting sections after each cutting increment.Then the displacement profiles are used to directly calculate the residual stress distributions up to the slot tip and hence,a stress distribution can be obtained after a cutting increment.Finally,all of the stress distributions are averaged to ultimately determine the original residual stress field.This method does not include any complex experimental operations or tedious derivation,and the resolution of stress variation is greatly improved by the continuous measurement of the released displacements.The presented method has been preliminarily verified by a specimen with residual stress introduced by a four-point bending test.The results show that residual stresses determined by the DIC-aided slitting technique agree well with those from finite element (FE) prediction.The residual stress in a friction stir welded aluminum specimen obtained by the presented technique is also consistent with the evaluations given by X-ray diffraction.Fur thermore,the residual stresses obtained by the DIC-aided slitting technique demonstrate higher accuracy and stability than the evaluations derived by the DIC-aided contour method.     

基于非局部理论的一维损伤模型的有限元算法

为了解决损伤局部化问题的网格依赖性,采用非局部理论,从损伤本构关系中引入特征长度入手,建立非局部损伤本构模型。提出了一种新的单元,将位移、非局部应变和非局部应变梯度作为基本的未知量,引入了C-1连续的形函数,并推导了非局部损伤本构模型的有限元格式。在ABAQUS中进行了二次开发,嵌入了本文提出的单元和本构模型,并用ABAQUS软件分析了一维损伤结构的受力情况。计算结果表明,数值模拟与理论推导结果的总体趋势吻合得较好,证明了本文方法的正确性与可行性。     

Hoff型夹层板振动分析的边界积分方程法

 本文利用静力问题的基本解建立了Hoff型夹层板振动分析的边界积分方程。在数值实现过程中,边界上和区域内的未知函数均釆用线性插值以提高计算精度。数值结果表明本文方法具有未知量少和计算精度高等优点。对于低阶频率,本文结果与精确解的误差不大于1％。     

计算机翼跨声速压力分布的格林函数方法

本文给出一种用格林定理将跨声速小扰动速势方程化成积分方程式,然后进行数值求解的方法,文中对场元强度进行松弛迭代,一般只需20次左右即达收敛。算例表明,该方法结果与差分法结果一致。     

旋翼桨叶动力分析及内力计算中有限元法的改进

提高用有限元法进行桨叶动力分析、特别是动内力计算的精度是本文工作的主要目的。本文建立了变剖面旋转梁协调单元族;提出利用单元动刚度矩阵直接由节点位移计算桨叶内力（简称为计算内力的“动刚度法”）。作为例子,用所导出的不同精度的有限单元和建议的内力计算方法计算了具有不连续结构特性的变剖面桨叶弯曲振动固有频率、振型和模态内力。计算结果表明,工作可以达到预期的目的。     

基于逆向有限元法的变形机翼鱼骨的变形重构

变形监测技术能够为自适应变形机翼的变形控制系统提供参考信息，是保证结构安全性以及优化结构的运行性能的重要手段。传统的基于光学成像的变形测量方法已经不能满足自适应智能结构的实时变形监测的要求。由于变形机翼表面受气动载荷影响，不便于直接在变形机翼蒙皮表面布置应变传感系统，目前还没有针对鱼骨结构这种真实复杂机翼结构的变形重构研究，大多针对机翼翼型的变形重构研究是将整个机翼简化成简单的翼形板、梁结构。针对上述问题，本文首次以真实复杂变形机翼主承力结构——鱼骨为研究对象，提出了一种基于逆向有限元（iFEM）算法与位移分段叠加思想结合的变形监测方法，根据Mindlin板变形理论建立四节点逆向壳单元，采用应变传感系统测得鱼骨结构表面应变分布作为算法输入，然后基于最小二乘变分方程求解结构应变场和位移场之间的传递函数，重构鱼骨结构的变形形状，为反演机翼翼型的变形形状提供方法。针对真实自适应变形机翼的主要承力构件开展了变形实验，实验结果表明，机翼鱼骨在分别偏转5°、10°、15°的情况下，逆向有限元法能准确重构鱼骨变形形状，验证了基于逆向有限元法的变形重构方法在真实自适应变形机翼结构变形重构研究中的有效性和准确性。     

大型工程结构的有限元动力模型优化法

研究以振动试验测得的飞机结构固有频率为目标值的有限元动力模型的优化问题。应用了一种单元组合法,该方法基于动力平衡方程以及正交性。通过修改质量矩阵的元素来优化,使计算出的固有频率与试验测得的频率值基本一致,并成功地应用于无人机水平尾翼模型上,效果良好     

夹层板精密成形的数值模拟

为了实现基于“点阵钉模、真空负压、蜂窝夹层”的双曲度面板精密成形过程的有限元数值模拟,分析了成形 2个阶段的变形过程,建立了简化的力学模形,以一种广义协调薄板弯曲四边形位移单元为基础构造了具有公共横向位移自由度的几何非线性四边形四节点分别考虑和不考虑横向剪切的多层板单元,通过牛顿 -拉斐逊法结合载荷、位移约束量的增量法求解几何非线性方程,较好地再现了真实成形过程