Next Article in Journal
Stress Suppression Design for Radiofrequency Microelectromechanical System Switch Based on a Flexible Substrate
Previous Article in Journal
Experimental and Numerical Investigations on the Seismic Performance of High-Strength Exterior Beam-Column Joints with Steel Fibers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Forced-Vibration Characteristics of Bowtie-Shaped Honeycomb Composite Sandwich Panel with Viscoelastic Damping Layer

1
School of Civil Engineering, Chongqing University, Chongqing 400045, China
2
Key Laboratory of New Technology for Construction of Cities in Mountain Area, Chongqing University, Chongqing 400045, China
*
Author to whom correspondence should be addressed.
Materials 2024, 17(16), 4067; https://doi.org/10.3390/ma17164067
Submission received: 17 July 2024 / Revised: 13 August 2024 / Accepted: 14 August 2024 / Published: 16 August 2024
(This article belongs to the Section Advanced Composites)

Abstract

:
The incorporation of viscoelastic layers in laminates can markedly enhance the damped dynamic characteristics. This study focuses on integrating viscoelastic layers into the composite facesheet of the bowtie-shaped honeycomb core composite sandwich panel (BHC-CSP). The homogenization of the damped BHC-CSP is performed by employing the variational asymptotic method. Based on the generalized total energy equation, the energy functional of the representative unit cell of the damped BHC-CSP is asymptotically analyzed. The warping function, derived following the principle of minimum potential energy, provides a basis for obtaining the corresponding Euler–Lagrange equation to ascertain the equivalent elastic properties of the damped BHC-CSP. Utilizing the developed two-dimensional equivalent model, the free-vibration characteristics of the damped BHC-CSP are examined across diverse boundary conditions while delving into the impact of an external viscous damping layer on the natural frequency of the damped BHC-CSP. The results reveal that intensified boundary constraints effectively diminish the effective vibration region of the damped BHC-CSP, thereby enhancing its overall stability. The introduction of a PMI foam layer proves effective in adjusting the stiffness and mass distribution of the damped BHC-CSP. Resonance characteristics are explored through frequency and time-domain analyses, highlighting the pivotal roles of the excitation position and receiver point in influencing the displacement and velocity responses. Although the stiffness is improved by incorporating a PMI foam layer, its effect on the damping performance of the damped BHC-CSP is minimal when compared to the T-SW308 foam layer.

1. Introduction

The mechanic performance of sandwich structures is significantly influenced by both the structural configuration and the core material [1,2]. In recent years, there have been notable advancements in honeycomb sandwich panel design, including more intricate core shapes and innovative combinations with high-performance composite facesheets. Among different types of honeycomb sandwich panels, the re-entrant hexagonal honeycomb sandwich panel is notable for possessing a negative Poisson’s ratio, displaying unconventional behavior under loading, expanding under tension, and contracting under compression [3,4]. This unique mechanical response enhances the deformability and flexibility, leading to improved energy absorption and dissipation capabilities. Consequently, it exhibits an exceptional performance in vibration damping and impact resistance [5,6].
Numerous works have been dedicated to exploring the vibration characteristics of honeycomb sandwich panels. Jiang et al. [7] explored the free-vibration characteristics of composite sandwich plates with a re-entrant honeycomb core using theory, experiments, and the FE method. Sun et al. [8] conducted numerical and experimental investigations on the impact damage of aluminum honeycomb sandwich panels, focusing on the effect of interface layer geometry. Zhang et al. [9] explored the low-speed impact behavior of CFRP composite panels and Nomex honeycomb sandwich panels using numerical and experimental methods. Li et al. [10] investigated the anti-knock performance of hexagonal aluminum honeycomb sandwich panels through explosion testing. Current research on the dynamic performance of honeycomb sandwich panels primarily concentrates on specific structural parameters, material combinations, or restricted impact conditions, thereby constraining the generalizability and universality of research findings.
Re-entrant honeycombs typically exhibit high porosity with a relatively low effective stiffness, thereby restricting their practical applications in fields demanding high stiffness. To enhance the stiffness of regular re-entrant honeycombs, numerous novel variants have been devised, as illustrated in Figure 1. For instance, Lu et al. [11] introduced narrow ribs within the re-entrant cellular structure, observing a clear linear correlation between Poisson’s ratio and Young’s modulus with varying rib thickness. Chen et al. [12] incorporated strengthening ribs into traditional auxetic unit cells perpendicular to the re-entrant direction and noted a substantial enhancement of Young’s modulus by approximately 200% in the strengthening direction without significant compromise to auxecitical properties.
Fu et al. [13] introduced a novel honeycomb design by incorporating a rhombic pattern into the traditional re-entrant hexagonal honeycomb structure to improve its in-plane mechanical properties. Zhang et al. [14] created a novel re-entrant honeycomb structure by incorporating wedge-shaped components into the traditional design. Yu et al. [15] studied the mechanical properties of a novel foam-filled 2D re-entrant hexagonal honeycomb structure subjected to planar compression. Luo et al. [16] filled slow recovery foam into re-entrant honeycombs to enhance its mechanical characteristics. Qi et al. [17] introduced a new re-entrant circular honeycomb design that features double-circle arc cell walls instead of sloped walls. Zhang et al. [18] introduced a re-entrant arc-shaped honeycomb model based on the bio-inspired structural design.
Inspired by recent advancements in the field, a novel composite sandwich panel, called the bowtie-shaped honeycomb core composite sandwich panel (BHC-CSP), has been developed by integrating vertical walls into the core layer of a re-entrant hexagonal honeycomb sandwich panel, as depicted in Figure 2. This modification enhances the contact area between adjacent cells and enhances the connections within the cellular structure, leading to an improved overall performance. Nonetheless, enhancing the structural damping characteristics is imperative to attain efficient vibration mitigation and noise management in line with engineering specifications. Inserting viscoelastic layers in laminates can significantly improve the damped dynamic properties. Moreover, the concept of interlaminar damping is also well-suited to the manufacturing processes involved in creating laminated structures.
The literature offers limited studies on the vibration and damping characteristics of composite laminates with viscoelastic layers [19,20,21,22]. Cupial et al. [20] applied the first-order shear deformation theory to evaluate the loss factors of three-layered rectangular plates with a viscoelastic core layer and laminated facesheets. Plagianakos et al. [22] presented a higher-order theory and utilized finite element analysis to assess the modal damping of laminates with ply-angle variations and interply damping layers. Yang et al. [23] investigated the vibration and damping characteristics of hybrid carbon fiber composite pyramidal truss sandwich panels with viscoelastic layers integrated into the facesheets. Jean et al. [24] investigated the damping characteristics of unidirectional glass fiber composites with interleaved viscoelastic layers. The results indicate that hybrid composite sandwich structures with viscoelastic layers provide an effective method for improving damping performance.
Due to its innovation and potential applications, a multiscale study on the vibration and damping characteristics of the BHC-CSP is essential. Zhong and Yu [25,26,27] introduced a novel multiscale modeling approach for composite structures utilizing the variational asymptotic method (VAM) developed by Hodges and Cesnik [28]. On the microscopic scale, the effective properties of composite structures are obtained through asymptotic homogenization analysis and assigned to equivalent models. The strain energy of representative cells is obtained through homogenization analysis by applying the principle of minimum total potential energy, and the necessary performance parameters for macroscopic structure calculation (such as equivalent stiffness) are derived. Subsequently, the stress and strain responses of the macroscopic scale structure are obtained through macroscopic analysis. The introduction of a small parameter (such as the width-to-thickness ratio) allows for the gradual expansion of the energy functional, achieving a balance between calculation accuracy and efficiency [29,30,31].
To the best of our knowledge, the forced-vibration characteristics of the BHC-CSP, considering the damping layer, have not been reported. Therefore, in this work, a more precise and reliable VAM-based model is established to analyze the vibration characteristics of the damped BHC-CSP. The paper proceeds as follows: Section 2 and Section 3 present the theoretical formula and constitutive relationship of the two-dimensional equivalent Reissner-Mindlin model (2D-ERM) using the VAM. In Section 4, the forced-vibration analysis of the 2D-ERM considering a viscoelastic damping layer is conducted. Section 5 validates the effectiveness and precision of the 2D-ERM in analyses of free and forced vibrations of the damped BHC-CSP under different cases. In Section 6, the impact of geometric parameters on the vibration characteristics of the damped BHC-CSP is investigated. Section 7 compares the vibration characteristics of various damped honeycomb sandwich panels. Finally, Section 8 summarizes the primary conclusions of the study.

