Next Article in Journal
Weather Impact on Airport Performance
Next Article in Special Issue
Exploitation of a Multifunctional Twistable Wing Trailing-Edge for Performance Improvement of a Turboprop 90-Seats Regional Aircraft
Previous Article in Journal
Flexural Mechanical Properties of Hybrid Epoxy Composites Reinforced with Nonwoven Made of Flax Fibres and Recycled Carbon Fibres
Previous Article in Special Issue
Pneumatically Powered Drilling of Carbon Fibre Composites Using Synthetic Biodegradable Lubricating Oil: An Experimental Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Free and Forced Vibration of Laminated and Sandwich Plates by Zig-Zag Theories Differently Accounting for Transverse Shear and Normal Deformability

Dipartimento di Ingegneria Meccanica e Aerospaziale, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129 Torino, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Aerospace 2018, 5(4), 108; https://doi.org/10.3390/aerospace5040108
Submission received: 13 September 2018 / Revised: 4 October 2018 / Accepted: 8 October 2018 / Published: 11 October 2018
(This article belongs to the Special Issue Adaptive/Smart Structures and Multifunctional Materials in Aerospace)

Abstract

:
A number of mixed and displacement-based zig-zag theories are derived from the zig-zag adaptive theory (ZZA). As a consequence of their different assumptions on displacement, strain, and stress fields, and layerwise functions, these theories account for the transverse shear and normal deformability in different ways, but their unknowns are independent of the number of layers. Some have features that are reminiscent of ones that have been published in the literature for the sake of comparison. Benchmarks with different length-to-thickness ratios, lay-ups, material properties, and simply supported or clamped edges are studied with the intended aim of contributing toward better understanding the influence of transverse anisotropy on free vibration and the response of blast-loaded, multilayered, and sandwich plates, as well as enhancing the existing database. The results show that only theories whose layerwise contributions identically satisfy interfacial stress constrains and whose displacement fields are redefined for each layer provide results that are in agreement with elasticity solutions and three-dimensional (3D) finite element analysis (FEA) (mixed solid elements with displacements and out-of-plane stresses as nodal degrees of freedom (d.o.f.)) with a low expansion order of polynomials in the in-plane and out-of-plane directions. The choice of their layerwise functions is shown to be immaterial, while theories with fixed kinematics are shown to be strongly case-sensitive and often inadequate (even for slender components).

Graphical Abstract

1. Introduction

Fibre-reinforced laminated and sandwich composites are fundamental to obtain a faster speed, longer range, larger payloads, less engine power, and a better operating economy of land, sea, and aerospace vehicles. Indeed, they possess excellent specific strength and stiffness properties along with many other advantages. However, their behaviour may be strongly influenced by local effects in their multi-phase structure, which can cause a relevant loss of strength and stiffness, and even lead to a possible catastrophic failure in service. To avert this, a three-dimensional description of their displacement and stress fields is required in the computer simulations. Displacements are no longer C1-continuous as in homogeneous non-layered materials, but instead they should be C0-continuous (zig-zag effect), this being the only way to guarantee the continuity of transverse shear and normal stresses and of the transverse normal stress gradient that is necessary to satisfy the fulfillment of local equilibrium equations in layered materials.
A multitude of theories is currently available for laminated and sandwich composites, which account for layerwise effects at varying degrees of accuracy and computational costs. The papers by Carrera et al. [1,2] and by Demasi [3] are cited wherein a broad discussion of this matter can be found. Theories can be summarily categorized into equivalent single-layer (ESL), discrete-layer (DL), and zig-zag (ZZ) formulations, and further into displacement-based and mixed theories, because displacements, strains, and stress fields can be chosen separately from one another.
Given the limited accuracy offered by ESL even to predict overall response quantities for certain loading, material properties, and stack-up [4,5], they still remain widespread by virtue of their low computational cost (see Burlayenko et al. [6] and by Jun et al. [7]), but nowadays, DL and ZZ are becoming more and more widespread in their applications. Indeed, they have the merit of accurately predicting the displacement and stress fields irrespective of lay-up, layer properties, and loading or boundary conditions. However, DL could overwhelm the computational capacity when structures of industrial interest are analysed (unless their use is limited to critical areas) due to having too many variables. Instead, ZZ having intermediate characteristics between ESL and DL can strike the right balance between accuracy and cost-saving, allowing designers’ demand of theories to be met in a simple, already accurate form. In particular, ZZ accounting for the transverse normal deformability effect has been proven to be suitable for carrying out ply-level stress analyses at a lower cost than DL.
Relevant examples of dynamic studies carried out by DL and ZZ that highlight their superior performance include the papers by [1,2,3,6,7] Boscolo and Banerjee [8], Khdeir and Aldraihem [9], Sayyad and Ghugal [10], Kazanci [11], Lin and Zhang [12], Vescovini et al. [13], and in addition to the previous ones, [14,15,16,17,18,19] are also cited. Moreover, papers [1,2,3,4,5,8,13,20] are examples that prove the limitations of ESL and the importance of the transverse normal deformability under static localised loading for certain boundary conditions and lay-ups [21] and even for accurately predicting the first free-vibration mode frequencies. Although ZZ theories are ultimately finding an ever-increasing number of applications thanks to their accuracy at an affordable cost, the literature shows that undeservedly Di Sciuva’s zig-zag-like fuction DZZ is less used than Murakami’s zig-zag-like fuction MZZ (for a definition of acronyms, see Table 1), despite their better accuracy with a low order of expansion of analytical solutions, as also demonstrated by the numerical results of this paper.
To contribute to the dissemination of DZZs and their further development, hereafter the typologies and specific characteristics of available ZZs are briefly reviewed. The intended aim is to explain the different behaviour of DZZ and MZZ as a consequence of their different layerwise contributions, namely Di Sciuva’s [22] (DZZ) or Murakami’s [23] (MZZ) zig-zag functions.
DZZ incorporates layerwise contributions as the product of assumed zig-zag functions and zig-zag amplitudes, the latter being computed to a priori satisfy the interfacial compatibility of out-of-plane stresses (for this reason, they are referred in literature as physically-based ZZ). Instead, MZZ assumes zig-zag functions that a priori feature a periodic change of the slope of displacements at interfaces, irrespective of the stack-up. As a consequence, stresses are assumed apart from the kinematics within the framework of Hellinger–Reissner variational theorem. For this reason, they are called kinematic-based ZZ. So, the development of MZZ is easier, and an efficient C0 formulation can more easily be generated, but they could not appropriately represent the displacement fields when the orientation angle of layers is aperiodic, as shown in the literature. According to findings by [1,2,5,24,25,26], it can be stated in general that MZZ are accurate only with a rather high expansion order of variables across the thickness and in the in-plane directions. As shown in [27,28,29,30,31,32,33], instead, DZZ accurately predict displacements and (very often but not in all cases) stresses with a lower expansion order.
As shown in the papers by Carrera et al. [1,2,24,25,26], Demasi [3], Mattei and Bardella [27], Icardi et al. [28,29], Li and Liu [30], Zhen and Wanji [31,32], and Shariyat [33], the transverse normal deformability, which is usually neglected because it is erroneously considered of little importance, on the contrary should be accounted for. Indeed, his contribution becomes of primary importance when boundary conditions other than simply-supported edges are considered, under a strong variation of mechanical properties across the thickness and for localised loading. However, a rather complex formulation of theories, i.e., with a sufficient number of parameters that can be defined in order to satisfy all of the relevant physical constraints, is required for these cases. Although very often power and Taylor’s series expansions are used across the thickness, hierarchic polynomials, and trigonometric and exponential functions, a combination of both or radial basis functions could be required to more efficiently represent variables [34,35,36,37,38,39].
The idea of separately hiring kinematics and stresses, so as to account for the transverse normal deformation without having to represent in a piecewise way the transverse displacement, was applied by Barut et al. to DZZ [40]. Other mixed DZZ that are based on a similar idea have a rather simple kinematics, but yet, the capability to quite accurately predict stresses has been formerly developed in HR form by Kim and Cho [41] and Tessler et al. [42]
Despite the notable contributions mentioned up to now, further research is required on ZZ because recent studies by Zhen and Wanji [43], Gherlone [44], and Groh and Weaver [45] have shown that MZZ could be less accurate than RZT [42], assuming the same degree of representation across the thickness, while [43] shows just the opposite for different cases. So, it is necessary to ascertain even whether three-dimensional (3D) DZZ can be as accurate as MZZ with a minor computational burden, so as to be effectively suitable for industrial purposes, considering also other theories in addition to those already examined in [44,45], especially those with a piecewise representation of the transverse displacement.
A massive amount of free vibration studies is currently found in the literature for cross-ply laminates, but other boundary conditions and lay-up that are equally interesting for practical applications need to be considered when testing the accuracy of theories (see, e.g., Li et al. [46]). It should be noted that for clamped edges, an incorrect vanishing transverse shear force resultant is predicted by ESL plate models; therefore, it is necessary to check on an adequate number of cases as to whether DZZ are immune from such mistakes. Also, lay-ups with a quite large variation of properties of layers, and in particular soft-core sandwiches, should be investigated to enhance the existing database. Another subject that needs further studies is what happens when zig-zag functions that are different from those commonly used [22,23] are assumed. It is also necessary to better clarify what effects have an a priori assumption of certain zig-zag amplitudes at certain interfaces, or even to all of the interfaces such as in MZZ, in order to save costs.
To contribute to this matter, in this paper, zig-zag theories in displacement-based and mixed form are particularized from the ZZA 3D zig-zag theory [29], assuming different layerwise functions. Consider that ZZA, here schematically defined as the adaptive zig-zag theory, assumes a variable kinematics with fixed degrees of freedom (d.o.f.) that consist of in-plane and transverse displacements u α 0 , w 0 and transverse shear rotations Γ α 0 at a middle reference plane. Moreover, it should be noticed that ZZA is capable of great accuracy, even with strong anisotropy at a cost that is still comparable to that of ESL.
The purpose of this paper is to test the accuracy of theories derived from ZZA under simplified hypotheses on displacement, deformations, and stress fields, so as to further reduce the computational cost, in order to assess whether and when the accuracy of ZZA can be preserved. The accuracy and the efficiency of these theories, which for the most part allow a redefinition of the coefficients layer-by-layer so as to satisfy all of the physical constraints and account for the normal transverse deformability effect (although in different ways), are assessed considering the free vibration and the response behaviour of blast-loaded multilayered and sandwich plates.
A new theory referred to as ZZA* is developed from ZZA in order to prove that once the coefficients of displacements are redefined across the thickness, the same accuracy degree can be achieved irrespective of the layerwise functions that are chosen. It will also be shown that as the coefficient is redefined, the same results are obtained even when zig-zag contributions do not explicitly appear. Other theories have been previously developed by the authors [47] in displacement-based and mixed form, either with Di Sciuva’s or Murakami’s zig-zag-like functions, and which also have characteristics similar to other theories that are already known in the literature, are considered for comparison.
Numerical applications aim to show the superiority of ZZA* over ZZA, including its lower computational costs with the same accuracy, for the quite large range of variation of lay-ups and boundary conditions that are considered. Benchmarks that are either retaken from the literature or new are considered, the former in order to enable a comparison of available theories, and the latter in order to enhance the existing database. The results by the present theories are compared to exact solutions whenever available or to 3D finite element analysis (FEA) [48]. The findings can be categorized as the new results of known benchmarks provided by the new theories and entirely new results for the new benchmarks, which could serve as test beds for future analytical and finite element models.
This study also aims to show the superiority of ZZA* over all of the other theories considered for most of the cases, since the use of simplified assumptions implies a loss of accuracy. It will be shown that the accuracy of the simplified theories, in particular those of Murakami, is strongly case-sensitive; therefore, they cannot be used interchangeably. Studies are carried out in closed form considering the free and forced vibrations of laminated and sandwich beams and plates, simply-supported and clamped edges, different lay-ups, and distinctly different material properties of constituent layers, so as to give rise to relevant 3D effects.

2. Preliminaries and Notations

Constituent layers are assumed to have a uniform arbitrary thickness h k and linear elastic properties. They are assumed to be perfectly bonded to each other and, as usual when the global-scale response is examined, the existing bonding resin interlayer is disregarded. For the same reason, sandwiches are described in homogenised form as multilayered structures with one or more thick soft intermediate layers as the cores, with the cell-scale effects being disregarded.
A rectangular, right-handed Cartesian coordinate reference system ( x , y , z ) is assumed as the reference frame, having ( x , y ) on the middle reference plane Ω of the laminated plate (origin in the lower left edge) and z as the thickness coordinate ( z [ 0.5 h , 0.5 h ] , with h being the overall thickness). Lx and Ly symbolise the plate side-length in the x and y-directions, while ( k ) z + and ( k ) z represent the upper and lower positions of the layer interfaces, respectively. Subscripts k and superscripts k are used to indicate that a quantity belongs to the layer k , while u and l mark the upper and lower faces of the laminate, and a comma is used to indicate spatial derivatives, e.g., ( . ) , x = / x , ( . ) , z = / z . The elastic in-plane and transverse displacement components are indicated as u α and u ς , respectively. Strains are assumed to be infinitesimal and, to distinguish their origin, they are specified as ε i j u = 1 / 2 ( u i , j + u j , i ) , ε i j σ = ( E i j k l ) 1 σ k l , respectively when they come from kinematic [ . ] u or stress-strain [ . ] σ relations. Once assumed as primary variables, they are indicated as ε i j , ε α ς , ε ς ς ( i , j 1 , 2 ; α x , β y ; 3 ς z ; γ α β = 2 ε α β ). Stresses from stress–strain constitutive relations are indicated as σ i j ε = E i j k l ε k l ( C i j k l = E i j k l 1 ) , while when they are assumed as primary variables, they are indicated as σ i j , σ α ς , σ ς ς .

2.1. Recalls on Mixed Variational Theorems

Here a generalized version of Hu-Washizu theorem is used, whose primary displacement boundary condition link is weakened as S u ( u i u ˜ i ) n j σ i j d S = 0 and which generates the following variational statement:
Π H W g = V [ ( ε i j u ε i j ) σ i j + ( σ i j ε σ i j ) ε i j + σ i j ε i j u b i u i ] d v S t t ˜ i u i d s S u [ ( u i u ˜ i ) n j σ i j + σ i j n j u i ] d s = 0
Here H W g represents the canonical functional, n j are the components of the external unit normal to the volume bounding surface, b i are the components of body forces (they will contain inertial forces, as explained in (4)) and V ( ε i j u ε i j ) σ i j d v , V ( σ i j ε σ i j ) ε i j d v ensure the consistency of assumed strain and stress fields with their counterparts obtained from stress-strain and strain-displacement relations. As usual, the prismatic volume V of the laminated plates is assumed to be bounded by a surface S that is split into a surface S t on which surface tractions are prescribed and a surface S u on which surface displacements are prescribed. Body forces b i on V , prescribed surface tractions t ˜ i on S t and prescribed displacements u ˜ i on S u are assumed to act.
Theories with displacements and out of plane stresses assumed separately are developed from the HR variational theorem:
Π H R = V [ σ i j ε i j u + ( γ i 3 u γ i 3 σ ) σ i 3 + ( ε 33 u ε ^ 33 σ ) σ 33 b i u i ] d v S t t ˜ i u i d s S u [ ( u i u ˜ i ) n j σ i j + σ i j n j u i ] d s = 0
where δ Π H R is the first variation of the HR canonical functional, and the following assumptions are made: γ i j ( i 3 ) ( 33 ) = 2 ε i j ( i 3 ) ( 33 ) ; ε ^ 33 σ = 1 / C 3333 ( σ 33 C 33 i j ε i j u ) , being i , j = 1 , 2 x , y ; 3 z .
Dynamic governing equations are obtained from the previous variational statements (1), (2), or from the principle of virtual work (PVW), for theories in displacement-based form, accounting for the work of inertial forces.

2.2. Construction of Analytical Solutions

