Next Article in Journal / Special Issue
Comparison of Tensile and Fatigue Properties of Copper Thin Film Depending on Process Method
Previous Article in Journal
A Novel Intelligent Method for Bearing Fault Diagnosis Based on EEMD Permutation Entropy and GG Clustering
Previous Article in Special Issue
Bulletproof Performance of Composite Plate Fabricated Using Shear Thickening Fluid and Natural Fiber Paper
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Active Vibration Suppression of Stiffened Composite Panels with Piezoelectric Materials under Blast Loads

1
Department of Aerospace Engineering, Chosun University, Gwangju 61452, Korea
2
National Research Foundation of Korea, Daejeon 34113, Korea
3
Department of Naval Architecture and Ocean Engineering, Chosun University, Gwangju 61452, Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(1), 387; https://doi.org/10.3390/app10010387
Submission received: 30 November 2019 / Revised: 19 December 2019 / Accepted: 26 December 2019 / Published: 4 January 2020
(This article belongs to the Special Issue Selected Papers from the ICMR 2019)

Abstract

:
Transient responses of stiffened panels with piezoelectric sensors and actuators are studied under normal blast loads. The air vehicles could be exposed to blast pulses generated by an explosion or shock-wave disturbances. Thus, active vibration suppression of the vehicles is important under blast loadings. The structural model is designed as a laminated composite panel with lead zirconate titanate (PZT) piezoceramic layers embedded on both top and bottom surfaces. A uniformly distributed blast load is assumed over the whole of the panel surface. The first-order shear deformation theory of plate is adopted, and the extended Hamilton’s principle is applied to derive the equations of motions. The numerical model is verified by the comparison with previous data. Using linear quadratic regulator (LQR) control algorithm, vibration characteristics and dynamic responses are compared. As piezoelectric patches are attached on the whole of the surface, the effect of the stiffener’s location is studied. Furthermore, the influences of the patch’s positions are also investigated through subjection to the blast wave. From various results, in order to get the best control performances, the research aims to find the optimum position of sensor and actuator pairs that is most effective under blast load environments.

1. Introduction

Recently, the research relevant to the piezoelectric materials to control the vibration of structures is actively reported in the literature. The applications of the piezoelectric materials are vibration suppression, shape control, active damping, size control, and energy harvesting, etc. The materials can be used as sensors and actuators. At first, direct piezoelectric effect means that the materials generate an electric charge while the model is subjected to a mechanical deformation. Conversely, they could be deformed when some electric charges are applied to structures. This is called the converse piezoelectric effect.
Many engineers have studied composite structures with piezoelectric sensors and actuators. Reddy [1] investigated theoretical formulations and finite element models using the shear deformation plate theories. They analyzed the laminated plates with sensors and actuators under mechanical and electrical loads. Liu et al. [2] presented the active vibration suppression of the laminated composite plates. Lam and Ng [3] also simulated theoretical formulation for the plates with piezoelectric materials. Balamurugan and Narayanan [4] showed the mechanics analysis and performances of piezolaminated plate and shell structures. In addition, Tolliver et al. [5] described the finite element analysis of the piezoelectric transducer. The model of multilayer piezoelectric actuators is derived based on the physical analysis by Zhang et al. [6]. The proposed methods are easy to handle and apply the piezo models. Furthermore, Qian [7] presented the design, optimization, and testing of piezoelectric harvesters for energy scavenging from human walking. Also, Xie et al. [8] designed a self-powered and wireless radio frequency (RF) transmission device using by combination of triboelectric nanogenerator (TENG) and off-the-shelf piezoelectric element.
On the other hand, the dynamic responses of structures under air blast loads have been studied for many years. Luccioni et al. [9] showed the evaluation of pressures and impulses produced by blast loads with the aid of hydrocodes. The numerical and experimental research of a stiffened laminated composite panel under a shock wave was represented by Turkmen and Mecitoglu [10]. Furthermore, Jacinto et al. [11] analyzed numerical responses of metallic panels subjected to explosive loadings.
Recently, research and development of stiffened structures is one of the important subjects in aircraft and aerospace industries. Globally, many engineers have tried to increase the stiffness by using composite stiffened structures in the space launch vehicle. Nowadays, Ma et al. [12,13] studied active vibration control and sound radiation from a rib stiffened panel using velocity feedback method. Also, Liu et al. [14] presented the dynamic analytical solution of piezoelectric stack using piezo-elasticity theory. In addition, Chung et al. [15] studied the dynamic behaviors of stiffened composite model with piezoelectric patches under airflow. The blast loads applied a sharp pressure to the composite structure in a moment.
Though numerous studies on piezoelectric stiffened models have been performed widely up until now, research to suppress the vibrations of stiffened models caused by aerodynamic flows such as blast loads has been studied in a limited range. In other words, many studies have been focused on damages to the structures. Therefore, active vibration suppression of piezoelectric stiffened laminated panels subjected to normal blast waves is studied in this paper. Due to the characteristics of piezoelectric materials as fast responders, the vibration behaviors of stiffened composite panels could be controlled easily by piezoelectric actuators under a blast loading. A stiffened panel with piezoelectric effect is formulated using finite element method. Numerical formulations and results are verified with previous data. Then, vibration behaviors and transient responses are controlled using linear quadratic regulator (LQR) algorithms. Furthermore, numerical results for vibration suppression of stiffened panels under blast load are discussed in detail. The effect of the stiffeners location of the cantilevered model is considered like as general wings of aircrafts. Furthermore, the influences of the piezoelectric patches position also investigated. Consequently, these studies could be easily applied to various types of structures.

