A BEM solution to transverse shear loading of beams

In this paper a boundary element method is developed for the solution of the general transverse shear loading problem of beams of arbitrary simply or multiply connected constant cross section. 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 center 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 at any interior point 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.     

Shear deformation effect in nonlinear analysis of spatial composite beams subjected to variable axial loading by bem

Summary In this paper a boundary element method is developed for the nonlinear analysis of composite beams of arbitrary doubly symmetric constant cross section, taking into account the shear deformation effect. 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 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. Numerical examples with great practical interest are worked out to illustrate the efficiency and the range of applications of the developed method. The influence of both the shear deformation effect and the variableness of the axial loading are remarkable.     

Shear deformation effect in second-order analysis of frames subjected to variable axial loading

In this paper a boundary element method is developed for the second-order analysis of frames consisting of beams of arbitrary simply or multiply connected constant cross section, taking into account shear deformation effect. Each beam is subjected to 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. Three boundary value problems are formulated with respect to the beam deflection, the axial displacement and to a stress function and solved employing a BEM approach. The evaluation of the shear deformation coefficients is accomplished from the aforementioned stress function 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 influence of both the shear deformation effect and the variableness of the axial loading are remarkable.     

Flexural - Torsional Nonlinear Analysis of Timoshenko Beam-Column of Arbitrary Cross Section by BEM

In this paper a boundary element method is developed for the nonlinear flexural - torsional analysis of Timoshenko beam-columns of arbitrary simply or multiply connected constant cross section, undergoing moderate large deflections under general boundary conditions. The beam-column is subjected to the combined action of an arbitrarily distributed or concentrated axial and transverse loading as well as to bending and twisting moments. To account for shear deformations, the concept of shear deformation coefficients is used. Seven boundary value problems are formulated with respect to the transverse displacements, to the axial displacement, to the angle of twist (which is assumed to be small), to the primary warping function 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 as well as the angle of twist 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 with great practical interest are worked out to illustrate the efficiency, the accuracy and the range of applications of the developed method. The influence of both the shear deformation effect and the variableness of the axial loading are remarkable.     

Warping shear stresses in nonuniform torsion by BEM

 In this paper a boundary element method is developed for the nonuniform torsion of simply or multiply connected prismatic bars of arbitrary cross section. The bar is subjected to an arbitrarily distributed twisting moment, while its edges are restrained by the most general linear torsional boundary conditions. Since warping is prevented, beside the Saint–Venant torsional shear stresses, the warping normal and shear stresses are also computed. Three boundary value problems with respect to the variable along the beam angle of twist and to the primary and secondary warping functions are formulated and solved employing a BEM approach. Both the warping and the torsion constants using only boundary discretization together with the torsional shear stresses and the warping normal and shear stresses are computed. Numerical results are presented to illustrate the method and demonstrate its efficiency and accuracy. The magnitude of the warping shear stresses due to restrained warping is investigated by numerical examples with great practical interest. Received: 13 November 2001 / Accepted: 2 October 2002     

Shear deformation effect in flexural–torsional vibrations of beams by BEM

In this paper, a boundary element method is developed for the general flexural–torsional vibration problem of Timoshenko beams of arbitrarily shaped cross section taking into account the effects of warping stiffness, warping and rotary inertia and shear deformation. The beam is subjected to arbitrarily transverse and/or torsional distributed or concentrated loading, while its edges are restrained by the most general linear boundary conditions. The resulting initial boundary value problem, described by three coupled partial differential equations, is solved employing a boundary integral equation approach. Besides the effectiveness and accuracy of the developed method, a significant advantage is that the displacements as well as the stress resultants are computed at any cross-section of the beam using the respective integral representations as mathematical formulae. All basic equations are formulated with respect to the principal shear axes coordinate system, which does not coincide with the principal bending one in a nonsymmetric cross section. To account for shear deformations, the concept of shear deformation coefficients is used. Six boundary value problems are formulated with respect to the transverse displacements, to the angle of twist, to the primary warping function and to two stress functions and solved using the Analog Equation Method, a BEM based method. Both free and forced vibrations are examined. Several beams are analysed to illustrate the method and demonstrate its efficiency and wherever possible its accuracy.     

Advanced 3D beam element of arbitrary composite cross section including generalized warping effects

