Next Article in Journal
Influence of Different Build Orientations and Heat Treatments on the Creep Properties of Inconel 718 Produced by PBF-LB
Next Article in Special Issue
Computational Analysis of Sandwich Panels with Graded Foam Cores Subjected to Combined Blast and Fragment Impact Loading
Previous Article in Journal
Friction Stir Welding of Dissimilar Al 6061-T6 to AISI 316 Stainless Steel: Microstructure and Mechanical Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Three-Dimensional Vibration Theory for Ultralight Cellular Sandwich Plates Subjected to Linearly Varying In-Plane Distributed Loads

1
The Institute of Xi’an Aerospace Solid Propulsion Technology, Xi’an 710025, China
2
National Key Laboratory of Solid Rocket Propulsion, Xi’an 710025, China
3
School of Mechanical Engineering, Xi’an Jiaotong University, Xi’an 710049, China
4
State Key Laboratory of Mechanics and Control of Mechanical Structures, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
5
Nanjing Center for Multifunctional Lightweight Materials and Structures, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
*
Authors to whom correspondence should be addressed.
Materials 2023, 16(11), 4086; https://doi.org/10.3390/ma16114086
Submission received: 16 April 2023 / Revised: 19 May 2023 / Accepted: 24 May 2023 / Published: 31 May 2023
(This article belongs to the Special Issue Mechanical Behavior and Numerical Simulation of Sandwich Composites)

Abstract

:
Thin structural elements such as large-scale covering plates of aerospace protection structures and vertical stabilizers of aircraft are strongly influenced by gravity (and/or acceleration); thus, exploring how the mechanical behaviors of such structures are affected by gravitational field is necessary. Built upon a zigzag displacement model, this study establishes a three-dimensional vibration theory for ultralight cellular-cored sandwich plates subjected to linearly varying in-plane distributed loads (due to, e.g., hyper gravity or acceleration), with the cross-section rotation angle induced by face sheet shearing accounted for. For selected boundary conditions, the theory enables quantifying the influence of core type (e.g., close-celled metal foams, triangular corrugated metal plates, and metal hexagonal honeycombs) on fundamental frequencies of the sandwich plates. For validation, three-dimensional finite element simulations are carried out, with good agreement achieved between theoretical predictions and simulation results. The validated theory is subsequently employed to evaluate how the geometric parameters of metal sandwich core and the mixture of metal cores and composite face sheets influence the fundamental frequencies. Triangular corrugated sandwich plate possesses the highest fundamental frequency, irrespective of boundary conditions. For each type of sandwich plate considered, the presence of in-plane distributed loads significantly affects its fundamental frequencies and modal shapes.

1. Introduction

With the development of modern flight vehicles, studies on advanced materials, topological structures, and biological processes in hyper gravity environments have evolved into an important research hotspot [1,2,3]. It has been demonstrated that gravity has noticeable effects on the precision of positioning systems, inertial navigators, and guided systems as well as the mechanical performance of thin structural elements applicated in aerospace vehicles [4,5,6]. Typically, such thin structural elements include large-scale covering shells/plates of aerospace protection structures, vertical stabilizers of accelerated aircraft or missiles, and structural components of specific systems strongly influenced by varying gravity [7,8]. These structures, relatively thin, are subjected to hyper gravity as the flight vehicle accelerates. It is thus necessary to explore how the mechanical behaviors (e.g., buckling and vibration) of thin-walled structures are affected by gravitational field.
As an elemental thin structure component, the mechanical performance of a monolithic plate in gravitational field has been extensively studied. At the early stage, the stability and vibration of a rectangular plate subjected to self-weight, or acceleration-induced body force, were characterized by treating the body force as linearly varying in-plane distributed loads (IPDLs) [9,10,11]. Afterwards, the elastic stability and vibration of a standing monolithic plate with different boundary conditions (BCs) were investigated by applying either the classical or the Reissner–Mindlin plate theory [12,13,14,15,16,17]. It has been demonstrated that IPDL significantly affects the mechanical properties of a monolithic plate. For instance, IPDL affects the potential energy of the plate, thus also changing its stability and vibration behaviors. The critical buckling top load decreases with increasing IPDL, while the fundamental frequency can be either greater or less than that of a no-load plate, depending on the type of BCs and the direction of accelerations.
In lieu of monolithic plates, sandwich structures are widely utilized in aircraft, aerospace vehicles, satellites, missiles, and the like. In particular, ultralight sandwich plates with thin face sheets and periodic lattice truss cores, for example, triangular corrugated (TCOR) plates and hexagonal honeycombs (HHON), receive increasing attention due to superior load-bearing capability and additional multifunctional attributes such as energy absorption, active cooling, and noise attenuation [18,19,20]. For instance, it has been demonstrated that TCOR and HHON sandwich plates possess excellent bending strength, blast resistance, and impact energy absorption [21,22]. However, while the influence of hyper gravity or acceleration on either vibration or stability of sandwich plates deserves much attention as in the case of monolithic plates, only the present authors managed to investigate how a corrugated sandwich plate would buckle when subjected to distributed body force under a variety of BCs [23]. At present, the vibration behaviors of ultralight cellular-cored sandwich plates subjected to IPDLs remain elusive.
This study aims to establish a theory to analyze the vibration performance of ultralight cellular-cored sandwich plates subjected to linearly varying IPDLs, with the effects of BCs duly accounted for. In previous theoretical analysis of sandwich plates, they are often regarded as a special case of laminated structure [24]. Additionally, there are four main kinds of theoretical displacement model: three-dimensional elasticity, equivalent single layer, layer wise, and zigzag [25,26]. The three-dimensional elasticity theory regards the laminated structure as a generalized three-dimensional solid model, and does not consider the special layered configuration of the laminates, and the calculation cost is generally relatively high. At present, the equivalent single-layer model has classical theories and various high-order theories, but it ignores the continuity conditions of interlayer displacement and transverse shear stress. Layer wise model and zigzag model can overcome the limitations of the previous two theoretical models by independently modeling the displacement field of each single layer, and achieve a balance between computational accuracy and computational efficiency. However, the unknown displacement functions of the layer wise model increase with the increase of the number of laminate layers. The zigzag model is actually a superposition form of the equivalent single-layer model and interlayer continuous terms. The number of unknown variables is certain, and the calculation amount is less than the layer wise model [27,28]. To sum up, this paper intends to use the zigzag displacement model. Zigzag theory was pioneered by Lekhnitskii for multilayered beams, and Di Sciuva is the first to apply the zigzag model to the vibration analysis of the sandwich structure [29,30]. The bending, buckling, and vibration of simply supported multilayer anisotropic plates are analyzed, and the in-plane displacement is distributed linearly along the thickness of the structure under the condition of continuity between layers. The defect of Di Sciuva’s model is that it cannot be used to analyze the global response of the structure under clamped boundary conditions. Tessler et al. then proposed a refined zigzag model based on Timoshenko’s beam theory [31]. The model is actually an extension of the first-order shear deformation theory, on which the piecewise linear zigzag functions are superimposed. The in-plane displacement and transverse displacement are consistent in the thickness direction, and the shear correction coefficient and shear stress continuity conditions are not required. After that, Iurlaro performed the bending and vibration behavior of the sandwich beam and plate through the first- and high- order zigzag functions [32,33]. It is found that the predicted results of the zigzag displacement model are closer to the experimental results than those of the Timoshenko’s beam model considering the shear correction factor. Furthermore, the researchers developed the refined high-order zigzag theory by means of the cubic and trigonometric functions to analyze the mechanical behaviors of laminated structures such as bending, buckling, stress, and viscoelasticity [34,35,36,37]. It is worth mentioning that the refined zigzag theory can be reduced to the classical straight line displacement field hypothesis in Allen’s works when the in-plane displacements are not considered and the rotations of the cross-section are expressed by -λw/∂x or -λw/∂y [38]. Allen and other related models for solving sandwich structure problems are only used to consider the sandwich structural characteristics of soft core and hard face sheets, and the assumption of the cross-section angle of face sheets and core is the same, and the shear effect of the face sheets is ignored [39].
Unlike most existing theories that ignore the shear effect of face sheets, we demonstrated in previous study [23] the necessity to consider such shear effects, especially for the case of composite face sheets under IPDL. Thus, the displacement field introduced in [23] on the basis of a zigzag hypothesis is adopted to consider the cross-section rotation angle induced by face sheet shearing, so that sufficiently accurate results can be obtained at low computational costs. Subsequently, the Hamilton principle and the Ritz method are utilized to obtain the dynamic governing equations. To compare how different cellular cores for sandwich construction perform dynamically, the geometric parameter optimizations of the metal configuration of TCOR and HHON are performed under the constraint of equal mass. Further, in the presence of IPDL, the influence of composite face sheets on fundamental frequency is explored under three different types of BCs.

2. Formulation

2.1. Problem Description