Closed-form solutions to dynamic governing equations are obtained, irrespective of the theory examined, expressing functional d.o.f. as a truncated series expansion of unknown amplitudes A Δ i and trial functions i ( x , y ) that individually satisfy the prescribed boundary conditions:
Δ = i = 1 m Δ A Δ i i ( x , y )
Then, we substitute these expressions within PVW, HR, or HW functionals and operating as specified immediately below. Here, Δ symbolises in turns u 0 , v 0 , w 0 , θ x , θ y , because middle plane displacements and rotations of the normal are assumed as the only functional d.o.f. for each theory of this paper. Mechanical boundary conditions are accounted for by determining a number of unknown amplitudes A Δ i , in proportion to the number of boundary conditions enforced, using the Lagrange multipliers method to account for the relationships resulting from each mechanical condition. The remaining amplitudes are determined deriving the governing functional with respect to still unknown amplitudes and equating to zero, having considered the work of inertial forces:
V [ b i u i ] d v = V [ ρ u ¨ i u i ] d v
within the functionals. In this way, an algebraic system is obtained whose solution provides the numerical value of each amplitude; then, the displacement, strain, and stress fields can be computed.
Since symbolic calculus being used to construct the theories, the applied distributed loading can be conveniently defined as a continuous or discontinuous general function χ ( x , y ) acting on upper and/or lower faces (or just on a part of them), so that energy contributions can be constructed in exact closed form. In this way, a series expansion representation with a large number of components is not necessary to represent discontinuous or otherwise complex loading distributions; then, the construction of the structural model in the numerical applications can be simplified, and at the same time made more accurate.
The trial functions have been adopted in individual applications are explicitly defined below and in Table 2 along with the expansion order used in each case. Note that although it is different for each benchmark, the order of expansion used is the same for all theories when a specific one is examined being equal to the minimum one that makes ZZA accurate. This is done to homogenise results and in order to compare theories under the same conditions.
At the clamped edge of cantilever beams, hereafter assumed at x = 0 by way of example, the following conditions are enforced
u 0 ( 0 , 0 ) = 0 ;   w 0 ( 0 , 0 ) = 0 ;   w 0 ( 0 , 0 ) , x = 0 ;   Γ x 0 ( 0 , 0 ) = 0 ;
the following further conditions are enforced:
u α ( 0 , z ) , z = 0 ;   u ς ( 0 , z ) , z = 0 ;   u ς ( 0 , z ) , x z = 0
to simulate that (5) holds identically across the thickness. To ensure that the transverse shear stress resultant force equals the constraint force, the following additional constraint
h / 2 h / 2 σ x z ( 0 , z ) d z = T
is enforced, while at the free edge x = L such resultant force is enforced to vanish
h / 2 h / 2 σ x z ( L , z ) d z = 0
It should be noted that the boundary conditions in Equations (6) to (8) are enforced using Lagrange multipliers methods; no further conditions are enforced on bending moments, but if necessary, they could be coerced choosing the sufficient expansion order in Equation (3).
At the supported edge of the propped cantilever beams at x = L, the following support condition:
w 0 ( L , h / 2 ) = 0
is enforced at the lower face z = h / 2 , while the condition in Equation (8) is reformulated as:
h / 2 h / 2 σ x z ( L , z ) d z = T L
At the simply supported edges, the following boundary conditions are enforced
w 0 ( 0 , y ) = 0 ; w 0 ( L x , y ) = 0 ; w 0 ( 0 , y ) , x x = 0 ; w 0 ( L x , y ) , x x = 0 w 0 ( x , 0 ) = 0 ; w 0 ( x , L y ) = 0 ; w 0 ( x , 0 ) , y y = 0 ; w 0 ( x , L y ) , y y = 0
on the reference mid-plane of plates at x = 0 , x = L x , and y = 0 , y = L y . Appropriate changes corresponding to the ones for simply supported beams are obtained. Table 2 provides the trial functions and expansion order that is assumed for each benchmark.
Any other boundary condition could be enforced in the same way, namely by choosing the trial functions that individually satisfy the prescribed boundary conditions. In the event of this not being possible and in order to satisfy the mechanical conditions, Lagrange multipliers method should be applied. In the numerical applications for cases g and h, the boundary condition to be further enforced is u ς ( x , h / 2 ) = u ς ( x , h / 2 ) , because some of their first free vibration modes exhibit anti-symmetric displacements in the transverse direction. In this way, a closed-form solution can be obtained with a reduced computational effort, instead of resorting to FEA, as is often done for such cases.

3. Higher-Order Theories

The laminated plate theories that are used as structural models are thoroughly examined below. Here, they are grouped into higher-order and lower-order ones. The latter represent simplified versions of the former, having mostly features that are similar to those of the theories that have been previously proposed in the literature. In the following section, it will be specified which are new and which have been previously developed by the authors. Governing equations will not be reported in explicit form, as they can be obtained in a straightforward way using standard techniques.

3.1. Features of the ZZA Theory and the Higher Order Theories Derived from It

As it forms the basis of all theories considered in this paper, the theoretical framework of ZZA below is expounded first. The through-thickness displacement field is postulated as [29]:
u α ( x , y , z ) = [ u α 0 ( x , y ) + z ( Γ α 0 ( x , y ) w 0 ( x , y ) , α ) ] 0 + [ F α u ( z ) ] i [ k = 1 n i Φ α k ( x , y ) ( z z k ) H k ( z ) + j = 1 n α C u j ( x , y ) H j ( z ) ] c
u ς ( x , y , z ) = [ w 0 ( x , y ) ] 0 + [ F ζ ( z ) ] i + [ k = 1 n i Ψ k ( x , y ) ( z z k ) H k ( z ) + k = 1 n i Ω k ( x , y ) ( z z k ) 2 H k ( z ) + j = 1 n C ς j ( x , y ) H j ( z ) ] c
Three kinds of contributions that are distinctly separated into lower- [ ] 0 , higher- [ ] i and layerwise [ ] c are incorporated, whose specific purpose is described below.
  • [ ] 0 is a linear contribution and contains the five functional degrees of freedom of the theory.
  • [ ] i contains higher-order terms. Any combination of independent functions could be assumed to represent [ F α u ( z ) ] i and [ F ζ ( z ) ] i ; however to include theory [28] as a particularization of ZZA [29], the following form is chosen:
    [ F α u ( z ) ] i = [ C α i ( x , y ) z 2 + D α i ( x , y ) z 3 + ( O z 4 ) ] i = [ 3 ( . ) ˜ α ] i + [ ( O z 4 ) ] i
    [ F ζ ( z ) ] i = [ b i ( x , y ) z + c i ( x , y ) z 2 + d i ( x , y ) z 3 + e i ( x , y ) z 4 + ( O z 5 ) ] i = [ 4 ( . ) ˜ ζ ] i + [ ( O z 5 ) ] i
    Higher-order contributions [ ( O z 4 ) ] i , + [ ( O z 5 ) ] i are characteristic of ZZA, while the terms [ 3 ( . ) ˜ α ] i , [ 4 ( . ) ˜ ζ ] i are the same as in the previous theory [28]. The closed-form expressions of coefficients C α i , D α i , b i to e i are obtained using symbolic calculus from enforcing the fulfilment of stress boundary conditions
    σ α ς = σ ς ς , ς = 0 ; σ ς ς = p 0 ( ± )
    Here, p 0 ( ± ) represents the distributed transverse loading acting on the upper p 0 ( + ) and lower p 0 ( ) faces. Of course, also, non-homogeneous conditions σ α ς ; σ β ς 0 could be enforced without any additional difficulty. For clarity, contributions [ ] i from Equations (14) and (15) are rearranged in the following way:
    U α i ( x , y , z ) = [ A α 2 z 2 + A α 3 z 3 ] + A α 4 z 4 + + A α n z n U ζ i ( x , y , z ) = [ A ς 1 z + A ς 2 z 2 + A ς 3 z 3 + A ς 4 z 4 ] + A ς 5 z 5 + + A ς n z n
    The lower-order contributions under the square brackets in Equation (17) are determined by enforcing the fulfilment of the boundary conditions in Equation (16) and the local equilibrium equations:
    σ α α , α + σ α β , β + σ α z , z = b α ;   σ α ς , α + σ β ς , ς + σ ς ς , ς = b ς
    at selected points across the thickness. Higher-order contributions A α 4 z 4 + + A α n z n , A ς 5 z 5 + + A ς n z n , which enable a variable-kinematics representation across the thickness, are computed for each fictitious computational layer i in which the laminate is subdivided by imposing the fulfilment of Equation (18). However, except otherwise stated, a third/fourth-order representation that embraces the whole laminate is used in the applications which is adequate to obtain accurate results. Indeed, as shown in [21], this choice is a valuable combination of accuracy and cost-saving also when extreme variations of material properties give rise to very strong layerwise effects. A single computational layer is used for laminates, while two or three layers are used for sandwiches. Note that the in-plane position of equilibrium points can be chosen appropriately for each case.
  • [ ] c represents the layerwise contributions; the expressions of zig-zag amplitudes Φ α k , Ψ k and Ω k are determined so that the continuity of out-of-plane stresses and the transverse normal stress gradient σ ς ς , ς at the layer interfaces is satisfied, as prescribed by the elasticity theory through the following stress compatibility conditions:
    σ α ς ( z + ( k ) ) = σ α ς ( z ( k ) ) ;   σ ς ς ( z + ( k ) ) = σ ς ς ( z ( k ) ) ;   σ ς ς , ς ( z + ( k ) ) = σ ς ς , ς ( z ( k ) )
    at the physical and mathematical layer interfaces. Layerwise contributions c u j α and C ς j restore the continuity of displacements at the mathematical layer interfaces. The symbols n i and n in the summation of Equations (12) and (13) are used to distinguish the number of physical interfaces from that of the mathematical layer interfaces, respectively. In detail, Φ α k enables the continuity of transverse shear stresses, while Ψ k , Ω k enable the continuity of the transverse normal stress and of its gradient at physical and mathematical layer interfaces. All together, these provide the right slope changes of displacements at the interfaces of layers with different material properties and/or orientations.
Elsewhere in this paper, the symbols − and + indicate the position just before and just after the interface, respectively. The term ( z z k ) H k appearing in in Equations (12) and (13) is Di Sciuva’s zig-zag function [22], while ( z z k ) 2 H k is Icardi’s parabolic zig-zag function [28], with Hk being the Heaviside unit step function (Hk = 0 for z < zk, while Hk = 1 for zzk).
The enforcement of the equilibrium and stress compatibility conditions in Equations (18) and (19) yields to a system of algebraic equations at each interface that is solved in closed form once and for all using a symbolic calculus tool. Notice that if just the material properties and/or the orientation of layers change, but not their number, symbolic expressions representing the solution will remain the same.
It will be shown forward that ( z z k ) H k and ( z z k ) 2 H k can be replaced with any other layerwise function keeping the results unchanged, provided that coefficients inside [ ] i are re-computed for each computational layer, as outlined above. Layerwise functions can even be omitted if a sufficient number of coefficients are incorporated in [ ] i whose expressions are determined by enforcing the fulfilment of the interfacial stress compatibility conditions in Equation (19). To show this, a new theory ZZA* is developed in Section 3.1.1, which is devoid of [ ] c , but incorporates new unknown coefficients within [ ] i , which will be indicated as [ ] i + c . Since this choice speeds up the computations of the coefficients of each layer, it turns into a computational advantage that grows with the number of computational layers, so it is worth taking into consideration.
The expressions of C u j α and C w j are determined in a straightforward way by enforcing the continuity of displacements at the mathematical layer interfaces
u α ( z + ( k ) ) = u α ( z ( k ) ) ;   u ς ( z + ( k ) ) = u ς ( z ( k ) )
Notice that as no d.o.f. derivatives are involved, the computation of C u j α and C w j is much easier and faster than those of Φ α k , Ψ k , and Ω k , which instead involve such derivatives. However, it must be considered that the computation of all of the above indicated zig-zag terms takes only an infinitesimal fraction of the overall calculation cost, and so remains compatible with that of ESL.
The strain energy updating technique (SEUPT) [29] can be used to obtain a C0 formulation of the ZZA theory, since derivatives of the d.o.f. are involved in the displacement field as a consequence of the enforcement of Equations (18) and (19).

3.1.1. ZZA* Displacement-Based Theory

This theory, which represents the new theoretical contribution brought by this paper, is developed with the intended aim to demonstrate that the choice of zig-zag functions is immaterial, provided that the coefficients of displacements [ ] i are recomputed as indicated above.
For this purpose the displacement field of ZZA* is assumed to be similar to that of ZZA except for the layerwise functions. In numerical applications it will be shown that the same accuracy of ZZA can be achieved without explicitly incorporating contributions [ ] c , and consequently obtaining a reduction of the computational burden.
It will be shown by the numerical results that the choice of zig-zag functions is immaterial if the coefficients [ ] i are redefined layer-by-layer, so these functions do not even need to be explicitly incorporated in the displacement field. The same result wouldn’t be achieved by keeping the coefficients fixed across the thickness, in which case the accuracy depends on the choice of zig-zag functions.
The displacement field of ZZA* is conceived in the following way:
u α ( x , y , z ) = [ u α 0 ( x , y ) + z ( Γ α 0 ( x , y ) w 0 ( x , y ) , α ) ] 0 + { k = 1 n B ˜ α i k ( x , y ) z + + [ C α i ( x , y ) z 2 ] + [ D α i ( x , y ) z 3 ] + k = 1 n i C ˜ α i k ( x , y ) } i + c
u ς ( x , y , z ) = [ w 0 ( x , y ) ] 0 + { [ b i ( x , y ) z + k = 1 n b ˜ i k ( x , y ) z ] + [ c i ( x , y ) z 2 + k = 1 n i c ˜ i k ( x , y ) z 2 ] + [ d i ( x , y ) z 3 ] + e i ( x , y ) z 4 + k = 1 n i d ˜ i k ( x , y ) } i + c
wherein the terms B ˜ α i k and C ˜ α i k serve the same purpose as Φ α k and α C u j in Equation (12) inside ZZA, while b ˜ i k and c ˜ i k have the same function of Ω k and Ψ k , and d ˜ i k has the function of C ς k in Equation (13). In the same way of ZZA, C α i , D α i , b i , c i , d i and e i allow the stress boundary conditions in Equation (16) and the local equilibrium equations in Equation (18) to be met. More specifically, b i and c i enable the fulfillment of the stress boundary conditions concerning σ ς ς and σ ς ς , ς over the lower bounding face, while they are cancelled in subsequent layers. Instead, C α i , D α i , d i , and e i allow satisfying the three equilibria (18) at two points for each intermediate layer. In the lower layer, such coefficients enable the two boundary conditions on σ α ς to be enforced along with the three equilibrium equations at a single point. In this way, free variables still remain that allow meeting three equilibrium equations at a single point across the upper layer and the boundary conditions at the upper bounding surface. When more equilibrium points are desired, each layer can be subdivided into two or more computational layers. However, this doesn’t mean an increased expansion order or number of variables, since coefficients can be recomputed using the same order of representation for each computational layer, and the d.o.f. remain fixed. As for ZZA, closed-form expressions of all of the coefficients of ZZA* within contributions [ ] i + c are determined once and for all for a specific lay-up using a symbolic calculus tool.

3.2. HWZZ Mixed Theory

Such a theory was developed in [47] in order to reduce the computational effort of ZZA by keeping only the essential contributions to the displacement, strain, and stress fields, and hopefully preserving its accuracy. Developed within the framework of the Hu–Washizu theorem, HWZZ forms the basis of another mixed theory here used and called HWZZM, which considers a different zig-zag layerwise function.

3.2.1. Master Displacements

The displacements of HWZZ are derived from those of ZZA neglecting the contributions of Ω k , because they are supposed to give imperceptible slope variations compared to those by Φ α k , Ψ k . Also higher-order and adaptive contributions A α 4 z 4 + + A α n z n and A ς 5 z 5 + + A ς n z n are neglected and no decomposition into mathematical layers is allowed. As a consequence, contributions by α C u j and C ς j are omitted. So, the displacement field of HWZZ is written:
u α ( x , y , z ) = [ u α 0 ( x , y ) + z ( Γ α 0 ( x , y ) w 0 ( x , y ) , α ) ] 0 + [ C α i ( x , y ) z 2 + D α i ( x , y ) z 3 ] i + [ k = 1 n i Φ α k ( x , y ) ( z z k ) H k ( z ) ] c u ς ( x , y , z ) = [ w 0 ( x , y ) ] 0 + [ b i ( x , y ) z + c i ( x , y ) z 2 + d i ( x , y ) z 3 + e i ( x , y ) z 4 ] i + [ k = 1 n i Ψ k ( x , y ) ( z z k ) H k ( z ) ] c

3.2.2. Master Strains

