张拉结构非线性分析两节点曲线单元有限元法

提出适合张拉结构几何非线性分析的两节点曲线单元有限元方法．假定索元的初始形状呈二次抛物线,根据单索的平衡条件、几何和物理关系建立了索元的位移函数;由拉格朗日应变的定义建立了可以考虑任意次高阶位移影响的索元轴向应变的精确表达式,并基于拉格明日描述方法和虚功原理得到了索元的非线性平衡方程与切线刚度矩阵．采用荷载增量法与Newton－Raphson法相结合的混合法进行了实例计算,结果表明：本文方法的精度明显优于两节点直线索单元,适合于大跨度索阿、索穹顶等张拉结构的几何非线性分析．     

A MESHLESS LOCAL PETROV-GALERKIN METHOD FOR GEOMETRICALLY NONLINEAR PROBLEMS

Nonlinear formulations of the meshless local Petrov-Galerkin （MLPG） method are presented for geometrically nonlinear problems. The method requires no mesh in computation and therefore avoids mesh distortion difficulties in the large deformation analysis. The essential boundary conditions in the present formulation axe imposed by a penalty method. An incremental and iterative solution procedure is used to solve geometrically nonlinear problems. Several examples are presented to demonstrate the effectiveness of the method in geometrically nonlinear problems analysis. Numerical results show that the MLPG method is an effective one and that the values of the unknown variable are quite accurate.     

一种基于基线力的平面几何非线性余能原理有限元模型

以基线力为状态变量,构造几何非线性问题单元的余能.该余能包括变形部分和转动部分.利用高玉臣提出的弹性大变形余能原理,以Lagrange乘子法松弛单元域内的平衡条件,得到了一个以基线力表达的修正的余能原理.在假设平面单元每一个边上的力为均匀分布的条件下,推导出一种几何非线性平面4节点有限元模型.运用MATLAB语言编制出相应的非线性有限元分析程序.数值算例结果表明该模型具有良好的性能.     

基面力概念在几何非线性余能有限元中的应用

以基面力为基本未知量描述一个弹性系统的应力状态并表征单元的余能,将大变形的余能分解为变形余能部分和转动余能部分,采用Lagrange乘子法放松单元的平衡方程,利用已有的弹性大变形余能原理建立了一种几何非线性显式有限元模型,编制了相应的几何非线性余能原理有限元程序. 数值算例表明:该方法具有较好的收敛性和计算精度,可进行大载荷步的大位移、大转动计算.      

非线性有限元的若干基本问题

本文介绍了非线性有限元中的若干基本问题。其中包括有关应变、应力和非线性平衡方程的一些基本概念,基于不同非线性广义变分原理的位移模式、杂交模式和拟协调模式几何非线性有限元及其在壳体屈曲问题中的应用等。      

数学规划加权残值法分析薄板与薄壳的几何非线性问题

使用数学规划加权残值涯分析了薄板大挠度及薄壳的非线性稳定问题。在两者应变协调方程严格满足的情况下,论证并利用平衡方程的单调性,建立数学规划问题,首先得到无布载荷四边简支方板中心挠度的最小上界及最大下界,经典Ｌｅｖｙ解位基其间。     

非线性动力学方程的精细积分算法

介绍用精细积分法求解动力学问题的原则和方法,通过实例证明用这种方法求非线性问题数值解的有效性．     

圆柱壳的非线性有限元分析及其并行求解

采取16节点曲线边等参元对圆柱壳的几何非线性进行有限元分析. 分析考虑了完全非线性运动关系以便预测在非线性区域的稳定平衡路径. 建立了基于广义非线性位移的有限元公式. 提出了基于全Lagrangian格式的非线性有限元分析的并行计算策略. 在集群环境下, 对圆柱壳的几何非线性分析进行了并行计算. 计算结果表明: 在集群环境下,所提并行算法具有良好的加速比和效率.     

基于欧拉描述的两节点索单元非线性有限元法

本文针对柔性悬索结构几何非线性分析的特点,提出了一种用欧拉描述来表示的两节点索单元非线性有限元模型,在索元变形后的位置上由虚功能建立了非线有限元基本方程及切线刚度矩阵。这样建立的非线性有限元分析方法可充分考虑拉索的几何非线性特性的影响并给悬索结构的初始平衡分析带来方便,算例结果表明,本文方法是精确有效的。     

非线性广义杂交元

选择新的三类变量建立了放松单元间连续性条件的用于非线性有限元分析的泛函，并由此建立了基于不协调模式的非线性广义杂交元方法。     

用复变量表示的环壳几何非线性方程及其摄动解

本文从小应变,中小转动问题的弹性薄壳一般方程出发,导出用复未知函数表示的非线性的环壳轴对称变形的基本方程,并用摄动方法分別求解了整环壳和开口环壳问题.所得到的结果与其他作者用数值方法得到的结果吻合很好。     

薄板弯曲问题的P型杂交解析有限元方法