Consider the vibration of a cellular sandwich plate with a gravitational field g, as shown in Figure 1a. The origin of the Cartesian coordinate system is placed at the geometric center of the plate. Let the x-y-z directions be parallel to the width a L , height L and thickness h of the plate, respectively. Cellular foams, TCOR plates, and HHONs are successively employed as the core to construct the sandwich plate, as shown in Figure 1b. Relevant geometric parameters of the sandwich core are labeled in Figure 1b: core thickness h c ; TCOR member thickness t t and corrugation angle θ ; HHON member thickness t h ; length l h , and angle between horizontal line and inclined cell wall α . In the present study, only regular HHONs are discussed so that α = 30 o .

2.2. Kinematics and Constitutive Equations

Similar to our previous study [23], each type of cellular core displayed in Figure 1 is viewed as an equivalent orthogonal layer. The displacement model considering the shear deformation of face sheet layers is shown in Figure 2. Upon assuming that the points on the cross-section of the core layer and the face sheet layers have the same rotations, the blue line representing these points moves to the new position marked as the red line. The two face sheets are taken as symmetrical with identical thickness (i.e., h f 1 = h f 2 = h f ) and material make. Then, the rotation angles of the cross-sections in the x-z plane are denoted by θ c and θ f for the equivalent core and face sheet layers, respectively, as shown in Figure 2. Due to consideration of shear deformation, the rotation angles are all less than the pure bending angle d w / d x , wherein w represents the consistent transverse displacement of the equivalent core and face sheet layers.
With interface displacement continuity considered, the three constituent layers of the sandwich plate have the following displacement components in the x and y directions:
a.
Face sheet layer 1 ( h c 2 z h 2 ) :
u f 1 ( x , y , z ) = ( z + h c 2 ) θ f ( x , y ) + h c 2 θ c ( x , y ) v f 1 ( x , y , z ) = ( z + h c 2 ) λ f ( x , y ) + h c 2 λ c ( x , y )
b.
Sandwich core layer ( h c 2 z h c 2 ) :
u c ( x , y , z ) = z θ c ( x , y ) v c ( x , y , z ) = z λ c ( x , y )
c.
Face sheet layer 2 ( h c 2 z h 2 ) :
u f 2 ( x , y , z ) = ( z h c 2 ) θ f ( x , y ) h c 2 θ c ( x , y ) v f 2 ( x , y , z ) = ( z h c 2 ) λ f ( x , y ) h c 2 λ c ( x , y )
where λ c and λ f denote the cross-section rotation angles of sandwich core layer and face sheet layers in y-z plane (Figure 2), respectively. It follows that, with small deformation assumed, the linear strain vectors of face sheets and core layer can be obtained as:
ε ( f 1 ,   f 2 ,   c ) = { ε x ε y γ x y γ y z γ x z } ( f 1 ,   f 2 ,   c )   T = { u x v y u y + v x v z + w y u z + w x } ( f 1 ,   f 2 ,   c )   T
For both the face sheets and core, the linear stress–strain relationship can be expressed as:
σ ( f 1 ,   f 2 ,   c ) = { σ x σ y τ x y τ y z τ x z } ( f 1 ,   f 2 ,   c )   T = [ C 11 C 12 0 0 0 C 12 C 22 0 0 0 0 0 C 44 0 0 0 0 0 C 55 0 0 0 0 0 C 66 ] ( f 1 ,   f 2 ,   c ) ε ( f 1 ,   f 2 ,   c )
where C i j ( f 1 ,   f 2 ,   c ) are the elastic coefficients. As the two face sheet layers are completely symmetrical, the expressions of C i j f 1 ,   f 2 can be expressed uniformly as: C 11 f = E 1 f 1 ν 12 f ν 21 f , C 12 f = ν 21 f E 1 f 1 ν 12 f ν 21 f , C 22 f = E 2 f 1 ν 12 f ν 21 f , C 44 f = G 12 f , C 55 f = G 23 f , C 66 f = G 13 f . C i j c refer to the equivalent elastic constants of sandwich core. Detailed expressions of C i j c for TCOR plates and HHONs are given in Appendix A.

2.3. Energy Formulation

To investigate the vibration behavior of a cellular sandwich plate, Hamilton’s principle is adopted: δ t ( U + V T ) d t = 0 . Therefore, the strain energy U , potential energy V , and kinetic energy T of the sandwich are given specifically as follows:
  • Strain energy:
U = 1 2 f 1 ,   f 2 ,   c V ( σ x ε x + σ y ε y + τ x y γ x y + τ y z γ y z + τ x z γ x z ) d V = 1 2 A z ( σ f 1 T ε f 1 + σ c T ε c + σ f 2 T ε f 2 ) d z d A
  • Potential energy generated from IPDLs:
V g = ( ρ m f 1 g h f 1 2 + ρ s g h c 2 + ρ m f 2 g h f 2 2 ) A ( L y ) ( w y ) 2 d A = ( ρ m f g h f + ρ s g h c 2 ) A ( L y ) ( w y ) 2 d A
where ρ m f is the density of the parent material of face sheet layers; ρ s is the equivalent density of sandwich core ( ρ s HB for HHON and ρ s CR for TCOR) given by:
ρ s HB = 2 ρ m c t h 3 l h ρ s CR = ρ m c t t 2 h c sin θ h c 2 h c tan θ = ρ m c t t h c cos θ
where ρ m c is the density of the parent material of cellular sandwich core.
  • Kinetic energy:
T = 1 2 f 1 ,   f 2 ,   c V ρ ( z ) ( u , t 2 + v , t 2 + w , t 2 ) d V = 1 2 A z ( d , t f 1   T ρ m f 1 ( z ) d , t f 1 + d , t c   T ρ s ( z ) d , t c + d , t f 2   T ρ m f 1 ( z ) d , t f 2 ) d z d A
where:
d ( f 1 ,   f 2 ) = { u f 1 ,   f 2 v f 1 ,   f 2 w } T d c = { u c v c w } T
The subscript ‘ , t ’ of displacement vector d ( f 1 ,   f 2 ,   c ) represents the first derivative with respect to time t .

3. Solution Procedures and Validation

3.1. Solution Procedure

Based on the formulation of Section 2, the Ritz method is employed to determine dynamic equations of the cellular sandwich plate and the corresponding eigenvalues [23,40]. Firstly, the Ritz functions are introduced as:
w ¯ ( x ¯ , y ¯ ) = r = 0 p i = 0 r w ¯ ξ ψ ξ w ( x ¯ , y ¯ ) sin ( ω t + φ ) θ c ( x ¯ , y ¯ ) = r = 0 p i = 0 r θ ξ c ψ ξ x c ( x ¯ , y ¯ ) sin ( ω t + φ ) λ c ( x ¯ , y ¯ ) = r = 0 p i = 0 r λ ξ c ψ ξ y c ( x ¯ , y ¯ ) sin ( ω t + φ ) θ f ( x ¯ , y ¯ ) = r = 0 p i = 0 r θ ξ f ψ ξ x f ( x ¯ , y ¯ ) sin ( ω t + φ ) λ f ( x ¯ , y ¯ ) = r = 0 p i = 0 r λ ξ f ψ ξ y f ( x ¯ , y ¯ ) sin ( ω t + φ )
where w ¯ is the non-dimensional transverse displacement, defined as w ¯ = 2 w L ; x ¯ and y ¯ are the non-dimensional coordinates, defined as x ¯ = 2 x L , y ¯ = 2 y L ; p is the degree of the complete polynomial space; w ¯ i ,   θ i c ,   λ i c ,   θ i f ,   λ i f are unknown coefficients to be varied with subscript ξ , defined by ξ = ( r + 1 ) ( r + 2 ) / 2 i ; ω is the natural frequency of sandwich plate, and ψ ξ w ,   ψ ξ x c ,   ψ ξ y c ,   ψ ξ x f ,   ψ ξ y f are the polynomial functions, which contain the basic functions ϕ w ,   ϕ x c ,   ϕ y c ,   ϕ x f ,   ϕ y f set to satisfy the geometric BCs as:
ψ ξ w , x c , y c , x f , y f = ϕ w , x c , y c , x f , y f ( x ¯ , y ¯ ) x ¯ i y ¯ r i = ( x ¯ a ) n 1 ( x ¯ + a ) n 2 ( y ¯ 1 ) n 3 ( y ¯ + 1 ) n 4 x ¯ i y ¯ r i
Here, depending upon the BC type and the direction of cross-section rotation angle, the value of n k   ( k = 1 , 2 , 3 , 4 ) is 0 or 1.
Four different BC types for the four edges of the plate displayed in Figure 3 are involved in subsequent verification study. For example, as listed in Table 1, the boundary condition SCFF means that the 1st edge is simply supported, the 2nd edge is clamped, and the 3rd and 4th edges are free. For different combinations of BC and cross-section rotation angle direction, Table 1 lists the values of n k in order.
Secondly, upon substituting the Ritz functions (Equation (11)) into the displacement field (Equations (1)–(3)) and then applying the strain vectors (Equation (4)) and constitutive equations (Equation (5)), the energy formulas of the cellular sandwich plate can be rewritten as:
U = 1 2 f ,   c A ¯ z ( Λ T Ν ( f ,   c )   T C ( f ,   c ) Ν ( f ,   c ) Λ ) d z d x ¯ d y ¯ V g = 1 2 A ¯ ( 1 y ¯ ) Λ T Ν T q ¯ Ν Λ d x ¯ d y ¯ T = 1 2 A ¯ z ( Λ , t T D ( f ,   c )   T ρ ( z ) D ( f ,   c ) Λ , t ) d z d x ¯ d y ¯
where Ν ( f ,   c ) , Ν , D ( f ,   c ) are the non-dimensional coefficient matrices of strain vectors and displacement vectors. These matrices can be conveniently obtained from Equations (1)–(4) and (11), thus their specific expressions are not presented herein for brevity. Λ is the vector of unknown coefficients, expressed as [ θ ξ c θ ξ f λ ξ c λ ξ f w ¯ ξ ] T sin ( ω t + φ ) . q ¯ is the non-dimensional parameter quantifying IPDLs, defined as ρ g L 3 D eq ; ρ represents the mass per area of sandwich plate in the x-y plane, and D eq is the parameter quantifying the ability of the sandwich plate to resist bending in the y-direction, given by:
ρ = h c / 2 h c / 2 ρ s d z + 2 h c / 2 h / 2 ρ m f d z D eq = h c / 2 h c / 2 z 2 C 22 c d z + 2 h c / 2 h / 2 z 2 C 22 f d z
Finally, by substituting the new energy formulas of Equation (13) into Hamilton’s principle and performing variational calculations on unknown coefficients Λ , the dynamic equations of the cellular sandwich plate expressed in the form of an eigenvalue problem are obtained as:
( K + q ¯ K g ω ¯ 2 M ) Λ = 0
where ω ¯ is the non-dimensional natural frequency, defined as ρ ω 2 L 4 D eq . The structural stiffness matrix K , the geometric stiffness matrix K g induced by IPDLs, and the mass matrix M are all 5th order symmetric matrices, while detailed expressions of K and M are presented in Appendix B and Appendix C, respectively. The geometric stiffness matrix is given by:
K g = A ¯ ( 1 y ¯ ) B T B d x ¯ d y ¯
In the current study, only changes in the first-order natural frequency and modal shape of a cellular sandwich plate are concerned, which can be obtained by solving the eigenvalue problem of Equation (15). With the mass matrix omitted, the critical buckling IPDL can be obtained by calculating the minimum eigenvalue of ( K ,   K g ) .