2. Stiffened Laminated Panel with Piezoelectric Materials

In aircraft vehicles, there are many skins with an array of stiffening ribs in the wings and fuselage. These structures have been widely used because of their light weight and ability to withstand dynamic loads. In most cases, the catastrophic disaster of the air vehicles was due to the failure of the fuselage under blast loads. The overpressure of the blast loading generated typically acts during few milliseconds. Thus, active vibration control of stiffened composite structures is an important problem. In this regard, a piezo-laminated composite model is designed with stiffeners (shown in Figure 1).
The stiffeners are placed parallel to the geometric coordinates x and y. In this figure, the thickness of the model, the thickness of stiffener, the width of the x-stiffener, and the width of the y-stiffener are t p , t y s , b x s , and b y s , respectively. The model is a laminated composite structure with lead zirconate titanate (PZT) piezoceramic layers embedded on both top and bottom surfaces to act as a sensor and an actuator, respectively. In addition, the sensor and actuator layers are designed as a single layer.

2.1. Finite Element Formulations

First of all, the k t h layer’s piezoelectric coupling relations between the elastic and the electrical fields could be presented, and the equations can be written as follows [15]:
D k = e k ε k + ε k E k
σ = Q ε e k T E k
where ε , σ , D , and E are the strain, stress component, electric displacement, and electrical field vectors, respectively. In addition, e , ε , and Q are piezoelectric constants, permittivity coefficients, and reduced elastic constants matrices, respectively. Equation (1) means the direct piezoelectric effect, and Equation (2) indicates the converse piezoelectric effect.
The stiffened plate elements are derived with the plate element and stiffener elements. The first-order shear deformation theory (FSDT) applied, the displacement fields are defined as follows:
u ( x , y , z ) = u p o ( x , y ) + z θ x p ( x , y ) , v ( x , y , z ) = v p o ( x , y ) + z θ y p ( x , y ) , w ( x , y , z ) = w p o ( x , y , z ) .
where u p o , v p o , and w p o are the mid-plane displacements, and θ x p and θ y p are the rotations of transverse normal about x and y , respectively. The displacement fields of the x -directional stiffener and y -directional stiffener are expressed as follows:
u x s ( x , z ) = u x s o ( x ) + z θ x s ( x ) , w x s ( x , z ) = w x s o ( x ) v y s ( y , z ) = v y s o ( y ) + z θ y s ( y ) , w y s ( y , z ) = w y s o ( y )
where x -axis is taken along the x -directional stiffener center line and z -axis is its upward normal line. The detailed formulations can be found in reference [16].
The strains are presented as a function of the nodal displacement variables.
ε p = { ε x x ε y y γ x y } p = { u p , x o v p , y o u p , y o + v p , x o } + z { θ x p , x θ y p , y θ x p , y + θ y p , x } = ε p o + z κ p = B m δ p + z B b δ p γ p = { γ x z γ y z } p = { w p , x o θ p , x w p , y o θ p , y } = B s δ p
where B m , B b , and B s are the derivative operators between the strain and nodal displacements.
The stress resultant can be written as
F = A i j ε p o + B i j κ p t p e T E d z M = B i j ε p o + D i j κ p t p z e T E d z V p = S γ p
where ( A i j , B i j , D i j ) = k = 1 n z k 1 z k Q i j ( 1 , z , z 2 ) d z and S i j = k = 1 n κ p z k 1 z k Q i j d z with ( i , j = 4 , 5 ) . κ p = 5 / 6 is shear correction factor. A i j , B i j , D i j , and S i j are extensional stiffness, bending-extension coupling stiffness, bending stiffness, and shear stiffness matrix, respectively.
Then, the electric field vector considering the piezoelectric actuator layers can be derived as follows:
E = { E x E y E z } = { 0 0 1 / t a } ϕ a = B ϕ ϕ a
where t a and ϕ a are the thickness of the actuator and the electric potential, respectively. The present element has elastic degrees of freedom u p o , v p o , θ x p , θ y p , and w p o per node and electrical degree of freedom ϕ per piezoelectric layer.
On the other hand, the strains of the x -stiffener are indicated as follows [17]:
ε x s = { ε x x γ x z } x s = [ 1 z 0 0 0 0 1 1 ] { u x s , x o θ x s , x θ x s w x s , x o } = H x s ε x s o
The stress resultant can be derived as follows:
{ F x s M x s V x s } = [ A x s , 11 B x s , 11 0 0 B x s , 11 D x s , 11 0 0 0 0 S x s , 55 S x s , 55 ] ε x s o = D x s ε x s o
where A x s , i j , B x s , i j , D x s , i j , S x s , i j = k = 1 n z k 1 z k Q i j ( 1 , z , z 2 , κ x s ) d z .
To obtain the energy equations, the strain energy function is derived using the variational principles as follows:
V e = V p e + V x s e = 1 2 δ p e T A ( B m T A i j B m + B m T B i j B b + B b T B i j B m + B b T D i j B b + B s T S i j B s ) d A δ p e + 1 2 δ p e T A ( B m T E 1 B ϕ + B b T E 2 B ϕ ) d A ϕ e + 1 2 δ p e T L T x s T B x s T D x s B x s T x s d x δ p e = 1 2 δ p e T K p e δ p e + 1 2 δ p e T K u ϕ e ϕ e + 1 2 δ p e T K x s e δ p e
where E 1 , E 2 = z p e i j T ( 1 , z ) d z , D x s = b x s t x s H x s T Q x s H x s d z . K p e is the element stiffness matrix of plate and K u ϕ e is the element elastic-electric coupling stiffness matrix. K x s e is the element stiffness matrix of the x -stiffeners.
The element kinetic energy can be written as
T e = T p e + T x s e = 1 2 δ ˙ p e T A [ P ( N ¯ p u T N ¯ p u + N ¯ p v T N ¯ p v + N ¯ p w T N ¯ p w ) + I ( N ¯ p θ x T N ¯ p θ x + N ¯ p θ y T N ¯ p θ y ) ] d A δ ˙ p e + 1 2 δ ˙ p e T T x s T L b [ P ( N ¯ x s u T N ¯ x s u + N ¯ x s w T N ¯ x s w ) + I N ¯ x s θ x T N ¯ x s θ x ] d x T x s δ ˙ p e = 1 2 δ ˙ p e T M p e δ ˙ p e + 1 2 δ ˙ p e T M x s e δ ˙ p e
where P , I = k = 1 n z k 1 z k ρ ( 1 , z 2 ) d z and N ¯ is the shape function matrices relating the primary variables.
M p e is the element mass matrix of plate and M x s e is the element mass matrix of the stiffeners.
The element electrical energy is calculated as follows:
W E e = 1 2 V p E T ( e ε p + ε E ) d V = 1 2 ϕ e T A ( B ϕ T E 1 B m + B ϕ T E 2 B b ) d A δ p e + 1 2 ϕ e T A B ϕ T E 3 B ϕ d A ϕ e = 1 2 ϕ e T K ϕ u δ p e + 1 2 ϕ e T K ϕ ϕ ϕ e
where K ϕ u is element elastic-electric coupling stiffness matrix and K ϕ ϕ is element electric stiffness matrix. E 3 = z p ε i j d z and V p is the volume of the piezoelectric layer.
The virtual work and the applied electrical charge density are derived as follows:
Δ W S e = Δ δ p e S 1 B f T f s d A Δ ϕ e T S 2 B ϕ T q s d A = Δ δ p e F s e + Δ ϕ e T F q e
where S 1 and S 2 are the surface areas that the surface forces and electrical charges are applied, respectively. F s e is the applied mechanical force due to surface forces, and F q e is the applied electrical charge.
The governing equation can be described using the Hamilton’s principle.
M u u e δ ¨ p e + K u u e δ p e + K u ϕ e ϕ e = F s e
K ϕ u e δ p e + K ϕ ϕ e ϕ e = F q e
where M u u e = M p e + M x s e , K u u e = K p e + K x s e .
Then, the global dynamic equation is obtained as follows:
M u u δ ¨ p + C u u δ ˙ p + ( K u u + K u ϕ K ϕ ϕ 1 K ϕ u ) δ p = F s + K u ϕ ϕ a
where damping matrix C u u is assumed as C u u = α M u u + β K u u .
Next, modal transformations and state-space formulations are derived. The nodal displacement is represented by
δ p = Φ η
The global dynamic equation can be transformed to the reduced modal space form as follows:
Φ T M u u Φ η + Φ T C u u Φ η ˙ + Φ T ( K u u + K u ϕ K ϕ ϕ 1 K ϕ u ) Φ η = Φ T F s + Φ T K u ϕ ϕ a
M ¯ η ¨ + C ¯ η ˙ + K ¯ η = F ¯ s + K ¯ u ϕ ϕ a
Introducing the state-space variable ξ = { η η ˙ } T , state-space form of the global dynamic equation is given by
ξ ˙ = A s t ξ + B s t ϕ a + U f
where ϕ a is the control input and some valuables are defined as follows:
A s t = [ 0 I M ¯ 1 K ¯ M ¯ 1 C ¯ ] , B s t = { 0 M ¯ 1 K ¯ u ϕ } , U f = { 0 M ¯ 1 F ¯ s }
The output equations can be written in the physical coordinates and state space coordinates as
y = C o δ p , y ¯ = { C o Φ 0 } ξ = C ¯ o ξ
where C o is the output matrix.
Finally, the linear quadratic regulator (LQR) control algorithm [4,18] is applied. The control gain can be obtained.
J = 1 2 0 ( ξ T Q ξ + ϕ a T R ϕ a ) d t
where R and Q are positive definite weighting matrices and semi-positive definite, respectively. The actuating voltage can be derived as follows:
ϕ a ( t ) = G c ξ = R 1 B s t T P ξ
Which is determined by the constant solution P of the Riccati equation.
A s t T P + P A s t P B s t R 1 B s t T P + Q = 0

