Post buckling analysis of shear deformable shallow shells by the boundary element method

This paper presents four boundary element formulations for post buckling analysis of shear deformable shallow shells. The main differences between the formulations rely on the way non‐linear terms are treated and on the number of degrees of freedom in the domain. Boundary integral equations are obtained by coupling boundary element formulation of shear deformable plate and two‐dimensional plane stress elasticity. Four different sets of non‐linear integral equations are presented. Some domain integrals are treated directly with domain discretization whereas others are dealt indirectly with the dual reciprocity method. Each set of non‐linear boundary integral equations are solved using an incremental approach, where loads and prescribed boundary conditions are applied in small but finite increments. The resulting systems of equations are solved using a purely incremental technique and the Newton–Raphson technique with the Arc length method. Finally, the effect of imperfections (obtained from a linear buckling analysis) on the post‐buckling behaviour of axially compressed shallow shells is investigated. Results of several benchmark examples are compared with the published work and good agreement is obtained. Copyright © 2010 John Wiley & Sons, Ltd.     

Buckling analysis of shear deformable plates by boundary element method

In this paper, the derivation and numerical implementation of boundary integral equations for the buckling analysis of shear deformable plates are presented. Plate buckling equations are derived as a standard eigenvalue problem. The formulation is formed by coupling boundary element formulations of shear deformable plate and two dimensional plane stress elasticity. The eigenvalue problem of plate buckling yields the critical load factor and buckling modes. The domain integrals which appear in this formulation are treated in two different ways: initially the integrals are evaluated using constant cells, and next, they are transformed into equivalent boundary integrals using the dual reciprocity method (DRM). Several examples with different geometry, loading and boundary conditions are presented to demonstrate the accuracy of the formulation. Copyright © 2004 John Wiley & Sons, Ltd.     

弹性地基板广义边值问题的边界元法

本文利用Ｈａｎｋｅｌ变换导出了弹性地基板弯曲问题的基本解，该基本解对于Ｗｉｎｋｌｅｒ地基、Ｐａｓｔｅｒｎａｋ地基和弹性半空间地基模型具有统一的表达形式。在此基础上，建立了适用于弹性地基板广义边值问题的边界积分方程组，最后文中给出了若干数值算例。     

矩形薄板面内非线性分布载荷下的辛弹性力学解

用辛弹性力学理论,根据平面矩形域本征向量展开解法,得到了对应于零本征值和非零本征值的本征向量解,以及含待定常数的面内应力分布通解,依据必须满足的应力边界条件,利用符号运算软件Maple,导出了矩形薄板在半余弦分布载荷作用下的面内应力表达式。为了验证方法的有效性和所得到的公式的正确性,具体分析了正方形薄板在两种非线性形式载荷——半余弦和抛物线分布载荷作用下的例子。算例结果与微分求积法及其有限元法得到的数值结果极其相近。基于所给出的结果,可望为工程应用中的屈曲分析提供合理的前期准备。     

A robust methodology for the simulation of postbuckling response of composite plates

The application of a robust numerical method, namely the differential quadrature method (DQM) for the analysis of buckling and postbuckling of laminated composite plates is introduced. The method is combined with an arc-length strategy to solve the resulting system of nonlinear equations. The treatment accounts for the effect of large deformation by including the von Karman strains. Imperfections in the form of global deflection, conforming to the preferred buckling modes are introduced throughout the plate. Case studies are used to evaluate the geometrically nonlinear response of composite plates under the set conditions, and the results are compared with those obtained by finite element solution and the results obtained from a published literature.     

Buckling analysis of shear deformable shallow shells by the boundary element method

In this work a boundary element (BE) formulation for buckling problem of shear deformable shallow shells is presented. A set of five boundary integral equations are obtained by coupling two-dimensional plane stress elasticity with shear deformable plate bending (Reissner). The domain integrals appearing in the formulation (due to the curvature and due to the domain load) are transferred into equivalent boundary integrals. The BE formulation is presented as an eigenvalue problem, to provide direct evaluation of critical load factors and buckling modes. Several examples are presented. The BE results for a cylindrical shallow shell with different curvatures are compared with other numerical solutions and good agreements are obtained.     

刚性基底上弹性约束矩形板的屈曲行为分析