3.2. Validation

In this section, the three-dimensional (3D) theory of free vibration established above for cellular-cored sandwich plates subjected to linearly varying in-plane distributed loads are validated.
Consider first the limiting case of sandwich plates, that is, the monolithic plate. For a monolithic rectangular plate (a = 0.5; Figure 1), Table 2 compares its critical buckling IPDL q ¯ and natural fundamental frequency ω ¯ calculated from the present method with existing theoretical predictions by Wang et al. [13,14] and Yu et al. [16,17]. Additionally included in Table 2 are 3D finite element (FE) simulation results obtained using the commercially available FE code ABAQUS. As shown in Table 2, overall, the present results agree well with existing theoretical predictions and FE calculations, irrespective of the BCs. Nonetheless, the values obtained for both q ¯ and ω ¯ using the proposed theory are consistently smaller than those of existing theoretical predictions that did not consider shear effect of the monolithic plate, thus more consistent with 3D FE simulation results. Further, while the present predictions are sensitive to the slenderness ratio L/h of the monolithic plate, such sensitivity is absent in the theoretical predictions of Wang et al. [13,14] and Yu et al. [16,17].
To further validate the proposed theory, direct FE simulations of cellular sandwich plates are performed, with 3D deformable four-node shell elements and eight-node solid elements (S4R, C3D8R) selected to model the sandwich core and the face sheets, respectively. Interactions of Tie are exerted between the face sheets and core layer. For simply supported sandwich plates, the freedoms of reference point located at the geometric center of the BC face are constrained, while the points in the BC face are coupled with the reference point. For sandwich plates having either clamped or free edges, the BCs can be straightforwardly implemented in ABAQUS.
For each FE simulation, a linear frequency analysis step is used to extract the natural frequencies and modal shapes. To consider the influence of IPDL, a prior static mechanical analysis before Frequency step is performed. At the static analysis step, a gravity load is applied and the option of large deformation is turned on. For calculating the critical buckling IPDL, only a linear buckle analysis step is required [10,23]. For both face sheets and core, aluminum alloy is selected as the material make, with density ρ Al = 2700   kg / m 3 , Poisson ratio ν = 0.3 , and Young’s modulus E Al = 70   GPa . Irrespective of core type, the sandwich plate has the following geometrical parameters: L = 4   m , h c = 0.09   m , h = 0.1   m , and a = 1 . It follows that the sandwich core has a relative density of ρ ¯ s = ρ s / ρ m = 0.058 .
Figure 4 compares the theoretical and numerical results obtained for both TCOR and HHON sandwich plates in terms of ω ¯ versus q ¯ curves with excellent agreement achieved, irrespective of BCs considered, SSSS, FSFS, and SCFF. With the value of q ¯ fixed, the sandwich plate with SSSS exhibits the largest ω ¯ , followed in order by those with FSFS and SCFF.

4. Results and Discussion

The proposed 3D vibration theory, validated in the previous section, is applied in this section to investigate how the BC and key geometrical/material parameters (e.g., core thickness h c , core density ρ s ¯ , corrugation angle θ , and face sheet material) of a cellular sandwich plate affects vibration behaviors. Firstly, how varying the BC of the sandwich plate affects its fundamental frequency under IPDL is quantified. Secondly, how the fundamental frequency depends upon geometric parameters is evaluated, with equal mass assumed. Thirdly, the frequency versus IPDL curves for four different kinds of plates under varying BCs are compared: monolithic plate, foam sandwich plate, TCOR sandwich plate, and HHON sandwich plate. Further, the material make of the face sheets is varied from metal to fiber-reinforced composite to see how this would affect its vibration performance. Lastly, first-order vibration modal shapes of sandwich plates are discussed.

4.1. Effect of Boundary Condition (BC)

For TCOR and HHON sandwich plates, the results of Figure 4 reveal that different types of BCs have similar influence on plate vibration under IPDL. With the increase of dimensionless critical buckling IPDL q ¯ , the dimensionless natural fundamental frequency ω ¯ exhibits a different downward variation trend as the BC is varied. A sandwich plate with SSSS BC is the most stable and possesses the highest frequency as most degrees of freedom of its four edges are constrained. As for the sandwich plate with two free edges (i.e., either FSFS or SCFF BC), the plate is more stable with opposite edges constrained (FSFS), even though the clamped bottom edge in the constrained adjacent edges (SCFF) restricts more degrees of freedom than the simply supported edge in FSFS. However, as q ¯ is increased, ω ¯ drops faster under FSFS than that under SCFF because the simply supported edge in SCFF is parallel to the direction of IPDL.

4.2. Effect of Geometric Parameters

