基于ANSYS二次开发的几何非线性平面梁单元

共旋坐标法(C.R法)具有在局部坐标系考虑材料非线性,通过局部坐标系与结构坐标系之间内力和切线刚度矩阵的转换矩阵来考虑几何非线性,从而实现两种非线性脱耦的优点,C.R法相对于其它非线性有限元列式而言较少运用于商业程序。本文利用ANSYS平台提供的单元二次开发工具——用户可编程特性(UPFs),开发了基于C.R法的几何非线性平面梁单元,给出了详细的算法及流程,通过多个算例对本文算法及程序进行了验证。研究表明:该方法能有效利用共旋法非线性单元和ANSYS商业程序的优点,对于方框架在对边中点受一对集中力的算例1,采用二次开发的用户梁单元与beam3梁单元在每个荷载步下收敛所需的迭代次数分别为3次和6次;对于预应力钢筋混凝土悬臂梁发生大变形的算例2,上述两种单元模型进行非线性计算能收敛的预应力加载系数分别为132和128,可知本文基于共旋法得到的ANSYS二次开发用户梁单元提高了非线性计算的效率和能力,可用于平面梁结构的几何非线性分析。     

模拟结构倒塌破坏的有限质点法

有限质点法是以向量式力学为基础的新兴结构分析方法,本文将其应用于冲击荷载作用下的网壳结构的倒塌破坏模拟中。以空间杆单元为例建立了有限质点法的基本方程,推导了求解几何和材料非线性问题的基本公式。为计算断裂问题,建立了空间杆单元的断裂准则和断裂模型,发展了有限质点法进行断裂分析的基本算法。通过对某双层网壳冲击荷载下破坏过程的模拟和分析,验证了该方法在结构倒塌破坏过程模拟中的有效性和适用性。     

空间杆系结构实用几何非线性分析

从简单实用的角度论述了空间杆系结构的几何非线性分析方法。文中分析了非线性有限元方法的求解过程,特别强调决定几何非线性收敛结果的关键问题,即由节点位移增量计算单元的内力增量。通过引入转随转坐标系,论述了平面和空间梁单元小应变时单元内力增量的计算问题。针对杆系结构的大应变问题,从有限应变理论出发进行分析,提出了对该问题的有效处理方法,并且用实例进行了验证。计算结果表明,该实用几何非线性分析方法是可靠和有效的。     

柔性多体动力学的共旋坐标法

共旋坐标法作为柔性多体动力学的建模方法已经有数十年的发展历程.作为解决大转动小应变问题的可靠方法,共旋坐标法建模简单,意义明确,并且能借用现有的有限元理论,求解时效率高,精度好.因此柔性多体动力学中,共旋坐标法不仅在非线性分析中有良好的表现,在工程实际中也有较好的应用潜力.柔性多体系统动力学中的大变形问题大多涉及梁、板壳等柔性结构,有大量的大转动小应变问题需要解决.本文对采用共旋坐标法来分析梁、板壳类结构的研究中存在的若干关键问题进行了相关的讨论,如有限转动理论、局部坐标系的定义、单元类型和相应的动力学问题等.比较了各种方法的优势与局限,并且展望了共旋坐标法的发展趋势.     

一种复合材料旋翼桨叶结构大变形分析方法的改进及验证

提出了一种适用于直升机旋翼复合材料桨叶大变形分析的改进方法。将旋翼桨叶变形分析分解为一维非线性分析和二维剖面特性分析,并考虑横向剪切、翘曲对剖面刚度及弹性耦合的影响;为使方法适用于旋翼气动弹性分析,将应变能中的广义应变用参考轴线处的弹性运动表示,保留所有非线性项,推导出计算复合材料桨叶大变形的公式;采用有限元法处理方程,对梁结构进行了分析,并将大变形状态下的位移计算结果与Princeton梁实验值、Minguet复合材料梁实验值以及中等变形梁理论计算结果进行了比较,验证了大变形状态下本文计算方法的正确性;此外与中等变形梁模型计算结果的对比,验证了本文方法在计算精度上的提高。     

超静定梁弹塑性加载过程分析的教学内容构建