对刚性基底上非加载边弹性约束矩形板在固支边均匀面内压力作用下的屈曲行为进行了分析研究。用里兹能量变分法建立确定板屈曲强度的本征值问题,利用满足边界条件的屈曲位移函数导出板屈曲系数计算公式。将该方法用于矩形钢管混凝土钢板局部屈曲的分析和计算,提出了目标板边界约束系数计算公式。试验结果表明,方形钢管混凝土柱钢板局部屈曲强度计算值与试验值吻合良好,钢板边界弹性约束的处理和约束系数的计算方法是有效的。     

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.     

Buckling of the CCFF orthotropic rectangular plates under in-plane pure bending

Solution of the buckling problem for the CCFF orthotropic plate subjected to in-plane pure bending is presented. The two parallel clamped edges of the plate are loaded by linearly distributed in-plane loads statically equivalent to the in-plane bending moments. The problem is solved using method of lines for partial differential equations and Galerkin’s method. The buckling problems are solved for isotropic, orthotropic and multilayered CFRP composite plates with various aspect ratios. Results of calculations of critical loads are compared with those based on finite-element modelling and analyses. The comparisons demonstrate efficiency of the proposed approach to the buckling analysis of composite CCFF plates with various dimensional and stiffness parameters.     

Exact buckling solutions for composite laminates: proper free edge conditions under in-plane loadings

Summary This paper considers the elastic buckling of symmetric cross-ply laminated rectangular plates with two parallel edges simply supported, one edge free and the remaining edge free, simply supported or clamped. The first-order shear deformation plate theory is used in the analysis. An error apparently made by previous researchers for boundary conditions at free edges subjected to in-plane loads is corrected. Closed-form buckling factors are obtained using a generalised Levy-type solution method to solve the differential equations which govern the buckling behaviour of the laminates. Comparisons are made with previously published results, and the differences between buckling factors obtained with the appropriate and inappropriate free edge conditions are examined. The variation of buckling factors with plate aspect ratio, thickness ratio and the number of layers is investigated. Sets of first-known buckling solutions for cross-ply laminates are reported in design charts and tables.     

Multiple-parameter reduced basis technique for bifurcation and post-buckling analyses of composite plates

A multiple-parameter reduced basis technique and a problem-adaptive computational algorithm are presented for the bifurcation and post-buckling analyses of composite plates subjected to combined loadings. The computational algorithm can be conveniently divided into three distinct stages. The first stage is that of determining the stability boundary. The plate is discretized by using displacement finite element models and the analysis region is reduced by exploiting the special symmetries exhibited by the response of the plate. The vector of unknown nodal displacements is expressed as a linear combination of a small number of path derivatives (derivatives of the nodal displacements with respect to path parameters), and a Rayleigh-Ritz technique is used to approximate the finite element equations by a small system of algebraic equations. The reduced equations are used to determine the stability boundary of the plate. In the second stage, a nonllnear solution in the vicinity of the stability boundary is obtained by using a bifurcation buckling mode as a predictor, and a set of reduced equations is generated. In the third stage, the reduced equations are used to trace post-buckling paths corresponding to various combinations of the load parameters. The potential of the proposed approach is discussed and its effectiveness is demonstrated by means of a numerical example of laminated composite plate subjected to combined compressive and shear loadings.     

考虑剪切效应的轴对称脱层圆板的后屈曲分析

基于Von Karman板理论,考虑横向剪切变形,建立了具脱层的轴对称层合圆板的后屈曲控制方程。应用正交配点法,将后屈曲控制方程、边界条件、以及连续条件转化为非线性方程组,然后进行迭代求解。讨论了不同脱层深度和脱层半径对层合圆板的屈曲及后屈曲特性影响,且与有关文献的结果进行了比较。     

Thermal buckling analysis of moderately thick functionally graded annular sector plates

In this article, thermal buckling analysis of moderately thick functionally graded annular sector plate is studied. The equilibrium and stability equations are derived using first order shear deformation plate theory. These equations are five highly coupled partial differential equations. By using an analytical method, the coupled stability equations are replaced by four decoupled equations. Solving the decoupled equations and satisfying the boundary conditions, the critical buckling temperature is found analytically. To this end, it is assumed that the annular sector plate is simply supported in radial edges and it has arbitrary boundary conditions along the circular edges. Thermal buckling of functionally graded annular sector plate for two types of thermal loading, uniform temperature rise and gradient through the thickness, are investigated. Finally, the effects of boundary conditions, power law index, plate thickness, annularity and sector angle on the critical buckling temperature of functionally graded annular sector plates are discussed in details.     