2. VAM-Based Homogenization of Composite Facesheets with Viscous Damping Layer

First, the VAM is employed to homogenize the laminated facesheet with the viscous damping layer, as shown in Figure 3d. The material kinematics are characterized by both the volume-averaged value and the corresponding difference:
u i ( y ; z ) = v i ( y ) + η w i ( y ; z ) ,
where u i and v i represent the displacements of the heterogeneous and homogeneous materials, respectively. y = y 1 , y 2 , y 3 and z = z 1 , z 2 , z 3 represent the meso- and microcoordinate systems, respectively. w i is the difference function signifying the variance between u i and v i ( i , j , k , l = 1 , 2 , 3 , and repeated subscripts imply summation).
The strain field within the heterogeneous material can be expressed as
ε i j ( y ; z ) = 1 2 u i ( y ; z ) y j + u j ( y ; z ) y i = ε ¯ i j + η w i j ,
where ε ¯ i j = 1 2 v i y j + v j y i represents the strain field of the homogeneous material, w i j = 1 2 w i z j + w j z i .
The kinematic variables of homogeneous materials can be expressed as
v i = 1 Ω Ω u i d Ω u i , ε ¯ ij = ε ij ,
where · represents the integration over the volume domain of the unit cell.
The constraint on the difference function can be obtained from Equation (3), such as
w i = 0 , w i j = 0 ,
The energy difference between the deformed inhomogeneous material and the homogeneous material can be specified as
δ J = δ Π micro Π macro λ k l w k l λ i w i = 0 ,
where Π micro = C i j k l ε i j ε k l = C i j k l ε ¯ i j + w i j ε ¯ k l + w k l ; Π macro = C ¯ i j k l c ¯ i j ε ¯ k l , C i j k l is the fourth-order elastic tensor; and λ k l and λ i are the Lagrange multipliers to introduce the constraints in Equation (4).
By minimizing the energy difference subject to the constraints in Equation (4), one obtains
min w i Eq · ( 4 ) 1 2 C i j k l ε i j ε k l = min w i Eq · ( 4 ) 1 2 C i j k l ε ¯ i j + w i j ε ¯ k l + w k l ,
The variational expression can be formulated as
U = 1 2 Γ h w + ϵ ¯ T C Γ h w + ϵ ¯ ,
where C is the 6 × 6 material matrix and w = w 1 , w 2 , w 3 T is the global strain array; ϵ ¯ = ε ¯ 11 ε ¯ 22 ε ¯ 33 2 ε ¯ 23 2 ε ¯ 13 2 ε ¯ 12 T , Γ h = / z 1 0 0 0 / z 3 / z 2 0 / z 2 0 / z 3 0 / z 3 0 0 / z 3 / z 2 / z 1 0 T .
The discrepancy function can be discretized utilizing shape functions that are defined on the unit cell:
w y i ; z i = S z i N y i ,
where the shape function S depends on the element type employed and N represents the nodal values of the difference function that needs to be solved.
Substituting Equation (8) into Equation (7), one can express the discrete form of the strain energy functional as
U = 1 2 N T D hh N + 2 N T D h ε e ¯ + ϵ ¯ T D ε ε ε ¯ ¯ 1 2 N T D ¯ N λ D h λ T N ,
where
D hh = Γ h S T C Γ h S , D h e = Γ h S T C , D c e = C , D h λ = Γ h S T .
Minimizing Equation (9) leads to the linear system
D hh N = D h s ϵ ¯ .
It is evident that N is directly proportional to ϵ , and the solution can be expressed as
N = N 0 ϵ ¯ , w = S N 0 ϵ ¯ .
By substituting Equation (12) into Equation (9), one can calculate the energy stored in the unit cell as
U = 1 2 ϵ ¯ T 2 N 0 T D h ε + N 0 T D h h N 0 + D ε E ϵ ¯ = Ω 2 ϵ ¯ T C ¯ ϵ ¯ ,
where C ¯ is the equivalent material matrix.

3. VAM-Based Equivalent Reissner-Mindlin Model for BHC-CSP

The process of establishing the 2D-ERM of the BHC-CSP using the VAM is depicted in Figure 3b. The VAM-based equivalent Reissner-Mindlin model involves representing the displacement field ( u i ) of the BHC-CSP using the displacements ( u ¯ i ) in the 2D-ERM and warping functions w i , such as
u 1 x 1 , x 2 , y 1 , y 2 , y 3 , t = u ¯ 1 x 1 , x 2 , t ζ y 3 u ¯ 3 , 1 x 1 , x 2 , t + ζ w 1 x 1 , x 2 , y 1 , y 2 , y 3 , t , u 2 x 1 , x 2 , y 1 , y 2 , y 3 , t = u ¯ 2 x 1 , x 2 , t ζ y 3 u ¯ 3 , 2 x 1 , x 2 , t + ζ w 2 x 1 , x 2 , y 1 , y 2 , y 3 , t , u ¯ 3 x 1 , x 2 , y 1 , y 2 , y 3 , t = u ¯ 3 x 1 , x 2 , t + ζ w 3 x 1 , x 2 , y 1 , y 2 , y 3 , t .
The explicit representation of u ¯ i can be derived from Equation (14),
u ¯ 1 = u 1 + ζ y 3 u ¯ 3 , 1 , u ¯ 2 = u 2 + ζ y 3 u ¯ 3 , 2 , u ¯ 3 = u 3 ,
where · represents the volume average over the unit cell.
Since the meso-coordinate y i originates at the geometric center of the unit cell, it follows that y 3 = 0 . This results in three constraints on the warping functions:
w i = 0 .
The concept of rotation tensor decomposition can be used to express the 3D strain components with small local rotation as
ε i j = 1 2 u i x j + u j x i .
The 3D strain field Γ can be represented in matrix form as
Γ e = ε 11 ε 22 2 ε 12 T = ϵ + ζ y 3 κ + I α w , α , 2 Γ s = 2 ε 13 2 ε 23 T = w , 3 + e α w 3 , α , Γ t = ε 33 = w 3 , 3 ,
where ( ) | | = ( ) 1 ( ) 2 T , ϵ = ϵ 11 2 ϵ 12 ϵ 22 T , κ = κ 11 κ 12 + κ 21 κ 22 T , and
I 1 = 1 0 0 1 0 0 ,   I 2 = 0 0 1 0 0 1 ,   e 1 = 1 0 ,   e 2 = 0 1 .
The panel’s stored strain energy can be represented as
U = s 1 Ω U Ω d x 1 d x 2 ,
where s is the reference surface, with an integral domain of a / 2 < x 1 < a / 2 , b / 2 < x 2 < b / 2 , where a and b are the length and width of the sandwich panel, respectively; U Ω / Ω represents the strain energy density.
Based on the principle of conservation energy, the strain energy is equivalent to the total energy stored by all its components:
U Ω = V b Γ b T D b Γ b d y 1 d y 2 d y 3 + V t Γ t T D t Γ t d y 1 d y 2 d y 3 + 4 × V A Γ A T D A Γ A d y 1 d y 2 d y 3 + V B Γ B T D B Γ B d y 1 d y 2 d y 3 + V c Γ C T D C Γ C d y 1 d y 2 d y 3 + 2 × V D Γ D T D D Γ D d y 1 d y 2 d y 3 ,
where subscripts b and t represent the bottom and top facesheets, respectively; A , B , C , and D represent different parts within the 1/4 core, as shown in Figure 4.
Since the coordinates originate from the geometric center of the core cell, the integration domains of the top and bottom facesheets in Equation (21) are identical, with the exception of V b , h c / 2 < y 3 < h c / 2 h f and V t , h c / 2 < y 3 < h c / 2 + h f . In the same way, the integral domain of the core layer is the same along the y 3 axis, and the differences are
V A , 0 < y 1 < l 1 / 2 + t c , l 2 / 2 t c < y 2 < l 2 / 2 V B , l 1 / 2 + l 2 / 2 l 3 / tan α < y 1 < l 1 / 2 , l 3 < y 2 < l 2 / 2 l 3 V C , l 1 / 2 + l 2 / 2 l 3 / tan α t c < y 1 < l 4 / 2 , 0 < y 2 < t c V D , l 1 / 2 + t c < y 1 < l 1 / 2 , l 2 / 2 l 3 < y 2 < l 2 / 2
Equation (21) can be written as
U = 1 2 Γ T D Γ = 1 2 Γ e 2 Γ s Γ t T C e C e s C e t C e s T C s C s t C e t T C s t T C t Γ e 2 Γ s Γ t ,
where C e , C e s , C e t , C s , C s t , and C t represent the respective submatrices of the 3D 6 × 6 material matrix.
The kinetic energy of the BHC-CSP can be obtained by
K = 1 2 V ρ v T v d V = K 2 D + K ,
where ρ denotes the mass density, v represents the absolute velocity of a generic point in the panel, and
K 2 D = 1 2 Ω μ ¯ V T V + 2 ω T μ ξ ¯ ˜ V + ω T j ω d Ω ,
K = 1 2 V ρ ( ω ˜ w + w ˙ ) T ( ω ˜ w + w ˙ ) + 2 ( V + ω ˜ ξ ) T ( ω ˜ w + w ˙ ) d V
where μ ¯ , μ ξ ¯ , and j are inertial constants commonly used in plate dynamics.
The virtual work resulting from applied loads can be calculated as
δ W ¯ = δ W ¯ 2 D + δ W ¯ ,
where
δ W ¯ 2 D = s f i δ q ¯ i + m α δ ψ ¯ α d s + s Q i δ q ¯ i + x 3 Q α δ ψ ¯ α d s ,
δ W ¯ = s P i δ w i + τ i δ w i + + β i δ w i d s + s Q i δ w i d s ,
with m α and f i denoting the generalized moments and forces, respectively; δ ψ ¯ i and δ q ¯ i are virtual rotation and displacement, respectively.
For simplicity, the virtual work conducted by the warping function along the lateral boundaries of the panel can be neglected to derive the extended Hamilton’s principle as
t 1 t 2 δ K 2 D + K U + δ W ¯ 2 D + δ W ¯ d t = 0 ,

