Shear deformation effect in non-linear analysis of composite beams of variable cross section

In this paper the non-linear analysis of a composite Timoshenko beam with arbitrary variable cross section undergoing moderate large deflections under general boundary conditions is presented employing the analog equation method (AEM), a BEM-based method. The composite beam consists of materials in contact, each of which can surround a finite number of inclusions. The materials have different elasticity and shear moduli with same Poisson's ratio and are firmly bonded together. The beam is subjected in an arbitrarily concentrated or distributed variable axial loading, while the shear loading is applied at the shear center of the cross section, avoiding in this way the induction of a twisting moment. To account for shear deformations, the concept of shear deformation coefficients is used. Five boundary value problems are formulated with respect to the transverse displacements, the axial displacement and to two stress functions and solved using the AEM. Application of the boundary element technique yields a system of non-linear equations from which the transverse and axial displacements are computed by an iterative process. The evaluation of the shear deformation coefficients is accomplished from the aforementioned stress functions using only boundary integration. Numerical examples are worked out to illustrate the efficiency, the accuracy, the range of applications of the developed method and the influence of the shear deformation effect.     

Nonlinear dynamic analysis of Timoshenko beams by BEM. Part I: Theory and numerical implementation

In this two-part contribution, a boundary element method is developed for the nonlinear dynamic analysis of beams of arbitrary doubly symmetric simply or multiply connected constant cross section, undergoing moderate large displacements and small deformations under general boundary conditions, taking into account the effects of shear deformation and rotary inertia. Part I is devoted to the theoretical developments and their numerical implementation and Part II discusses analytical and numerical results obtained from both analytical or numerical research efforts from the literature and the proposed method. The beam is subjected to the combined action of arbitrarily distributed or concentrated transverse loading and bending moments in both directions as well as to axial loading. To account for shear deformations, the concept of shear deformation coefficients is used. Five boundary value problems are formulated with respect to the transverse displacements, to the axial displacement and to two stress functions and solved using the Analog Equation Method, a BEM based method. Application of the boundary element technique yields a nonlinear coupled system of equations of motion. The solution of this system is accomplished iteratively by employing the average acceleration method in combination with the modified Newton–Raphson method. The evaluation of the shear deformation coefficients is accomplished from the aforementioned stress functions using only boundary integration. The proposed model takes into account the coupling effects of bending and shear deformations along the member, as well as the shear forces along the span induced by the applied axial loading.

     

Nonlinear analysis of shear deformable beam-columns partially supported on tensionless three-parameter foundation

In this paper, a boundary element method is developed for the nonlinear analysis of shear deformable beam-columns of arbitrary doubly symmetric simply or multiply connected constant cross-section, partially supported on tensionless three-parameter foundation, undergoing moderate large deflections under general boundary conditions. The beam-column is subjected to the combined action of arbitrarily distributed or concentrated transverse loading and bending moments in both directions as well as to axial loading. To account for shear deformations, the concept of shear deformation coefficients is used. Five boundary value problems are formulated with respect to the transverse displacements, to the axial displacement and to two stress functions and solved using the Analog Equation Method, a BEM-based method. Application of the boundary element technique yields a system of nonlinear equations from which the transverse and axial displacements are computed by an iterative process. The evaluation of the shear deformation coefficients is accomplished from the aforementioned stress functions using only boundary integration. The proposed model takes into account the coupling effects of bending and shear deformations along the member as well as the shear forces along the span induced by the applied axial loading. Numerical examples are worked out to illustrate the efficiency, wherever possible, the accuracy and the range of applications of the developed method.     

Nonlinear dynamic analysis of Timoshenko beams by BEM. Part II: applications and validation

In this two-part contribution, a boundary element method is developed for the nonlinear dynamic analysis of beams of arbitrary doubly symmetric simply or multiply connected constant cross section, undergoing moderate large displacements and small deformations under general boundary conditions, taking into account the effects of shear deformation and rotary inertia. In Part I the governing equations of the aforementioned problem have been derived, leading to the formulation of five boundary value problems with respect to the transverse displacements, to the axial displacement and to two stress functions. These problems are numerically solved using the Analog Equation Method, a BEM based method. In this Part II, numerical examples are worked out to illustrate the efficiency, the accuracy and the range of applications of the developed method. Thus, the results obtained from the proposed method are presented as compared with those from both analytical and numerical research efforts from the literature. More specifically, the shear deformation effect in nonlinear free vibration analysis, the influence of geometric nonlinearities in forced vibration analysis, the shear deformation effect in nonlinear forced vibration analysis, the nonlinear dynamic analysis of Timoshenko beams subjected to arbitrary axial and transverse in both directions loading, the free vibration analysis of Timoshenko beams with very flexible boundary conditions and the stability under axial loading (Mathieu problem) are presented and discussed through examples of practical interest.

     

