Limit analysis of multi-layered plates. Part II: Shear effects

In the first part of this work [Dallot, J., Sab, K., 2007. Limit analysis of multi-layered plates. Part I: the homogenized Love-Kirchhoff model. J. Mech. Phys. Solids, in press, doi:10.1016/j.jmps.2007.05.005], the limit analysis of a multi-layered plastic plate submitted to out-of-plane loads was studied. The authors have shown that a homogeneous equivalent Love-Kirchhoff plate can be substituted for the heterogeneous multi-layered plate, as the slenderness (length-to-thickness) ratio goes to infinity. In fact, the out-of-plane shear stresses are shown to become asymptotically negligible when compared to in-plane stresses, as the slenderness ratio goes to infinity. Actually, failure of thick multi-layered structures often occurs by shearing in the core layers and sliding at the interfaces between the layers. Both shearing and sliding are caused by the out-of-plane shear stresses. The purpose of the present paper is to build an enhanced Multi-particular Model for Multi-layered Material (M4) taking into account shear stress effects. In this model, each layer is seen as a Reissner-Mindlin plate interacting with its neighboring layers through interfaces. The proposed model is asymptotically consistent with the homogenized Love-Kirchhoff model described in the first part of the work, as the slenderness ratio goes to infinity. Kinematic and static methods for the determination of the limit load of a thick multi-layered plate which is submitted to out-of-plane distributed forces are described. The special case of multi-layered plates under cylindrical bending conditions is studied. These conditions lead to simplifications which often allow for the analytical resolution of the Love-Kirchhoff and the M4 limit analysis problems. The benefit of the proposed M4 model is demonstrated on an example. A comparison between the heterogeneous 3D model, the Love-Kirchhoff model and the M4 model is performed on a three-layer sandwich plate under cylindrical bending conditions. Finite element calculations are used to solve the 3D problem, while both the Love-Kirchhoff and the M4 problems are analytically solved. It is shown that, when the contrast between the core and the skins strengths is high, the Love-Kirchhoff model fails to capture the plastic collapse modes that cause the ruin of the sandwich plate. These modes are well captured by the M4 model which predicts limit loads that are very consistent with the limit loads predicted by the heterogeneous 3D model (the relative error is found to be smaller than 1%).     

A computational homogenization approach for the yield design of periodic thin plates. Part II : Upper bound yield design calculation of the homogenized structure

In the first part of this work (Bleyer and de Buhan, 2014), the determination of the macroscopic strength criterion of periodic thin plates has been addressed by means of the yield design homogenization theory and its associated numerical procedures. The present paper aims at using such numerically computed homogenized strength criteria in order to evaluate limit load estimates of global plate structures. The yield line method being a common kinematic approach for the yield design of plates, which enables to obtain upper bound estimates quite efficiently, it is first shown that its extension to the case of complex strength criteria as those calculated from the homogenization method, necessitates the computation of a function depending on one single parameter. A simple analytical example on a reinforced rectangular plate illustrates the simplicity of the method. The case of numerical yield line method being also rapidly mentioned, a more refined finite element-based upper bound approach is also proposed, taking dissipation through curvature as well as angular jumps into account. In this case, an approximation procedure is proposed to treat the curvature term, based upon an algorithm approximating the original macroscopic strength criterion by a convex hull of ellipsoids. Numerical examples are presented to assess the efficiency of the different methods.     

A computational homogenization approach for the yield design of periodic thin plates. Part I: Construction of the macroscopic strength criterion

The purpose of this paper is to propose numerical methods to determine the macroscopic bending strength criterion of periodically heterogeneous thin plates in the framework of yield design (or limit analysis) theory. The macroscopic strength criterion of the heterogeneous plate is obtained by solving an auxiliary yield design problem formulated on the unit cell, that is the elementary domain reproducing the plate strength properties by periodicity. In the present work, it is assumed that the plate thickness is small compared to the unit cell characteristic length, so that the unit cell can still be considered as a thin plate itself. Yield design static and kinematic approaches for solving the auxiliary problem are, therefore, formulated with a Love–Kirchhoff plate model. Finite elements consistent with this model are proposed to solve both approaches and it is shown that the corresponding optimization problems belong to the class of second-order cone programming (SOCP), for which very efficient solvers are available. Macroscopic strength criteria are computed for different type of heterogeneous plates (reinforced, perforated plates,…) by comparing the results of the static and the kinematic approaches. Information on the unit cell failure modes can also be obtained by representing the optimal failure mechanisms. In a companion paper, the so-obtained homogenized strength criteria will be used to compute ultimate loads of global plate structures.     

On the stability of multi-layered fluids. Part I: analysis