Out-of-plane strains ε z z , γ x z , γ y z are constructed on the basis of the ZZA layer-by-layer representation of displacements as:
u α ( x , y , z ) = [ u α 0 ( x , y ) + z ( Γ α 0 ( x , y ) w 0 ( x , y ) , α ) ] 0 + [ C α i ( x , y ) z 2 + D α i ( x , y ) z 3 ] i + [ k = 1 n i Φ α k ( x , y ) ( z z k ) H k ( z ) + j = 1 C u j α ( x , y ) H j ] c w ς ( x , y , z ) = [ w 0 ( x , y ) ] 0 + [ b i ( x , y ) z + c i ( x , y ) z 2 + d i ( x , y ) z 3 + e i ( x , y ) z 4 ] i + [ k = 1 n i Ψ k ( x , y ) ( z z k ) H k ( z ) + j = 1 C ς j ( x , y ) H j ] c
The symbol ( . ) states that they refer to the computational layer . On the contrary, no decomposition into mathematical layers is allowed for the in-plane strains:
ε x x ( x , y , z ) = U ˜ ( x , y , z ) , x + k = 1 s Φ x , x k ( z z k ) H k ε y y ( x , y , z ) = V ˜ ( x , y , z ) , y + k = 1 s Φ y , y k ( z z k ) H k ε z z ( x , y , z ) = W ˜ ( x , y , z ) , z + k = 1 s Ψ k H k γ x z ( x , y , z ) = [ U ˜ ( x , y , z ) , z + k = 1 s Φ x k H k + W ˜ ( x , y , z ) , x + k = 1 s Ψ , x k ( z z k ) H k ] γ y z ( x , y , z ) = [ V ˜ ( x , y , z ) , z + k = 1 s Φ y k H k + W ˜ ( x , y , z ) , y + k = 1 s Ψ , y k ( z z k ) H k ] γ x y ( x , y , z ) = [ U ˜ ( x , y , z ) , y + k = 1 s Φ x , y k ( z z k ) H k + V ˜ ( x , y , z ) , x + k = 1 s Φ y , x k ( z z k ) H k ]
The expressions of membrane stresses σ x x , σ y y , and σ x y are obtained in a straightforward way from stress–strain relations and previous strains, while those of the out-of-plane counterparts are assumed as specified next.

3.2.3. Master Stresses

Out-of-plane stresses are obtained from membrane stresses by integrating local equilibrium equations:
σ x z = h / 2 h / 2 ( b x σ x x , x σ x y , y ) d z ; σ y z = h / 2 h / 2 ( b y σ x y , x σ y y , y ) d z ; σ z z = h / 2 h / 2 ( b z σ x z , x σ y z , y ) d z
and then
σ z z = h / 2 h / 2 { b z [ h / 2 h / 2 ( b x , x σ x x , x x σ x y , x y ) d z ] [ h / 2 h / 2 ( b y , y σ x y , x y σ y y , y y ) d z ] } d z
In this way, stress jumps resulting from the omission of contributions by Ω k are recovered. As a consequence of simplifying the assumptions that are made, HWZZ needs to be post-processed by ZZA in order to accurately represent displacements when very strong layerwise effects rise.

3.3. Other Mixed Theories

3.3.1. HWZZM Theory

Murakami’s zig-zag contributions are assumed [51], which are variations of the canonical form
M k ( z ) = ( 1 ) k ζ k
where:
ζ k = a k z b k , a k = 2 z k + 1 z k , b k = z k + 1 + z k z k + 1 z k
In fact, the displacement field of HWZZM is assumed as:
u α ( x , y , z ) = [ u α 0 ( x , y ) + z ( Γ α 0 ( x , y ) w 0 ( x , y ) , α ) ] 0 + [ F α u ( z ) ] i +     [ A k u α ( z ) [ 2 z z k + 1 z k z k + 1 + z k z k + 1 z k ] + C α k ( x , y ) ] c
u ς ( x , y , z ) = [ w 0 ( x , y ) ] 0 + [ F ζ ( z ) ] i +     [ A k u ς ( z ) [ 2 z z k + 1 z k z k + 1 + z k z k + 1 z k ] + B k u ς ( z ) [ ( 2 z ) 2 z k + 1 z k ] + C ς k ( x , y ) ] c
Differently to Murakami’s theories proposed in the literature, here, multiplier coefficients A k u ε , A k u ς , and B k u ς are incorporated, which can be defined differently across the thickness to improve the accuracy. As has been done previously, their expressions are a priori determined by enforcing the continuity of the transverse shear and normal stress, and of the transverse normal stress gradient at layer interfaces, so that they are no longer uniform across the thickness. As before, C α k and C ς k restore the continuity of displacements at interfaces of mathematical layers.
HWZZM is developed starting from the displacement field in Equations (29) and (30) in the same way as HWZZ; namely, the decomposition into fictitious computational layers is not allowed for the displacement and in-plane strain fields, while it is allowed for out-of-plane strains. Contributions to in-plane displacements over the third-order and to the transverse displacement over the fourth-order are neglected. Stress–boundary conditions are enforced at the top and bottom laminate faces, while local equilibrium equations are enforced at the inner layers. Membrane stresses σ x x , σ y y , and σ x y come from stress–strain relations, while out-of-plane master stresses are derived integrating local equilibrium equations. The numerical findings will show that although different zig-zag functions contradistinguish ZZA, HWZZ, and HWZZM, their results are indistinguishable, so there will be evidence that the choice of such functions is immaterial, provided that [ ] i and [ ] c are recomputed at each interface.

3.3.2. HWZZM* New Theory and Theories Derived from It

In order to assess the effect of the choice of different zig-zag functions, a new theory called HWZZM* is derived from ZZA*, similar to how HWZZ was derived from ZZA. Then, other theories are particularized from HWZZM, assuming A k u α ( z ) , A k u ς ( z ) , and B k u ς ( z ) are equal to those at a specific interface, which differ from one another.

3.3.2.1. HWZZM* Theory

The displacement field of HWZZM* is assumed in the following form
u α ( x , y , z ) = [ u α 0 ( x , y ) + z ( Γ α 0 ( x , y ) w 0 ( x , y ) , α ) ] 0 + { k = 1 n B ˜ α i k ( x , y ) z + + [ C α i ( x , y ) z 2 ] + [ D α i ( x , y ) z 3 ] } i + c
u ς ( x , y , z ) = [ w 0 ( x , y ) ] 0 + { [ b i ( x , y ) z + k = 1 n b ˜ i k ( x , y ) z ] + [ c i ( x , y ) z 2 ] + + [ d i ( x , y ) z 3 ] + e i ( x , y ) z 4 } i + c
whose amplitudes are computed at each interface from the enforcement of the stress compatibility conditions in Equation (19). In this case, terms [ ] i only serve to satisfy the stress–boundary conditions and local equilibrium equations, while the satisfaction of stress compatibility conditions is demanded to [ ] c .
It should be noted that Equations (31) and (32) imply a little reduction of the processing time, which is equal to 10% per each layer with respect to ZZA and to a somewhat lesser extent equal to 6% with respect to HWZZ. Of course, such an advantage will become more consistent as soon as the number of computational/constituent layers increase.

3.3.2.2. Theories with a Priori Chosen Zig-Zag Amplitudes

Several theories are particularized assuming a priori zig-zag amplitudes as the ones by HWZZM competing at a specific interface and then keeping them unchanged across the thickness, which is something similar to what is done using Murakami’s zig-zag functions, whose amplitudes are assumed a priori.
The amplitudes of HWZZMA are assumed to be coincident with those of HWZZM at the first interface from below, while the A k u ς ( z ) and B k u ς ( z ) of HWZZMB are assumed to be the same as those of HWZZMA at the first interface from below, but instead, A k u α ( z ) is calculated as in HWZZM theory. In HWZZMC, only B k u ς ( z ) is assumed to be uniform across the thickness and coincident with that of HWZZM at the first interface from below, while the remaining are computed at each interface. B k u ς ( z ) are neglected in HWZZM0 theory, while A k u α ( z ) and A k u ς ( z ) are assumed in the same way as HWZZMB.
The HWZZMB2 and HWZZMC2 theories are the same as those of HWZZMB and HWZZMC, respectively, but currently, the amplitudes are assumed to be coincident with those of HWZZM at the first interface from above. Since zig-zag amplitudes are assumed, discontinuous out-of-plane stress may result; hence, the integration of local equilibrium equations is required to obviate this discontinuity, with a corresponding increase in costs by 0.9%. However, it must be considered that zig-zag amplitudes are not being computed at each interface, so a processing time saving of 10% is obtained, and in the end, a positive balance is achieved. However, applications will show a loss of accuracy, making it vain.

4. Lower-Order Theories

Lower-order theories are particularized from ZZA through limiting assumptions that take features that are reminiscent to those of the theories in the literature.

4.1. MHR Theory

MHR considers piecewise cubic in-plane displacements, wherein Murakami’s zig-zag function in Equation (28) is incorporated as the layerwise function, alongside a fourth-order polynomial transverse displacement [47]:
u α ( x , y , z ) = [ u α 0 ( x , y ) + z ( Γ α 0 ( x , y ) w 0 ( x , y ) , α ) ] 0 + [ C α ( x , y ) z 2 + D α ( x , y ) z 3 ] i + u α z ( x , y ) M k ( z ) u ς ( x , y , z ) = [ w 0 ( x , y ) ] 0 + [ a ( x , y ) z + b ( x , y ) z 2 + c ( x , y ) z 3 + d ( x , y ) z 4 ] i
Coefficients C α , D α , a , b , c and d are still calculated by enforcing the fulfilment of the stress boundaries conditions in Equation (16), while u α z is calculated by enforcing the fulfilment of the first and second equilibrium in Equation (18) at the middle plane of the laminate. Since out-of-plane stresses may be still discontinuous at layer interfaces, and since u α z is assumed to be uniform across the thickness, the expressions of out-of-plane stresses are derived integrating local equilibrium equations within the framework of an HR variational theorem (2). A refined version that is referred as MHR±, is obtained by determining the right sign of Murakami’s zig-zag function (28) on a physical basis, instead of being forced to reverse at interfaces by the coefficient ( 1 ) k . The right slope is determined at each interface without a cost burden evaluating what sign ( ± 1 ) k attain the lowest residual force norm from three local equilibrium equations.

4.2. MHR4, MHWZZA and MHWZZA4 Theories

A refined variant MHR4 of MHR is obtained assuming the in-plane displacement field u α by Equation (31), and a fourth-order piecewise variation of the transverse displacement [47]:
u ς ( x , y , z ) = [ w 0 ( x , y ) ] 0 + [ a ( x , y ) z + b ( x , y ) z 2 + c ( x , y ) z 3 + d ( x , y ) z 4 ] i + w z ( x , y ) M k ( z )
Coefficients a to d are determined by enforcing the fulfillment of the stress boundary conditions in Equation (16), whereas w z is calculated by enforcing the fulfilment of the third local equilibrium equation at the middle plane.
MHWZZA is developed assuming the same master displacement field as that of MHR in Equation (33) and the same master strain and stress fields as those of the HWZZ model, in Equations (25) and (26), respectively. To improve the accuracy, the displacement, strain, and stress fields of MHWZZA are recovered using ZZA as the post-processor.
Similarly to MHR±, a refined theory MHR4± is obtained from MHR4 determining the sign of Murakami’s zig-zag functions on a physical basis.
A further theory MHWZZA4 is derived assuming the in-plane displacement field in Equation (31) by MHR, the transverse displacement in Equation (13) by ZZA, and as the master strain and stress fields those by HWZZ in Equations (25) and (26), respectively. The only substantial difference of MHWZZA4 with respect to HWZZ and ZZA is a different zig-zag function; the just-mentioned theories along with ZZA* and HWZZM* make it possible to verify whether the accuracy is sensitive to the choice of zig-zag functions.

4.3. HRZZ and HRZZ4

HRZZ theory is developed postulating a uniform transverse displacement and a third-order zig-zag representation of in-plane displacements [47]:
u α ( x , y , z ) = [ u α 0 ( x , y ) + z ( Γ α 0 ( x , y ) w 0 ( x , y ) , α ) ] 0 + [ C α i ( x , y ) z 2 + D α i ( x , y ) z 3 ] i + [ k = 1 n i Φ α k ( x , y ) ( z z k ) H k ( z ) + k = 1 C u k α ( x , y ) H k ] c u ς ( x , y , z ) = w 0 ( x , y )
Within the framework of the HR theorem, the transverse normal stress σ 33 is assumed to be the same as that of the ZZA model, while the transverse shear stresses σ i 3 are derived from the equilibrium equations assuming kinematic relations in Equation (25) to define membrane stresses. However, currently, second and higher-order derivatives of the d.o.f. are neglected, and a unique computational layer is assumed. Since a uniform transverse displacement is chosen and transformed, reduced stiffness properties are assumed. Then, ε i j σ = C i j k l σ k l , C i j k l is the inverse of E i j k l ; i , j = 1 , 2 x , y ; 3 z ; ε i j u = 1 / 2 ( u i , j + u j , i ) ; γ i j ( i 3 ) ( 33 ) = 2 ε i j ( i 3 ) ( 33 ) , and ε ^ 33 σ = 1 / C 3333 ( σ 33 C 33 i j ε i j u ) .
In order to increase the accuracy of HRZZ, the ZZA theory will be used as the post-processor and the results obtained in this way will be indicated as HRZZ PP in the figures.
HRZZ4 assumes the same in-plane representation of HRZZ, the following fourth-order polynomial approximation of the transverse displacement [47]:
u ς ( x , y , z ) = [ w 0 ( x , y ) ] 0 + [ b ( x , y ) z + c ( x , y ) z 2 + d ( x , y ) z 3 + e ( x , y ) z 4 ] i
and the same stress fields of HRZZ. In this case, ε 33 u being no longer null, it is unnecessary to use transformed, reduced stiffness properties. Similar to in the previous theories, coefficients b to e of Equation (36) are determined by enforcing the stress boundary conditions at the upper + and lower faces (16). The out-of-plane master shear stresses σ i 3 are obtained through integrating local equilibrium equations, while the σ 33 appearing in Equation (2) is assumed to be the same as that of ZZA.

5. Numerical Assessments and Discussion

The accuracy of previous theories is assessed, analysing free vibration modes and the transient dynamic behaviour under impulsive loading of simply supported and clamped, laminated, and sandwich beams and plates.
The lay-up, geometric, and material properties of each case are grouped in Table 3 and Table 4; the trial functions and the expansion order that are used in each study are reported in Table 2; the normalization of quantities is brought in Table 5; while the results are carried in Table 6, Table 7, Table 8, Table 9, Table 10, Table 11, Table 12, Table 13, Table 14, Table 15, Table 16, Table 17, Table 18, Table 19 and Table 20 and the computational cost is reported in Table 21 and Table 22.

5.1. Propped Cantilever Sandwich Plate in Cylindrical Bending under Uniform Static Loading

Firstly, a static analysis of two propped cantilever sandwich plates in cylindrical bending under a uniform load p 0 ( + ) on the top face are considered, in order to preliminary assess the accuracy of the theories in predicting the displacement and stress fields. The structure is clamped at x = 0 and restrained on the lower face at x = Lx; the elastic moduli of faces and core are assumed to be Eu/El = 1.6 and Eu/Ec = 166.6, respectively, and with a Poisson’s ratio υ = 0.3 ; the lower (l) and upper (u) faces are tl = c/2 and tu = c/4 thick, respectively, where c is the core thickness. The results of the theories are compared to the ones by 3D FEA [48], whose elements are formulated in order to fulfil Equations (16) and (18), and the displacements and out-of-plane stresses are assumed as nodal d.o.f. The length-to-thickness ratios Lx/h = 5.714 (case a1) and 20 (case a2) are considered.
Regarding case a1, according to Mattei and Bardella [27] and Icardi and Sola [21], its through-thickness displacement and stress fields are very challenging, and require a very accurate description of the transverse shear force resultant at the clamped (7) and at the simply supported edges (10), and a very accurate description of the transverse normal stress; otherwise, incorrect stress predictions are obtained. Moreover, what makes this a tough case is the opposite sign that is assumed at the upper and lower faces by the transverse shear stress at the supported edge. Figure 1 shows the results for the thickest case (Lx/h = 5.714). It can be seen that, contrary to what was postulated by Murakami’s zig-zag function, at the supported edge, the slope of u ς never reverses, while u α doesn’t reverse near the upper interface, so MHR and MHR4 obtain inaccurate results. Even MHR± and MHR4±, whose interfacial displacement slope is determined on a physical basis, are not adequate, because their kinematics are too poor. MHWZZA4 and MHWZZA, incorporating strain and stress fields from HWZZ, obtain better results than MHR and MHR4, while HRZZ calculates an incorrect null transverse displacement at the supported edge. Summarizing, all of the lower-order theories obtain inaccurate results. Since the results of HWZZM are in good agreement with those of the adaptive theories ZZA, ZZA*, HWZZ, and 3D FEA, it is proven that the choice of layerwise functions is immaterial as long as higher-order coefficients are defined as in Section 3.1 and Section 3.3. HWZZMB, HWZZMB2, and HWZZM0 assuming arbitrarily zig-zag amplitudes can accurately predict the axial displacement but not the transverse one, while HWZZMC and HWZZMC2 calculate precise displacements, whereas HWZZMA always provides wrong results (so they are not reported in the figures). It should be noticed that the behaviour of theories arbitrarily assuming zig-zag amplitudes is very case-dependent, and only some are quite accurate.
The results for the slender case Lx/h = 20 are reported in Figure 2; the calculations show that a still significant difference between the predictions of theories exists, and σ α ς differs from that which is expected for the rather thin sandwich structures. Again, only higher-order theories obtain a good degree of accuracy, and HWZZMA is so inaccurate that results cannot be reported. The obvious conclusion is that even for a length-to-thickness ratio of 20, an accurate description of kinematics is of primary importance in the present case.
On the basis of the processing time required (see Table 21 and Table 22 and Section 5.10), it can be concluded that the adaptive theories ZZA, ZZA*, HWZZ, HWZZM, and HWZZM* are the most efficient ones, as they combine accuracy and low cost. It should be remembered that the same representation with the same low expansion order (Section 2.2) dictated by industrial needs is used for all of the theories. Since the results in the literature for MZZ have shown that greatly increasing the expansion order along with that of the representation across the thickness obviously produces more accurate results, prior statements are valid only for the conditions that are examined here. Below, dynamic tests are carried out that consider progressively more challenging benchmarks that highlight the need for sophisticated theories.

