A spatially curved‐beam element with warping and Wagner effects

This paper presents a new spatially curved‐beam element with warping and Wagner effects that can be used for the non‐linear large displacement analysis of members that are curved in space. The non‐linear behaviour of members curved in space shows that the Wagner effects are substantial in the large twist rotation analysis. Most existing finite beam element models, such as ABAQUS and ANSYS cannot predict the non‐linear large displacement response of members curved in space correctly because the Wagner effects, viz. the Wagner moment and the corresponding finite strain terms, have not been considered in these finite beam elements. As a consequence, these finite beam elements do not provide correct predictions for the out‐of‐plane buckling and postbuckling behaviour of arches as well. In this paper, the symmetric tangent stiffness matrix has been derived based on the finite rotations parameterized by the conventional displacements. The warping and Wagner effects: both the Wagner moment and the corresponding finite strain terms and their constitutive relationship, are included in the spatially curved‐beam element. Two components of the initial curvature, the initial twist and their interactions with the displacements are also considered in the spatially curved‐beam element. This ensures that the large twist rotation analysis for the members curved in space is accurate. Comparisons with existing experimental, analytical and numerical results show that the spatially curved‐beam element is accurate and efficient for the non‐linear elastic analysis of curved members, buckling and postbuckling analysis of arches, and in its ability to predict large deflections and twist rotations in more arbitrarily curved members. Copyright © 2005 John Wiley & Sons, Ltd.     

Exact static element stiffness matrices of non‐symmetric thin‐walled curved beams

An evaluation procedure of exact static stiffness matrices for curved beams with non‐symmetric thin‐walled cross section are rigorously presented for the static analysis. Higher‐order differential equations for a uniform curved beam element are first transformed into a set of the first‐order simultaneous ordinary differential equations by introducing 14 displacement parameters where displacement modes corresponding to zero eigenvalues are suitably taken into account. This numerical technique is then accomplished via a generalized linear eigenvalue problem with non‐symmetric matrices. Next, the displacement functions of displacement parameters are exactly calculated by determining general solutions of simultaneous non‐homogeneous differential equations. Finally an exact stiffness matrix is evaluated using force–deformation relationships. In order to demonstrate the validity and effectiveness of this method, displacements and normal stresses of cantilever thin‐walled curved beams subjected to tip loads are evaluated and compared with those by thin‐walled curved beam elements as well as shell elements. Copyright © 2004 John Wiley & Sons, Ltd.     

Spatial stability and free vibration of shear flexible thin-walled elastic beams. II: Numerical approach

In a companion paper,¹ equations of motion and closed-form solutions for spatial stability and free vibration analysis of shear flexible thin-walled elastic beams were analytically derived from the linearized Hellinger–Reissner principle. In this paper, elastic and geometric stiffness matrices and consistent mass matrix for finite element analysis are evaluated by using isoparametric and Hermitian interpolation polynomials. Isoparametric interpolation functions with 2, 3 and 4 nodes per element are utilized in isoparametric beam elements, and in Hermitian beam elements, the third- and fifth-order Hermitian polynomials including shear deformation effects are newly derived and applied for the calculation of element matrices. In order to verify the validity of the finite element formulation, both analytic and numerical solutions for spatial buckling and free vibration problems including shear effects are presented and compared.     

Exact stiffness matrix of mono-symmetric composite I-beams with arbitrary lamination

The exact stiffness matrix, based on the simultaneous solution of the ordinary differential equations, for the static analysis of mono-symmetric arbitrarily laminated composite I-beams is presented herein. For this, a general thin-walled composite beam theory with arbitrary lamination including torsional warping is developed by introducing Vlasov’s assumption. The equilibrium equations and force–deformation relations are derived from energy principles. The explicit expressions for displacement parameters are then derived using the displacement state vector consisting of 14 displacement parameters, and the exact stiffness matrix is determined using the force–deformation relations. In addition, the analytical solutions for symmetrically laminated composite beams with various boundary conditions are derived as a special case. Finally, a finite element procedure based on Hermitian interpolation polynomial is developed. To demonstrate the validity and the accuracy of this study, the numerical solutions are presented and compared with the analytical solutions and the finite element results using the Hermitian beam elements and ABAQUS’s shell element.     

一种新的集成非线性杆件单元刚度矩阵的方法

对于非线性杆件单元,本文提出一种新的简便有效的集成单元刚度矩阵的算法。该方法直接从结构力学中的位移法的概念出发,通过解析积分或数值积分求解积分算子,由积分算子线性组合,能快速求解考虑弯、剪、扭、轴压等各种非线形刚度的杆件单元的刚度矩阵。该方法具有广泛的普适性,能适用于所有多项式、插值多项式、解析式、离散点描述的变刚度、变截面直杆的单元刚度矩阵集成计算。文中通过求解线性直杆单元刚度矩阵验证了该方法的正确性。     

