Next Article in Journal
Analysis of the Temperature Distribution on the Surface of Saddle-Shaped Briquettes Consolidated in the Roller Press
Next Article in Special Issue
Tie Rod-Equivalent Non-Linear Constitutive Law for Uniformly Loaded Cables
Previous Article in Journal
Reduction of Load Capacity of Fiber Cement Board Facade Cladding under the Influence of Fire
Previous Article in Special Issue
Shear Wave Splitting and Polarization in Anisotropic Fluid-Infiltrating Porous Media: A Numerical Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Third-Order Theory for the Bending Analysis of Laminated Thin and Thick Plates Including the Strain Gradient Effect

by
Michele Bacciocchi
1,2,* and
Angelo Marcello Tarantino
2,3
1
Dipartimento di Economia, Scienze e Diritto (DESD), University of San Marino, Via Consiglio dei Sessanta, 47891 Dogana, San Marino
2
Centro di Ricerca Interdipartimentale Costruzioni e del Territorio (CRICT), Via Vivarelli, 41125 Modena, Italy
3
Department of Engineering “Enzo Ferrari” (DIEF), University of Modena and Reggio Emilia, Via Vivarelli, 41125 Modena, Italy
*
Author to whom correspondence should be addressed.
Materials 2021, 14(7), 1771; https://doi.org/10.3390/ma14071771
Submission received: 6 March 2021 / Revised: 26 March 2021 / Accepted: 31 March 2021 / Published: 3 April 2021

Abstract

:
The aim of the paper is the development of a third-order theory for laminated composite plates that is able to accurately investigate their bending behavior in terms of displacements and stresses. The starting point is given by the corresponding Reddy’s Third-order Shear Deformation Theory (TSDT). This model is then generalized to consider simultaneously the Classical Laminated Plate Theory (CLPT), as well as the First-order Shear Deformation Theory (FSDT). The constitutive laws are modified according to the principles of the nonlocal strain gradient approach. The fundamental equations are solved analytically by means of the Navier methodology taking into account cross-ply and angle-ply lamination schemes. The numerical applications are presented to highlight the nonlocal effects on static behavior.

1. Introduction

Higher-order plates theories for laminates have been introduced in the last decades to avoid some issues related to the use of lower-order and simpler approaches, such as the Classical Laminated Plate Theory (CLPT) and First-order Shear Deformation Theory (FSDT) [1,2]. In particular, higher-order equivalent single-layer theories allow to obtain a more accurate interlaminar stress analysis and do not need to introduce the shear correction factor [3,4,5,6]. The displacement field that characterizes the Third-order Shear Deformation Theory (TSDT), for instance, determines a quadratic profile of shear strains and stresses along the thickness [7,8], due to its cubic expansion in the thickness coordinate. Consequently, there is no need for the shear correction factor [9,10,11,12]. The importance of these cubic terms in the analysis of laminates has been recently highlighted in the paper by Petrolo and Carrera [13], in which the best theory diagrams for multilayered structures have been widely discussed.
In general, higher-order approaches have been justified by the use of more and more advanced materials [14,15,16,17] and the need of innovative configurations for the optimal design of structures [18,19]. In particular, their introduction could be essential when these innovative constituents are included in the stacking sequences of multilayered or sandwich structures [20,21,22,23,24,25]. These aspects have been clearly emphasized in the works by Carrera [26,27,28], Carrera and Giunta [29], and Carrera et al. [30]. In these works, accurate and effective higher-order structural models based on a unified formulation have been presented.
The increasing number of applications involving micro- and nanostructures [31,32,33,34] has proven that the size-dependent features of the advanced constituents could have not negligible effects on their mechanical behavior [35,36,37]. These aspects have been highlighted by so-called multiscale analysis [38,39,40,41]. In these circumstances, structural theories based on classical elasticity could turn out to be inadequate to model such innovative mediums [42,43,44,45,46,47]. Nonlocal theories have been developed to overcome these issues, as illustrated in the first papers by Eringen [48,49]. Some other subsequent contributions could be mentioned to this aim. For instance, Luciano and Willis investigated the nonlocal constitutive behavior of an infinite composite laminate [50]. Nanorods made of functionally graded materials were studied by Barretta et al. [51], who included the gradient Eringen model in their framework. The same approach was used for Timoshenko nanobeams in bending [52]. Analogously, the effects of nonlocal elasticity were emphasized in the papers by Apuzzo et al. [53,54] for the investigation of torsional behavior and dynamic response of Bernoulli–Euler nanobeams. More recently, similar remarks and observations have been drawn in [55,56,57]. In general, different nonlocal theoretical frameworks can be developed. Examples of nonlocal approaches that should be mentioned for completeness purposes are the strain and stress gradient theories [58,59,60], the modified couple stress theory [61,62,63,64], and the ones based on micropolar formulations [65,66,67]. A comprehensive literature review concerning nonlocal elasticity can be found in the paper by Zhao et al. [68] for completeness purposes.
The current paper is developed within the strain gradient theoretical framework [69]. According to this approach, the size effects of the materials are modeled by means of a length-scale parameter, which controls the second-order derivatives of the strain components [70]. Consequently, the stresses are functions not only of the strains in an evaluation point, but also depend on the divergence of the gradient of the strains. These aspects have been emphasized in the papers by Aifantis [71] and by Askes and Aifantis [72], where an overview of gradient elasticity formulations in statics and dynamics has been presented. Therefore, higher-order derivatives of displacements are involved in the three-dimensional constitutive laws [73,74,75,76]. These aspects have been also recently emphasized in papers [77,78] where peculiar Finite Element (FE) formulations have been developed to take into account the strain gradient effect. In particular, the use of higher-order Hermite interpolating polynomials for the approximation of both membrane and bending degrees of freedom have been disclosed.
The plate theory that this nonlocal effect is included in is based on the TSDT for laminates [1], in which the various layers that define the stacking sequence have orthotropic features [79,80]. The kinematic model is written to take into account simultaneously lower-order approaches, such as the CLPT and FSDT, for comparison purposes. After the brief introduction of the main topics of the paper (Section 1), the theory is presented by using a matrix compact notation in Section 2. Here, the strong form of the governing equations is obtained, once the definitions of both strains and stress resultants are carried out. The fundamental system is solved analytically by means of the Navier approach as shown in [73]. The algebraic system of equations is written in Section 3, highlighting the contributions related to the strain gradient effect for cross-ply and angle-ply simply supported laminated composite plates. As remarked in [1], the Navier methodology can be efficiently applied to deal with these configurations. Section 4 is focused on the numerical applications. The proposed approach is validated for both classical and nonlocal elasticity through comparison with the results available in the literature, taking into account the three different structural models. Then, the results are extended to emphasize the influence of the nonlocal parameter on the static behavior, which is expressed in terms of displacements and stress components. The through-the-thickness stress profiles are also provided. Finally, the main achievements are drawn in Section 5. Appendix A collects instead the analytical expressions of the terms involved in the algebraic formulation for cross-ply and angle-ply laminates.

2. Nonlocal Structural Model