To compare the mechanical properties (i.e., fundamental frequency and critical buckling IPDL) between TCOR and HHON sandwich plates having equal mass, geometric parameter optimization is performed first. With the cost of computation in mind, only SSSS and FSFS BCs are considered in the optimization. The face sheet thickness is fixed at L / 2 h f = 400 . To visualize the discrepancy among the calculated curves more clearly, non-dimensional natural fundamental frequency 100 ω L ( ρ / E ) Al and non-dimensional IPDL 1000 g L ( ρ / E ) Al are introduced, with the subscript ‘Al’ denoting aluminum alloy. As the mass of face sheets is fixed, using the expressions of equivalent density of sandwich core as shown in Equation (8) leads to the non-dimensional mass of HHON core M ¯ HB and that of TCOR core M ¯ CR as:
M ¯ HB = M HB M f = ρ s HB h c L 2 ρ m h f L 2 = ρ ¯ s HB h c h f = 2 t h 3 l h h c h f
M ¯ CR = M CR M f = ρ s CR h c L 2 ρ m h f L 2 = ρ ¯ s CR h c h f = t t h c cos θ h c h f
where M HB and M CR represent separately the mass of HHON core and TCOR core, and M f represents the mass of a single face sheet. Equation (17) implies that the geometry and mass of HHON core with specified M ¯ HB can be determined by two independent ratios, that is, ρ ¯ s HB and h c / h f . The variations of ρ ¯ s HB correspond to changes in t h / l h . The thickness t h or length l h of the honeycomb cell wall can be determined as one of that which is given. h c / h f determines the thickness of the sandwich core.
Unlike the HHON core, two independent ratios, that is, ρ ¯ s CR and h c / h f , determine the mass of a TCOR core, but not its geometry, because ρ ¯ s CR contains two separate variables: t t / h c and θ . Hence, for the TCOR core, the influence of t t / h c and θ on the frequency versus IPDL curve is firstly discussed under the constraint of equal mass, as shown in Figure 5, with M ¯ CR = 1.04 , ρ ¯ s CR = 0.058 , and h c / h f = 18 . For both the SSSS and FSFS BCs, representative angles of π / 6 , π / 4 and π / 3 are discussed. The results of Figure 5 show that the TCOR sandwich plate exhibits slightly higher frequencies at the combination of bigger inclination angle and thinner corrugation member thickness. However, a larger corrugation angle makes it more difficult to fabricate the sandwich. Hence, in the rest of this study, the corrugation angle is set to π / 3 in consideration of processing.
Figure 6 compares the frequency versus IPDL curves of TCOR and HHON sandwich plates with equal mass ( M ¯ CR = M ¯ HB = 1.04 ) at different combinations of core density ρ ¯ s and thickness h c / h f . Both the fundamental frequency and critical buckling IPDL decrease as core equivalent density is increased (or, equivalently, core thickness is decreased). To a certain extent, as the sandwich plate will degenerate to a monolithic plate as its core density is increased or core thickness is decreased, these results suggest the superiority of sandwich plate over its monolithic counterpart.
The results of Figure 6 reveal that both the critical buckling IPDL and fundamental frequency of TCOR sandwich plates are higher than those of HHON sandwich plates under either SSSS or FSFS BC. The discrepancy between the two sandwich types becomes negligible as core density is increased or core thickness is decreased, since both TCOR and HHON sandwich plates approach a monolithic plate in the limit. Moreover, the superiority of TCOR core to HHON core is more obvious under FSFS BC, compared with SSSS BC, for the TCOR core possesses stronger elastic constants in the vertical direction, resulting in its much higher structural stiffness of the sandwich plate under FSFS BC.
Essentially, changes in the curves of Figure 5 and Figure 6 are attributed to variations of elastic constants when the core type of cellular sandwich plate is varied. As compared in Figure 7, the great in-plane elastic anisotropy exhibited by a TCOR core can be seen intuitively from the discrepancy between C 22 c and C 11 c (or C 12 c ), the former over ten times larger than the latter. Under either FSFS or SSSS BC, this great anisotropy enables a TCOR sandwich plate to become superior relative to an HHON sandwich. Moreover, the in-plane shear elastic constant C 44 c and the transverse shear constant C 55 c of the HHON core are much smaller than those of the TOCR core. This also contributes to the relatively poor performance of an HHON sandwich plate. Overall, the better stability of a TCOR sandwich plate is be attributed to its higher elastic constants, that is, C 22 c , C 44 c , and C 55 c .

4.3. Comparison between Metal and Composite Face Sheets

In nature and engineering applications, HHON sandwiches are more commonly found than TCOR sandwiches. Therefore, in this section, for enhanced in-plane mechanical performance, the metallic face sheets are replaced by composite face sheets to construct HHON sandwich plates. The composite face sheets are made of carbon fiber-reinforced composite T700/3234, with its mechanical properties given by Ref. [41]: E 1 = 110   GPa , E 2 = E 3 = 8.7   GPa , G 12 = G 13 = G 23 = 4   GPa , ν 13 = 0.3 , and ν 21 = ν 23 = 0.32 . Figure 8 illustrates three different fiber stacking types of composite face sheets for HHON sandwich construction.
As shown in Figure 9a for SSSS BC, under the constraint of equal mass, an HHON sandwich plate with composite face sheets has lower fundamental frequency than its counterpart with Al face sheets. Further, varying fiber stacking changes the anisotropy of the HHON sandwich plate, thus also changing its fundamental frequency. Simultaneously, the critical buckling IPDL is dependent upon the degree of in-plane anisotropy. In particular, when the fiber direction is parallel to the IPDL, that is, the Y-direction fiber, the sandwich plate is more stable.
In contrast, under FSFS BC as shown in Figure 9b, both the frequency and critical buckling IPDL are higher for a sandwich plate having Y- or XY-direction fibers laying in composite face sheets, compared to its counterpart with Al face sheets. This is because under the FSFS BC, the advantage of a stronger stiffness (i.e., Y-direction fiber stacking) in the y-direction can be fully exploited.
Figure 9c compares the frequency versus IPDL curves of different sandwich plates with equal mass under SFSF BC, including monolithic plate, all-metallic cellular sandwich plates with different cores, and cellular sandwich plates with composite face sheets. For all the plates considered here, the fundamental frequency remains nearly constant as the IPDL is increased (called Stage I), and then decreases when the IPDL exceeds a specific value (called Stage II). (This transition from Stage I to Stage II may be attributed to the shift of modal shape as the IPDL is increased under SFSF BC, as discussed in Section 4.5.) When the dimensionless IPDL, g / g H O M , is less than 10, the dimensionless frequency of the HHON sandwich plate, ω / ω H O M , is almost the same as that of the TCOR sandwich. As the g / g H O M exceeds 10, the ω / ω H O M of the TCOR sandwich becomes larger than that of the HHON sandwich. However, upon enhancing the x-direction structural stiffness via X-direction and XY-direction fiber stacking in composite face sheets, the fundamental frequency of the HHON sandwich plate remarkably increases, although its critical buckling IPDL decreases.
It may be concluded from the results of Figure 9 that increased frequency and critical buckling IPDL can be realized only under partial BCs for sandwich plates with composite face sheets that have specific fiber stacking types.
Next, to demonstrate the necessity of considering the shear effect of face sheets for both vibration and buckling analysis, Figure 10 compares the results between TCOR sandwich plates with and without considering such shear effect. When shear effect is ignored, the cross-section rotation angles θ f and λ f in Equations (1) and (3) representing the shear effect of face sheets are replaced with w / x and w / y . From Figure 10a, it can be seen that the difference of frequency between cases with and without considering shear effect is small when the face sheets are thin, even less than 2%. However, the difference increases significantly as the face sheets become thicker, which means that shear effect becomes important for relatively thick face sheets. At the same time, the critical buckling IPDLs of Figure 10b reveal that the difference between the two groups of results is about twice as large as that in frequency. This means that the shear effect of face sheets becomes more important when the IPDL is considered.

4.4. Effect of Sandwich Core Type

Figure 9 compares the frequencies and critical buckling IPDLs among the four different types of plate having equal mass: homogeneous (monolithic) plate, aluminum foam-cored sandwich plate, TCOR-cored sandwich plate, and HHON-cored sandwich plate. For simplicity, the face sheets of all the plates are made of Al alloy. The aluminum foam has closed cells and the following material properties: ρ Al - f = 540   kg / m 3 , υ Al - f = 0.3 , and E Al - f = 405   MPa . From Figure 9, it is seen that the frequency of foam-cored sandwich plate is almost thrice than that of the homogeneous plate. As for the TCOR and HHON sandwich plates, the deviation is even larger: the frequencies are more than four times that of the homogeneous plate. Further, the sandwich plates have much larger critical buckling IPDLs than the homogeneous plate, confirming that the former is significantly more stable than the latter. Among the three cellular core types, the TCOR and HHON are more efficient in enhancing the structural stability of a sandwich construction than the foam. In addition, with the highest frequency and largest critical buckling IPDL, the TCOR sandwich plate outperforms the others.
Comparing the structural characteristics of the sandwich plates with different core types, the TCOR and HHON sandwich cores both have the higher thickness than the foam core under the equal mass constraint, since the equivalent core density of either TCOR or HHON is much smaller than that of foam. Therefore, the increase of the core thickness is a benefit to the improvement of the stability of the sandwich plate because it can increase the overall flexural stiffness of the sandwich structure. The flexural stiffness is proportional to the cubic thickness of the structure. However, the premise is that when the loads act, the core plate of TCOR and HHON cannot exhibit local buckling. In other words, the core should not be too thin. Moreover, if the core thickness is kept constant and only the core density is increased, the critical buckling IPDLs of the structure will increase, but the fundamental frequency will first increase and then decrease. As the core density increases the overall stiffness of the structure, the mass of the structure also increases. A similar structural parametric study can be referred to in our previous study [37,42].

4.5. Comparison of Modal Shapes

Figure 11 and Figure 12 present the FE simulated first-order vibration modal shapes of TCOR and HHON sandwich plates subjected to IPDL under different BCs. Overall, the presence of IPDL shifts the maximum modal displacement downward from the plate center. Especially for the cases under SFSF BC, the presence of IPDL leads to an obvious change in modal shape.
Figure 13 displays the theoretically predicted contours of first-order vibration modal shapes considering IPDLs, which are consistent with the FE simulation results of Figure 11 and Figure 12. This again confirms the viability of the three-dimensional vibration theory developed in the present study for cellular-cored sandwich plates. Specifically, as the IPDLs increased, the modal shift of a sandwich plate under SFSF is presented in Figure 14. It is seen that one half-wave along the y-direction evolves into two half-waves.
As previously mentioned, this modal shift corresponds to the transition of the curves displayed in Figure 9c. At the same time, this phenomenon did not occur under SSSS and FSFS BCs when the IPDLs increased. This is because the out-plane displacements of the top and bottom edge are simultaneously constrained, and a wave parallel to the direction of the IPDLs exists in their modal shapes. This will cause their fundamental frequency to drop faster than the SFSF BC under IPDLs. The presence of a wave perpendicular to the direction of the IPDLs can reduce the fundamental frequency affected by the IPDLs. The results in Figure 4 and Figure 9 can prove this phenomenon as the fundamental frequency under SSSS BC drops slower than FSFS BC, and that under SFSF BC is the slowest when the IPDLs are not particularly strong. Moreover, the maximum modal displacement of the modal shapes all move downward from the plate center under the considered BCs. However, the constraints of the out-plane displacement on the top or bottom edge constrained the infinite deformation (an ideal situation) when IPDLs keep increasing. However, the ‘infinite deformation’ can occur under SFSF BC when the IPDLs act on the whole plate area and keep increasing, which looks like a plate rotates around the x axis. Then, an ‘S’ wave paralleled to the direction of IPDLs appeared, and the fundamental frequency drops faster as shown in Figure 9c.

