A BEM formulation based on Reissner's hypothesis for analysing the coupled stretching–bending problem of building floor structures

In this work, a plate bending formulation of the boundary element method (BEM) based on the Reissner's hypothesis to perform linear analysis of plates reinforced by rectangular beams is extended to consider the beams not displayed over their middle surface. Therefore eccentricity effects are taken into account. The building floor structure is modelled as a stiffened plate which is treated as a single body without dividing it into beam and plate elements. Moreover the equilibrium and compatibility conditions are automatically imposed by the integral equations. In the proposed model the final system of equation is obtained by coupling the bending problem to the stretching problem. Besides, in order to reduce the number of degrees of freedom, both the displacements and tractions are approximated along the beam width, leading to a model where the values are defined on the beams axis. In order to validate the proposed formulation, the numerical results are compared to a well know finite element code.     

一种分析AGS结构的三角形加筋板壳单元

基于精细三角形Mindlin板单元构造了21个自由度三角形复合材料加筋板壳单元。在该单元构造过程中,考虑了肋骨弯曲、扭转、面内剪切和横向剪切变形的影响;由于肋骨和蒙皮的位移插值函数采用了相同的形函数,保证了两者变形的协调性,同时又放松了肋骨转动的约束,故与传统的板单元相比,能较好地反映了蒙皮和肋骨的变形特征。在此单元中,肋骨放置的数量、位置和角度可以任意,为结构的单元网格剖分带来了很大的便利。算例结果验证了单元的有效性,特别是在分析高肋结构时,显示出其比传统加筋板单元具有更高精度的优点。还以一典型的先进复合材料等格栅加筋板(AGS)为例,讨论了该结构弯曲变形的力学特性。     

Boundary element analysis of reinforced shear deformable shells

In this paper, stiffened shear‐deformable shells are analysed using the boundary element method. Coupled boundary integral equations are presented for describing curved shells under general loading conditions. The equations are based on boundary integral equations for plane stress and plate bending, with coupling terms arising from the curvature of the shell. Domain integrals are transformed into boundary integrals using the dual reciprocity technique. Stiffeners are modelled as curved beams, continuously attached to the shell. Numerical solutions calculated using the present method are compared with finite element results in two examples. Copyright © 2002 John Wiley & Sons, Ltd.     

Stiffened cracked plates analysis by dual boundary element method

In this paper a boundary element formulation for analysis of shear deformable stiffened cracked plates is presented. By coupling boundary element formulation of shear deformable plate and two dimensional plane stress elasticity, dual boundary integral equations are presented. The interaction forces between stiffeners and the plate are treated as line distributed body forces along the attachment. Both concentric and eccentric stiffeners have been considered. Rectangular stiffened plate containing a single crack and double cracks subjected to uniform distributed moment on the crack surface and uniform shear load on the plate are analysed by the proposed method. Good agreement has been achieved compared with analytical solutions.     

复合材料加筋板的动力分析

本文构造了九自由度三角形拟协调罚函数复合材料板单元与六自由度复合材料梁单元,考虑了剪切变形与转动惯量的影响。用这两种单元对复合材料加筋板的自由振动、阻尼特性和瞬态响应问题进行了研究,给出一些有益结果。      

复合材料大变形任意加筋板单元

本文构造了用于复合材料偏心加筋板状结构大变形分析的任意加筋板单元。在此模型中,肋骨可放置在板单元内任意位置和任意方向,因而在网格划分时,可不必顾及肋骨的具体位置。在板和肋骨的基本方程中,引用Von-Karman大变形理论计及大变形的影响,按照 Mindlin-Reissoner一阶剪切变形理论考虑横向剪切变形。任意加筋板单元具有和常规板单元相同的自由度。数值计算表明本文的有限元模型具有很高的精度。     

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.     

联肢加劲钢板剪力墙滞回性能试验研究与数值分析