5.2. Free Vibration Modes of Simply Supported Laminated and Sandwich Plates in Cylindrical Bending

First, the analysis of simply supported cross-ply [0/90/0] (case b1) and [0/90]2 (case b2) plates in cylindrical bending is presented, which primarily serves as a preliminary test of the accuracy of 3D FEA [48] in solving dynamics problems, this having yet been tested. Then, a [0/core/0] simply supported sandwich plate in cylindrical bending (cases b3) is studied, which has more marked layerwise effects.
The first fundamental frequency predicted by the present theories for case b1 is given in Table 6 considering a length-to-thickness ratio of 10. Comparisons are carried out with the exact solution and the results of EFSDT and EHSDT theories by Kim [49].
It can be seen that only ZZA, HWZZ, HWZZM, ZZA*, and HWZZM* higher-order theories provide very accurate results, but a sufficient accuracy is obtained by the HRZZ, HRZZ4, MHWZZA, MHWZZA4, MHR, MHR4, MHR±, MHR4±, and HWZZMC simplified theories. Instead, HWZZMA, HWZZMB, HWZZMB2, and HWZZM0 give an inaccurate prediction of it. However, case b1 does not appear to be severe, as many theories prove adequate. Although they have dissimilar characteristics, it is noted that HSDT and FSDT (shear correction factor of 5/6 chosen to minimize the error) calculate even the fundamental frequency inaccurately. The same considerations apply for case b2, with the same length-to-thickness ratio of 10, as shown in Table 6; therefore, the previous considerations are not repeated.
The fundamental frequency for a [0/core/0] simply supported sandwich plate constituting case b3 [49] is reported in Table 7 for three different length-to-thickness ratios (four, 10 and 20). Contrary to what one would expect, errors don’t dramatically decrease as Lx/h increases, and the behaviour of theories remains quite diversified. In particular, MHWZZA, MHWZZA4, and HWZZMA, HSDT, and FSDT (shear correction factor is 5/6) prove to be inaccurate irrespective of the length-to-thickness ratio.
Modal displacements and stresses are reported In Figure 3, Figure 4 and Figure 5, where the results are normalized as shown in Table 5, and show that the axial displacement u α and the in-plane stress σ α α are accurately calculated except by MHWZZA and MHWZZA4, while the transverse shear stress σ α ς is inaccurately calculated by MHR, MHR4, MHR±, MHR4±, MHWZZA, and MHWZZA4 for the thickest case. An even bigger scattering is shown for u ς , which is accurately calculated only by ZZA, HWZZ, HWZZM, ZZA*, and HWZZM*, while on the contrary, the transverse normal stress σ ς ς is accurately calculated by all of the theories. Note that the results by HWZZMA are never reported in Figure 3 and Figure 4 as being too wrong. For the same reason, the stress and displacement distributions are not reported for FSDT and HSDT. Such big errors highlight the inapplicability of the most simplified theories even for slender cases, which is contrary to what is often claimed in the literature.

5.3. Fundamental Frequency and Modal Displacements and Stresses of a Simply Supported, Anti-Symmetric Cross-Ply Plate

The fundamental frequency of [0/90/0/90], simply supported, cross-ply square plates with a length-to-thickness ratio of five [4], is examined in Table 8 for increasing the values of the degree of orthotropy E1/E2 ranging from three to 40 (case c1).
Comparisons are given with the results by a global–local higher-order theory [4] and by Matsunaga [56], which were obtained using a higher-order ESL. It can be seen that all of the theories of this paper provide fairly accurate predictions of the fundamental frequency in correspondence with the lowest orthotropy ratios, while quite dispersed results are shown as this ratio reaches the value of 20. In particular, very incorrect results are given by ESL, FSDT (shear correction factor π 2 / 12 ), and HSDT for E1/E2 values greater than 30 and by theories MHR, MHR4, MHWZZA, MHWZZA4, HWZZM0, HWZZMB2, and HWZZMC2. Adaptive theories ZZA, ZZA*, HWZZ, HWZZM, and HWZZM* instead always provide results that are in very well agreement with the exact solution, irrespective of the orthotropy ratio.
The through-thickness variation of modal displacement u ς , transverse shear σ α ς , and transverse normal σ ς ς modal stresses are reported in Figure 6 and Figure 7 for ratios E1/E2 of 3 and 40, respectively. Note that u α is correctly captured by all of the theories only for E1/E2 = 3, so just the results for this highest ratio are reported. The results by HWZZMB2 and by HWZZM0 are only reported for E1/E2 = 40, as they are inaccurate for all of the other cases. All of the other theories provide a quite accurate prediction of the in-plane displacement as well as of σ α ς , with the only exceptions being MHR and MHR4, while there is a bigger scattering of results regarding σ ς ς .
The fundamental frequencies for case c2 with an orthotropy ratio E1/E2 of 30 and a length-to-thickness ratio of 10/3 are reported in Table 9. Modal displacements and modal in-plane and transverse shear and normal stresses for this case are shown in Figure 8, along with the results by 3D FEA, Zhen and Wanji [4], and Matsunaga [56]. These results show that MHWZZA, MHWZZA4, MHR, and MHR4 overestimate the fundamental frequency, while MHWZZA, MHWZZA4, MHR, and MHR4 incorrectly predict the in-plane displacement and stress, while MHR± and MHR4±, whose slope is defined on a physical basis, obtain better results than their counterparts MHR and MHR4 with the slope assumed a priori.
A rather large dispersion of results is shown for u ς and for the transverse normal stress, while the transverse shear stress is erroneously provided by HRZZ, HRZZ4, MHWZZA, MHWZZA4, MHR, and MHR4 only across the first layer. Again, in this case, with the slope being defined on a physical basis, HSDT and FSDT (the latter uses π2/12 as the shear correction factor in order to improve the accuracy of the results) appear inadequate to perform the analysis.

5.4. Fundamental Frequency of a Cross-Ply Plate with Different Thickness Ratios and Boundary Conditions

The fundamental frequency of a [0/90/0] cross-ply square plate with a thickness ratio ranging from 4 to 100 and simply supported (case d1), all clamped, or simply supported edges and clamped on opposite sides (case d2), are reported in Table 10 and Table 11, respectively. The results are compared to those by Di Sciuva and Icardi [50] (RFSDT, RHSDT, theories and RHQ40 elements) and to 3D FEA [48], as well as to those by FSDT and HSDT.
The results of case d1 show a very well agreement among the theories each other, as well as with 3D FEA reference solutions when the length-to-thickness ratio increases, because the layerwise effects become less important, while remarkable differences are shown for the thickest cases. All of the theories appear to be accurate except for HWZZM0, MHWZZA, and MHWZZA4 (even for moderately thin plates), FSDT (shear correction factor 5/6), and HSDT (the latter two are accurate only for the lengths-to-thickness greater than 10).
The results for the different boundary conditions considered (case d2), namely clamped and simply supported edges (CSCS) and all-clamped edges (CCCC), are reported in Table 11 for a length-to-thickness ratio Lx/h = 10, along with the results for simply supported edges (SSSS).
It is noted that only ZZA, ZZA*, HWZZ, HWZZM and HWZZM* give predictions of the fundamental frequency that are always very well in agreement with 3D FEA, while all of the other theories are less accurate.
It is also noted that the errors are greater for CSCS and CCCC, which therefore prove more problematic than SSSS because the mechanical boundary constraints are more difficult to satisfy identically for clamped edges. In this case, FSDT (shear correction factor 5/6) and HSDT do not provide valid results for any boundary condition. Therefore, it is deduced that for cases d1 and d2 FSDT and HSDT are unsuitable, similarly to other theories with a fixed representation.
So, it is confirmed also in this case that only adaptive theories whose coefficients can be redefined across the thicknesses (ZZA, ZZA*, HWZZ, HWZZM, and HWZZM*) obtain always accurate results with low computational cost (see Section 5.10), and then, they should be preferred in the applications.

5.5. First and Higher-Order Free Vibration Frequencies of Potpourri Cases

Here, the vibration behaviour of plates made of different materials and with different boundary conditions is discussed. Table 12, Table 13 and Table 14 report the first five frequencies for a [0/90/0] square plate (the intermediate layer has a thickness of h/2, while the outer ones have a thickness of h/4) here referred to as case e1, which is retaken from Kapuria et al. [5] (ZIGT and TOT theories). This plate is simply supported along the edges parallel to the y-axis and free at the other two edges (SFSF). Length-to-thickness ratios of 5, 10 and 20 are considered.
Again, adaptive theories ZZA, HWZZ, ZZA*, HWZZM, ZZA*, and HWZZM* appear as the most accurate among those considered, irrespective of the mode examined. On the contrary, the other theories exhibit errors that grow with the frequency order and with the increasing thickness. Only ZZA, HWZZ, ZZA*, HWZZM, ZZA*, HWZZM*, HRZZ, HRZZ4, and MHWZZA4 obtain quite accurate results for length-to-thickness ratios of 5, while MHR, MHR4, MHR±, MHR4±, HWZZMA, and HWZZM0 can’t get the fourth and the fifth frequencies for the intermediate length-to-thickness ratio of 10. Instead, for Lx/h = 20, all of the theories except FSDT (shear correction factor π 2 / 12 ) and HSDT give accurate predictions of frequencies, as the layerwise effects tend to wears off, even if contrary to what is claimed in the literature, this does not always occur in all of the examined cases.
The first eight frequencies for a [0/90/0] square plate with clamped edges (CCCC), which is retaken from Zhen and Wanji [4] and Liew [52], and is here referred to as case e2, are reported in Table 15. In this case, in addition to 3D FEA, comparisons can be carried out with results obtained assuming a linear in-plane displacement across the thickness and a uniform transverse displacement. Errors less than 2% are shown by adaptive theories ZZA, ZZA*, HWZZ, HWZZM, and HWZZM* with respect to 3D FEA, whereas the other theories exhibit larger errors that increase with the order of frequency that is considered. However, guessing the shear correction factor, accurate results can be achieved by FSDT (in the current case with a shear correction factor of 5/6), as also shown in [56], despite the very simple kinematics. This gives a reason at least for the case that is currently examined, to those who consider ESL suitable for dynamic analysis. However, it could be argued that an appropriate shear correction factor could not easily be chosen in the industrial applications.
Table 16 reports the first six frequencies of case e3, which concerns a simply supported (SSSS) [0/90/core/0/90] soft-core sandwich plate with a length-to-thickness ratio of 10. It should be noted that for this case, it is necessary to enforce two additional conditions across the upper layer in the adaptive theories, which means the first two local equilibrium equations; otherwise, the accuracy drops. Having done this, again, all of the higher-order adaptive theories obtain accurate results for all of the frequencies.
Also HRZZ, HRZZ4, HWZZMB, HWZZMC, and HWZZMC2 provide quite accurate results, while other lower-order theories are inaccurate. As in many other cases, MHR and MHR4 are not adequate, because Murakami’s rule is not respected, so MHR± and MHR4±, whose slope sign is determined on a physical basis, obtain better results. Also, in this case, FSDT (shear correction factor 5/6) and HSDT (no need of shear correction factor) obtain inaccurate results.

5.6. Through-Thickness Mode of a Simply Supported, Cross-Ply Plate in Cylindrical Bending

A simply supported cylindrically bent [0/90] plate (case f) is now studied, which is retaken from Pagani et al. [53], who analysed it via finite elements. The interesting aspect of this case is that great cross-section deformations are already involved by the first four frequencies, which could undermine the concept of the plate on which the theories that are considered in this paper are based. Indeed, deformations from the first to the fourth frequency are respectively: a bending mode, a bending/torsional mode, a torsional mode, and finally an axial/shear mode.
The structure is 200-mm thick, while its length-to-thickness and length-to-side ratios are 10. In this case, a suited choice of even and odd trial functions should be made to obtain accurate results by the present theories, and a sufficiently high expansion order of the representation must also be considered, as indicated in Table 2.
The results for this case, which are reported in Table 17, where they are compared to the present FEA results and the ones from [53], show that adaptive ZZA, ZZA*, HWZZ, HWZZM, and HWZZM* theories only commit an error of the order of 2% or less, albeit their representation order across the thickness and their number of unknowns are lower of those of the theory that was used to construct finite elements in [53]. Since there is only one interface, HWZZMA, HWZZMB, HWZZMB2, HWZZMC, and HWZZMC2 achieve the same accuracy of HWZZM, whereas all of the other lower-order theories predict wrong frequencies, which is a sign that a refined kinematics is required. In this case, the results by FSDT and HSDT aren’t reported as being totally wrong.

5.7. First Five Free Vibration Frequencies of a Soft-Core Sandwich with Strong Transverse Normal and Shear Deformability Effects

For the purpose of checking which theory is able to effectively capture strong 3D effects related to transverse normal and shear deformability, case g is examined. It concerns a simply supported cylindrically bent sandwich plate with stiff faces and a compliant core that has not yet been considered in the literature. Each face has a thickness of 0.2 h, and is made of three layers having different thickness and properties, as indicated in Table 3, while the core is 0.6-h thick. The first and the third face layers proceeding from the outside toward the inside are made of the same very stiff material (m2), while the interposed layer is made of a more compliant material in tension, compression, and shear (m1); finally, the core is made of the most compliant material (m3), as indicated in Table 4. A length-to-thickness ratio of five is considered.
Table 18 reports the first five free vibration frequencies for this case, as predicted by the theories of this paper and by 3D FEA. The results indicate that the first, second, and fifth frequencies represent the bending modes, while the third and fourth ones represent more interesting motions that occur across the thickness in a symmetrical manner with respect to the mid-plane. So, specific boundary conditions should be enforced to get rid of these modes.
Moreover, the fulfillment of local equilibrium equations should be enforced near the core interfaces in the adaptive theories, the only ones where it is possible to do it; otherwise, poor results similar to those of other theories are obtained. The results demonstrate the superior accuracy of the adaptive theories ZZA, ZZA*, HWZZA, HWZZM, and HWZZM*, which were obtained thanks to the imposition of these constraints, irrespective of whether they directly considered a piecewise transverse displacement or recovered the effects of the normal transverse deformation differently. Naturally, in this case, the gap with other theories is much more marked than in the previous cases because of the greater importance assumed by the transverse normal deformation, it being the one that was less accurately reproduced. In this case, HWZZMC and HWZZMC2 obtain inaccurate results, with their errors increasing with the frequency order, i.e., the third and fourth frequencies are progressively more inaccurate.
From Figure 9, which reports the in-plane modal displacement, in-plane and transverse shear modal stresses for the first mode, and a minor subject to errors, it can be seen that all of the theories except HWZZMA, MHR, MHR4, MHR±, and MHR4± obtain quite accurate results, which is a sign that Murakami’s theories are inappropriate for this case. However, it is noted that the greatest errors occur on the transverse modal displacement and the normal modal stress. Figure 10 reports the modal stresses and displacements for the third frequency by ZZA, ZZA*, HWZZA, HWZZM, HWZZM*, HWZZMC, and HWZZMC2, while the results by other theories are not reported as being too wrong.