2.2. Blast Loads

The air vehicles could be exposed to blast pulses generated by an explosion or shock-wave disturbances. Figure 2 presents the expected form of an ideal blast wave. It is increased at the shock front by an abrupt pressure. A negative phase follows, and oscillations quickly die away. These oscillations are not dominant compared with the first positive phase. The pressure is also assumed to be distributed over the whole of the surface. The total pressure is defined by Friedländer decay function [19] as follows:
P ( t ) b l a s t = p m ( 1 t t p ) e a t / t p
where p m , t p , and a are peak pressure, positive phase duration, and the waveform parameter, respectively. The parameters of the Friedländer decay equation are chosen as p m = 28.9 kPa, t p = 1.8 ms, and a = 0.35 [10].

3. Numerical Results and Discussions

3.1. Code Verifications

To verify the accuracy of the present numerical results, three cases are compared with the previous data in the literature. The first case is the cantilevered composite model with the upper and lower surfaces by piezoelectric ceramics. The plates have four composite layers [ 45 / 45 / 45 / 45 ] . The total thickness of the composite model is 1 mm and each layer has the same thickness of 0.25 mm. The two outer piezoelectric layers are of 0.1 mm thickness. In addition, the adhesive layers are neglected. The panel consists of T300/976 graphite-epoxy composites, and the piezoelectric material is PZT F1195N, as presented in Table 1 [21]. The centerline deflections of the composite model with various input voltages are presented in Figure 3. The results are in good agreement with the data obtained by Lam et al. [21].
Secondly, the deflection responses of the model under blast loads. The material properties of E , ν , and ρ are 206.84 GPa, 0.3, and 7900 kg/m3, respectively. The dynamic response is calculated using mode superposition method, and time step is 0.05 ms. Figure 4 shows the deflection history at the center of the model. The present results are almost same with the previous works by Sheikh et al. [22].
The last case is the polyester model with stiffeners placed at the center. Two opposite edges are free, and the others are clamped. The structure and stiffeners are made with the isotropic materials as follows: E = 7.0 GPa, G = 2.6 GPa, ρ = 1504.2 kg/m3 and ν = 0.345. The natural frequencies of the model are well-matched with the previous data [16,23], as shown in Table 2.