完成了3个1/3比例的3层联肢钢板剪力墙试件的低周反复加载试验。3个试件的钢板剪力墙分别采用非加劲、槽钢竖向加劲和井字加劲的形式,钢板剪力墙的竖向边缘构件采用方钢管混凝土。得到了联肢钢板剪力墙试件的荷载-位移滞回曲线和破坏形态,对试件的骨架曲线、应力发展、延性及耗能能力等进行了分析。采用有限元软件ABAQUS对试件进行了数值模拟。结果表明:非加劲和槽钢竖向加劲墙板先屈曲后屈服,井字加劲墙板先屈服后屈曲,墙板屈服后连梁与钢板剪力墙边框梁相继屈服。方钢管混凝土柱脚屈服较早,屈服后仍具有良好的承载力和弹塑性变形能力。采用非加劲墙板的试件承载力最低,滞回环捏缩效应最严重,其次是采用槽钢竖向加劲墙板的试件。采用井字加劲墙板的试件滞回环较饱满。井字加劲和槽钢竖向加劲试件的峰值荷载分别比非加劲试件的峰值荷载提高了11.7%和6.9%,井字加劲和槽钢竖向加劲试件的等效黏滞阻尼系数分别比非加劲试件的等效黏滞阻尼系数提高了65.9%和19.9%。各试件的延性系数均大于4.5,表明不同加劲形式的联肢钢板剪力墙均具有良好的延性。数值分析与试验结果吻合较好,可充分地反映试件的滞回性能和破坏过程。加劲肋对连梁和边缘构件的内力影响较小,但可显著提高剪力墙板的抗剪承载力。相较于两片单肢钢板剪力墙,联肢钢板剪力墙的承载力和耗能能力均有大于20%的提高。     

Building floor analysis by the boundary element method

In this work, a boundary element formulation to analyse slabs reinforced by beams, combined or not to define a grid sub-system, is proposed. Kirchhoff s hypothesis is assumed for the plate elements. The beams elements are not required to be displayed over the plate surface, therefore eccentricity effects are taken into account. The formulation is derived by assuming a zoned body where beam elements are introduced by degenerating plate sub-regions. After finding properly a single reciprocity for the whole body, the required integral representations are derived. The integral representations derived for this complex structural element take into account the bending and stretching effects of booth structural elements working together. The standard equilibrium and compatibility conditions along interface are naturally imposed. Moreover, the amount of degrees of freedom required along the interfaces is substantially reduced, leading therefore to small and more accurate algebraic system of linear equations. Several examples are then shown to illustrate the accuracy of the formulation, comparing the obtained results with analytical and other numerical solutions.The authors wish to thank FAPESP (São Paulo State Foundation for Scientific Research) for the financial support and Prof. Selma Hissae Shimura da Nobrega by the collaboration to run the ADINA code.     

B-Spline Wavelet on Interval Finite Element Method for Static and Vibration Analysis of Stiffened Flexible Thin Plate

A new wavelet finite element method (WFEM) is constructed in this paper and two elements for bending and free vibration problems of a stiffened plate are analyzed. By means of generalized potential energy function and virtual work principle, the formulations of the bending and free vibration problems of the stiffened plate are derived separately. Then, the scaling functions of the B-spline wavelet on the interval (BSWI) are introduced to discrete the solving field variables instead of conventional polynomial interpolation. Finally, the corresponding two problems can be resolved following the traditional finite element frame. There are some advantages of the constructed elements in structural analysis. Due to the excellent features of the wavelet, such as multi-scale and localization characteristics, and the excellent numerical approximation property of the BSWI, the precise and efficient analysis can be achieved. Besides, transformation matrix is used to translate the meaningless wavelet coefficients into physical space, thus the resolving process is simplified. In order to verify the superiority of the constructed method in stiffened plate analysis, several numerical examples are given in the end.     

A BEM formulation based on Reissner’s theory to perform simple bending analysis of plates reinforced by rectangular beams

