A quadrilateral element based on refined global-local higher-order theory for coupling bending and extension thermo-elastic multilayered plates

In this paper a refined higher-order global-local theory is presented to analyze the laminated plates coupled bending and extension under thermo-mechanical loading. The in-plane displacement fields are composed of a third-order polynomial of global coordinate z in the thickness direction and 1,2–3 order power series of local coordinate ζ_k in the thickness direction of each layer, which is identical to the 1,2–3 global-local higher-order theory by Li and Liu [Li, X.Y., Liu, D., 1997. Generalized laminate theories based on double superposition hypothesis. Int. J. Numer. Methods Eng. 40, 1197–1212] Moreover, a second-order polynomial of global coordinate z in the thickness direction is chosen as transverse displacement field. The transverse shear stresses can satisfy continuity at interfaces, and the number of unknowns does not depend on the layer numbers of the laminate.Based on this theory, a quadrilateral laminated plate element satisfying the requirement of C¹ continuity is presented. By solving both bending and thermal expansion problems of laminates, it can be found that the present refined theory is very accurate and obviously superior to the existing 1,2–3 global-local higher-order theory. The most attractive feature of this theory is that the transverse shear stresses can be accurately predicted from direct use of constitutive equations without any post-processing method. It is also shown that the present quadrilateral element possesses higher accuracy.     

Stochastic finite element nonlinear free vibration analysis of piezoelectric functionally graded materials beam subjected to thermo-piezoelectric loadings with material uncertainties

In this paper, second order statistics of large amplitude free flexural vibration of shear deformable functionally graded materials (FGMs) beams with surface-bonded piezoelectric layers subjected to thermopiezoelectric loadings with random material properties are studied. The material properties such as Young’s modulus, shear modulus, Poisson’s ratio and thermal expansion coefficients of FGMs and piezoelectric materials with volume fraction exponent are modeled as independent random variables. The temperature field considered is assumed to be uniform and non-uniform distribution over the plate thickness and electric field is assumed to be the transverse components E _z only. The mechanical properties are assumed to be temperature dependent (TD) and temperature independent (TID). The basic formulation is based on higher order shear deformation theory (HSDT) with von-Karman nonlinear strain kinematics. A C ⁰ nonlinear finite element method (FEM) based on direct iterative approach combined with mean centered first order perturbation technique (FOPT) is developed for the solution of random eigenvalue problem. Comparison studies have been carried out with those results available in the literature and Monte Carlo simulation (MCS) through normal Gaussian probability density function.     

Bending solution of third-order orthotropic Reddy plates with asymmetric interfacial crack

In this paper Reddy’s third-order shear deformable plate theory is applied to asymmetrically delaminated orthotropic composite plates under antiplane–inplane shear fracture mode. A double-plate system is utilized to capture the mechanical behavior of the uncracked plate portion. An assumed displacement field is used and modified in order to satisfy the traction-free conditions at the top and bottom plate boundaries. Moreover, the system of exact kinematic conditions was also implemented into the novel plate model. An important improvement of this work compared to previous papers is the continuity condition of the shear strains at the interface of the double-plate system. Applying these conditions it is shown that the nineteen parameters of the third-order displacement field can be reduced to nine. Using the simplified displacement field the governing equations are derived, as well. The solution of a simply-supported delaminated plate is presented using the state-space model and the displacement, strain and stress fields are determined, respectively. The energy release rate and mode mixity distributions are calculated using the 3D J-integral. The analytical results are compared to those by finite element computations and it is concluded that the present model is the most accurate one among the previous plate theory-based approaches.     

Bending analysis of functionally graded CNT reinforced doubly curved singly ruled truncated rhombic cone