3.2. The Evaluation of Stiffeners’ Locations

Figure 5 shows the cantilevered laminated composite panel with two x-stiffeners. Piezoelectric ceramics are boded at the upper and lower surfaces of the structure. The material properties of the composite panel (T300/976 graphite-epoxy) and the piezoceramic (PZT G1195N) are already presented in Table 1. The structural model is consist of four layers ( 15 / 15 / 15 / 15 ) with ( 0.4 / 0.4 / 0.4 / 0.4 ) mm layers. The stiffeners with width 6 mm and depth 12 mm are made of the same materials. The thickness of each piezoceramic layer is 0.2 mm.
As mentioned in advance, the upper piezo layer acts as a sensor and the lower layer is assumed as an actuator. Jacinto et al. [11] investigated the experimental and computational analysis of a panel under blast loads, and the data are used in this research. The LQR optimal control algorithm is adopted to the structure. Modal superposition method is applied to the analysis. In addition, the initial damping ratio is assumed to be 0.7%. The transient responses are calculated with the Newmark-direct integration method.
At first, the effect of the stiffener’s location on the suppression of the model is studied. The two kinds of cases are presented in Figure 5b,c. The stiffeners are placed at both end sides in Case 1, while the stiffeners are located in the middle part in Case 2. Then, the first six natural frequencies of the stiffened panel with different stiffener location listed in Table 3. It can be seen from Table 3 that the natural frequencies of the Stiffener Case 1 are higher than that of the Stiffener Case 2. This means that the stiffened panel in Case 1 is stiffer than Case 2.
Figure 6 shows magnitude-frequency responses for the two cases of stiffened panel using LQR method. The weighting matrices Q and R are as follows:
Q = 1 7 × 10 35 [ Ω 2 0 0 I ] , R = 1 × 10 2   and   1 × 10 3
where Ω 2 is a diagonal matrix that stores the squares of the natural frequencies ω i 2 and I is the unit matrix with suitable size. When the location of the stiffeners is centered (Case 2), the magnitude of the model is the higher than in the other case (Case 1). In addition, the peak point appeared earlier than Case 1.
Figure 7 describes phase-frequency responses for the two cases of stiffened panel using LQR method. The results represent that controlled responses have the lower magnitude of the first, second, and third bending modes. When the stiffened panel is under uniformly distributed blast pressures, the vertical displacement ( W 2 ) at the tip (0.4, 0.1 m), the transverse bending and the lateral twisting for various R are given in Figure 8, Figure 9 and Figure 10. The transverse bending ( W T ) and the lateral twisting ( W R ) are defined as in Ref. [4], and the equations are as follows:
W T = W 2 ( W 1 + W 3 ) / 2 C ,   W R = W 1 W 3 C
where W 1 , W 2 , and W 3 are the vertical deflections at locations already shown in Figure 5a.
Figure 8 represents vertical displacement responses with two locations of stiffener using LQR control method. The vertical displacement is initially high in Case 1, however the displacement quickly decreased according to time. This is due to the increased stiffness with stiffeners location. On the other hand, both cases showed effective responses to the LQR controller. In other words, it can be seen that the piezoelectric actuator works effectively.
Figure 9 shows transverse bending responses with two locations of stiffener using the LQR control method. The results of the Stiffener Cases 1 and 2 show a similar tendency of vertical displacement in Figure 8. The transverse bending is also initially high in Case 1, however the bending quickly decreased according to time. The vibration oscillations are damped out more quickly with higher control parameters.
Figure 10 shows lateral twisting responses with two locations of stiffeners using the LQR control method. Unlike the previous results in Figure 8 and Figure 9, the lateral twisting was of higher initial value in Case 2. Overall, it could be mentioned that the bending modes can be controlled effectively by the piezoelectric actuators, but the control of torsion is weak for two stiffener cases.