In this work, the plate bending formulation of the boundary element method (BEM), based on the Reissner’s hypothesis, is extended to the analysis of plates reinforced by rectangular beams. This composed structure is modelled by a zoned plate, being the beams represented by narrow sub-regions with larger thickness. The integral equations are derived by applying the weighted residual method to each sub-region, and summing them to get the equation for the whole plate. Equilibrium and compatibility conditions are automatically imposed by the integral equations, which treat this composed structure as a single body. In order to decrease the number of degrees of freedom, some approximations are considered for both displacements and tractions along the beam width. The accuracy of the proposed model is illustrated by simple examples whose exact solution are known as well as by more complex examples whose numerical results are compared with a well-known finite element code.     

Finite element solution to an elastica problem of beams

In this study, finite element solution procedures are developed for an elastica problem of inextensible beams. The element stiffness matrices are obtained by using Galerkin's method. Results of a numerical example compare reasonably well with those obtained by using the elliptical integral. Extensions of this research are currently being made to include membrane effects of beams and to develop plate element stiffness matrices for elastica problems of plate structures.     

A new finite element for buckling analysis of laminated stiffened plates

A new stiffened plate element for stability analysis of laminated stiffened plates has been presented. The basic plate element is a combination of Allman's plane stress triangular element and a Discrete Kirchhoff–Mindlin plate bending element. The element includes transverse shear effects. The model accommodates any number of arbitrarily oriented stiffeners within the plate element and eliminates constraints on the mesh division of the plate. The element has no problem associated with shear locking – a phenomenon usually encountered in isoparametric elements. The stability analysis of laminated stiffened plates has been carried out under different loading conditions with the present element.     

Simple and efficient shear flexible two-node arch/beam and four-node cylindrical shell/plate finite elements

Several simple and accurate C° two-node arch/beam and four-node cylindrical shell/plate finite elements are presented in this paper. The formulation used here is based on the refined theory of thick cylindrical shells and the quasi-conforming element technique. Unlike most C° elements, the element stiffness matrix presented here is given explicitly. In spite of their simplicity, these C° finite elements posseses linear bending strains and are free from the deficiencies existing in curved C° elements such as shear and membrane locking, spurious kinematic modes and numerical ill-conditioning. These finite elements are valid not only for thick/thin beams and plates, but also for arches/straight beams and cylindrical shells/plates. Furthermore, these C° elements can automatically reduce to the corresponding C¹ beam and plate elements and give the C° beam element obtained by the reduced integration as a special case. Several numerical examples indicate that the simple two-node arch/beam and four-node cylindrical shell/plate elements given in this paper are superior to the existing C° elements with the same element degrees of freedom. Only the formulation of the rectangular cylindrical shell and plate element is presented in this paper. The formulation of an arbitrarily quadrilateral plate element will be presented in a follow-up paper³².     

用耦合有限元／边界元方法研究加筋板的声传输

采用耦合的有限元／边界元方法,建立了考虑流体结构耦合的加筋板结构声传输计算模型,并利用该模型计算研究了加筋板的声传输特性。数值计算表明：结构的传声损失与结构的固有频率和激励频率密切相关,避免共振仍是产要的;改变板厚或肋骨惯性矩后,由于结构固有频率发生了变化会使得结构传声损失曲线上极小值的位置发生改变;当声波从空气经结构向空气中传递时,增大板厚或肋骨惯性矩可明显增大结构的传声损失;而当声波从空气经结构向水中传递时,增大板厚或肋骨惯性矩也增大了结构的传声损失,但效果明显不如当声波从空气经结构向空气中传递时显著。     

Buckling analysis of corrugated plates using a mesh-free Galerkin method based on the first-order shear deformation theory

This paper deals with elastic buckling analysis of stiffened and un-stiffened corrugated plates via a mesh-free Galerkin method based on the first-order shear deformation theory (FSDT). The corrugated plates are approximated by orthotropic plates of uniform thickness that have different elastic properties along the two perpendicular directions of the plates. The key to the approximation is that the equivalaent elastic properties of the orthotropic plates are derived by applying constant curvature conditions to the corrugated sheet. The stiffened corrugated plates are analyzed as stiffened orthotropic plates. The stiffeners are modelled as beams. The stiffness matrix of the stiffened corrugated plate is obtained by superimposing the strain energy of the equivalent orthotropic plate and the beams after implementing the displacement compatibility conditions between the plate and the beams. The mesh free characteristic of the proposed method guarantee that the stiffeners can be placed anywhere on the plate, and that remeshing is avoided when the stiffener positions change. A few selected examples are studied to demonstrate the accuracy and convergence of the proposed method. The results obtained for these examples, when possible, are compared with the ANSYS solutions or other available solutions in literature. Good agreement is evident for all cases. Some new results for both trapezoidally and sinusoidally corrugated plates are then reported.     