"超静定梁的塑性极限分析" 作为塑性力学教材中的一节内容,阐述了如何用"机动法" 和"静力法" 求最终的塑性极限破坏载荷,却没有分析超静定梁的弹塑性加载变形过程. 通过把结构力学中计算弹性位移的单位载荷法扩展应用到超静定梁的弹塑性加载过程,以均布载荷作用下两端固支超静定梁的弹塑性加载和变形全过程分析为例,构建了超静定梁弹塑性加载过程分析的教学内容,给出了两端固支超静定梁在均布载荷加载过程中弯矩内力和挠度随外载荷而变化的解析公式. 主要目的是引导学生掌握超静定梁复杂的非线性弹塑性加载变形全过程的分析方法,可供塑性力学教材改编时参考引用.     

基于六面体网格的向量式有限元分析及应用

基于向量式有限元基本原理,给出了八节点六面体等参实体单元的基本公式,通过投影方式将空间曲面六面体转换为投影六面体,采用参考面的逆向运动求解节点纯变形,通过单元形函的虚功方程计算节点内力;针对坐标模式和内力积分模式等关键问题提出了有效的处理方案。编制了六面体实体单元的数值计算程序,并进行工程结构算例分析。结果表明,所编制程序可有效模拟实体结构的静力、动力及大变形大位移行为分析,验证了本文理论和程序的正确性和实用性。     

基于六面体网格的向量式有限元分析及应用

基于向量式有限元基本原理,给出了八节点六面体等参实体单元的基本公式,通过投影方式将空间曲面六面体转换为投影六面体,采用参考面的逆向运动求解节点纯变形,通过单元形函的虚功方程计算节点内力;针对坐标模式和内力积分模式等关键问题提出了有效的处理方案。编制了六面体实体单元的数值计算程序,并进行工程结构算例分析。结果表明,所编制程序可有效模拟实体结构的静力、动力及大变形大位移行为分析,验证了本文理论和程序的正确性和实用性。     

固体力学中的统一分析梁与有限节线法

针对现有梁理论不能解决的杆状类工程结构构件和结构体系,本文提出了一种新的梁模型——统一分析梁与一种新的数值分析方法——有限节线法。利用统一分析梁和有限节线法不仅可以解决任意杆状类结构构件或结构体系的力学分析问题,而且当问题的性质与传统梁理论的前提条件一致时,会得出同样精度的解答。算例计算结果证明了统一分析梁的合理性与有限节线法的正确性。     

基于共旋列式的空间梁的稳定性分析

基于独立于单元的共旋列式(EICR),将一种几何线性的无剪切锁死的Timoshenko梁单元扩展用于空间梁结构的几何非线性分析。考虑到三维分析中发生大转动时转动顺序的不可交换性,也即转动自由度不能作为向量采用加法规则更新,采用了四元变量来存储和更新转动自由度,使得共旋列式适用于位移任意大和转动任意大但应变很小的几何非线性分析。同时改进了Riks弧长法使之适用于带有大转动的三维几何非线性分析。给出了几个数值算例,结果表明本文方法具有较高的计算精度和效率。     

Geometrical nonlinear elasto-plastic analysis of tensegrity systems via the co-rotational method

An efficient finite element formulation is presented for geometrical nonlinear elasto-plastic analyses of tensegrity systems based on the co-rotational method. Large displacement of a space rod element is decomposed into a rigid body motion in the global coordinate system and a pure small deformation in the local coordinate system. A new form of tangent stiffness matrix, including elastic and elasto-plastic stages is derived based on the proposed approach. An incremental-iterative solution strategy in conjunction with the Newton-Raphson method is employed to obtain the geometrical nonlinear elasto-plastic behavior of tensegrities. Several numerical examples are given to illustrate the validity and efficiency of the proposed algorithm for geometrical nonlinear elasto-plastic analyses of tensegrity structures.     

有限转动张量的保角参量及在多柔体系统中的应用

采用保角转动参数描述了多体系统中的大转动张量．该方法消除了传统的欧拉参数描述所必需的约束方程，并且适于大变形部件的建模需要．利用以上结果建立了含大变形梁状部件的多体系统的力学模型．     

浅曲梁大挠度分析的余能有限元方法