3.3. Evaluation of Positions for Piezoelectric Patches

The effect of positions is studied for the piezoelectric patches on the transient behaviors of the stiffened model under blast load. To investigate the effect of sensors and actuators locations on the control performances, the stiffened model is performed for two different positions of piezoelectric pairs that are attached on top and bottom surfaces, as shown in Figure 11. Each model has eight piezoelectric patches, and the dimensions of each sensor and actuator patch are 0.05 m × 0.05 m. The center points of the Models I and II are 0.05 m and 0.35 m away from the clamped edge, respectively. To investigate the control effects, two parameters are chosen. One parameter is settling time that oscillation does not exceed the limit value of each transient responses. The limit values of the vertical displacement, the transverse bending, and the lateral twisting are 0.4 mm, 0.002, and 0.002, respectively. For other parameter, the limit is the maximum response of the models.
Table 4 shows that the settling time of the two models is obtained using the LQR control algorithm under uniformly distributed blast loads. The models I and II with Stiffener Case 1 have longer settling time than that of the others in the vertical displacement and the transverse bending. The vertical displacement and the transverse bending of the stiffener’s location at Stiffener Case 2 is damped out more quickly for the Model I and II. Also, lateral twisting closed to the free end is damped out more quickly.
Table 5 presents the comparison of maximum responses. These results indicate that Model II has the lower peak response for many cases of stiffeners’ locations. From Table 4 and Table 5, different control performance results are obtained in various positions of piezoelectric pairs. It is important to obtain the optimal design with maximum control effect.