Dynamic analysis of 3-D beam elements including warping and shear deformation effects

In this paper the dynamic analysis of 3-D beam elements restrained at their edges by the most general linear torsional, transverse or longitudinal boundary conditions and subjected in arbitrarily distributed dynamic twisting, bending, transverse or longitudinal loading is presented. For the solution of the problem at hand, a boundary element method is developed for the construction of the 14 × 14 stiffness matrix and the corresponding nodal load vector of a member of an arbitrarily shaped simply or multiply connected cross section, taking into account both warping and shear deformation effects, which together with the respective mass and damping matrices lead to the formulation of the equation of motion. To account for shear deformations, the concept of shear deformation coefficients is used, defining these factors using a strain energy approach. Eight boundary value problems with respect to the variable along the bar angle of twist, to the primary warping function, to a fictitious function, to the beam transverse and longitudinal displacements and to two stress functions are formulated and solved employing a pure BEM approach that is only boundary discretization is used. Both free and forced transverse, longitudinal or torsional vibrations are considered, taking also into account effects of transverse, longitudinal, rotatory, torsional and warping inertia and damping resistance. Numerical examples are presented to illustrate the method and demonstrate its efficiency and accuracy. The influence of the warping effect especially in members of open form cross section is analyzed through examples demonstrating the importance of the inclusion of the warping degrees of freedom in the dynamic analysis of a space frame. Moreover, the discrepancy in the dynamic analysis of a member of a spatial structure arising from the ignorance of the shear deformation effect necessitates the inclusion of this additional effect, especially in thick walled cross section members.     

Nonlinear flexural–torsional dynamic analysis of beams of variable doubly symmetric cross section—application to wind turbine towers

In this paper, a boundary element solution is developed for the nonlinear flexural–torsional dynamic analysis of beams of arbitrary doubly symmetric variable cross section, undergoing moderate large displacements, and twisting rotations under general boundary conditions, taking into account the effect of rotary and warping inertia. The beam is subjected to the combined action of arbitrarily distributed or concentrated transverse loading in both directions and to twisting and/or axial loading. Four boundary-value problems are formulated with respect to the transverse displacements, to the axial displacement, and to the angle of twist and solved using the Analog Equation Method, a Boundary Element Method (BEM) based technique. Application of the boundary element technique yields a system of nonlinear coupled Differential–Algebraic Equations (DAE) of motion, which is solved iteratively using the Petzold–Gear Backward Differentiation Formula (BDF), a linear multistep method for differential equations coupled with algebraic equations. Numerical examples of great practical interest including wind turbine towers are worked out, while the influence of the nonlinear effects to the response of beams of variable cross section is investigated.     

A BEM solution to transverse shear loading of composite beams

In this paper a boundary element method is developed for the solution of the general transverse shear loading problem of composite beams of arbitrary constant cross-section. The composite beam consists of materials in contact, each of which can surround a finite number of inclusions. The materials have different elasticity and shear moduli with same Poisson’s ratio and are firmly bonded together. The analysis of the beam is accomplished with respect to a coordinate system that has its origin at the centroid of the cross-section, while its axes are not necessarily the principal ones. The transverse shear loading is applied at the shear centre of the cross-section, avoiding in this way the induction of a twisting moment. Two boundary value problems that take into account the effect of Poisson’s ratio are formulated with respect to stress functions and solved employing a pure BEM approach, that is only boundary discretization is used. The evaluation of the transverse shear stresses is accomplished by direct differentiation of these stress functions, while both the coordinates of the shear center and the shear deformation coefficients are obtained from these functions using only boundary integration. Numerical examples with great practical interest are worked out to illustrate the efficiency, the accuracy and the range of applications of the developed method. The accuracy of the proposed shear deformation coefficients compared with those obtained from a 3-D FEM solution of the ‘exact’ elastic beam theory is remarkable.     