5. Conclusions

With cross-section rotation angles of both the face sheets and core layer considered, a three-dimensional (3D) theory is established to characterize the vibration performance of ultralight cellular-cored sandwich plates subjected to linearly varying in-plane distributed loads (IPDLs) under different boundary conditions (BCs). For validation, 3D finite element simulations are carried out, with good agreement between theory and simulation achieved. The influence of sandwich core type on the fundamental frequency of the sandwich plate is quantified, including close-celled foams, triangular corrugated (TCOR) metal plates, and metal hexagonal honeycombs (HHON). Systematic parametric study is conducted to explore the influence of the key geometric parameters (e.g., core thickness h c , core density ρ s ¯ , and corrugation angle θ ) of the metal cellular sandwich core, the mixture application of composite face sheets, and the BCs on fundamental frequency and critical buckling IPDL. The main conclusions are summarized as follows.
(1)
The fundamental frequency and critical buckling IPDL of a sandwich plate are both much higher than a monolithic plate with equal mass. As the sandwich core, TCOR and HHON are more efficient in enhancing the structural stability than the foam.
(2)
For TCOR metal sandwich plates, the frequency and critical buckling IPDL are not sensitive to the inclination angle of corrugation. However, the frequency and critical buckling IPDL of both TCOR and HHON sandwich plates are quite sensitive to either core density or core thickness.
(3)
The frequency versus IPDL curves and the vibration modal shapes are quite different for sandwich plates under different BCs (i.e., SSSS, SFSF, and FSFS). Especially for the case under SFSF BC, the vibration modal may shift as the IPDL is increased.
(4)
Using fiber-reinforced composite face sheets in lieu of metal face sheets enhances the performance of the HHON sandwich plate in terms of fundamental frequency and critical buckling IPDL. However, such enhancement occurs only under partial BCs with specific types of fiber stacking.

Author Contributions

Formal analysis, F.-H.L. and Y.W.; Writing—original draft, F.-H.L.; Writing—review & editing, B.H., A.-H.Z., K.L. and T.-J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key R&D Program of China (2022YFB4601804, 2022YFB4603103).

Data Availability Statement

The data presented in this study are available on request from the references and the corresponding author. The data are not publicly available due to the privacy of program data.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Equivalent Elastic Constants of Sandwich Cores

Equivalent elastic constants of an empty TCOR core [42]:
C 11 c = E 1 t t cos 3 θ ( 1 ν 12 ν 21 ) h c + E 1 1 ν 12 ν 21 ( t t sin θ h c ) 3 sin θ cos θ C 22 c = 2 E 2 t t sin θ ( 1 ν 12 ν 21 ) h c sin 2 θ C 12 c = ν 12 E 2 t t cos θ ( 1 ν 12 ν 21 ) h c C 44 c = E 1 1 ν 13 ν 31 t t sin 2 θ cos θ ( 1 ν 12 ν 21 ) h c + E 1 1 ν 13 ν 31 ( t t sin θ h c ) 3 ( c o s 2 θ s i n 2 θ ) 2 4 sin θ cos θ C 55 c = G 23 c t t sin 2 θ h c cos θ C 66 c = G 13 c t t sin 2 θ h c cos θ
Equivalent elastic constants of an empty HHON core [43]:
C 11 c = 1 ν 23 ν 32 E 2 E 3 Δ ; C 12 c = ν 21 + ν 31 ν 23 E 2 E 3 Δ ; C 22 c = 1 ν 13 ν 31 E 1 E 3 Δ C 44 c = E ( t h l b ) 3 × ( 1 + sin α ) N cos α C 55 c = G t h / l h ( 1 + sin α ) cos α × [ l b cos 2 α l h + l b 2 l h + 3 t h tan α 4 l h t h sin 2 α ( 2 sin α 1 ) 2 l h cos α ] C 66 c = G t h / l h ( 1 + sin α ) cos α × [ l b cos 2 α l h + 3 t h tan α 4 l h t h cos α ( 2 sin α 1 ) 2 l h ]
E 1 = E ( t h l b ) 3 × cos α ( 1 + sin α ) sin 2 α [ l b 2 l b 2 + ( 2.4 + 1.5 ν + cot 2 α ) t h 2 ] E 2 = E ( t h l b ) 3 × 1 + sin α cos 3 α [ l b 2 l b 2 + ( 2.4 + 1.5 ν + tan 2 α + 2 h b l b cos 2 α ) t h 2 ] E 3 = E ( 1 l b 2 l h 2 ) ν 12 = cos 2 α ( 1 + sin α ) sin α × l b 2 + ( 1.4 + 1.5 ν ) t h 2 l b 2 + ( 2.4 + 1.5 ν + cot 2 α ) t h 2 ν 21 = ( 1 + sin α ) sin α cos 2 α × l b 2 + ( 1.4 + 1.5 ν ) t h 2 l b 2 + ( 2.4 + 1.5 ν + tan 2 α + 2 h b l b cos 2 α ) t h 2 ν 31 = ν 32 = ν ; ν 13 = E 1 ν 31 E 3 ; ν 23 = E 2 ν 32 E 3 Δ = 1 ν 12 ν 12 ν 23 ν 32 ν 13 ν 31 2 ν 21 ν 32 ν 13 E 1 E 2 E 3 l b = l h t h 2 cos α N = l b 2 3 l b 2 + t h 2 [ ( 2.4 + 1.5 ν ) ( 3 + sin α ) + ( 1 + sin α ) ( ( 1 + sin α ) tan 2 α + sin α ) ]
where E and ν are the Young’s modulus and Poisson ratio of the parent material.

Appendix B. Elements of Structural Matrix K

