直接基于单元平衡的变截面梁反应分析方法

固定形状的单元位移插值函数不能合理地近似变截面梁内部的位移变化,从而影响了传统梁单元用于计算变截面梁的精度.采用直接基于单元平衡的思想给出了计算变截面梁反应的有限元方法,解决了单元位移插值函数局限性所带来的问题.导出了变截面梁单元的单元刚度矩阵、单元等效节点荷载和单元一致质量矩阵.在此基础上,利用编制的程序进行了算例验证与分析.算例验证了本文理论的正确性,表明本文方法具有很高的计算精度.     

一种考虑剪切变形的平行四边形厚/薄板通用单元

根据Timoshenko二广义位移梁理论，构造了深梁位移场的插值函数。利用斜坐标系与直角坐标系的变换关系、有限条带思想和深梁位移插值函数，构造了一种考虑剪切变形的平行四边形厚／薄板弯曲通用单元的位移(曲率、剪应变、转角、横向位移)插值函数，导出了刚度矩阵和非结点荷载等效力。并对简支阍支方板、Razzaque斜板、四边简支斜交板弯曲进行了数值计算。算例表明此单元有较好的精度，对于薄板不出现剪切闭锁，可适应于目前桥梁建设中大量采用的斜交板桥结构分析。     

一个改进的平面梁单元

根据有限单元法基本原理 ,提出了一个变截面平面梁单元 ,推导了其单元钢度矩阵。这一改进的梁单元用于分析梁高呈线性变化及二次抛物线变化的矩形截面梁 ,将得到准确解。文中给出了一个变截面悬臂梁算例 ,计算表明 ,这一改进的梁单元使变截面梁的分析大大简化     

变截面Timoshenko梁的单元刚度矩阵

变截面构件在工程中应用广泛,在对变截面梁进行数值计算时,需要建立变截面梁单元的刚度矩阵。该文采用势能驻值原理,考虑了轴力引起的几何非线性和剪切变形的影响,将梁截面刚度的变化率作为小量,得到了近似到二阶的单元刚度矩阵。在构造位移模式时,从梁的微分平衡方程出发,得到同样近似到二阶、分别以三次和五次多项式表示的剪切和弯曲位移模式。该文还证明了单元刚度矩阵的奇异性,给出了轴压刚度的表达式,定量论证了与某些精确解的误差,表明在一定范围内,该文的结果具有足够的精度。最后以一个计算实例说明该文的单元刚度矩阵具有较快的收敛性。     

变截面梁热变形位移的计算

变截面梁是工程中广泛应用的结构元件,温度的变化将使梁产生应力和位移。随着变截面梁弯曲问题的解决,该类问题也就不难处理。本文将介绍变截面梁热变形位移的计算公 ...     

中性线修正型变截面梁类构件压电控制

传统绝对节点坐标法(absolute nodal coordinate formulation,ANCF)在变截面梁类构件建模过程中常以几何中位线等效构造单元中性线,难以对变截面单元位移场状态进行精确描述.为解决此类问题,本文以中细型变截面梁类构件为研究对象,深入考虑变截面结构几何因素及复合材料属性对变截面梁类构件中性...     

等截面梁有限变形的传递函数解法

本文应用传递函数方法对等截面梁的有限变形进行了分析,对于等截面梁的有限变形问题,该方法从变分方程出发把问题表述为状态空间的形式,然后利用Gauss积分对轴力进行加速迭代求解,不需要进行增量迭代即可取得具有良好计算精度的数值结果,对简单受力的等截面梁情况该解可以看作是所讨论问题的精确解。对于受力比较复杂或者阶梯变截面梁情况,为减少运算量,可以和有限元法类似,采用多个单元进行拼接,从而得到问题的解。数值算例表明,本方法具有半解析、精度高、收敛快等特点。     

考虑变形耦合的几何非线性空间梁单元

以"精确几何模型梁单元"为代表的很多几何非线性梁单元,在构造过程中分别对描述截面转动的转角和描述截面形心位置的位移进行了独立插值,由此引起了诸如运动学描述冗余和剪切闭锁等困难。其根本原因在于单元形函数没能体现细长梁中的变形耦合关系。本文对这类传统单元进行了改造,通过深入研究单元变形之间的内在联系,提出了一种变形场完全满足Bernoulli梁变形耦合关系的新单元,避免了构造过程中对转动矢量的插值,并通过数值算例检验了单元的有效性。     

