Next Article in Journal
On the Microstructure and Properties of Nb-12Ti-18Si-6Ta-5Al-5Cr-2.5W-1Hf (at.%) Silicide-Based Alloys with Ge and Sn Additions
Previous Article in Journal
Influence of the Size and Type of Pores on Brick Resistance to Freeze-Thaw Cycles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geometric Non-Linear Analysis of Auxetic Hybrid Laminated Beams Containing CNT Reinforced Composite Materials

1
State Key Laboratory of Ocean Engineering, Shanghai Jiao Tong University, Shanghai 200240, China
2
Shanghai Key Laboratory for Digital Maintenance of Buildings and Infrastructure, School of Naval Architecture, Ocean and Civil Engineering, Shanghai Jiao Tong University, Shanghai 200240, China
3
School of Civil Engineering, University of Birmingham, Birmingham B15 2TT, UK
4
College of Civil Engineering, Nanjing Tech University (Nanjing Tech), Nanjing 211816, China
*
Author to whom correspondence should be addressed.
Materials 2020, 13(17), 3718; https://doi.org/10.3390/ma13173718
Submission received: 10 July 2020 / Revised: 12 August 2020 / Accepted: 17 August 2020 / Published: 22 August 2020
(This article belongs to the Section Construction and Building Materials)

Abstract

:
In the current work, a novel hybrid laminate with negative Poisson’s ratio (NPR) is developed by considering auxetic laminate which is composed of carbon nanotube-reinforced composite (CNTRC) and fiber-reinforced composite (FRC) materials. The maximum magnitude of out-of-plane NPR is identified in the case of (20 F/20 C/−20 C/20 C) S laminate as well. Meanwhile, a method for the geometric non-linear analysis of hybrid laminated beam with NPR including the non-linear bending, free, and forced vibrations is proposed. The beam deformation is modeled by combining higher-order shear-deformation theory (HSDT) and large deflection theory. Based on a two-step perturbation approach, the asymptotic solutions of the governing equations are obtained to capture the linear and non-linear frequencies and load-deflection curves. Moreover, a two-step perturbation methodology in conjunction with fourth-order Runge–Kutta method is employed to solve the forced-vibration problem. Several key factors, such as CNT distribution, variations in the elastic foundation, and thermal stress, are considered in the exhaustive analysis. Theoretical results for some particular cases are given to examine the geometric non-linearity behavior of hybrid beam with NPR as well as positive Poisson’s ratio (PPR).

1. Introduction

Materials and structures with negative Poisson’s ratio (NPR) behave in a counter-intuitive manner: when compressed (stretched) in the axial direction, they contract (expand) transversely. The materials and structures that exhibit this feature are also termed as “auxetics” [1]. Lakes [2] first reported NPR behavior in polyurethane (PU) foam with re-entrant structures. Wojciechowski [3] presented the first thermodynamically stable molecular model to study the mechanisms for generating auxetic behavior of solid. NPR was analyzed in systems containing rigid rotating hexamers [4]. Thereafter, analysis on cellular auxetics [5,6], multi-material auxetics [7,8], and auxetic composites [9,10] was carried out. A more systematic review of the development and application for the investigation of auxetic material and structures was reported by Ren et al. [11] and Lakes [12].
In recent years, laminated beams or plates with NPR are developed with more applications as primary structural elements in many fields. Fascinating properties of composite materials with NPR have led researchers to find outstanding applications in the field of aerospace [13,14], automobile industry [15,16], and civil engineering [17,18]. Recently, studies on the design of auxetic structures have propelled to achieve improvements in the mechanical properties such as impact resistance [19,20] and energy absorption [21,22]. On the other hand, with the growing applications of auxetic materials in different industries, developing multi-scale models for design and capturing the non-linear static and dynamic responses of the auxetic laminates under various loading scenarios is critical, when related to the structural design. A series of work has been dedicated to the study of composite beams, in particular, investigations have been carried out on their static and dynamic behavior [23,24,25,26]. However, the value of NPRs (ranging from −0.2 to −1) is not for any particular composite material in real engineering, but is a hypothetical value in a model analysis.
Fiber-reinforced composite (FRC) materials have attracted the attention of researchers, primarily due to the strong anisotropy they offer when designing laminate with NPR. Zhang et al. [27] showed that both the particular stacking sequence and the individual ply material (strongly anisotropic) are essential for a laminate to exhibit NPR. The authors also presented an optimal angle for ply and particular stacking sequences were presented. Evans et al. [28] specially designed a software to predict the effective engineering constants. It was reported that the NPR property can be obtained by designing stacking sequences in the laminated plates. Lempriere [29] measured that the effective Poisson’s ratio (EPR) in orthotropic materials is −0.4, occurring at θ = 45° orientation. Clarke et al. [30] reported that the EPR of the laminates show negative values for lay angle in a range between 15° and 30° for a ( ± θ ). It should be noted that when layers are oriented at angles equal in magnitude but with opposite signs, the appropriate + or – sign is used. Herakovich [31] investigated the auxetic characteristics of laminated structure made of graphite-epoxy to determine the value of v 13 e . Such laminates exhibited a very wide range of NPR ranging from a peak of 0.49 for a laminate with ply angle of 90° to as low as −0.21 for laminate with ply angle of (±25)S. Hine et al. [32] reported that the out-of-plane Poisson’s ratio reached −1/2 when a high modulus of elasticity carbon fiber is used in the laminates. Matsuda et al. [33] observed that peak values of NPR in carbon fiber-reinforced plastic laminates were around −0.7 when the axis is oriented at 25 o. The influence of Young’s modulus ratio (E1/E2), the type of resin and the volume of fraction on the ERP ( v 13 e ) of an angle-ply [±θ]2s plate was investigated by Harkati et al. [34,35]. It was shown that the NPR of Kevlar and carbon reinforced composite plate is −0.746 at θ = 20°. Therefore, it can be concluded form the above discussion that the maximum value of NPR of the laminates greatly depends on the ply orientations and stacking sequences in laminates [36,37,38]. One of the practical applications of such composite structures with NPR is the low-velocity impact resistance [39,40] and static indentation resistance [41]. Alderson and Coenen [39] found that the auxetic FRC laminates showed increased load and energy absorption to failure at low levels of impact energy. Zhou et al. [40] carried out experimental studies on the low-velocity impact response of 3D auxetic composites. They found that 3D auxetic textile composite exhibited excellent impact protective performance compared to 3D non-auxetic textile composite. Coenen and Alderson [41] manufactured auxetic laminates and evaluated their static indentation resistance in comparison with two laminates having near zero and larger positive Poisson’s ratio. The results showed enhancement in load sustained and energy absorption by auxetic laminates.
Carbon nanotubes (CNTs) are extensively used in different fields of industry and research with the development in fabrication technology [42,43,44]. A composite laminate structure made from CNTRCs is prone to exhibit an auxetic feature due to inherent special properties such as strong anisotropy. Furthermore, functionally graded (FG) materials are frequently employed in engineering applications and the laminae with various CNT volume fractions are used to achieve excellent mechanical properties [45]. Based on the studies mentioned above, the auxetic concept of FRC laminates is introduced for designing FRC/CNTRC hybrid laminate with NPR. This hybrid laminate is modeled by placing FRCs in the outermost layers and FG-CNTRC core. Shen was the first investigator to analyze the particular characteristics and behavior of the FG-CNTRC structures at different scales [45]. Following this landmark research, many studies on the forced vibration of FG-CNTRC beams [46,47] and plates [48] were carried out. The free vibration of FG-CNTRC plate were examined by Huu Quoc et al. [49] with a new refined plate theory. Based on Timoshenko beam theory (TBT), the non-linear vibration, forced vibration, and bending behavior of FG-CNTRC beams were performed by Mohammadimehr et al. [50] and Ansari et al. [51]. Considering the matrix cracks, Fan and Wang [52,53] presented an investigation on the non-linear static and dynamic responses of hybrid laminated structures made either from FRCs or CNTRCs. Based on an element-free numerical approach in conjunction with self-consistent model, the effect of matrix cracks on bending and vibration characteristics of hybrid laminated plates were reported by Lei et al. [54,55]. In addition to the applications of FRC/CNTRC hybrid beam or plate, the adoptions of the FRC/CNTRC hybrid concept were later extended to the design of composite blades by Zhang et al. [56,57]. They carried out several investigations exploring the non-linear vibration characteristics of hybrid composited blades. In addition, the design and analysis of the hybrid laminates comprising multi-materials with different strengths have been carried out by [58,59,60,61].
From the literature survey, it can be noted that no research has been carried out to analyze of hybrid beams using FRC/CNTRC materials with NPR at different external conditions. Therefore, the primary objective of the current work is to consider both hybrid configuration and NPR and analyze the static and dynamic characteristics of these auxetic structures. Furthermore, a non-linear model is developed based on the higher-order shear-deformation theory (HSDT) and the Von Kármán large deflection assumption. The motion equations due to thermal stress and reaction force due to foundation are derived. The effects of the CNT material distribution, environment condition, foundation type on the free and forced-vibration feature, and non-linear bending behavior of hybrid laminated beams are investigated in detail.

2. Design of Hybrid Laminate with NPR