K = [ K α α cc K α α cf K α β cc K α β cf K α w c K α α ff K α β fc K α β ff K α w f K β β cc K β β cf K β w c Sym K β β ff K β w f K w w ]
where:
K α α i j cc = 1 D eq { ( h c / 2 h c / 2 z 2 C 11 c d z + h c / 2 h / 2 ( h c 2 ) 2 C 11 f d z ) A ¯ ψ ξ x c x ¯ ψ ξ x c x ¯ d x ¯ d y ¯ + ( h c / 2 h c / 2 z 2 C 44 c d z + h c / 2 h / 2 ( h c 2 ) 2 C 44 f d z ) A ¯ ψ ξ x c y ¯ ψ ξ x c y ¯ d x ¯ d y ¯ + L 2 h c / 2 h c / 2 C 66 c d z A ¯ ψ ξ x c ψ ξ x c d x ¯ d y ¯ }
K α α i j cf = 1 D eq { h c / 2 h / 2 [ ( h c 2 ) 2 + z h c ] C 11 f d z A ¯ ψ ξ x c x ¯ ψ ξ x f x ¯ d x ¯ d y ¯ + h c / 2 h / 2 [ ( h c 2 ) 2 + z h c ] C 44 f d z A ¯ ψ ξ x c y ¯ ψ ξ x f y ¯ d x ¯ d y ¯ }
K α β i j cc = 1 D eq { h c / 2 h c / 2 z 2 C 12 c d z A ¯ ψ ξ x c x ¯ ψ ξ y c y ¯ d x ¯ d y ¯ + h c / 2 h c / 2 z 2 C 12 c d z A ¯ ψ ξ x c x ¯ ψ ξ x f y ¯ d x ¯ d y ¯ + h c / 2 h / 2 ( h c 2 ) 2 C 12 f d z A ¯ ψ ξ x c x ¯ ψ ξ y c y ¯ d x ¯ d y ¯ + h c / 2 h / 2 ( h c 2 ) 2 C 44 f d z A ¯ ψ ξ x c y ¯ ψ ξ y c x ¯ d x ¯ d y ¯ }
K α β i j cf = 1 D eq { h c / 2 h / 2 [ ( h c 2 ) 2 + z h c ] C 12 f d z A ¯ ψ ξ x c x ¯ ψ ξ y f y ¯ d x ¯ d y ¯ + h c / 2 h / 2 [ ( h c 2 ) 2 + z h c ] C 44 f d z A ¯ ψ ξ x c y ¯ ψ ξ y f x ¯ d x ¯ d y ¯ }
K α w i j c = L 2 D eq h c / 2 h c / 2 C 66 c d z A ¯ ψ ξ x c ψ ξ w x ¯ d x ¯ d y ¯
K α α i j ff = 1 D eq { h c / 2 h / 2 ( z + h c 2 ) 2 C 11 f d z A ¯ ψ ξ x f x ¯ ψ ξ x f x ¯ d x ¯ d y ¯ + h c / 2 h / 2 ( z + h c 2 ) 2 C 44 f d z A ¯ ψ ξ x f y ¯ ψ ξ x f y ¯ d x ¯ d y ¯ + L 2 h c / 2 h / 2 C 66 f d z A ¯ ψ ξ x f ψ ξ x f d x ¯ d y ¯ }
K α β i j fc = 1 D eq { [ h c / 2 h / 2 ( h c 2 ) 2 C 12 f d z + h c / 2 h / 2 z h c 2 C 12 f d z ] A ¯ ψ ξ x f x ¯ ψ ξ y c y ¯ d x ¯ d y ¯ + [ h c / 2 h / 2 ( h c 2 ) 2 C 44 f d z + h c / 2 h / 2 z h c 2 C 44 f d z ] A ¯ ψ ξ x f y ¯ ψ ξ y c x ¯ d x ¯ d y ¯ }
K α β i j ff = 1 D eq { h c / 2 h / 2 ( z + h c 2 ) 2 C 12 f d z A ¯ ψ ξ x f x ¯ ψ ξ y f y ¯ d x ¯ d y ¯ + h c / 2 h / 2 ( z + h c 2 ) 2 C 44 f d z A ¯ ψ ξ x f y ¯ ψ ξ y f x ¯ d x ¯ d y ¯ }
K α w i j f = L 2 D eq h c / 2 h / 2 C 66 f d z A ¯ ψ ξ x f ψ ξ w x ¯ d x ¯ d y ¯
K β β i j cc = 1 D eq { [ h c / 2 h c / 2 z 2 C 22 c d z + h c / 2 h / 2 ( h c 2 ) 2 C 22 f d z ] A ¯ ψ ξ y c y ¯ ψ ξ y c y ¯ d x ¯ d y ¯ + [ h c / 2 h c / 2 z 2 C 44 c d z + h c / 2 h / 2 ( h c 2 ) 2 C 44 f d z ] A ¯ ϕ i y c x ¯ ϕ j y c x ¯ d x ¯ d y ¯ + L 2 h c / 2 h c / 2 C 55 c d z A ¯ ψ ξ y c ψ ξ y c d x ¯ d y ¯ }
K β β i j cf = 1 D eq { [ h c / 2 h / 2 ( h c 2 ) 2 C 22 f d z + h c / 2 h / 2 z h c 2 C 22 f d z ] A ¯ ψ ξ y c y ¯ ψ ξ y f y ¯ d x ¯ d y ¯ + [ h c / 2 h / 2 ( h c 2 ) 2 C 44 f d z + h c / 2 h / 2 z h c 2 C 44 f d z ] A ¯ ψ ξ y c x ¯ ψ ξ y f x ¯ d x ¯ d y ¯ }
K β w i j c = L 2 D eq h c / 2 h c / 2 C 55 c d z A ¯ ψ ξ y c ψ ξ w y ¯ d x ¯ d y ¯
K β β i j ff = 1 D eq { h c / 2 h / 2 ( z + h c 2 ) 2 C 22 f d z A ¯ ψ ξ y f y ¯ ψ ξ y f y ¯ d x ¯ d y ¯ + h c / 2 h / 2 ( z + h c 2 ) 2 C 44 f d z A ¯ ψ ξ y f x ¯ ψ ξ y f x ¯ d x ¯ d y ¯ + L 2 h c / 2 h / 2 C 66 f d z A ¯ ψ ξ y f ψ ξ y f d x ¯ d y ¯ }
K β w i j f = L 2 D eq h c / 2 h / 2 C 55 f d z A ¯ ψ ξ y f ψ ξ w y ¯ d x ¯ d y ¯
K w w i j = 1 D eq { [ h c / 2 h c / 2 C 55 c d z + h c / 2 h / 2 C 55 f d z ] A ¯ ψ ξ w y ¯ ψ ξ w y ¯ d x ¯ d y ¯ + [ h c / 2 h c / 2 C 66 c d z + h c / 2 h / 2 C 66 f d z ] A ¯ ψ ξ w x ¯ ψ x w x ¯ d x ¯ d y ¯ }

Appendix C. Elements of Mass Matrix M

M = [ m 11 m 12 m 13 m 14 m 15 m 22 m 23 m 24 m 25 m 33 m 34 m 35 Sym m 44 m 45 m 55 ]
where:
m 11 = 1 ρ L 2 [ h c / 2 h c / 2 ρ s z 2 d z + h / 2 h c / 2 ρ f ( h c / 2 ) 2 d z + h c / 2 h / 2 ρ f ( h c / 2 ) 2 d z ] A ¯ ψ ξ x c ψ ξ x c d x ¯ d y ¯ m 12 = 1 ρ L 2 [ h / 2 h c / 2 ρ f h c 2 ( z + h c 2 ) d z + h c / 2 h / 2 ρ f h c 2 ( z h c 2 ) d z ] A ¯ ψ ξ x c ψ ξ x f d x ¯ d y ¯ m 13 = m 14 = m 15 = 0 m 22 = 1 ρ L 2 [ h / 2 h c / 2 ρ f ( z + h c 2 ) 2 d z + h c / 2 h / 2 ρ f ( z h c 2 ) 2 d z ] A ¯ ψ ξ x f ψ ξ x f d x ¯ d y ¯ m 23 = m 24 = m 25 = 0 m 33 = 1 ρ L 2 [ h c / 2 h c / 2 ρ s z 2 d z + h / 2 h c / 2 ρ f ( h c / 2 ) 2 d z + h c / 2 h / 2 ρ f ( h c / 2 ) 2 d z ] A ¯ ψ ξ y c ψ ξ y c d x ¯ d y ¯ m 34 = 1 ρ L 2 [ h / 2 h c / 2 ρ f h c 2 ( z + h c 2 ) d z + h c / 2 h / 2 ρ f h c 2 ( z h c 2 ) d z ] A ¯ ψ ξ y c ψ ξ y f d x ¯ d y ¯ m 35 = 0 m 44 = 1 ρ L 2 [ h / 2 h c / 2 ρ f ( z + h c 2 ) 2 d z + h c / 2 h / 2 ρ f ( z h c 2 ) 2 d z ] A ¯ ψ ξ y f ψ ξ y f d x ¯ d y ¯ m 45 = 0 m 55 = A ¯ ψ ξ w ψ ξ w d x ¯ d y ¯