Nonlinear dynamic analysis of circular plates with varying thickness

The BEM is developed for nonlinear free and forced vibrations of circular plates with variable thickness undergoing large deflections. General boundary conditions are considered, which may be also nonlinear. The problem is formulated in terms of displacements. The solution is based on the concept of the analog equation, according to which the two coupled nonlinear differential equations with variable coefficients pertaining to the in-plane radial and transverse deformation are converted to two uncoupled linear ones of a substitute beam with unit axial and unit bending stiffness, respectively, under fictitious quasi-static load distributions. Numerical examples are presented which illustrate the method and demonstrate its accuracy.     

Post-critical behavior of damped beam columns with variable cross section subjected to distributed follower forces

In this paper the post-critical behavior of beam columns with variable mass and stiffness properties subjected to follower forces arbitrarily distributed along their length in the presence of damping (both internal and external) is investigated using a complete nonlinear dynamic analysis. Although the static nonlinear analysis is more economical in computational cost, it is associated only with the loss of local stability via flutter or divergence. Thus, the nonlinear dynamic analysis is adopted in order to examine the global stability of the system. The governing equations of hyperbolic type are derived in terms of the displacements by considering (a) nonlinear response including the axial deformation, (b) nonlinear response excluding the axial deformation and (c) linear response. Moreover, as the cross-sectional properties of the beam vary along its axis, the resulting coupled nonlinear differential equations have variable coefficients. Their solution is achieved using the analog equation method (AEM) of Katsikadelis. Besides its accuracy and effectiveness, this method overcomes the shortcoming of a possible FEM solution which may experience a lack of convergence. The problems treated in this investigation include beam columns with various load distributions, such as constant, linear and parabolic. Some of the conclusions detected in studying the nonlinear dynamic stability of Beck’s column with variable cross section (Katsikadelis and Tsiatas, Nonlinear dynamic stability of damped Beck’s column with variable cross section. Int. J. Non-linear Mech. 42, 164–171, 2007), are also valid for the case of distributed loads. The important, however, finding is that the post-critical response under distributed loads depends on the law of distribution of mass and stiffness properties, which may lead also to explosive flutter (unbounded amplitude), in contrast to Beck’s column (end-tip load) where the motion is always bounded.     

Shear deformation effect in plates stiffened by parallel beams

In this paper a general solution for the analysis of shear deformable stiffened plates subjected to arbitrary loading is presented. According to the proposed model, the arbitrarily placed parallel stiffening beams of arbitrary doubly symmetric cross section are isolated from the plate by sections in the lower outer surface of the plate, taking into account the arising tractions in all directions at the fictitious interfaces. These tractions are integrated with respect to each half of the interface width resulting two interface lines, along which the loading of the beams as well as the additional loading of the plate is defined. Their unknown distribution is established by applying continuity conditions in all directions at the interfaces. The utilization of two interface lines for each beam enables the nonuniform distribution of the interface transverse shear forces and the nonuniform torsional response of the beams to be taken into account. The analysis of both the plate and the beams is accomplished on their deformed shape taking into account second-order effects. The analysis of the plate is based on Reissner’s theory, which may be considered as the standard thick plate theory with which all others are compared, while the analysis of the beams is performed employing the linearized second order theory taking into account shear deformation effect. Six boundary value problems are formulated and solved using the analog equation method (AEM), a BEM based method. The solution of the aforementioned plate and beam problems, which are nonlinearly coupled, is achieved using iterative numerical methods. The adopted model permits the evaluation of the shear forces at the interfaces in both directions, the knowledge of which is very important in the design of prefabricated ribbed plates. The effectiveness, the range of applications of the proposed method and the influence of shear deformation effect are illustrated by working out numerical examples with great practical interest.     

弹性曲梁在机械和热载荷共同作用下的几何非线性模型及其数值解

基于Euler-Bernoulli梁的几何非线性理论,建立了弹性曲梁在任意分布机械载荷和热载荷共同作用下的几何非线性静平衡控制方程。该模型不仅计及了轴线伸长,同时也精确地考虑了梁的初始曲率对变形的影响以及轴向变形与弯曲变形之间的相互耦合效应。应用打靶法数值求解了半圆形曲梁在横向均匀升温作用下的非线性弯曲问题,数值比较了轴向伸长对曲梁变形的影响。     