5.8. Free Vibration of a Thick Simply Supported [0/90/Core/0/90] Sandwich Plate

The first six free vibration modes of a simply supported sandwich plate having the same thickness of faces and core of case e3 and the same orientation is studied assuming a thicker length-to-thickness ratio of five. With the same purpose, stiffer faces are considered, as indicated in Table 3, while the core has the same properties as case e3. These choices, which are distinctive of a benchmark that has never been investigated before, here called case h, enhance the layerwise effects, and consequently they should highlight the quite different behaviours of the theories.
Table 19 reports the first six free vibration frequencies predicted by the present theories and by 3D FEA for this new case. Similar to the previous case g, the first five frequencies are bending modes, while the sixth represents a motion that occurs in a symmetrical manner in the thickness direction. So, the same considerations about the constraints that must be imposed apply again, and consequently, it is still demonstrated that the material properties and thickness of constituent layers constitutes a strong discriminatory effect on the accuracy of the theories.
Considerations that are similar to those of the previous case g apply, because only adaptive theories are always in well agreement with 3D FEA results. All of the other theories, except HWZZMC and HWZZMC2, which obtain results accurate enough, give inaccurate predictions, especially for the sixth mode.

5.9. Blast Pulse Loading

In this section, two square plates with a different lay-up and subject to step and exponential blast pulse loadings are analysed. The first is a sandwich plate that is subject to a step pulse (case i1), which has a length Lx of 609.6 mm, its core is 25.4-mm thick, and its two faces, each one being a five-layer laminate, have a total thickness of 1.905 mm. The second is a laminated [0/90/0] square plate whose central layer is two times thicker than the outer ones, whose length is 2540 mm and whose overall thickness is 170 mm (case i2), which is subject to an exponential blast pulse. Such case studies are retaken respectively from Hause and Librescu [54], where the step blast pulse overpressure loading is described as:
p = { p 0 ( + )   i f   t < 5   ms 0   i f   t 5   ms
and from Librescu and Noisier [55], where the exponential blast pulse overpressure is considered as:
p = p 0 ( + ) ( 1 t / t p ) e 2 t / t p
where tp is 0.1 s.
The results for the first case i1 are reported in Figure 11 as the central plate deflection at the mid-plane normalized to the plate thickness. It can be seen that also in this case, only theories ZZA, ZZA*, HRZZ, HRZZ4, and HRZZ4* provide a correct time-variation of the deflection, which is in good agreement with 3D FEA. Notice that the results by [54] differ from the FEM results because of the lower-order model used therein. However, we see a similar behaviour where errors tend to disappear with respect to the present theories if a single half-wave in the x and y directions is considered instead of the expansion order that is reported in Table 2.
As for the other cases considered previously, again, MHR and MHR4 appear to be inadequate because Murakami’s rule is not respected, while their counterparts MHR± and MHR4±, whose slopes are computed on a physical basis, appear to be more accurate.
The results for the case i2 are reported in Figure 12, where again the transverse displacement is reported at z = 0 and at the center of the plate, and is still normalized to the plate thickness. Since now the layer has the same material properties and a symmetric stack-up, the layerwise effects fade, so the discrepancies between the results by theories also fade.
Anyway, MHWZZA4, MHR4, MHR4±, HWZZMA, and HWZZM0 still appear inaccurate, whereas the adaptive theories once again obtain results that are in very good agreement with 3D FEA. Table 20, which reports results for specific instants (0.9 ms, 4.5 ms, 5.6 ms, 6.5 ms, and 7.4 ms for i1 and 3.9 ms, 8.0 ms, 19.96 ms, and 24.08 ms for i2) show that the erroneous predictions (both frequency and amplitude of response are wrong) given by HRZZ, HRZZ4, HWZZMB, HWZZMB2, and HWZZMC2 are erroneous, which had not been noticed for the initial instant (0.9 ms for i1 and 3.9 ms for i2, which corresponded to the occurrence of the first peak of ZZA in each of the two cases) considered in Figure 11 and Figure 12).
It is worth noting that in both cases i1 and i2, there is no detectable difference when the transverse displacement is evaluated in points across the thickness other than at the middle plane as reported in the tables and figures, because both structures are thin. However, splitting the core into two halves whose upper half is much more compliant than the lower one (case i3), similar as to when a face is damaged, and assuming Lx/h = 10 and a different orientation of layers of faces, the results of Figure 13 are obtained, which show visible differences between the results of the theories. Those assuming a uniform or a polynomial transverse displacement in this case don’t account properly for this effect, so less accurate results are obtained, see e.g., MHR±, MHWZZA, MHWZZA4, HRZZ, HRZZ4, FSDT, and HSDT.
It is specified that a Newmark implicit time integration scheme was adopted for solving transient dynamic equations, since alternative explicit time integration schemes need extremely small time steps to be stable. However, for reasons of the stability of the algorithm, small time-steps are still required (30 μ s ) to limit the convergence and rounding errors. Since geometrical and material non-linearity are disregarded, the system to solve is a linear system, and the computational burden isn’t adversely affected by such a small time step.

5.10. Computational Effort of Theories

Table 21 and Table 22 reports the calculation times that are necessary to solve each of the benchmarks considered by the examined theories, which being based on the same five d.o.f., therefore have a memory storage occupation that is practically indistinguishable from one another. As closed-form solutions are considered even when other researchers recoursed to FEA due to the complexity of the solutions in the cases examined, calculation times are very short for all of the theories; that is, they remain comparable to those of FSDT and HSDT.
This testifies to the efficiency of the present adaptive and higher-order, because they require just a reduced expansion order, both with regard to the in-plane and through-thickness representation to achieve accurate results for all of the challenging benchmarks examined. So, it can be said that a level of accuracy that is comparable to that of the FEA has been obtained with a lower computational burden. However, FEA remains indispensable for solving the problems of industrial complexity, while the preliminary parametric studies can be performed as in this paper.
It appears that the MHR and MRH4 theories have the lowest processing time out of all of the theories, but this advantage is totally negated because they provide inaccurate results whenever strong layerwise effects rise, because in these cases, the slope varies differently from what is expected by Murakami’s rule. Although slightly more expensive, MHR± and MHR4± obtain often rather accurate results, since their slope sign of displacements is decided on a physical basis. HRZZ and HRZZ4 result in slower processing times than the adaptive theories for static cases whenever stresses must be computed through time-consuming procedures, but if this is not required, that is only global quantities are required, they result in faster processing times than the adaptive theories. However, this advantage nullified the results of HRZZ and HRZZ4, which were inaccurate in almost all of the cases that were considered.
Higher-order theories HWZZ, HWZZM, ZZA*, and HWZZM* provided rather more accurate results in all of the cases examined, and required a little longer processing time than HRZZ and HRZZ4. However, it is noted that HWZZM* with a priori assumed zig-zag amplitudes requires 20% less processing time than HWZZ, but they don’t appear to be the most accurate theories. In particular, HWZZMA, HWZZMB, HWZZMB2, HWZZMC, and HWZZMC2 appear to be inadequate in many cases. ZZA, HWZZ, HWZZM, ZZM*, and HWZZM* don’t qualify for the lowest calculation time between all of the theories, appear to be the most efficient theories, and thus are preferred in the applications, which are always very accurate and still have affordable costs. However, the best of such adaptive theories from this point of view turns out to be ZZA*, which has a slightly lower calculation cost.

6. Concluding Remarks

Various displacement-based and mixed zig-zag theories, which differ in the layerwise functions and in the scheme of the through-thickness representation of the displacements that are used, have been applied to investigate the free vibration behaviour and the response of blast-loaded laminated and sandwich plates with different length-to-thickness ratios, lay-ups, constituent materials, and boundary conditions. To homogenise the results, they are compared using the same type and order of representation as closed-form solutions, with the appropriate trial functions being selected for each benchmark to minimize the expansion order. The intended aim is to evaluate the merits and drawbacks of theories in order to establish which are significantly much more accurate and efficient.
The numerical applications show the importance of very accurately accounting for the transverse normal deformability whenever the layers have non-uniform mechanical properties and a different thickness. Indeed, adaptive zig-zag theories whose layerwise contributions identically satisfy interfacial stress constrains and whose displacement fields are redefined for each layer prove superiority. ZZA* theory shows that the choice of zig-zag functions is immaterial whenever the coefficients of displacements are recomputed across the computational layers. In this context, zig-zag functions can even be omitted, as the stress continuity constraints can be enforced in order to define the coefficients of displacement fields in a more computationally efficient way.
The accuracy of results is shown to be independent of the choice of zig-zag functions for ZZA*, but this result is extensible to all of the theories that in the same way provide a redefinition of the coefficients of displacements across the thickness, so as to satisfy the physical constraints. Vice versa, the theories whose coefficients of displacements are fixed fail to be accurate whenever strong layerwise effects rise or there is a strong transverse anisotropy, since finding a kind of fixed representation that is always suitable is impossible, unless a very high order of representation is used. That is the opposite of what this paper sets out, which is wanting to find accurate solutions at a low cost. Indeed, the accuracy of theories with a fixed representation appears to be largely case-dependent. Mixed theories such as MHWZZA and MHWZZA4 based on Murakami’s zig-zag function (as well as all those for which zig-zag amplitudes are a priori assumed) are often proven inaccurate, although not in all cases, even though they benefit from strain and stress fields by adaptive theories. The same happens even when the slope sign of displacements at interfaces is established on a physical basis, at least for the low orders of the in-plane and through-thickness representation that are considered in this paper, which however allow the adaptive theories to be already very accurate. Anyhow, it is not easy to discern for which cases the limiting assumptions of such theories do not have weight. Therefore, it is not possible to establish a general rule, although the results undoubtedly show that the theories accounting for layerwise effects without the determination of zig-zag amplitudes on a physical basis cannot provide an adequate level of accuracy with the low expansion orders that are considered in this paper.
A simplified uniform or polynomial representation of the transverse displacement is shown to be ineffective even when the strain and stress fields are retaken from other more accurate structural models, such as for MHWZZA. In particular, FSDT and HSDT theories are proven to be inaccurate in the majority of the examined cases.
Although the adaptive theories whose coefficients of displacements are redefined across the thickness do not get the lowest processing time, they were proven to be the efficient ones, since they always achieve the best accuracy with a processing time that is still short, while the other theories have lower calculation times, but also a much lower accuracy.

Author Contributions