In this paper, an advanced 20 × 20 stiffness matrix and the corresponding nodal load vector of a member of arbitrary composite cross section is developed taking into account shear lag effects due to both flexure and torsion (the case of the three‐dimensional beam element of arbitrary homogeneous cross section is treated as a special one). The composite member consists of materials in contact each of which can surround a finite number of inclusions. Nonuniform warping distributions are taken into account by employing four independent warping parameters multiplying a shear warping function in each direction and two torsional warping functions. Ten boundary value problems with respect to the kinematical components are formulated and solved employing the analog equation method, a BEM‐based technique. The aforementioned boundary value problems are formulated employing either an improved stress field arising from the correction of shear stress components or the stress field arising directly from displacement considerations. The warping functions and the geometric constants including the additional ones due to warping are evaluated employing a pure BEM approach, that is, only boundary discretization of the cross section is used. Numerical results are presented to illustrate the method and demonstrate its efficiency and accuracy. The deviations arising from the use of the advanced 20 × 20 stiffness matrix and the classical 12 × 12 or 14 × 14 ones employed in commercial software packages are illustrated through examples of great practical interest. Moreover, the influence of nonuniform warping effects necessitating the use of the aforementioned ‘advanced’ stiffness matrix is also demonstrated. Copyright © 2015 John Wiley & Sons, Ltd.     

考虑滑移、剪力滞后和剪切变形的钢-混凝土组合梁解析解

通过在Newmark 模型中引入(1)描述横向非均匀分布的纵向位移的翘曲形函数和(2)描述钢梁腹板剪切变形的Timoshenko 梁假定,建立了一个能考虑滑移、剪力滞后和剪切变形的钢-混凝土组合梁模型,并推导了均布荷载作用下的解析解。最后通过4 个算例验证了模型和解析解的正确性和适用性,并显示了考虑组合梁剪切变形的必要性。另外,算例还表明,在组合梁的三维有限元建模中采用Timoshenko 梁单元来考虑钢梁的剪切变形会导致显著的误差。     

复合材料梁的几何精确非线性建模及非经典效应

基于Hodges的广义Timoshenko梁理论对具有任意剖面形状、任意材料分布及大变形的复合材料梁进行几何精确非线性建模,采用旋转张量分解法计算梁内任意一点的应变,采用变分渐近法确定梁剖面的任意翘曲,采用平衡方程由二次渐近精确的应变能导出广义Timoshenko应变能,采用广义Hamilton原理建立梁的几何精确非线性运动方程。将所建模型用于复合材料梁的静动力分析,通过与实验数据的对比,验证了建模方法的准确性,并进一步研究了剖面翘曲及横向剪切变形非经典效应对复合材料梁的影响。研究表明,剖面翘曲对复合材料梁的静变形和固有频率有显著影响,横向剪切变形对复合材料梁的静变形和固有频率的影响与梁的长度/剖面高度比有关。     

A C1 finite element for flexural and torsional analysis of rectangular piezoelectric laminated/sandwich composite beams

This work deals with the development of a new C¹ finite element for analysing the bending and torsional behaviour of rectangular piezoelectric laminated/sandwich composite beams. The formulation includes transverse shear, warping due to torsion, and elastic–electric coupling effects. It also accounts for the inter-layer continuity condition at the interfaces between layers, and the boundary conditions at the upper and lower surfaces of the beam. The shear strain is represented by a cosine function of a higher order in nature and thus avoiding shear correction factors. The warping function obtained from a three-dimensional elasticity solution is incorporated in the present model. An exact integration is employed in evaluating various energy terms due to the application of field consistency approach while interpolating the transverse shear and torsional strains. The variation of the electric potential through the thickness is taken care of in the formulation based on the observation of three-dimensional solution. The performance of the laminated piezoelectric element is tested comparing with analytical results as well as with the reference solutions evaluated using three-dimensional finite element procedure. A detailed study is conducted to highlight the influence of length-to-thickness ratio on the displacements, stresses and electric potential field of piezoelectric laminated beam structures subjected to flexural and torsional loadings. Copyright © 2004 John Wiley & Sons, Ltd.     

A C1 finite element including transverse shear and torsion warping for rectangular sandwich beams