以高玉臣提出的弹性大变形余能原理为基础,利用Lagrange乘子,放松平衡方程和力边界条件对余能泛函的约束,推导出广义的余能原理.根据极分解定理,将变形分为刚性转动和纯变形两部分,则余能也包含相应的两部分,一部分与刚性转动有关,而另一部分与纯变形有关.使用线弹性本构关系,建立了可用于几何非线性计算的有限元模型.应用更新的Lagrange列式法,给出了增量形式的有限元公式.数值计算结果表明,该方法可用于浅曲粱的几何大变形计算.     

Discrete element-embedded finite element model for simulation of soft particle motion and deformation

The motion and deformation of soft particles are commonly encountered and important in many applications. A discrete element-embedded finite element model (DEFEM) is proposed to solve soft particle motion and deformation, which combines discrete element and finite element methods. The collisional surface of soft particles is covered by several dynamical embedded discrete elements (EDEs) to model the collisional external forces of the particles. The particle deformation, motion, and rotation are independent of each other in the DEFEM. The deformation and internal forces are simulated using the finite element model, whereas the particle rotation and motion calculations are based on the discrete element model. By inheriting the advantages of existing coupling methods, the contact force and contact search between soft particles are improved with the aid of the EDE. Soft particle packing is simulated using the DEFEM for two cases: particle accumulation along a rectangular straight wall and a wall with an inclined angle. The large particle deformation in the lower layers can be simulated using current methods, where the deformed particle shape is either irregular in the marginal region or nearly hexagonal in the tightly packed central region. This method can also be used to simulate the deformation, motion, and heat transfer of non-spherical soft particles.     

Geometric nonlinear dynamic analysis of curved beams using curved beam element

Instead of using the previous straight beam element to approximate the curved beam,in this paper,a curvilinear coordinate is employed to describe the deformations,and a new curved beam element is proposed to model the curved beam.Based on exact nonlinear strain-displacement relation,virtual work principle is used to derive dynamic equations for a rotating curved beam,with the effects of axial extensibility,shear deformation and rotary inertia taken into account.The constant matrices are solved numerically utilizing the Gauss quadrature integration method.Newmark and Newton-Raphson iteration methods are adopted to solve the differential equations of the rigid-flexible coupling system.The present results are compared with those obtained by commercial programs to validate the present finite method.In order to further illustrate the convergence and efficiency characteristics of the present modeling and computation formulation,comparison of the results of the present formulation with those of the ADAMS software are made.Furthermore,the present results obtained from linear formulation are compared with those from nonlinear formulation,and the special dynamic characteristics of the curved beam are concluded by comparison with those of the straight beam.     

Large displacement analysis of shear deformable composite beams with interlayer slips

This paper presents a novel geometric non-linear finite element formulation for the analysis of shear deformable two-layer beams with interlayer slips. We adopt the co-rotational approach where the motion of the element is decomposed into two parts: a rigid body motion which defines a local coordinate system and a small deformational motion of the element relative to this local coordinate system. The main advantage of this approach is that the transformation matrices relating local and global quantities are independent to the choice of the geometrical linear local element. The effect of transverse shear deformation of the layers is taken into account by assuming that each layer behaves as a Timoshenko beam element. The layers are assumed to be continuously connected and partial interaction is considered by considering a continuous relationship between the interface shear flow and the corresponding slip. In order to avoid curvature and shear locking phenomena, the local linear element is formulated using “exact” displacement shape functions derived from the closed-form solution of the governing equations of a two-layer beam element. Finally, three numerical applications are presented in order to assess the performance of the proposed formulation.     

Geometrically exact curved beam element using internal force field defined in deformed configuration

This paper presents a study on the development of high-performance finite elements for geometrically nonlinear analysis of frame structures with curved members. Based on the geometrically exact beam theory, a highly efficient and accurate mixed finite element is developed. A new approach is proposed for constructing the independent internal force field by including major terms satisfying equilibrium conditions in the deformed configuration. An element-level equilibrium iteration procedure is employed for the condensation of element internal degrees of freedom during the nonlinear solution. Numerical results are presented to demonstrate the excellent performance of the element developed, and it is shown that even when each structural member is modelled with just one element, accurate solutions can still be achieved.