The bending analysis of functionally graded carbon nanotube (CNT) reinforced doubly curved singly ruled truncated rhombic cone is investigated. In this study, a simple C⁰ isoparametric finite element formulation based on third order shear deformation theory is presented. To characterize the membrane-flexure behavior observed in a CNT reinforced truncated rhombic cone, a displacement field involving higher-order terms in in-plane fields is considered. The proposed kinematics field incorporates for transverse shear deformation and nonlinear variation of the in-plane displacement field through the thickness to predict the overall response of the CNT reinforced truncated rhombic cone in an accurate sense. The material properties of the CNT reinforced truncated rhombic cone are estimated according to the rule of mixture. The present model eliminates the need of shear correction factor and imposed zero-transverse shear strain at upper and lower surface of the truncated rhombic cone. The new feature in present model is simultaneous inclusion of twist curvature in strain field as well as curvature in displacement field that makes it suitable for moderately thick and deep truncated rhombic cone. The proposed new mathematical model is implemented in finite element code written in FORTRAN. The proposed model has been validated with analytical, experimental, and finite element results from the literature. This is first attempt to study bending response of CNT reinforced doubly curved singly ruled truncated rhombic cone. The effect of CNT distribution, boundary condition, loading pattern, and other geometric parameters are also examined.     

A nth-order shear deformation theory for the free vibration analysis on the isotropic plates

In the present paper, a nth-order shear deformation theory is used to perform the free vibration analysis of the isotropic plates. The present nth-order shear deformation theory satisfies the zero transverse shear stress boundary conditions on the top and bottom surface of the plate. Reddy??s third order theory can be considered as a special case of present nth-order theory (n=3). The governing equations and boundary conditions are derived by the principle of virtual work. The governing differential equations of the isotropic plates are solved by the meshless radial point collocation method based on the thin plate spline radial basis function. The effectiveness of the present theory is demonstrated by applying it to free vibration problem of the square and circular isotropic plate.     

An advanced higher-order theory for laminated composite plates with general lamination angles

This paper proposes a higher-order shear deformation theory to predict the bending response of the laminated composite and sandwich plates with general lamination configurations.The proposed theory a priori satisfies the continuity conditions of transverse shear stresses at interfaces.Moreover,the number of unknown variables is independent of the number of layers.The first derivatives of transverse displacements have been taken out from the inplane displacement fields,so that the C 0 shape functions are only required during its finite element implementation.Due to C 0 continuity requirements,the proposed model can be conveniently extended for implementation in commercial finite element codes.To verify the proposed theory,the fournode C 0 quadrilateral element is employed for the interpolation of all the displacement parameters defined at each nodal point on the composite plate.Numerical results show that following the proposed theory,simple C 0 finite elements could accurately predict the interlaminar stresses of laminated composite and sandwich plates directly from a constitutive equation,which has caused difficulty for the other global higher order theories.     

A general nonlinear theory of elastic shells

This paper presents a general nonlinear theory of elastic shells for large deflections and finite strains in reference to a certain natural state. By expanding the displacement components into power series in the coordinate θ³ normal to the undeformed middle surface of shells, the expansions of the Cauchy-Green strain tensors are expressed in terms of these expanded displacement components. Through the modified Hellinger-Reissner variational principle for a three-dimensional elastic continuum, a set of the fundamental shell equations is derived in terms of the expanded Cauchy-Green strain tensors and Kirchhoff stress resultants. The Love-Kirchhoff hypothesis is not assumed and higher order stretching and bending are taken into consideration. For elastic shells of isotropic materials, assuming the strain-energy to be an analytic function of the strain measures, general nonlinear constitutive equations are then derived. Thus, a complete and consistent two-dimensional shell theory incorporating the geometrical and physical nonlinearities is established. The classical theories of shells are directly derivable from the present results by proper truncations of the series.     

Large Strain Shear Compression Test of Sheet Metal Specimens

A hybrid experimental-computational procedure to establish accurate true stress-plastic strain curve of sheet metal specimen covering the large plastic strain region using shear compression test data is described. A new shear compression jig assembly with a machined gage slot inclined at 35° to the horizontal plane of the assembly is designed and fabricated. The novel design of the shear compression jig assembly fulfills the requirement to maintain a uniform volume of yielded material with characteristic maximum plastic strain level across the gage region of the Shear Compression Metal Sheet (SCMS) specimen. The approach relies on a one-to-one correlation between measured global load–displacement response of the shear compression jig assembly with SCMS specimen to the local stress-plastic strain behavior of the material. Such correlations have been demonstrated using finite element (FE) simulation of the shear compression test. Coefficients of the proposed correlations and their dependency on relative plastic modulus were determined. The procedure has been established for materials with relative plastic modulus in the range 5?×?10^?4?<?(E _p /E)?<?0.01. It can be readily extended to materials with relative plastic modulus values beyond the range considered in this study. Nonlinear characteristic hardening of the material could be established through piecewise linear consideration of the measured load–displacement curve. Validity of the procedure is established by close comparison of measured and FE-predicted load–displacement curve when the provisional hardening curve is employed as input material data in the simulation. The procedure has successfully been demonstrated in establishing the true stress-plastic strain curve of a demonstrator 0.0627C steel SCMS specimen to a plastic strain level of 49.2 pct.     