The linear stability of the Poiseuille flow of multi-layered different fluids, described mathematically by a system of Orr-Somerfeld differential equations, is investigated. A spectral method is used to rewrite this system into a generalized eigenvalue problem, which can be solved with the QZ-algorithm. Special attention is paid to the tractibility of the interfacial conditions of the stability problem. Since we will limit ourselves to a linear stability analysis, the analytical treatment of the interfacial conditions is simplified. Some results related to simple flow configurations are presented. The origin of certain regions of interfacial instability is explained by simple analytical reasoning.     

Nonlinear analysis of skewed sandwich plates

Nonlinear static and dynamic behaviours of freely supported Rhombic sandwich plates have been studied following Banerjee's hypothesis. Numerical results for 0° skew angle are compared with other known results. Results for other skew angles are believed to be new.

---

Sommario Si studia, seguendo l'ipotesi di Banerjee, il comportamento nonlineare statico e dinamico di piastre rombiche semplicemente appoggiate. Si presentano risultati numerici relativi a piastre rombiche e rettangolari: questi ultimi vengono paragonati a risultati già noti, mentre i primi si ritengono nuovi.

     

A Bending-Gradient model for thick plates. Part I: Theory

This is the first part of a two-part paper dedicated to a new plate theory for out-of-plane loaded thick plates where the static unknowns are those of the Kirchhoff–Love theory (3 in-plane stresses and 3 bending moments), to which six components are added representing the gradient of the bending moment. The new theory, called the Bending-Gradient plate theory is described in the present paper. It is an extension to arbitrarily layered plates of the Reissner–Mindlin plate theory which appears as a special case of the Bending-Gradient plate theory when the plate is homogeneous. However, we demonstrate also that, in the general case, the Bending-Gradient model cannot be reduced to a Reissner–Mindlin model. In part two (Lebée and Sab, 2011), the Bending-Gradient theory is applied to multilayered plates and its predictions are compared to those of the Reissner–Mindlin theory and to full 3D Pagano’s exact solutions. The main conclusion of the second part is that the Bending-Gradient gives good predictions of both deflection and shear stress distributions in any material configuration. Moreover, under some symmetry conditions, the Bending-Gradient model coincides with the second-order approximation of the exact solution as the slenderness ratio L/h goes to infinity.     

Influence of skin/core debonding on free vibration behavior of foam and honeycomb cored sandwich plates

The dynamic behavior of partially delaminated at the skin/core interface sandwich plates with flexible cores is studied. The commercial finite element code ABAQUS is used to calculate natural frequencies and mode shapes of the sandwich plates containing a debonding zone. The influence of the debonding size, debonding location and types of debonding on the modal parameters of damaged sandwich plates with various boundary conditions is investigated. The results of dynamic analysis illustrated that they can be useful for analyzing practical problems related to the non-destructive damage detection of partially debonded sandwich plates.     

On the finite displacement analysis of quadrangular plates

In this paper, a numerical method for the linear and geometrically non-linear static analysis of thin plates is presented. The method begins with the elasticity equations pertaining to strain components, stresses, displacement components, strain energy and work due to externally applied loads. The plate geometry is defined by a quadrangular boundary with four straight edges and the natural coordinates in conjunction with the Cartesian coordinates are used to map the geometry. The matrix equation of equilibrium is derived using the work-energy principle with the displacement fields expressed by algebraic polynomials, the coefficients of which are then manipulated to satisfy the kinematic boundary conditions. To validate the results from the present method, square plates having all sides fully fixed and all sides simply supported without in-plane movement are analysed. Comparison is made for the uniformly loaded square plate with the results obtained by Levy who solved the non-linear plate bending problem using the Th.von Ka’rma’n’s equations. Rhombic plates are examined and numerical results corresponding to these cases are presented in this paper. Very good comparison of the results regarding deflection and bending stresses with other sources available in the literature is found.     

A new approach to large deflection analysis of heated sandwich plates

This paper represents non-linear analysis of sandwich plates under thermal loading in which a new approach is followed. Analysis of rectangular sandwich plates has been carried out in detail. Numerical results have been computed and compared with other known results.

---

Sommario Questo lavoro illustra l'analisi non lineare delle piastre sandwich sotto carico termico seguendo un nuovo approccio. Si è svolta in dettaglio l'analisi delle piastre rettangolari calcolando valori numerici e paragonandoli con risultati precedentemente noti.

     

Cell size effects for vibration analysis and design of sandwich beams

In this work, sandwich beams are studied to reveal the underlying size effects of the periodic core cells for the first time within the framework of free vibration analysis of such an advanced lightweight structure. The energy equivalence method is formulated as a theoretical approach that takes into account the cell size effect. It is compared with the asymptotic homogenization method and direct finite element method systematically to show their consistence and applicability. The accuracy of free vibration responses predicted by the detailed finite element model is used as the standard of comparison. It is shown that the cell size is an important parameter characterizing the cellular core rigidities that influence vibration responses. The homogenization model agrees exactly with the asymptotic solution of the analytical expression of the beam model only whenever the cell size tends to be infinitely small.     