Buckling Analysis of Plates Stiffened by Parallel Beams

In this paper a general solution for the elastic buckling analysis of plates stiffened by arbitrarily placed parallel beams of arbitrary doubly symmetric cross section subjected to an arbitrary inplane 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. 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. The unknown distribution of the aforementioned integrated tractions is established by applying continuity conditions in all directions at the two interface lines, while the analysis of both the plate and the beams is accomplished on their deformed shape. The method of analysis is based on the capability to establish the elastic and the corresponding geometric stiffness matrices of the stiffened plate with respect to a set of nodal points. Thus, the original eigenvalue problem for the differential equation of buckling is converted into a typical linear eigenvalue problem, from which the buckling loads are established numerically. For the calculation of the elastic and geometric stiffness matrices six boundary value problems are formulated and solved using the Analog Equation Method (AEM), a BEM-based method. Numerical examples with practical interest are presented. The accuracy of the results of the proposed model compared with those obtained from a 3-D FEM solution is remarkable.     

A BEM formulation for analysing the coupled stretching–bending problem of plates reinforced by rectangular beams with columns defined in the domain

In this work, a boundary element formulation to analyse plates reinforced by rectangular beams, with columns defined in the domain is proposed. The model is based on Kirchhoff hypothesis and the beams are not required to be displayed over the plate surface, therefore eccentricity effects are taken into account. The presented boundary element method formulation is derived by applying the reciprocity theorem to zoned plates, where beams are treated as thin sub-regions with larger rigidities. The integral representations derived for this complex structural element consider the bending and stretching effects of both structural elements working together. The standard equilibrium and compatibility conditions along interface are naturally imposed, being the bending tractions eliminated along interfaces. The in-plane tractions and the bending and in-plane displacements are approximated along the beam width, reducing the number of degrees of freedom. The columns are introduced into the formulation by considering domain points where tractions can be prescribed. Some examples are then shown to illustrate the accuracy of the formulation, comparing the obtained results with other numerical solutions.     

带有加强筋的Mindlin板动态刚度阵法

以加筋中厚矩形板为研究对象，推导了加筋板的动态刚度阵，为动态刚度阵法提供一种新单元。板的运动微分方程由Mindlin厚板理论给出，同时还考虑了板平面内的振动。对于板上加强筋的处理，则通过Hamilton原理对板的运动方程作相应的修正，最终得到加筋板的运动微分方程。而方程的解析解直接用于单元刚度阵的推导，所得加筋板单元的动态刚度阵结合传统有限元方法的单元组装和求解方法即可用于计算整个结构的动力响应。此外，还给出了加筋板单元的均方响应计算公式，可用来计算结构的平均振动能量。最后通过数值算例验证本文方法，计算结果与传统有限元方法进行分析比较。     

A super element model for non-linear analysis of stiffened box structures

A numerical model for non-linear static and dynamic analysis of stiffened box structures is presented. The model is based on a new super element formulation which provides complete C¹ continuity for both plate and beam elements. Geometric and material non-linearities are included and the temporal equations are solved by the implicit Newmark-β method with Newton-Raphson subiteration. The new formulation has been applied to the static, vibration and transient analysis of various structures such as flat plates, folded plates and rectangular boxes. Both isotropic and beam stiffened structures are considered and the results obtained are compared with other available solutions. It is observed that the new super element formulation can provide reasonable solutions to both linear and non-linear problems of stiffened box structures. The mathematical formulation of the model is presented in this paper, while the numerical verifications are given in the companion paper.¹