4. Conclusions

This paper presents vibration suppression control of a stiffened model with piezoelectric sensors and actuators under blast loads. The air vehicles could be exposed to blast pulses generated by an explosion or shock-wave disturbances. The blast wave is assumed to be uniformly distributed on the panel surface. The linear quadratic regulator (LQR) method is chosen as a control algorithm. As piezoelectric patches are distributed on the whole surfaces of the model, the effect of stiffener’s location is investigated. Furthermore, the effect of piezoelectric patches position is also studied on the transient responses of the stiffened structure under blast load. The most effective position of patches can be determined with settling time and maximum responses. Consequently, this presented analysis could be easily applied to various types of models. In order to get the best control performances, the research aimed to find the optimum position of sensor and actuator pairs that is most effective under blast load environments.

Author Contributions

C.-Y.L. participated in the design, computational analysis and writing. J.-Y.J. played an important role in initial overall research. S.-M.J. participated in the supervision of analysis and paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Reddy, J.N. On laminated composite plates with integrated sensors and actuators. Eng. Struct. 1999, 21, 569–593. [Google Scholar] [CrossRef]
  2. Liu, G.R.; Peng, X.Q.; Lam, K.Y.; Tani, J. Vibration control simulation of laminated composite plates with integrated piezoelectrics. J. Sound Vibr. 1999, 220, 827–846. [Google Scholar] [CrossRef]
  3. Lam, K.Y.; Ng, T.Y. Active control of composite plates with integrated piezoelectric sensors and actuators under various dynamic loading conditions. Smart Mater. Struct. 1999, 8, 223–237. [Google Scholar] [CrossRef]
  4. Balamurugan, V.; Narayanan, S. Shell finite element for smart piezoelectric composite plate/shell structures and its application to the study of active vibration control. Finite Elem. Anal. Des. 2001, 37, 713–738. [Google Scholar] [CrossRef]
  5. Tolliver, L.; Xu, T.B.; Jiang, X. Finite element analysis of the piezoelectric stacked-hybats transducer. Smart Mater. Struct. 2013, 22, 035015. [Google Scholar] [CrossRef]
  6. Zhang, Y.K.; Lu, T.F.; Al-Sarawi, S. Formulation of a simple distributed-parameter model of multilayer piezoelectric actuators. J. Intell. Mater. Syst. Struct. 2016, 27, 1485–1491. [Google Scholar] [CrossRef]
  7. Qian, F.; Xu, T.B.; Zuo, L. Design, optimization, modeling and testing of a piezoelectric footwear energy harvester. Energy Convers. Manag. 2018, 171, 1352–1364. [Google Scholar] [CrossRef]
  8. Xie, Y.; Long, J.; Zhao, P.; Chen, J.; Luo, J.; Zhang, Z.; Li, K.; Han, Y.; Hao, X.; Qu, Z.; et al. A self-powered radio frequency (RF) transmission system based on the combination of triboelectric nanogenerator (TENG) and piezoelectric element for disaster rescue/relief. Nano Energy 2018, 54, 331–340. [Google Scholar] [CrossRef]
  9. Luccioni, B.; Ambrosini, R.D.; Danesi, R. Blast load assessment using hydrocodes. Eng. Struct. 2006, 28, 1736–1744. [Google Scholar] [CrossRef] [Green Version]
  10. Turkmen, H.S.; Mecitoglu, Z. Dynamic response of a stiffened laminated composite plate subjected to blast load. J. Sound Vibr. 1999, 221, 371–389. [Google Scholar] [CrossRef]
  11. Jacinto, A.C.; Ambrosini, R.K.; Danesi, R.F. Experimental and computational analysis of plates under air blast loading. Int. J. Impact Eng. 2001, 25, 927–947. [Google Scholar] [CrossRef]
  12. Ma, X.; Wang, L.; Xu, J. Active Vibration control of rib stiffened plate by using decentralized velocity feedback controllers with inertial actuators. Appl. Sci. 2019, 9, 3188. [Google Scholar] [CrossRef] [Green Version]
  13. Ma, X.; Chen, K.; Xu, J. Active control of sound transmission through orthogonally rib stiffened double-panel structure: Mechanism analysis. Appl. Sci. 2019, 9, 3286. [Google Scholar] [CrossRef] [Green Version]
  14. Liu, X.; Wang, J.; Li, W. Dynamic analytical solution of a piezoelectric stack utilized in an actuator and a generator. Appl. Sci. 2018, 8, 1779. [Google Scholar] [CrossRef] [Green Version]
  15. Chung, N.T.; Thuy, N.N.; Thu, D.T.N.; Chau, L.H. Numerical and experimental analysis of the dynamic behavior of piezoelectric stiffened composite plates subjected to airflow. Math. Probl. Eng. 2019, 2019, 10. [Google Scholar] [CrossRef]
  16. Attaf, B.; Hollaway, L. Vibration analyses of stiffened and unstiffened composite plates subjected to in-plane loads. Composites 1990, 21, 117–126. [Google Scholar] [CrossRef]
  17. Sadek, E.A.; Tawfik, S.A. A finite element model for the analysis of stiffened laminated plates. Comput. Struct. 2000, 75, 369–383. [Google Scholar] [CrossRef]
  18. Kumar, E.V.; Raaja, G.S.; Jerome, J. Adaptive PSO for optimal LQR tracking control of 2 DoF laboratory helicopter. Appl. Soft. Comput. 2016, 41, 77–90. [Google Scholar] [CrossRef]
  19. Gupta, A.D.; Gregory, F.H.; Bitting, R.L.; Bhattacharya, S. Dynamic analysis of an explosively loaded hinged rectangular plate. Comput. Struct. 1987, 26, 339–344. [Google Scholar] [CrossRef]
  20. Hetherington, J.G.; Smith, P.D. Blast and Ballistic Loading of Structures; Butterworth-Heinemann Ltd.: London, UK, 1994; pp. 30–57. [Google Scholar]
  21. Lam, K.Y.; Peng, X.Q.; Liu, G.R.; Reddy, J.N. A finite-element model for piezoelectric composite laminates. Smart Mater. Struct. 1997, 6, 583–591. [Google Scholar] [CrossRef]
  22. Sheikh, A.H.; Mukhopadhyay, M. Linear and nonlinear transient vibration analysis of stiffened plate structures. Finite Elem. Anal. Des. 2002, 38, 477–502. [Google Scholar] [CrossRef]
  23. Mukhergee, A.; Joshi, S.P.; Ganguli, A. Active vibration control of piezolaminated stiffened plates. Compos. Struct. 2002, 55, 435–443. [Google Scholar] [CrossRef]