Primitive variable dual reciprocity boundary element method solution of incompressible Navier–Stokes equations

This paper describes the solution to transient incompressible two-dimensional Navier–Stokes equations in primitive variables by the dual reciprocity boundary element method. The coupled set of mass and momentum equations is structured by the fundamental solution of the Laplace equation. The dual reciprocity method is based on the augmented thin plate splines. All derivatives involved are calculated through integral representation formulas. Numerical example include convergence studies with different mesh size for the classical lid-driven cavity problem at Re=100 and comparison with the results obtained through calculation of the derivatives from global interpolation formulas. The accuracy of the solution is assessed by comparison with the Ghia–Ghia–Shin finite difference solution as a reference.     

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.     

Evaluation of crack development in ribbed panels

The problem of elastic equilibrium for a finite anisotropic ribbed panel weakened by a rectilinear crack is considered. Within the framework of a hypothesis of continuous contact of a rib and a plate along a line special presentations are given for obtaining and solving a set of singular integral equations. Absence of an unknown function in the notch contour makes it possible to reduce the order of the set of singular integral equations by one and to propose an effective algorithm for numerical solution of the problem. The effect of material anisotropy, plate edges, and reinforcement stiffness on the stress intensity factor at the crack tip are studied. Calculated results are compared with experimental data for the case of an isotropic plate.Translated from Problemy Prochnosti, No. 3, pp. 63–68, March, 1991.     

A BEM-based meshless solution to buckling and vibration problems of orthotropicplates

In this paper a BEM-based meshless solution is presented to buckling and vibration problems of Kirchhoff orthotropic plates with arbitrary shape. The plate is subjected to compressive centrally applied load together with arbitrarily transverse distributed or concentrated loading, while its edges are restrained by the most general linear boundary conditions. The resulting buckling and vibration problems are described by partial differential equations in terms of the deflection. Both problems are solved employing the Analog Equation Method (AEM). According to this method the fourth-order partial differential equation describing the response of the orthotropic plate is converted to an equivalent linear problem for an isotropic plate subjected only to a fictitious load under the same boundary conditions. The AEM is applied to the problem at hand as a boundary-only method by approximating the fictitious load with a radial basis function series. Thus, the method retains all the advantages of the pure BEM using a known simple fundamental solution. Example problems are presented for orthotropic plates, subjected to compressive or vibratory loading, to illustrate the method and demonstrate its efficiency and its accuracy.     

Exact solution for stability analysis of moderately thick functionally graded sector plates on elastic foundation

In this article, an exact analytical solution for buckling analysis of moderately thick functionally graded (FG) sector plates resting on Winkler elastic foundation is presented. The equilibrium equations are derived according to the first order shear deformation plate theory. Because of the coupling between the bending and stretching equilibrium equations of FG plates, these plates have deflection under in-plane loads lower than the critical buckling load acting on the mid-plane. The conditions under which FG plates remain flat in the pre-buckling configuration are investigated and the stability equations are obtained based on the flat plate assumption in the pre-buckling state. The stability equations are simplified into decoupled equations and solved analytically for plates having simply supported boundary condition on the straight edges. The critical buckling load is obtained and the effects of geometrical parameters and power law index on the stability of functionally graded sector plates are studded. The results for the critical buckling load of moderately thick functionally graded sector plates resting on elastic foundation are reported for the first time.     

模糊随机结构屈曲问题的区间有限元法

利用可靠度理论中的置信区间和模糊集理论中的λ截集方法,可以把结构中的随机参数和模糊参数转化为区间数.这时模糊随机有限元平衡方程转化为区间方程组,因此利用区间运算方法或蒙特卡洛直接抽样法可以求解模糊随机结构屈曲问题,并且得到一个区间解.如果结合结构稳定性理论,在某些特定的情况下,其计算量与求解一个相应的确定性问题有限元法的计算量相当.