THE GENERATION OF NON－LINEAR STIFFNESS MATRIX OF TRIANGLE ELEMENT WHENCONSIDERING BOTH THE BENDING AND IN－PLANEME MBRANE FORCES

Using Stricklin Melhod￣[5],we have this paper has derived the formulas for the ge-neration of non-linear element stiffness matrix of a triangle element when considering both the bending and the in-plane membrane forces. A computer programme for the calculation of large deflection and inner forces of shallow shells is designed on theseformulas. The central deflection curve computed by this programme is compared with other pertaining results.

Finite element model for arch bridge vibration dynamics considering effect of suspender length adjustment on geometry stiffness matrix

In this paper, we established a finite element (FEM) model to analyze the dynamic characteristics of arch bridges. In this model, the effects of adjustment to the length of a suspender on its geometry stiffness matrix are stressed. The FEM equations of mechanics characteristics, natural frequency and main mode are set up based on the first order matrix perturbation theory. Applicantion of the proposed model to analyze a real arch bridge proved the improvement in the simulation precision of dynamical characteristics of the arch bridge by considering the effects of suspender length variation.

THE DYNAMIC STIFFNESS MATRIX OF THE FINITE ANNULAR PLATE ELEMENT

The dynamic deformation of harmonic vibration is used as the shape functions of the finite annular plate element, and sonic integration difficulties related to the Bessel's functions are solved in this paper. Then the dynamic stiffness matrix of the finite annular plate element is established in closed form and checked by the direct stiffness method. The paper has given wide convcrage for decomposing the dynamic matrix into the power series of frequency square. By utilizing the axial symmetry of annular elements, the modes with different numbers of nodal diameters at s separately treated. Thus some terse and complete results are obtained as the foundation of structural characteristic analysis and dynamic response compulation.

New tangent stiffness matrix for geometrically nonlinear analysis of space frames

A three-dimensional beam element is derived based on the principle of stationary total potential energy for geometrically nonlinear analysis of space frames. A new tangent stiffness matrix, which allows for high order effects of element deformations, replaces the conventional incremental secant stiffness matrix. Two deformation stiffness matrices due to the variation of axial force and bending moments are included in the tangent stiffness. They are functions of element deformations and incorporate the coupling among axial, lateral and torsional deformations. A correction matrix is added to the tangent stiffness matrix to make displacement derivatives equivalent to the commutative rotational degrees of freedom. Numerical examples show that the proposed dement is accurate and efficient in predicting the nonlinear behavior, such as axial-torsional and flexural-torsional buckling, of space frames even when fewer elements are used to model a member.

New Formula for Geometric Stiffness Matrix Calculation

The standard formula for geometric stiffness matrix calculation, which is convenient for most engineering applications, is seen to be unsatisfactory for large strains because of poor accuracy, low convergence rate, and stability. For very large compressions, the tangent stiffness in the direction of the compression can even become negative, which can be regarded as physical nonsense. So in many cases rubber materials exposed to great compression cannot be analyzed, or the analysis could lead to very poor convergence. Problems with the standard geometric stiffness matrix can even occur with a small strain in the case of plastic yielding, which eventuates even greater practical problems. The authors demonstrate that amore precisional approach would not lead to such strange and theoretically unjustified results. An improved formula that would eliminate the disadvantages mentioned above and leads to higher convergence rate and more robust computations is suggested in this paper. The new formula can be derived from the principle of virtual work using a modified Green-Lagrange strain tensor, or from equilibrium conditions where in the choice of a specific strain measure is not needed for the geometric stiffness derivation (which can also be used for derivation of geometric stiffness of a rigid truss member). The new formula has been verified in practice with many calculations and implemented in the RFEM and SCIA Engineer programs. The advantages of the new formula in comparison with the standard formula are shown using several examples.

Systematic Elastostatic Stiffness Model of Over-Constrained Parallel Manipulators Without Additional Constraint Equations

The establishment of an elastostatic stiffness model for over constrained parallel manipulators(PMs),particularly those with over constrained subclosed loops,poses a challenge while ensuring numerical stability.This study addresses this issue by proposing a systematic elastostatic stiffness model based on matrix structural analysis(MSA)and independent displacement coordinates(IDCs)extraction techniques.To begin,the closed-loop PM is transformed into an open-loop PM by eliminating constraints.A subassembly element is then introduced,which considers the flexibility of both rods and joints.This approach helps circumvent the numerical instability typically encountered with traditional constraint equations.The IDCs and analytical constraint equations of nodes constrained by various joints are summarized in the appendix,utilizing multipoint constraint theory and singularity analysis,all unified within a single coordinate frame.Subsequently,the open-loop mechanism is efficiently closed by referencing the constraint equations presented in the appendix,alongside its elastostatic model.The proposed method proves to be both modeling and computationally efficient due to the comprehensive summary of the constraint equations in the Appendix,eliminating the need for additional equations.An example utilizing an over constrained subclosed loops demonstrate the application of the proposed method.In conclusion,the model proposed in this study enriches the theory of elastostatic stiffness modeling of PMs and provides an effective solution for stiffness modeling challenges they present.