References

  1. Siddhardha, K. Autonomous reduced-gravity enabling quadrotor test-bed: Design, modelling and flight test analysis. Aerosp. Sci. Technol. 2019, 86, 64–77. [Google Scholar] [CrossRef]
  2. Zhao, Z.; Zhao, H.; Chang, Z.; Feng, X. Analysis of bending and buckling of pre-twisted beams: A bioinspired study. Acta Mech. Sin. 2014, 30, 507–515. [Google Scholar] [CrossRef]
  3. Buffone, C.; Grishaev, V.; Glushchuk, A. Experimental investigation of liquid retention in a cyclone evaporator under variable gravity conditions. Appl. Therm. Eng. 2016, 99, 235–243. [Google Scholar] [CrossRef]
  4. Yao, Z.; Zhao, R.; Zega, V.; Corigliano, A. A metaplate for complete 3D vibration isolation. Eur. J. Mech. A Solids 2020, 84, 104016. [Google Scholar] [CrossRef]
  5. Bahaadini, R.; Saidi, A.R. Stability analysis of thin-walled spinning reinforced pipes conveying fluid in thermal environment. Eur. J. Mech. A Solids 2018, 72, 298–309. [Google Scholar] [CrossRef]
  6. Simons, D.A.; Leissa, A.W. Vibration of rectangular cantilever plates subjected to in-plane acceleration loads. J. Sound Vib. 1971, 17, 407–422. [Google Scholar] [CrossRef]
  7. Saxena, I.F.; Guzman, N.; Hui, K.; Mal, A.K. Disbond detection in a composite honeycomb structure of an aircraft vertical stabilizer by fiber Bragg grating detecting guided ultrasound waves. J. Mech. Eng. Sci. 2017, 231, 3001–3010. [Google Scholar] [CrossRef]
  8. Xu, Q.; Liu, J. An improved dynamic model for silicone material beam with large deformation. Acta Mech. Sin. 2018, 34, 744–753. [Google Scholar] [CrossRef]
  9. Wang, T.M.; Sussaman, J.M. Elastic stability of a simply supported plate under linearly variable compressive stresses. AIAA J. 1967, 5, 1362–1364. [Google Scholar] [CrossRef]
  10. Han, B.; Li, F.H.; Zhang, Q.C.; Chen, C.Q.; Lu, T.J. Stability and initial post-buckling analysis of a standing sandwich beam under terminal force and self-weight. Arch. Appl. Mech. 2016, 86, 1063–1082. [Google Scholar] [CrossRef]
  11. Fauconneau, G.; Marangoni, R.D. Natural frequencies and elastic stability of a simply-supported rectangular plate under linearly varying compressive loads. Int. J. Solids Struct. 1971, 7, 473–493. [Google Scholar] [CrossRef]
  12. Brown, C.J. Elastic buckling of plates subjected to distributed tangential loads. Comput. Struct. 1991, 41, 151–1565. [Google Scholar] [CrossRef]
  13. Wang, C.M.; Xiang, Y.; Wang, C.Y. Buckling of standing vertical plates under body forces. Int. J. Struct. Stab. Dyn. 2002, 2, 151–161. [Google Scholar] [CrossRef]
  14. Wang, C.Y. Buckling of a heavy standing plate with top load. Thin Walled Struct. 2010, 48, 127–133. [Google Scholar] [CrossRef]
  15. Bodaghi, M.; Saidi, A.R. Buckling behavior of standing laminated Mindlin plates subjected to body force and vertical loading. Compos. Struct. 2011, 93, 538–547. [Google Scholar] [CrossRef]
  16. Yu, L.H.; Wang, C.Y. Fundamental frequency of a standing heavy plate with vertical simply-supported edges. J. Sound Vib. 2009, 321, 1–7. [Google Scholar] [CrossRef]
  17. Yu, L.H.; Wang, C.Y. Vibration of a standing plate with simply supported vertical sides and weakened by a horizontal hinge. Thin Walled Struct. 2011, 49, 899–901. [Google Scholar] [CrossRef]
  18. Wang, W.B.; Yang, X.H.; Han, B.; Zhang, Q.C.; Wang, X.F.; Lu, T.J. Analytical design of effective thermal conductivity for fluid-saturated prismatic cellular metal honeycombs. Theor. Appl. Mech. Lett. 2016, 6, 69–75. [Google Scholar] [CrossRef]
  19. Han, B.; Yue, Z.S.; Wu, H.; Zhang, Q.; Lu, T.J. Superior compressive performance of hierarchical origami-corrugation metallic sandwich structures based on selective laser melting. Compos. Struct. 2022, 300, 116181. [Google Scholar] [CrossRef]
  20. Tang, Y.F.; Li, F.H.; Xin, F.; Lu, T.J. Heterogeneously perforated honeycomb-corrugation hybrid sandwich panel as sound absorber. Mater. Des. 2017, 134, 502–512. [Google Scholar] [CrossRef]
  21. Wu, H.X.; Liu, Y.; Zhang, X.C. In-plane crushing behavior and energy absorption design of composite honeycombs. Acta Mech. Sin. 2018, 34, 1108–1123. [Google Scholar] [CrossRef]
  22. Wang, X.; Li, X.; Yue, Z.S.; Zhang, Q.C.; Du, S.F.; Yang, Z.K.; Han, B.; Lu, T.J. Optimal design of metallic corrugated sandwich panels with polyurea-metal laminate face sheets for simultaneous vibration attenuation and structural stiffness. Compos. Struct. 2021, 256, 112994. [Google Scholar] [CrossRef]
  23. Li, F.H.; Han, B.; Zhang, Q.C.; Jin, F.; Lu, T.J. Buckling of a standing corrugated sandwich plate subjected to body force and terminal load. Thin Walled Struct. 2018, 127, 688–699. [Google Scholar] [CrossRef]
  24. Volmir, A.S. Stability of sandwich plates and shells. In Stability of Deformable S.; Nauka Publishing House: Moscow, Russia, 1967; Chapter 19. (In Russian) [Google Scholar]
  25. Liew, K.M.; Pan, Z.Z.; Zhang, L.W. An overview of layerwise theories for composite laminates and structures: Development, numerical implementation and application. Compos. Struct. 2019, 216, 240–259. [Google Scholar] [CrossRef]
  26. Sayyad, A.S.; Ghugal, Y.M. Bending, buckling and free vibration of laminated composite and sandwich beams: Acritical review of literature. Compos. Struct. 2017, 171, 486–504. [Google Scholar] [CrossRef]
  27. Iurlaro, L.; Gherlone, M.; Di Sciuva, M.; Tessler, A. Assessment of the refined zigzag theory for bending, vibration, and buckling of sandwich plates: A comparative study of different theories. Compos. Struct. 2013, 106, 777–792. [Google Scholar] [CrossRef]
  28. Akif, K.; Mehmet, D.; Timon, R. A novel mixed finite element formulation based on the refined zigzag theory for the stress analysis of laminated composite plates. Compos. Struct. 2021, 267, 113886. [Google Scholar]
  29. Lekhnitskii, S.G. Strength calculation of composite beams. Vestn. Inzhen Tekhnikov 1935, 9, 137–148. [Google Scholar]
  30. Di Sciuva, M. Bending, vibration and buckling of simply supported thick multilayered orthortopic plates. J. Sound Vib. 1986, 105, 425–442. [Google Scholar] [CrossRef]
  31. Tessler, A.; Di Sciuva, M.; Gherlone, M. A refined zigzag beam theory for composite and sandwich beams. J. Compos. Mater. 2009, 43, 1051–1081. [Google Scholar] [CrossRef]
  32. Iurlaro, L.; Ascione, A.; Gherlone, M. Free vibration analysis of sandwich beams using the refined zigzag theory: An experimental assessment. Meccanica 2015, 50, 2525–2535. [Google Scholar] [CrossRef]
  33. Iurlaro, L.; Gherlone, M.; Di Sciuva, M. Bending and free vibration analysis of functionally graded sandwich plates using the refined zigzag theory. J. Sandw. Struct. Mater. 2014, 16, 669–699. [Google Scholar] [CrossRef]
  34. Nguyen, S.N.; Lee, J.; Cho, M. Atriangular finite element using Laplace transform for viscoelastic laminated composite plates based on efficient higher-order zigzag theory. Compos. Struct. 2016, 155, 223–244. [Google Scholar] [CrossRef]
  35. Dorduncu, M. Stress analysis of sandwich plates with functionally graded cores using peridynamic differential operator and refined zigzag theory. Thin Walled Struct. 2020, 146, 106468. [Google Scholar] [CrossRef]
  36. Di Sciuva, M.; Sorrenti, M. Bending, free vibration and buckling of functionally graded carbon nanotube-reinforced sandwich plates, using the extended refined zigzag theory. Compos. Struct. 2019, 227, 111324. [Google Scholar] [CrossRef]
  37. Han, B.; Hui, W.W.; Zhang, Q.C.; Zhao, Z.Y.; Jin, F.; Zhang, Q.; Lu, T.J.; Lu, B.H. A refined quasi-3D zigzag beam theory for free vibration and stability analysis of multilayered composite beams subjected to thermomechanical loading. Compos. Struct. 2018, 204, 620–633. [Google Scholar] [CrossRef]
  38. Allen, H.G. Analysis and Design of Structural Sandwich Panels; Pergamon Press: Oxford, UK, 1969. [Google Scholar]
  39. Lou, J.; Ma, L.; Wu, L.Z. Free vibration analysis of simply supported sandwich beams with lattice truss core. Mater. Sci. Eng. B 2012, 177, 1712–1716. [Google Scholar] [CrossRef]
  40. Wang, C.M.; Aung, T.M. Plastic buckling analysis of thick plates using p-Ritz method. Int. J. Solids Struct. 2007, 44, 6239–6255. [Google Scholar] [CrossRef]
  41. Li, M.; Wu, L.Z.; Ma, L. Structural response of all-composite pyramidal truss core sandwich columns in end compression. Compos. Struct. 2011, 93, 1964–1972. [Google Scholar] [CrossRef]
  42. Han, B.; Qin, K.K.; Zhang, Q.C.; Zhang, Q.; Lu, T.J.; Lu, B.H. Free vibration and buckling of foam-filled composite corrugated sandwich plates under thermal loading. Compos. Struct. 2017, 172, 173–189. [Google Scholar] [CrossRef]
  43. Malek, S.; Gibson, L. Effective elastic properties of periodic HHONs. Mech. Mater. 2015, 91, 226–240. [Google Scholar] [CrossRef]