As we know that both the ply orientations and stacking sequence have a significant effect on the mechanical response of laminated structures. Therefore, the effective engineering constants are usually used for the convenience of engineers in describing the mechanical behavior of the laminates. Sun et al. [36] and Chen et al. [62] presented a model for EPR for general thick laminates. However, only the extensional response was taken into consideration while the bending and bending-extension coupling characteristics of the laminates were neglected. Therefore, their model fails to provide accurate solutions for EPR of an asymmetric angle-ply laminated plate. Considering the effects of bending and bending-extension coupling, the general solutions of the EPRs for an arbitrary angle-ply laminate are derived in the previous works [63,64,65].
For a laminate where the material parameters of the layers are not distributed symmetrically along the section, the general formula of v 13 e can be expressed as follows:
ν 13 e = | A 13 B 11 B 5 - 3 D | / | A 11 B 11 B 5 - 1 D | ,
where the matrixes (A11, A13, B11, B5-1, B5-3, D) are defined in [63,64,65].
For symmetric laminates, the bending-extension coupling stiffnesses Bij (i,j = 1,2…6) are zero. The preceding expression simplifies as follows:
v 13 e = A 16 ( A 22 A 36 A 23 A 26 ) + A 13 ( A 26 2 A 22 A 66 ) + A 12 ( A 23 A 66 A 26 A 36 ) A 26 2 A 33 2 A 23 A 26 A 36 + A 23 2 A 66 + A 22 ( A 36 2 A 33 A 66 ) ,
where Aij represents the beam stiffnesses, which are defined in terms of the transformed elastic coefficients ( C ¯ i j )k as:
A i j = k = 1 N h k 1 h k ( C ¯ i j ) k d Z , ( i , j = 1 ~ 6 ) ,
[ C ¯ i j ] 1 = [ S ¯ i j ] ,
where the compliance constants ( S ¯ i j ) of a laminate whose fibers direction makes an angle θ with the direction of X-axis (see Figure 1) and c = cosθ, s = sinθ.
[ S ¯ 11 S ¯ 12 S ¯ 22 S ¯ 16 S ¯ 26 S ¯ 66 ] = [ c 4 2 c 2 s 2 s 4 c 2 s 2 c 2 s 2 c 4 + s 4 c 2 s 2 c 2 s 2 s 4 2 c 2 s 2 c 4 c 2 s 2 2 c 3 s 2 ( c s 3 c 3 s ) 2 c s 3 c s ( c 2 s 2 ) 2 c s 3 2 ( c 3 s c s 3 ) 2 c 3 s c s ( c 2 s 2 ) 4 c 2 s 2 8 c 2 s 2 4 c 2 s 2 ( c 2 s 2 ) 2 ]   [ S 11 S 12 S 22 S 66 ] ,
[ S ¯ 13 S ¯ 23 S ¯ 33 ] = [ c 2 s 2 0 s 2 c 2 0 0 0 1 ] [ S 13 S 23 S 33 ] ,
where the compliance constants of laminate Sij are given as follow:
[ S 11 S 12 S 13 S 22 S 23 S 33 S 44 S 55 S 66 ] = [ 1 / E 11 ν 12 / E 11 ν 13 / E 11 1 / E 22 ν 23 / E 22 1 / E 33 1 / G 23 1 / G 13 1 / G 12 ]
where the basic material parameters of each layer of laminate are introduced as follow by referring to Figure 1.
Eii, (I = 1,2,3) = Young’s moduli in i, (i=1,2,3) directions
Gij, (ij = 12, 13, 23) = Shear moduli in i-j planes, respectively
vij, ij = 12, 13, 23) = Poisson’s ratios (the subscripts i and j represent the loading and strain directions, respectively).
The theoretical solution presented above can predict the EPR of an arbitrary angle-ply laminates. A systematic investigation of the symmetric hybrid laminates has been carried out. The hybrid laminated beam is designed by placing FRC in the outermost layers and CNTRC in the rest of the layers which gives excellent performance. The lay-ups of the hybrid laminated beam are considered eight-layered and angle-ply (θF/θC/-θC/θC)S. The subscript/index “S” indicates that the laminate is symmetric and θ is the angle of CNT or fiber orientation. Unless otherwise stated, the superscript F represents the layer of FRC while the superscript C represents the other layers with CNTRC. These configurations are for an identical material and have constant thickness of 0.125 mm and 0.25 mm for FRC and CNT layers, respectively. Several types of CNT distributions are taken into consideration and the type of CNT volume fraction (VCN) is given. The temperature-related material properties of CNTRCs are predicted by the extended micromechanical model [45] as summarized in Table 1. It should be noted that the thermal expansion coefficients (α11, α22) are calculated by Equation (15), which is obtained by the Schapery model [66].
To assess how the distribution and volume fractions of CNT influence the static and dynamic behavior of the beam, we considered five configurations summarized in Table 2. The volume fraction of fiber is fixed as Vf = 0.6. It differ in the following characteristics: FG-Λ: 0.6/(0.11)2/(0.14)2/(0.17)2/0.6, FG-V: 0.6/(0.17)2/(0.14)2/(0.11)2/0.6, FG-X:(0.6/0.11/0.14/0.17)S, FG-O: (0.6/0.17/0.14/0.11)S, and uniform distribution (UD): (0.6/0.14/0.14/0.14)S. It can be assumed that G13 = G23 = G12.
The EPR ( v 13 e ) of the hybrid FRC/CNTRC laminated beam is obtained by adjusting specific ply orientations and stacking sequence. The results in terms of NPR and PPR of hybrid laminated beam and the corresponding ply orientations are presented in Table 3.

3. Theoretical Modeling of Hybrid Beams