Eigenvalue analysis of compatible and incompatible rectangular four-node quadrilateral elements

Exact analytical expressions for the eigenvalues of the elastic stiffness matrix are obtained for the four-node, rectangular, quadrilateral element. A procedure is given for identifying alternative hourglass modes and eigenvalues which render the element incompatible but with non-monotonic convergence assured. A convergence study confirms that for the special case of when the hourglass modes coincide with beam bending the element can serve as a beam element. Analytical expressions are given for the resulting element stiffness matrix.     

A new approach for displacement functions of a curved Timoshenko beam element in motions normal to its initial plane

In existing literature, either analytical methods or numerical methods, the formulations for free vibration analysis of circularly curved beams normal to its initial plane are somewhat complicated, particularly if the effects of both shear deformation (SD) and rotary inertia (RI) are considered. It is hoped that the simple approach presented in this paper may improve the above‐mentioned drawback of the existing techniques. First, the three functions for axial (or normal to plane) displacement and rotational angles about radial and circumferential (or tangential) axes of a curved beam element were assumed. Since each function consists of six integration constants, one has 18 unknown constants for the three assumed displacement functions. Next, from the last three displacement functions, the three force–displacement differential equations and the three static equilibrium equations for the arc element, one obtained three polynomial expressions. Equating to zero the coefficients of the terms in each of the last three expressions, respectively, one obtained 17 simultaneous equations as functions of the 18 unknown constants. Excluding the five dependent ones among the last 17 equations, one obtained 12 independent simultaneous equations. Solving the last 12 independent equations, one obtained a unique solution in terms of six unknown constants. Finally, imposing the six boundary conditions at the two ends of an arc element, one determined the last six unknown constants and completely defined the three displacement functions. By means of the last displacement functions, one may calculate the shape functions, stiffness matrix, mass matrix and external loading vector for each arc element and then perform the free and forced vibration analyses of the entire curved beam. Good agreement between the results of this paper and those of the existing literature confirms the reliability of the presented theory. Copyright © 2005 John Wiley & Sons, Ltd.     

一种有效的分析任意空间形状曲杆单元的位移函数

利用经典弹性理论和微分几何、矩阵方法等数学理论,基于空间自然坐标系和随体坐标系,通过求解应变与位移之间关系的微分方程,得到了一种能完全反映任意空间形状圆截面曲杆单元刚体位移和常应变等模式的位移函数。给出了两个算例,通过将采用导出的位移函数建立的有限元解与解析解进行了比较,验证了它的正确性,同时通过将基于位移函数导出的有限元解与应用软件得到的解进行了比较,其计算精度,特别是应力计算精度大为改善,验证了它的有效性。由于计算效率高,提出的曲杆单元可望在三维大曲率井的钻柱非线性分析等工程实际中发挥作用。     

A semi‐analytical curved element for linear elasticity based on the scaled boundary finite element method

This work introduces a semi‐analytical formulation for the simulation and modeling of curved structures based on the scaled boundary finite element method (SBFEM). This approach adapts the fundamental idea of the SBFEM concept to scale a boundary to describe a geometry. Until now, scaling in SBFEM has exclusively been performed along a straight coordinate that enlarges, shrinks, or shifts a given boundary. In this novel approach, scaling is based on a polar or cylindrical coordinate system such that a boundary is shifted along a curved scaling direction. The derived formulations are used to compute the static and dynamic stiffness matrices of homogeneous curved structures. The resulting elements can be coupled to general SBFEM or FEM domains. For elastodynamic problems, computations are performed in the frequency domain. Results of this work are validated using the global matrix method and standard finite element analysis. Copyright © 2016 John Wiley & Sons, Ltd.     

Curved beam elements with penalty relaxation