3.1. First Approximation

Utilizing the VAM allows one to derive the leading terms in Equation (30) by disregarding smaller terms, such as
t 1 t 2 δ K 2 D Ω U 0 d Ω + δ W ¯ 2 D d t = 0 ,
where the first approximation of the strain energy can be obtained by eliminating smaller energy contributions due to w i , α , such as
2 U 0 = ϵ + ζ y 3 κ T C e ϵ + ζ y 3 κ + 2 ϵ + ζ y 3 κ T C e s w , 3 + C e t w 3 , 3 + w , 3 T C s w , 3 + 2 w , 3 T C s t w 3 , 3 + w 3 , 3 T C t w 3 , 3 .
Introducing a Lagrange multiplier λ i allows one to derive the associated Euler–Lagrange equation as
ϵ + ζ y 3 κ T C e s + w , 3 T C s + w , 3 T C s t , 3 = λ , ϵ + ζ y 3 κ T C e t + w , 3 T C s t + w 3 , 3 C t , 3 = λ 3 ,
where λ | | = λ 1 λ 2 T .
The free boundary conditions at the top and bottom surfaces of the sandwich panel can be determined by
ϵ + ζ y 3 κ T C e s + w , 3 T C s + w , 3 T C s t + / = 0 , ϵ + ζ y 3 κ T C e t + w , 3 T C s t + w 3 , 3 C t + / = 0 ,
where the superscript “+/−” indicates the quantity on the top and bottom surfaces of the panel, respectively.
Given these conditions, one can represent the solutions for w | | and w 3 as
w = ϵ + ζ y 3 κ C ¯ e s C s 1 T , w 3 = ϵ + ζ y 3 κ C ¯ e t C ¯ t 1 ,
where
C ¯ e s = C e s C ¯ e t C s t T C ¯ t 1 , C ¯ e t = C e t C e s C s 1 C s t , C ¯ t = C t C s t T C s 1 C s t .
The stain energy of the 2D-ERM can be determined by substituting Equation (35) in Equation (32):
U 2 D = 1 2 ϵ + ζ y 3 κ T K ϵ + ζ y 3 κ = 1 2 ϵ κ T A B B D ϵ κ ,
with
A = K , B = ζ y 3 K , D = ζ y 3 2 K , K = C e C ¯ e s C s 1 C e s T C ¯ e t C e t T / C ¯ t .

3.2. Second Approximation

The first approximation aligns with classical plate theory, allowing for the calculation of in-plane stresses. However, to accurately consider out-of-plane stresses, the second approximation is necessary. This involves perturbing the warping function as
w = v ¯ , w 3 = v ¯ 3 + D w ,
where w = [ ε κ ] T , D = D e t T D t x 3 D e t T D t .
Substituting Equation (39) into Equation (31) and then into Equation (32), the expression for the second approximate energy can be obtained as
2 Π 1 = v ¯ , 3 T D s v ¯ , 3 + D t v ¯ 3 , 3 2 + 2 v ¯ T C , 3 w , α + 2 v ¯ T D s t D w , α 2 v ¯ T p 2 v ¯ T τ 2 v ¯ T β .
The related Euler–Lagrange equation can be derived as
D s v ¯ , 3 + D s t D w , α , 3 = C , 3 t D w , α + g , 3 + λ ,
where C = e T D x 3 D , g , 3 = p .
Since v ¯ 3 does not depend on v ¯ , v ¯ 3 has a trivial solution. Consequently, the solution for v ¯ can be determined as
v ¯ = C ¯ + L α w , α + g ¯ ,
where
C ¯ , 3 = D s 1 C , C ¯ = 0 , g ¯ 3 = D s 1 g ¯ , g ¯ = 0 , L α w , α = C ¯ / h , C ¯ = C + x 3 h D α 1 2 D ± D s e α D , g ¯ = g + x 3 h g 1 2 g ± ,
and ( · ) = ( · ) ( · ) + , ( · ) ± = ( · ) + + ( · ) .
The second-order approximate energy can be formulated in the form of the Reissner-Mindlin model as
2 Π 1 = w T A w + w , α T B α β w , β 2 w T F ,
where B α β is additional shear stiffness and F is a load-related item:
A = D x 3 D x 3 D x 3 2 D , B α β = D S α β D T D C ¯ α T D s 1 C ¯ β + L α T C β , 3 , F = D T p 3 C ¯ T D s 1 g ¯ , a L α p ¯ + p , α .
The constitutive relationship for the 2D-ERM can be derived as
N 11 N 22 N 12 M 11 M 22 M 12 Q 1 Q 2 = A 11 A 12 A 16 B 11 B 12 B 16 Y 11 Y 12 A 12 A 22 A 26 B 12 B 22 B 26 Y 21 Y 22 A 16 A 26 A 66 B 16 B 26 B 66 Y 31 Y 32 B 11 B 12 B 16 D 11 D 12 D 16 Y 41 Y 42 B 12 B 22 B 26 D 12 D 22 D 26 Y 51 Y 52 B 16 B 26 B 66 D 16 D 26 D 66 Y 61 Y 62 Y 11 Y 21 Y 31 Y 41 Y 51 Y 61 G 11 G 12 Y 12 Y 22 Y 32 Y 42 Y 22 Y 62 G 12 G 22 ϵ 11 ϵ 22 2 ϵ 12 κ 11 κ 22 2 κ 12 2 γ 13 2 γ 23
where N α β , M α β , and Q α are the in-plane force, moment, and transverse shear stress resultants, respectively.

4. Free- and Forced-Vibration Analysis of the Damped 2D-ERM