Design sensitivity analysis for the homogenized elasticity tensor of a polymer filled with rubber particles

The main purpose of this work is the computational simulation of the sensitivity coefficients of the homogenized tensor for a polymer filled with rubber particles with respect to the material parameters of the constituents. The Representative Volume Element (RVE) of this composite contains a single spherical particle, and the composite components are treated as homogeneous isotropic media, resulting in an isotropic effective homogenized material. The sensitivity analysis presented in this paper is performed via the provided semi-analytical technique using the commercial FEM code ABAQUS and the symbolic computation package MAPLE. The analytical method applied for comparison uses the additional algebraic formulas derived for the homogenized tensor for a medium filled with spherical inclusions, while the FEM-based technique employs the polynomial response functions recovered from the Weighted Least-Squares Method. The homogenization technique consists of equating the strain energies for the real composite and the artificial isotropic material characterized by the effective elasticity tensor. The homogenization problem is solved using ABAQUS by the application of uniform deformations on specific outer surfaces of the composite RVE and the use of tetrahedral finite elements C3D4. The energy approach will allow for the future application of more realistic constitutive models of rubber-filled polymers such as that of Mullins and for RVEs of larger size that contain an agglomeration of rubber particles.     

Static, stability and dynamic analysis of gradient elastic flexural Kirchhoff plates

The governing equation of motion of gradient elastic flexural Kirchhoff plates, including the effect of in-plane constant forces on bending, is explicitly derived. This is accomplished by appropriately combining the equations of flexural motion in terms of moments, shear and in-plane forces, the moment–stress relations and the stress–strain equations of a simple strain gradient elastic theory with just one constant (the internal length squared), in addition to the two classical elastic moduli. The resulting partial differential equation in terms of the lateral deflection of the plate is of the sixth order instead of the fourth, which is the case for the classical elastic case. Three boundary value problems dealing with static, stability and dynamic analysis of a rectangular simply supported all-around gradient elastic flexural plate are solved analytically. Non-classical boundary conditions, in additional to the classical ones, have to be utilized. An assessment of the effect of the gradient coefficient on the static or dynamic response of the plate, its buckling load and natural frequencies is also made by comparing the gradient type of solutions against the classical ones.     

Asymptotic analysis of a thin linearly elastic plate equipped with a periodic distribution of stiffeners

We derive several models of thin plates equipped with a periodic distribution of stiffeners. Depending on the orders of magnitude of the different parameters involved, diverse situations arise, from classical Kirchhoff–Love behaviour with additional energy term to full rigidification.     

Wave propagation in periodic buckled beams. Part I: Analytical models and numerical simulations

Periodic buckled beams possess a geometrically nonlinear, load–deformation relationship and intrinsic length scales such that stable, nonlinear waves are possible. Modeling buckled beams as a chain of masses and nonlinear springs which account for transverse and coupling effects, homogenization of the discretized system leads to the Boussinesq equation. Since the sign of the dispersive and nonlinear terms depends on the level of buckling and support type (guided or pinned), compressive supersonic, rarefaction supersonic, compressive subsonic and rarefaction subsonic solitary waves are predicted, and their existence is validated using finite element simulations of the structure. Large dynamic deformations, which cannot be approximated with a polynomial of degree two, lead to strongly nonlinear equations for which closed-form solutions are proposed.     

Finite element analysis of post-buckling dynamics in plates—Part I: An asymptotic approach

Various static and dynamic aspects of post-buckled thin plates, including the transition of buckled patterns, post-buckling dynamics, secondary bifurcation, and dynamic snapping (mode jumping phenomenon), are investigated systematically using asymptotical and non-stationary finite element methods. In part I, the secondary dynamic instability and the local post-secondary buckling behavior of thin rectangular plates under generalized (mechanical and thermal) loading is investigated using an asymptotic numerical method which combines Koiter’s nonlinear instability theory with the finite element technique. A dynamic multi-mode reduction method—similar to its static single-mode counterpart: Liapunov–Schmidt reduction—is developed in this perturbation approach. Post-secondary buckling equilibrium branches are obtained by solving the reduced low-dimensional parametric equations and their stability properties are determined directly by checking the eigenvalues of the resulting Jacobian matrix. Typical post-secondary buckling forms—transcritical, supercritical and subcritical bifurcations are observed according to different combinations of boundary conditions and load types. Geometric imperfection analysis shows that not only the secondary bifurcation load but also changes in the fundamental post-secondary buckling behavior are affected. The post-buckling dynamics and the global analysis of mode jumping of the plates are addressed in part II.     