本文采用两套变量构造有限元试函数空间，在单元内部要求试函数精确满足平衡微分方程，在单元边界上对位移和转角分别用Ｐｅａｎｏ升阶函数插值，然后利用广义变分原理建立了一种薄板弯曲问题的Ｐ型杂交解析有限方法，与常规有限元法相比，该方法不心进行过细的网格剖分，通过增加单元插值多项式的阶数Ｐ来提高精度，此外，该方法还具有积分计算只需在单元边界上进行、单元钢度矩阵和载荷向量具有嵌入结构、协调程度可以自动控制等优     

GEOMETRICALLY NONLINEAR FINITE ELEMENT MODEL OF SPATIAL THIN-WALLED BEAMS WITH GENERAL OPEN CROSS SECTION

Based on the theory of Timoshenko and thin-walled beams, a new finite element model of spatial thin-walled beams with general open cross sections is presented in the paper, in which several factors are included such as lateral shear deformation, warp generated by nonuni- form torsion and second-order shear stress, coupling of flexure and torsion, and large displacement with small strain. With an additional internal node in the element, the element stiffness matrix is deduced by incremental virtual work in updated Lagrangian （UL） formulation. Numerical examples demonstrate that the presented model well describes the geometrically nonlinear property of spatial thin-walled beams.     

非线性复合材料杂交应力有限元的有效迭代方法

推导了面内剪应力应变关系非线性的复合材料的杂交应力有限元列式,给出了位移迭代和应力迭代的策略和步骤.提出一种非线性应力场迭代格式的改进方案,不仅提高了收敛速度,而且克服了大载荷下简单迭代法循环迭代而无法收敛的关键问题,使得所提出的非线性杂交应力元方法几乎对任意大载荷都能够收敛.数值算例表明该方法是确实可行的.     

不同本构模型对橡胶制品有限元法适应性研究

为确定不同本构模型对橡胶制品进行非线性有限元分析的适应性,采用HyperMesh 和ABAQUS软件建立了橡胶标准试件和防尘罩制品的非线性有限元模型. 根据标准试件单轴拉伸试验数据,利用ABAQUS软件拟合了Mooney——Rivlin, Ogden 及Yeoh 3 种本构模型的特征参数,并给出了试件和制品轴向拉伸和压缩工况有限元分析结果. 实测试验与有限元分析结果对比发现:在应变不大于100% 时,两参数的Mooney——Rivlin本构模型计算误差较小;在应变大于100% 时, 用Yeoh 和Ogden 本构模型计算误差较小. 该结论对大变形橡胶元件的非线性特性理论分析具有实际意义.     

黏弹性功能梯度材料裂纹问题的有限元方法

针对组分材料体积含量任意分布的黏弹性功能梯度材料裂纹问题建立有限元分析途径. 通过Laplace变换,将黏弹性问题转化到象空间中求解,基于反映材料非均匀的梯度单元和裂纹尖端奇异特性的奇异单元计算象空间中的位移、应力和应变场,应用虚拟裂纹闭合方法得到应变能释放率,分别由应力和应变能释放率确定应力强度因子. 给出这些断裂参量在物理空间和象空间之间的对应关系,由数值逆变换求出其在物理空间的相应值. 文中分析两端均匀受拉的黏弹性边裂纹板条,首先针对松弛模量表示为空间函数和时间函数乘积的特殊梯度材料进行计算,结合对应原理验证方法的有效性. 然后分析组分材料体积含量具有任意梯度分布的情形,由Mori-Tanaka方法预测象空间中的等效松弛模量. 计算结果表明,蠕变加载条件下,应变能释放率随时间增加,其增大程度与黏弹性组分材料体积含量相关. 由于梯度材料的非均匀黏弹性性质,产生应力重新分布,导致应力强度因子随时间变化,其变化范围与组分材料的体积含量分布方式有关.     

求解固体力学问题的强?弱耦合形式单元微分法

单元微分法是一种新型强形式有限单元法. 与弱形式算法相比, 该算法直接对控制方程进行离散, 不需要用到数值积分. 因此该算法有较简单的形式, 并且其在计算系数矩阵时具有极高的效率. 但作为一种强形式算法, 单元微分法往往需要较多网格或者更高阶单元才能达到满意的计算精度. 与此同时, 对于一些包含奇异点的模型, 如在多材料界面、间断边界条件、裂纹尖端等处, 传统单元微分法往往得不到较精确的计算结果. 为了克服这些缺点, 本文提出了将伽辽金有限元法与单元微分法相结合的强?弱耦合算法, 即整体模型采用单元微分法的同时, 在奇异点附近或某些关键部件采用有限元法. 该策略在保留单元微分法高效率与简洁形式等优点的同时, 确保了求解奇异问题的精度. 在处理大规模问题时, 针对关键部件采用有限元法, 其他部件采用单元微分法, 可以在得到较精确结果的同时, 极大提高整体计算效率. 在本文中, 给出了两个典型算例, 一个是具有切口的二维问题, 一个是复杂的三维发动机问题. 针对这两个问题, 分析了该耦合算法在求二维奇异问题和三维大规模问题时的精度与效率.      

用于几何非线性分析的内参型非协调元法

本文把内参型非协调元的理论和方法推广应用于几何非线性分析。对方法的合理性、可行性、单元收敛条件及实用简化方案等诸方面作了初步探讨。文章以大位移平面问题为例进行了详细讨论,给出了有限元非线性列式的全过程和数值结果。