The theoretical framework is developed in this Section for a rectangular plate whose planar size is a × b , considering a Cartesian reference system x y z . The plate is made of a sequence of N L orthotropic layers. The generic k-th ply has a thickness h k = z k + 1 z k , with z k + 1 , z k being its upper and lower thickness coordinates. The overall plate thickness is given by h = k = 1 N L h k . These geometric features are all shown and specified in Figure 1.
The three-dimensional displacements collected in U x , y , z = U V W T can be written in terms of the five degrees of freedom, which are three translations u , v , w and two rotations ϕ x , ϕ y . The vector u x , y = u v w ϕ x ϕ y T is conveniently introduced. The displacement field assumes the following compact aspect:
U = I ( 0 ) u + z I ( 1 ) u c 1 F I ( 1 ) u + D ( 0 ) I ( 3 ) u ,
where D ( 0 ) is a differential operator given by
D ( 0 ) = 0 0 x 0 0 y 0 0 0 ,
whereas the matrices I ( 0 ) , I ( 1 ) , I ( 3 ) are defined as
I ( 0 ) = 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 , I ( 1 ) = 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 , I ( 3 ) = 0 0 1 0 0 0 0 1 0 0 0 0 1 0 0 .
The choice of the structural theory defines the values of the constant c 1 and the thickness function F z . In particular, the TSDT is obtained for c 1 = 4 3 h 2 and F = z 3 . By setting c 1 = 0 , the FSDT is achieved instead. On the other hand, the CLPT can be defined by using c 1 = 1 and F = z . The fundamental assumptions of each theory, therefore, are different according to this choice [1]. It should be specified that only the FSDT requires the shear correction factor in the definitions of the shear stresses. The value of 5/6 is considered to this aim.
The membrane strain components ε = ε x x ε x y γ x y T and the transverse shear strains γ = γ x z γ y z T can be written as follows:
ε = ε ( 0 ) + z ε ( 1 ) + F ε ( 3 ) , γ = γ ( 0 ) + F γ ( 2 ) ,
where F = F z . The terms introduced in (4) are discussed below. According to the notation employed by Reddy [1], the membrane strains ε ( 0 ) assume the following aspect:
ε ( 0 ) = ε x x ( 0 ) ε y y ( 0 ) γ x y ( 0 ) T = D ( m ) I ( 0 ) u .
On the other hand, the curvatures ε ( 1 ) are computed as follows:
ε ( 1 ) = ε x x ( 1 ) ε y y ( 1 ) γ x y ( 1 ) T = D ( m ) I ( 1 ) u ,
where the differential operator D ( m ) is given by
D ( m ) = x 0 0 0 y 0 y x 0 .
Higher-order membrane strains vector ε ( 3 ) is defined as
ε ( 3 ) = ε x x ( 3 ) ε y y ( 3 ) γ x y ( 3 ) T = c 1 D ( m ) I ( 1 ) + D ( b ) I ( 3 ) u ,
in which the following definition is used for the differential operator D ( b ) :
D ( b ) = 0 0 2 x 2 0 0 2 y 2 0 0 2 2 x y .
The shear strains are now discussed. The components collected in γ ( 0 ) can be written as
γ ( 0 ) = γ x z ( 0 ) γ y z ( 0 ) T = I s ( 1 ) u + D s ( 0 ) I s ( 3 ) u ,
where
I s ( 1 ) = 0 0 0 1 0 0 0 0 0 1 , I s ( 3 ) = 0 0 1 0 0 0 0 1 0 0 ,
whereas the differential operator D s ( 0 ) is given by
D s ( 0 ) = 0 x 0 y .
Likewise, higher-order shear terms included in γ ( 2 ) are defined as
γ ( 2 ) = γ x z ( 2 ) γ y z ( 2 ) T = c 1 I s ( 1 ) u + D s ( 0 ) I s ( 3 ) u .
The governing equations are derived by means of the principle of virtual displacements [1]
δ Φ + δ L = 0 ,
in which δ Φ is the strain energy variation, whereas δ L represents the work done by applied external forces. If a laminated composite plate made of N L layers is considered, the variation δ Φ can be defined as
δ Φ = k = 1 N L A z k z k + 1 δ ε T σ ( k ) + δ γ T τ ( k ) d z d A ,
in which A denotes the plate middle surface. The membrane stresses in the k-th layer are specified by σ ( k ) = σ x x ( k ) σ y y ( k ) σ x y ( k ) T and assume the following aspect [77,78]:
σ ( k ) = 1 2 2 Q ¯ ( k ) ε ,
where is the nonlocal parameter linked to the influence of the micro/macroscale interactions, 2 = 2 x 2 + 2 y 2 is the Laplacian, and Q ¯ ( k ) is the plane-stress-reduced stiffness coefficients matrix given by
Q ¯ ( k ) = Q ¯ 11 Q ¯ 12 Q ¯ 16 Q ¯ 12 Q ¯ 22 Q ¯ 26 Q ¯ 16 Q ¯ 26 Q ¯ 66 ( k ) .
Its terms are computed as a function of the Young’s moduli E 1 , E 2 ; Poisson’s ratio ν 12 ; and shear modulus G 12 of the orthotropic medium [1]. The constitutive law (16) takes into account the strain gradient effect. Analogously, a similar relation can be written for the shear stresses τ ( k ) = τ x z ( k ) τ y z ( k ) T
τ ( k ) = 1 2 2 Q ¯ s ( k ) γ ,
in which Q ¯ s ( k ) is the stiffness coefficients matrix related to the shear shown below
Q ¯ s ( k ) = Q ¯ 44 Q ¯ 45 Q ¯ 45 Q ¯ 55 ( k ) ,
whose terms are computed as a function of the shear moduli G 13 , G 23 [1]. It is important to highlight that both membrane and shear stresses are characterized by a classical and a nonlocal part, which is the one multiplied by .
The external loads can be collected into the vector q = q x q y q z M x M y T , which includes five load components. Consequently, the work done by external forces δ L can be written as
δ L = A δ u T q d A .
The governing equations can be obtained by conveniently introducing the stress resultants as the integrals of the stress components along the thickness of the layer. The following quantities are defined:
N = N x x N y y N x y T = k = 1 N L z k z k + 1 σ ( k ) d z , M = M x x M y y M x y T = k = 1 N L z k z k + 1 σ ( k ) z d z , P = P x x P y y P x y T = k = 1 N L z k z k + 1 σ ( k ) F d z , Q = Q x Q y T = k = 1 N L z k z k + 1 τ ( k ) d z , R = R x R y T = k = 1 N L z k z k + 1 τ ( k ) F d z .
The system of five differential equations in terms of stress resultants that governs the static behavior of the plates is carried out by performing the proper manipulations starting from the principle of virtual displacements [77]. By using a compact matrix form, it becomes
I ( 0 ) T D ( m ) T N + I ( 1 ) T D ( m ) T M * + c 1 I ( 3 ) T D ( b ) T P I s ( 1 ) T Q * + I s ( 3 ) T D s ( 0 ) T Q * + q = 0 ,
where
M * = M c 1 P , Q * = Q c 1 R .
Boundary conditions in terms of primary and secondary variables are obtained as shown in the book by Reddy [1], since they are not affected by the strain gradient effect.
It is convenient at this point to express the stress resultants in the fundamental Equation (22) as a function of the displacement vector u , recalling the constitutive laws (16) and (18), as well as the definitions of the strains shown in (4). The following relations are carried out, as far as the membrane stress resultants are concerned:
N = { A D ( m ) I ( 0 ) + B D ( m ) I ( 1 ) c 1 E D ( m ) I ( 1 ) + D ( b ) I ( 3 ) 2 [ A D x x ( m ) I ( 0 ) + A D y y ( m ) I ( 0 ) + B D x x ( m ) I ( 1 ) + B D y y ( m ) I ( 1 ) c 1 E D x x ( m ) I ( 1 ) + D y y ( m ) I ( 1 ) + D x x ( b ) I ( 3 ) + D y y ( b ) I ( 3 ) ] } u ,
M = { B D ( m ) I ( 0 ) + D D ( m ) I ( 1 ) c 1 F D ( m ) I ( 1 ) + D ( b ) I ( 3 ) 2 [ B D x x ( m ) I ( 0 ) + B D y y ( m ) I ( 0 ) + D D x x ( m ) I ( 1 ) + D D y y ( m ) I ( 1 ) c 1 F D x x ( m ) I ( 1 ) + D y y ( m ) I ( 1 ) + D x x ( b ) I ( 3 ) + D y y ( b ) I ( 3 ) ] } u ,
P = { E D ( m ) I ( 0 ) + F D ( m ) I ( 1 ) c 1 H D ( m ) I ( 1 ) + D ( b ) I ( 3 ) 2 [ E D x x ( m ) I ( 0 ) + E D y y ( m ) I ( 0 ) + F D x x ( m ) I ( 1 ) + F D y y ( m ) I ( 1 ) c 1 H D x x ( m ) I ( 1 ) + D y y ( m ) I ( 1 ) + D x x ( b ) I ( 3 ) + D y y ( b ) I ( 3 ) ] } u ,
where the terms into the constitutive matrices A , B , D , E , F , H are given by
A , B , D , E , F , H i j = k = 1 N L z k z k + 1 Q ¯ i j ( k ) 1 , z , z 2 , F , z F , F 2 d z
for i , j = 1 , 2 , 6 . The following differential operators that appear due to the Laplacian are also introduced: D x x ( m ) = P x y ( 20 ) D ( m ) , D y y ( m ) = P x y ( 02 ) D ( m ) and D x x ( b ) = P x y ( 20 ) D ( b ) , D y y ( b ) = P x y ( 02 ) D ( b ) , in which P x y ( p q ) is given by
P x y ( p q ) = p + q x p y q .
On the other hand, the stress resultants related to shear forces are defined as follows:
Q = { A s I s ( 1 ) + A s D s ( 0 ) I s ( 3 ) c 1 L s I s ( 1 ) + D s ( 0 ) I s ( 3 ) 2 [ A s P x y ( 20 ) I s ( 1 ) + A s P x y ( 02 ) I s ( 1 ) + A s D s x x ( 0 ) I s ( 3 ) + A s D s y y ( 0 ) I s ( 3 ) c 1 L s P x y ( 20 ) I s ( 1 ) + P x y ( 02 ) I s ( 1 ) + D s x x ( 0 ) I s ( 3 ) + D s y y ( 0 ) I s ( 3 ) ] } u ,
R = { L s I s ( 1 ) + L s D s ( 0 ) I s ( 3 ) c 1 N s I s ( 1 ) + D s ( 0 ) I s ( 3 ) 2 [ L s P x y ( 20 ) I s ( 1 ) + L s P x y ( 02 ) I s ( 1 ) + L s D s x x ( 0 ) I s ( 3 ) + L s D s y y ( 0 ) I s ( 3 ) c 1 N s P x y ( 20 ) I s ( 1 ) + P x y ( 02 ) I s ( 1 ) + D s x x ( 0 ) I s ( 3 ) + D s y y ( 0 ) I s ( 3 ) ] } u ,
in which the terms collected in the constitutive matrices A s , L s , N s , for i , j = 4 , 5 , are given by
A s , L s , N s i j = k = 1 N L z k z k + 1 κ s Q s i j ( k ) 1 , F , F 2 d z ,
where κ s is the shear correction factor. Its value is different from the unity only for the FSDT, in which it is assumed equal to 5/6. The differential operators D s x x ( 0 ) , D s y y ( 0 ) are computed as D s x x ( 0 ) = P x y ( 20 ) D s ( 0 ) , D s y y ( 0 ) = P x y ( 02 ) D s ( 0 ) .

3. Solution Procedure

Once the nonlocal governing equations are written in terms of the displacements collected in u , they can be solved analytically by means of the Navier methodology. As illustrated in [1], this approach can be applied only for some peculiar lamination schemes, which are antisymmetric cross-ply and antisymmetric angle-ply, respectively. These two cases are analyzed separately in the following, assuming simply supported boundary conditions for both circumstances. The solution to the current static problem is provided by the algebraic linear system shown below
K Δ = F ,
in which Δ = U m n V m n W m n X m n Y m n T is the vector that collects the unknown coefficients that determine the displacement amplitudes. On the other hand, K and F are the stiffness matrix and the load vector, respectively. The terms included in the symmetric matrix K will be specified for each lamination scheme. On the other hand, the load vector has the following definition, assuming that only transverse surface forces are applied: F = 0 0 Q m n 0 0 T , where Q m n is equal to q z for a sinusoidally distributed load, having introduced the same expansion also for the applied external forces [1].

3.1. Cross-Ply Laminates