柔性梁大变形计算的一种新方法

研究了柔性梁大变形问题。常规Lagrangian有限元格式在处理大变形问题时,由于其单元插值函数不满足位移场的协调性要求,从而需要划分较多的单元,才能得到较好的结果。本文首先推导了Lagrangian坐标描述下的位移场变量满足的协调关系式,利用此关系式给出了位移场协调的非线性单元插值函数。基于虚功原理导出了梁大变形问题的非线性控制方程,数值计算结果证明了本文方法的正确性和有效性。     

钢框架高等分析中考虑局部屈曲的子结构法

将相对自由度概念引入三维实体等参数单元,得到了16结点相对自由度壳元的插值表达式,并基于更新拉格朗日构型的虚位移原理推导了该壳单元的刚度矩阵。在单元刚度矩阵形成以后采用子结构法对内部结点自由度进行聚缩,根据薄壁构件相关假定并考虑截面翘曲因素的影响,通过边界结点的位移协调条件把壳单元组装成每端7个自由度的薄壁梁-柱单元参与整体分析,可以考虑板件局部屈曲的影响;突破了传统分析的紧凑型截面限制,并通过算例证明了本文方法的精确度和有效性。     

Refined beam elements with only displacement variables and plate/shell capabilities

This paper proposes a refined beam formulation with displacement variables only. Lagrange-type polynomials, in fact, are used to interpolate the displacement field over the beam cross-section. Three- (L3), four- (L4), and nine-point (L9) polynomials are considered which lead to linear, quasi-linear (bilinear), and quadratic displacement field approximations over the beam cross-section. Finite elements are obtained by employing the principle of virtual displacements in conjunction with the Unified Formulation (UF). With UF application the finite element matrices and vectors are expressed in terms of fundamental nuclei whose forms do not depend on the assumptions made (L3, L4, or L9). Additional refined beam models are implemented by introducing further discretizations over the beam cross-section in terms of the implemented L3, L4, and L9 elements. A number of numerical problems have been solved and compared with results given by classical beam theories (Euler-Bernoulli and Timoshenko), refined beam theories based on the use of Taylor-type expansions in the neighborhood of the beam axis, and solid element models from commercial codes. Poisson locking correction is analyzed. Applications to compact, thin-walled open/closed sections are discussed. The investigation conducted shows that: (1) the proposed formulation is very suitable to increase accuracy when localized effects have to be detected; (2) it leads to shell-like results in case of thin-walled closed cross-section analysis as well as in open cross-section analysis; (3) it allows us to modify the boundary conditions over the cross-section easily by introducing localized constraints; (4) it allows us to introduce geometrical boundary conditions along the beam axis which lead to plate/shell-like cases.     

埃尔米特梁单元的分块对角与高阶质量矩阵

埃尔米特梁单元常用的集中质量矩阵,是由挠度自由度对应的一致质量矩阵元素通过行求和或节点积分构造。然而,数值结果表明该集中质量矩阵在求解包含自由端的梁振动问题时,会出现频率精度掉阶现象。本文首先从保障质量矩阵最优收敛性的数值积分精度出发,分别针对三次和五次梁单元,发展了质量矩阵的梯度增强节点积分方案。利用梯度增强节点积分方案,可以得到具有分块对角形式的单元质量矩阵,而其组装的整体质量矩阵除边界节点外仍然呈现对角形式。对于两种单元,其分块对角质量矩阵分别具有4阶最优精度和6阶次优精度。再者,将标准一致质量矩阵和具有同阶精度的梯度增强节点积分质量矩阵进行优化组合,建立了具有超收敛特性的高阶质量矩阵。最后,通过数值算例系统验证了三次和五次单元的分块对角与高阶质量矩阵的频率计算精度。     

Stiffness matrices of thin-walled composite beam with mono-symmetric I- and channel-sections on two-parameter elastic foundation

For the coupled analysis of thin-walled composite beam under the initial axial force and on two-parameter elastic foundation with mono-symmetric I- and channel-sections, the stiffness matrices are derived. The stiffness matrices developed by this study are based on the homogeneous forms of simultaneous ordinary differential equations using the eigen-problem. For this, from the elastic strain energy, the potential energy due to the initial axial force and the strain energy considering the foundation effects, the equilibrium equations and force–displacement relationships are derived. The exact displacement functions for displacement parameters are evaluated by determining the eigenmodes corresponding to multiple non-zero and zero eigenvalues. Then the element stiffness matrix is determined using the force–displacement relationships. For the purpose of comparison, the finite element model based on the classical Hermitian interpolation polynomial is presented. In order to verify the accuracy and the superiority of the beam elements developed herein, the numerical solutions are presented and compared with results from the Hermitian beam elements and the ABAQUS’s shell elements. Particularly, the influence of the initial compressive and tensile forces, the fiber orientation, and the boundary conditions on the coupled behavior of composite beam with mono-symmetric I- and channel-sections is parametrically investigated.     