A Bending-Gradient model for thick plates,Part II: Closed-form solutions for cylindrical bending of laminates

In the first part (Lebée and Sab, 2010a) of this two-part paper we have presented 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 Bending-Gradient plate theory is an extension to arbitrarily layered plates of the Reissner–Mindlin plate theory which appears as a special case when the plate is homogeneous. Moreover, we demonstrated that, in the general case, the Bending-Gradient model cannot be reduced to a Reissner–Mindlin model. In this paper, the Bending-Gradient theory is applied to laminated plates and its predictions are compared to those of Reissner–Mindlin theory and to full 3D (Pagano, 1969) exact solutions. The main conclusion is that the Bending-Gradient gives good predictions of deflection, shear stress distributions and in-plane displacement 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.     

A Refined Shear Deformation Theory of Bending of Isotropic Plates*

ABSTRACT

A nonlinear, in-plane displacement assumption is proposed, based on an undetermined variation df/dz of transverse shear strains through the plate thickness. A second-order ordinary differential equation for f(z) and two surface conditions, as well as a set of eighth-order partial differential equations and four associated boundary conditions, are derived from the principle of minimum potential energy. Coupling exists between the partial and ordinary differential equations. In the homogeneous solutions for the former, in addition to an interior solution contribution, there exist two edge-zone solution contributions, one of which induces self-equilibrated (in the thickness direction) boundary stresses. Three examples are calculated using the present theory. The last gives the stress couple and maximum-stress concentration factors at the free edge of a circular hole in a large bent plate. Numerical results for the examples are compared with those given by three-dimensional elasticity theory and several two-dimensional theories. It is found that the present theory can accurately predict nonlinear variations of in-plane stresses through the thickness of a plate.     

Existence of Solutions for a Mathematical Model Related to Solid–Solid Phase Transitions in Shape Memory Alloys

We consider a strongly nonlinear PDE system describing solid–solid phase transitions in shape memory alloys. The system accounts for the evolution of an order parameter χ (related to different symmetries of the crystal lattice in the phase configurations), of the stress (and the displacement u), and of the absolute temperature ?. The resulting equations present several technical difficulties to be tackled; in particular, we emphasize the presence of nonlinear coupling terms, higher order dissipative contributions, possibly multivalued operators. As for the evolution of temperature, a highly nonlinear parabolic equation has to be solved for a right hand side that is controlled only in L ¹. We prove the existence of a solution for a regularized version by use of a time discretization technique. Then, we perform suitable a priori estimates which allow us pass to the limit and find a weak global-in-time solution to the system.     

Non-linear theories of beams and plates accounting for moderate rotations and material length scales

The primary objective of this paper is to formulate the governing equations of shear deformable beams and plates that account for moderate rotations and microstructural material length scales. This is done using two different approaches: (1) a modified von Kármán non-linear theory with modified couple stress model and (2) a gradient elasticity theory of fully constrained finitely deforming hyperelastic cosserat continuum where the directors are constrained to rotate with the body rotation. Such theories would be useful in determining the response of elastic continua, for example, consisting of embedded stiff short fibers or inclusions and that accounts for certain longer range interactions. Unlike a conventional approach based on postulating additional balance laws or ad hoc addition of terms to the strain energy functional, the approaches presented here extend existing ideas to thermodynamically consistent models. Two major ideas introduced are: (1) inclusion of the same order terms in the strain–displacement relations as those in the conventional von Kármán non-linear strains and (2) the use of the polar decomposition theorem as a constraint and a representation for finite rotations in terms of displacement gradients for large deformation beam and plate theories. Classical couple stress theory is recovered for small strains from the ideas expressed in (1) and (2). As a part of this development, an overview of Eringen׳s non-local, Mindlin׳s modified couple stress theory, and the gradient elasticity theory of Srinivasa–Reddy is presented.     