Figure 1. (a) Cellular sandwich plate subject to linearly varying in-plane distributed load (IPDL) and (b) cellular foam, triangular corrugated (TCOR) metal plate, and hexagonal metal honeycomb (HHON) as sandwich core.
Figure 1. (a) Cellular sandwich plate subject to linearly varying in-plane distributed load (IPDL) and (b) cellular foam, triangular corrugated (TCOR) metal plate, and hexagonal metal honeycomb (HHON) as sandwich core.
Materials 16 04086 g001
Figure 2. Displacement model considering cross-sectional rotation due to shear of face sheet layers (x-z plane) [23].
Figure 2. Displacement model considering cross-sectional rotation due to shear of face sheet layers (x-z plane) [23].
Materials 16 04086 g002
Figure 3. Order of plate edges for boundary conditions specified in Table 1.
Figure 3. Order of plate edges for boundary conditions specified in Table 1.
Materials 16 04086 g003
Figure 4. Comparison between theoretical predictions (TH) and FE simulations (FE): natural fundamental frequency ω ¯ plotted as a function of critical buckling IPDL parameter q ¯ for (a) TCOR and (b) HHON sandwich plates under selected BCs.
Figure 4. Comparison between theoretical predictions (TH) and FE simulations (FE): natural fundamental frequency ω ¯ plotted as a function of critical buckling IPDL parameter q ¯ for (a) TCOR and (b) HHON sandwich plates under selected BCs.
Materials 16 04086 g004
Figure 5. Effect of inclination angle on frequency versus IPDL curves for TCOR sandwich plates having equal mass ( M ¯ CR = 1.04 , ρ ¯ s CR = 0.058 , h c / h f = 18 ) with (a) SSSS BC and (b) FSFS BC.
Figure 5. Effect of inclination angle on frequency versus IPDL curves for TCOR sandwich plates having equal mass ( M ¯ CR = 1.04 , ρ ¯ s CR = 0.058 , h c / h f = 18 ) with (a) SSSS BC and (b) FSFS BC.
Materials 16 04086 g005
Figure 6. Effect of core thickness h c / h f or equivalent core density ρ ¯ s on frequency versus IPDL curves of TCOR and HHON sandwich plates under equal mass ( M ¯ CR = M ¯ HB = 1.04 , θ = π / 3 ) with (a) SSSS BC and (b) FSFS BC.
Figure 6. Effect of core thickness h c / h f or equivalent core density ρ ¯ s on frequency versus IPDL curves of TCOR and HHON sandwich plates under equal mass ( M ¯ CR = M ¯ HB = 1.04 , θ = π / 3 ) with (a) SSSS BC and (b) FSFS BC.
Materials 16 04086 g006
Figure 7. (a) In-plane and (b) transverse elastic constants plotted as functions of equivalent density: comparison between TCOR and HHON cores.
Figure 7. (a) In-plane and (b) transverse elastic constants plotted as functions of equivalent density: comparison between TCOR and HHON cores.
Materials 16 04086 g007
Figure 8. Fiber stacking types of composite face sheets for sandwich plate construction: X-direction and Y-direction refer to a single layer of fiber laying in x- and y-direction, respectively; XY-direction denotes two stacking layers of fibers laying in x- and y-direction, respectively, with x-direction fibers placed inside and closer to sandwich core.
Figure 8. Fiber stacking types of composite face sheets for sandwich plate construction: X-direction and Y-direction refer to a single layer of fiber laying in x- and y-direction, respectively; XY-direction denotes two stacking layers of fibers laying in x- and y-direction, respectively, with x-direction fibers placed inside and closer to sandwich core.
Materials 16 04086 g008
Figure 9. Comparison of dimensionless frequencies of different sandwich plates having equal mass (with M ¯ CR = 1.04 , ρ ¯ s = 0.127 , h c / h f = 8.2 ). ω HOM and g HOM denote the fundamental frequency and critical buckling IPDL of a homogeneous plate, respectively. The face sheets are made of either Al alloy (solid lines) or fiber-reinforced composite laminates (T700/3234, dotted lines).
Figure 9. Comparison of dimensionless frequencies of different sandwich plates having equal mass (with M ¯ CR = 1.04 , ρ ¯ s = 0.127 , h c / h f = 8.2 ). ω HOM and g HOM denote the fundamental frequency and critical buckling IPDL of a homogeneous plate, respectively. The face sheets are made of either Al alloy (solid lines) or fiber-reinforced composite laminates (T700/3234, dotted lines).
Materials 16 04086 g009
Figure 10. The calculated (a) natural frequency and (b) critical buckling IPDL for TCOR sandwich plates under SCFF BC, considering both the case with the shear effect of face sheets and that without the shear effect of face sheets. The face sheets are made by the XY-direction fiber composite laminates.
Figure 10. The calculated (a) natural frequency and (b) critical buckling IPDL for TCOR sandwich plates under SCFF BC, considering both the case with the shear effect of face sheets and that without the shear effect of face sheets. The face sheets are made by the XY-direction fiber composite laminates.
Materials 16 04086 g010
Figure 11. FE calculated first-order vibration modal shapes of TCOR sandwich plates under three types of BC: (a1c1) IPDL not considered; (a2c2) IPDL considered.
Figure 11. FE calculated first-order vibration modal shapes of TCOR sandwich plates under three types of BC: (a1c1) IPDL not considered; (a2c2) IPDL considered.
Materials 16 04086 g011
Figure 12. FE calculated first-order vibration modal shapes of HHON sandwich plates under three types of BC: (a1c1) IPDL not considered; (a2c2) IPDL considered.
Figure 12. FE calculated first-order vibration modal shapes of HHON sandwich plates under three types of BC: (a1c1) IPDL not considered; (a2c2) IPDL considered.
Materials 16 04086 g012
Figure 13. Theoretically predicted contours of first-order vibration modal shapes for sandwich plates subject to IPDL under three types of BC.
Figure 13. Theoretically predicted contours of first-order vibration modal shapes for sandwich plates subject to IPDL under three types of BC.
Materials 16 04086 g013
Figure 14. First-order vibration modal shape of a sandwich plate with SFSF BC subjected to: (1) small IPDL and (2) large IPDL.
Figure 14. First-order vibration modal shape of a sandwich plate with SFSF BC subjected to: (1) small IPDL and (2) large IPDL.
Materials 16 04086 g014
Table 1. Values of n k   ( k = 1 , 2 , 3 , 4 ) for selected combinations of cross-section rotation angle direction and BC.
Table 1. Values of n k   ( k = 1 , 2 , 3 , 4 ) for selected combinations of cross-section rotation angle direction and BC.
BC
ϕ x c , ϕ x f
ϕ y c , ϕ y f
ϕ w
SSSS (1234)0,0,1,11,1,0,01,1,1,1
FSFS (1234)0,0,1,10,0,0,00,0,1,1
SFSF (1234)0,0,0,01,1,0,01,1,0,0
SCFF a (1234 b)0,0,0,10,1,0,10,1,0,1
a S—simply supported (rotational degree of freedom around the edge-axis unconstrained); C— clamped (six degrees of freedom constrained); F—free (no degree of freedom constrained). b Order of plate edges are specified in Figure 3.
Table 2. Dimensionless natural fundamental frequency ( ω ¯ ) and critical buckling IPDL ( q ¯ ) of a monolithic plate (a = 0.5): comparison between the present predictions and existing theoretical and FE results.
Table 2. Dimensionless natural fundamental frequency ( ω ¯ ) and critical buckling IPDL ( q ¯ ) of a monolithic plate (a = 0.5): comparison between the present predictions and existing theoretical and FE results.
CaseMethodSSSFSFSFSCSFSSSS
ω ¯ ( q ¯ = 0 ) FE* (L/h = 40)40.20838.16840.69247.859
FE (L/h = 100)40.52638.35341.01148.414
Present study (L/h = 40)40.99838.77641.48649.098
Present study (L/h = 100)41.16538.91341.66749.308
Yu et al. [16,17]41.20438.95041.70649.351
q ¯ ( ω ¯ = 0 ) FE (L/h = 40)195.64687.064280.818195.858
FE (L/h = 100)203.09293.052287.700203.372
Present study (L/h = 40)210.30499.528299.924210.583
Present study (L/h = 100)213.167101.303305.095213.468
Wang et al. [13,14]213.72--306.09214.02
* FE results calculated with FE code ABAQUS (element type: C3D8R); existing theoretical predictions [13,14,16,17] did not account for shear effect of the monolithic plate and are independent of plate slenderness ratio L/h [10].
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, F.-H.; Han, B.; Zhang, A.-H.; Liu, K.; Wang, Y.; Lu, T.-J. A Three-Dimensional Vibration Theory for Ultralight Cellular Sandwich Plates Subjected to Linearly Varying In-Plane Distributed Loads. Materials 2023, 16, 4086. https://doi.org/10.3390/ma16114086

AMA Style

Li F-H, Han B, Zhang A-H, Liu K, Wang Y, Lu T-J. A Three-Dimensional Vibration Theory for Ultralight Cellular Sandwich Plates Subjected to Linearly Varying In-Plane Distributed Loads. Materials. 2023; 16(11):4086. https://doi.org/10.3390/ma16114086

Chicago/Turabian Style

Li, Fei-Hao, Bin Han, Ai-Hua Zhang, Kai Liu, Ying Wang, and Tian-Jian Lu. 2023. "A Three-Dimensional Vibration Theory for Ultralight Cellular Sandwich Plates Subjected to Linearly Varying In-Plane Distributed Loads" Materials 16, no. 11: 4086. https://doi.org/10.3390/ma16114086

APA Style

Li, F. -H., Han, B., Zhang, A. -H., Liu, K., Wang, Y., & Lu, T. -J. (2023). A Three-Dimensional Vibration Theory for Ultralight Cellular Sandwich Plates Subjected to Linearly Varying In-Plane Distributed Loads. Materials, 16(11), 4086. https://doi.org/10.3390/ma16114086

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