Stiffness Matrix Derivation of Space Beam Element at Elevated Temperature

Element stiffness equation is very important in structural analysis, and directly influences the accuracy of the results. At present, derivation method of element stiffness equation is relatively mature under ambient temperature, and the elastic phrase of material stress-strain curve is generally adopted as physical equation in derivation. However, the material stress-strain relationship is very complicated at elevated temperature, and its form is not unique, which brings great difficulty to the derivation of element stiffness equation. Referring to the derivation method of element stiffness equation at ambient temperature, by using the continuous function of stress-strain-temperature at elevated temperature, and based on the principle of virtual work, the stiffness equation of space beam element and the formulas of stiffness matrix are derived in this paper, which provide basis for finite element analysis on structures at elevated temperature.

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.

A Finite Element Method for Cracked Components of Structures

In this paper, a method is developed for determining the effective stiffness of the cracked component. The stiffness matrix of the cracked component is integrated into the global stiffness matrix of the finite element model of the global platform for the FE calculation of the structure in any environmental conditions. The stiffness matrix equation of the cracked component is derived by use of the finite variation principle and fracture mechanics. The equivalent parameters defining the element that simulates the cracked component are mathematically presented, and can be easily used for the FE calculation of large scale cracked structures together with any finite element program. The theories developed are validated by both lab tests and numerical calculations, and applied to the evaluation of crack effect on the strength of a fixed platform and a self-elevating drilling rig.

GEOMETRICALLY NONLINEAR FE FORMULATIONS FOR THE MACRO-ELEMENT UNIPLET OF FOLDABLE STRUCTURES

Geometrically nonlinear stiffness matrix due to large displacement small strain was firstly formulated explicitly for the basic components of pantographic foldable structures,namely, the uniplet, derived from a three node beam element.The formulation of the uniplet stiffness matrix is based on the precise nonlinear finite element theory and the displacement harmonized and internal force constraints are applied directly to the deformation modes of the three node beam element. The formulations were derived in general form, and can be simplified for particular foldable structures, such as flat, cylindrical and spherical structures.Finally, two examples were presented to illustrate the applications of the stiffness matrix evolved.

A CATENARY ELEMENT FOR THE ANALYSIS OF CABLE STRUCTURES

Based on analytical equations, a cat ena ry element is presented for the finite element analysis of cable structures. Com pared with usually used element(3_node element, 5_node element), a program with the proposed element is of less computer time and better accuracy.

Time-varying stiffness analysis of double-row tapered roller bearing based on the mapping structure of bearing stiffness matrix

Time-varying stiffness is one of the most important dynamic characteristics of rolling element bearings.The method of analyzing the elements in the bearing stiffness matrix is usually adopted to investigate the characteristics of bearing stiffness.Linear mapping structure of the bearing stiffness matrix is helpful to understand the varying compliance excitation and its influence on vibration transmission.In this study,a method to analyze the mapping structure of bearing stiffness matrix is proposed based on the singular value decomposition of block matrices in the stiffness matrix.Not only does this method have the advantages of coordinate transformation independence and unit independence,but also the analysis procedure involved is geometrically intuitive.The time-varying stiffness matrix of double-row tapered bearing is calculated and analyzed using the proposed method under two representative load cases.The principal stiffnesses and principal axes defined in the method together indicate the dominant and insignificant stiffness properties with the corresponding directions,and the vibration transmission properties are also revealed.Besides,the coupling behaviors between different shaft motions are found during the analysis of mapping structure.The mechanism of the generation of varying compliance excitation is also revealed.

Gershgorin and Rayleigh Bounds on the Eigenvalues of the Finite-Element Global Matrices via Optimal Similarity Transformations

The large finite element global stiffness matrix is an algebraic, discreet, even-order, differential operator of zero row sums. Direct application of the, practically convenient, readily applied, Gershgorin’s eigenvalue bounding theorem to this matrix inherently fails to foresee its positive definiteness, predictably, and routinely failing to produce a nontrivial lower bound on the least eigenvalue of this, theoretically assured to be positive definite, matrix. Considered here are practical methods for producing an optimal similarity transformation for the finite-elements global stiffness matrix, following which non trivial, realistic, lower bounds on the least eigenvalue can be located, then further improved. The technique is restricted here to the common case of a global stiffness matrix having only non-positive off-diagonal entries. For such a matrix application of the Gershgorin bounding method may be carried out by a mere matrix vector multiplication.

基于共旋坐标法的几何非线性裂纹梁单元