Linear and higher order displacement theories for adhesively bonded lap joints

This work presents an adhesive model for stress analysis of bonded lap joints, which can be applied to model thin and thick adhesive layers. In this theory, linear variations of displacement components along the adhesive thickness are firstly assumed, and the longitudinal strain and the Poisson's effect of the adhesive are modeled. A differential form of the equilibrium equations for the adherends is analytically solved by means of compatible relations of the adhesive deformation. The derived shear and peel stresses are compared with the classical adhesive model of continuous springs with constant shear and peel stresses, and validated with two-dimensional finite element results of the geometrically nonlinear analysis using a commercial package. The numerical results show that the present linear displacement theory can be applied to both thin and moderately thick adhesive layers. The present formulation of the linear displacement theory is then extended to the higher order displacement theory for stress analysis of a thick adhesive, whose numerical results are also compared with those of the finite element computation.     

An efficient higher order zigzag theory for composite and sandwich beams subjected to thermal loading

A new efficient higher order zigzag theory is presented for thermal stress analysis of laminated beams under thermal loads, with modification of the third order zigzag model by inclusion of the explicit contribution of the thermal expansion coefficient α₃ in the approximation of the transverse displacement w. The thermal field is approximated as piecewise linear across the thickness. The displacement field is expressed in terms of the thermal field and only three primary displacement variables by satisfying exactly the conditions of zero transverse shear stress at the top and the bottom and its continuity at the layer interfaces. The governing equations are derived using the principle of virtual work. Fourier series solutions are obtained for simply-supported beams. Comparison with the exact thermo-elasticity solution for thermal stress analysis under two kinds of thermal loads establishes that the present zigzag theory is generally very accurate and superior to the existing zigzag theory for composite and sandwich beams.     

A new trigonometric shear deformation theory for isotropic,laminated composite and sandwich plates

A new trigonometric shear deformation theory for isotropic and composite laminated and sandwich plates, is developed. The new displacement field depends on a parameter “m”, whose value is determined so as to give results closest to the 3D elasticity bending solutions. The theory accounts for adequate distribution of the transverse shear strains through the plate thickness and tangential stress-free boundary conditions on the plate boundary surface, thus a shear correction factor is not required. Plate governing equations and boundary conditions are derived by employing the principle of virtual work. The Navier-type exact solutions for static bending analysis are presented for sinusoidally and uniformly distributed loads. The accuracy of the present theory is ascertained by comparing it with various available results in the literature. The results show that the present model performs as good as the Reddy’s and Touratier’s shear deformation theories for analyzing the static behavior of isotropic and composite laminated and sandwich plates.     

Asymptotic construction of Reissner-like composite plate theory with accurate strain recovery

The focus of this paper is to develop an asymptotically correct theory for composite laminated plates when each lamina exhibits monoclinic material symmetry. The development starts with formulation of the three-dimensional (3-D), anisotropic elasticity problem in which the deformation of the reference surface is expressed in terms of intrinsic two-dimensional (2-D) variables. The variational asymptotic method is then used to rigorously split this 3-D problem into a linear one-dimensional normal-line analysis and a nonlinear 2-D plate analysis accounting for classical as well as transverse shear deformation. The normal-line analysis provides a constitutive law between the generalized, 2-D strains and stress resultants as well as recovering relations to approximately but accurately express the 3-D displacement, strain and stress fields in terms of plate variables calculated in the plate analysis. It is known that more than one theory may exist that is asymptotically correct to a given order. This nonuniqueness is used to cast a strain energy functional that is asymptotically correct through the second order into a simple “Reissner-like” plate theory. Although it is not possible in general to construct an asymptotically correct Reissner-like composite plate theory, an optimization procedure is used to drive the present theory as close to being asymptotically correct as possible while maintaining the beauty of the Reissner-like formulation. Numerical results are presented to compare with the exact solution as well as a previous similar yet very different theory. The present theory has excellent agreement with the previous theory and exact results.     