Conceptualization, U.I. and A.U.; Data curation, U.I. and A.U.; Formal analysis, U.I. and A.U.; Methodology, U.I. and A.U.; Project administration, U.I. and A.U.; Software, U.I. and A.U.; Supervision, U.I.; Validation, U.I. and A.U.; Visualization, U.I. and A.U.; Writing—original draft, U.I. and A.U.; Writing—review & editing, U.I. and A.U.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Carrera, E. A study of transverse normal stress effects on vibration of multilayered plates and shells. J. Sound Vib. 1999, 225, 803–829. [Google Scholar] [CrossRef]
  2. Carrera, E.; Ciuffreda, A. Bending of composites and sandwich plates subjected to localized lateral loadings: A comparison of various theories. Compos. Struct. 2005, 68, 185–202. [Google Scholar] [CrossRef]
  3. Demasi, L. Refined multilayered plate elements based on Murakami zig-zag functions. Compos. Struct. 2004, 70, 308–316. [Google Scholar] [CrossRef]
  4. Zhen, W.; Wanji, C. Free vibration of laminated composite and sandwich plates using global–local higher-order theory. J. Sound Vib. 2006, 298, 333–349. [Google Scholar] [CrossRef]
  5. Kapuria, S.; Dumir, P.C.; Jain, N.K. Assessment of zig-zag theory for static loading, buckling, free and forced response of composite and sandwich beams. Compos. Struct. 2004, 64, 317–327. [Google Scholar] [CrossRef]
  6. Burlayenko, V.N.; Altenbach, H.; Sadowski, T. An evaluation of displacement-based finite element models used for free vibration analysis of homogeneous and composite plates. J. Sound Vib. 2015, 358, 152–175. [Google Scholar] [CrossRef]
  7. Li, J.; Hu, X.; Li, X. Free vibration analyses of axially loaded laminated composite beams using a unified higher-order shear deformation theory and dynamic stiffness method. Compos. Struct. 2016, 158, 308–322. [Google Scholar] [CrossRef]
  8. Boscolo, M.; Banerjee, J.R. Layer-wise dynamic stiffness solution for free vibration analysis of laminated composite plates. J. Sound Vib. 2014, 333, 200–227. [Google Scholar] [CrossRef]
  9. Khdeir, A.A.; Aldraihem, O.J. Free vibration of sandwich beams with soft core. Compos. Struct. 2016, 154, 179–189. [Google Scholar] [CrossRef]
  10. Sayyad, A.S.; Ghugal, Y.M. On the free vibration analysis of laminated composite and sandwich plates: A review of recent literature with some numerical results. Compos. Struct. 2015, 129, 177–201. [Google Scholar] [CrossRef]
  11. Kazancı, Z. A review on the response of blast loaded laminated composite plates. Prog. Aerosp. Sci. 2016, 81, 49–59. [Google Scholar] [CrossRef]
  12. Lin, T.R.; Zhang, K. An analytical study of the free and forced vibration response of a ribbed plate with free boundary conditions. J. Sound Vib. 2018, 422, 15–33. [Google Scholar] [CrossRef]
  13. Vescovini, R.; Dozio, L.; D’Ottavio, M.; Polit, O. On the application of the Ritz method to free vibration and buckling analysis of highly anisotropic plates. Compos. Struct. 2018, 192, 460–474. [Google Scholar] [CrossRef]
  14. Rekatsinas, C.S.; Nastos, C.V.; Theodosiou, T.C.; Saravanos, D.A. A time-domain high-order spectral finite element for the simulation of symmetric and anti-symmetric guided waves in laminated composite strips. Wave Mot. 2015, 53, 1–19. [Google Scholar] [CrossRef]
  15. Valisetty, R.R.; Rehfield, L.W. Application of ply level analysis to flexural wave vibration. J. Sound Vib. 1988, 126, 183–194. [Google Scholar] [CrossRef]
  16. Cho, K.N.; Bert, C.W.; Striz, A.G. Free vibration of laminated rectangular plates analyzed by higher-order individual-layer theory. J. Sound Vib. 1991, 145, 429–442. [Google Scholar] [CrossRef]
  17. Noiser, A.; Kapania, R.K.; Reddy, J.N. Free vibration analysis of laminated plates using a layer-wise theory. AIAA J. 1993, 31, 2335–2346. [Google Scholar] [CrossRef]
  18. Wu, C.P.; Chen, W.Y. Vibration and stability of laminated plates based on a local higher-order plate theory. Sound Vib. 1994, 177, 503–520. [Google Scholar] [CrossRef]
  19. Carrera, E. Layerwise mixed models for accurate vibration analysis of multilayered plates. Int. J. Appl. Mech. 1998, 65, 820–828. [Google Scholar] [CrossRef]
  20. Rao, M.K.; Desai, Y.M.; Chitnis, M.R. Free vibrations of laminated beams using mixed theory. Compos. Struct. 2001, 52, 149–160. [Google Scholar] [CrossRef] [Green Version]
  21. Icardi, U.; Sola, F. Assessment of recent zig-zag theories for laminated and sandwich structures. Compos. Part B 2016, 97, 26–52. [Google Scholar] [CrossRef]
  22. Di Sciuva, M. A refinement of the transverse shear deformation theory for multilayered orthotropic plates. In Proceedings of the AIDAA National Congress, Napoli, Italy, 25–28 October 1983; pp. 84–92. [Google Scholar]
  23. Murakami, H. Laminated composite plate theory with improved in-plane responses. ASME Appl. Mech. 1986, 53, 661–666. [Google Scholar] [CrossRef]
  24. Carrera, E. An assessment of mixed and classical theories on global and local response of multilayered orthotropic plates. Compos. Struct. 2000, 50, 183–198. [Google Scholar] [CrossRef]
  25. Carrera, E. On the use of the Murakami’s zig-zag function in the modeling of layered plates and shells. Compos. Struct. 2004, 82, 541–554. [Google Scholar] [CrossRef]
  26. Brischetto, S.; Carrera, E.; Demasi, L. Improved response of asymmetrically laminated sandwich plates by using Zig-Zag functions. J. Sand. Struct. Mat. 2009, 11, 257–267. [Google Scholar] [CrossRef]
  27. Mattei, O.; Bardella, L. A structural model for plane sandwich beams including transverse core deformability and arbitrary boundary conditions. Eur. J. Mech. Part A Solids 2016, 58, 172–186. [Google Scholar] [CrossRef]
  28. Icardi, U. Higher-order zig-zag model for analysis of thick composite beams with inclusion of transverse normal stress and sublaminates approximations. Compos. Part B 2001, 32, 343–354. [Google Scholar] [CrossRef]
  29. Icardi, U.; Sola, F. Development of an efficient zig-zag model with variable representation of displacements across the thickness. J. Eng. Mech. 2014, 140, 531–541. [Google Scholar] [CrossRef]
  30. Li, X.Y.; Liu, D. Generalized laminate theories based on double superposition hypothesis. Int. J. Numer. Meth. Eng. 1997, 40, 197–212. [Google Scholar] [CrossRef]
  31. Zhen, W.; Wanji, C. An efficient higher-order theory and finite element for laminated plates subjected to thermal loading. Compos. Struct. 2006, 73, 99–109. [Google Scholar] [CrossRef]
  32. Zhen, W.; Wanji, C. A study of global–local higher-order theories for laminated composite plates. Compos. Struct. 2007, 79, 44–54. [Google Scholar] [CrossRef]
  33. Shariyat, M. A generalized global-local high-order theory for bending and vibration analyses of sandwich plates subjected to thermo-mecanical loads. Compos. Struct. 2010, 92, 130–143. [Google Scholar] [CrossRef]
  34. Surana, K.; Nguyen, S. Two-dimensional curved beam element with higher order hierarchical transverse approximation for laminated composites. Compos. Struct. 1990, 36, 499–511. [Google Scholar] [CrossRef]
  35. Karama, M.; Afaq, K.; Mistou, S. Mechanical behaviour of laminated composite beam by the new multi-layered laminated composite structures model with transverse shear stress continuity. Int. J. Solids Struct. 2003, 40, 1525–1546. [Google Scholar] [CrossRef]
  36. Jun, L.; Hongxing, H. Dynamic stiffness analysis of laminated composite beams using trigonometric shear deformation theory. Compos. Struct. 2009, 89, 433–442. [Google Scholar] [CrossRef]
  37. Vidal, P.; Polit, O. Assessment of the refined sinus model for the non-linear analysis of composite beams. Compos. Struct. 2009, 87, 370–381. [Google Scholar] [CrossRef] [Green Version]
  38. Rodrigues, J.D.; Roque, C.M.C.; Ferreira, A.J.M.; Carrera, E.; Cinefra, M. Radial basis functions-finite differences collocation and a unified formulation for bending, vibration and buckling analysis of laminated plates, according to Murakami’s zig-zag theory. Compos. Struct. 2011, 93, 1613–1620. [Google Scholar] [CrossRef]
  39. Mantari, J.; Oktem, A.; Soares, C.G. A new higher order shear deformation theory for sandwich and composite laminated plates. Compos. Part B 2012, 43, 1489–1499. [Google Scholar] [CrossRef]
  40. Barut, A.; Madenci, E.; Tessler, A. A refined zigzag theory for laminated composite and sandwich plates incorporating thickness stretch deformation. In Proceedings of the 53rd AIAA/ASME/ASCE/AHS/ACS Structures, Structural Dynamics and Materials Conference, Honolulu, HI, USA, 23–26 April 2012; p. 1705. [Google Scholar]
  41. Kim, J.S.; Cho, M. Enhanced first-order theory based on mixed formulation and transverse normal effect. Int. J. Solids Struct. 2007, 44, 1256–1276. [Google Scholar] [CrossRef]
  42. Tessler, A.; Di Sciuva, M.; Gherlone, M. A refined zigzag beam theory for composite and sandwich beams. J. Compos. Mat. 2009, 43, 1051–1081. [Google Scholar] [CrossRef]
  43. Zhen, W.; Wanji, C. A global higher-order zig-zag model in terms of the HW variational theorem for multi-layered composite beams. Compos. Struct. 2016, 158, 128–136. [Google Scholar] [CrossRef]
  44. Gherlone, M. On the use of zigzag functions in equivalent single layer theories for laminated composite and sandwich beams: A comparative study and some observations on external weak layers. ASME Appl. Mech. 2013, 80, 1–19. [Google Scholar] [CrossRef]
  45. Groh, R.M.J.; Weaver, P.M. On displacement-based and mixed-variational equivalent single layer theories for modeling highly heterogeneous laminated beams. Int. J. Solids Struct. 2015, 59, 147–170. [Google Scholar] [CrossRef]
  46. Li, R.; Wang, P.; Zheng, X.; Wang, B. New benchmark solutions for free vibration of clamped rectangular thick plates and their variants. Appl. Math. Lett. 2018, 78, 88–94. [Google Scholar] [CrossRef]
  47. Icardi, U.; Urraci, A. Novel HW mixed zig-zag theory accounting for transverse normal deformability and lower-order counterparts assessed by old and new elastostatic benchmarks. Aerosp. Sci. Technol. 2018, 80, 541–571. [Google Scholar] [CrossRef]
  48. Icardi, U.; Atzori, A. Simple, efficient mixed solid element for accurate analysis of local effects in laminated and sandwich composites. Adv. Eng. Soft. 2004, 35, 843–859. [Google Scholar] [CrossRef]
  49. Kim, J.S. Free vibration of laminated and sandwich plates using enhanced plate theories. J. Sound Vib. 2007, 308, 268–286. [Google Scholar] [CrossRef]
  50. Di Sciuva, M.; Icardi, U. Numerical studies on bending, free vibration and buckling of multilayered anisotropic plates. L’Aerot. Miss. Spaz. 1993, 72, 1–2. [Google Scholar]
  51. Icardi, U.; Urraci, A. Dipartimento di Ingegneria Meccanica e Aerospaziale, Politecnico di Torino, Torino, Italy. Unpublished work. 2018. [Google Scholar]
  52. Liew, K.M. Solving the vibration of thick symmetric laminates by Reissner/Mindlin plate theory and the P-Ritz method. J. Sound Vib. 1996, 198, 343–360. [Google Scholar] [CrossRef]
  53. Pagani, A.; Carrera, E.; Boscolo, M.; Banerjee, J.R. Refined dynamic stiffness elements applied to free vibration analysis of generally laminated composite beams with arbitrary boundary conditions. Compos. Struct. 2014, 110, 305–316. [Google Scholar] [CrossRef]
  54. Hause, T.; Librescu, L. Dynamic response of anisotropic sandwich flat panels to explosive pressure pulses. Int. J. Imp. Eng. 2005, 31, 607–628. [Google Scholar] [CrossRef]
  55. Librescu, L.; Nosier, A. Response of laminated composite flat panels to sonic boom and explosive blast loadings. AIAA J. 1990, 28, 345–352. [Google Scholar] [CrossRef]
  56. Matsunaga, H. Vibration and stability of cross-ply laminated composite plates according to a global higher-order plate theory. Compos. Struct. 2000, 48, 231–244. [Google Scholar] [CrossRef]