Non-linear flexural–torsional dynamic analysis of beams of arbitrary cross section by BEM

In this paper, a boundary element method is developed for the non-linear flexural–torsional dynamic analysis of beams of arbitrary, simply or multiply connected, constant cross section, undergoing moderately large deflections and twisting rotations under general boundary conditions, taking into account the effects of rotary and warping inertia. The beam is subjected to the combined action of arbitrarily distributed or concentrated transverse loading in both directions as well as to twisting and/or axial loading. Four boundary value problems are formulated with respect to the transverse displacements, to the axial displacement and to the angle of twist and solved using the Analog Equation Method, a BEM based method. Application of the boundary element technique leads to a system of non-linear coupled Differential–Algebraic Equations (DAE) of motion, which is solved iteratively using the Petzold–Gear Backward Differentiation Formula (BDF), a linear multistep method for differential equations coupled to algebraic equations. The geometric, inertia, torsion and warping constants are evaluated employing the Boundary Element Method. The proposed model takes into account, both the Wagner's coefficients and the shortening effect. Numerical examples are worked out to illustrate the efficiency, wherever possible the accuracy, the range of applications of the developed method as well as the influence of the non-linear effects to the response of the beam.     

用材料力学公式计算任意截面短梁正应力时的修正项研究

对材料力学中梁的弯曲应力公式增加一修正项,以反映短梁弯剪翘曲变形对应力分布的影响。提出一种根据短梁横截面边界形状及艾瑞应力函数求解应力修正项的方法,应用弹性力学空间问题的一般理论,通过应力平衡方程、应变相容方程及应力边界条件,建立了关于任意截面短梁的应力修正项及剪应力的基本方程。在所建立的基本方程基础上,导出了矩形截面和圆形截面短梁修正应力的具体计算公式,该修正应力与均布荷载大小及弹性模量与剪切模量之比均成正比,但与截面惯性矩成反比。数值算例表明,本文方法计算的应力与通用有限元软件ANSYS计算的结果吻合良好,从而验证了本文方法及其基本公式的正确性。     

压电复合材料层合梁的分岔、混沌动力学与控制

研究了简支压电复合材料层合梁在轴向、横向载荷共同作用下的非线性动力学、分岔和混沌动力学响应. 基于vonKarman理论和Reddy高阶剪切变形理论,推导出了压电复合层合梁的动力学方程. 利用Galerkin法离散偏微分方程,得到两个自由度非线性控制方程,并且利用多尺度法得到了平均方程. 基于平均方程,研究了压电层合梁系统的动态分岔,分析了系统各种参数对倍周期分岔的影响及变化规律. 结果表明,压电复合材料层合梁周期运动的稳定性和混沌运动对外激励的变化非常敏感,通过控制压电激励,可以控制压电复合材料层合梁的振动,保持系统的稳定性,即控制系统产生倍周期分岔解,从而阻止系统通过倍周期分岔进入混沌运动,并给出了控制分岔图.      

考虑自平衡条件T形曲梁的静力学特性分析

为了准确分析T 形曲梁的静力学特性,该文考虑了剪滞翘曲应力自平衡条件、剪力滞后和剪切变形等因素的影响. 同时为了更好地反映T 形曲梁翼板的位移变化,4 个广义位移函数被引入,分析中以能量变分原理为基础建立了T 形曲梁静力学特性的控制微分方程和自然边界条件. 算例中,分析了剪滞翘曲应力自平衡条件、不同载荷形式和曲梁半径R 等因素对T 形曲梁静力学特性的影响,该文解析解与有限元数值解吻合更好,说明了该文方法的有效性.     

Shear correction factors and an energy-consistent beam theory

Presented here is a new derivation of shear correction factors for isotropic beams by matching the exact shear stress resultants and shear strain energy with those of the equivalent first-order shear deformation theory. Moreover, a new method of deriving in-plane and shear warping functions from available elasticity solutions is shown. The derived exact warping functions can be used to check the accuracy of a two-dimensional sectional finite-element analysis of central solutions. The physical meaning of a shear correction factor is shown to be the ratio of the geometric average to the energy average of the transverse shear strain on a cross section. Examples are shown for circular and rectangular cross sections, and the obtained shear correction factors are compared with those of Cowper (1966) . The energy-averaged shear representative is also used to derive Timoshenkos beam theory.     

计及剪切影响的RC梁的非线性分析