基于B-样条小波的复合材料层合板灵敏度分析

基于BSWI样条小波有限元方法,将响应和响应的灵敏度系数同时看作状态变量,在Hamilton体系下推导了复合材料层合板响应和响应灵敏度系数的混合控制方程。基于该混合控制方程研究了对称铺层四边固支复合材料层合板的位移响应灵敏度系数在z方向上的分布情况,并将计算结果与有限差分法进行比较。结果表明:材料参数E₁,E₂和G₁₂对位移响应影响明显,其中以E₁对位移响应的影响最大,E₂次之。与经验法相比,半解析法计算结果相对误差小于10^-4,这说明本文的计算方法是正确的,且降低了计算成本和程序实现难度,提高了计算精度。     

Damping for large-amplitude vibrations of plates and curved panels,Part 1: Modeling and experiments

Theoretical and experimental non-linear vibrations of thin rectangular plates and curved panels subjected to out-of-plane harmonic excitation are investigated. Experiments have been performed on isotropic and laminated sandwich plates and panels with supported and free boundary conditions. A sophisticated measuring technique has been developed to characterize the non-linear behavior experimentally by using a Laser Doppler Vibrometer and a stepped-sine testing procedure. The theoretical approach is based on Donnell's non-linear shell theory (since the tested plates are very thin) but retaining in-plane inertia, taking into account the effect of geometric imperfections. A unified energy approach has been utilized to obtain the discretized non-linear equations of motion by using the linear natural modes of vibration. Moreover, a pseudo arc-length continuation and collocation scheme has been used to obtain the periodic solutions and perform bifurcation analysis. Comparisons between numerical simulations and the experiments show good qualitative and quantitative agreement. It is found that, in order to simulate large-amplitude vibrations, a damping value much larger than the linear modal damping should be considered. This indicates a very large and non-linear increase of damping with the increase of the excitation and vibration amplitude for plates and curved panels with different shape, boundary conditions and materials.     

中厚板统计分析的随机边界元法

本文建立了分析含随机材料参数并具厚度不均匀性的中厚板问题的随机边界元法,基于Taylor级数展开技术,分析和到广义位移的均值和一阶偏差的积分方程,其中将材料参数的随机性和厚度的不均匀性作为等效荷载处理,从而得到广义边界位移或面力的均值和协方差,并进一步求出部点广义位移和内力的均值和协方差,最后用本文方法计算了两个数例,并对所得结果进行了分析,探讨。     

The effective elastic moduli of columnar composites made of cylindrically anisotropic phases with rough interfaces

The present work aims to determine the effective elastic moduli of a composite having a columnar microstructure and made of two cylindrically anisotropic phases perfectly bonded at their interface oscillating quickly and periodically along the circular circumferential direction. To achieve this objective, a two-scale homogenization method is elaborated. First, the micro-to-meso upscaling is carried out by applying an asymptotic analysis, and the zone in which the interface oscillates is correspondingly homogenized as an equivalent interphase whose elastic properties are analytically and exactly determined. Second, the meso-to-macro upscaling is accomplished by using the composite cylinder assemblage model, and closed-form solutions are derived for the effective elastic moduli of the composite. Two important cases in which rough interfaces exhibit comb and saw-tooth profiles are studied in detail. The analytical results given by the two-scale homogenization procedure are shown to agree well with the numerical ones provided by the finite element method and to verify the universal relations existing between the effective elastic moduli of a two-phase columnar composite.     

Fracture mechanics of plates and shells applied to fail-safe analysis of fuselage Part I: Theory

In this paper, a general theory on the asymptotic field near the crack tip for plates and shells with and without shear deformation effect is established. It is found that four stress intensity factors, two for symmetrical and antisymmetrical stretching and two for symmetrical and antisymmetrical bending, are required to describe arbitrary asymptotic fields near the crack tip for plates without shear deformation. An additional stress intensity factor is required for the transverse shearing force induced by antisymmetrical bending when the shear deformation is included in the analysis. It is also proven by means of the complex variable technique that for problems of plates with shear deformation, there exist similarities in the asymptotic expressions of moments and membrane forces and also in the asymptotic expressions of in-plane displacements and rotations of the mid-surface. The energy release rate associated with crack growth in the direction of the crack line can be expressed in terms of stress intensity factors by means of Irwin's method of work and energy associated with a virtual crack extension. A combined stress intensity factor can be defined through the total energy release rate. The theory of the fracture of plates is generalized and applied to the study of problems in the fracture of shells. An example of an infinitely long cylindrical shell with a circumferential crack subjected to remote axial tension is given to demonstrate the application of the theory and to test the accuracy of the numerical analysis used for the problem.