A layerwise C0-type higher order shear deformation theory for laminated composite and sandwich plates

A novel layerwise C⁰-type higher order shear deformation theory (layerwise C⁰-type HSDT) for the analysis of laminated composite and sandwich plates is proposed. A C⁰-type HSDT is used in each lamina layer and the continuity of in-plane displacements and transverse shear stresses at inner-laminar layer is consolidated. The present layerwise theory retains only seven variables without increasing the number of variables when the number of lamina layers are intensified. The shear stresses through the plate thickness derived from the constitutive equation of the present theory have the same shape as those calculated from the equilibrium equation. In addition, the artificial constraints are added in the principle of virtual displacements (PVD) and are certainly fulfilled through a penalty approach. In this paper, two C⁰-continuity numerical methods, such as the Finite Element Method (FEM) and Bézier isogeometric element (BIEM) are utilized to solve a discrete system of equations derived from the PVD. Several numerical examples with various geometries, aspect ratios, stiffness ratios, and boundary conditions are investigated and compared with the 3D elasticity solution, the analytical, as well as, numerical solutions based on various plate theories.     

Stochastic nonlinear bending response of laminated composite plates with system randomness under lateral pressure and thermal loading

In this paper, the effect of sensitivity of randomness in system parameters on the nonlinear transverse central deflection response of laminated composite plates subjected to transverse uniform lateral pressure and thermal loading is examined. System parameters such as the lamina material properties, expansion of thermal coefficients, lamina plate thickness and lateral load are modelled as basic random variables. A higher order shear deformation theory in the von-Karman sense is used to model the system behavior of the laminated plate. A direct iterative-based C ⁰ nonlinear finite element method in conjunction with the first-order perturbation technique developed by the authors is extended for thermal problem to obtain the second-order response statistics, i.e., mean and variance of the nonlinear transverse deflection of the plate. Typical numerical results of composite plates with temperature independent and dependent material properties subjected to uniform temperature and combination of uniform and transverse temperature are obtained for various combinations of geometric parameters, uniform lateral pressures, staking sequences and boundary conditions. The results have been compared with those available in the literature and an independent Monte Carlo simulation.     

Interface crack between isotropic Kirchhoff plates

In this work Kirchhoff plate theory is used to calculate the energy release rate function in delaminated isotropic plates. The approximation is based on the consideration of the equilibrium equations and the displacement continuity between the interface plane of a double-plate model. It is shown that the interface shear stresses are governed by a fourth order partial differential equation system. As an example, a simply supported delaminated plate subjected to a point force is analyzed adopting Lévy plate formulation and the mode-II and mode-III energy release rate distributions along the crack front were calculated by the J-integral. To confirm the analytical results the 3D finite element model of the delaminated plate was created, the energy release rates were calculated by the virtual crack-closure technique and the J-integral. The results indicate a good agreement between analysis and numerical computation.     

Nonlinear bending behavior of orthotropic Mindlin plate resting on orthotropic Pasternak foundation using GDQM

Rectangular plates resting on elastic foundations are operational activities of large transportation aircraft on runways, footings, foundation of spillway dam, civil building in cold regions, and bridge structures. Hence, in the present work, nonlinear bending analysis of embedded rectangular plates is investigated based on orthotropic Mindlin plate theory. The elastic medium is simulated by orthotropic Pasternak foundation. Adopting the nonlinear strain–displacement relation, the governing equations are derived based on energy method and Hamilton’s principle. The generalized differential quadrature method is performed for the case when all four ends are clamped supported. The influences of the plate thickness, shear-locking, elastic medium constants, and applied force on the nonlinear bending of the rectangular plate are studied. Results indicate that increasing the plate thickness decreases the deflection of the plate. It is also observed that increasing the applied force increases the deflection of the plate. Furthermore, considering elastic medium decreases deflection of the plate, and the effect of the Pasternak-type is higher than the Winkler-type on the maximum deflection of the plate. Also, it is found that the present results have good agreement with previous researches.