Shallowly curved beam elements, including shear deformation and rotary inertia effects, are derived from Hamilton's variational principle. Different degree polynomials, labelled ‘anisoparametric’, are used to interpolate the kinematic variables, instead of uniform interpolations as in the conventional isoparametric procedure. This approach yields a correct representation of the bending strain and, importantly, the membrane and transverse shear strains. Consequently, the severe shortcomings of the exactly integrated isoparametric elements, characterized by excessively stiff solutions in the thin regime (a phenomenon often referred to as membrane and shear locking), are overcome. Uniform (isoparametric-like) nodal patterns are achieved by explicitly enforcing higher-degree penalty modes in the membrane and shear strains. This procedure preserves the compatibility of the kinematic field and the capability of the element to move rigidly without straining. Exact quadratures are used on all element matrices, producing a correct rank stiffness matrix, a consistent load vector and a consistent mass matrix. The elements suffer no limitations over the entire theoretical range of the slenderness ratio. For further enhancement and, particularly, in coarse-mesh situations, an effective relaxation of penalty constraints at the local element level is introduced. This technique ensures a well-conditioned stiffness matrix. Although the element penalty constraints are relaxed, the corresponding global structure constraints are enforced as is required by the analytic theory. Particular attention is given to the simplest element—a two-node, six degree-of-freedom beam in which all strains are constant. Solutions to static and free vibration arch and ring problems are presented, demonstrating the exceptional modelling capabilities of this element.     

考虑构件弹性的电动Stewart并联机构刚度建模与仿真

基于集总解析建模方法和构件有限元分析建立包含驱动副、被动万向铰链和运动杆件弹性变形以及预载作用下的Stewart机构刚度矩阵模型。采用添加虚拟铰链等效构件弹性的方式,将分支等效为一系列刚性构件经由主、被动副以及虚拟铰链连接的形式,给出了运动关节和虚拟铰链变量对机构末端位姿的运动学Jacobian矩阵的数值计算方法,应用虚功原理得到静平衡方程,最终建立了机构无预载以及预载下的刚度矩阵模型。该模型不仅考虑了控制环路刚度,还将构件柔性的有限元分析结果与解析建模相结合,在降低计算成本的同时保证了精度。通过一机构分析实例,考察了两种模型下刚度分布的差异。     

Free vibration analysis of arches using curved beam elements

The natural frequencies and mode shapes for the radial (in‐plane) bending vibrations of the uniform circular arches were investigated by means of the finite arch (curved beam) elements. Instead of the complicated explicit shape functions of the arch element given by the existing literature, the simple implicit shape functions associated with the tangential, radial (or normal) and rotational displacements of the arch element were derived and presented in matrix form. Based on the relationship between the nodal forces and the nodal displacements of a two‐node six‐degree‐of‐freedom arch element, the elemental stiffness matrix was derived, and based on the equation of kinetic energy and the implicit shape functions of an arch element the elemental consistent mass matrix with rotary inertia effect considered was obtained. Assembly of the foregoing elemental property matrices yields the overall stiffness and mass matrices of the complete curved beam. The standard techniques were used to determine the natural frequencies and mode shapes for the curved beam with various boundary conditions and subtended angles. In addition to the typical circular arches with constant curvatures, a hybrid beam constructed by using an arch segment connected with a straight beam segment at each of its two ends was also studied. For simplicity, a lumped mass model for the arch element was also presented. All numerical results were compared with the existing literature or those obtained from the finite element method based on the conventional straight beam element and good agreements were achieved. Copyright © 2003 John Wiley & Sons, Ltd.     

Efficient laminated composite beam element subjected to variable axial force for coupled stability analysis

A numerically efficient laminated composite beam element subjected to a variable axial force is presented for a coupled stability analysis. The analytical technique is used to present the thin-walled laminated composite beam theory considering the transverse shear and the restrained warping-induced shear deformation based on an orthogonal Cartesian coordinate system. The elastic strain energy and the potential energy due to the variable axial force are introduced. The equilibrium equations are derived from the energy principle, and explicit expressions for the displacement parameters are presented using the power series expansions of displacement components. Finally, the member stiffness matrix is determined using the force–displacement relations. In order to verify accuracy and efficiency of the beam element developed in this study, numerical results are presented and compared with results from other researchers and the finite beam element results, and the detailed finite shell element analysis results using ABAQUS; especially, the influence of variable axial forces, the fiber orientation, and boundary conditions on the buckling behavior of the laminated composite beams is parametrically investigated.     

Power series expansion of the general stiffness matrix for beam elements

The general stiffness matrix for a beam element is derived from the Bernoulli–Euler differential equation with the inclusion of axial forces. The terms of this matrix are expanded into a power series as a function of the two variables: the axial force, and; the vibrating frequency. It is shown that the first three terms of the resulting series, which are derived in the technical literature from assumed static displacement functions, correspond respectively to the elastic stiffness matrix, the consistent mass matrix, and the geometric matrix. Higher order terms up to the second order terms of the series expansion are obtained explicitly. Also a discussion is presented for establishing the region of convergence of the series expansion for the dynamic stiffness matrix, the stability matrix, and the general stiffness matrix.     

基于静力凝聚的高精度变截面梁单元及其几何非线性分析研究

