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.     

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.     

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 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.     

A Displacement Solution to Transverse Shear Loading of Composite Beams by BEM

In this paper the boundary element method is employed to develop a displacement solution for the general transverse shear loading problem of composite beams of arbitrary constant cross section. The composite beam (thin or thick walled) consists of materials in contact, each of which can surround a finite number of inclusions. The materials have different elasticity and shear moduli 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 bending 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. The evaluation of the transverse shear stresses at any interior point is accomplished by direct differentiation of a warping function. The shear deformation coefficients are obtained from the solution of two boundary value problems with respect to warping functions appropriately arising from the aforementioned one using only boundary integration, while the coordinates of the shear center are obtained from these functions using again only boundary integration. Three boundary value problems are formulated with respect to corresponding warping functions and solved employing a pure BEM approach. Numerical examples are worked out to illustrate the efficiency, the accuracy and the range of applications of the developed method. The accuracy of the obtained values of the resultant transverse shear stresses compared with those obtained from an exact solution is remarkable.     

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.     

Static behavior of composite beams using various refined shear deformation theories

Static behavior of composite beams with arbitrary lay-ups using various refined shear deformation theories is presented. The developed theories, which do not require shear correction factor, account for parabolical variation of shear strains and consequently shear stresses through the depth of the beam. In addition, they have strong similarity with Euler–Bernoulli beam theory in some aspects such as governing equations, boundary conditions, and stress resultant expressions. A two-noded C¹ finite element with six degree-of-freedom per node which accounts for shear deformation effects and all coupling coming from the material anisotropy is developed to solve the problem. Numerical results are performed for symmetric and anti-symmetric cross-ply composite beams under the uniformly distributed load and concentrated load. The effects of fiber angle and lay-ups on the shear deformation parameter and extension-bending-shear-torsion response are investigated.     

Analysis of shear-deformable composite beams in postbuckling

The nonlinear response of composite beams modeled according to higher-order shear deformation theories in postbuckling is investigated. The beam ends are restrained from axial movement, and as a result the contribution of the midplane stretching is considered. The equations of motion and the boundary conditions are derived using Hamilton’s principle. The shear deformation effect on the critical buckling load and static postbuckling response is introduced using classical, first-order, and higher-order shear deformation theories. This paper presents an exact solution for the static postbuckling response of a symmetrically laminated simply supported shear-deformable composite beam. The shear effect is shown to have a significant contribution to both the buckling and postbuckling behaviors. Results of this analysis show that classical and first-order theories underestimate the amplitude of buckling while all higher-order theories, considered in this study, yield very close results for the static postbuckling response.     

Elastic properties of wood with rectangular cross section under combined static axial force and torque

It is rare for a component-member of a structure to be subjected to a simple stress state. Usually, it is subjected to a multi-axial stress state in many cases. Therefore, in order to more efficiently design a structure, it is necessary to fully understand the mechanical properties of constituent materials under such a state. In this report, the effect of combined axial force and torque loading on the elastic behavior of wood (Japanese beech and Japanese cypress) was examined. As the elastic behavior, the initial slopes of the stress-strain relationships obtained from combined loading tests are estimated. The specimen had a rectangular cross section with one of its major axes lying in the fiber (longitudinal) direction. The axial force and torque were applied in the fiber direction (along L) and about an axis lying in the L direction, respectively. Combined loading tests were performed using the proportional deformation loading method and the initial constant loading method. The results obtained were summarized as follows: (1) The effect of differences in loading methods on the relationships between shear stiffnesses and the states of combined stresses was confirmed, in particular, for Japanese cypress. (2) Differences in axial stiffness were observed between the two species under compression-shear combined stress state. While the axial stiffness of Japanese beech was not affected under the combined stress state, that of Japanese cypress tended to increase under compression-shear combined stress state. (3) The difference in shear or axial stiffness between the two planes was considered to be almost constant; however, when the axial or shear stress component of the combined stresses became dominant, the difference between the two planes tended to show a larger variation.     

基于拓扑优化的直升机旋翼桨叶剖面设计