A new three‐noded C¹ beam finite element is derived for the analysis of sandwich beams. The formulation includes transverse shear and warping due to torsion. It also accounts for the interlaminar continuity conditions at the interfaces between the layers, and the boundary conditions at the upper and lower surfaces of the beam. The transverse shear deformation is represented by a cosine function of a higher order. This allows us to avoid using shear correction factors. A warping function obtained from a three‐dimensional elasticity solution is used in the present model. Since the field consistency approach is accounted for interpolating the transverse strain and torsional strain, an exact integration scheme is employed in evaluating the strain energy terms. Performance of the element is tested by comparing the present results with exact three‐dimensional solu‐tions available for laminates under bending, and the elasticity three‐dimensional solution deduced from the de Saint‐Venant solution including both torsion with warping and bending. In addition, three‐dimensional solid finite elements using 27 noded‐brick elements have been used to bring out a reference solution not available for sandwich structures having high shear modular ratio between skins and core. A detailed parametric study is carried out to show the effects of various parameters such as length‐to‐thickness ratio, shear modular ratio, boundary conditions, free (de Saint‐Venant) and constrained torsion. Copyright © 1999 John Wiley & Sons, Ltd.     

Nonlinear analysis of concrete shells including effects of normal and transverse shear stresses

The previously developed numerical model of the authors for the analysis of conventional reinforced and prestressed concrete shells under short‐term and long‐term loading was improved by including the effects of transverse shear stresses on the shell failure. The 9‐node degenerated shell element with the layered material model through the thickness of the shell was used. The reinforcement was modelled as a separate layer. To include the effect of transverse shear stresses on the shell failure, the failure criterion for concrete and longitudinal reinforcement was defined by a relation of transverse shear stresses and normal stresses in two mutually perpendicular vertical planes. The total transverse shear bearing capacity of the shell cross‐section is obtained by summing up the concrete and reinforcement contributions. The developed numerical model and appropriate software were verified based on experimental tests.     

Higher order refined computational model with 12 degrees of freedom for the stress analysis of antisymmetric angle-ply plates – analytical solutions

Analytical formulations and solutions for the stress analysis of simply supported antisymmetric angle-ply composite and sandwich plates hitherto not reported in the literature based on a higher order refined computational model with twelve degrees of freedom already reported in the literature are presented. The theoretical model presented herein incorporates laminate deformations which account for the effects of transverse shear deformation, transverse normal strain/stress and a nonlinear variation of in-plane displacements with respect to the thickness coordinate thus modelling the warping of transverse cross sections more accurately and eliminating the need for shear correction coefficients. In addition, two higher order computational models, one with nine and the other with five degrees of freedom already available in the literature are also considered for comparison. The equations of equilibrium are obtained using Principle of Minimum Potential Energy (PMPE). Solutions are obtained in closed form using Navier’s technique by solving the boundary value problem. Accuracy of the theoretical formulations and the solution method is first ascertained by comparing the results with that already available in the literature. After establishing the accuracy of the solutions, numerical results with real properties using all the computational models are presented for the stress analysis of multilayer antisymmetric angle-ply composite and sandwich plates, which will serve as a benchmark for future investigations.     

Flexural-Torsional Buckling and Vibration Analysis of Composite Beams

In this paper the general flexural-torsional buckling and vibration problems of composite Euler-Bernoulli beams of arbitrarily shaped cross section are solved using a boundary element method. The general character of the proposed method is verified from the formulation of all basic equations with respect to an arbitrary coordinate system, which is not restricted to the principal one. The composite beam consists of materials in contact each of which can surround a finite number of inclusions. It is subjected to a compressive centrally applied load together with arbitrarily transverse and/or torsional distributed or concentrated loading, while its edges are restrained by the most general linear boundary conditions. The resulting problems are (i) the flexural-torsional buckling problem, which is described by three coupled ordinary differential equations and (ii) the flexural-torsional vibration problem, which is described by three coupled partial differential equations. Both problems are solved employing a boundary integral equation approach. Besides the effectiveness and accuracy of the developed method, a significant advantage is that the method can treat composite beams of both thin and thick walled cross sections taking into account the warping along the thickness of the walls. The proposed method overcomes the shortcoming of possible thin tube theory (TTT) solution, which its utilization has been proven to be prohibitive even in thin walled homogeneous sections. Example problems of composite beams are analysed, subjected to compressive or vibratory loading, to illustrate the method and demonstrate its efficiency and wherever possible its accuracy. Moreover, useful conclusions are drawn from the buckling and dynamic response of the beam.     

复合材料层板层间缺陷分析——剪切滑移

根据三维弹性平衡方程和层间剪切滑移条件,导出了一个复合材料层板层间剪切滑移模型。本文模型具有一般形式的二维板壳理论的位移场及其平衡方程,但因引入了能反映层板界面粘贴情况及板面条件的剪切变形函数,模型因而简单又精确。层板的弯曲问题和屈曲问题被考虑,层间弱粘贴的影响被讨论。数值结果与精确解比较,表明了本文模型的高精度。      

A refined higher-order composite box beam theory

A three-dimensional theory is developed to model composite box beams with arbitrary wall thicknesses. The theory, which is based on a refined displacement field, approximates the three-dimensional elasticity solution so that the beam cross-sectional properties are not reduced to one-dimensional beam parameters. Both in-plane and out-of-plane warping are included automatically in the formulation. The model can accurately capture the transverse shear stresses through the thickness of each wall while satisfying stress-free boundary conditions on the inner and outer surfaces of the beam. Numerical results are presented for beams with varying wall thicknesses and aspect ratios. The static results are correlated with available experimental data and show excellent agreement. Results presented for thick-walled box beams show the importance of including transverse shear in the formulation and the difficulty of defining a ‘beam’ twist for the entire cross-section.     

自然弯扭梁的耦合振动特性分析

以自然弯扭梁理论为基础对具有一般横截面形状空间曲梁的耦合振动特性进行了研究。 在该梁的运动控制方程中,位移函数和广义翘曲坐标均被定义在形心轴上,且在分析中包括了转动惯量、横向剪切变形以及和扭转有关的翘曲对振动的影响。通过对数学计算软件MATHEMATICA的精确运用可以得到该梁振型的解析表达式,精确的固有频率则可用搜索的方法来确定。为了证明理论的有效性,对两端固支椭圆截面曲梁的固有频率和振型进行了求解,并把数值计算结果同使用PATRAN梁单元的有限元结果进行了比较。     

ㄇ形梁剪力滞效应的解析解

精确分析ㄇ形梁在纵横向荷载共同作用下,其横断面上正应力分布规律对于计算其有效宽度有重要意义,应用力法原理,先将ㄇ形梁和翼板截开成矩形截面梁和平面应力板,在截面上代之以赘余的分布剪力,对于平面应力板,通过利用板变形的对称性来简化其边界条件,然后假设一个满足板的控制方程的Airy应力函数求得板的应力和位移,再利用Timoshenko梁理论求得梁的挠度和转角,根据截面上梁与板的纵向位移相等的变形协调条件便可最终确定截面上的分布剪力,给出的数值算例验证了方法的有效性,并与铜陵长江公路大桥主梁的模型有限元结果和试验结果作了对比,解析解法还可用来检验其他各种数值计算方法的精度,并可推广到其他多跨薄壁结构梁桥的膜应力分析中。     

A generalized layerwise higher-order shear deformation theory for laminated composite and sandwich plates based on isogeometric analysis

This paper presents a generalized layerwise higher-order shear deformation theory for laminated composite and sandwich plates. We exploit a higher-order shear deformation theory in each layer such that the continuity of the displacement and transverse shear stresses at the layer interfaces is ensured. Thanks for enforcing the continuity of the displacement and transverse shear stresses at an inner-laminar layer, the minimum number of variables is retained from the present theory in comparison with other layerwise theories. The method requires only five variables, the same as what obtained from the first- and higher-order shear deformation theories. In comparison with the shear deformation theories based on the equivalent single layer, the present theory is capable of producing a higher accuracy for inner-laminar layer shear stresses. The free boundary conditions of transverse shear stresses at the top and bottom surfaces of the plate are fulfilled without any shear correction factors. The discrete system equations are derived from the Galerkin weak form, and the solution is obtained by isogeometric analysis (IGA). The discrete form requires the C¹ continuity of the transverse displacement, and hence NURBS basis functions in IGA naturally ensure this condition. The laminated composite and sandwich plates with various geometries, aspect ratios, stiffness ratios and boundary conditions are studied. The obtained results are compared with the 3D elasticity solution, the analytical as well as numerical solutions based on various plate theories.