基于Euler-Bernoulli梁单元基本假定,通过静力凝聚获得截面特性沿单元轴向连续变化的变截面梁单元高精度刚度矩阵,并提出一种基于随动坐标法求解变截面梁杆结构大位移、大转动、小应变问题的新思路。首先依据插值理论和非线性有限元理论推导出三节点变截面梁单元的切线刚度矩阵,然后使用静力凝聚方法消除中间节点自由度,从而得到一种新型非线性两节点变截面梁单元。结合随动坐标法,在变形后位形上建立随动坐标系,得到变截面梁单元的大位移全量平衡方程。实例计算表明,该新型变截面梁单元具有较高的计算精度,可应用于变截面梁杆系统大位移几何非线性分析。     

An accurate shear-deformable six-node triangular plate element for laminated composite structures

A six-node triangle plate/shell element is developed for the analysis of laminated composite structures. This model is formulated using Hamilton's principle along with a first-order (Reissner/Mindlin) shear deformation theory. The element is based upon an isoparametric representation along with an interdependent interpolation strategy; bicubic polynomials for the transverse displacement and biquadratic polynomials for the element geometry, in-plane displacements and rotations. The resulting element, which is evaluated using exact numerical integration, has correct rank and is free of shear ‘locking’. Numerical results are presented that validate the new element and prove its outstanding convergence capabilities in comparison to existing triangular elements using standardized test problems (elastic eigenvalue analysis, patch test, static simply supported square-plate solutions) and experimentally measured vibration data of cantilevered isotropic and composite plates.     

Instability of laminated composite thin-walled open-section beams

Instability of thin-walled open-section laminated composite beams is studied using the finite element method. A two-noded, 8 df per node thin-walled open-section laminated composite beam finite element has been used. The displacements of the element reference axis are expressed in terms of one-dimensional first order Hermite interpolation polynomials, and line member assumptions are invoked in formulation of the elastic stiffness matrix and geometric stiffness matrix. The nonlinear expressions for the strains occurring in thin-walled open-section beams, when subjected to axial, flexural and torsional loads, are incorporated in a general instability analysis. Several problems for which continuum solutions (exact/approximate) are possible have been solved in order to evaluate the performance of finite element. Next its applicability is demonstrated by predicting the buckling loads for the following problems of laminated composites: (i) two layer (45°/−45°) composite Z section cantilever beam and (ii) three layer (0°/45°/0°) composite Z section cantilever beam.     

几何非线性高性能复合材料筋混凝土梁Heterosis组合壳单元

对于高性能碳纤维增强聚合物复合材料(CFRP)筋混凝土梁,研究几何非线性组合壳单元模型,对预应力CFRP筋混凝土梁进行了全过程分析。引入Von Karman理论,推导了局部坐标系下Piola2Kirchhoff 应力矩阵和几何刚度矩阵;分别采用组合壳单元和分层壳单元模拟预应力CFRP 筋和玻璃纤维增强聚合物复合材料(GFRP)筋,并推导了CFRP筋对组合壳单元刚度矩阵的贡献,同时采用Heterosis选择积分技术以避免剪切锁定和零能量模式,研制了相应的非线性计算程序。计算结果与试验数据对比可知,挠度发展规律和预应力CFRP筋应变发展规律均吻合良好,说明了研究单元的有效性及研制程序的正确性;CFRP筋具有高强度性能,梁试件破坏时CFRP筋均未失效;利用预应力CFRP筋应变重分布系数研究了梁的刚度退化规律,表明采用GFRP筋代替普通钢筋在加载后期会使梁的刚度退化减小。      

Exact solutions for extensible circular curved timoshenko beams with nonhomogeneous elastic boundary conditions

Summary A generalized Green function ofnth-order ordinary differential equation with forcing function composed of the delta function and its derivatives is obtained. The generalized Green function can be easily and effectively applied to both the boundary value problems and the initial value problems. The generalized Green function is expressed in terms ofn linearly independent normalized homogeneous solutions. It is the generalization of those given by Pan and Hohenstein, and Kanwal. Accordingly, the exact solution for static analysis of an extensible circular curved Timoshenko beam with general nonhomogeneous elastic boundary conditions, subjected to any transverse, tangential and moment loads is obtained. The three coupled governing differential equations are uncoupled into one complete sixth-order ordinary differential characteristic equation in the tangential displacement. The explicit relations between the angle of rotation due to bending, the transverse displacement and the tangential displacement are obtained. The deflection curves due to a unit generalized displacement at nodal coordinate, and the exact element stiffness matrix are derived based on the solution for the general system. A finite element method can be developed based on the results for the dynamic analysis. Meanwhile, the stiffness locking phenomena accompanied in some other curved beam element methods does not exist in the proposed method.