Variational base and solution strategies for non-linear force-based beam finite elements

This paper presents the variational bases for the non-linear force-based beam elements. The element state determination of these elements is obtained exactly from a two-field functional with independent stress and strain fields. The variational base of the non-linear force-based beam elements implemented in a general purpose displacement-based finite element program requires the inclusion of independent displacement field in the formulation. For this purpose, a three-field functional is considered with independent displacement, stress, and strain fields. Various local and global solution strategies come out from the mixed formulation of the beam element, and these are shown to yield the algorithms presented for non-linear force formulation beam elements in literature; thus removing any doubts on their variational bases. The presented numerical examples demonstrate the accuracy and robustness of the solution algorithms adapted for mixed formulation elements over popularly used displacement-based beam finite elements even for large structural systems.     

DISPLACEMENT CALCULATION FOR STRAIGHT BAR WITH VARIABLE STIFFNESS BASED ON GRAPH MULTIPLICATION METHOD1)

为解决图乘法计算变截面杆件位移时分割图形数目多,计算容易出错这一问题,提出了一种刚度折算的方法,把长度已知的变截面杆件等价转换为两个一定长度的等截面杆件,通过刚度与位移的关系进行折算,然后叠加即可得到正确结果。可显著减少图乘次数,降低计算量,提高计算效率和正确率,值得推广。     

General analytic solution of dynamic response of beams with nonhomogeneity and variable cross-section

In this paper, a new method, the step-reduction method, is proposed to investigate the dynamic response of the Bernoulli-Euler beams with arbitrary nonhomogeneity and arbitrary variable cross-section under arbitrary loads. Both free vibration and forced vibration of such beams are studied. The new method requires to discretize the space domain into a number of elements. Each element can be treated as a homogeneous one with uniform thickness. Therefore, the general analytical solution of homogeneous beams with uniform cross-section can be used in each element. Then, the general analytic solution of the whole beam in terms of initial parameters can be obtained by satisfying the physical and geometric continuity conditions at the adjacent elements. In the case of free vibration, the frequency equation in analytic form can be obtained, and in the case of forced vibration, a final solution in analytical form can also be obtained which is involved in solving a set of simultaneous algebraic equations with only     

The influence of the shear deformations on the local stress state of pultruded composite profiles

Due to the relevance of shear deformability, the practical use of composite profiles still conflicts with the serviceability requirements related to the stiffness demand for civil applications. Moreover, when dealing with shear deformable beams, it is also well recognized that displacement-based 1-D models can lead to inaccurate stress predictions.Hence, a relatively simple beam model allowing to evaluate both strains and stresses accurately may represent a useful tool.The main motivation of the present paper is precisely to investigate these features by presenting relevant numerical results dealing with the mechanical response of pultruded composite profiles with thin-walled open cross-section made of both Glass (GFRP) and Carbon Fiber Reinforced Plastic (CFRP).Comparisons with solutions given via classical 1-D/2-D mechanical models are also provided, which highlight the accuracy of the proposed kinematics, especially with the aim of a reliable stress evaluation.     

Analytical test functions for free vibration analysis of rotating non-homogeneous Timoshenko beams

In this paper, the governing equations for free vibration of a non-homogeneous rotating Timoshenko beam, having uniform cross-section, is studied using an inverse problem approach, for both cantilever and pinned-free boundary conditions. The bending displacement and the rotation due to bending are assumed to be simple polynomials which satisfy all four boundary conditions. It is found that for certain polynomial variations of the material mass density, elastic modulus and shear modulus, along the length of the beam, the assumed polynomials serve as simple closed form solutions to the coupled second order governing differential equations with variable coefficients. It is found that there are an infinite number of analytical polynomial functions possible for material mass density, shear modulus and elastic modulus distributions, which share the same frequency and mode shape for a particular mode. The derived results are intended to serve as benchmark solutions for testing approximate or numerical methods used for the vibration analysis of rotating non-homogeneous Timoshenko beams.