Laminates consist of layers of composites reinforced with CNTRC and FRC. Consider a hybrid laminated beam composed of eight layers with lamination scheme (θF/θC/-θC/θC)S resting on a continuous visco-Pasternak foundation as shown in Figure 1. Figure 1a defines the coordinate system used in development of the hybrid laminated beam analysis. The XYZ coordinate system is assumed to have its origin on the middle face of the beam so that the middle surface lies in the XY-plane. Q(X,t) is the out-of-plane static or dynamic load. The displacement at a point on the X, Y, and Z directions are U ¯ , V ¯ , and W ¯ , respectively.
The simply supported beam is resting on a three-parameter foundation including the Winkler foundation ( K ¯ 1 ), shearing layer stiffness ( K ¯ 2 ), and damping parameter C ¯ d . The reaction force from the visco-Pasternak foundation P0 (X, t ¯ ) is given by:
p 0 ( X , t ¯ ) = K ¯ 1 W ¯ ( X , t ¯ ) K ¯ 2 2 W ¯ ( X , t ¯ ) X 2 + C ¯ d 2 W ¯ ( X , t ¯ ) t ¯ 2
The method of analysis is based on the HSDT [67] for the laminated beam undergoing large deflection. The effect of the elevated temperature is considered by introducing thermal stress resultants N ¯ T ,   M ¯ T , and P ¯ T as shown in Appendix A. The motion equations are given as follow:
S 11 4 W ¯ X 4 + S 12 3 Ψ ¯ x X 3 + B ¯ 11 A ¯ 11 2 N ¯ T X 2 + 2 M ¯ T X 2 + N ¯ x 2 W ¯ X 2 + Q ( X , t ¯ )
= p 0 ( X , t ¯ ) + I 1 2 W ¯ t ¯ 2 + I ^ 5 3 Ψ ¯ x X t ¯ 2 4 3 h 2 I ^ 7 4 W ¯ X 2 t ¯ 2
S 21 3 W ¯ X 3 + S 22 2 Ψ ¯ x X 2 S 23 ( W ¯ X + Ψ ¯ x ) S 26 N ¯ T X + S ¯ T X = I ˜ 3 2 Ψ ¯ x t ¯ 2 4 I ˜ 5 3 h 2 3 W ¯ X t ¯ 2
N ¯ x = L 1 0 L { [ A ¯ 11 2 ( W ¯ X ) 2 4 E ¯ 11 3 h 2 ( 2 W ¯ X 2 + Ψ ¯ x X ) + B ¯ 11 Ψ ¯ x X ] N ¯ T } d X
where Ψ ¯ x denotes rotation about the longitudinal axes. The reduced stiffnesses of the beam ( A ¯ 11 , B ¯ 11 ,   E ¯ 11 ) and the coefficients Sij and inertias Ii are defined and given in Appendix A.
The formulae for the forces owing to thermal stress are given as:
( N ¯ T , M ¯ T , P ¯ T ) = k = 1 N b t k 1 t k [ A x ] k ( 1 , Z , Z 3 ) Δ T   d Z
S ¯ T = M ¯ T 4 P ¯ T 3 h 2
where ΔT is the temperature increment from a reference state (T0 = 300 K), ΔT = TT0. T is set as 400 K or 500 K.
The coefficient Ax is known in terms of the thermal expansion coefficients (α11, α22):
A x = Q ¯ 11 ( c 2 α 11 + s 2 α 22 ) + Q ¯ 12 ( s 2 α 11 + c 2 α 22 ) + Q ¯ 16 2 c s ( α 11 α 22 )
α i j = { ( V C N E 11 C N α 22 C N + V m E m α m ) / ( V C N E 11 C N + V m E m )   i j = 11   ( 1 + ν 12 C N ) V C N α 22 C N + ( 1 + ν m ) V m α m ν 12 α 11   i j = 22
where αij (ij = 11,22) are the thermal expansion coefficients and subscripts 11 and 22 denote the longitudinal and transverse directions, respectively. The formulae for calculating the α11 and α22 are obtained from [66].
It should be noted that Equation (11) represent restricted boundary conditions. It is suitable for a beam with immovable boundary condition i.e., the longitudinal displacement is equal to zero at both ends of the beam. It should be noted that this immovable boundary condition is unacceptable for the compressive post-buckling analysis of beam.
The non-linear motion equations for the vibration can be solved by a two-step perturbation approach proposed by Shen [68,69]. By deriving the dimensionless forms of the motion Equations (9)–(11), we get:
γ 11 4 W x 4 γ 12 3 Ψ x x 3 π { 0 π [ γ 13 2 ( W x ) 2 + γ 14 Ψ x x γ 15 2 W x 2 ] d x } 2 W x 2 + C 1 2 W x 2 γ 16 2 N T x 2 C 2 2 M T x 2 = λ q ( K 1 W K 2 2 W x 2 + C d W t ) + γ 17 2 W t 2 + γ 18 3 Ψ x x t 2 + γ 19 4 W x 2 t 2
γ 21 3 W x 3 γ 22 2 Ψ x x 2 + γ 23 ( W x + Ψ x ) γ 26 N T x C 3 S T x = γ 28 2 Ψ x t 2 + γ 29 3 W x t 2
where
( C 1 , C 2 , C 3 ) = ( γ T 1 , γ T 3 , γ T 3 γ T 6 ) Δ T
The non-dimensional parameters mentioned above can be written as follow:
( x , W , Ψ x ) = ( π X L , W ¯ L , Ψ ¯ x π ) ,   ( N x , M x ,   P x ) = L 2 π 2 h D ¯ 11 ( h N ¯ x , M ¯ x ,   4 3 h 2 P ¯ x ) , ( K 1 , K 2 ) = L 2 π 2 D ¯ 11 ( L 2 π 2 K ¯ 1 ,   K ¯ 2 ) ,   ( k 1 , k 2 ) = L 2 E 0 I ( L 2 K ¯ 1 ,   K ¯ 2 ) , ( γ 11 ,   γ 12 ,   γ 21 ,   γ 22 ) = 1 D ¯ 11 ( S 11 ,   S 12 ,   S 21 ,   S 22 ) ,   γ 23 = L 2 π 2 D ¯ 11 S 23 , ( γ 13 ,   γ 14 ,   γ 15 ) = L π 2 D ¯ 11 ( L A ¯ 11 , B ¯ 11 4 3 h 2 E ¯ 11 ,   4 3 h 2 E ¯ 11 ) ,   ( γ 16 ,   γ 26 ) = 1 A ¯ 11 L ( B ¯ 11 ,   B ¯ 11 4 3 h 2 E ¯ 11 ) , ( γ 17 , γ 18 , γ 19 , γ 28 ,   γ 29 ) =   ( I 1 L 2 π 2 I ^ 5 ,   4 3 h 2 I ^ 7 ,   I ˜ 3 ,   4 3 h 2 I ˜ 5 ) b E 0 ρ 0 D ¯ 11 ,   t = π   t ¯ L E 0 ρ 0 ,   γ T 1 = L 2 A x T π 2 D ¯ 11 , ω L = Ω L L π ρ 0 E 0 ,   ( γ T 3 ,   γ T 6 ) = L 2 π 2 h D ¯ 11 ( D x T ,   4 3 h 2 F x T ) ,   λ q = Q L 3 π 4 D ¯ 11 ,   C d = C ¯ d a 3 π 3 D 11 * E 0 ρ 0
in which E0 = Em and ρ0 = ρm are the elastic modulus and density of the matrix. A x T , D x T , F x T are the functions of thickness and Ax can be defined by:
( A x T ,   D x T ,   F x T ) = b k = 1 N t k 1 t k [ A x ] k ( 1 , Z , Z 3 ) d Z

4. Free Vibration Analysis

For non-linear vibration problem, the solutions for Equations (16) and (17) consist of an additional displacement term and initial displacement term as a result of the varying temperature. The initial deflection under the thermal loading for these structures are reported from Shen [63]. By considering, τ = εt, the solution equations can be expanded as a function with a small perturbation parameter εj (j = 1,2,3,….) as given below:
Ψ x ( x , τ , ε ) = ε ψ x 1 ( x , τ ) + ε 2 ψ x 2 ( x , τ ) + ε 3 ψ x 3 ( x , τ ) + Ψ x ( x , τ , ε ) = ε ψ x 1 ( x , τ ) + ε 2 ψ x 2 ( x , τ ) + ε 3 ψ x 3 ( x , τ ) + λ q ( x , τ , ε ) = ε λ q 1 ( x , τ ) + ε 2 λ q 2 ( x , τ ) + ε 3 λ q 3 ( x , τ ) +
The boundary conditions for τ = 0 can be expressed as:
W | τ = 0 = Ψ x | τ = 0 = 0 ,   W τ | τ = 0 = Ψ x τ | τ = 0 = 0 , ( simply   supported ) .
If we substitute τ = εt and Equation (21) into Equations (16) and (17) and collect all the terms of the same order of εi (i = 1,2,3,…), we will get a set of perturbation equations. The solution of Euler-Bernoulli beam with simply supported conditions can be used for perturbation equations with εi (i = 1)
w 1 ( x , τ ) = A 10 ( 1 ) ( τ ) sin m x
Then, the equations with εj (j = 1,2,3,...) can be solved step by step, we will get:
W ( x , t , ε ) = ε A 10 ( 1 ) ( t ) sin m x + O ( ε 4 )
Ψ x ( x , t , ε ) = ε B 10 ( 1 ) ( t ) cos m x + ε 3 B 10 ( 3 ) ( t ) cos m x + O ( ε 4 )
λ q ( x , t , ε ) = ( ε A ¨ 10 ( 1 ) ( t ) g 30 + ε A 10 ( 1 ) ( t ) g 31 ) sin m x + ( ( ε A 10 ( 1 ) ( t ) ) 2 g 321 + ( ε A 10 ( 1 ) ) ( ε A 10 * ) g 322 ) sin m x + ( ε A 10 ( 1 ) ( t ) ) 3 g 33 sin m x + O ( ε 4 )
Introducing expression τ = t and ε A 10 ( 1 ) =Wm in Equation (26) and applying Galerkin procedure yielded Equation (26) which can be re-written as:
g 30 d 2 ( W m ) d t 2 + g c d ( W m ) d t + g 31 ( W m ) + g 32 ( W m ) 2 + g 33 ( W m ) 3 λ ¯ q ( t ) = 0
The numerical solution of Equation (27) corresponds to the forced-vibration response of the structure. It can be converted to a free vibration problem by setting the uniform load to zero. Therefore, the solution of the differential equation with g c = λ ¯ q (t) = 0 can be obtained as follow:
ω N L = ω L ( 1 + 9 g 31 g 33 10 g 32 2 12 g 31 2 A 2 ) 1 / 2
ω L = ( g 31 / g 30 ) 1 / 2
where ωNL and ωL represent the dimensionless non-linear and linear frequency, respectively. A = Wm = W ¯ m /L. W ¯ m is the amplitude of deflection. According to Equation (19), the corresponding linear frequency can be expressed as Ω = ωL(π/L)(E0/ρ0)1/2. The details of g30, g31, g32, and g33 are described in Appendix B.
For the purpose of verification, Table 4 and Figure 2 show the comparison of the solutions obtained by the present method and the method proposed by Fan & Wang [70] with the present solutions of a hybrid laminated beam. The reinforced material of the hybrid beam consists of two parts: (1) for FRC, E 11 F = 233.05 GPa, v 12 f = 0.2,   E 22 F = 23.1 GPa, ρf = 1750 kg/m³, G 12 F = 8.96 GPa, (2) for CNTRC, E 11 C N = 5646.6 GPa, v 12 C N = 0.175,   E 22 C N = 7080.0 GPa, ρCN = 1400 kg/m³ and   G 12 C N = 1944.5 GPa. The material properties for the same matrix are Em = 2.5 GPa, ρm = 1150 kg/m3, and νm = 0.34. The value of geometric parameters used are: L = 20 h, h = 0.5mm.
Table 4 compares the fundamental frequencies Ω ˜   = Ω(L2/h) ρ 0 / E 0 for (θC/90F)S hybrid beam with different values of VCN and the angle of the CNT orientation with θC = 15°, 30°, and 45°. It can be found that the predicted frequencies tally reasonably well with those results obtained from Fan & Wang [70]. Meanwhile, the non-linear vibration of a (0C/90F)S hybrid laminated beam is validated with Fan &Wang study, as plotted in Figure 2. It can be observed from Figure 2 that the results in terms of frequency ratio obtained from the presented model are very close to the predictions made by Fan & Wang [70]. Based on the above discussion, it can be concluded that the proposed model is accurate to analyze free vibration of the hybrid laminated beams.
In the following, a detailed investigation of the vibration of hybrid laminated beams has been carried out. The effect of temperature T, distribution pattern, and the foundation constants (k1, k2) on the free vibration of hybrid laminated beams are scrutinized. In the following study, the angle-ply (25 F/20 C/-20 C/20 C) S and (25 F/90 C/-90 C/90 C) S hybrid laminated beams are adopted.
The effect of temperature T on the first four frequencies Ω ˜ i (i = 1, 2, 3, 4) of hybrid laminated beams with five different distribution pattern is investigated in Table 5. It is observed that increase in temperature cause decrease in the frequencies of both the angle-ply hybrid laminated beams. In addition, the temperature effect of the first-order frequency is more significant than on the other order frequencies. Consequently, the influence of degradation caused by thermal stress should be taken into consideration in the design and service of hybrid laminated beams.
Table 6 shows the effects of foundation constants on the hybrid laminated beam with length/thickness ratio of 20. Three different foundation coefficients (k1, k2) are (103, 10), (103, 0), and (0, 0), respectively. It should be noted that structures with (k1, k2) = (0, 0) are selected for comparison. As can be seen from Table 6, the natural frequencies of the hybrid laminated beam for higher (k1, k2) are more than the other cases.
The frequency ratio–deflection curves for all the five types of hybrid laminated beams are described in Figure 1. As expected, FG-X has the highest frequency ratio while FG-O has the lowest frequency ratio. At the room temperature, FG-Λ and FG-V show similar higher amplitude frequency (see Figure 3).
The influence of temperature field on the non-linear vibration behavior of UD and FG-X hybrid laminated beam with NPR is examined by considering three different temperature values (T = 300 K,400 K,500 K). Due to the presence of thermal stress, increase in the temperature leads to a higher frequency ratio, as shown in Figure 4. It is found that (25 F/20 C/−20 C/20 C) S laminated beams are more sensitive to the temperature effect than those with (25 F/90 C−90 C/90 C) S.
The effect of the foundation coefficients on the resulting lager amplitude vibration frequency of UD and FG-X are examined in Figure 5, for a member with L = 10 h. From Figure 5, it is confirmed that the beam becomes more rigid as the foundation coefficients of the beam increases. Hence, the frequency ratio of the UD and FG-X increases.

5. Forced-Vibration Analysis

In this section, the dynamic response of hybrid laminated beams is investigated by considering uniform load Q which is dependent on time as plotted in Figure 1. We need to determine the relationship between central deflection and time. Hence, a numerical procedure for solving the second order differential Equation (27) is employed here. Given the initial value Wm (t0) and W ˜ ˙ m   (t0) at initial time t0 = 0, Equation (27) can be solved to obtain the central deflection-time relationship for the beam under applied load by employing the fourth-order Runge–Kutta method. Here, it should also be noted that the initial deflection for FG-Λ or FG-V will be triggered because of thermal stress.
In the following, the time response history of the verification analyses is presented. Figure 6 shows the comparison of the solution of a cross-ply (0/90/0/90/0)S graphene-reinforced composite (GRC) laminated beam under a sudden load Q = 12 Mpa by the method presented by Fan et al. [71]. and the method presented in this study. The dimensions and material properties of the GRC layer are as follows: L = 20 h = 60 mm, E 11 G = 1.812 TPa, E 22 G = 1.87 TPa, G 12 G = 0.683 TPa, v 12 G = 0.177, ρG = 4118 kg/m³, Em = 2.5 GPa, ρm = 1150 kg/m³. vm = 0.34. The subscripts, G and m, represent the graphene and matrix, respectively. In the forced-vibration stage (time range 0 to 0.2 ms), the central deflection-time curve is described in Figure 6. It is found from Figure 6 that the forced response predicted from the proposed model follows the same trend and agrees well with that provided by Fan et al. [71].
Let us examine the forced-vibration behavior by applying the proposed method to the cases (25 F/20 C/−20 C/20 C) S. Figure 7 exhibits the forced-vibration curves of the uniform distribution (UD) beam under the reference temperature (T = 300 K). Four different dynamic loads are taken into consideration. The values of the parameters used are L = 20 h, h = 1.75 mm. In the forced-vibration region (time range 0 to 0.8ms), it can be observed that the step load leads to highest deflection among the four while lowest deflection occurs at the exponential load. At time range of 0.8 to 1.6ms, the free vibration behavior of hybrid laminated beams with an initial deflection and velocity is triggered resulting of the dynamic load is released. Please note that a sudden uniform load is taken as the applied load in the parametric study.
Figure 8 depicts central deflection versus time for UD and FG-CNTRC beam under the reference temperature (T = 300 K). Five different distribution patterns designed above are considered. The value of load amplitude is taken as Q = 0.4 MPa. Figure 8 shows deflection-time curves among the above designed distributing patterns. The curves of the UD pattern are used at the same temperature for reference. It can be observed that the lowest deflection is obtained by FG-Λ beam as the UD plate has the highest deflection. Thus, for the next parametric analysis, UD and FG-Λ are considered/studied.
The above approach is also used to investigate the influence of changing the temperature field on UD and FG-Λ beams. Figure 9 shows the influence of the thermal stress on UD and FG-Λ beams. The applied load is 45 MPa for (25 F/20 C/−20 C/20 C) S and (25 F/90 C/−90 C/90 C) S with L = 5 h. Incremental temperatures are marked on these curves. The curves show that the amplitude and period of response increase with increase in temperature. This is again attributed to the fact that the development of thermal resultants reduces the overall plate stiffness.
Figure 10 shows the maximum deflection with time variation of hybrid beams (25 F/20 C/−20 C/20 C) S and (25 F/90 C/−90 C/90 C)S for different foundations stiffnesses ((k1, k2 Cd) = (0, 0, 0), (102, 0, 0), (102, 10, 1), (102, 10,2)). Beams without foundation (k1 = k2 = Cd = 0) are selected as a reference. The results confirm that the foundation stiffness decrease the transverse deflection. Moreover, the deflection of beams on the viscoelastic foundation decreases with increasing time.

6. Non-Linear Bending Analysis

In this section, non-linear bending response of hybrid laminated beams with NPR and PPR will be discussed in detail. For this purpose, we will determine relationship between applied pressure and deflection of the beam. For the static analysis, this relationship is independent of the change in time (t). Therefore, the transverse load is modified to be uniformly distributed, and Q (X, t) = Q (X) = q0. W is independent of time. Under these assumptions, Equations (16) and (17) can be simplified as:
γ 11 4 W x 4 γ 12 3 Ψ x x 3 π { 0 π [ γ 13 2 ( W x ) 2 + γ 14 Ψ x x γ 15 2 W x 2 ] d x } 2 W x 2 + C 1 2 W x 2
γ 16 2 N T x 2 C 2 2 M T x 2 = λ q ( K 1 W K 2 2 W x 2 )
γ 21 3 W x 3 γ 22 2 Ψ x x 2 + γ 23 ( W x + Ψ x ) γ 26 N T x C 3 S T x = 0
A perturbation approach is also used to derive the solutions of Equations (30) and (31). The solution equations can be expanded as a function with a small perturbation parameter εj (j = 1,2,3,…).
W ( x , ε ) = j = 1 ε j w j ( x ) Ψ x ( x , ε ) = j = 1 ε j ψ x j ( x ) λ q ( x , ε ) = j = 1 ε j λ j ( x )
Following the perturbation solutions procedure, the solutions for the equations with the first order ε can be assumed as:
w 1 ( x ) = A 10 ( 1 ) sin m x
The thermal load can be expanded as the Fourier modes as:
[ M x T S x T ] = k = 1 , 3 , 1 k sin k x [ M x ( 1 ) S x ( 1 ) ]
( M x ( 1 ) , S x ( 1 ) ) = 4 π ( γ T 3 , γ T 3 γ T 6 ) Δ T
For the free vibration analysis, the half-wavelength (m) can be assumed as m = 1,2,3…, whereas, for the bending study m is adopted as 1. Substituting Equations (32) and (35) in the motion equations, we get the asymptotic solution as:
W ( x , ε ) = ε A 10 ( 1 ) sin m x + O ( ε 4 )
Ψ ˜ x ( x , ε ) = ε B 10 ( 1 ) cos m x + O ( ε 4 )
and
λ q ( x , ε ) = λ q 1 ( A 10 ( 1 ) ε ) + λ q 2 ( A 10 ( 1 ) ε ) 2 + λ q 3 ( A 10 ( 1 ) ε ) 3 + O ( ε 4 )
Equations (36) and (37) are the functions of A 10 ( 1 ) . In Equations (36) and (37), (ε A 10 ( 1 ) ) is taken as the second perturbation parameter related to the dimensionless maximum deflection Wm. From Equations (36a) and (37), the relationships between the load and central displacement are obtained as:
q 0 L 3 π 4 D ¯ 11 = A W ( 0 ) + A W ( 1 ) ( W ¯ m L ) + A W ( 2 ) ( W ¯ m L ) 2 + A W ( 3 ) ( W ¯ m L ) 3 +
in which
A W ( 0 ) =   [ γ T 3 ( γ T 3 γ T 6 ) γ 12 γ 22 + γ 23 ] Δ T , A w ( 1 ) = π 4 [ ( γ 11 γ 12 γ 21 γ 23 γ 22 + γ 23 ) + ( K 1 + K 2 ) γ T 1 Δ T ] , A W ( 2 ) = π 2 2 ( γ 15 γ 14 γ 21 γ 23 γ 22 + γ 23 ) , A W ( 3 ) = π 3 16 γ 13 ,
The non-linear bending of a hybrid laminated beam with cross-ply (0 C/90 F/0 C/90 F/0 C) under a uniformly distributed load is studied in the current analysis. Both the ends of the beam are taken as simply supported. The same analysis has been reported by Fan & Wang [72]. Constituent materials of the beam are similar to the cases of hybrid beam listed in Table 4. Present dimensionless deflections are compared with those of Fan & Wang [72] for the hybrid laminated beam with L/h = 10, 20, 40. For the sake of consistency, we have used qoL3/E0I, which is adopted by Fan & Wang [72]. Nevertheless, qoL4/E0I is non-dimensional load which is adopted in this paper. According to curves plotted in Figure 11, similar trend and good accordance may be observed between the results.
Figure 12 illustrates the dimensionless load-deflection curves for UD and FG beam under a uniform pressure. Five different distribution patterns designed above are considered. W ¯ m /h and qoL4/E0I represent the dimensionless defection and load, respectively. Figure 12 shows the deflection-time curves among the above-mentioned distributing patterns. The curve of the UD pattern is used at the same temperature for reference. UD beam has the maximum central deflection as compared with the other distribution patterns while FG-Λ beam has the minimum central displacement. Thus, for the next comprehensive studies, UD and FG-Λ are taken as the case studies.
Figure 13 demonstrates the influence of the changing temperature field on the bending response of UD and FG-Λ with the above-mentioned conditions and L = 5 h. It may be concluded from Figure 13 that the central deflection increases with a rise in temperature. It is also observed that (25 F/20 C/−20 C/20 C) S laminated beams are less sensitive to variations in the temperature than (25 F/90 C/−90 C/90 C) S laminated beam.
Figure 14 plots the effect of the foundation stiffnesses on the bending characteristics of UD and FG-Λ with reference environmental conditions and L = 20 h. The predicted load-deflection curves are given for (k1, k2) = (0, 0), (103, 0), (103, 102). In Figure 14, it is evident that the dimensionless central deflection decreases by increasing the foundation stiffnesses. Meanwhile, the effect of reinforcement of foundation on (25 F/90 C/−90 C/90 C) S laminated beam is more obvious than the other beams.

7. Conclusions

In this study, a model for free vibration, non-linear forced vibration, and bending analyses of hybrid laminated beams with NPR and PPR under various external conditions is proposed. The model developed accounts for the FG configurations and NPRs in hybrid laminated beams. Five types of volume fractions of CNTs are considered which include: UD, FG-Λ, FG-V, and FG-X. For a configuration of (25 F/20 C/−20 C/20 C) S and (25 F/90 C/−90 C/90 C) S, the out-of-plane Poisson’s ratios exhibit NPR and PPR, respectively. The solutions for the static and dynamic responses are derived by using a two-step perturbation approach by using a fourth-order Runge–Kutta method. The theoretical model predicts reliable values of frequencies and deflections in comparison with the results available in the existing model in the literature. The influence of the temperature variations and elastic foundation on the static and dynamic behavior of hybrid laminated beams is investigated in detail.
Following comprehensive results are concluded:
  • The distribution type of CNTs show considerable effect on the static and dynamic behavior of FRC/CNTRC plate with NPR or PPR. It is concluded that the FG-X beams have a lower frequency ratio. Moreover, for the non-linear bending and dynamic response, the FG-Λ beam showed better performance than that with UD under different external conditions.
  • An increment in the temperature make the frequencies ratio and central defection considerably larger, whereas the increasing foundation stiffness will result in an opposite effect.
  • The dynamic responses considering viscosity of foundation show obvious differences from those obtained with Pasternak foundation.
  • The comparative studies reveal that the non-linear vibration of (25 F/20 C/−20 C/20 C) S laminated beams are more sensitive to the temperature field, whereas, the non-linear bending in the (25 F/20 C/−20 C/20 C) S laminated beams are less sensitive to the changes in temperature and foundation coefficient than those with (25 F/90 C/−90 C/90 C) S.
The results of the current work demonstrate that the NPR has a significant effect on the non-linear bending, free vibration, and forced-vibration characteristics of FRC/CNTRC hybrid laminated beams. Furthermore, the results show that the thermal-mechanical characteristics of FRC/CNTRC hybrid laminated beam can be remarkably improved by the determination of proper distribution of CNTs. For example, to have satisfactory performance in high temperature environments from hybrid laminated beams, beams with FG-X or FG-Λ arrangement can be preferred. Meanwhile, it is hoped that the current work will give an insight into the non-linear behavior of hybrid laminated beams with NPR and will be helpful for the further investigations.

Author Contributions

Conceptualization: X.-h.H., J.Y. Methodology: X.-h.H., J.Y. Writing review and editing, X.-h.H., J.Y., I.A. Writing original draft X.-h.H., I.A. Project administration J.Y., X.-e.W., X.R. Investigation X.-e.W., X.R. Validation X.-e.W., X.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Science Research Plan of Shanghai Municipal Science and Technology Committee, grant number 18DZ1205603, National Natural Science Foundation of China [Grant No. 51908352, 51978330, 51808286] and Natural Science Foundation of Jiangsu Province [Grant No. BK20180710].

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In Equations (9) and (10)
S 11 = 4 3 h 2 ( B ¯ 11 E ¯ 11 A ¯ 11 F ¯ 11 ) ,   S 12 = D ¯ 11 4 F ¯ 11 3 h 2 + B ¯ 11 A ¯ 11 ( 4 3 h 2 E ¯ 11 B ¯ 11 ) , S 21 = 4 3 h 2 [ 4 3 h 2 H ¯ 11 + E ¯ 11 A ¯ 11 ( B ¯ 11 4 3 h 2 E ¯ 11 ) F ¯ 11 ] , S 22 = D ¯ 11 4 F ¯ 11 3 h 2 + 4 3 h 2 ( 4 H ¯ 11 3 h 2 F ¯ 11 ) 1 A ¯ 11 ( B ¯ 11 4 E ¯ 11 3 h 2 ) 2 , S 23 = A ¯ 55 4 D ¯ 55 h 2 4 h 2 ( D ¯ 55 4 F ¯ 55 h 2 ) ,   S 26 = 1 A ¯ 11 ( B ¯ 11 4 E ¯ 11 3 h 2 ) ,
( A ¯ 11 ,   B ¯ 11 ,   D ¯ 11 ,   E ¯ 11 ,   F ¯ 11 ,   H ¯ 11 ) = k = 1 N b h k 1 h k ( Q ˜ 11 ) k ( 1 , Z , Z 2 , Z 3 , Z 4 , Z 6 ) d Z
( A ¯ 55 , D ¯ 55 , F ¯ 55 ) = k = 1 N b h k 1 h k ( Q ˜ 55 ) k ( 1 , Z 2 , Z 4 ) d Z
in which the detailed definition of Q ˜ 11 and Q ˜ 55 are found in [73]
Q ˜ 11 = Q ¯ 11 + ( Q ¯ 16 ) 2 Q ¯ 22 2 Q ¯ 12 Q ¯ 16 Q ¯ 26 + ( Q ¯ 12 ) 2 Q ¯ 66 ( Q ¯ 26 ) 2 Q ¯ 22 Q ¯ 66 , Q ˜ 55 = Q ¯ 55 ( Q ¯ 45 ) 2 Q ¯ 44 ,
and ( Q ¯ i j ) k are the components of the transformed stiffness of kth lamina with orientation θ ( k ) , defined by
[ Q ¯ 11 Q ¯ 12 Q ¯ 22 Q ¯ 16 Q ¯ 26 Q ¯ 66 ] = [ c 4 2 c 2 s 2 s 4 4 c 2 s 2 c 2 s 2 c 4 + s 4 c 2 s 2 4 c 2 s 2 s 4 2 c 2 s 2 c 4 4 c 2 s 2 c 3 s c s 3 c 3 s c s 3 2 c s ( c 2 s 2 ) c s 3 c 3 s c s 3 c 3 s 2 c s ( c 2 s 2 ) c 2 s 2 2 c 2 s 2 c 2 s 2 ( c 2 s 2 ) 2 ]   [ Q 11 Q 12 Q 22 Q 66 ]
[ Q ¯ 44 Q ¯ 45 Q ¯ 55 ] = [ c 2 c s s 2 s 2 c s c 2 ] [ Q 44 Q 55 ]
where
Q 11 = E 11 ( 1 v 12 v 21 ) 1 ,   Q 22 = E 22 ( 1 v 12 v 21 ) 1 ,   Q 12 = v 21 E 11 ( 1 v 12 v 21 ) 1 Q 44 = G 23 ,   Q 55 = G 13 ,   Q 66 = G 12
E11, E22, v12 etc… are summarized in Table 1.
In Equations (9) and (10), the coefficient Ii (i = 1,2,3,4,5,7) can be calculated as:
( I 1 , I 2 , I 3 , I 4 , I 5 , I 7 ) = k = 1 N ρ k h k 1 h k ( 1 , Z , Z 2 , Z 3 , Z 4 , Z 6 ) d Z
and
I ¯ 2 = I 2 4 I 4 3 h 2 ,   I ¯ 3 = I 3 8 3 h 2 ( I 5 2 3 h 2 I 7 ) ,   I ˜ 3 = I ¯ 3 I ¯ 2 I ¯ 2 I 1 ,   I ¯ 5 = I 5 4 I 7 3 h 2 , I ˜ 5 = I ¯ 5 I ¯ 2 I 4 I 1 ,   I ^ 5 = I ˜ 3 + 4 I ˜ 5 3 h 2 ,   I ˜ 7 = I 7 I 4 I 4 I 1 ,   I ^ 7 = I ˜ 5 + 4 I ˜ 7 3 h 2 ,

Appendix B

In Equation (27)
g 30 = γ 17 + m 2 ( γ 19 + γ 18 ( γ 21 m 2 γ 23 ) γ 22 m 2 + γ 23 ) + γ 21 m 4 γ 22 m 2 + γ 23 ( γ 29 + γ 28 ( γ 21 m 2 γ 23 ) γ 22 m 2 + γ 23 ) , g 31 = m 4 ( γ 11 γ 12 ( γ 21 m 2 γ 23 ) γ 22 m 2 + γ 23 ) + ( K 1 + K 2 m 2 ) γ T 1 m 2 Δ T + 2 π m ( γ 15 γ 14 ( γ 21 m 2 γ 23 ) γ 22 m 2 + γ 23 ) Φ + 2 π m 2 ( γ 15 γ 14 ( γ 21 γ 23 ) γ 22 + γ 23 ) Φ + 3 π 2 8 m 2 γ 13 Φ 2 , g 32 = { 2 π m 3 ( γ 15 γ 14 γ 21 m 2 γ 23 γ 22 m 2 + γ 23 ) π 2 4 3 γ 13 Φ   if   m = 1 2 π m 3 ( γ 15 γ 14 γ 21 m 2 γ 23 γ 22 m 2 + γ 23 ) π 2 4 m 2 γ 13 Φ   if   m 1 g 33 = π 2 m 4 γ 13 4 ,
in which
Φ = ε A 10 * = ( λ * G 10 ) G 20 G 10 ( λ * G 10 ) 2 C 33 ( λ * G 10 ) 3 +
where
λ * = 4 π ( γ T 3 γ 12 ( γ T 3 γ T 6 ) γ 22 + γ 23 ) Δ T , G 10 = ( γ 11 γ 12 ( γ 21 γ 23 ) γ 22 + γ 23 ) + ( K 1 + K 2 ) + 2 π γ 14 ( γ T 3 γ T 6 ) Δ T γ 22 + γ 23 , G 20 = 2 π ( γ 15 γ 14 ( γ 21 γ 23 ) γ 22 + γ 23 ) ,   G 30 = π 2 γ 13 4 ,   C 33 = G 30 G 10 ( 2 G 20 G 10 ) 2

References

  1. Evans, K.E.; Nkansah, M.A.; Hutchinson, I.J. Molecular Network Design. Nature 1991, 353, 124. [Google Scholar] [CrossRef]
  2. Lakes, R. Foam Structures with a Negative Poisson’s Ratio. Science 1987, 235, 1038–1040. [Google Scholar] [CrossRef] [PubMed]
  3. Wojciechowski, K.W. Constant Thermodynamic Tension Monte Carlo Studies of Elastic Properties of a Two-Dimensional System of Hard Cyclic Hexamers. Mol. Phys. 1987, 61, 1247–1258. [Google Scholar] [CrossRef]
  4. Wojciechowski, K.W. Two-Dimensional Isotropic System with a Negative Poisson Ratio. Phys. Lett. A 1989, 137, 60–64. [Google Scholar] [CrossRef]
  5. Grima, J.N.; Evans, K.E. Auxetic Behavior from Rotating Squares. J. Mater. Sci. Lett. 2000, 19, 1563–1565. [Google Scholar] [CrossRef]
  6. Hoover, W.G.; Hoover, C.G. Searching for Auxetics with DYNA3D and ParaDyn. Phys. Status Solidi B 2005, 242, 585–594. [Google Scholar] [CrossRef] [Green Version]
  7. Hiller, J.; Lipson, H. Tunable Digital Material Properties for 3D Voxel Printers. Rapid Prototyping J. 2010, 16, 241–247. [Google Scholar] [CrossRef] [Green Version]
  8. Wang, K.; Chang, Y.; Chen, Y.; Zhang, C.; Wang, B. Designable Dual-Material Auxetic Metamaterials Using Three-Dimensional Printing. Mater. Des. 2015, 67, 159–164. [Google Scholar] [CrossRef]
  9. Chen, Y.J.; Scarpa, F.; Farrow, I.R.; Liu, Y.J.; Leng, J.S. Composite Flexible Skin with Large Negative Poisson’s Ratio Range: Numerical and Experimental Analysis. Smart Mater. Struct. 2013, 22, 45005. [Google Scholar] [CrossRef]
  10. Fan, Y.; Xiang, Y.; Shen, H.-S. Temperature-Dependent Mechanical Properties of Graphene/Cu Nanocomposites with In-Plane Negative Poisson’s Ratios. Research 2020, 2020, 5618021. [Google Scholar] [CrossRef] [Green Version]
  11. Ren, X.; Das, R.; Tran, P.; Ngo, T.D.; Xie, Y.M. Auxetic Metamaterials and Structures: A Review. Smart Mater. Struct. 2018, 27, 023001. [Google Scholar] [CrossRef]
  12. Lakes, R.S. Negative-Poisson’s-Ratio Materials: Auxetic Solids. Annu. Rev. Mater. Res. 2017, 47, 63–81. [Google Scholar] [CrossRef]
  13. Siguier, J.; Mandeville, J. Test Procedures to Evaluate Spacecraft Materials Ejecta upon Hypervelocity Impact. Proc. Inst. Mech. Eng. Part. G-J. Aerosp. Eng. 2007, 221, 969–974. [Google Scholar] [CrossRef]
  14. Yu, Y.; Shen, H.-S. A Comparison of Nonlinear Vibration and Bending of Hybrid CNTRC/metal Laminated Plates with Positive and Negative Poisson’s Ratios. Int. J. Mech. Sci. 2020, 183, 105790. [Google Scholar] [CrossRef]
  15. Ju, J.; Kim, D.; Kim, K. Flexible Cellular Solid Spokes of a Non-Pneumatic Tire. Compos. Struct. 2012, 94, 2285–2295. [Google Scholar] [CrossRef]
  16. Wu, T.; Li, M.; Zhu, X.; Lu, X. Research on Non-Pneumatic Tire with Gradient Anti-Tetrachiral Structures. Mech. Adv. Mater. Struc. 2020. [Google Scholar] [CrossRef]
  17. Zahra, T.; Dhanasekar, M. Characterisation of Cementitious Polymer Mortar—Auxetic Foam Composites. Constr. Build. Mater. 2017, 147, 143–159. [Google Scholar] [CrossRef] [Green Version]
  18. Asad, M.; Dhanasekar, M.; Zahra, T.; Thambiratnam, D. Characterisation of Polymer Cement Mortar Composites Containing Carbon Fibre Or Auxetic Fabric Overlays and Inserts Under Flexure. Constr. Build. Mater. 2019, 224, 863–879. [Google Scholar] [CrossRef]
  19. Qiao, J.X.; Chen, C.Q. Impact Resistance of Uniform and Functionally Graded Auxetic Double Arrowhead Honeycombs. Int. J. Impact Eng. 2015, 83, 47–58. [Google Scholar] [CrossRef]
  20. Jiang, L.; Hu, H. Low-Velocity Impact Response of Multilayer Orthogonal Structural Composite with Auxetic Effect. Compos. Struct. 2017, 169, 62–68. [Google Scholar] [CrossRef]
  21. Chen, Z.; Wang, Z.; Zhou, S.; Shao, J.; Wu, X. Novel Negative Poisson’s Ratio Lattice Structures with Enhanced Stiffness and Energy Absorption Capacity. Materials 2018, 11, 1095. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Jiang, L.; Hu, H. Finite Element Modeling of Multilayer Orthogonal Auxetic Composites under Low-Velocity Impact. Materials 2017, 10, 908. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Chen, X.; Feng, Z. Dynamic Behaviour of a Thin Laminated Plate Embedded with Auxetic Layers Subject to In-Plane Excitation. Mech. Res. Commun. 2017, 85, 45–52. [Google Scholar] [CrossRef]
  24. Mansouri, M.H.; Shariyat, M. Biaxial Thermo-Mechanical Buckling of Orthotropic Auxetic FGM Plates with Temperature and Moisture Dependent Material Properties on Elastic Foundations. Compos. Pt. B-Eng. 2015, 83, 88–104. [Google Scholar] [CrossRef]
  25. Asemi, K.; Shariyat, M. Three-Dimensional Biaxial Post-Buckling Analysis of Heterogeneous Auxetic Rectangular Plates on Elastic Foundations by New Criteria. Comput. Method. Appl. M. 2016, 302, 1–26. [Google Scholar] [CrossRef]
  26. Lim, T. Shear Deformation in Beams with Negative Poisson’s Ratio. Proc. Inst. Mech. Eng. Part L J. Mater. Des. Appl. 2014, 229, 447–454. [Google Scholar] [CrossRef]
  27. Zhang, R.G.; Yeh, H.-L.; Yeh, H.-Y. A Preliminary Study of Negative Poisson’s Ratio of Laminated Fiber Reinforced Composites. J. Reinf. Plast. Comp. 1998, 17, 1651–1664. [Google Scholar] [CrossRef]
  28. Evans, K.E.; Donoghue, J.P.; Alderson, K.L. The Design, Matching and Manufacture of Auxetic Carbon Fibre Laminates. J. Compos. Mater. 2004, 38, 95–106. [Google Scholar] [CrossRef]
  29. Lempriere, M.B. Poisson’s Ratio in Orthotropic Materials. Aiaa J. 1968, 6, 2226–2227. [Google Scholar] [CrossRef]
  30. Clarke, J.F.; Duckett, R.A.; Hine, P.J.; Hutchinson, I.J.; Ward, I.M. Negative Poisson’s Ratios in Angle-Ply Laminates: Theory and Experiment. Composites 1994, 25, 863–868. [Google Scholar] [CrossRef]
  31. Herakovich, C.T. Composite Laminates with Negative Through-the-Thickness. J. Compos. Mater. 1984, 18, 447–455. [Google Scholar] [CrossRef] [Green Version]
  32. Hine, P.J.; Duckett, R.A.; Ward, I.M. Negative Poisson’s Ratios in Angle-Ply Laminates. J. Mater. Sci. Lett. 1997, 16, 541–544. [Google Scholar] [CrossRef]
  33. Matsuda, T.; Goto, K.; Kubota, N.; Ohno, N. Negative Through-The-Thickness Poisson’s Ratio of Elastic–Viscoplastic Angle-Ply Carbon Fiber-Reinforced Plastic Laminates: Homogenization Analysis. Int. J. Plast. 2014, 63, 152–169. [Google Scholar] [CrossRef]
  34. Hadi Harkati, E.; Bezazi, A.; Scarpa, F.; Alderson, K.; Alderson, A. Modelling the Influence of the Orientation and Fibre Reinforcement on the Negative Poisson’s Ratio in Composite Laminates. Phys. Status Solidi B 2007, 244, 883–892. [Google Scholar] [CrossRef]
  35. Hadi Harkati, E.; Bezazi, A.; Boukharouba, W.; Scarpa, F. Influence of Carbon Fibre on the Through-The-Thickness NPR Behaviour of Composite Laminates. Phys. Status Solidi B 2009, 246, 2111–2117. [Google Scholar] [CrossRef]
  36. Sun, C.T.; Li, S.J. Three-Dimensional Effective Elastic Constants for Thick Laminates. J. Compos. Mater. 1988, 22, 629–639. [Google Scholar] [CrossRef]
  37. Yeh, H.-L.; Yeh, H.-Y. Dilatation and Through-the-Thickness Poisson’s Ratio of Composite Laminates. J. Reinf. Plast. Comp. 2000, 19, 966–991. [Google Scholar] [CrossRef]
  38. Peel, L.D. Exploration of High and Negative Poisson’s Ratio Elastomer-Matrix Laminates. Phys. Status Solidi B 2007, 244, 988–1003. [Google Scholar] [CrossRef]
  39. Alderson, K.L.; Coenen, V.L. The Low Velocity Impact Response of Auxetic Carbon Fibre Laminates. Phys. Status Solidi B 2008, 245, 489–496. [Google Scholar] [CrossRef]
  40. Zhou, L.; Zeng, J.; Jiang, L.; Hu, H. Low-Velocity Impact Properties of 3D Auxetic Textile Composite. J. Mater. Sci. 2018, 53, 3899–3914. [Google Scholar] [CrossRef]
  41. Coenen, V.L.; Alderson, K.L. Mechanisms of Failure in the Static Indentation Resistance of Auxetic Carbon Fibre Laminates. Phys. Status Solidi B 2011, 248, 66–72. [Google Scholar] [CrossRef]
  42. Wang, R.; Xie, L.; Hameed, S.; Wang, C.; Ying, Y. Mechanisms and Applications of Carbon Nanotubes in Terahertz Devices: A Review. Carbon 2018, 132, 42–58. [Google Scholar] [CrossRef]
  43. Beitollahi, H.; Movahedifar, F.; Tajik, S.; Jahani, S. A Review on the Effects of Introducing CNTs in the Modification Process of Electrochemical Sensors. Electroanalysis 2018, 31, 1195–1203. [Google Scholar] [CrossRef]
  44. De Menezes, B.R.C.; Rodrigues, K.F.; da Silva Fonseca, B.C.; Ribas, R.G.; do Amaral Montanheiro, T.L.; Thim, G.P. Recent Advances in the Use of Carbon Nanotubes as Smart Biomaterials. J. Mater. Chen. B 2019, 7, 1343–1360. [Google Scholar] [CrossRef]
  45. Shen, H.-S. Nonlinear Bending of Functionally Graded Carbon Nanotube-Reinforced Composite Plates in Thermal Environments. Compos. Struct. 2009, 91, 9–19. [Google Scholar] [CrossRef]
  46. Asadi, H.; Beheshti, A.R. On the Nonlinear Dynamic Responses of FG-CNTRC Beams Exposed to Aerothermal Loads Using Third-Order Piston Theory. Acta Mech. 2018, 229, 2413–2430. [Google Scholar] [CrossRef]
  47. Mohammadimehr, M.; Monajemi, A.A.; Afshari, H. Free and Forced Vibration Analysis of Viscoelastic Damped FG-CNT Reinforced Micro Composite Beams. Microsyst. Technol. 2017. [Google Scholar] [CrossRef]
  48. Fu, T.; Chen, Z.; Yu, H.; Wang, Z.; Liu, X. Mechanical Behavior of Laminated Functionally Graded Carbon Nanotube Reinforced Composite Plates Resting on Elastic Foundations in Thermal Environments. J. Compos. Mater. 2018, 53, 1159–1179. [Google Scholar] [CrossRef]
  49. Huu Quoc, T.; Minh Tu, T.; Van Tham, V. Free Vibration Analysis of Smart Laminated Functionally Graded CNT Reinforced Composite Plates via New Four-Variable Refined Plate Theory. Materials 2019, 12, 3675. [Google Scholar] [CrossRef] [Green Version]
  50. Mohammadimehr, M.; Shahedi, S.; Rousta Navi, B. Nonlinear Vibration Analysis of FG-CNTRC Sandwich Timoshenko Beam Based on Modified Couple Stress Theory Subjected to Longitudinal Magnetic Field Using Generalized Differential Quadrature Method. Proc. Inst. Mech. Eng. Part. C J. Eng. Mech. Eng. Sci. 2016, 231, 3866–3885. [Google Scholar] [CrossRef]
  51. Ansari, R.; Faghih Shojaei, M.; Mohammadi, V.; Gholami, R.; Sadeghi, F. Nonlinear Forced Vibration Analysis of Functionally Graded Carbon Nanotube-Reinforced Composite Timoshenko Beams. Compos. Struct. 2014, 113, 316–327. [Google Scholar] [CrossRef]
  52. Fan, Y.; Wang, H. Nonlinear Vibration of Matrix Cracked Laminated Beams Containing Carbon Nanotube Reinforced Composite Layers in Thermal Environments. Compos. Struct. 2015, 124, 35–43. [Google Scholar] [CrossRef]
  53. Fan, Y.; Wang, H. Nonlinear Low-Velocity Impact on Damped and Matrix-Cracked Hybrid Laminated Beams Containing Carbon Nanotube Reinforced Composite Layers. Nonlinear Dynam. 2017, 89, 1863–1876. [Google Scholar] [CrossRef]
  54. Lei, Z.X.; Yin, B.B.; Liew, K.M. Bending and Vibration Behaviors of Matrix Cracked Hybrid Laminated Plates Containing CNTR-FG Layers and FRC Layers. Compos. Struct. 2018, 184, 314–326. [Google Scholar] [CrossRef]
  55. Lei, Z.X.; Zhang, L.W.; Liew, K.M. Modeling Large Amplitude Vibration of Matrix Cracked Hybrid Laminated Plates Containing CNTR-FG Layers. Appl. Math. Model. 2018, 55, 33–48. [Google Scholar] [CrossRef]
  56. Pan, Z.Z.; Chen, X.H.; Zhang, L.W. Modeling Large Amplitude Vibration of Pretwisted Hybrid Composite Blades Containing CNTRC Layers and Matrix Cracked FRC Layers. Appl. Math. Model. 2020, 83, 640–659. [Google Scholar] [CrossRef]
  57. Zhang, L.W.; Pan, Z.Z.; Chen, X.H. Vibration Characteristics of Matrix Cracked Pretwisted Hybrid Composite Blades Containing CNTRC Layers. J. Sound Vib. 2020, 473, 115–135. [Google Scholar] [CrossRef]
  58. Kuo, S. Vibration of Thermally Post-Buckled Hybrid Laminates with Two Non-Uniformly Distributed Fibers. Thin Wall. Struct. 2017, 114, 70–79. [Google Scholar] [CrossRef]
  59. Papa, I.; Boccarusso, L.; Langella, A.; Lopresto, V. Carbon/Glass Hybrid Composite Laminates in Vinylester Resin: Bending and Low Velocity Impact Tests. Compos. Struct. 2020, 232, 111571. [Google Scholar] [CrossRef]
  60. Asadi, H.; Bodaghi, M.; Shakeri, M.; Aghdam, M.M. An Analytical Approach for Nonlinear Vibration and Thermal Stability of Shape Memory Alloy Hybrid Laminated Composite Beams. Eur. J. Mech. A-Solids 2013, 42, 454–468. [Google Scholar] [CrossRef]
  61. Zhang, J.; Chaisombat, K.; He, S.; Wang, C.H. Hybrid Composite Laminates Reinforced with Glass/Carbon Woven Fabrics for Lightweight Load Bearing Structures. Mater. Des. 2012, 36, 75–80. [Google Scholar] [CrossRef]
  62. Chen, J.H.; Tsai, S.W. Three-Dimensional Effective Moduli of Symmetric Laminates. J. Compos. Mater. 1996, 30, 906–907. [Google Scholar] [CrossRef]
  63. Yang, J.; Huang, X.-H.; Shen, H.-S. Nonlinear Flexural Behavior of Temperature-Dependent FG-CNTRC Laminated Beams with Negative Poisson’s Ratio Resting on the Pasternak Foundation. Eng. Struct. 2020, 207, 110250. [Google Scholar] [CrossRef]
  64. Yang, J.; Huang, X.-H.; Shen, H.-S. Nonlinear Vibration of Temperature-Dependent FG-CNTRC Laminated Plates with Negative Poisson’s Ratio. Thin Wall. Struct. 2020, 148, 106514. [Google Scholar] [CrossRef]
  65. Huang, X.-H.; Yang, J.; Wang, X.-E.; Azim, I. Combined Analytical and Numerical Approach for Auxetic FG-CNTRC Plate Subjected to a Sudden Load. Eng. Comput. 2020. [Google Scholar] [CrossRef]
  66. Schapery, R.A. Thermal Expansion Coefficients of Composite Materials Based on Energy Principles. J. Compos. Mater. 1968, 2, 380–404. [Google Scholar] [CrossRef]
  67. Reddy, J.N. A Simple Higher-Order Theory for Laminated Composite Plates. J. Appl. Mech. 1984, 51, 745–752. [Google Scholar] [CrossRef]
  68. Shen, H.-S. A Two-Step Perturbation Method in Nonlinear Analysis of Beams, Plates and Shells; John Wiley & Sons Inc.: Singapore, 2013. [Google Scholar]
  69. Yang, J.; Huang, X.-H.; Shen, H.-S. Nonlinear Vibration of Temperature-Dependent FG-CNTRC Laminated Beams with Negative Poisson’s Ratio. Int. J. Struct. Stab. Dy. 2020, 20, 2050043. [Google Scholar] [CrossRef]
  70. Fan, Y.; Wang, H. The Effects of Matrix Cracks on the Nonlinear Vibration Characteristics of Shear Deformable Laminated Beams Containing Carbon Nanotube Reinforced Composite Layers. Int. J. Mech. Sci. 2017, 124–125, 216–228. [Google Scholar] [CrossRef]
  71. Fan, Y.; Xiang, Y.; Shen, H.-S. Nonlinear Dynamics of Temperature-Dependent FG-GRC Laminated Beams Resting on Visco-Pasternak Foundations. Int. J. Struct. Stab. Dyn. 2020, 20, 2050012. [Google Scholar] [CrossRef]
  72. Fan, Y.; Wang, H. The Effects of Matrix Cracks on the Nonlinear Bending and Thermal Postbuckling of Shear Deformable Laminated Beams Containing Carbon Nanotube Reinforced Composite Layers and Piezoelectric Fiber Reinforced Composite Layers. Compos. Pt. B Eng. 2016, 106, 28–41. [Google Scholar] [CrossRef]
  73. Li, J.; Wu, Z.; Kong, X.; Li, X.; Wu, W. Comparison of Various Shear Deformation Theories for Free Vibration of Laminated Composite Beams with General Lay-Ups. Compos. Struct. 2014, 108, 767–778. [Google Scholar] [CrossRef]
Figure 1. Various types of hybrid laminated beams and reference system: (a) Geometry and coordinate system of a hybrid laminated beam on a visco-elastic foundation; (b) Four types of CNT distribution in the cross section of a hybrid laminated beam.
Figure 1. Various types of hybrid laminated beams and reference system: (a) Geometry and coordinate system of a hybrid laminated beam on a visco-elastic foundation; (b) Four types of CNT distribution in the cross section of a hybrid laminated beam.
Materials 13 03718 g001
Figure 2. Comparison of ωNL/ωL for angle-ply (0 C/90 F) S hybrid laminated beam.
Figure 2. Comparison of ωNL/ωL for angle-ply (0 C/90 F) S hybrid laminated beam.
Materials 13 03718 g002
Figure 3. Frequency ratio (ωNL/ωL) versus dimensionless deflection ( W ¯ / h ) of hybrid laminated beams: (a) (25 F/20 C/−20 C/20 C) S; (b) (25 F/90 C/−90 C/90 C) S.
Figure 3. Frequency ratio (ωNL/ωL) versus dimensionless deflection ( W ¯ / h ) of hybrid laminated beams: (a) (25 F/20 C/−20 C/20 C) S; (b) (25 F/90 C/−90 C/90 C) S.
Materials 13 03718 g003
Figure 4. Frequency ratio (ωNL/ωL) versus dimensionless deflection ( W ¯ / h ) of hybrid beams with various temperature conditions (T = 300 K, 400 K, 500 K): (a) (25 F/20 C/−20 C/20 C)S; (b) (25 F/90 C/−90 C/90 C) S.
Figure 4. Frequency ratio (ωNL/ωL) versus dimensionless deflection ( W ¯ / h ) of hybrid beams with various temperature conditions (T = 300 K, 400 K, 500 K): (a) (25 F/20 C/−20 C/20 C)S; (b) (25 F/90 C/−90 C/90 C) S.
Materials 13 03718 g004
Figure 5. Frequency ratio versus dimensionless deflection of hybrid beams with various elastic foundation constants at reference temperature (T = 300 K): (a) (25 F/20 C/−20 C/20 C) S; (b) (25 F/90 C/−90 C/90 C)S.
Figure 5. Frequency ratio versus dimensionless deflection of hybrid beams with various elastic foundation constants at reference temperature (T = 300 K): (a) (25 F/20 C/−20 C/20 C) S; (b) (25 F/90 C/−90 C/90 C)S.
Materials 13 03718 g005
Figure 6. Comparison of forced-vibration curves for angle-ply (0/90/0/90/0)S GRC laminated beam under a sudden uniform load.
Figure 6. Comparison of forced-vibration curves for angle-ply (0/90/0/90/0)S GRC laminated beam under a sudden uniform load.
Materials 13 03718 g006
Figure 7. Central deflection versus time of hybrid beams under different dynamic loading scenarios (T = 300 K): (a) Type of dynamic loads; (b) Vibration curve.
Figure 7. Central deflection versus time of hybrid beams under different dynamic loading scenarios (T = 300 K): (a) Type of dynamic loads; (b) Vibration curve.
Materials 13 03718 g007
Figure 8. Central deflection versus time (t) of hybrid laminated beams: (a) (25 F/20 C/−20 C/20 C) S; (b) (25 F/90 C/−90 C/90 C) S.
Figure 8. Central deflection versus time (t) of hybrid laminated beams: (a) (25 F/20 C/−20 C/20 C) S; (b) (25 F/90 C/−90 C/90 C) S.
Materials 13 03718 g008
Figure 9. Deflection versus time (t) of hybrid laminated beams with various temperature conditions (T = 300 K, 350 K, 400 K): (a)(25 F/20 C/−20 C/20 C) S; (b) (25 F/90 C/−90 C/90 C) S.
Figure 9. Deflection versus time (t) of hybrid laminated beams with various temperature conditions (T = 300 K, 350 K, 400 K): (a)(25 F/20 C/−20 C/20 C) S; (b) (25 F/90 C/−90 C/90 C) S.
Materials 13 03718 g009
Figure 10. Deflection versus time (t) of hybrid laminated beams with various foundation constants at the reference temperature (T = 300 K): (a) (25 F/20 C/−20 C/20 C) S; (b) (25 F/90 C/−90 C/90 C) S.
Figure 10. Deflection versus time (t) of hybrid laminated beams with various foundation constants at the reference temperature (T = 300 K): (a) (25 F/20 C/−20 C/20 C) S; (b) (25 F/90 C/−90 C/90 C) S.
Materials 13 03718 g010
Figure 11. Comparisons of load-deflection relationships of a hybrid laminated beam.
Figure 11. Comparisons of load-deflection relationships of a hybrid laminated beam.
Materials 13 03718 g011
Figure 12. Load-deflection relationships of hybrid laminated beams: (a) (25 F/20 C/−20 C/20 C) S; (b) (25 F/90 C/−90 C/90 C) S.
Figure 12. Load-deflection relationships of hybrid laminated beams: (a) (25 F/20 C/−20 C/20 C) S; (b) (25 F/90 C/−90 C/90 C) S.
Materials 13 03718 g012
Figure 13. Load-deflections relationships of hybrid laminated beams at various temperatures: (a) (25 F/20 C/−20 C/20 C) S; (b) (25 F/90 C/−90 C/90 C) S.
Figure 13. Load-deflections relationships of hybrid laminated beams at various temperatures: (a) (25 F/20 C/−20 C/20 C) S; (b) (25 F/90 C/−90 C/90 C) S.
Materials 13 03718 g013
Figure 14. Load-deflection relationships of hybrid laminated beams with various foundation constants: (a) (25 F/20 C/−20 C/20 C) S; (b) (25 F/90 C/−90 C/90 C) S.
Figure 14. Load-deflection relationships of hybrid laminated beams with various foundation constants: (a) (25 F/20 C/−20 C/20 C) S; (b) (25 F/90 C/−90 C/90 C) S.
Materials 13 03718 g014
Table 1. Mechanical parameters with different temperature for each ply with CNTRC or FRC.
Table 1. Mechanical parameters with different temperature for each ply with CNTRC or FRC.
Volume FractionTE11E22G12v12α11α22
(K)(MPa)(MPa)(MPa) (×10−6/K)(×10−5/K)
VCN = 0.11300 K94,416.772203.74822.280.32193.58305.3182
400 K92,708.581710.53638.2530.32194.25145.5640
500 K91,682.211217.32454.2220.32194.61235.8198
VCN = 0.14300 K120,384.602297.68857.330.31693.55315.1582
400 K118,327.71783.46665.460.31694.22695.3949
500 K117,144.401269.22473.5860.31694.59405.6414
VCN = 0.17300 K144,771.383493.881303.660.31203.53374.9979
400 K142,387.802711.951011.910.31204.21115.2255
500 K141,057.961930.01720.1440.31204.58215.4627
300 K140,670.005279.791731.780.2560−0.26813.1446
Vf = 0.6400 K140,482.004242.641380.010.2560−0.31823.2665
500 K140,294.003130.951008.980.256−0.37453.3885
Table 2. Distribution and volume fractions of CNT and fiber for hybrid beams.
Table 2. Distribution and volume fractions of CNT and fiber for hybrid beams.
TypesUDFG-XFG-OFG-ΛFG-VSymmetric Laminates in the Case Study
Each Ply
Ply10.60.60.60.60.6 Materials 13 03718 i001
Ply20.140.170.110.110.17
Ply30.140.140.140.110.17
Ply40.140.110.170.140.14
Ply50.140.110.170.140.14
Ply60.140.140.140.170.11
Ply70.140.170.110.170.11
Ply80.60.60.60.60.6
Table 3. EPRs ( v 13 e ) of hybrid laminated beams for various temperature conditions.
Table 3. EPRs ( v 13 e ) of hybrid laminated beams for various temperature conditions.
FG-UDFG-VFG-ΛFG-XFG-O
(25F/20C/−20C/20C)S
300K−0.4566−0.3972−0.3972−0.4213−0.4213
400K−0.5318−0.5021−0.5021−0.5167−0.5167
500K−0.7263−0.6583−0.6583−0.6743−0.6743
(25F/90C/−90C/90C)S
300K0.30100.30160.30160.30170.3017
400K0.30060.30120.30120.30130.3013
500K0.30050.30090.30090.30100.3010
Table 4. Comparison of the fundamental frequencies Ω ˜   = Ω(L2/h) ρ 0 / E 0 for (θC/90F)S hybrid beam with different values of the lamination angle of the CNTRC layer.
Table 4. Comparison of the fundamental frequencies Ω ˜   = Ω(L2/h) ρ 0 / E 0 for (θC/90F)S hybrid beam with different values of the lamination angle of the CNTRC layer.
VCN θC
153045
0.12Fan & Wang [70]5.27645 3.29216 2.79561
Proposed method5.284653.305502.81136
0.17Fan & Wang [70]6.669204.092263.42397
Proposed method6.675564.102903.43674
0.28Fan & Wang [70]7.121074.296873.57689
Proposed method7.126934.306893.58898
Table 5. Natural frequencies Ω ˜ i   = Ωi(a2/h) ρ 0 / E 0 for hybrid laminated beams in thermal environments (L/h = 5, h = 1.75 mm).
Table 5. Natural frequencies Ω ˜ i   = Ωi(a2/h) ρ 0 / E 0 for hybrid laminated beams in thermal environments (L/h = 5, h = 1.75 mm).
T (K)Lay-up Ω ˜ 1 Ω ˜ 2 Ω ˜ 3 Ω ˜ 4
300(25 F/20 C/-20 C/20 C) SUD4.691113.694823.262232.9778
FG-V& Λ 4.923514.368324.388834.5382
FG-X5.021914.347724.120334.0164
FG-O4.807114.510825.035735.7004
(25 F/90 C/-90 C/90 C) SUD3.929012.432022.045131.9209
FG-V& Λ 4.051612.885522.913633.2090
FG-X4.095212.850122.665232.7112
FG-O3.994912.993923.439534.2276
400(25 F/20 C/-20 C/20 C) SUD3.596111.403719.645028.0129
FG-V& Λ 3.902012.097220.758229.5337
FG-X3.924712.000720.415128.9379
FG-O3.714412.178821.300030.5349
(25 F/90 C/-90 C/90 C) SUD3.365210.881819.350828.0456
FG-V& Λ 3.466911.272820.107129.1715
FG-X3.495111.225419.864328.7019
FG-O3.407411.373320.587330.0945
500(25 F/20 C/-20 C/20 C) SUD1.82588.378715.013021.7286
FG-V& Λ 2.49749.276216.323623.4585
FG-X2.25488.942415.708822.5508
FG-O1.99299.145816.582424.0827
(25 F/90 C/-90 C/90 C) SUD2.75469.120216.264423.5940
FG-V& Λ 2.83689.445116.897924.5401
FG-X2.84789.387016.667424.1110
FG-O2.77469.531117.313625.3358
Table 6. Natural frequencies Ω ˜ i   = Ωi(L2/h) ρ 0 / E 0 for hybrid laminated beams resting on elastic foundations (L/h = 20, h=1.75 mm).
Table 6. Natural frequencies Ω ˜ i   = Ωi(L2/h) ρ 0 / E 0 for hybrid laminated beams resting on elastic foundations (L/h = 20, h=1.75 mm).
(k1, k2)Lay-Up Ω ˜ 1 Ω ˜ 2 Ω ˜ 3 Ω ˜ 4
(0, 0)(25 F/20 C/-20 C/20 C) SUD5.517421.269245.245575.0582
FG-V& Λ 5.791222.324447.488878.7764
FG-X6.016723.082348.779880.3508
FG-O5.505421.366145.882676.9132
(25 F/90 C/-90 C/90 C) SUD4.538217.045337.017462.8646
FG-V& Λ 4.479717.533738.121064.8256
FG-X4.566517.837638.677165.5239
FG-O4.358517.112237.377763.9190
(103, 0)(25 F/20 C/-20 C/20 C) SUD10.397023.012646.083775.5615
FG-V& Λ 10.544923.991148.288079.2561
FG-X10.670424.698149.558480.8215
FG-O10.390723.102046.708977.4038
(25 F/90 C/-90 C/90 C) SUD9.831119.176038.036163.4629
FG-V& Λ 9.885519.611439.110965.4059
FG-X9.925119.883639.643566.0983
FG-O9.831219.235438.386664.5070
(103, 102)(25 F/20 C/-20 C/20 C) SUD13.592028.885052.949281.062
FG-V& Λ 13.705429.670554.878586.4789
FG-X13.802230.245456.001687.9216
FG-O13.587128.955853.491284.7756
(25 F/90 C/-90 C/90 C) SUD13.164125.930846.107372.2569
FG-V& Λ 13.204826.254346.997473.9677
FG-X13.234426.458347.442474.5835
FG-O13.164225.974446.394773.1687

Share and Cite

MDPI and ACS Style

Huang, X.-h.; Yang, J.; Azim, I.; Wang, X.-e.; Ren, X. Geometric Non-Linear Analysis of Auxetic Hybrid Laminated Beams Containing CNT Reinforced Composite Materials. Materials 2020, 13, 3718. https://doi.org/10.3390/ma13173718

AMA Style

Huang X-h, Yang J, Azim I, Wang X-e, Ren X. Geometric Non-Linear Analysis of Auxetic Hybrid Laminated Beams Containing CNT Reinforced Composite Materials. Materials. 2020; 13(17):3718. https://doi.org/10.3390/ma13173718

Chicago/Turabian Style

Huang, Xu-hao, Jian Yang, Iftikhar Azim, Xing-er Wang, and Xin Ren. 2020. "Geometric Non-Linear Analysis of Auxetic Hybrid Laminated Beams Containing CNT Reinforced Composite Materials" Materials 13, no. 17: 3718. https://doi.org/10.3390/ma13173718

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