The free-vibration equation of the 2D-ERM is
D 1 2 ψ 1 x 1 2 + 2 ψ 1 x 2 2 + D 3 2 ψ 2 x 1 x 2 + C 1 u ¯ 3 x 1 ψ 1 + γ 4 ψ 1 = 0
D 2 2 ψ 2 x 2 2 + 2 ψ y x 2 + D 3 2 ψ 1 x 1 x 2 + C 2 u ¯ 3 x 2 ψ 2 + γ 4 ψ 2 = 0
C 1 2 u ¯ 3 x 1 2 ψ 1 x 1 + C 2 2 u ¯ 3 x 2 2 ψ 2 x 2 + β 4 u ¯ 3 = 0
where
D 1 = D 11 D 66 , D 2 = D 22 D 66 , D 3 = D 12 + D 66 D 66 , C 1 = C 55 D 66 , C 2 = C 44 D 66 , γ 4 = ρ J ω 2 D 66 , β 4 = ρ h ω 2 D 66
Differentiating Equation (49) with respect to x 1 and x 2 , and substituting the results into Equations (47) and (48), respectively, to eliminate ψ 1 and ψ 2 , one obtains
1 1 C 1 D 1 2 x 1 2 D 3 C 1 C 2 2 x 1 2 + 2 x 2 2 + γ 4 ψ 1 = 1 + D 3 C 1 2 x 2 2 + C 1 C 2 2 x 1 2 + β 4 C 2 u ¯ 3 x 1
1 1 C 2 D 2 2 x 2 2 D 3 C 2 C 1 2 x 2 2 + 2 x 1 2 + γ 4 ψ 2 = 1 + D 3 C 2 2 x 1 2 + C 2 C 1 2 x 2 2 + β 4 C 1 u ¯ 3 x 2
Equation (49) can be rewritten as
C 1 2 x 1 2 + C 2 2 x 2 2 + β 4 u ¯ 3 = C 1 ψ 1 x 1 + C 2 ψ 2 x 2
The damped forced-vibration equilibrium equation for the 2D-ERM is
M 11 x 1 M 12 x 2 + Q 1 ρ J 2 ψ 1 t 2 = 0 ,
M 12 x 1 M 22 x 2 + Q 2 ρ J 2 ψ 2 t 2 = 0 ,
Q 1 x 1 + Q 2 x 2 + p = ρ h 2 u ¯ 3 t 2 + c u ¯ 3 t .
where the load p is only a function of time t, and c is the damping coefficient.
By substituting M 11 ,   M 22 ,   M 12 ,   Q 1 , and Q 2 into Equations (54)–(56), one obtains the differential equation of damped forced vibration as
D 11 2 ψ 1 x 1 2 + D 66 2 ψ 1 x 2 2 + D 12 + D 66 2 ψ 2 x 1 x 2 + C 55 u ¯ 3 x 1 ψ 1 ρ J 2 ψ 1 t 2 = 0 ,
D 22 2 ψ 2 x 2 2 + D 66 2 ψ 2 x 1 2 + D 12 + D 66 2 ψ 1 x 1 x 2 + C 44 u ¯ 3 x 2 ψ 2 ρ J 2 ψ 2 t 2 = 0 ,
C 55 2 u ¯ 3 x 1 2 ψ 1 x 1 + C 44 2 u ¯ 3 x 2 2 ψ 2 x 2 + p = ρ h 2 u ¯ 3 t 2 + c u ¯ 3 t .
where
C 44 = κ G 11 h , C 55 = κ G 22 h .
Based on the vibration modes and natural frequencies obtained from free-vibration analysis, the locations of receiver and excitation points under dynamic loading, such as harmonic loading, are determined based on the method outlined in Ref. [32]. Specifically, the obtained effective properties of the damped BHC-CSP are incorporated into the equivalent model through the ’general shell’ feature in the S4R element, serving as an enhancement to the commercial program. Moreover, the ’frequency’ procedure in ABAQUS software (version 6.14) is utilized to determine the natural frequencies and vibration modes. Additionally, the modal/steady-state dynamics solver in ABAQUS software is employed for time- and frequency-domain forced-vibration analysis, respectively.

5. Free- and Forced-Vibration Analysis of Damped BHC-CSP

The dynamic performance of the damped BHC-CSP, particularly regarding free- and forced-vibration characteristics, plays a pivotal role in ensuring structural reliability and durability. The addition of damping layers to the sandwich panel introduces significant complexity, rendering analytical solutions and experimental data from the existing literature unavailable. Therefore, this section relies on numerical simulations as a valuable substitute to investigate the free and forced vibration of the damped BHC-CSP, providing insights that may inform future experimental endeavors and potentially result in substantial time and resource savings. Six typical boundary conditions (CCCC, CCCS, CCSS, SSSS, CCFF, and SSSF), as outlined in Figure 5, are used to investigate its free-vibration characteristics, with the letters C, S, and F denoting the clamped, simply supported, and free sides, respectively.

5.1. Model Information

Figure 6a illustrates the core cell of the BHC-CSP, featuring the following structural specifications: t c = 0.8 mm, l 1 = 16.8 mm, l 2 = 16 mm, l 3 = 2.4 mm, and l 4 = 24 mm. In addition, the core height is h c = 6 mm while the face panel thickness is h f = 0.6 mm. The BHC-CSP consists of 15 unit cells in both x 1 and x 2 directions, with dimensions of L = 360 mm in length and W = 240 mm in width. The top and bottom facesheets are composed of CFRP material with a layup configuration of [45/−45/0/90]s and a fiber volume fraction of 60%. The core layer is constructed from aluminum metal. Detailed properties of each material can be referenced in Table 1.
The viscoelastic materials are Polyimide (PMI) foam with a higher modulus and ternary ethylene propylene rubber foam (T-SW308) with a lower modulus, with their material properties detailed in Table 1. It can be anticipated that the developed method is suitable for other damping materials with differing moduli. A previous study indicates that incorporating viscoelastic damping layers on the facesheet and core layer individually had similar effects on the overall structural stiffness. However, positioning the damping layer in the facesheet led to superior enhancement of the damping performance of the sandwich panel. This superiority can be attributed to the higher strains experienced by the top and bottom facesheets of the sandwich panel during torsional or bending deformations. Therefore, the viscoelastic damping layer is integrated into the facesheet, as illustrated in Figure 7. It is mentioned that the 3D-FEM utilizes 340,653 C3D10 elements and 600,301 nodes, employing a structured meshing approach with a mesh size set to 0.2 mm after a mesh convergence study.

5.2. Free-Vibration Analysis of BHC-CSP with Viscoelastic Damping Layer

This section investigates the impact of PMI damping layers on the natural frequencies of composite sandwich panels. Figure 8 compares the first six natural frequencies of the BHC-CSP with varying thicknesses of the PMI damping layers under six different cases, as predicted by the 2D-ERM and 3D-FEM. Specifically, the red dots correspond to the 3D-FEM data, while the blue dots correspond to the 2D-ERM data. The results indicate that the VAM-based 2D-ERM accurately captures the impact of the viscoelastic damping layer on the natural frequencies. The natural frequencies exhibit a gradual increase with variations in the PMI foam thickness from 0.1 to 0.6 mm. Notably, the increase in higher-order natural frequencies is proportionally comparable to that of the lower-order natural frequencies. This behavior stems from PMI foam with a higher modulus; the addition of a PMI layer effectively increases the panel thickness, thereby augmenting the overall structural stiffness. Despite the rise in the sandwich panel’s total mass due to the inclusion of PMI foam, this increase in mass fails to counterbalance the amplified natural frequencies resulting from enhanced stiffness. Particularly when the PMI damping layer is thin, the effect of stiffness enhancement on natural frequencies becomes more discernible.
In addition, the relatively uniform distribution of PMI foam within the facesheet of the sandwich panel results in a uniform effect on the stiffness and mass of the entire panel, thereby exerting comparable influences on both higher- and lower-order frequencies. To sum up, integrating a constrained viscoelastic damping layer allows for a precise adjustment of the stiffness and mass distribution within the sandwich panel, laying the groundwork for the tailored design of panels with distinct vibration characteristics.

5.3. Frequency-Domain Analysis of BHC-CSP with Viscoelastic Damping Layer

Forced-vibration analysis assesses the panel’s response to external dynamic loads, a critical step in evaluating and optimizing the panel’s performance in its operational environment. This section investigates the impact of various viscoelastic damping layers on the resonant response of the BHC-CSP. The modal dynamic analysis involves single-point excitation and output to capture additional resonance responses for the first six frequencies under SSSS BCs.
The excitation force, set as a harmonic load with a 10 kN amplitude, sweeps across frequencies ranging from 0 to 3000 Hz, with a damping ratio of 0.02 for each mode. To explore the impact of excitation location on the structural frequency-domain response, three receiving points—A, B, and C—at coordinates (45, 60), (90, 60), and (180, 60), respectively, are selected for investigation, as shown in Figure 9. Through the steady-state modal dynamic analysis step in the ABAQUS software, a frequency-domain response analysis is conducted on the sandwich panel with or without these viscoelastic layers.
Figure 10a,b compare the frequency-domain response predicted by the 3D-FEM and 2D-ERM under Cases 8 and 9, where points B and C are excited. The frequency-domain response curves are obtained at three distinct receiving points for six models (BM, BM-2D, PMI, PMI-2D, T-SW308, and T-SW308-2D). In this context, BM, PMI, and T-SW308 represent the BHC-CSP without the viscoelastic damping layer and with the PMI and T-SW308 viscoelastic layers, respectively. On the other hand, BM-2D, PMI-2D, and T-SW308-2D represent the 2D equivalent models of the corresponding sandwich panels.
It is evident that the equivalent model is applicable for studying the frequency-domain response of the damped BHC-CSP. Specifically, the PMI-2D and T-SW308-2D models can effectively predict the resonant response of the corresponding viscoelastic-damped sandwich panels while meeting engineering accuracy requirements. An examination of the first-order resonance amplitudes reveals that the addition of the PMI foam viscoelastic layer has an insignificant impact on the panel’s resonance amplitude, while incorporating the T-SW308 foam viscoelastic layer effectively reduces the panel’s resonance amplitude. Adding the PMI foam viscoelastic layer results in an increase of 19.98% in the fundamental frequency, whereas adding the T-SW308 foam viscoelastic layer leads to a 5.00% decrease in the fundamental frequency. This is attributed to the soft texture of T-SW308 foam, which to some extent reduces the panel’s stiffness, consequently lowering the natural frequency.
In summary, a comprehensive comparison shows that the PMI damping layer can effectively increase the panel’s stiffness with minimal impact on the damping performance, while the T-SW308 damping layer can significantly enhance the damping performance without significantly altering the stiffness of the sandwich panel. This indicates that T-SW308 foam is an excellent damping layer material with a superior damping performance.

5.4. Time-Domain Analysis of BHC-CSP with Viscoelastic Damping Layer

The time-domain responses of the sandwich panel with the PMI damping layer, T-SW308 damping layer, and without a viscoelastic damping layer are examined by applying excitations of F = 10 sin 2 π f i t at points B and C, as shown in Figure 11. The results indicate that the inclusion of the T-SW308 damping layer leads to a significant reduction in the displacement and velocity amplitudes of the sandwich panel. This reduction suggests a notable decrease in the overall system energy, indicating that the viscoelastic damping layer effectively dampens vibration amplitudes, reducing the likelihood of resonant amplification and structural damage.
Moreover, following the addition of the viscoelastic damping layer, the displacement–time and velocity–time histories exhibit a tendency to stabilize more rapidly. This illustrates that the viscoelastic damping layer aids in dissipating vibrational energy more efficiently, enabling the system to achieve an equilibrium state in a shorter duration. However, excessive damping can result in the over-dissipation of vibrational energy, potentially slowing down the free decay process of vibrations. Consequently, in engineering practice, it is critical to comprehensively account for other factors when enhancing the damping performance. This approach ensures a balanced consideration of both stability and sensitivity in the panel, ultimately facilitating an optimal system performance in vibrational settings.

5.5. Comparison of Dynamic Analysis Efficiency

The numerical simulations were carried out on the Dell Precision 3660 workstation. Table 2 delineates a thorough comparison of the computational efficiency between the two models utilized for dynamic analysis. It illustrates that the time needed for free vibration and forced vibration in the frequency and time domain using the 2D-ERM is notably lower in comparison to the 3D-FEM, showcasing reductions of 1.9%, 2.6%, and 1.6%, respectively. These results emphasize the significant improvement in dynamic analysis efficiency through the VAM-based 2D-ERM, presenting a more efficient tool for dynamic performance evaluation in engineering fields.
The comparison highlights the superior efficiency and accuracy of the 2D-ERM in quantitatively assessing the free- and forced-vibration behavior of the BHC-CSP with the viscoelastic damping layer. However, in three-dimensional models with intricate geometric and physical characteristics, the limitations of the 2D-ERM become apparent, as it fails to capture the detailed geometric intricacies of complex structures as comprehensively as the 3D-FEM. Consequently, when selecting a model for dynamic analysis, it is essential to strike a balance between model accuracy and computational efficiency.

6. Parameters-Influenced Analysis

This section examines the influence of four crucial geometric parameters on the natural frequency, equivalent density, and displacement–time history of the damped BHC-CSP utilizing the 2D-ERM. The alterations in geometric parameters are detailed in Table 3 for reference.

6.1. Height Ratio ( h c / h f )

Figure 12 compares the trends in the variation in the first four natural frequencies, equivalent densities, and displacement–time history of the damped BHC-CSP with different height ratios. The examination of natural frequencies illustrates a nonlinear rise in the first four natural frequencies as the core layer height increases, with the growth rate gradually declining. This phenomenon stems from a more pronounced escalation in equivalent stiffness relative to equivalent density, resulting in a significant upsurge in the natural frequency. The declining growth rate indicates reduced benefits of further heightening the core layer beyond a certain threshold.
Figure 12b indicates that with an increase in the core layer height, the panel exhibits decreased sensitivity to external excitation, leading to a gradual decline in the resonance amplitude. This decline is attributed to the marked augmentation in overall structural stiffness with an increased core layer height, resulting in a more dispersed vibration energy and amplified damping effects, consequently reducing the resonance amplitude. These findings underscore the effectiveness of adjusting the core layer height in enhancing the resonance response of sandwich panels.

6.2. Wall Thickness-to-Height Ratio ( t c / h c )

Figure 13 compares the trends in the variations in the first four natural frequencies, equivalent densities, and displacement–time history curves of the damped BHC-CSP with different wall thickness-to-height ratios. An examination of the natural frequencies indicates that as the wall thickness increases, there is a decreasing trend in natural frequency, with a slight mutation observed at the ratio of 0.25. This decrease is attributed to a sudden increase in equivalent stiffness at this ratio, which surpasses the increase in equivalent density, resulting in an increase in natural frequency (the equivalent density, equivalent tensile stiffness, and equivalent bending stiffness increase by 15%, 105%, and 66%, respectively). Figure 13b illustrates that the resonance amplitude gradually decreases with an increase in the cell wall thickness. This reduction can be attributed to the corresponding increase in the structure’s equivalent stiffness, particularly the shear stiffness and torsional stiffness, both crucial for resisting shear and torsional deformations, which commonly constitute significant components of the vibration modes.

6.3. Re-Entrant Angle ( θ )

Figure 14 compares the trends in the variations in the first four natural frequencies, equivalent densities, and displacement–time history curves of the damped BHC-CSP with different re-entrant angles. An examination of the natural frequencies reveals a gradual decrease in equivalent density as the re-entrant angle increases, while the changes in natural frequency align with the variations in directional equivalent stiffness. An investigation of the resonance response in Figure 14b indicates that the resonance amplitude of the panel is more sensitive to changes in directional equivalent stiffness. A detailed comparison of resonance response and natural frequencies indicates that the panel demonstrates higher natural frequencies and relatively smaller resonance amplitudes at α = 40 , mitigating the risk of external excitation-induced resonance failure.

6.4. Aspect Ratio ( l 4 / l 2 )

Figure 15 compares the trends in the variations in the first four natural frequencies, equivalent densities, and displacement–time history curves of the damped BHC-CSP with different aspect ratios. An examination of the natural frequencies unveils that the first four natural frequencies rise as the aspect ratio increases. This increase is attributed to the smaller decrease in structural equivalent stiffness (8.88%) compared to the decrease in equivalent density (14.09%). An examination of the resonance response in Figure 15b reveals that the resonance amplitude of the panel escalates with the aspect ratio increment. This phenomenon can be explained by the panel becoming more flexible as the overall stiffness decreases, allowing for increased energy absorption during vibration and consequently leading to a larger dynamic response during resonance. Notably, as the l 4 / l 2 ratio shifts from 1.1 to 1.6, the fundamental frequency of the panel increases by 5.72%, while the resonance amplitude experiences a 4.24% increase. These findings underscore the importance of considering the aspect ratio’s influence on the dynamic performance of the damped BHC-CSP.

6.5. Summary

To visually demonstrate the influence of geometric parameters on the dynamic performance of the damped BHC-CSP, Table 4 lists the unit changes in the equivalent density ( Δ ρ ), fundamental frequency ( Δ ω ), resonance amplitude ( Δ R A ), and equivalent tensile and bending stiffness ( Δ A 22 and Δ D 22 , respectively) with the variation in different parameters. It is evident that the core layer height exerts the most significant influence on the equivalent bending stiffness, fundamental frequency, and resonance amplitude. This observation suggests that adjusting the core layer height can effectively enhance the panel’s dynamic response to external excitations. Moreover, variations in the re-entrant angle, which impact the volume of the core layer per unit length, play a crucial role in adjusting the panel’s tensile performance. By adjusting the equivalent stiffness, the panel’s response efficiency to external excitations can be improved, leading to reduced vibration amplitudes and the effective mitigation of structural fatigue and resonance damage risks. Therefore, optimizing the equivalent stiffness and overall dynamic response can be achieved through adjustments to the re-entrant angle of the bowtie-shaped honeycomb.

7. Comparison of Vibration Characteristics between Different Honeycomb Sandwich Panels

Distinct geometric shapes of various honeycomb structures result in unique energy distribution, load transfer paths, and deformation modes. Analyzing the differences in the vibration characteristics of diverse honeycomb sandwich panels aids in enhancing the structural design rationale. This, in turn, facilitates the development of honeycomb structures with superior adaptability and mechanical properties.
In this section, apart from the BHC-CSP, three typical honeycomb sandwich panels are investigated, including the conventional hexagonal core sandwich panel (CHC-CSP), re-entrant hexagonal core sandwich panel (RHC-CSP), and star-shaped honeycomb core sandwich panel (SHC-CSP), as illustrated in Figure 16. The geometric parameters are as follows: CHC-CSP, l 1 = 7.6 mm ,   α = 120 ; RHC-CSP, l 1 = 18.3 mm ,   α = 60 , l 3 = 8.32 mm ; and SHC-CSP, α = 30 ,   l 1 = 12.6 mm . The geometric parameters for the bowtie honeycomb remain consistent with those outlined in Section 5.1. The layup mode of CFRP is [±45/0/90]s, with a fiber volume fraction of 60%, while the core layer material remains unchanged. The core layer weights for the hexagonal honeycomb, re-entrant honeycomb, bowtie-shaped honeycomb, and star-shaped honeycomb were 1.29 g, 1.27 g, 0.84 g, and 2.03 g, respectively. The core volumes were 477.25 mm3, 469.93 mm3, 310.57 mm3, and 750.54 mm3, respectively.

7.1. Comparison of Free-Vibration Characteristics

Figure 17a–f compare the first six natural frequencies of four damped honeycomb sandwich panels under SSSS BCs predicted by the 3D-FEM and 2D-ERM. It is observed that the VAM-based equivalent model can accurately predict the vibration modes of the respective damped honeycomb sandwich panel. While the relative errors of natural frequencies predicted by the equivalent model increase with the modal order, the relative errors fall between 0.701% and 4.576%. This finding confirms the precision of utilizing the VAM for analyzing the free-vibration characteristics of damped honeycomb sandwich panels with distinct honeycomb cores. Comparing the natural frequencies of the different damped honeycomb sandwich panels reveals that the fundamental frequency of the CHC-CSP, SHC-CSP, and BHC-CSP increases by 14.286%, 12.936%, and 1.209%, respectively, relative to the RHC-CSP. This suggests that the RHC-CSP is more susceptible to external excitation and resonance. In contrast, the CHC-CSP and SHC-CSP exhibit lower sensitivity to external excitation and demonstrate a more stable performance.

7.2. Comparison of Forced-Vibration Characteristics

The displacement–time and velocity–time history curves at point C of different damped honeycomb sandwich panels (the CHC-CSP, RHC-CSP, and BHC-CSP) were compared in Figure 18. Upon analyzing the displacement–time history curves, it is evident that the steady-state maximum displacement amplitudes of the CHC-CSP, RHC-CSP, and BHC-CSP in relation to the SHC-CSP increased by 30.719%, 19.350%, and 6.025%, respectively. This suggests that the CHC-CSP exhibits the highest increase in displacement amplitude, followed by the RHC-CSP and then the BHC-CSP.
Similarly, the steady-state maximum velocity amplitudes of the CHC-CSP, RHC-CSP, and BHC-CSP relative to the SHC-CSP increased by 39.627%, 27.733%, and 8.683%, respectively. This indicates that the CHC-CSP also shows the largest increase in velocity amplitude, followed by the RHC-CSP and then the BHC-CSP. These findings further elucidate the differences in vibration characteristics among the different damped honeycomb sandwich panels, with the CHC-CSP showcasing the highest displacement and velocity amplitudes, followed by the RHC-CSP and BHC-CSP.
It is important to highlight that the displacement–time history predicted by the VAM-based 2D-ERM aligns closely with those derived from the corresponding 3D-FEM, meeting the required engineering precision. The CHC-CSP is relatively insensitive to external excitation. However, it exhibits the largest displacement and velocity amplitudes during resonance, making the panel highly susceptible to damage induced by resonance. On the other hand, the BHC-CSP not only displays low sensitivity to external excitation but also demonstrates smaller displacement and velocity amplitudes during resonance. This suggests that the bowtie-shaped honeycomb design can effectively mitigate resonance effects and possesses a favorable anti-vibration performance.
In summary, while the CHC-CSP shows resilience to external excitations, its susceptibility to resonance-induced damage is a concern. In contrast, the BHC-CSP not only resists external excitations but also minimizes the impact of resonance, highlighting its superior anti-vibration capabilities.

8. Conclusions

This work aims to utilize the VAM to decompose the BHC-CSP with a viscoelastic damping layer into a unit-cell-based constitutive model and two-dimensional equivalent Reissner-Mindlin model (2D-ERM), enabling the analysis of the free- and forced-vibration characteristics of the damped BHC-CSP. The efficiency and accuracy of the 2D-ERM were assessed by comparing the results from the 3D-FEM. The key conclusions drawn from the study are as follows:
(1)
Reinforcing the boundary constraints diminishes the effective vibration region of the panel, necessitating higher frequencies to attain equivalent amplitudes. The increase in the panel’s mass resulting from the inclusion of a PMI damping layer is less pronounced compared to the augmentation in stiffness. Given the relatively uniform distribution of the PMI damping layer within the panel, its impact on both high-order and low-order vibration modes is similar. The maximum relative error in resonance amplitude prediction by the 2D-ERM stands at 8.204%, highlighting the model’s capability to accurately predict the frequency-domain response of the panel. In contrast to the PMI damping layer, the T-SW308 damping layer significantly bolsters the panel’s damping efficiency while maintaining a relatively low stiffness level. An analysis of the time-domain data reveals the panel’s relatively low sensitivity to external excitations within the frequency band of 150∼300 Hz.
(2)
An appreciable increase in the equivalent stiffness of the panel occurs with the elevation of the core layer height, where the rise in equivalent stiffness surpasses the increase in equivalent density, consequently leading to a notable enhancement in natural frequencies. Given that shear and torsional deformation modes significantly contribute to the vibration analysis, the overall structural vibration response is particularly responsive to alterations in shear and torsional stiffness. At a re-entrant angle of 40°, the sandwich panel demonstrates a heightened natural frequency and comparatively reduced resonance amplitude. Notably, adjustments in the aspect ratio exhibit minimal impact relative to other geometric parameters on the natural frequency.
(3)
The hexagonal honeycomb sandwich panel and star-shaped honeycomb sandwich panel exhibit reduced sensitivity to external excitation. Furthermore, the star-shaped honeycomb sandwich panel not only displays lower sensitivity to external excitation but also features smaller displacement and velocity amplitudes during resonance. The 2D-ERM requires only a fraction of the time—1.9%, 2.6%, and 1.6%—compared to what the 3D-FEM necessitates for free-vibration and frequency- and time-domain analysis, respectively. Therefore, when dealing with three-dimensional spatial models that exhibit substantial geometric and physical complexity and necessitate a trade-off between model precision and computational efficiency, the VAM-based 2D equivalent model emerges as a more effective tool for dynamic analysis.

Author Contributions

Conceptualization, S.M.; investigation, writing, and analysis, Y.Z. and M.Z.; revision, R.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant numbers 51778088 and 52073036).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available upon request due to restrictions, e.g., privacy or ethical. The data presented in this study are available upon request from the corresponding author. The data are not publicly available due to subsequent analyses and publications.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Sun, G.; Huo, X.; Wang, H.; Hazell, P.J.; Li, Q. On the structural parameters of honeycomb-core sandwich panels against low-velocity impact. Compos. Part B-Eng. 2021, 216, 108881. [Google Scholar] [CrossRef]
  2. Ma, Q.; Rejab, M.R.M.; Siregar, J.P.; Guan, Z. A review of the recent trends on core structures and impact response of sandwich panels. J. Compos. Mater. 2021, 55, 2513–2555. [Google Scholar] [CrossRef]
  3. Qi, C.; Pei, L.Z.; Remennikov, A.; Yang, S.; Liu, J.; Wang, J.S.; Liao, X.W. Parametric study and optimization of the protect system containing a re-entrant hexagon cored sandwich panel under blast impact. Compos. Struct. 2020, 252, 112711. [Google Scholar] [CrossRef]
  4. Chen, G.; Zhang, P.; Deng, N.; Cai, S.; Cheng, Y.; Liu, J. Paper tube-guided blast response of sandwich panels with auxetic re-entrant and regular hexagonal honeycomb cores—An experimental study. Eng. Struct. 2022, 253, 113790. [Google Scholar] [CrossRef]
  5. Usta, F.; Türkmen, H.S.; Scarpa, F. High-velocity impact resistance of doubly curved sandwich panels with re-entrant honeycomb and foam core. Int. J. Impact Eng. 2022, 165, 104230. [Google Scholar] [CrossRef]
  6. Ma, N.; Han, Q.; Han, S.; Li, C. Hierarchical re-entrant honeycomb metamaterial for energy absorption and vibration insulation. Int. J. Mech. Sci. 2023, 250, 108307. [Google Scholar] [CrossRef]
  7. Razgordanisharahi, A.; Ghassabi, A.A.; Hellmich, C. Free vibration analysis of cylindrical honeycomb sandwich panels using state-space Levy method. Thin-Walled Struct. 2023, 182, 110308. [Google Scholar] [CrossRef]
  8. Sun, M.; Wowk, D.; Mechefske, C.; Alexander, E.; Kim, I.Y. Surface and honeycomb core damage in adhesively bonded aluminum sandwich panels subjected to low-velocity impact. Compos. Part B-Eng. 2022, 230, 109506. [Google Scholar] [CrossRef]
  9. Zhang, X.; Xu, F.; Zang, Y.; Feng, W. Experimental and numerical investigation on damage behavior of honeycomb sandwich panel subjected to low-velocity impact. Compos. Struct. 2020, 236, 111882. [Google Scholar] [CrossRef]
  10. Li, X.; Zhang, P.; Wang, Z.; Wu, G.; Zhao, L. Dynamic behavior of aluminum honeycomb sandwich panels under air blast, Experiment and numerical analysis. Compos. Struct. 2014, 108, 1001–1008. [Google Scholar] [CrossRef]
  11. Lu, Z.X.; Li, X.; Yang, Z.Y.; Xie, F. Novel structure with negative Poisson’s ratio and enhanced Young’s modulus. Compos. Struct. 2016, 138, 243–252. [Google Scholar] [CrossRef]
  12. Chen, Z.; Wu, X.; Xie, Y.M.; Wang, Z.; Zhou, S. Re-entrant auxetic lattices with enhanced stiffness, A numerical study. Int. J. Mech. Sci. 2020, 178, 105619. [Google Scholar] [CrossRef]
  13. Fu, M.H.; Chen, Y.; Hu, L.L. A novel auxetic honeycomb with enhanced in-plane stiffness and buckling strength. Compos. Struct. 2017, 160, 574–585. [Google Scholar] [CrossRef]
  14. Zhang, X.Y.; Ren, X.; Zhang, Y.; Xie, Y.M. A novel auxetic metamaterial with enhanced mechanical properties and tunable auxeticity. Thin-Walled Struct. 2022, 174, 109162. [Google Scholar] [CrossRef]
  15. Yu, R.; Luo, W.; Yuan, H.; Liu, J.; He, W.; Yu, Z. Experimental and numerical research on foam filled re-entrant cellular structure with negative Poisson’s ratio. Thin-Walled Struct. 2020, 153, 106679. [Google Scholar] [CrossRef]
  16. Luo, H.C.; Ren, X.; Zhang, Y.; Zhang, X.Y.; Zhang, X.G.; Luo, C.; Cheng, X.; Xie, Y.M. Mechanical properties of foam-filled hexagonal and re-entrant honeycombs under uniaxial compression. Compos. Struct. 2022, 280, 114922. [Google Scholar] [CrossRef]
  17. Qi, C.; Jiang, F.; Remennikov, A.; Pei, L.Z.; Liu, J.; Wang, J.S.; Liao, X.W.; Yang, S. Quasi-static crushing behavior of novel re-entrant circular auxetic honeycombs. Compos. Part B-Eng. 2020, 197, 108117. [Google Scholar] [CrossRef]
  18. Zhang, X.C.; An, C.C.; Shen, Z.F.; Wu, H.X.; Yang, W.G.; Bai, J.P. Dynamic crushing responses of bio-inspired re-entrant auxetic honeycombs under in-plane impact loading. Mater. Today Commun. 2020, 23, 100918. [Google Scholar] [CrossRef]
  19. Saravanos, D.A.; Pereira, J.M. Effects of interply damping layers on the dynamic characteristics ofcomposite plates. AIAA J. 1992, 30, 2906–2913. [Google Scholar] [CrossRef]
  20. Wang, H.J.; Chen, L.W. Vibration and damping analysis of a three-layered composite annular plate with a viscoelastic mid-layer. Compos. Struct. 2002, 58, 563–570. [Google Scholar] [CrossRef]
  21. Yim, J.H.; Cho, S.Y.; Seo, Y.J.; Jang, B.Z. A study on material damping of 0 laminated composite sandwich cantilever beams with a viscoelastic layer. Compos. Struct. 2003, 60, 367–374. [Google Scholar] [CrossRef]
  22. Plagianakos, T.S.; Saravanos, D.A. High-order layerwise mechanics and finite element for the damped dynamic characteristics of sandwich composite beams. Int. J. Solids Struct. 2004, 41, 6853–6871. [Google Scholar] [CrossRef]
  23. Yang, J.S.; Xiong, J.; Ma, L.; Wang, B.; Zhang, G.; Wu, L. Vibration and damping characteristics of hybrid carbon fiber composite pyramidal truss sandwich panels with viscoelastic layers. Compos. Struct. 2013, 106, 570–580. [Google Scholar] [CrossRef]
  24. Berthelot, J.M.; Sefrani, Y. Damping analysis of unidirectional glass fiber composites with interleaved viscoelastic layers, experimental investigation and discussion. J. Compos. Mater. 2006, 40, 1911–1932. [Google Scholar] [CrossRef]
  25. Zhong, Y.F.; Qin, W.; Yu, W. Variational asymptotic homogenization of magnetoelectro-elastic materials with coated fibers. Compos. Struc. 2015, 133, 300–311. [Google Scholar] [CrossRef]
  26. Peng, X.; Zhong, Y.F.; Shi, Z. Global buckling analysis of composite honeycomb sandwich plate with negative Poisson’s ratio (CSP-RHC) using variational asymptotic equivalent model. Compos. Struct. 2021, 264, 113721. [Google Scholar]
  27. Chen, J.; Zhong, Y.; Luo, Q.; Shi, Z. Static and dynamic analysis of Isogrid Stiffened Composite Plates (ISCP) using equivalent model based on variational asymptotic method. Thin-Walled Struct. 2021, 163, 107671. [Google Scholar] [CrossRef]
  28. Cesnik, C.E.; Hodges, D.H. VABS, A new concept for composite rotor blade cross-sectional modeling. J. Am. Helicopter Soc. 1997, 42, 27–38. [Google Scholar] [CrossRef]
  29. Peng, X.; Zhong, Y.F.; Shi, Z. Free flexural vibration analysis of composite sandwich plate with reentrant honeycomb cores using homogenized plate model. J. Sound Vib. 2022, 529, 116955. [Google Scholar] [CrossRef]
  30. Grima, J.N.; Oliveri, L.; Attard, D.; Ellul, B.; Gatt, R.; Cicala, G.; Recca, G. Hexagonal honeycombs with zero Poisson’s ratios and enhanced stiffness. Adv. Eng. Mater. 2010, 12, 855–862. [Google Scholar] [CrossRef]
  31. Li, H.; Dong, B.; Cao, J.; Zhao, J.; Xiong, J.; Yang, Y.; Du, D.; Sun, W.; Wang, X.; Wu, H. Vibration behaviours of foam-filled grille composite sandwich cylindrical shells. Int. J. Mech. Sci. 2023, 256, 108533. [Google Scholar] [CrossRef]
  32. Wattanasakulpong, N.; Eiadtrong, S. Transient responses of sandwich plates with a functionally graded porous core, Jacobi–Ritz method. Int. J. Struct. Stab. Dyn. 2023, 23, 2350039. [Google Scholar] [CrossRef]
Figure 1. Recent studies focusing on enhancing the stiffness of re-entrant honeycomb structures (the yellow ribs serve as additional stiffeners): (a) introducing strengthening ribs [11,12], (b) rhombic [13] and wedge-shaped reinforcement [14], (c) filled foam reinforcement [15,16], and (d) modifying the shape of the inclined ribs [17,18].
Figure 1. Recent studies focusing on enhancing the stiffness of re-entrant honeycomb structures (the yellow ribs serve as additional stiffeners): (a) introducing strengthening ribs [11,12], (b) rhombic [13] and wedge-shaped reinforcement [14], (c) filled foam reinforcement [15,16], and (d) modifying the shape of the inclined ribs [17,18].
Materials 17 04067 g001
Figure 2. Regularre-entrant honeycomb and its variant bowtie-shaped honeycomb achieved by adding vertical walls highlighted in red boxes in (b).
Figure 2. Regularre-entrant honeycomb and its variant bowtie-shaped honeycomb achieved by adding vertical walls highlighted in red boxes in (b).
Materials 17 04067 g002
Figure 3. Schematic diagram of equivalent analysis of bowtie-shaped honeycomb core composite sandwich panel (BHC-CSP).
Figure 3. Schematic diagram of equivalent analysis of bowtie-shaped honeycomb core composite sandwich panel (BHC-CSP).
Materials 17 04067 g003
Figure 4. Core cell decomposition for strain energy integral.
Figure 4. Core cell decomposition for strain energy integral.
Materials 17 04067 g004
Figure 5. Schematic diagram of different boundary conditions for free-vibration analysis.
Figure 5. Schematic diagram of different boundary conditions for free-vibration analysis.
Materials 17 04067 g005
Figure 6. Composition of BHC-CSP.
Figure 6. Composition of BHC-CSP.
Materials 17 04067 g006
Figure 7. Laminated structure of viscous layers and composite materials.
Figure 7. Laminated structure of viscous layers and composite materials.
Materials 17 04067 g007
Figure 8. Effect of varying thicknesses of PMI damping layers on the natural frequency of BHC-CSP.
Figure 8. Effect of varying thicknesses of PMI damping layers on the natural frequency of BHC-CSP.
Materials 17 04067 g008
Figure 9. Schematic diagram of excitation location on the sandwich panel under SSSS BCs.
Figure 9. Schematic diagram of excitation location on the sandwich panel under SSSS BCs.
Materials 17 04067 g009
Figure 10. Effect of viscous damping layers on frequency-domain response of damped BHC-CSP.
Figure 10. Effect of viscous damping layers on frequency-domain response of damped BHC-CSP.
Materials 17 04067 g010
Figure 11. Effect of viscous damping layers on the time-domain response of damped BHC-CSP.
Figure 11. Effect of viscous damping layers on the time-domain response of damped BHC-CSP.
Materials 17 04067 g011
Figure 12. Influence of height ratio of core to facesheet on the vibration characteristics of damped BHC-CSP.
Figure 12. Influence of height ratio of core to facesheet on the vibration characteristics of damped BHC-CSP.
Materials 17 04067 g012
Figure 13. Influence of wall thickness-to-height ratio on the vibration characteristics of damped BHC-CSP.
Figure 13. Influence of wall thickness-to-height ratio on the vibration characteristics of damped BHC-CSP.
Materials 17 04067 g013
Figure 14. Influence of re-entrant angle on the vibration characteristics of damped BHC-CSP.
Figure 14. Influence of re-entrant angle on the vibration characteristics of damped BHC-CSP.
Materials 17 04067 g014
Figure 15. Influence of aspect ratio on the vibration characteristics of damped BHC-CSP.
Figure 15. Influence of aspect ratio on the vibration characteristics of damped BHC-CSP.
Materials 17 04067 g015
Figure 16. Different honeycomb sandwich panels.
Figure 16. Different honeycomb sandwich panels.
Materials 17 04067 g016
Figure 17. Comparison of natural frequencies for different honeycomb sandwich panels.
Figure 17. Comparison of natural frequencies for different honeycomb sandwich panels.
Materials 17 04067 g017
Figure 18. Comparison of time-domain response for different honeycomb sandwich panels.
Figure 18. Comparison of time-domain response for different honeycomb sandwich panels.
Materials 17 04067 g018
Table 1. Material properties.
Table 1. Material properties.
Material ParameterCarbon FiberEpoxy ResinAluminiumPMIT-SW308
Elasticity modulus E 1 / GPa 105.54.570263430
E 2 = E 3 / GPa 11.34.5
Shear modulus G 12 = G 13 / GPa 3.231.727113200
G 23 / GPa 3.181.7
Poisson’s ratio ν 12 = ν 13 0.280.380.3--
ν 23 0.530.38
Density ρ / g · cm 3 1.491.62.70.2050.184
Table 2. Comparison of dynamic analysis efficiency between 2D-ERM and 3D-FEM.
Table 2. Comparison of dynamic analysis efficiency between 2D-ERM and 3D-FEM.
ModelsElement
Type
Element
Number
Node
Number
Calculation Time (sec)
Free
Vibration
Frequency-Domain
Analysis
Time-Domain
Analysis
Unit cell
homogenization
C3D1086,636140,257Homogenization time, 14
2D-ERMS4R5400555126220423
3D-FEMC3D10340,653600,3011362848325,920
Table 3. Variation range of geometric parameters.
Table 3. Variation range of geometric parameters.
Geometrical ParameterLayup Mode
h c / h f t c / h c α l 4 / l 2
2.5∼150.2 45 1.5[45/−45/0/90]s
100.05∼0.30 45 1.5[45/−45/0/90]s
100.2 10 60 1.5[45/−45/0/90]s
100.2 45 1.1∼1.6[45/−45/0/90]s
Table 4. Influence of geometric parameters on vibration characteristics of damped BHC-CSP.
Table 4. Influence of geometric parameters on vibration characteristics of damped BHC-CSP.
ParametersRanges | Δ ρ |
g / cm 3
Δ A 22 Δ ρ
KN · cm 3 / g
Δ D 22 Δ ρ
KN · cm 4 / g
Δ ω Δ ρ
Hz · cm 3 / g
Δ R A Δ ρ
cm 4 / g
h c / h f 2.5∼154.59051.48375.06490.24923.603
t c / h c 0.05∼0.304.286120.88510.17450.2460.446
α 10 60 2.509261.15814.70547.4690.522
l 4 / l 2 1.1∼1.60.87413.3658.42331.7320.499
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

Miao, S.; Zhong, Y.; Zhang, M.; Liu, R. Forced-Vibration Characteristics of Bowtie-Shaped Honeycomb Composite Sandwich Panel with Viscoelastic Damping Layer. Materials 2024, 17, 4067. https://doi.org/10.3390/ma17164067

AMA Style

Miao S, Zhong Y, Zhang M, Liu R. Forced-Vibration Characteristics of Bowtie-Shaped Honeycomb Composite Sandwich Panel with Viscoelastic Damping Layer. Materials. 2024; 17(16):4067. https://doi.org/10.3390/ma17164067

Chicago/Turabian Style

Miao, Siqi, Yifeng Zhong, Mingtao Zhang, and Rong Liu. 2024. "Forced-Vibration Characteristics of Bowtie-Shaped Honeycomb Composite Sandwich Panel with Viscoelastic Damping Layer" Materials 17, no. 16: 4067. https://doi.org/10.3390/ma17164067

APA Style

Miao, S., Zhong, Y., Zhang, M., & Liu, R. (2024). Forced-Vibration Characteristics of Bowtie-Shaped Honeycomb Composite Sandwich Panel with Viscoelastic Damping Layer. Materials, 17(16), 4067. https://doi.org/10.3390/ma17164067

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