现有有限元软件模拟裂纹结构主要采用指派裂纹的实体和壳单元,但缺少梁单元。因此,在弹性梁单元的曲率函数中引入Dirac-δ函数表征局部损伤的曲率变化,得到位移分布函数;推导出应变矩阵和单元刚度矩阵,提出了裂纹梁单元。基于共旋坐标法理论进一步发展了考虑裂纹损伤的几何非线性裂纹梁单元,给出了单元的具体算法及流程。利用课题组开发的有限元求解框架,在程序中实现了裂纹梁单元及共旋列式裂纹梁单元。经过典型算例,对比分析了裂纹位置、裂纹深度及外荷载对裂纹梁弯曲性能的影响。结果表明:该单元能够准确描述裂纹结构的挠度变化,同时与现有模拟裂纹结构弯曲变形的单元相比,提出的裂纹梁和共旋列式裂纹梁在精度分别提高了约1%和4%,并且能够退化成无损伤的梁单元和几何非线性梁单元。

基于Feon-Abaqus的杆梁混合结构的有限元分析与应用

由于不同类型单元的结合处其结点自由度不一致,在杆梁混合结构的内力分析过程中整体刚度矩阵的集成存在效率较低、存储空间较大等问题.Feon框架采用的集成方式在运算速度和存储经济性等方面仍存在改进空间,且Feon目前仅能处理线性情形.Abaqus虽能处理各种非线性问题,但缺少针对混合结构的参数化建模解决方法.将Feon与Abaqus进行结合,可以为杆梁混合结构的线性及非线性有限元分析提供更为有效的方法.该文首先改进Feon中整体刚度矩阵的集成方法,新方法显著减小了计算中的过渡矩阵和最终刚度矩阵规模,提高了处理速度和存储经济性.其次基于对Feon框架的改进,提供了在Python平台上直接对杆梁混合结构进行有限元分析的方法,将处理范围由线性拓展到非线性,并编写了相应的Feon-Abaqus转换程序,提供了Abaqus处理混合结构参数化建模解决方法.最后通过算例验证了改进方法在确保计算精度的同时能提高计算速度,验证了基于Feon-Abaqus在杆梁混合结构线性和非线性分析中的有效性.

有限元总体刚度矩阵集成的教学

单元刚度矩阵组装到总体刚度矩阵的原理和方法是有限元教学中基本而重要的内容,因涉及的节点平衡方程问题简单而讲解不够深入,容易导致留下认知隐患。文中介绍了有限元总体刚度矩阵集成的两种教学设计,一种涉及节点平衡方程,另一种不涉及节点平衡方程。针对由节点平衡方程引出的似是而非的认知问题展开教学。教学过程中采用提问、讨论、推导、动画演示等方式以加深学生的认知和理解。

一种复合驱动并联机构刚度分析

针对直升机旋翼支撑平台,设计了一种复合驱动四自由度并联机构,利用复合驱动分支增大了并联机构工作空间。推导了并联机构中各驱动分支和动平台之间的运动关系式,考虑约束力、驱动力和驱动力矩综合作用下复合驱动分支的横向变形、轴向变形和扭转变形,结合机构的运动学方程推导了复合驱动并联机构的刚度矩阵。利用有限元软件验证了并联机构整体刚度矩阵的正确性。仿真结果表明该并联机构的刚度满足使用需要,为含复合驱动并联机构的应用和推广奠定了理论基础。

考虑螺栓预紧力作用范围不均影响的结合面刚度识别

螺栓连接结合面的刚度和阻尼的变化会引起整个螺栓连接结构的动态特性变化,准确得到结合面动态参数在工程中具有极大的实用价值。基于螺栓连接结构,提出螺栓预紧力作用范围内结合面刚度分布不均的结合面等效改进模型,并对有限元建模中不同数量刚度矩阵单元分布作出分析。采用实验与有限元分析相结合的方法对螺栓连接结合面进行刚度参数辨识。结果表明:在一定程度上,增加刚度矩阵单元数量可以提高固有频率解的精度。同时考虑螺栓预紧力作用范围内结合面刚度分布不均,可以有效地提高螺栓连接结构的等效建模精度。

三维杆系结构的几何非线性有限元分析

为了更准确地描述杆系结构的几何非线性性能,建立了一种基于三维梁单元有限元分析的计算方法。引入了考虑两方向曲率和扭转角变化的坐标转换矩阵来描述任意增量下的单元平移和转动;采用了包括轴向变形和扭转的非线性项的刚度矩阵来考虑高阶非线性项的影响。应用广义位移控制法进行增量迭代,编制了相应的三维梁单元非线性计算程序NL_Beam3D。通过对几个例子进行的分析,验证了该方法可较好地考虑结构几何非线性。

求解动荷载作用下多层粘弹性半空间轴对称问题的精确刚度矩阵法

利用Laplace-Hankel联合积分变换,推导出了单层粘弹性半空间轴对称问题在动荷载作用下,层间完全接触情况的刚度矩阵,然后按传统的有限元方法组成总体刚度矩阵。通过求解由总体刚度矩阵所构成的代数方程就可解出动荷载作用下多层粘弹性半空间轴对称问题的矩阵。由于刚度矩阵的元素中只含有负指数项,计算时不会出现溢出的现象。本文还成功地应用了Durbin的Laplace逆变换的数值方法,求解出了多层粘弹性体的时域解。最后,文中还给出了路面动弯沉的计算结果与实测结果的对比来证明推导结果的准确性。