Figure 1. Stiffened model with piezoelectric sensors and actuators.
Figure 1. Stiffened model with piezoelectric sensors and actuators.
Applsci 10 00387 g001
Figure 2. Variations of blast loading according to time [20].
Figure 2. Variations of blast loading according to time [20].
Applsci 10 00387 g002
Figure 3. The centerline deflections with various input voltages.
Figure 3. The centerline deflections with various input voltages.
Applsci 10 00387 g003
Figure 4. Displacement-time responses of the center of the square model subjected to air blast.
Figure 4. Displacement-time responses of the center of the square model subjected to air blast.
Applsci 10 00387 g004
Figure 5. Cantilevered piezoelectric composite model and two cases of the stiffener’s location: (a) Stiffened composite model; (b) Stiffener Case 1; (c) Stiffener Case 2.
Figure 5. Cantilevered piezoelectric composite model and two cases of the stiffener’s location: (a) Stiffened composite model; (b) Stiffener Case 1; (c) Stiffener Case 2.
Applsci 10 00387 g005
Figure 6. Magnitude-frequency responses with two locations of stiffeners using linear quadratic regulator (LQR) control: (a) Stiffener Case 1; (b) Stiffener Case 2.
Figure 6. Magnitude-frequency responses with two locations of stiffeners using linear quadratic regulator (LQR) control: (a) Stiffener Case 1; (b) Stiffener Case 2.
Applsci 10 00387 g006
Figure 7. Phase-frequency responses with two locations of stiffener using LQR control: (a) Stiffener Case 1; (b) Stiffener Case 2.
Figure 7. Phase-frequency responses with two locations of stiffener using LQR control: (a) Stiffener Case 1; (b) Stiffener Case 2.
Applsci 10 00387 g007
Figure 8. Vertical displacement responses with two locations of stiffener using LQR control: (a) Stiffener Case 1; (b) Stiffener Case 2.
Figure 8. Vertical displacement responses with two locations of stiffener using LQR control: (a) Stiffener Case 1; (b) Stiffener Case 2.
Applsci 10 00387 g008
Figure 9. Transverse bending responses with two locations of stiffener using LQR control: (a) Stiffener Case 1; (b) Stiffener Case 2.
Figure 9. Transverse bending responses with two locations of stiffener using LQR control: (a) Stiffener Case 1; (b) Stiffener Case 2.
Applsci 10 00387 g009
Figure 10. Lateral twisting responses with two locations of stiffener using LQR control: (a) Stiffener Case 1; (b) Stiffener Case 2.
Figure 10. Lateral twisting responses with two locations of stiffener using LQR control: (a) Stiffener Case 1; (b) Stiffener Case 2.
Applsci 10 00387 g010
Figure 11. Two models of piezoelectric patch position.
Figure 11. Two models of piezoelectric patch position.
Applsci 10 00387 g011
Table 1. Material properties [21].
Table 1. Material properties [21].
PropertiesPZT PiezoceramicT300/976
Young’s modulus (GPa): E 11 63.0150
E 22 = E 33 63.09.0
Poisson’s ratio: ν 12 = ν 13 0.30.3
ν 23 0.30.3
Shear modulus (GPa): G 12 = G 13 24.27.10
G 23 24.22.50
Density (kg/m3): ρ 76001600
Piezoelectric constants (m/V): d 31 = d 32 254 × 10−12-
Electrical permittivity (F/m): 11 = 22 15.3 × 10−9-
33 15.0 × 10−9-
Table 2. Natural frequencies of stiffened models.
Table 2. Natural frequencies of stiffened models.
Natural Frequency (Hz)Ref. [16]Ref. [23]Present
1st68.6168.4768.13
2nd71.2068.6668.77
3rd124.70119.59117.55
4th150.40162.16160.128
5th183.20177.11177.45
6th184.80177.39178.31
Table 3. Natural frequencies of different stiffener cases.
Table 3. Natural frequencies of different stiffener cases.
Natural Frequency (Hz)Stiffener Case 1Stiffener Case 2
1st67.609839.6451
2nd117.727378.6680
3rd155.3218133.6110
4th212.1017178.2795
5th290.4140214.3758
6th379.4205273.2984
Table 4. Settling time of two models for LQR control.
Table 4. Settling time of two models for LQR control.
Stiffener LocationModel IModel II
Vertical displacement (s)case 10.28480.2204
case 20.11390.1844
Transverse bending (s)case 10.28570.1883
case 20.25690.1910
Lateral twisting (s)case 10.20270.0592
case 21.04720.3362
Table 5. Maximum responses of two models for LQR control.
Table 5. Maximum responses of two models for LQR control.
Stiffener LocationModel IModel II
Maximum vertical displacement (mm)case 18.06965.6853
case 21.83782.6561
Maximum transverse bending (mm)case 10.03590.0183
case 20.04530.0127
Maximum lateral twisting (mm)case 10.00700.0035
case 20.01950.0056

Share and Cite

MDPI and ACS Style

Lee, C.-Y.; Jung, J.-Y.; Jeong, S.-M. Active Vibration Suppression of Stiffened Composite Panels with Piezoelectric Materials under Blast Loads. Appl. Sci. 2020, 10, 387. https://doi.org/10.3390/app10010387

AMA Style

Lee C-Y, Jung J-Y, Jeong S-M. Active Vibration Suppression of Stiffened Composite Panels with Piezoelectric Materials under Blast Loads. Applied Sciences. 2020; 10(1):387. https://doi.org/10.3390/app10010387

Chicago/Turabian Style

Lee, Chang-Yull, Jin-Young Jung, and Se-Min Jeong. 2020. "Active Vibration Suppression of Stiffened Composite Panels with Piezoelectric Materials under Blast Loads" Applied Sciences 10, no. 1: 387. https://doi.org/10.3390/app10010387

APA Style

Lee, C. -Y., Jung, J. -Y., & Jeong, S. -M. (2020). Active Vibration Suppression of Stiffened Composite Panels with Piezoelectric Materials under Blast Loads. Applied Sciences, 10(1), 387. https://doi.org/10.3390/app10010387

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