Figure 1. Case a1: Propped cantilever sandwich plate in cylindrical bending (Lx/h = 5.714) under a uniform loading on the top layer.
Figure 1. Case a1: Propped cantilever sandwich plate in cylindrical bending (Lx/h = 5.714) under a uniform loading on the top layer.
Aerospace 05 00108 g001
Figure 2. Case a2: Propped cantilever sandwich plate in cylindrical bending (Lx/h = 20) under a uniform loading on the top layer.
Figure 2. Case a2: Propped cantilever sandwich plate in cylindrical bending (Lx/h = 20) under a uniform loading on the top layer.
Aerospace 05 00108 g002
Figure 3. Case b3: Normalized modal displacements and stresses of a simply supported sandwich plate in cylindrical bending (Lx/h = 4).
Figure 3. Case b3: Normalized modal displacements and stresses of a simply supported sandwich plate in cylindrical bending (Lx/h = 4).
Aerospace 05 00108 g003
Figure 4. Case b3: Normalized modal displacements and stresses of a simply supported sandwich plate in cylindrical bending (Lx/h = 10).
Figure 4. Case b3: Normalized modal displacements and stresses of a simply supported sandwich plate in cylindrical bending (Lx/h = 10).
Aerospace 05 00108 g004
Figure 5. Case b3: Normalized modal displacements and stresses of a simply supported sandwich plate in cylindrical bending (Lx/h = 20).
Figure 5. Case b3: Normalized modal displacements and stresses of a simply supported sandwich plate in cylindrical bending (Lx/h = 20).
Aerospace 05 00108 g005aAerospace 05 00108 g005b
Figure 6. Case c1: Normalized modal displacements and stresses of a simply supported laminated plate (Lx/h = 5, E1/E2 = 3).
Figure 6. Case c1: Normalized modal displacements and stresses of a simply supported laminated plate (Lx/h = 5, E1/E2 = 3).
Aerospace 05 00108 g006
Figure 7. Case c1: Normalized modal displacements and stresses of a simply supported laminated plate (Lx/h = 5, E1/E2 = 40).
Figure 7. Case c1: Normalized modal displacements and stresses of a simply supported laminated plate (Lx/h = 5, E1/E2 = 40).
Aerospace 05 00108 g007aAerospace 05 00108 g007b
Figure 8. Case c2: Normalized modal displacements and stresses of a simply supported laminated plate (Lx/h = 10/3).
Figure 8. Case c2: Normalized modal displacements and stresses of a simply supported laminated plate (Lx/h = 10/3).
Aerospace 05 00108 g008
Figure 9. Case g: Normalized modal displacements and stresses of a simply supported sandwich plate in cylindrical bending (Lx/h = 5, first mode).
Figure 9. Case g: Normalized modal displacements and stresses of a simply supported sandwich plate in cylindrical bending (Lx/h = 5, first mode).
Aerospace 05 00108 g009aAerospace 05 00108 g009b
Figure 10. Case g: Normalized modal displacements and stresses of a simply supported sandwich plate in cylindrical bending (Lx/h = 5, third mode).
Figure 10. Case g: Normalized modal displacements and stresses of a simply supported sandwich plate in cylindrical bending (Lx/h = 5, third mode).
Aerospace 05 00108 g010
Figure 11. Case i1: Normalized transverse displacement of a simply supported sandwich plate under a step blast pulse loading (Lx/h = 20.8696).
Figure 11. Case i1: Normalized transverse displacement of a simply supported sandwich plate under a step blast pulse loading (Lx/h = 20.8696).
Aerospace 05 00108 g011
Figure 12. Case i2: Normalized transverse displacement of a simply supported laminated plate under an exponential blast pulse loading (Lx/h = 14.941).
Figure 12. Case i2: Normalized transverse displacement of a simply supported laminated plate under an exponential blast pulse loading (Lx/h = 14.941).
Aerospace 05 00108 g012
Figure 13. Case i3: Normalized transverse displacement of a simply supported sandwich plate under a step blast pulse loading (Lx/h = 10).
Figure 13. Case i3: Normalized transverse displacement of a simply supported sandwich plate under a step blast pulse loading (Lx/h = 10).
Aerospace 05 00108 g013
Table 1. Acronyms.
Table 1. Acronyms.
AcronymMeaningAcronymMeaning
3D FEAMixed solid 3D elements (Ref. [48]).MHRHR mixed theory with Murakami’s zig-zag function (see Section 4.1, Ref. [47]).
DLDiscrete-layer theories.MHR±MHR with slope defined on a physical basis (see Section 4.1).
DZZDi Sciuva’s zig-zag theory (Ref. [22]).MHR4MHR with fourth-order polynomial transverse displacement (see Section 4.2, Ref. [47]).
EFSDTTheory from Ref. [49].MHR4±MHR4 with slope defined on a physical basis (see Section 4.2).
EHSDTTheory from Ref. [49].MHWZZAModified HWZZ theory, type A (see Section 4.2, Ref. [47]).
ESLEquivalent single-layer.MHWZZA4Modified MHWZZA theory, with fourth-order piecewise polynomial transverse displacement (see Section 4.2, Ref. [47]).
FSDTFirst-order shear deformation theory.MZZMurakami’s like zig-zag theory (Ref. [23]).
HRHellinger–Reissner variational theorem.PVWPrinciple of virtual work.
HRZZHR zig-zag theory (see Section 4.3, Ref. [47]).RFSDTTheory from Ref. [50].
HRZZ PPPost-processed HRZZ (see Section 4.3).RHQ40Finite element from Ref. [50].
HRZZ4HRZZ with fourth-order polynomial transverse displacement (see Section 4.3, Ref. [47]).RHSDTTheory from Ref. [50].
HSDTHigher-order shear deformation theory.RZTTheory from Ref. [42].
HWHu–Washizu variational theorem.SEUPTStrain energy update technique (Ref. [29]).
HWgHu–Washizu canonical functional (see Section 2.1).TOTTheory from Ref. [5].
HWZZHW zig-zag mixed theory (see Section 3.2, Ref. [47]).ZIGTTheory from Ref. [5].
HWZZMModified HWZZ theory, type M, (see Section 3.3.1).ZZPlate theories with zig-zag contributions.
HWZZM*Modified HWZZ theory, type M*, (see Section 3.3.2).ZZAZig-zag adaptive theory (see Section 3.1, Ref. [29]).
HWZZM(†)Modified HWZZM theories, type † (see Section 3.3.2.2).ZZA*Modified ZZA theory (see Section 3.1.1).
† = A, B, B2, C, C2, 0 (variants of theory).
Table 2. Trial functions.
Table 2. Trial functions.
CaseExpansion OrderMesh (xa·yb·zh) (+)Trial Function
b1116·2·60 u 0 ( x , y ) = m = 1 M A m cos ( m π x L x ) ; w 0 ( x , y ) = m = 1 M C m sin ( m π x L x ) ; γ x 0 ( x , y ) = m = 1 M D m cos ( m π x L x ) ;
b2116·2·60
b3116·2·60
e1516·2·60
g316·2·60
c1110·10·28 u 0 ( x , y ) = m = 1 M n = 1 N A m n cos ( m π L x x ) sin ( n π L y y ) ; v 0 ( x , y ) = m = 1 M n = 1 N B m n sin ( m π L x x ) cos ( n π L y y ) ;
w 0 ( x , y ) = m = 1 M n = 1 N C m n sin ( m π L x x ) sin ( n π L y y ) ;
γ x 0 ( x , y ) = m = 1 M n = 1 N D m n cos ( m π L x x ) sin ( n π L y y ) ; γ y 0 ( x , y ) = m = 1 M n = 1 N E m n sin ( m π L x x ) cos ( n π L y y ) ;
c2110·10·28
d1110·10·28
d2110·10·28
e3610·10·28
h510·10·28
i11110·10·28
i21110·10·28
i31110·10·28
a1916·2·60 u 0 ( x , y ) = i = 1 I A i ( x L ) i ; w 0 ( x , y ) = i = 1 I C i ( x L ) i ; γ x 0 ( x , y ) = i = 1 I D i ( x L ) i
a2916·2·60
d2410·10·28 u α 0 ( x , y ) = j = 1 J i = 1 I A α i ( x L x ) i ( y L y ) j ; γ α 0 ( x , y ) = j = 1 J i = 1 I D α i ( x L x ) i ( y L y ) j ;
e21010·10·28
f610·10·28
(+) A uniform mesh is used; xa and yb represent the number of elements in x and y directions, respectively, zh is the number of elements across the thickness.
Table 3. Data of cases.
Table 3. Data of cases.
CaseLay-UpLayer ThicknessMaterialBCSLx/hLy/Lx
a1 [27][0/0/0][(2h/7)/(4h/7)/(h/7)][n/n/n]CS5.714-
a2[0/0/0][(2h/7)/(4h/7)/(h/7)][n/n/n]CS20-
b1 [49][0/90/0][h/3]3[p]3SS10-
b2 [49][0/90/0/90][0.25h]4[p]4SS10-
b3 [49] ()[0/0/0][0.1h/0.8h/0.1h][p/mc/p]SS4, 10, 20-
c1 [4][0/90/0/90][0.25h]4[a]4SSSS51
c2 [4][90/0/90/0][0.25h]4[b]4SSSS10/31
d1 [5][0/90/0][h/3]3[c]3SSSS4, 10, 20, 30, 50, 1001
d2 [50][0/90/0][h/3]3[c]3SSSS, CCCC, CSCS101
e1 [5][0/90/0][0.25h/0.5h/0.25h][d]3SFSF5, 10, 20-
e2 [52][0/90/0][h/3]3[e]3CCCC101
e3 [4][0/90/0/0/90][(h/24)2/(5h/12)]S[f2/g]SSSSS101
f [53][0/90][h/2]2[h]2SS100.1
g (*)[0]8[0.025h/0.05h/0.125h/0.3h]S[m2/m1/m2/m3]SSS5-
h (*)[0/90/0/0/90][(h/24)2/(5h/12)]S[l1/l2/g/l1/l2]SSSS51
i1 [54][(45/−45)2/45/0]S[(0.381mm)5/(12.7mm)]S[o15/o2]SSSSS20.86961
i2 [55][0/90/0][(h/4)/(h/2)/(h/4)][q]3SSSS14.9411
i3 (*†)[(0/90)2/02]S[(0.381mm)5/(12.7mm)]S[o15/o2/o3/o15]SSSS101
* Transverse anisotropy; Strong layerwise effects; (BCS: boundary conditions name; SS simply supported, C clamped).
Table 4. Mechanical Properties.
Table 4. Mechanical Properties.
Material Nameabcdefghl1l2m1m2m3mcn [iso]o1o2o3pq
E1 [GPa]M130E225E218140E21316.89 × 10−325E233.51391330.050.1-206.840.1380.0138172.4132.4
E2 [GPa]---10.3-10.346.89 × 10−3-83.475110.050.1-5.1710.1380.01386.910.8
E3 [GPa]E2E2E210.3E210.346.89 × 10−3E283.475110.020.1M25.1710. 1380. 01386.910.8
G12 [GPa]0.6E20.6E20.5E27.170.6E26.2053.45 × 10−30.5E22.261.73750.0280.02170.04-2.5510.10270.010273.455.6
G13 [GPa]0.6E20.6E20.5E27.170.6E26.8953.45 × 10−30.2E22.261.73750.0280.02170.04-2.5510.10270.010273.455.6
G23 [GPa]0.5E20.5E20.2E22.870.5E26.8953.45 × 10−30.2E230.6950.0280.02170.04-2.5510.062050.0062051.385.6
υ120.250.250.250.250.250.2200.250.350.250.250.250.150.250.330.250.350.350.250.24
υ130.250.250.250.250.250.2200.250.350.250.250.250.150.250.330.250.350.350.250.24
υ230.250.250.250.330.250.4900.250.330.250.250.250.150.250.330.250.020.020.250.24
Density [kg/m3]ρρρ1587ρ162797ρ162716271558.351558.3516.3136ρρ1558.3516.313616.3136ρ1443
M1 E1/E2 = 3, 10, 20, 30, 40; M2 Eu/El = 1.6, Eu/Ec = 166.66; [iso] = isotropic E1 = E2 = E3 G12 = G13 = G23.
Table 5. Normalization of displacements, stresses and frequencies.
Table 5. Normalization of displacements, stresses and frequencies.
CaseNormalization
a1 u α ¯ = u ( L x , z ) h p 0   u ς ¯ = u ς ( L x , z ) h p 0   σ α α ¯ = σ α α ( L x , z ) p 0 ( L x / h ) 2   σ α ς ¯ = A σ α ς ( L x , z ) P 0 L x   σ ς ς ¯ = σ ς ς ( L x , z ) p 0
a2 u α ¯ = u ( L x , z ) h p 0   u ς ¯ = u ς ( L x , z ) h p 0   σ α α ¯ = σ α α ( L x , z ) p 0 ( L x / h ) 2   σ α ς ¯ = A σ α ς ( L x , z ) P 0 L x   σ ς ς ¯ = σ ς ς ( L x , z ) p 0
b1 b2 b3 ω ¯ = ω h ρ _ M A T p G 12 _ M A T p   u i ¯ = u i | u i | max   σ i j ¯ = σ i j | σ i j | max
c1 ω ¯ = 10 ω h ρ _ M A T a E 2 _ M A T a   u i ¯ = u i | u i | max   σ i j ¯ = σ i j | σ i j | max
c2 ω ¯ = 10 ω h ρ _ M A T b E 2 _ M A T b   u i ¯ = u i | u i | max   σ i j ¯ = σ i j | σ i j | max
d1 d2 ω ¯ = ω L x h L x 2   ρ _ M A T c E 2 _ M A T c
e1 ω ¯ = ω L x 2 h ρ _ M A T d E 2 _ M A T d
e2 ω ¯ = ω L y 2 π 2 ρ _ M A T e h   D 0   D 0 = E 2 _ M A T e h 3 12 ( 1 ν 21 _ M A T e ν 12 _ M A T e )   ν 21 _ M A T e = 0.00625
e3 ω ¯ = ω L x 2 h ρ _ M A T f E 2 _ M A T f
f ω ¯ = ω L x 2 L y ρ _ M A T h E 2 _ M A T h
g ω ¯ = ω L x 2 h ρ _ M A T m 2 E 2 _ M A T m 2   u i ¯ = u i | u i | max   σ i j ¯ = σ i j | σ i j | max
h ω ¯ = ω L x 2 h ρ _ M A T l 2 E 2 _ M A T l 2
i1 i2 u ς ¯ = u ς h
i3 u ς ¯ = u ς w , w   static   response
Table 6. Normalized fundamental frequencies, cases b1 and b2.
Table 6. Normalized fundamental frequencies, cases b1 and b2.
Theory0/90/0 (b1)0/90/0/90 (b2)Theory0/90/0 (b1)0/90/0/90 (b2)
3D [49]0.14620.1095HWZZMC0.14630.1095
EFSDT [49]0.14480.1125HWZZMC20.14650.1095
EHSDT [49]0.14600.1090HWZZM00.15500.1220
HSDT [49]0.15050.1139HRZZ0.14620.1094
FSDT [49]0.16160.1186HRZZ40.14620.1095
3D FEA [48]0.14640.1096MHWZZA0.14600.1095
ZZA0.14630.1095MHWZZA40.14610.1094
ZZA*0.14630.1095MHR0.14630.1115
HWZZ0.14630.1095MHR40.14630.1162
HWZZM0.14630.1095MHR±0.14630.1115
HWZZM*0.14630.1095MHR4±0.14630.1162
HWZZMA0.24990.1130HSDT0.15110.1157
HWZZMB0.21990.1095FSDT0.15650.1162
HWZZMB20.15180.1124
Table 7. Normalized fundamental frequencies, case b3.
Table 7. Normalized fundamental frequencies, case b3.
Lx/h41020Lx/h41020
3D [49]0.10110.03430.0155HWZZMC0.10170.03430.0156
EFSDT [49]0.08480.03320.0153HWZZMC20.10170.03430.0155
EHSDT [49]0.09720.03410.0155HWZZM00.10160.03430.0155
FSDT [49]0.33250.10050.0316HRZZ0.09920.03420.0155
3D FEA [48]0.10110.03430.0155HRZZ40.09950.03420.0155
ZZA0.10150.03430.0155MHWZZA0.11070.03730.0170
ZZA*0.10150.03430.0155MHWZZA40.10830.03740.0170
HWZZ0.10150.03430.0155MHR0.10150.03430.0155
HWZZM0.10150.03430.0155MHR40.10150.03430.0155
HWZZM*0.10150.03430.0155MHR±0.10150.03430.0155
HWZZMA0.12500.03570.0159MHR4±0.10150.03430.0155
HWZZMB0.10190.03430.0155HSDT0.13840.05030.0213
HWZZMB20.10160.03430.0155FSDT0.30550.09560.0310
Table 8. Normalized fundamental frequencies, case c1.
Table 8. Normalized fundamental frequencies, case c1.
TheoryE1/E2 = 310203040
Exact [4]2.61823.25783.76224.06604.2719
Zhen and Wanji [4]2.61273.25133.75234.05324.2568
Matsunaga [56]2.60213.23803.74004.04254.2477
3D FEA [48]2.61183.24923.75004.05094.2546
ZZA2.60263.23873.74064.04304.2481
ZZA*2.60263.23873.74064.04304.2481
HWZZ2.60263.23873.74064.04304.2481
HWZZM2.60263.23873.74064.04304.2481
HWZZM*2.60263.23873.74064.04304.2481
HWZZMA2.61583.25913.76644.07164.2783
HWZZMB2.61213.25293.75724.06034.2653
HWZZMB23.11574.39805.72656.71867.4660
HWZZMC2.60643.24683.75374.06014.2681
HWZZMC22.62773.30043.85414.19704.4252
HWZZM02.88903.45503.91884.20474.4024
HRZZ2.59773.22893.72504.02304.2247
HRZZ42.59803.22933.72554.02384.2256
MHWZZA2.08444.04683.30203.13472.9820
MHWZZA42.06394.00523.26263.11242.9785
MHR2.61013.36964.12434.69405.1627
MHR42.62653.39704.16164.73895.2142
MHR±2.60383.29153.91654.34404.6645
MHR4±2.63123.37744.01874.42544.7118
HSDT2.60033.27813.85054.21394.4686
FSDT2.59863.28363.86514.23424.4919
Table 9. Normalized fundamental frequencies (NFQ), case c2.
Table 9. Normalized fundamental frequencies (NFQ), case c2.
c2NFQTheoryNFQTheoryNFQ
3D FEA [48]0.7041HWZZMB20.7268MHWZZA40.7652
ZZA0.7044HWZZMC0.7043MHR0.9394
ZZA*0.7044HWZZMC20.7199MHR40.7596
HWZZ0.7044HWZZM00.7047MHR±0.7609
HWZZM0.7044HRZZ0.6860MHR4±0.7674
HWZZM*0.7044HRZZ40.6855HSDT0.7432
HWZZMA0.7094MHWZZA0.7505FSDT0.7451
HWZZMB0.7084
Table 10. Normalized fundamental frequencies, case d1.
Table 10. Normalized fundamental frequencies, case d1.
Lx/h410203050100
RFSDT [50]7.163011.700014.012014.645015.009015.1720
RHSDT [50]7.027011.462013.888014.579014.983015.1650
3D FEA [48]6.843611.430613.923114.659315.069615.2550
ZZA6.925411.458313.888914.580014.983315.1654
ZZA*6.925411.458313.888914.580014.983315.1654
HWZZ6.925411.458313.888914.580014.983315.1654
HWZZM6.925411.458313.888914.580014.983315.1654
HWZZM*6.925411.458313.888914.580014.983315.1654
HWZZMA9.059714.281415.081515.037115.074815.1730
HWZZMB8.300715.940318.080215.843315.155415.1760
HWZZMB27.625112.683613.924814.587114.984315.1655
HWZZMC6.999911.466013.889714.580214.983415.1654
HWZZMC27.140911.482313.891514.580614.983415.1654
HWZZM07.285912.535617.732621.723627.429833.3973
HRZZ6.389611.450213.886214.578714.982815.1653
HRZZ46.910411.456913.888814.579914.983315.1654
MHWZZA9.96329.100012.309913.895214.738515.1043
MHWZZA49.98159.105412.310013.895214.738515.1043
MHR6.991311.464713.889614.580214.983415.1654
MHR47.472511.664414.096114.743715.071215.1924
MHR±6.991311.464713.889614.580214.983415.1654
MHR4±7.472511.664414.096114.743715.071215.1924
HSDT7.116011.790014.060014.670015.019015.1750
FSDT7.413012.163014.230014.757015.053015.1830
Table 11. Normalized fundamental frequencies, case d2.
Table 11. Normalized fundamental frequencies, case d2.
TheorySSSSCSCSCCCCTheorySSSSCSCSCCCC
RHQ40 [50]11.400015.955017.4810HWZZM012.535617.142618.2727
3D FEA [48]11.430615.389516.6658HRZZ11.450216.130417.5659
ZZA11.458315.187516.4575HRZZ411.456916.148117.6134
ZZA*11.458315.187516.4575MHWZZA9.10006.031320.8865
HWZZ11.458315.187516.4575MHWZZA49.10546.336821.0304
HWZZM11.458315.187516.4575MHR11.464716.614518.0505
HWZZM*11.458315.187516.4575MHR411.664416.971318.4802
HWZZMA14.281417.439519.2037MHR±11.464716.614518.0505
HWZZMB15.940320.803120.9520MHR4±11.664416.971318.4802
HWZZMB212.683616.840917.8968HSDT11.790017.415718.5237
HWZZMC11.466014.566614.7679FSDT12.163016.443617.5603
HWZZMC211.482316.503318.0240
Table 12. Normalized fundamental frequencies, case e1.
Table 12. Normalized fundamental frequencies, case e1.
Lx/h = 5Mode 1Mode 2Mode 3Mode 4Mode 5
Exact [5]6.806016.515026.688037.255048.0350
Zhen and Wanji [4]6.816116.615427.421239.823054.1984
ZIGT [5]6.812816.729727.862340.906056.1529
TOT [5]6.976216.878327.141738.186450.3407
FSDT [5]7.445818.249128.769741.166849.2839
3D FEA [48]6.811516.534326.747537.374648.3752
ZZA6.811616.525326.712337.280149.1224
ZZA*6.811616.525326.712337.280149.1224
HWZZ6.811616.525326.712337.280149.1224
HWZZM6.811616.525326.712337.280149.1224
HWZZM*6.811616.525326.712337.280149.1224
HWZZMA6.835117.543527.226039.087655.4790
HWZZMB6.813616.550426.815037.566649.9111
HWZZMB26.816516.587626.968037.974750.6078
HWZZMC6.815316.579526.944637.923250.5189
HWZZMC26.815016.576926.940037.921750.5219
HWZZM06.813616.690727.127438.019650.2967
HRZZ6.814016.502626.606237.032248.7296
HRZZ46.808016.487026.605737.109148.9903
MHWZZA6.417515.573225.202135.871146.9424
MHWZZA46.827116.567326.819937.470249.4131
MHR6.819816.646327.206638.518051.6381
MHR46.820216.661427.270638.649751.8742
MHR±6.819816.646327.206638.518051.6381
MHR4±6.820216.661427.270638.649751.8742
HSDT6.982816.882827.148538.171351.2504
FSDT7.423718.156528.588838.785650.1428
Table 13. Normalized fundamental frequencies, case e1.
Table 13. Normalized fundamental frequencies, case e1.
Lx/h = 10Mode 1Mode 2Mode 3Mode 4Mode 5
Exact [5]9.343427.224046.419066.058086.1690
Zhen and Wanji [4]9.374827.326746.660966.689887.7990
ZIGT [5]9.343427.251246.651166.916888.4094
TOT [5]9.474227.904647.625967.511387.7200
FSDT [5]9.763929.783151.478772.994194.1827
3D FEA [48]9.353427.196246.474766.302486.8487
ZZA9.355627.246346.449266.054487.6999
ZZA*9.355627.246346.449266.054487.6999
HWZZ9.355627.246346.449266.054487.6999
HWZZM9.355627.246346.449266.054487.6999
HWZZM*9.355627.246346.449266.054487.6999
HWZZMA9.363427.340346.950170.120289.1944
HWZZMB9.356127.254446.485666.155088.2045
HWZZMB29.356827.266146.538966.304288.5122
HWZZMC9.356127.261246.523966.271888.4563
HWZZMC29.356027.259846.518966.261588.4411
HWZZM09.356227.436246.835466.738289.0772
HRZZ9.358927.256046.434665.963587.7726
HRZZ49.354927.232046.386665.901287.7148
MHWZZA7.498720.791842.663865.006887.5899
MHWZZA49.373327.306746.548166.221488.3184
MHR9.357727.279246.613466.536489.1480
MHR49.357727.280846.629366.596289.3332
MHR±9.357727.279246.613466.536489.1480
MHR4±9.357727.280846.629366.596289.3332
HSDT9.483227.931247.636867.482289.5987
FSDT9.760829.695051.264872.572395.7078
Table 14. Normalized fundamental frequencies, case e1.
Table 14. Normalized fundamental frequencies, case e1.
Lx/h = 20Mode 1Mode 2Mode 3Mode 4Mode 5
Exact [5]10.640037.374071.7440108.8900147.0400
Zhen and Wanji [4]10.689737.681072.6193110.5376149.6296
ZIGT [5]10.640037.374071.7440108.9989147.4811
TOT [5]10.693237.897273.2506111.6123150.8630
FSDT [5]10.789039.055877.0531119.1257162.4792
3D FEA [48]10.652737.375571.5020108.9389147.0909
ZZA10.657537.422271.8164108.9263149.6777
ZZA*10.657537.422271.8164108.9263149.6777
HWZZ10.657537.422271.8164108.9263149.6777
HWZZM10.657537.422271.8164108.9263149.6777
HWZZM*10.657537.422271.8164108.9263149.6777
HWZZMA10.660137.453671.9498109.2996150.7192
HWZZMB10.657737.424371.8267108.9586149.7543
HWZZMB210.658037.427271.8418109.0057149.8559
HWZZMC10.657537.424671.8334108.9862149.8207
HWZZMC210.657537.424271.8312108.9804149.8098
HWZZM010.658337.635572.2589109.6109150.6415
HRZZ10.658937.435771.8491108.9650149.2023
HRZZ410.657437.419771.7993108.8690149.2281
MHWZZA10.565237.120271.2611108.0743146.5888
MHWZZA410.671937.495271.9809109.1660148.0695
MHR10.658137.430671.8579109.0565150.0877
MHR410.658137.430771.8588109.0624150.1300
MHR±10.658137.430671.8579109.0565150.0877
MHR4±10.658137.430771.8588109.0624150.1300
HSDT10.706137.932773.3318111.6654153.6012
FSDT10.804739.043376.9153118.7217164.5034
Table 15. Normalized fundamental frequencies, case e2.
Table 15. Normalized fundamental frequencies, case e2.
TheoryMode 1Mode 2Mode 3Mode 4Mode 5Mode 6Mode 7Mode 8
Liew [54]4.44706.64207.70009.18509.738011.399011.644012.4660
Zhen and Wanji [4]4.54006.52408.17809.47309.492011.769012.395012.9040
3D FEA [48]4.48156.46377.99389.30259.350711.372211.532912.3221
ZZA4.46826.55618.10659.32719.337111.175611.757412.1556
ZZA*4.46826.55618.10659.32719.337111.175611.757412.1556
HWZZ4.46826.55618.10659.32719.337111.175611.757412.1556
HWZZM4.46826.55618.10659.32719.337111.175611.757412.1556
HWZZM*4.46826.55618.10659.32719.337111.175611.757412.1556
HWZZMA5.85207.52088.172410.489411.172211.256912.035313.1691
HWZZMB4.95706.92028.19599.612910.702811.167511.980013.0014
HWZZMB26.08087.76358.873910.049210.894911.211412.685113.2650
HWZZMC4.48596.62438.16069.64889.853811.188411.834013.2442
HWZZMC24.50016.66688.09399.646410.348911.183411.828513.6387
HWZZM04.43806.46928.62689.635610.215711.167011.985812.6024
HRZZ4.47806.57568.65939.07359.525810.981911.633512.1995
HRZZ44.49236.56228.61209.44189.441810.911411.641012.2389
MHWZZA5.770412.557512.829214.256815.368018.804921.744128.1598
MHWZZA45.891711.122712.912614.115715.215818.618721.528827.8810
MHR4.62366.49969.02789.56439.738012.289312.338114.5644
MHR45.757810.832812.832814.395614.415619.156619.943920.5045
MHR±4.62366.49969.02789.56439.738012.289312.338114.5644
MHR4±5.757810.832812.832814.395614.415619.156619.943920.5045
HSDT4.24246.48519.08189.608410.351511.588014.393113.3191
FSDT4.18166.59798.51479.317110.294811.377011.530613.1945
Table 16. Normalized fundamental frequencies, case e3.
Table 16. Normalized fundamental frequencies, case e3.
TheoryMode 1Mode 2Mode 3Mode 4Mode 5Mode 6
Zhen and Wanji [4]1.94453.37964.59145.52686.51288.4311
3D FEA [48]1.84973.23084.30915.28106.15837.7796
ZZA1.85663.23414.32155.27786.17827.8258
ZZA*1.85663.23414.32155.27786.17827.8258
HWZZ1.85663.23414.32155.27786.17827.8258
HWZZM1.85663.23414.32155.27786.17827.8258
HWZZM*1.85663.23414.32155.27786.17827.8258
HWZZMA2.99458.246413.002023.189327.644829.3959
HWZZMB1.86993.22384.47806.46647.09738.4051
HWZZMB22.55943.39024.47656.74777.117219.0397
HWZZMC4.32683.25724.36695.32036.23117.8900
HWZZMC21.86033.32974.47685.44476.33948.0418
HWZZM03.01337.492010.391011.321515.134621.7804
HRZZ1.83593.15464.13424.94885.63266.5177
HRZZ41.83753.16284.15414.98325.69316.6864
MHWZZA4.59916.39562.56022.985955.003235.4337
MHWZZA44.55366.33232.53492.956354.458635.0828
MHR13.945629.503941.610968.742665.861771.6472
MHR413.945629.503941.610968.742665.861771.6472
MHR±1.85423.24534.32815.26746.15647.7888
MHR4±11.069014.113014.972315.274416.532418.4374
HSDT6.579712.796814.566122.932326.754230.9263
FSDT13.975231.013642.094651.720959.247772.5151
Table 17. Normalized fundamental frequencies, case f.
Table 17. Normalized fundamental frequencies, case f.
TheoryMode 1Mode 2Mode 3Mode 4
3D FEA [53]5.77208.722017.821095.3190
N = 3 (30 DOFS) [53]5.79768.842819.538098.1640
3D FEA [48]5.77698.791417.870595.3828
ZZA5.77708.842417.796297.1215
ZZA*5.77708.842417.796297.1215
HWZZ5.77708.842417.796297.1215
HWZZM5.77708.842417.796297.1215
HWZZM*5.77708.842417.796297.1215
HWZZMA5.77708.842417.796297.1215
HWZZMB5.77708.842417.796297.1215
HWZZMB25.77708.842417.796297.1215
HWZZMC5.77708.842417.796297.1215
HWZZMC25.77708.842417.796297.1215
HWZZM05.77918.852217.7800101.5429
HRZZ5.77228.838918.4970101.6298
HRZZ45,77678,808718,4986101,6402
MHWZZA6.11019.833522.217399.8945
MHWZZA45.80639.833022.207199.2816
MHR5.87669.649220.563997.0202
MHR45.89559.649220.2339105.1698
MHR±5.87669.649220.563997.0202
MHR4±5.89559.649220.2339105.1698
Table 18. Normalized fundamental frequencies, case g.
Table 18. Normalized fundamental frequencies, case g.
TheoryMode 1Mode 2Mode 3Mode 4Mode 5
3D FEA [48]4.178510.423814.658716.497919.0105
ZZA4.178810.426314.686416.527219.0239
ZZA*4.178810.426314.686416.527219.0239
HWZZ4.178810.426314.686416.527219.0239
HWZZM4.178810.426314.686416.527219.0239
HWZZM*4.178810.426314.686416.527219.0239
HWZZMA14.666351.829582.2952171.820889.6849
HWZZMB4.235033.590991.3964157.627968.3833
HWZZMB24.379024.830691.5540159.309158.0426
HWZZMC4.179010.429214.685716.527519.1303
HWZZMC24.179010.428814.686016.527619.0365
HWZZM04.230011.514091.2793157.433941.1479
HRZZ4.13739.920688.7309151.579216.6322
HRZZ44.13809.926688.6647151.253116.6531
MHWZZA3.07765.009920.313129.15206.3628
MHWZZA44.197711.334421.334324.626121.2708
MHR15.205042.468277.9566154.708270.0325
MHR415.201142.287265.227565.402766.8219
MHR±15.205042.468277.9566154.708270.0325
MHR4±15.062032.367817.925426.715757.4925
HSDT14.189237.7022******61.7475
FSDT16.388649.9499******87.0203
*** Not provided by the theory.
Table 19. Normalized fundamental frequencies, case h.
Table 19. Normalized fundamental frequencies, case h.
TheoryMode 1Mode 2Mode 3Mode 4Mode 5Mode 6
3D FEA [48]1.68822.87963.47234.30334.68995.7441
ZZA1.68982.88553.47774.31714.70305.7500
ZZA*1.68982.88553.47774.31714.70305.7500
HWZZ1.68982.88553.47774.31714.70305.7500
HWZZM1.68982.88553.47774.31714.70305.7500
HWZZM*1.68982.88553.47774.31714.70305.7500
HWZZMA5.31417.211011.333319.326421.711746.2461
HWZZMB1.70102.89133.56054.36934.710134.1888
HWZZMB21.76514.72103.54735.057810.654434.7002
HWZZMC1.68982.88553.47784.31714.70305.7500
HWZZMC21.68982.88553.47784.31714.70305.7500
HWZZM02.51695.63645.90198.882711.310634.3120
HRZZ1.68232.85173.39404.16484.590734.3046
HRZZ41.68212.85253.39654.17204.594834.1832
MHWZZA11.76542.71532.72643.75266.87371.4635
MHWZZA41.17763.93254.31654.39504.56565.6519
MHR12.714715.138016.428827.162627.600964.6322
MHR412.762616.612122.268927.777127.868775.2673
MHR±1.69592.90973.49194.34054.764361.7387
MHR4±5.15105.83566.77047.26187.267266.5689
HSDT16.561028.720637.728344.030144.2319***
FSDT11.078317.636120.978425.061925.3697***
*** Not provided by the theory.
Table 20. Normalized fundamental frequencies, cases i1 and i2.
Table 20. Normalized fundamental frequencies, cases i1 and i2.
Casei1i2
t [s]0.00090.00450.00560.00650.00740.00390.00800.019960.02408
3D FEA [48]2.02682.0165−1.32911.3278−1.31740.4064-0.04670.3316−0.1147
ZZA2.03992.0179−1.32981.3286−1.31810.4079−0.04680.3319−0.1147
ZZA*2.03992.0179−1.32981.3286−1.31810.4079−0.04680.3319−0.1147
HWZZ2.03992.0179−1.32981.3286−1.31810.4079−0.04680.3319−0.1147
HWZZM2.03992.0179−1.32981.3286−1.31810.4079−0.04680.3319−0.1147
HWZZM*2.03992.0179−1.32981.3286−1.31810.4079−0.04680.3319−0.1147
HWZZMA1.52261.2548−0.25190.2502−0.18670.4774−0.09630.3901−0.1635
HWZZMB0.86070.8663−0.57650.5915−0.59140.3870−0.04200.3162−0.0798
HWZZMB20.00000.0001−0.0001−0.0001−0.00010.3097−0.03570.2522−0.0870
HWZZMC0.67180.5765−0.46770.1563−0.22020.4080−0.04680.3320−0.1147
HWZZMC21.07410.9654−0.59350.9727−1.08300.4064−0.04660.3308−0.1142
HWZZM00.85620.8598−0.57460.5855−0.58220.29120.0232−0.04430.2659
HRZZ2.02191.9176−1.63411.4543−1.51710.4079−0.04680.3320−0.1147
HRZZ42.02121.9175−1.63321.4540−1.51690.4076−0.04680.3317−0.1146
MHWZZA1.65441.1687−0.09940.1028−0.03310.4080−0.04680.3320−0.1147
MHWZZA41.78471.7135−0.86180.8970−0.83440.4600−0.03620.24850.0311
MHR0.77740.9828−0.4954−0.07650.74640.4077−0.04660.3319−0.1146
MHR40.73201.0386−0.2750−0.22670.56040.37370.00190.11970.1580
MHR±1.97341.9659−1.25951.2568−1.25860.4077−0.04660.3319−0.1146
MHR4±1.25270.0008−0.62260.8942−1.07960.37370.00190.11970.1580
Table 21. Processing time[s].
Table 21. Processing time[s].
Casea1a2b1b2b3c1c2d1d2
ZZA15.067115.97195.38666.87905.119429.699230.673526.445727.3202
HRZZ18.231218.22614.61176.18884.898827.797728.437024.201925.0022
HRZZ418.223718.48915.01385.03875.030234.108736.442626.352627.2240
HWZZ12.427112.84904.96406.27614.867927.359128.573924.528625.3396
MHR6.95746.62582.81074.82882.791822.108723.642917.371217.9456
MHR46.49466.97022.90935.14522.685323.059924.098717.849318.4395
MHWZZA7.23597.69523.76065.26133.664025.695925.696020.530421.2093
MHWZZA47.83657.58613.76025.26083.663625.701225.841220.555321.2350
HWZZM11.534411.70594.18875.59544.001427.236827.160422.083122.8133
HWZZMA11.526511.60184.15955.40613.940126.792226.416121.638122.3536
HWZZMB11.530711.62894.18175.52163.981926.534926.869221.824222.5459
HWZZMC11.531411.64574.18695.59263.919826.560526.895121.828822.5506
HWZZMB211.531011.63894.16595.54903.990526.579726.914621.854422.5770
HWZZMC211.531711.64014.18495.59513.899626.590526.925521.818422.5399
HWZZM011.428711.59124.15545.50793.799426.261726.341121.459622.1692
MHR±6.95746.62583.03884.87702.819722.108723.642917.464818.0423
MHR4±6.49466.97023.03845.19672.712223.059924.098717.945118.5385
ZZA*11.495111.61253.87225.17223.837825.330225.259220.758121.6727
HWZZM*10.957711.00353.93745.31563.801324.063724.512319.510420.3723
FSDT--3.03973.81512.61008.76248.896811.709212.0963
HSDT--3.25074.18392.613411.560811.676413.181113.6169
Processing time of 3D FEA is always about 12 times longer than that of ZZA, excluding meshing preparation.
Table 22. Processing time[s].
Table 22. Processing time[s].
Casee1e2e3fghi1i2i3
ZZA15.214649.899852.378820.991620.441557.4363147.685976.1909143.1814
HRZZ13.923745.666057.033419.210618.707252.5631135.155569.7264130.4250
HRZZ415.161049.724153.085720.917720.369657.2341147.166175.9227150.2565
HWZZ14.111646.282337.595419.469818.959653.2725136.979570.6674138.1438
MHR9.993932.777238.630113.788613.427337.727797.009250.046893.0318
MHR410.268933.679344.432714.168113.796838.766099.679151.4242100.9251
MHWZZA11.811438.738344.486516.296315.869244.5891114.651959.1487117.1169
MHWZZA411.825738.785247.793116.316015.888444.6430114.790659.2202116.9716
HWZZM12.704741.668046.830117.528717.069447.9613123.322863.6220119.6847
HWZZMA12.448740.828447.232817.175516.725446.9949120.837962.3400117.8773
HWZZMB12.555841.179547.242717.323216.869347.3990121.877062.8761122.6692
HWZZMC12.558441.188147.298117.326916.872847.4089121.902562.8892121.4148
HWZZMB212.573141.236547.220217.347216.892647.4645122.045562.9630125.2797
HWZZMC212.552441.168646.443617.318616.864847.3864121.844662.8594125.6217
HWZZM012.346040.491537.797917.033816.587446.6071119.840761.8256114.0284
MHR±10.047732.953838.837613.862913.499637.930997.531750.316493.58166
MHR4±10.324133.860247.793114.244213.870938.9742100.214551.700498.46074
ZZA*12.069538.751344.488616.477016.215944.6040114.690259.8047112.7978
HWZZM*11.224637.209641.819315.486614.763342.3786106.661956.2164106.5141
FSDT6.736422.093825.3415-9.050825.430765.390033.734564.1452
HSDT7.583224.871128.5271-10.188528.627573.609937.975275.5128
Processing time reported include symbolic computation. For calculation, a laptop computer with quad-core [email protected] GHz. 64-bit operating system and 8.00 GB RAM was used.

Share and Cite

MDPI and ACS Style

Icardi, U.; Urraci, A. Free and Forced Vibration of Laminated and Sandwich Plates by Zig-Zag Theories Differently Accounting for Transverse Shear and Normal Deformability. Aerospace 2018, 5, 108. https://doi.org/10.3390/aerospace5040108

AMA Style

Icardi U, Urraci A. Free and Forced Vibration of Laminated and Sandwich Plates by Zig-Zag Theories Differently Accounting for Transverse Shear and Normal Deformability. Aerospace. 2018; 5(4):108. https://doi.org/10.3390/aerospace5040108

Chicago/Turabian Style

Icardi, Ugo, and Andrea Urraci. 2018. "Free and Forced Vibration of Laminated and Sandwich Plates by Zig-Zag Theories Differently Accounting for Transverse Shear and Normal Deformability" Aerospace 5, no. 4: 108. https://doi.org/10.3390/aerospace5040108

APA Style

Icardi, U., & Urraci, A. (2018). Free and Forced Vibration of Laminated and Sandwich Plates by Zig-Zag Theories Differently Accounting for Transverse Shear and Normal Deformability. Aerospace, 5(4), 108. https://doi.org/10.3390/aerospace5040108

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop