Abstract:The buckling and parametric resonance characteristics of laminated composite spherical sandwich shell panels with viscoelastic material (VEM) core are investigated in the present analysis considering full geometric nonlinearity in the Green–Lagrange sense. The study includes the longitudinal strain and normal strain in the transverse direction along with transverse shear deformation of the VEM core. The core displacements are considered to be varying linearly along the thickness and those of the face sheets follow first-order shear deformation theory. An eight-noded sandwich shell finite element of the serendipity family is adopted to discretize the sandwich shell panel domain. The finite element-based equation of motion is derived using Hamilton's principle in the form of the Mathieu–Hill-type equation. The dynamic instability regions are obtained by applying Hsu's criteria-based Saito–Otomi conditions to the transformed equation motion. An in-house finite element-based code is developed in the MATLAB platform to solve the stability problem and to establish the stability regions. A parametric study is carried out to investigate the influence of different system parameters on the critical buckling load and the parametric resonance of the sandwich shell panels. It is noted that an increase in core and constraining layer thicknesses increases the critical buckling load of the sandwich shell panels. The stability boundaries are observed to shift toward a higher-excitation-frequency region in the stability diagram with an increase in constraining layer thickness and a decrease in aspect ratio.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is about the static and dynamic stabilities of spherical sandwich shell plates with a viscoelastic material core layer and laminated composite panels under uniaxial and biaxial harmonic excitations. Specifically, the buckling and parametric resonance characteristics of these structures under geometrically nonlinear (in the Green - Lagrange sense) conditions are studied. By considering the longitudinal strain, normal strain and transverse shear deformation of the core layer, the paper uses the finite element method (FEM) and Hamilton's principle to establish the equations of motion, and uses Hsu's criteria and Saito - Otomi conditions to determine the regions of dynamic instability.
### Main research contents:
1. **Static stability analysis**:
- The influence of different system parameters (such as core layer thickness, constraint layer thickness, aspect ratio and boundary conditions) on the critical buckling load is studied.
- The equations of motion for static stability are established by the finite element method (FEM) and Hamilton's principle.
2. **Dynamic stability analysis**:
- The dynamic stability of spherical sandwich shell plates under uniaxial and biaxial harmonic excitations is studied.
- The boundaries of the dynamic instability regions are determined by Hsu's criteria and Saito - Otomi conditions.
### Research methods:
- **Finite element method (FEM)**: The sandwich shell plate domain is discretized by using eight - node sandwich shell elements (serendipity family).
- **Hamilton's principle**: It is used to derive the finite - element equations of motion based on the Mathieu - Hill equation.
- **Hsu's criteria and Saito - Otomi conditions**: They are used to determine the boundaries of the dynamic instability regions.
### Main conclusions:
- Increasing the thicknesses of the core layer and the constraint layer can improve the critical buckling load of the sandwich shell plate.
- As the thickness of the constraint layer increases and the aspect ratio decreases, the frequency range of the stability boundary moves towards the region of higher excitation frequencies.
- For symmetric and antisymmetric laminated composite panels, the buckling parameters show different trends with the change of the number of layers.
### Formula representation:
- **Geometrically nonlinear strain**:
\[
\epsilon_{xx}^{NL} = \frac{1}{2} \left[ \left( U_{i,x} + \frac{W_i}{R_x} \right)^2 + \left( V_{i,x} \right)^2 + \left( W_{i,x} - \frac{U_i}{R_x} \right)^2 \right]
\]
\[
\epsilon_{yy}^{NL} = \frac{1}{2} \left[ \left( U_{i,y} \right)^2 + \left( V_{i,y} + \frac{W_i}{R_y} \right)^2 + \left( W_{i,y} - \frac{V_i}{R_y} \right)^2 \right]
\]
\[
\epsilon_{xy}^{NL} = \left( U_{i,y} \right) \left( U_{i,x} + \frac{W_i}{R_x} \right) + \left( V_{i,y} + \frac{W_i}{R_y} \right) \left( V_{i,x} \right) + \left( W_{i,y} - \frac{V_i}{R_y} \right) \left( W_{i,x} - \frac{U_i}{R_x} \right)
\]
- **Dynamics equation**:
\[
M \ddot{u} + K u - N_0(t) K_G u = 0
\]
where \( M \) is the mass matrix, \( K \) is the elastic stiffness matrix, and \( K_G \) is the geometric stiffness.