In order to apply the Navier approach for the cross-ply sequence at issue, the stiffnesses written below are all equal to zero:
A 16 = A 26 = A s 45 = B 16 = B 26 = D 16 = D 26 = 0 , E 16 = E 26 = F 16 = F 26 = H 16 = H 26 = L s 45 = N s 45 = 0 .
In this circumstance, the solution can be sought assuming the following expansion:
u = n = 1 m = 1 Δ cos α x sin β y sin α x cos β y sin α x sin β y cos α x sin β y sin α x cos β y ,
where ○ stands for the elementwise product and α = m π / a , β = n π / b . The explicit expressions of the coefficients K i j included in the stiffness matrix K , for i , j = 1 , , 5 , whose size is clearly 5 × 5 are listed in Appendix A, recalling that K i j = K j i .

3.2. Angle-Ply Laminates

The Navier approach can be applied for angle-ply laminates if the following stiffnesses are all equal to zero:
A 16 = A 26 = A s 45 = B 11 = B 12 = B 22 = B 66 = D 16 = D 26 = 0 , E 11 = E 12 = E 22 = E 66 = F 16 = F 26 = H 16 = H 26 = L s 45 = N s 45 = 0 .
In this case, the solution is obtained assuming the following expansion:
u = n = 1 m = 1 Δ sin α x cos β y cos α x sin β y sin α x sin β y cos α x sin β y sin α x cos β y .
The expressions of the coefficients K i j , for i , j = 1 , , 5 , of the stiffness matrix K , are listed in Appendix A, presenting only those terms that are different from the ones valid for cross-ply sequences.
The solution in terms of Δ can be easily obtained from Equation (32), for both laminates under consideration. Once the amplitudes in Δ are computed, definitions (34) and (36) allow us to obtain the displacements within the reference domain. Consequently, the strains can be also deduced. The membrane stress components σ ( k ) can be evaluated as well, following the procedure illustrated in the book by Reddy [1] through the constitutive relations previously presented. On the other hand, the shear stress components τ ( k ) are determined from the three-dimensional equilibrium elasticity equations [1,73]. The complete procedure is omitted for conciseness purposes.

4. Numerical Results

The current Section aims to present the results of the static analyses. Due to the general features of the theoretical approach, the solutions are presented for different nonlocal theories, which are CLPT, FSDT and TSDT, setting properly the values of c 1 and F . As far as the mechanical features are concerned, the ratio between the longitudinal and transverse Young’s moduli E 1 / E 2 is specified in each application, whereas the other quantities are taken as G 12 = G 13 = 0.5 E 2 , G 23 = 0.2 E 2 , ν 12 = 0.25 for Material 1, or G 12 = G 13 = 0.6 E 2 , G 23 = 0.5 E 2 , ν 12 = 0.25 for Material 2. The lamination schemes, instead, are denoted by ( θ ( 1 ) / / θ ( k ) / / θ ( N L ) ) , where θ ( k ) stands for the orientation of the k-th layer. The results are presented in dimensionless form. In particular, the central deflection w ¯ is given by
w ¯ = w a 2 , b 2 E 2 h 3 a 4 q z .
On the other hand, the stress components are evaluated as follows, unless differently specified:
σ ¯ x x = σ x x a 2 , b 2 , h 2 h 2 b 2 q z , σ ¯ y y = σ y y a 2 , b 2 , h 4 h 2 b 2 q z , σ ¯ x y = σ x y 0 , 0 , h 2 h 2 b 2 q z , σ ¯ x z = σ x z 0 , b 2 , 0 h b q z , σ ¯ y z = σ y z a 2 , 0 , 0 h b q z .
It should be specified that the values of the stresses presented in this Section are all related to the classical component, which can be deducted from definitions (16) and (18) following the approach used in [73,77]. The analyses are carried out for increasing values of the dimensionless nonlocal parameter / a 2 , in order to show the effect of the strain gradient on the static solutions.
The first application aims to investigate the central deflection w ¯ as a function of side-to-thickness ratio a / h of a square plate for different lamination schemes: cross-ply ( 0 / 90 / 90 / 0 ) and angle-ply ( 45 / 45 ) . The results are shown in Figure 2, assuming E 1 / E 2 = 25 for the cross-ply (Material 1) and E 1 / E 2 = 40 for the angle-ply (Material 2).
The graphs include the structural response obtained by means of the three theories. Each model is related to a different line style: solid line for TSDT, dotted for FSDT, and dash-dotted for CLPT. The color, instead, is linked to the value of the nonlocal parameter. The same choice is also kept in the following figures. It can be observed that the greater is the nonlocal effect and the lower is the vertical deflection, independently from the theory. In other words, the central deflection is reduced by increasing the value of the nonlocal parameter ( / a ) 2 . The FSDT and TSDT, moreover, are characterized by a comparable behavior and are highly affected by plate thickness. By increasing the ratio a / h , their displacements tend to the results of the CLPT, which do not depend on that geometric ratio. The corresponding curves, in fact, are described by rectilinear functions. Similar behaviors are obtained in both lamination schemes.
In the next test, a ( 0 / 90 / 0 ) cross-ply square plate is considered. Material 1 is taken into account in this circumstance, assuming E 1 / E 2 = 25 . The results are presented in Table 1 for different values of a / h , varying the structural theory. Where available, the analytical solutions by Reddy [1] are provided, in terms of both displacement and stress components. The reference results are clearly presented only for classical elasticity, assuming ( / a ) 2 = 0 . The comparison proves a very good agreement between the current approach and the reference one.
Table 2, instead, aims to extend these results to the nonlocal elasticity framework, varying the ratio ( / a ) 2 . The same reduction of the displacement values observed in Figure 2 can also be seen numerically in this circumstance. The three theories provide really close results for thin plates even if the strain gradient effect is taken into account. It should be noted from these tables that noticeable differences can be seen especially in terms of membrane stresses in thick configurations (defined by a / h = 4 ), if the results related to the different theories are compared. As it will be highlighted in the following paragraphs, this is due to the fact that the TSDT is characterized by nonlinear stress profiles. This nonlinearity is particularly emphasized for thicker plates, whereas it is reduced for lower values of a / h .
A ( 0 / 90 / 90 / 0 ) cross-ply lamination scheme is considered in the next application. Even in this case, Material 1 is taken into account to describe the orthotropic features of the layers, with E 1 / E 2 = 25 . The results are shown in Table 3 for the classical elasticity. In the same Table, the solutions shown in [1] are presented for comparison purposes.
As in the previous test, the present configuration is analyzed also in the nonlocal framework. Results for increasing values of ( / a ) 2 are included in Table 4 for the sake of completeness.
The reference solutions are all related to classical elasticity. In the following application, the current methodology is also verified with regard to nonlocal elasticity, taking into account the solutions presented in [77]. To this aim, only CLPT is considered in accordance with the reference paper, due to the availability of the literature. Several cross-ply lamination schemes are considered, assuming E 1 / E 2 = 25 and Material 1 as orthotropic features. The results are presented in Table 5 for various stacking sequences, varying the value of ( / a ) 2 , in terms of central deflections and membrane stress components. In this circumstance, the following dimensionless quantities are considered for the stresses:
σ ¯ y y = σ y y a 2 , b 2 , h 2 h 2 b 2 q z , σ ¯ x y = σ x y a , b , h 2 h 2 b 2 q z ,
whereas the same expression introduced before is used for σ ¯ x x . The results shown in Table 5 prove that the present solutions are in very good agreement with the reference ones, also in the framework of nonlocal elasticity.
The next application is focused on the bending analysis of antisymmetric angle ply laminates, characterized by ( θ / θ / ) as the lamination scheme and θ being the arbitrary orientation of each layer. Firstly, the analyses are presented in terms of dimensionless central deflection w ¯ , considering square plates made of Material 2 and E 1 / E 2 = 40 as orthotropic ratio. The results are presented for several side-to-thickness values a / h , orthotropic angles θ , and different structural theories in Table 6, as far as classical elasticity is concerned, which means ( / a ) 2 = 0 . The values are in good agreement with the ones taken as references [1].
As in the previous cases, the analyses are extended to nonlocal elasticity by increasing the value of ( / a ) 2 but keeping the same geometric ratios and mechanical features. The results are presented in Table 7. Even in these circumstances, the differences between TSDT and FSDT decrease for lower values of thickness, and the displacements tend to the ones of the CLPT.
A ( 45 / 45 ) 4 laminated plate square plate is considered in the next application. Each orthotropic layer is made of Material 1, with E 1 / E 2 = 25 . This test aims to evaluate the stress components in angle-ply configurations. With respect to the dimensionless values shown in (38), only the transverse shear stress σ ¯ x z is specified in a different thickness coordinate as specified below
σ ¯ x z = σ x z 0 , b 2 , h 2 h b q z .
It should be recalled that in the following application σ ¯ x x = σ ¯ y y and σ ¯ x z = σ ¯ y z , therefore, their values are not repeated twice. The results related to the classical elasticity are shown in Table 8, varying a / h and the structural theory. Even if the displacements become closer independently from the considered approach, the stresses assume different values especially if the thickness is greater. In this circumstance, the reference solutions are available only for the FSDT [1].
On the other hand, the nonlocal counterpart of the current application is shown in Table 9. The increasing values of ( / a ) 2 emphasize the differences in terms of stresses, especially if the plates are thicker.
The last tests aim to present the stress analysis in terms of the through-the-thickness distributions of the various components, highlighting the differences that could arise by varying the structural approach (TSDT, FSDT, and CLPT) for different geometric ratios. The effect of ( / a ) 2 is also investigated. Firstly, a ( 0 / 90 / 90 / 0 ) cross-ply square plate is analyzed, for a / h equal to 4 and 10, respectively. The orthotropic features of the layers are obtained by setting E 1 / E 2 = 25 and selecting Material 1 as constituent. The membrane and shear dimensionless stresses are given by
σ ¯ x x = σ x x a 2 , b 2 , z ¯ h 2 b 2 q z , σ ¯ y y = σ y y a 2 , b 2 , z ¯ h 2 b 2 q z , σ ¯ x y = σ x y 0 , 0 , z ¯ h 2 b 2 q z , σ ¯ x z = σ x z 0 , b 2 , z ¯ h b q z , σ ¯ y z = σ y z a 2 , 0 , z ¯ h b q z ,
where z ¯ = 2 z / h stands for the dimensionless thickness coordinate. The through-the-thickness stress distributions in Figure 3 are related to the case a / h = 4 , which is representative for thick plates. It can be observed that in both classical and nonlocal elasticity, the TSDT is characterized by significantly different profiles due to the higher-order features of the displacement field. This aspect is more evident in the membrane stress components, which are obtained by means of the application of constitutive laws. In fact, it should be recalled that the TSDT allows a cubic representation of the stress profiles, whereas a linear variation is associated with the FSDT and CLPT. This feature gives rise to noticeable differences if thick plates are investigated, as it can be seen from the plot of σ ¯ x x in Figure 3. The variation of the transverse stresses, instead, is always characterized by curved and continuous profiles since they are equilibrium-derived, according to the procedure shown in [1]. A lower value of thickness is considered in the graphs shown in Figure 4, assuming a / h = 10 . Due to this choice, the stress distributions tend to the same value, for both ( / a ) 2 = 0 and ( / a ) 2 = 0.10 , independently from the theory. Therefore, the differences among the theories is practically negligible starting from a / h = 10 , which is typically the geometric ratio that characterizes moderately thick plates.
A similar analysis is carried out for a ( 45 / 45 ) angle-ply square plate, for a / h equal to 4 and 10, respectively. In this case, the orthotropic features of the layers are given by E 1 / E 2 = 40 . Material 2 is taken into account as constituent. With respect to the dimensionless expressions introduced in (41), a different value of σ ¯ y z is considered, which is defined below:
σ ¯ y z = σ y z 0 , b 2 , z ¯ h b q z .
The graphical plots are shown in Figure 5 and Figure 6 for a / h = 4 and a / h = 10 , respectively. It can be observed that in these circumstances, the FSDT and CLPT are overlapped. The nonlinear behavior of TSDT is clearly more evident for thicker plates, for both classical and nonlocal elasticity.

5. Conclusions

A theoretical framework able to simultaneously and accurately study thick and thin laminated composite plates has been presented in the paper. In particular, the chosen kinematic model is able to describe several theories if properly set, which are the CLPT, FSDT, and TSDT. The theories have been modified to include the strain gradient effect, in order to take into account nonlocal contributions in the evaluation of stresses. The proposed approach is general, and the fundamental static equations have been written for arbitrary configurations by using a compact matrix notation. In order to provide an analytical solution, the Navier methodology has been applied to deal with cross-ply and angle-ply lamination schemes. The explicit definitions needed in the solution procedure have been presented, highlighting also the contributions linked to the strain gradient effect. The approach has been validated numerically by means of comparison with the results accessible in the literature, if available, for both classical and nonlocal elasticity. The results have been presented in terms of displacements and stresses, varying the geometric features, mechanical properties, lamination schemes, as well as the influence of the nonlocal effect. The differences that have been observed among the theories have been emphasized. In particular, the results have proven that the stress components are noticeably different by changing the plate theory, especially for higher values of thickness. Finally, it should be specified that the proposed solutions could be used for further advancements of this topic in nonlocal elasticity and could be taken as benchmarks for future comparisons, especially if numerical methods are developed.

Author Contributions

Conceptualization, M.B. and A.M.T.; methodology, M.B. and A.M.T.; software, M.B. and A.M.T.; validation, M.B. and A.M.T.; formal analysis, M.B. and A.M.T.; investigation, M.B. and A.M.T.; resources, M.B. and A.M.T.; data curation, M.B. and A.M.T.; writing—original draft preparation, M.B. and A.M.T.; writing—review and editing, M.B. and A.M.T.; visualization, M.B. and A.M.T.; supervision, M.B. and A.M.T.; project administration, M.B. and A.M.T.; funding acquisition, M.B. and A.M.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Acknowledgments

Financial support from the Italian Ministry of University and Research (MIUR) in the framework of the Project FISR 2019: “Eco_Earth” (code 00245) is gratefully acknowledged. Financial support from the Italian Ministry of Education, University and Research (MIUR) in the framework of the Project PRIN “Modelling of constitutive laws for traditional and innovative building materials” (code 2017HFPKZY) is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Algebraic Form of the Stiffness Matrix

As specified in (32), the static problems at issue are governed by the algebraic linear system K Δ = F . The terms included in the symmetric matrix K are presented in this Appendix, for the two cases investigated in the paper.

Appendix A.1. Cross-Ply Laminates

As far as cross-ply lamination schemes are concerned, the terms of the first row of the matrix K are given by
K 11 = A 11 α 2 + A 66 β 2 + 2 [ A 11 ( α 4 + α 2 β 2 ) + A 66 β 4 + α 2 β 2 ] ,
K 12 = A 12 + A 66 α β + 2 [ A 12 + A 66 ( α 3 β + α β 3 ) ] ,
K 13 = c 1 { E 11 α 3 + E 12 + 2 E 66 α β 2 + 2 [ E 11 ( α 5 + α 3 β 2 ) + E 12 + 2 E 66 ( α 3 β 2 + α β 4 ) ] } ,
K 14 = B 11 c 1 E 11 α 2 + B 66 c 1 E 66 β 2 + 2 [ B 11 c 1 E 11 ( α 4 + α 2 β 2 ) + B 66 c 1 E 66 ( β 4 + α 2 β 2 ) } ,
K 15 = B 12 c 1 E 12 + B 66 c 1 E 66 α β + 2 [ B 12 c 1 E 12 + B 66 c 1 E 66 ( α 3 β + α β 3 ) ] .
The terms of the second row are defined as
K 22 = A 66 α 2 + A 22 β 2 + 2 [ A 66 ( α 4 + α 2 β 2 ) + A 22 ( β 4 + α 2 β 2 ) ] ,
K 23 = c 1 { E 22 β 3 + E 12 + 2 E 66 α 2 β + 2 [ E 22 ( α β 3 + β 5 ) + E 12 + 2 E 66 ( α 4 β + α 2 β 3 ) ] } ,
K 24 = B 12 c 1 E 12 + B 66 c 1 E 66 α β + 2 [ B 12 c 1 E 12 + B 66 c 1 E 66 ( α 3 β + α β 3 ) ] ,
K 25 = B 22 c 1 E 22 β 2 + B 66 c 1 E 66 α 2 + 2 [ B 22 c 1 E 22 ( α 2 β 2 + β 4 ) + B 66 c 1 E 66 ( α 4 + α 2 β 2 ) ] .
As far as the terms of the third row are concerned, one gets
K 33 = ( A s 44 2 c 1 L s 44 + c 1 2 N s 44 ) α 2 + ( A s 55 2 c 1 L s 55 + c 1 2 N s 55 ) β 2 + 2 [ ( A s 44 2 c 1 L s 44 + c 1 2 N s 44 ) ( α 4 + α 2 β 2 ) + ( A s 55 2 c 1 L s 55 + c 1 2 N s 55 ) ( β 4 + α 2 β 2 ) ] + c 1 2 { H 11 α 4 + 2 H 12 + 2 H 66 α 2 β 2 + H 22 β 4 + 2 [ H 11 ( α 6 + α 4 β 2 ) + 2 H 12 + 2 H 66 ( α 4 β 2 + α 2 β 4 ) + H 22 ( β 6 + α 2 β 4 ) ] } ,
K 34 = ( A s 44 2 c 1 L s 44 + c 1 2 N s 44 ) α + 2 [ ( A s 44 2 c 1 L s 44 + c 1 2 N s 44 ) ( α 3 + α β 2 ) ] c 1 { F 11 c 1 H 11 α 3 + F 12 c 1 H 12 + 2 F 66 2 c 1 H 66 α β 2 + 2 [ F 11 c 1 H 11 ( α 5 + α 3 β 2 ) + F 12 c 1 H 12 + 2 F 66 2 c 1 H 66 ( α 3 β 2 + α β 4 ) ] } ,
K 35 = ( A s 55 2 c 1 L s 55 + c 1 2 N s 55 ) α + 2 [ ( A s 55 2 c 1 L s 55 + c 1 2 N s 55 ) ( α 3 + α β 2 ) ] c 1 { F 22 c 1 H 22 β 3 + F 12 c 1 H 12 + 2 F 66 2 c 1 H 66 α 2 β + 2 [ F 22 c 1 H 22 ( β 5 + α 2 β 3 ) + F 12 c 1 H 12 + 2 F 66 2 c 1 H 66 ( α 4 β + α 2 β 3 ) ] } .
Instead, the terms of the fourth row are given by
K 44 = A s 44 2 c 1 L s 44 + c 1 2 N s 44 + 2 [ ( A s 44 2 c 1 L s 44 + c 1 2 N s 44 ) ( α 2 + β 2 ) ] + ( D 11 2 c 1 F 11 + c 1 2 H 11 ) α 2 + ( D 66 2 c 1 F 66 + c 1 2 H 66 ) β 2 + ( D 11 2 c 1 F 11 + c 1 2 H 11 ) ( α 4 + α 2 β 2 ) + ( D 66 2 c 1 F 66 + c 1 2 H 66 ) ( β 4 + α 2 β 2 ) ] } ,
K 45 = ( D 12 2 c 1 F 12 + c 1 2 H 12 ) α β + ( D 66 2 c 1 F 66 + c 1 2 H 66 ) α β + 2 [ ( D 12 2 c 1 F 12 + c 1 2 H 12 ) ( α 3 β + α β 3 ) + ( D 66 2 c 1 F 66 + c 1 2 H 66 ) ( α 3 β + α β 3 ) ] } .
Finally, the terms included in the last row assume the following aspect:
K 55 = A s 55 2 c 1 L s 55 + c 1 2 N s 55 + 2 [ ( A s 55 2 c 1 L s 55 + c 1 2 N s 55 ) ( α 2 + β 2 ) ] + ( D 66 2 c 1 F 66 + c 1 2 H 66 ) α 2 + ( D 22 2 c 1 F 22 + c 1 2 H 22 ) β 2 + ( D 66 2 c 1 F 66 + c 1 2 H 66 ) ( α 4 + α 2 β 2 ) + ( D 22 2 c 1 F 22 + c 1 2 H 22 ) ( β 4 + α 2 β 2 ) ] } .

Appendix A.2. Angle-Ply Laminates

The coefficients included in the stiffness matrix are presented below for angle-ply laminates. When omitted, their definitions are the same as those presented for the cross-ply sequences. The terms of the first row are given by
K 13 = c 1 { 3 E 16 α 2 β + E 26 β 3 + 2 [ 3 E 16 ( α 4 β + α 2 β 3 ) + E 26 ( α 2 β 3 + β 5 ) ] } ,
K 14 = 2 B 16 c 1 E 16 α β + 2 [ 2 B 16 c 1 E 16 ( α 3 β + α β 3 ) ] ,
K 15 = B 16 c 1 E 16 α 2 + B 26 c 1 E 26 β 2 + 2 [ B 16 c 1 E 16 ( α 4 + α 2 β 2 ) + B 26 c 1 E 26 ( α 2 β 2 + β 4 ) ] .
The terms of the second row are defined as
K 23 = c 1 { E 16 α 3 + 3 E 26 α β 2 + 2 [ E 16 ( α 5 + α 3 β 2 ) + 3 E 26 ( α 3 β 2 + α β 4 ) ] } ,
K 24 = B 16 c 1 E 16 α 2 + B 26 c 1 E 26 β 2 + 2 [ B 16 c 1 E 16 ( α 4 + α 2 β 2 ) + B 26 c 1 E 26 ( α 2 β 2 + β 4 ) ] ,
K 25 = 2 B 26 c 1 E 26 α β + 2 [ 2 B 26 c 1 E 26 ( α 3 β + α β 3 ) ] .

References

  1. Reddy, J.N. Mechanics of Laminated Composite Plates and Shells: Theory and Analysis, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2004. [Google Scholar]
  2. Amabili, M.; Reddy, J.N. The nonlinear, third-order thickness and shear deformation theory for statics and dynamics of laminated composite shells. Compos. Struct. 2020, 244, 112265. [Google Scholar] [CrossRef]
  3. Reddy, J.N. A review of refined theories of laminated composite plates. Shock Vib. Dig. 1990, 22, 3–17. [Google Scholar] [CrossRef]
  4. Reddy, J.N. On refined theories of composite laminates. Meccanica 1990, 25, 230–238. [Google Scholar] [CrossRef]
  5. Patni, M.; Minera, S.; Groh, R.M.J.; Pirrera, A.; Weaver, P.M. Three-dimensional stress analysis for laminated composite and sandwich structures. Compos. Part B Eng. 2018, 155, 299–328. [Google Scholar] [CrossRef]
  6. Hii, A.K.W.; Minera, S.; Groh, R.M.J.; Pirrera, A.; Kawashita, L.F. Three-dimensional stress analyses of complex laminated shells with a variable-kinematics continuum shell element. Compos. Struct. 2019, 229, 111405. [Google Scholar] [CrossRef]
  7. Phan, N.D.; Reddy, J.N. Analysis of laminated composite plates using a higher-order shear deformation theory. Int. J. Numer. Methods Eng. 1985, 21, 2201–2219. [Google Scholar] [CrossRef]
  8. Gutierrez Rivera, M.; Reddy, J.N.; Amabili, M. A continuum eight-parameter shell finite element for large deformation analysis. Mech. Adv. Mater. Struct. 2020, 27, 551–560. [Google Scholar] [CrossRef]
  9. Nguyen, H.N.; Tan, T.C.; Luat, D.T.; Phan, V.D.; Thom, D.V.; Minh, P.V. Research on the buckling behavior of functionally graded plates with stiffeners based on the third-order shear deformation theory. Materials 2019, 12, 1262. [Google Scholar] [CrossRef] [Green Version]
  10. Qin, B.; Zhao, X.; Liu, H.; Yu, Y.; Wang, Q. Free vibration analysis of curved laminated composite beams with different shapes, lamination schemes, and boundary conditions. Materials 2020, 13, 1010. [Google Scholar] [CrossRef] [Green Version]
  11. Valencia Murillo, C.; Gutierrez Rivera, M.; Reddy, J.N. Linear Vibration Analysis of Shells Using a Seven-Parameter Spectral/hp Finite Element Model. Appl. Sci. 2020, 10, 5102. [Google Scholar] [CrossRef]
  12. Mishra, B.B.; Kumar, A.; Zaburko, J.; Sadowska-Buraczewska, B.; Barnat-Hunek, D. Dynamic Response of Angle Ply Laminates with Uncertainties Using MARS, ANN-PSO, GPR and ANFIS. Materials 2021, 14, 395. [Google Scholar] [CrossRef]
  13. Petrolo, M.; Carrera, E. Best theory diagrams for multilayered structures via shell finite elements. Adv. Model. Simul. Eng. Sci. 2019, 6, 1–23. [Google Scholar] [CrossRef]
  14. Nobili, A.; Falope, F.O. Impregnated carbon fabric–reinforced cementitious matrix composite for rehabilitation of the Finale Emilia hospital roofs: Case study. J. Compos. Constr. 2017, 21, 05017001. [Google Scholar] [CrossRef] [Green Version]
  15. Signorini, C.; Nobili, A.; Falope, F.O. Mechanical performance and crack pattern analysis of aged Carbon Fabric Cementitious Matrix (CFRCM) composites. Compos. Struct. 2018, 202, 1114–1120. [Google Scholar] [CrossRef]
  16. Falope, F.O.; Lanzoni, L.; Tarantino, A.M. Double lap shear test on steel fabric reinforced cementitious matrix (SFRCM). Compos. Struct. 2018, 201, 503–513. [Google Scholar] [CrossRef]
  17. Falope, F.O.; Lanzoni, L.; Tarantino, A.M. Modified hinged beam test on steel fabric reinforced cementitious matrix (SFRCM). Compos. Part B Eng. 2018, 146, 232–243. [Google Scholar] [CrossRef]
  18. Arbind, A.; Reddy, J.N. A general higher-order shell theory for compressible isotropic hyperelastic materials using orthonormal moving frame. Int. J. Numer. Methods Eng. 2021, 122, 235–269. [Google Scholar] [CrossRef]
  19. Liguori, F.S.; Zucco, G.; Madeo, A.; Garcea, G.; Leonetti, L.; Weaver, P.M. An isogeometric framework for the optimal design of variable stiffness shells undergoing large deformations. Int. J. Solids Struct. 2021, 210, 18–34. [Google Scholar] [CrossRef]
  20. Bacciocchi, M.; Tarantino, A.M. Natural Frequency Analysis of Functionally Graded Orthotropic Cross-Ply Plates Based on the Finite Element Method. Math. Comput. Appl. 2019, 24, 52. [Google Scholar] [CrossRef] [Green Version]
  21. Bacciocchi, M.; Luciano, R.; Majorana, C.; Tarantino, A.M. Free vibrations of sandwich plates with damaged soft-core and non-uniform mechanical properties: Modeling and finite element analysis. Materials 2019, 12, 2444. [Google Scholar] [CrossRef] [Green Version]
  22. Bacciocchi, M.; Tarantino, A.M. Critical buckling load of honeycomb sandwich panels reinforced by three-phase orthotropic skins enhanced by carbon nanotubes. Compos. Struct. 2020, 237, 111904. [Google Scholar] [CrossRef]
  23. Moleiro, F.; Carrera, E.; Li, G.; Cinefra, M.; Reddy, J.N. Hygro-thermo-mechanical modelling of multilayered plates: Hybrid composite laminates, fibre metal laminates and sandwich plates. Compos. Part B Eng. 2019, 177, 107388. [Google Scholar] [CrossRef]
  24. Moleiro, F.; Carrera, E.; Ferreira, A.; Reddy, J.N. Hygro-thermo-mechanical modelling and analysis of multilayered plates with embedded functionally graded material layers. Compos. Struct. 2020, 233, 111442. [Google Scholar] [CrossRef]
  25. Gorgeri, A.; Vescovini, R.; Dozio, L. Sublaminate variable kinematics shell models for functionally graded sandwich panels: Bending and free vibration response. Mech. Adv. Mater. Struct. 2020, 1–18. [Google Scholar] [CrossRef]
  26. Carrera, E. A refined multilayered finite-element model applied to linear and non-linear analysis of sandwich plates. Compos. Sci. Technol. 1998, 58, 1553–1569. [Google Scholar] [CrossRef]
  27. Carrera, E. Theories and finite elements for multilayered plates and shells: A unified compact formulation with numerical assessment and benchmarking. Arch. Comput. Methods Eng. 2003, 10, 215–296. [Google Scholar] [CrossRef]
  28. Carrera, E. Historical review of zig-zag theories for multilayered plates and shells. Appl. Mech. Rev. 2003, 56, 287–308. [Google Scholar] [CrossRef]
  29. Carrera, E.; Giunta, G. Refined beam theories based on a unified formulation. Int. J. Appl. Mech. 2010, 2, 117–143. [Google Scholar] [CrossRef]
  30. Carrera, E.; Kaleel, I.; Petrolo, M. Elastoplastic analysis of compact and thin-walled structures using classical and refined beam finite element models. Mech. Adv. Mater. Struct. 2019, 26, 274–286. [Google Scholar] [CrossRef]
  31. Di Maida, P.; Falope, F.O. Euler-Bernoulli nanobeam welded to a compressible semi-infinite substrate. Model. Simul. Eng. 2016, 2016. [Google Scholar] [CrossRef] [Green Version]
  32. Jankowski, P.; Żur, K.K.; Kim, J.; Reddy, J.N. On the bifurcation buckling and vibration of porous nanobeams. Compos. Struct. 2020, 250, 112632. [Google Scholar] [CrossRef]
  33. Farajpour, A.; Żur, K.K.; Kim, J.; Reddy, J.N. Nonlinear frequency behaviour of magneto-electromechanical mass nanosensors using vibrating MEE nanoplates with multiple nanoparticles. Compos. Struct. 2021, 260, 113458. [Google Scholar] [CrossRef]
  34. Żur, K.K.; Arefi, M.; Kim, J.; Reddy, J.N. Free vibration and buckling analyses of magneto-electro-elastic FGM nanoplates based on nonlocal modified higher-order sinusoidal shear deformation theory. Compos. Part B Eng. 2020, 182, 107601. [Google Scholar] [CrossRef]
  35. Wang, B.; Zhou, S.; Zhao, J.; Chen, X. A size-dependent Kirchhoff micro-plate model based on strain gradient elasticity theory. Eur. J. Mech. A/Solids 2011, 30, 517–524. [Google Scholar] [CrossRef]
  36. Lazopoulos, K. On bending of strain gradient elastic micro-plates. Mech. Res. Commun. 2009, 36, 777–783. [Google Scholar] [CrossRef]
  37. Lu, L.; Guo, X.; Zhao, J. A unified size-dependent plate model based on nonlocal strain gradient theory including surface effects. Appl. Math. Model. 2019, 68, 583–602. [Google Scholar] [CrossRef]
  38. Brands, B.; Davydov, D.; Mergheim, J.; Steinmann, P. Reduced-Order Modelling and Homogenisation in Magneto-Mechanics: A Numerical Comparison of Established Hyper-Reduction Methods. Math. Comput. Appl. 2019, 24, 20. [Google Scholar] [CrossRef] [Green Version]
  39. Huang, W.; Xu, R.; Yang, J.; Huang, Q.; Hu, H. Data-driven multiscale simulation of FRP based on material twins. Compos. Struct. 2021, 256, 113013. [Google Scholar] [CrossRef]
  40. Xu, M.; Gitman, I.M.; Wei, P.; Askes, H. Finite element implementation of a multi-scale dynamic piezomagnetic continuum model. Comput. Struct. 2020, 240, 106352. [Google Scholar] [CrossRef]
  41. Fantuzzi, N.; Trovalusci, P.; Luciano, R. Multiscale Analysis of Anisotropic Materials with Hexagonal Microstructure as Micropolar Continua. Int. J. Multiscale Comput. Eng. 2020, 18. [Google Scholar] [CrossRef]
  42. Tarantino, A.M. Crack propagation in finite elastodynamics. Math. Mech. Solids 2005, 10, 577–601. [Google Scholar] [CrossRef]
  43. Tarantino, A.M.; Lanzoni, L.; Falope, F.O. The Bending Theory of Fully Nonlinear Beams; Springer: New York, NY, USA, 2019. [Google Scholar]
  44. Falope, F.O.; Lanzoni, L.; Tarantino, A.M. The bending of fully nonlinear beams. Theoretical, numerical and experimental analyses. Int. J. Eng. Sci. 2019, 145, 103167. [Google Scholar] [CrossRef]
  45. Falope, F.O.; Lanzoni, L.; Tarantino, A.M. Bending device and anticlastic surface measurement of solids under large deformations and displacements. Mech. Res. Commun. 2019, 97, 52–56. [Google Scholar] [CrossRef]
  46. Falope, F.O.; Lanzoni, L.; Tarantino, A.M. FE Analyses of Hyperelastic Solids under Large Bending: The Role of the Searle Parameter and Eulerian Slenderness. Materials 2020, 13, 1597. [Google Scholar] [CrossRef] [Green Version]
  47. Falope, F.O.; Lanzoni, L.; Radi, E. Buckling of a Timoshenko beam bonded to an elastic half-plane: Effects of sharp and smooth beam edges. Int. J. Solids Struct. 2020, 185, 222–239. [Google Scholar] [CrossRef] [Green Version]
  48. Eringen, A.C. Nonlocal polar elastic continua. Int. J. Eng. Sci. 1972, 10, 1–16. [Google Scholar] [CrossRef]
  49. Eringen, A.C. On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves. J. Appl. Phys. 1983, 54, 4703–4710. [Google Scholar] [CrossRef]
  50. Luciano, R.; Willis, J. Non-local constitutive response of a random laminate subjected to configuration-dependent body force. J. Mech. Phys. Solids 2001, 49, 431–444. [Google Scholar] [CrossRef]
  51. Barretta, R.; Feo, L.; Luciano, R.; Marotti De Sciarra, F. A gradient Eringen model for functionally graded nanorods. Compos. Struct. 2015, 131, 1124–1131. [Google Scholar] [CrossRef]
  52. Barretta, R.; Feo, L.; Luciano, R.; Marotti de Sciarra, F.; Penna, R. Functionally graded Timoshenko nanobeams: A novel nonlocal gradient formulation. Compos. Part B Eng. 2016, 100, 208–219. [Google Scholar] [CrossRef]
  53. Apuzzo, A.; Barretta, R.; Canadija, M.; Feo, L.; Luciano, R.; Marotti De Sciarra, F. A closed-form model for torsion of nanobeams with an enhanced nonlocal formulation. Compos. Part B Eng. 2017, 108, 315–324. [Google Scholar] [CrossRef]
  54. Apuzzo, A.; Barretta, R.; Luciano, R.; Marotti De Sciarra, F.; Penna, R. Free vibrations of Bernoulli-Euler nano-beams by the stress-driven nonlocal integral model. Compos. Part B Eng. 2017, 123, 105–111. [Google Scholar] [CrossRef]
  55. Tuna, M.; Kirca, M.; Trovalusci, P. Deformation of atomic models and their equivalent continuum counterparts using Eringen’s two-phase local/nonlocal model. Mech. Res. Commun. 2019, 97, 26–32. [Google Scholar] [CrossRef] [Green Version]
  56. Tuna, M.; Trovalusci, P. Scale dependent continuum approaches for discontinuous assemblies: Explicit and implicit non-local models. Mech. Res. Commun. 2020, 103, 103461. [Google Scholar] [CrossRef]
  57. Tuna, M.; Trovalusci, P. Stress distribution around an elliptic hole in a plate with ‘implicit’ and ‘explicit’ non-local models. Compos. Struct. 2021, 256, 113003. [Google Scholar] [CrossRef]
  58. Beheshti, A. Large deformation analysis of strain-gradient elastic beams. Comput. Struct. 2016, 177, 162–175. [Google Scholar] [CrossRef]
  59. Thai, S.; Thai, H.T.; Vo, T.P.; Patel, V.I. Size-dependant behaviour of functionally graded microplates based on the modified strain gradient elasticity theory and isogeometric analysis. Comput. Struct. 2017, 190, 219–241. [Google Scholar] [CrossRef]
  60. Bleyer, J.; Hassen, G. Automated formulation and resolution of limit analysis problems. Comput. Struct. 2021, 243, 106341. [Google Scholar] [CrossRef]
  61. Ashoori, A.; Mahmoodi, M.J. A nonlinear thick plate formulation based on the modified strain gradient theory. Mech. Adv. Mater. Struct. 2018, 25, 813–819. [Google Scholar] [CrossRef]
  62. Kim, J.; Żur, K.K.; Reddy, J.N. Bending, free vibration, and buckling of modified couples stress-based functionally graded porous micro-plates. Compos. Struct. 2019, 209, 879–888. [Google Scholar] [CrossRef]
  63. Thanh, C.L.; Ferreira, A.; Wahab, M.A. A refined size-dependent couple stress theory for laminated composite micro-plates using isogeometric analysis. Thin-Walled Struct. 2019, 145, 106427. [Google Scholar] [CrossRef]
  64. Choi, J.H.; Lee, B.C.; Sim, G.D. A 10-node tetrahedral element with condensed Lagrange multipliers for the modified couple stress theory. Comput. Struct. 2021, 246, 106476. [Google Scholar] [CrossRef]
  65. Grbčić, S.; Ibrahimbegović, A.; Jelenić, G. Variational formulation of micropolar elasticity using 3D hexahedral finite-element interpolation with incompatible modes. Comput. Struct. 2018, 205, 1–14. [Google Scholar] [CrossRef] [Green Version]
  66. Fantuzzi, N.; Trovalusci, P.; Dharasura, S. Mechanical Behavior of Anisotropic Composite Materials as Micropolar Continua. Front. Mater. 2019, 6, 59. [Google Scholar] [CrossRef]
  67. Fantuzzi, N.; Trovalusci, P.; Luciano, R. Material Symmetries in Homogenized Hexagonal-Shaped Composites as Cosserat Continua. Symmetry 2020, 12, 441. [Google Scholar] [CrossRef] [Green Version]
  68. Zhao, J.; Chen, W.J.; Lo, S.H. A refined nonconforming quadrilateral element for couple stress/strain gradient elasticity. Int. J. Numer. Methods Eng. 2011, 85, 269–288. [Google Scholar] [CrossRef]
  69. Yang, F.; Chong, A.; Lam, D.; Tong, P. Couple stress based strain gradient theory for elasticity. Int. J. Solids Struct. 2002, 39, 2731–2743. [Google Scholar] [CrossRef]
  70. Eremeyev, V.A.; Altenbach, H. On the Direct Approach in the Theory of Second Gradient Plates. In Shell and Membrane Theories in Mechanics and Biology: From Macro- to Nanoscale Structures; Altenbach, H., Mikhasev, G.I., Eds.; Springer International Publishing: Cham, Switzerland, 2015; pp. 147–154. [Google Scholar]
  71. Aifantis, E.C. Update on a class of gradient theories. Mech. Mater. 2003, 35, 259–280. [Google Scholar] [CrossRef]
  72. Askes, H.; Aifantis, E.C. Gradient elasticity in statics and dynamics: An overview of formulations, length scale identification procedures, finite element implementations and new results. Int. J. Solids Struct. 2011, 48, 1962–1990. [Google Scholar] [CrossRef]
  73. Cornacchia, F.; Fantuzzi, N.; Luciano, R.; Penna, R. Solution for cross-and angle-ply laminated Kirchhoff nano plates in bending using strain gradient theory. Compos. Part B Eng. 2019, 173, 107006. [Google Scholar] [CrossRef]
  74. Cornacchia, F.; Fabbrocino, F.; Fantuzzi, N.; Luciano, R.; Penna, R. Analytical solution of cross-and angle-ply nano plates with strain gradient theory for linear vibrations and buckling. Mech. Adv. Mater. Struct. 2019, 1–15. [Google Scholar] [CrossRef]
  75. Tocci Monaco, G.; Fantuzzi, N.; Fabbrocino, F.; Luciano, R. Hygro-thermal vibrations and buckling of laminated nanoplates via nonlocal strain gradient theory. Compos. Struct. 2020, 262, 113337. [Google Scholar] [CrossRef]
  76. Tocci Monaco, G.; Fantuzzi, N.; Fabbrocino, F.; Luciano, R. Critical Temperatures for Vibrations and Buckling of Magneto-Electro-Elastic Nonlocal Strain Gradient Plates. Nanomaterials 2021, 11, 87. [Google Scholar] [CrossRef] [PubMed]
  77. Bacciocchi, M.; Fantuzzi, N.; Ferreira, A.J.M. Conforming and nonconforming laminated finite element Kirchhoff nanoplates in bending using strain gradient theory. Comput. Struct. 2020, 239, 106322. [Google Scholar] [CrossRef]
  78. Bacciocchi, M.; Fantuzzi, N.; Ferreira, A.J.M. Static finite element analysis of thin laminated strain gradient nanoplates in hygro-thermal environment. Contin. Mech. Thermodyn. 2020, 1–24. [Google Scholar] [CrossRef]
  79. Falope, F.O.; Lanzoni, L.; Radi, E.; Tarantino, A.M. Thin film bonded to elastic orthotropic substrate under thermal loading. J. Strain Anal. Eng. Des. 2016, 51, 256–269. [Google Scholar] [CrossRef] [Green Version]
  80. Falope, F.O.; Radi, E. Finite thin cover on an orthotropic elastic half plane. Model. Simul. Eng. 2016, 2016. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Laminated composite plate: geometry and stacking sequence.
Figure 1. Laminated composite plate: geometry and stacking sequence.
Materials 14 01771 g001
Figure 2. Dimensionless central deflection w ¯ versus side-to-thickness ratio a / h of a square plate subjected to a sinusoidally distributed load varying the nonlocal ratio ( / a ) 2 , for different lamination schemes: (a) cross-ply ( 0 / 90 / 90 / 0 ) ; (b) angle-ply ( 45 / 45 ) . TSDT—Third-order Shear Deformation Theory; FSDT—First-order Shear Deformation Theory; CLPT—Classical Laminated Plate Theory.
Figure 2. Dimensionless central deflection w ¯ versus side-to-thickness ratio a / h of a square plate subjected to a sinusoidally distributed load varying the nonlocal ratio ( / a ) 2 , for different lamination schemes: (a) cross-ply ( 0 / 90 / 90 / 0 ) ; (b) angle-ply ( 45 / 45 ) . TSDT—Third-order Shear Deformation Theory; FSDT—First-order Shear Deformation Theory; CLPT—Classical Laminated Plate Theory.
Materials 14 01771 g002
Figure 3. Through-the-thickness distributions of stresses of a ( 0 / 90 / 90 / 0 ) square plate characterized by a / h = 4 for different theories and values of ( / a ) 2 . The following stress components are considered: (a) σ ¯ x x ; (b) σ ¯ y y ; (c) σ ¯ x z ; (d) σ ¯ y z ; (e) σ ¯ x y .
Figure 3. Through-the-thickness distributions of stresses of a ( 0 / 90 / 90 / 0 ) square plate characterized by a / h = 4 for different theories and values of ( / a ) 2 . The following stress components are considered: (a) σ ¯ x x ; (b) σ ¯ y y ; (c) σ ¯ x z ; (d) σ ¯ y z ; (e) σ ¯ x y .
Materials 14 01771 g003
Figure 4. Through-the-thickness distributions of stresses of a ( 0 / 90 / 90 / 0 ) square plate characterized by a / h = 10 , for different theories and values of ( / a ) 2 . The following stress components are considered: (a) σ ¯ x x ; (b) σ ¯ y y ; (c) σ ¯ x z ; (d) σ ¯ y z ; (e) σ ¯ x y .
Figure 4. Through-the-thickness distributions of stresses of a ( 0 / 90 / 90 / 0 ) square plate characterized by a / h = 10 , for different theories and values of ( / a ) 2 . The following stress components are considered: (a) σ ¯ x x ; (b) σ ¯ y y ; (c) σ ¯ x z ; (d) σ ¯ y z ; (e) σ ¯ x y .
Materials 14 01771 g004
Figure 5. Through-the-thickness distributions of stresses of a ( 45 / 45 ) square plate characterized by a / h = 4 for different theories and values of ( / a ) 2 . The following stress components are considered: (a) σ ¯ x x ; (b) σ ¯ y y ; (c) σ ¯ x z ; (d) σ ¯ y z ; (e) σ ¯ x y .
Figure 5. Through-the-thickness distributions of stresses of a ( 45 / 45 ) square plate characterized by a / h = 4 for different theories and values of ( / a ) 2 . The following stress components are considered: (a) σ ¯ x x ; (b) σ ¯ y y ; (c) σ ¯ x z ; (d) σ ¯ y z ; (e) σ ¯ x y .
Materials 14 01771 g005
Figure 6. Through-the-thickness distributions of stresses of a ( 45 / 45 ) square plate characterized by a / h = 10 , for different theories and values of ( / a ) 2 . The following stress components are considered: (a) σ ¯ x x ; (b) σ ¯ y y ; (c) σ ¯ x z ; (d) σ ¯ y z ; (e) σ ¯ x y .
Figure 6. Through-the-thickness distributions of stresses of a ( 45 / 45 ) square plate characterized by a / h = 10 , for different theories and values of ( / a ) 2 . The following stress components are considered: (a) σ ¯ x x ; (b) σ ¯ y y ; (c) σ ¯ x z ; (d) σ ¯ y z ; (e) σ ¯ x y .
Materials 14 01771 g006
Table 1. Dimensionless central displacement and stress components of a cross-ply square plate with ( 0 / 90 / 0 ) as the lamination scheme, for ( / a ) 2 = 0 .
Table 1. Dimensionless central displacement and stress components of a cross-ply square plate with ( 0 / 90 / 0 ) as the lamination scheme, for ( / a ) 2 = 0 .
a / h Theory 10 2 w ¯ σ ¯ xx σ ¯ yy σ ¯ yz σ ¯ xz σ ¯ xy
4TSDT [1]1.92180.7345-0.2086--
TSDT1.92170.73440.03150.20860.28550.0497
FSDT [1]1.77580.4370-0.1968--
FSDT1.77570.43700.03070.19680.33680.0369
10TSDT [1]0.71250.5684-0.1167--
TSDT0.71250.56840.01830.11670.36930.0277
FSDT [1]0.66930.5134-0.1108--
FSDT0.66930.51340.01760.11080.38060.0252
100TSDT [1]0.43420.5390-0.0827--
TSDT0.43420.53900.01340.08270.39480.0214
FSDT [1]0.43370.5384-0.0827--
FSDT0.43370.53840.01340.08270.39500.0213
CLPT [1]0.43130.5387-0.0823--
CLPT0.43120.53870.01330.08230.39510.0213
Table 2. Dimensionless central displacement and stress components of a cross-ply square plate with ( 0 / 90 / 0 ) as the lamination scheme, varying the nonlocal ratio ( / a ) 2 .
Table 2. Dimensionless central displacement and stress components of a cross-ply square plate with ( 0 / 90 / 0 ) as the lamination scheme, varying the nonlocal ratio ( / a ) 2 .
( / a ) 2 a / h Theory 10 2 w ¯ σ ¯ xx σ ¯ yy σ ¯ yz σ ¯ xz σ ¯ xy
0.054TSDT0.96720.36960.01580.10500.14370.0250
FSDT0.89370.21990.01540.09910.16950.0186
10TSDT0.35860.28600.00920.05870.18590.0139
FSDT0.33680.25840.00890.05580.19160.0127
100TSDT0.21850.27130.00670.04160.19870.0107
FSDT0.21830.27100.00670.04160.19880.0107
CLPT0.21700.27110.00670.04140.19890.0107
0.14TSDT0.64620.24700.01060.07010.09600.0167
FSDT0.59710.14690.01030.06620.11320.0124
10TSDT0.23960.19110.00620.03920.12420.0093
FSDT0.22510.17260.00590.03730.12800.0085
100TSDT0.14600.18120.00450.02780.13280.0072
FSDT0.14580.18100.00450.02780.13280.0072
CLPT0.14500.18110.00450.02770.13290.0072
Table 3. Dimensionless central displacement and stress components of a cross-ply square plate with ( 0 / 90 / 90 / 0 ) as the lamination scheme, for ( / a ) 2 = 0 .
Table 3. Dimensionless central displacement and stress components of a cross-ply square plate with ( 0 / 90 / 90 / 0 ) as the lamination scheme, for ( / a ) 2 = 0 .
a / h Theory 10 2 w ¯ σ ¯ xx σ ¯ yy σ ¯ yz σ ¯ xz σ ¯ xy
4TSDT [1]1.89400.66500.63200.23900.20600.0440
TSDT1.89360.66510.63220.29850.23050.0440
FSDT [1]1.71000.40600.57600.19600.14000.0308
FSDT1.70950.40590.57640.27990.26860.0308
10TSDT [1]0.71500.54600.38900.15300.26400.0268
TSDT0.71470.54560.38880.19240.30690.0268
FSDT [1]0.66300.49890.36100.13000.16700.0241
FSDT0.66270.49890.36140.18070.31810.0241
20TSDT [1]0.50600.53900.30400.12300.28200.0228
TSDT0.50600.53930.30430.15410.32990.0228
FSDT [1]0.49100.52700.29600.10900.17500.0221
FSDT0.49120.52730.29560.15030.33320.0221
100TSDT [1]0.43400.53900.27100.11200.29000.0213
TSDT0.43430.53870.27080.13890.33890.0213
FSDT [1]0.43400.53800.27000.10100.17800.0213
FSDT0.43370.53820.27040.13870.33900.0213
CLPT [1]0.43100.53900.26900.13800.33900.0213
CLPT0.43120.53870.26930.13820.33930.0213
Table 4. Dimensionless central displacement and stress components of a cross-ply square plate with ( 0 / 90 / 90 / 0 ) as the lamination scheme, varying the nonlocal ratio ( / a ) 2 .
Table 4. Dimensionless central displacement and stress components of a cross-ply square plate with ( 0 / 90 / 90 / 0 ) as the lamination scheme, varying the nonlocal ratio ( / a ) 2 .
( / a ) 2 a / h Theory 10 2 w ¯ σ ¯ xx σ ¯ yy σ ¯ yz σ ¯ xz σ ¯ xy
0.054TSDT0.95300.33470.31820.15020.11600.0222
FSDT0.86040.20430.29010.14090.13520.0155
10TSDT0.35970.27460.19570.09680.15450.0135
FSDT0.33350.25110.18190.09090.16010.0121
20TSDT0.25470.27140.15310.07760.16610.0115
FSDT0.24720.26540.14880.07570.16770.0111
100TSDT0.21860.27110.13630.06990.17050.0107
FSDT0.21830.27090.13610.06980.17060.0107
CLPT0.21700.27110.13560.06960.17070.0107
0.14TSDT0.63670.22360.21260.10040.07750.0148
FSDT0.57480.13650.19380.09410.09030.0104
10TSDT0.24030.18350.13070.06470.10320.0090
FSDT0.22280.16780.12150.06080.10700.0081
20TSDT0.17020.18130.10230.05180.11090.0077
FSDT0.16520.17730.09940.05060.11200.0074
100TSDT0.14600.18110.09110.04670.11390.0072
FSDT0.14580.18100.09090.04660.11400.0072
CLPT0.14500.18110.09060.04650.11410.0072
Table 5. Dimensionless central displacement and membrane stress components of several cross-ply square plates, varying the nonlocal ratio ( / a ) 2 . The reference solutions (Ref.) are taken from [77].
Table 5. Dimensionless central displacement and membrane stress components of several cross-ply square plates, varying the nonlocal ratio ( / a ) 2 . The reference solutions (Ref.) are taken from [77].
10 2 w ¯ σ ¯ xx σ ¯ yy σ ¯ xy
( / a ) 2 SchemeRef. [77]PresentRef. [77]PresentRef. [77]PresentRef. [77]Present
0.05 ( 0 ) 0.21700.21700.27110.27110.01340.01340.01070.0107
( 0 / 90 ) 0.53530.53520.04240.04240.36020.36020.02640.0264
( 0 / 90 ) 2 0.25490.25490.01800.01800.24500.24500.01260.0126
( 0 / 90 ) 4 0.22540.22540.01490.01490.24910.24910.01110.0111
0.1 ( 0 ) 0.14500.14500.18110.18110.00900.00900.00720.0072
( 0 / 90 ) 0.35760.35760.02840.02840.24070.24060.01760.0176
( 0 / 90 ) 2 0.17030.17030.01200.01200.16370.16370.00840.0084
( 0 / 90 ) 4 0.15060.15060.01000.01000.16640.16640.00740.0074
Table 6. Dimensionless central displacement 10 2 w ¯ of a angle-ply square plate with ( θ / θ / ) as lamination scheme, for ( / a ) 2 = 0 . The number of plies is denoted by N L .
Table 6. Dimensionless central displacement 10 2 w ¯ of a angle-ply square plate with ( θ / θ / ) as lamination scheme, for ( / a ) 2 = 0 . The number of plies is denoted by N L .
θ = 5 ° θ = 30 ° θ = 45 °
a / h Theory N L = 2 N L = 6 N L = 2 N L = 6 N L = 2 N L = 6
4TSDT [1]1.26251.22821.08380.88511.02030.8375
TSDT1.26241.22821.08360.88511.02020.8375
FSDT [1]1.31651.26471.21550.89941.15760.8531
FSDT1.31651.26471.21540.89941.15750.8531
10TSDT [1]0.48480.44850.59160.30070.55810.2745
TSDT0.48480.44850.59150.30060.55800.2745
FSDT [1]0.48830.44910.60990.29890.57730.2728
FSDT0.48820.44910.60980.29890.57720.2728
20TSDT [1]0.35790.32090.51800.21270.48970.1905
TSDT0.35790.32090.51790.21270.48960.1905
FSDT [1]0.35860.32080.52240.21210.49440.1899
FSDT0.35850.32080.52240.21210.49430.1899
100TSDT [1]0.31620.27890.49420.18420.46760.1634
TSDT0.31620.27890.49410.18420.46760.1634
FSDT [1]0.31620.27890.49440.18420.46780.1633
FSDT0.31620.27890.49430.18420.46780.1633
CLPT [1]0.31450.27710.49320.18310.46670.1622
CLPT0.31440.27710.49310.18310.46670.1622
Table 7. Dimensionless central displacement 10 2 w ¯ of a angle-ply square plate with ( θ / θ / ) as lamination scheme, varying the nonlocal ratio ( / a ) 2 . The number of plies is denoted by N L .
Table 7. Dimensionless central displacement 10 2 w ¯ of a angle-ply square plate with ( θ / θ / ) as lamination scheme, varying the nonlocal ratio ( / a ) 2 . The number of plies is denoted by N L .
θ = 5 ° θ = 30 ° θ = 45 °
( / a ) 2 a / h Theory N L = 2 N L = 6 N L = 2 N L = 6 N L = 2 N L = 6
0.054TSDT0.63530.61810.54540.44550.51340.4215
FSDT0.66250.63650.61170.45270.58250.4293
10TSDT0.24400.22570.29770.15130.28080.1382
FSDT0.24570.22600.30690.15040.29050.1373
20TSDT0.18010.16150.26060.10700.24640.0959
FSDT0.18040.16140.26290.10670.24880.0956
100TSDT0.15910.14040.24870.09270.23530.0822
FSDT0.15910.14030.24880.09270.23540.0822
CLPT0.15830.13950.24820.09210.23490.0816
0.14TSDT0.42450.41300.36440.29760.34310.2816
FSDT0.44270.42530.40870.30240.38920.2868
10TSDT0.16300.15080.19890.10110.18760.0923
FSDT0.16420.15100.20510.10050.19410.0917
20TSDT0.12030.10790.17410.07150.16460.0640
FSDT0.12060.10790.17560.07130.16620.0638
100TSDT0.10630.09380.16620.06200.15720.0549
FSDT0.10630.09380.16620.06190.15730.0549
CLPT0.10570.09320.16580.06160.15690.0546
Table 8. Dimensionless central displacement and stress components of a cross-ply square plate with ( 45 / 45 ) 4 as lamination scheme, for ( / a ) 2 = 0 .
Table 8. Dimensionless central displacement and stress components of a cross-ply square plate with ( 45 / 45 ) 4 as lamination scheme, for ( / a ) 2 = 0 .
a / h Theory 10 2 w ¯ σ ¯ xx σ ¯ xy σ ¯ xz
4TSDT1.29310.25110.24060.2116
FSDT1.33170.14450.13840.2487
10TSDT0.42070.16230.15540.2425
FSDT [1]0.41980.14450.13840.2487
FSDT0.41980.14450.13840.2487
20TSDT0.29000.14900.14270.2471
FSDT [1]0.28960.14450.13840.2487
FSDT0.28960.14450.13840.2487
100TSDT0.24790.14470.13860.2486
FSDT [1]0.24790.14450.13840.2487
FSDT0.24790.14450.13840.2487
CLPT [1]0.24620.14450.13840.2487
CLPT0.24610.14450.13840.2487
Table 9. Dimensionless central displacement and stress components of a cross-ply square plate with ( 45 / 45 ) 4 as lamination scheme, varying the nonlocal ratio ( / a ) 2 .
Table 9. Dimensionless central displacement and stress components of a cross-ply square plate with ( 45 / 45 ) 4 as lamination scheme, varying the nonlocal ratio ( / a ) 2 .
( / a ) 2 a / h Theory 10 2 w ¯ σ ¯ xx σ ¯ xy σ ¯ xz
0.054TSDT0.65080.12640.12110.1065
FSDT0.67020.07270.06970.1251
10TSDT0.21170.08170.07820.1220
FSDT0.21130.07270.06970.1251
20TSDT0.14600.07500.07180.1244
FSDT0.14570.07270.06970.1251
100TSDT0.12480.07280.06970.1251
FSDT0.12480.07270.06970.1251
CLPT0.12390.07270.06970.1251
0.14TSDT0.43480.08440.08090.0711
FSDT0.44780.04860.04650.0836
10TSDT0.14150.05460.05230.0815
FSDT0.14120.04860.04650.0836
20TSDT0.09750.05010.04800.0831
FSDT0.09740.04860.04650.0836
100TSDT0.08340.04860.04660.0836
FSDT0.08340.04860.04650.0836
CLPT0.08280.04860.04650.0836
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bacciocchi, M.; Tarantino, A.M. Third-Order Theory for the Bending Analysis of Laminated Thin and Thick Plates Including the Strain Gradient Effect. Materials 2021, 14, 1771. https://doi.org/10.3390/ma14071771

AMA Style

Bacciocchi M, Tarantino AM. Third-Order Theory for the Bending Analysis of Laminated Thin and Thick Plates Including the Strain Gradient Effect. Materials. 2021; 14(7):1771. https://doi.org/10.3390/ma14071771

Chicago/Turabian Style

Bacciocchi, Michele, and Angelo Marcello Tarantino. 2021. "Third-Order Theory for the Bending Analysis of Laminated Thin and Thick Plates Including the Strain Gradient Effect" Materials 14, no. 7: 1771. https://doi.org/10.3390/ma14071771

APA Style

Bacciocchi, M., & Tarantino, A. M. (2021). Third-Order Theory for the Bending Analysis of Laminated Thin and Thick Plates Including the Strain Gradient Effect. Materials, 14(7), 1771. https://doi.org/10.3390/ma14071771

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