为提高钢筋混凝土RC梁的计算效率和精度,提出了一种基于梁截面弯矩-曲率关系的宏观有限元方法,可用于各种跨高比RC梁的材料非线性分析。首先假定了混凝土和钢筋的非线性应力-应变关系,然后引入经过修正的Rodriguez截面模型,根据边界顶点把截面划分成若干梯形单元,利用quasi-Newton法求解由两个变量耦合而成的截面非线性平衡方程,由此建立RC截面的弯矩-曲率关系。在此基础上利用Timoshenko梁弯曲理论建立考虑横向剪切变形影响的RC梁的有限元分析模型。通过对试验梁的分析对比验证了所提出的分析方法的适用性。     

A nonlinear planar beam formulation with stretch and shear deformations under end forces and moments

A new nonlinear planar beam formulation with stretch and shear deformations is developed in this work to study equilibria of a beam under arbitrary end forces and moments. The slope angle and stretch strain of the centroid line, and shear strain of cross-sections, are chosen as dependent variables in this formulation, and end forces and moments can be either prescribed or resultant forces and moments due to constraints. Static equations of equilibria are derived from the principle of virtual work, which consist of one second-order ordinary differential equation and two algebraic equations. These equations are discretized using the finite difference method, and equilibria of the beam can be accurately calculated. For practical, geometrically nonlinear beam problems, stretch and shear strains are usually small, and a good approximate solution of the equations can be derived from the solution of the corresponding Euler–Bernoulli beam problem. The bending deformation of the beam is the only important one in a slender beam, and stretch and shear strains can be derived from it, which give a theoretical validation of the accuracy and applicability of the nonlinear Euler–Bernoulli beam formulation. Relations between end forces and moments and relative displacements of two ends of the beam can be easily calculated. This formulation is powerful in the study of buckling of beams with various boundary conditions under compression, and can be used to calculate post-buckling equilibria of beams. Higher-order buckling modes of a long slender beam that have complex configurations are also studied using this formulation.     

THERMAL POST-BUCKLING OF FUNCTIONALLY GRADED MATERIAL TIMOSHENKO BEAMS

Analysis of thermal post-buckling of FGM (Functionally Graded Material) Timoshenko beams subjected to transversely non-uniform temperature rise is presented. By accurately considering the axial extension and transverse shear deformation in the sense of theory of Timoshenko beam, geometrical nonlinear governing equations including seven basic unknown functions for functionally graded beams subjected to mechanical and thermal loads were formulated. In the analysis, it was assumed that the material properties of the beam vary continuously as a power function of the thickness coordinate. By using a shooting method, the obtained nonlinear boundary value problem was numerically solved and thermal buckling and post-buckling response of transversely non-uniformly heated FGM Timoshenko beams with fixed-fixed edges were obtained. Characteristic curves of the buckling deformation of the beam varying with thermal load and the power law index are plotted. The effects of material gradient property on the buckling deformation and critical temperature of beam were discussed in details. The results show that there exists the tension-bend coupling deformation in the uniformly heated beam because of the transversely non-uniform characteristic of materials.     

A new efficient method to evaluate exact stiffness and mass matrices of non-uniform beams resting on an elastic foundation

In this paper, a new efficient method to evaluate the exact stiffness and mass matrices of a non-uniform Bernoulli–Euler beam resting on an elastic Winkler foundation is presented. The non-uniformity may result from variable cross-section and/or from inhomogeneous linearly elastic material. It is assumed that there is no abrupt variation in the cross-section of the beam so that the Euler–Bernoulli theory is valid. The method is based on the integration of the exact shape functions which are derived from the solution of the axial deformation problem of a non-uniform bar and the bending problem of a non-uniform beam which are both formulated in terms of the two displacement components. The governing differential equations are uncoupled with variable coefficients and are solved within the framework of the analog equation concept. According to this, the two differential equations with variable coefficients are replaced by two linear ones pertaining to the axial and transverse deformation of a substitute beam with unit axial and bending stiffness, respectively, under ideal load distributions. The key point of the method is the evaluation of the two ideal loads which in this work is achieved by approximating them by two polynomials. More specifically, the axial ideal load is approximated by a linear polynomial while the transverse one by a cubic polynomial. The numerical implementation of the method is simple, and the results are compared favorably to those obtained by exact solutions available in literature.