提出了一种基于拓扑优化的直升机旋翼桨叶剖面设计方法。采用了有限元方法计算直升机旋翼桨叶剖面刚度特性, 截面考虑了剪切和翘曲变形, 并消除了翘曲位移和刚体位移之间的耦合作用。基于SIMP拓扑优化算法, 以旋翼桨叶平均柔度或者剖面刚度为设计目标, 桨叶重量为约束函数, 建立了旋翼桨叶拓扑优化模型。提出的敏度求解算法具有较高的计算精度, 采用序列线性规划算法对旋翼桨叶剖面进行优化设计。结果表明在展长较小并且承受均布升力载荷情况下, Ⅱ型截面梁的柔度最小, 而当展长增大时, 工字梁截面具有最小的柔度。此外, 旋翼桨叶外载荷等对优化结果也有较大的影响。提出的拓扑优化方法适合于概念设计阶段的直升机旋翼桨叶剖面设计。     

基于n阶GBT初始轴向载荷影响下FGM梁的振动特性

研究了初始轴向载荷影响下弹性地基功能梯度材料(FGM)梁的振动特性。基于一种拓展的n阶广义剪切变形梁理论(n-GBT),以轴向位移、剪切变形挠度与弯曲变形挠度为基本未知函数,应用Hamilton原理,建立了该系统自由振动问题力学模型的控制方程。引入边界控制参数,采用一种改进型广义微分求积(MGDQ)法获得了FGM梁的静动态响应。通过算例验证并给出了GBT阶次n的理想取值,丰富梁理论的同时,可供验证或改进其它各种剪切变形梁理论;提供的数值分析方法切实可行,拓展了GDQ法的使用范围。最后,着重讨论并分析了初始轴向载荷、边界条件、梯度指标、地基刚度、跨厚比等参数对FGM梁振动特性的影响。     

Flexural–torsional postbuckling analysis of beams of arbitrary cross section

In this article, the postbuckling analysis of axially compressed elements of arbitrary cross section is presented taking into account moderately large displacements, moderately large angles of twist and employing nonlinear relationships between bending moments and curvatures. The elements are supported by the most general boundary conditions including elastic support or restraint. Based on Galerkin’s method and approximating the displacement field of the element by polynomial expressions the governing differential equations lead to a nonlinear algebraic system. The geometric, inertia, torsion and warping constants of the arbitrary beam cross section are evaluated employing the boundary element method. The proposed formulation does not stand on the assumption of a thin-walled structure and therefore the cross section’s torsional rigidity is evaluated exactly without using the so-called Saint–Venant’s torsional constant. Both the Wagner’s coefficients and the shortening effect are taken into account, while their influence is examined and discussed. Numerical examples are worked out to illustrate the efficiency, the accuracy and the range of applications of the developed method.     

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

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

Analytical and numerical investigation of the stiffness matrix for edge‐cracked circular shafts

This paper examines two engineering methods of evaluating the stress intensity factors for cracked beams and bars subjected to a combined loading and proposes innovative formulations, as far as the circular cross section is concerned. Based on the definition of the stress intensity factors, the compliance matrix is determined as the inverse of the stiffness matrix, modelling the cracked section of a beam through a line‐spring approximation with interactive forces computed within fracture mechanics. A comparative evaluation of numerical predictions based on the proposed methods is also performed with methods available from the literature. Results for free vibration analyses of beams with transverse non‐propagating open cracks are presented and compared in order to estimate the accuracy and efficiency of the proposed methods, where a good agreement is generally found. More specifically, two different coupling effects are herein analysed for circular beams subjected to a combined bending, axial and shear loading, first, and a combined bending, shear and torsion loading, subsequently.     

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.     

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     

一种考虑剪滞效应的斜支承连续箱梁挠曲扭转分析方法

通过在支承坐标系下考虑约束条件,提出一种适用于斜支承连续箱梁挠曲扭转分析的薄壁箱梁单元。该单元具有10 个自由度,可方便地考虑斜支承连续箱梁的剪滞效应和扭翘变形。选取挠曲剪滞微分方程和约束扭转微分方程的齐次解作为单元位移函数,推导出单元刚度矩阵各元素的具体表达式。从剪滞翘曲应力的轴向平衡条件出发,建立双室箱形断面的剪滞翘曲位移函数,并给出了剪滞翘曲几何特性的一般计算公式。用所编制的电算程序SSCBA 对一个3 跨斜支承双室连续箱梁模型进行计算,计算值与实测值和ANSYS 壳单元结果均吻合较好,证实该箱梁单元是可靠的。计算表明:在跨中偏心荷载作用下,斜支承连续箱梁的剪滞翘曲变形和约束扭转翘曲变形对应力分布具有显著影响。     

On approximate solutions for the deformation of plane anisotropic beams

Plane deformation of anisotropic beams with narrow rectangular cross sections exhibits coupling of stretching, bending and transverse shearing. For anisotropic cantilever beams with a stiff end-cap under end forces and an end couple, assessments were made for approximate solutions by comparing these with numerically exact finite element (FE) solutions. Specific attention is given to point-wise or approximate satisfaction of the end-fixity conditions. As approximate methodologies, (i) the elementary polynomial form of Airy's stress function for the plane stress problem in a rectangular region, (ii) a Timoshenko-type beam theory, and (iii) the Bernoulli-Euler beam theory were selected. Among these, only the polynomial form of Airy's stress function violates the point-wise end-fixity conditions. Both the polynomial Airy stress function and the Timoshenko-type beam theory successfully model the effects of transverse shear deformation and the coupling of stretching and transverse deflection. Analytical solutions demonstrate that the normal shear coupling effect increases linearly with the thickness-to-span ratios in axial normal stress and axial displacement, while the coupling manifests quadratically in transverse displacement. The comparison of end displacements with the numerically exact FE solutions indicates that the polynomial form of Airy's stress function is no better than the Timoshenko-type beam theory. Similar conclusions were reached for the problem of uniformly loaded cantilever beams. It has been found that the accurate prediction of the deformation of thick anisotropic beams with significant normal-shear coupling requires the use of higher order theories.     

A mixed finite element method for beam and frame problems

 In this work we consider solutions for the Euler-Bernoulli and Timoshenko theories of beams in which material behavior may be elastic or inelastic. The formulation relies on the integration of the local constitutive equation over the beam cross section to develop the relations for beam resultants. For this case we include axial, bending and shear effects. This permits consideration in a direct manner of elastic and inelastic behavior with or without shear deformation. A finite element solution method is presented from a three-field variational form based on an extension of the Hu–Washizu principle to permit inelastic material behavior. The approximation for beams uses equilibrium satisfying axial force and bending moments in each element combined with discontinuous strain approximations. Shear forces are computed as derivative of bending moment and, thus, also satisfy equilibrium. For quasi-static applications no interpolation is needed for the displacement fields, these are merely expressed in terms of nodal values. The development results in a straight forward, variationally consistent formulation which shares all the properties of so-called flexibility methods. Moreover, the approach leads to a shear deformable formulation which is free of locking effects – identical to the behavior of flexibility based elements. The advantages of the approach are illustrated with a few numerical examples. Dedicated to the memory of Prof. Mike Crisfield, for his cheerfulness and cooperation as a colleague and friend over many years.     

Analysis of plates stiffened by parallel beams

In this paper a general solution for the analysis of plates stiffened by parallel beams subjected to an arbitrary loading is presented. According to the proposed model, the stiffening beams 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. The aforementioned integrated tractions result in the loading of the beams as well as the additional loading of the plate. Their distribution is established by applying continuity conditions in all directions at the interfaces. The analysis of both the plate and the beams is accomplished on their deformed shape taking into account second‐order effects. Six boundary value problems with respect to the plate transverse deflection, to the plate inplane displacement components, to the beam transverse deflections, to the beam axial deformation and to the beam non‐uniform angle of twist are formulated and solved using the analog equation method (AEM), a boundary element method (BEM)‐based method employing a boundary integral equation approach. The solution of the aforementioned plate and beam problems, which are non‐linearly coupled, is achieved using iterative numerical methods. The adopted model describes better the actual response of the plate beams system and 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 evaluated lateral deflections of the plate–beams system are found to exhibit considerable discrepancy from those of other models, which neglect inplane and axial forces and deformations. Copyright © 2006 John Wiley & Sons, Ltd.