Next Article in Journal
An Optical Test Strip for the Detection of Benzoic Acid in Food
Previous Article in Journal
Visual Odometry Based on Structural Matching of Local Invariant Features Using Stereo Camera Sensor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electromechanical Impedance Response of a Cracked Timoshenko Beam

203 Office, Xi’an Research Institute of High-Technology, Xi’an 710025, China
*
Author to whom correspondence should be addressed.
Sensors 2011, 11(7), 7285-7301; https://doi.org/10.3390/s110707285
Submission received: 1 June 2011 / Revised: 16 July 2011 / Accepted: 17 July 2011 / Published: 22 July 2011
(This article belongs to the Section Physical Sensors)

Abstract

: Typically, the Electromechanical Impedance (EMI) technique does not use an analytical model for basic damage identification. However, an accurate model is necessary for getting more information about any damage. In this paper, an EMI model is presented for predicting the electromechanical impedance of a cracked beam structure quantitatively. A coupled system of a cracked Timoshenko beam with a pair of PZT patches bonded on the top and bottom surfaces has been considered, where the bonding layers are assumed as a Kelvin-Voigt material. The shear lag model is introduced to describe the load transfer between the PZT patches and the beam structure. The beam crack is simulated as a massless torsional spring; the dynamic equations of the coupled system are derived, which include the crack information and the inertial forces of both PZT patches and adhesive layers. According to the boundary conditions and continuity conditions, the analytical expression of the admittance of PZT patch is obtained. In the case study, the influences of crack and the inertial forces of PZT patches are analyzed. The results show that: (1) the inertial forces affects significantly in high frequency band; and (2) the use of appropriate frequency range can improve the accuracy of damage identification.

1. Introduction

In recent years, the Electromechanical Impedance (EMI) technique has emerged as a promising structural health monitoring (SHM) method. It has been successfully applied to various engineering disciplines, including aerospace and aircraft structures [16] and civil structures [79]. In this technique, a piezoelectric ceramic lead-zirconate-titanate (PZT) patch is surface bonded to the monitored structure or embedded into a new composite construction and excited by an alternating voltage sweeping signal. Any physical change in the structures will result in a change of the structural mechanical impedance. Due to the electromechanical coupling between the PZT transducers, the EMI signature extracted from PZT transducer is directly related to the mechanical impedance of the monitored structure. Consequently, any structural physical change will induce changes in the EMI signature of the PZT transducer. Hence, for SHM applications, PZT EMI spectra can be compared with a baseline measurement during the service period of the monitored structure. Any change in the spectra is an indication of a change of the structural integrity, which may be caused by the presence of damage.

Although the EMI-based SHM does not typically use an analytical model for basic damage identification, it is necessary to establish a model for more advanced features of SHM, such as damage evaluation and prognosis [10]. Many researchers have investigated the EMI model. As early as the 1990s, Liang et al. [11] proposed the first one-dimensional (1D) EMI-type model for a PZT-structure interaction system. In his model the electromechanical admittance (inverse of impedance) of PZT can be expressed analytically by structure impedance. Zhou et al. [12] extended the 1D impedance method to model a two-dimensional (2D) PZT element coupled to a 2D monitored structure. On the basis of the concept of “effective impedance”, Bhalla and Soh [13,14] improved Zhou’s model. Yang presented a generic model for predicting the electromechanical impedance of one-dimensional and two-dimensional PZT—structure interaction systems [15]. Based on the concept “sum impedance”, Annamdas and Soh [16,17] proposed a three-dimensional (3D) EMI model. In general, all the models above ignored the adhesive layer between the PZT transducer and the monitored structure. However, many experimental and theoretical analysis results [1822] have demonstrated that the mechanical interaction between the PZT transducer and the monitored structure occurs through interfacial shear stress. Hence, the adhesive layer has to be considered. Ong et al. presented an EMI model which considered the shear lag effect of the bond layer [23]. Suresh Bhalla et al. [24] incorporated the shear lag effect into the existing 1D and 2D EMI and obtained an improved model for them. Yan et al. [25,26] considered the adhesive layer, and presented an EMI model for Timoshenko beams and Mindlin-Herrmann rods. Although the factor of the adhesive layer was considered, the inertial terms of PZT and adhesive layer produced by the motion along with the monitored structure have not been taken into account in the above models. Pietrzakowski [22] noticed the inertial terms of PZT patches and studied the influence of bonding layer on the beam response, but his method is not suitable for high frequency EMI techniques due to the use of Euler beams.

For some simple pristine structure, the analytical EMI model can be obtained. However, when damages are induced in structures resulting in the possibly inhomogeneity of material properties, it is difficult to derive an analytical formulation. In order to quantitatively identify structural damages, numerical methods and approximate approaches have been adopted. Naidu and Soh [27] and Tseng and Wang [28] obtained the relationships between the EMI signatures and the structural changes by the Finite Element Method (FEM). However, for the purpose of predicting the response accurately in a high frequency range, a very small size element is needed. Therefore the FEM method is time-consuming. Yang [29] and Xu [30] applied the p-Ritz method to establish EMI models for health monitoring of beams, plates and cylindrical shell structures with various boundaries. Their models could calculate the EMI response effectively and accurately below 50 kHz, however, the accuracy of prediction was reduced with the increased frequency.

In order to quantitatively identify the structural damage in a more efficient way, a more accurate EMI model is needed. This study proposes an electromechanical impedance model for health monitoring of beam structures. Different to the existing EMI model, the proposed model not only contains physical parameters of damage and the shear lag effect of the bonding layer, but also includes the transversal inertial forces of PZT patches and adhesive induced by the transversal motion of beam structures.

In this paper, a pair of PZT patches bonded symmetrically onto the top and bottom surfaces of a rectangular beam with a crack was activated out of phase to create a pure bending excitation. Pure extensions in PZT patches were assumed. A shear lag model was applied to describe the behavior of bonding layer which is assumed to be a Kelvin-Voigt material. The PZT-adhesive layer-cracked beam coupled structural system was considered. The coupled system is divided into four sections due to the crack cross-section and the location of PZT patches. The crack is simulated by a massless torsional spring. Taking into account the inertial forces of PZT patches and bonding layer, the boundary problem is formulated by the dynamic equations of the coupled system. According to the boundary conditions and continuity conditions, the solution of the coupled system can be obtained, and then the analytical relationship between PZT admittance (inverse of impedance) and the damage parameters such as location and depth is derived. Finally, numerical results are presented and discussed to validate the proposed theoretical analysis.

2. Formulation of the Coupled Structure and Analysis

As shown in Figure 1, two identical PZT patches are bonded symmetrically onto the top and bottom surfaces of the Timoshenko beam. Two out-of-phase alternating electrical fields are applied to the PZT patches.

The beam was activated to create a pure bending vibration. A crack with a depth of hc is located at x1. For a better analysis, the beam is divided into four parts due to the crack cross-section and the location of PZT patch (Figure 1). The dynamic behavior of each part is governed by different equations. According to boundary conditions and continuity conditions, the dynamic equations of PZT-structure interaction system can be solved. The analysis method used here is similar to the method used by Pietrzakowski [22], but different from that work, the Timoshenko beam theory and the inertial of bonding layers are taken into account in this paper.

2.1. Dynamic Equations

The governing equations for the section with PZT patch (x2 < x < x3) are derived considering an infinitesimal element shown in Figure 2. The longitudinal motion of the PZT patches is defined by introducing inertial forces. The governing equations can be expressed as:

E p 2 u p x 2 ρ p 2 u p t 2 τ h p = 0 E p 2 u p x 2 ρ p 2 u p t 2 τ h p = 0
where Ep, ρp, hp, denote Young’s modulus, density and thickness of the two PZT patches, respectively; up and u p are longitudinal displacement of the upper PZT and the lower PZT, respectively; The shear stresses τ and τ′ transmitted by the piezoelements are determined by the following stress-strain relation:
τ = G a * ( u p u s ) h a τ = G a * ( u p u s ) h a
where us and u s are the upper and lower beam surface longitudinal displacement, respectively; ha is the bonding layer thickness; hb is the beam thickness; and w is the transverse displacement of beam; G a * is the linear function of differential operator which for the Kelvin-Voigt model of bonding layer material can be expressed as:
G a * = G a ( 1 + C b t )
where Ga is Kirchhoff’s modulus; Cb denotes the retardation time.

Because the two PZT patches are identical and are actived by two out-of-phase alternating electrical fields with same amplitude, so, the following relationship can be obtained:

τ = τ u p = u p u s = u s = h b 2 w x

Considering the dynamic coupling between the piezoelement and the beam, and including the mass of piezoelement and adhesive in the inertial force of intensity qb (shown in Figure 2), the transverse motion of Timoshenko beam is described. The governing equations of beam are expressed as follows:

M x + τ h b b b Q + γ ρ b I 2 ϕ t 2 = 0 ,     Q x γ ρ b A b 2 w t 2 = 0 ,     M = E b I ϕ x ,     Q = κ A b G b ( w x ϕ )
where Eb, Gb are the Young’s modulus and shear modulus, respectively; κ = π2/12 is the shear correction factor; I is the inertia moment; M, Q are the bending moment, transverse force of beam, respectively; ϕ is the rotation angle of beam due to pure bending; bb, ρb, Ab denote the width, density and cross-section area of beam, respectively; the mass ratio γ is determined by following relation:
γ = ρ b A b + 2 ρ p A p + 2 ρ a A a ρ b   A b A b = h b b b , A a = h a b a , A p = h p b p
where ρa is the density of adhesive layer.

From Equations (15), the following equations can be obtained:

E p 2 u p x 2 G a * h p h a ( u p + h b 2 w 2 ) = ρ p 2 u p t 2
E b I 4 w x 4 + h b b b G a * h a ( u p x + h p 2 2 w x 2 ) I γ ρ b ( E b κ G b + 1 ) 4 w x 2 t 2 γ ρ b A b 2 w t 2 I γ 2 ρ b 2 κ G b 4 w t 4 = 0

Assuming the high frequency harmonic voltage loaded on the upper and lower PZT patch are V(x,t) = V̄(x)eiωt and V′(x,t) = V̄(x)ei(ωt+π): then the steady state response of coupled structural system can be expressed as:

u p ( x , t ) = u ¯ p ( x ) e iωt , u s ( x , t ) = u ¯ s ( x ) e iωt , τ ( x , t ) = τ ¯ ( x ) e iωt , w ( x , t ) = w ¯ ( x ) e iωt

From Equation (8) and Equation (9), we can obtain:

ɛ ¯ p = u ¯ p x = E b I h a h b b b G ¯ a 4 w ¯ x 4 h b 2 2 w ¯ x 2 + I h a γ ρ b ω 2 h b b b G ¯ a ( E b κ G b + 1 ) 2 w ¯ x 2 + ( I h a γ 2 ρ b 2 ω 4 h b b b G ¯ a κ G b h a γ ρ b A b ω 2 h b b b G ¯ a ) w ¯
where, a = Ga (1 + iCbω).

Substituting Equation (10) into Equation (7), the following equation can be obtained:

d 6 w ¯ d x 6 + k 1 d 4 w ¯ d x 4 + k 2 d 2 w ¯ d x 2 + k 3 w ¯ = 0
where:
k 1 = ( h b 2 b b G ¯ a 2 E b I h a + G ¯ a h p h a E p ρ p ω 2 E p γ ρ b ω 2 κ G b γ ρ b ω 2 E b ) k 2 = ω 4 ρ p γ ρ b E p E b + ω 4 γ 2 ρ b 2 κ G b E b G ¯ a γ ρ b ω 2 E p h a E b h p G ¯ a γ ρ b ω 2 E p h a h p κ G b ω 2 γ ρ b A b I E b + ω 4 γ ρ b ρ b κ G b E p h b 2 b b G ¯ a ρ p ω 2 2 E p h a E b I k 3 = ( G ¯ a γ 2 ρ b 2 ω 4 h p h a E p E b κ G b ρ p γ 2 ρ b 2 ω 6 E p E b κ G b G ¯ a γ ρ b A b ω 2 h p h a E p E b I + γ ρ b A b ω 4 E p E b I )

The solution of Equation (11) can be expressed as:

w ¯ ( x ) = c 1 e β 1 x + c 2 e β 1 x + c 3 e β 2 x + c 4 e β 2 x + c 5 e β 3 x + c 6 e β 3 x , x [ x 2 + , x 3 ]
where, ci (i = 1,⋯,6) are the undetermined constant coefficients, ± βi (i = 1,2,3) are the characteristic roots of the following equation:
β 6 + k 1 β 4 + k 2 β 2 + k 3 = 0

Substituting the Equation (12), Equation (9) into Equation (10), we can obtain:

ɛ ¯ p = g 1 c 1 e β 1 x + g 1 c 2 e β 1 x + g 2 c 3 e β 2 x + g 2 c 4 e β 2 x + g 3 c 5 e β 3 x + g 3 c 6 e β 3 x
where
g i = E b I h a h b b b G ¯ a β i 4 h b 2 β i 2 + I h a γ ρ b ω 2 h b b b G ¯ a ( E b κ G b + 1 ) β i 2 + ( I h a γ 2 ρ b 2 ω 4 h b b b G ¯ a κ G b h a γ ρ b A b ω 2 h b b b G ¯ a ) , ( i = 1 , 2 , 3 )

Substituting the Equation (12), Equation (9) into Equation (7), we can obtain:

u ¯ p = g 1 β 1 c 1 e β 1 x + g 1 β 1 c 2 e β 1 x + g 2 β 2 c 3 e β 2 x + g 2 β 2 c 4 e β 2 x + g 3 β 3 c 5 e β 3 x + g 3 β 3 c 6 e β 3 x

By substituting the Equations (2), (11), (15) into Equation (5), respectively, the bending moment M, shear stress Q and rotational angle ϕ can be obtained. In order to solve the unknown constants ci (i = 1,....,6),the boundary conditions and the continuity conditions are required, so, the analytical expression of other beam section is also needed. The basic equations of a uniform beam section without bonded PZT patches are:

M x Q + ρ b I 2 ϕ t 2 = 0 , Q x ρ b A b 2 w t 2 = 0 , M = E b I ϕ x , Q = κ A b G b ( w x ϕ )

From Equation (16), the following relations can be obtained:

E b I d 4 w ¯ d x 4 + ( E b κ G b + 1 ) I ρ b ω 2 d 2 w ¯ d x 2 ( ρ b A b ω 2 ρ b 2 I ω 4 κ G b ) w ¯ = 0

The steady state solution of the Equation (17) can be expressed as:

w ¯ ( x ) = c 7 e 1 2 k 4 x + c 8 e 1 2 k 4 x + c 9 e 1 2 k 5 x + c 10 e 1 2 k 5 x , x [ 0 , x 1 ] w ¯ ( x ) = c 11 e 1 2 k 4 x + c 12 e 1 2 k 4 x + c 13 e 1 2 k 5 x + c 14 e 1 2 k 5 x , x [ x 1 + , x 2 ] w ¯ ( x ) = c 15 e 1 2 k 4 x + c 16 e 1 2 k 4 x + c 17 e 1 2 k 5 x + c 18 e 1 2 k 5 x , x [ x 3 + , l ]
where, ci (i = 7,…,18) are the undetermined constant coefficients; and k4, k5 can be expressed as:
k 4 = ( 2 k 6 2 k 6 2 4 k 7 ) 1 / 2 , k 5 = ( 2 k 6 + 2 k 6 2 4 k 7 ) 1 / 2 k 6 = ( 1 κ G b + 1 E b ) ρ b ω 2 , k 7 = ( ρ b A b ω 2 E b I ρ b 2 ω 4 E b κ G b )

Substituting the Equation (18) into Equation (16), we can obtain the bending moment M, shear stress Q and rotational angle ϕ of the beam without bonded PZT patches. In this study, the crack is assumed to be a fully open crack, and the depth of crack is relative small. Therefore, the crack on the beam can be modeled by a weightless rotational spring [31] as shown in Figure 3. The softness of the spring θ is a function of the beam thickness and the depth of crack [31,32] and can be expressed as:

θ = 6 π ( 1 μ 2 ) h b ( 0.6272 ξ 2 0.4533 ξ 3 + 4.5948 ξ 4 9.9736 ξ 5 + 20.2948 ξ 6 33.0351 ξ 7 + 47.1063 ξ 8 40.7556 ξ 9 + 19.6 ξ 10 )
where ξ = hc/hb; hc is the depth of crack, and μ is the Poisson’s ratio of the beam.

Considering the continuity of transverse placement, bending moment, shear stress and discontinuity of slope at the crack, we can obtain the following equation:

w + = w , M + = M , Q + = Q , ϕ + = ϕ + θ ϕ x
where, “+” denotes the right side of crack; “−” denotes the left side of crack.

2.2. Steady-State Solution

There are 18 undetermined constant coefficients ci (i = 1,⋯,18) in the steady analytical solution of coupled system. The undetermined coefficients can be determined through boundary condition and continuity condition. For the classical end of beam, one has the following equations:

hinged end

w = 0 , M = 0

clamped end

w = 0 , ϕ = 0

free end

M = 0 , Q = 0

For the cantilever beam, the boundary conditions are expressed as:

w ¯ ( 0 ) = ϕ ¯ ( 0 ) = 0
M ¯ ( 0 ) = Q ¯ ( 0 ) = 0

The boundary conditions of PZT patch can be expressed as:

ɛ ¯ p ( x 2 ) = d 31 E ¯ 3 = d 31 V ¯ h p , ɛ ¯ p ( x 3 ) = d 31 E ¯ 3 = d 31 V ¯ h p

According to the continuity of beam deflection, slope, bending moment and shear stress at the borders of the sections at x = x1, x = x2, x = x3 we can obtain the following equations:

ϕ ¯ ( x 1 ) = ϕ ¯ ( x 1 + ) + θ ( d ϕ ¯ d x | x = x 1 ) ,    ϕ ¯ ( x 2 ) = ϕ ¯ ( x 2 + ) ,    ϕ ¯ ( x 3 ) = ϕ ¯ ( x 3 + ) ,
M ¯ ( x 1 ) = M ¯ ( x 1 + ) ,    M ¯ ( x 2 ) = M ¯ ( x 2 + ) ,    M ¯ ( x 3 ) = M ¯ ( x 3 + ) ,
Q ¯ ( x 1 ) = Q ¯ ( x 1 + ) ,    Q ¯ ( x 2 ) = Q ¯ ( x 2 + ) ,    Q ¯ ( x 3 ) = Q ¯ ( x 3 + )
w ¯ ( x 1 ) = w ¯ ( x 1 + ) ,    w ¯ ( x 2 ) = w ¯ ( x 2 + ) ,    w ¯ ( x 3 ) = w ¯ ( x 3 + )

The 18 unknown coefficients ci = (i = 1,⋯,18) are obtained from the system of algebraic equations determined by the boundary and continuity conditions [Equations (2430)].

3. Electromechanical Signatures

Consider a pure extension of PZT patches, and then the corresponding constitutive equations of PZT patches can be expressed as [33]:

ɛ 1 = σ 1 / E ¯ p + d 31 E 3 D 3 = d 31 σ 1 + ɛ ¯ 33 T E 3
where ε1, σ1 are the strain and stress along x direction, respectively; D3, E3 are the electric flux density and electric field intensity along height direction, respectively; Ēp = Ep (1 + ) is the complex Young’s modulus of the PZT material at zero electric field with η denoting the mechanical loss factor; ɛ ¯ 33 T = ɛ 33 T ( 1 j δ ) is the complex dielectric constant at zero stress with δ denoting the dielectric loss factor of PZT patch.

From Equation (31), we can obtain the following equation:

D 3 = d 31 ( ɛ 1 d 31 E 3 ) E ¯ p + ɛ ¯ 33 T E 3 = d 31 E ¯ p ɛ 1 + ( ɛ ¯ 33 T d 31 2 E ¯ p ) E 3

The electric current passing through the upper PZT patch can be determined from the electric displacement as:

I p = i ω 0 b a x 2 x 3 D 3 d x d y = i ω b a d 31 E ¯ p [ u p ( x 3 ) u p ( x 2 ) ] + i ω E 3 b a ( x 3 x 2 ) ( ɛ ¯ 33 T d 31 2 E ¯ p )
where Ip is the electric current passing through PZT patch.

Electric admittance of the upper PZT patch can be determined as:

Y = I p V = i ω b a d 31 E ¯ p ( u ¯ p ( x 3 ) u ¯ p ( x 2 ) ) V ¯ + i ω ( ɛ ¯ 33 T d 31 2 E ¯ p ) b a ( x 3 x 2 ) h p = i ω b a d 31 E ¯ p ( u ¯ p ( x 3 ) u ¯ p ( x 2 ) ) V ¯ + i ω ( ɛ ¯ 33 T d 31 2 E ¯ p ) b a l p h p = i ω b a d 31 E ¯ p V ¯ i = 1 3 g i β i [ ( e β i x 3 e β i x 2 ) c 2 × i ( e β i x 3 e β i x 2 ) c i ] + i ω ( ɛ ¯ 33 T d 31 2 E ¯ p ) b a l p h p
where Y is the electric admittance of PZT patch.

For a piezoelectric system without damage on the PZT patches, the parameters of the PZT patches can be regarded as constants. Hence, the change of admittance of PZT patch is only determined by the first term on the right side of Equation (34). From Equation (34), we can see that any changes in beam structure will lead to a change of the admittance signature of the PZT patch.

4. Numerical Examples

In this paper, a cantilever beam with a pair of PZT patches bonded symmetrically onto its top and bottom surfaces is studied. The PZT patches are located at x2 = 200 mm from the left side of beam (Figure 1). Geometric parameters and material constants of the beam, PZT patches and adhesive layer are listed in Tables 13, respectively.

By considering the inertia terms of PZT patches and bonding layers caused by their motion with the beam, the values of γ against the beam thickness hb can be calculated and the result is listed in Table 4. If the inertia terms are not taken into account, then γ = 1. In order to study the effect of inertia, the differences of the admittances of PZT under two scenarios (considering/without considering the inertia effect) are compared. To verify the accuracy and the reliability of the proposed analytical model, the results are compared to the simulation results obtained by FEA (finite element analysis), which can predict well the experimental results [34]. Yang et al. developed a multi-physics simulation method of EMI modeling, which can effect direct acquisition of PZT electrical admittance [35]. Yang’s method was also used in this work. The test specimen is numerically modeled in the ANSYS 10.0 workspace as illustrated in Figure 4. The constitutive data, in accordance to PZT-5A, are assigned to the PZT patch as given in Table 5. The properties of beam and bonding layer are listed in Tables 1 and 3, respectively. The PZT patches are modeled with Solid 5 elements and the bonding layer and beam with Solid 45 elements. The sizes of the elements are less than 1.0 mm. An alternating (sinusoidal) voltage of 1 V was applied across the PZT patch for excitation. The detailed process can be found in the references [34] and [35].

The results shown in Figures 57 suggest that the proposed analytical model provides reasonable predictions of the the FEA results, as the major resonance peaks are well predicted. For any beam thickness hb, the corresponding frequencies to the peaks of admittance curves considering the inertia terms are smaller than those without considering the inertia terms; and the difference between the curves is increased with the increasing frequency. When hb = 10 mm (Figure 7), the difference between curves of γ = 1.58 and γ = 1.0 is small, however, it is increased with the increasing frequency, and the difference of peak frequency is about 0.45 kHz near 96 kHz. From Figure 7, It can be seen that the inertia effect cannot be ignored if the frequency is high (>80 kHz). When hb = 5 mm (Figure 6), the difference between curves of γ = 2.116 and γ = 1.0 is more obvious than that of hb = 10 mm, and the difference of peak frequency is about 1.5 kHz near 98 kHz. From Figure 6, it can be seen that the influence of inertia can’t be ignored when the frequency is higher than 40 kHz. When hb = 2.5 mm (Figure 5), the curves of γ = 3.321 have a distinct difference with the curve of γ = 1.0. The difference of peak frequency between γ = 3.321 and γ = 1.0 near 41 kHz is about 0.73 kHz.

As shown in Figures 57, when the mass ratio γ is large, the inertia forces of PZT patches and adhesive layers can influence the admittance signature greatly, even in the low frequency range, and the inertia term must be considered; when the mass ratio γ is small, the inertial term also need to be considered in the high frequency range. Because the EMI technique uses high-frequency alternating current, the proposed EMI model which takes inertial term into account will predict more accurate results.

Consider a beam with a crack located at x = 100. The parameters of the beam are listed in Table 1, and the beam thickness hb = 5 mm. the depth of crack is hc, the relative depth ξ = hc / hb. To verify the reliability of the proposed method, the admittance signatures of PZT patch from Ansys are compared with the data obtained by the analytical model. The crack is modeled in the ANSYS 10.0 as a slot with a width of 0.1 mm. the size of mesh near the slot is smaller than that of other areas (Figure 8). The results are shown in Figure 9. From Figure 9, it can be seen that the proposed method can predict admittance signatures of damaged structural as well as the FEA method. In order to study the influence of the crack on the admittance signature of the PZT patch, the admittance signatures are calculated when ξ = 0.05, ξ = 0.1, ξ = 0.2 and ξ = 0.4, respectively. The corresponding results are shown in Figure 10, Figure 11 and Figure 12. The resonant peaks of the admittance signatures shift towards the left with increasing crack depth. The influence of the crack increased with increasing frequency. These results are in good accord with the experimental phenomena in reference [30] and reference [33]. In fact, the change of admittance signature and the decrease of the peek frequencies reflect the decreasing local stiffness due to the crack, so the EMI technique can be used to identify the crack damage. From the changes of the admittance signature caused by the appearance of damage, both the location and quantity of the damage can be identified by using a certain back-calculation algorithm.

The crack has little influence on the admittance signature in some frequency bands, such as 45 kHz–47 kHz in Figure 10, 80 kHz–83 kHz in Figure 11 and 191 kHz–193 kHz in Figure 12. The frequency band is related to the crack location. The same phenomenon was also found by Youdi [33]. The reason is that if the crack is located at the node of a displacement modal, the crack has no influence on the modal frequency; otherwise, the modal frequency will reduce with the increasing crack depth. Hence, in order to improve the damage identification accuracy, the sweeping frequency band needs to contain enough peek frequencies to avoid the error due to the absence of influence of the crack around some frequencies. From the signatures of different damage extent, we can calculate the RMSD (root mean square deviation) value as follows:

RMSD = i = 1 n ( Y i d Y i u ) 2 i = 1 n ( Y i u ) 2
where Y is the admittance of PZT, the superscripts d and u denote the signature of the damaged structure and undamaged structure, respectively, and n denotes the number of sample points.

The RMSD values are shown in Figure 10, Figure 11 and Figure 12. From these figures, it can be seen that for the same extent of damage, the RMSD value increases with increasing frequencies. In other words, a higher frequency is more sensitive to damage than a lower frequency. However, when the driven frequency is higher than 200 kHz, the signature can be easily affected by temperature and bonding layer, hence, the driving frequencies used in EMI techniques are normally less than 200 kHz [35]. According to the engineering experience with the EMI technique, the minimum RMSD value which can be viewed as a reliable indicator of the existence of structural damage can be set as 1% [36]. In this work, the crack damage with a depth hc = 0.5 mm (ξ = 0.1) can be detected, as the RMSD value in the frequency range 180 kHz–200 kHz is greater than 5%.

5. Conclusions

An EMI model of a cracked Timoshenko beam with a pair of PZT patches system has been developed by considering the inertial forces of both the PZT patches and the bonding layers. The theoretical analysis and numerical tests are focused on the influence of crack and the influence of inertial forces because of the beam transverse motion. Through numerical tests results, the following conclusions can be drawn:

  • The inertial forces of PZT patches and bonding layers produced by the transverse motion of beam can be ignored in the low frequency band but should be considered in the high frequency band, especially for a thin beam structure. Because the EMI technique employs high frequency, taking the inertial forces into account is necessary when monitoring a thin beam structure.

  • The admittance signature of the PZT patch can reflect the crack damage very well, especially in the high frequency band. In some frequency bands, the crack has little influence on the admittance signature, while in other frequency bands the peak frequency of the admittance signature decreases with increasing frequency. In order to improve the accuracy of damage identification, the high frequency band which contains many peak frequencies should be chosen.

Based on the proposed EMI model, future work is planned to identify the crack damage quantitatively. To quantify the crack damage in plates and shells, further research is needed to establish the corresponding EMI models for cracked plates and cracked shells.

References

  1. Park, G; Sohn, H; Farrar, CR; Inman, DJ. Overview of piezoelectric impedance-based health monitoring and path forward. Shock Vib. Dig 2003, 35, 451–463. [Google Scholar]
  2. Giurgiutiu, V; Reynolds, A; Rogers, CA. Experimental investigation of E/M impedance health monitoring for spot welded structural joints. J. Intell. Mater. Syst. Struct 1999, 10, 802–812. [Google Scholar]
  3. Park, G; Rutherford, AC; Sohn, H; Farrar, CR. An outlier analysis framework for impedance-based structural health monitoring. J. Sound Vib 2005, 286, 229–250. [Google Scholar]
  4. Ayres, JW; Lalande, F; Chaudhry, Z; Rogers, CA. Qualitative impedance-based health monitoring of civil infrastructures. Smart Mater. Struct 1998, 7, 599–605. [Google Scholar]
  5. Park, G; Cudney, HH; Inman, DJ. Feasibility of using impedance-based damage assessment for pipeline structures. Earthquake Eng. Struc. Dynam 2001, 30, 1463–1474. [Google Scholar]
  6. Giurgiutiu, V; Zagrai, A. Damage detection in thin plates and aerospace structures with the electro-mechanical impedance method. Struct. Health Monit 2005, 4, 99–118. [Google Scholar]
  7. Xing, KJ; Fritzen, CP. Monitoring of growing fatigue damage using the E/M impedance method. Key Eng. Mater 2007, 347, 153–158. [Google Scholar]
  8. Shin, SW; Oh, TK. Application of electro-mechanical impedance sensing technique for online monitoring of strength development in concrete using smart PZT patches. Constr. Build. Mater 2009, 23, 1185–1188. [Google Scholar]
  9. Mourna, JD; Steffen, V. Impedance-based health monitoring for aeronautic structures using statistical meta-modeling. J. Intell. Mater. Syst. Struct 2006, 17, 1023–1036. [Google Scholar]
  10. Daniel, MP. High Frequency Modeling and Experimental Analysis for Implementation of Impedance-Based Structural Health Monitoring. PhD Thesis, The Virginia Polytechnic Institute and State University, Blacksburg, VA, USA, 2006.
  11. Liang, C; Sun, FP; Rogers, CA. Coupled electromechanical analysis of adaptive material systems-determination of the actuator power consumption and system energy transfer. J. Intell. Mater. Syst. Struct 1994, 5, 12–20. [Google Scholar]
  12. Zhou, SW; Liang, C; Rogers, CA. An impedance-based system modeling approach for induced strain actuator-driven structures. J. Vib. Acoust 1996, 118, 323–331. [Google Scholar]
  13. Bhalla, S; Soh, CK. Structural health monitoring by piezo-impedance transducers. I. Modeling. J. Aerosp. Eng 2004, 17, 154–165. [Google Scholar]
  14. Bhalla, S; Soh, CK. Structural health monitoring by piezo-impedance transducers. II. Applications. J. Aerosp. Eng 2004, 17, 166–175. [Google Scholar]
  15. Yang, YW; Xu, JF; Soh, CK. Generic impedance-based model for structure-piezoceramic interacting system. J. Aerosp. Eng 2005, 18, 93–101. [Google Scholar]
  16. Annamdas, VGM; Soh, CK. An electromechanical impedance model of a piezoceramic transducer-structure in the presence of thick adhesive bonding. Smart Mater. Struct 2007, 16, 673–686. [Google Scholar]
  17. Annamdas, VGM; Soh, CK. Three-dimensional electromechanical impedance model for multiple piezoceramic transducers-structure interaction. J. Aerosp. Eng 2008, 21, 35–44. [Google Scholar]
  18. Crawley, EF; de Luis, J. Use of piezoelectric actuators as elements of intelligent structures. AIAA J 1987, 25, 1373–1385. [Google Scholar]
  19. Qing, XP; Chan, HL; Beard, SJ; Ooi, TK; Marotta, SA. Effect of adhesive on the performance of piezoelectric elements used to monitor structural health. Int. J. Adhes. Adhes 2006, 26, 622–628. [Google Scholar]
  20. Wang, XD; Huang, GL. Wave propagation in electro-mechanical structures: induced by surface-bonded piezoelectric actuators. J. Intell. Mater. Syst. Struct 2001, 12, 105–115. [Google Scholar]
  21. Wang, XD; Huang, GL. The coupled dynamic behavior of piezoelectric sensors bonded to elastic media. J. Intell. Mater. Syst. Struct 2006, 17, 883–894. [Google Scholar]
  22. Pietrzakowski, M. Active damping of beams by piezoelectric system: Effects of bonding layer properties. Int. J. Solids Struct 2001, 38, 7885–7897. [Google Scholar]
  23. Ong, CW; Yang, YW; Wong, YT; Bhalla, S; Lu, Y; Soh, CK. The effects of adhesive on the electro-mechanical response of a piezoceramic transducer coupled smart system. Proc. SPIE 2003, 5062, 241–247. [Google Scholar]
  24. Bhalla, S; Kumar, P; Gupta, A; Datta, TK. Simplified impedance model for adhesively bonded piezo-impedance transducers. J. Aerosp. Eng 2009, 22, 373–382. [Google Scholar]
  25. Yan, W; Lim, CW; Cai, JB; Chen, WQ. An electromechanical impedance approach for quantitative damage detection in Timoshenko beams with piezoelectric patches. Smart Mater. Struct 2007, 16, 1390–1400. [Google Scholar]
  26. Yan, W; Lim, CW; Chen, WQ; Cai, JB. Modeling of EMI response of damaged Mindlin-Herrmann rod. Int. J. Mech. Sci 2007, 49, 1355–1365. [Google Scholar]
  27. Naidu, ASK; Soh, CK. Identifying damage location with admittance signatures of smart piezo-transducers. J Intell Mater Syst Struct 2004, 15, 627–642. [Google Scholar]
  28. Tseng, KKH; Wang, L. Impedance-based method for nondestructive damage identification. J. Aerosp. Eng 2005, 131, 58–64. [Google Scholar]
  29. Yang, YW; Hu, Y. Electromechanical impedance modeling of PZT transducers for health monitoring of cylindrical shell structures. Smart Mater. Struct 2008, 17, 015005. [Google Scholar]
  30. Xu, JF; Yang, YW; Soh, CK. Electromechanical impedance-based structural health monitoring with evolutionary programming. J. Aerosp. Eng 2004, 17, 182–193. [Google Scholar]
  31. Chondros, TG; Dimarogonas, AD; Yao, J. A continuous cracked beam vibration theory. J. Sound Vibra 1998, 215, 17–34. [Google Scholar]
  32. Wang, DS; Zhu, HP; Jin, K; Chen, XQ. Impedance analysis and damage detection on piezoelectric smart beam. Chin. J. Solid Mech 2008, 29, 402–407. [Google Scholar]
  33. Kuang, YD. Study on the Electro-Mechanical Characteristics of the Piezoelectric-Actuated Non-Cracked and Cracked Beams. Ph.D. Thesis, Huazhong University of Science and Technology, Wuhan, China, 2007.
  34. Liu, W; Giurgiutiu, V. Finite element simulation of piezoelectric wafer active sensors for structural health monitoring with coupled-field elements. Proc. SPIE 2007, 6529, 715–726. [Google Scholar]
  35. Yang, YW; Lim, YY; Soh, CK. Practical issues related to the application of the electromechanical impedance technique in the structural health monitoring of civil structures: II. Numerical verification. Smart Mater. Struct 2008, 17, 035009. [Google Scholar]
  36. Yang, YW; Hu, YH; Lu, Y. Sensitivity of PZT impedance sensors for damage detection of concrete structures. Sensors 2008, 8, 327–346. [Google Scholar]
Figure 1. Beam with a pair PZT patches bonded on its surface.
Figure 1. Beam with a pair PZT patches bonded on its surface.
Sensors 11 07285f1 1024
Figure 2. Infinitesimal beam element with piezoelements and bonding layers.
Figure 2. Infinitesimal beam element with piezoelements and bonding layers.
Sensors 11 07285f2 1024
Figure 3. The equivalent model of the cracked beam.
Figure 3. The equivalent model of the cracked beam.
Sensors 11 07285f3 1024
Figure 4. FE model of beam and PZT patches in Ansys.
Figure 4. FE model of beam and PZT patches in Ansys.
Sensors 11 07285f4 1024
Figure 5. The influence of inertia terms of PZT patches and bonding layers when hb = 2.5 mm. (a) 40 kHz–60 kHz; (b) 80 kHz–100 kHz.
Figure 5. The influence of inertia terms of PZT patches and bonding layers when hb = 2.5 mm. (a) 40 kHz–60 kHz; (b) 80 kHz–100 kHz.
Sensors 11 07285f5 1024
Figure 6. The influence of inertia terms of PZT patches and bonding layers when hb = 5 mm. (a) 40 kHz–60 kHz; (b) 80 kHz–100 kHz.
Figure 6. The influence of inertia terms of PZT patches and bonding layers when hb = 5 mm. (a) 40 kHz–60 kHz; (b) 80 kHz–100 kHz.
Sensors 11 07285f6 1024
Figure 7. The influence of inertia terms of PZT patches and bonding layers when hb = 10 mm. (a) 40 kHz–60 kHz; (b) 80 kHz–100 kHz.
Figure 7. The influence of inertia terms of PZT patches and bonding layers when hb = 10 mm. (a) 40 kHz–60 kHz; (b) 80 kHz–100 kHz.
Sensors 11 07285f7 1024
Figure 8. FE model of beam with a crack in Ansys.
Figure 8. FE model of beam with a crack in Ansys.
Sensors 11 07285f8 1024
Figure 9. Comparison of admittance signatures against frequency plot between the ANSYS simulation and the analytical results. (a) 40 kHz–60 kHz; (b) 80 kHz–100 kHz.
Figure 9. Comparison of admittance signatures against frequency plot between the ANSYS simulation and the analytical results. (a) 40 kHz–60 kHz; (b) 80 kHz–100 kHz.
Sensors 11 07285f9 1024
Figure 10. The influence of different extent crack damage in 40 kHz–60 kHz. (a) the admittance signatures of different damage extent; (b) RMSD value of different damage extent.
Figure 10. The influence of different extent crack damage in 40 kHz–60 kHz. (a) the admittance signatures of different damage extent; (b) RMSD value of different damage extent.
Sensors 11 07285f10 1024
Figure 11. The influence of different extent crack damage in 80 kHz–100 kHz. (a) the admittance signatures of different damage extent; (b) RMSD value of different damage extent.
Figure 11. The influence of different extent crack damage in 80 kHz–100 kHz. (a) the admittance signatures of different damage extent; (b) RMSD value of different damage extent.
Sensors 11 07285f11 1024
Figure 12. The influence of different extent crack damage in 180 kHz–200 kHz. (a) the admittance signatures of different damage extent; (b) RMSD value of different damage extent.
Figure 12. The influence of different extent crack damage in 180 kHz–200 kHz. (a) the admittance signatures of different damage extent; (b) RMSD value of different damage extent.
Sensors 11 07285f12 1024
Table 1. Properties of the beam.
Table 1. Properties of the beam.
lb (mm)bb (mm)Eb (Pa)Gb (Pa)Poisson’s ratioρb (kg/m3)Damping ratio
30086.60E102.33E100.3327000.01
Table 2. Properties of the PZT patch.
Table 2. Properties of the PZT patch.
lp (mm)hp (mm)bp (mm)Ep (Pa)ρp (kg/m3)d31 (m/V)ε33 (F/m)δη
100.586.1E107750−1.71E-101.53E-80.020.03
Table 3. Properties of the adhesive layer.
Table 3. Properties of the adhesive layer.
la (mm)ha (mm)ba (mm)Ga (Pa)ρa (kg/m3)CbPoisson’s ratio
100.181.0E917001.0E-70.38
Table 4. γ vs. hb.
Table 4. γ vs. hb.
hb (mm)2.551030
γ3.3212.1161.5801.193
Table 5. Constitutive data of the PZT-5A.
Table 5. Constitutive data of the PZT-5A.
S(m2/N)d(C/N)ε (F/m)
S11 = S22 = 16.4E-12d31 = −171E-12ε11 = 1.53E-8
S12 = −5.74E-12d32 = −171E-12ε22 = 1.53E-8
S13 = S23 = −7.2E-12d33 = 374E-12ε33 = 1.50E-8
S33 = 18.8E-12d42 = 584E-12
S44 = S55 = 47.5E-12d52 = 584E-12
S66 = 44.3E-12

Share and Cite

MDPI and ACS Style

Zhang, Y.; Xu, F.; Chen, J.; Wu, C.; Wen, D. Electromechanical Impedance Response of a Cracked Timoshenko Beam. Sensors 2011, 11, 7285-7301. https://doi.org/10.3390/s110707285

AMA Style

Zhang Y, Xu F, Chen J, Wu C, Wen D. Electromechanical Impedance Response of a Cracked Timoshenko Beam. Sensors. 2011; 11(7):7285-7301. https://doi.org/10.3390/s110707285

Chicago/Turabian Style

Zhang, Yuxiang, Fuhou Xu, Jiazhao Chen, Cuiqin Wu, and Dongdong Wen. 2011. "Electromechanical Impedance Response of a Cracked Timoshenko Beam" Sensors 11, no. 7: 7285-7301. https://doi.org/10.3390/s110707285

APA Style

Zhang, Y., Xu, F., Chen, J., Wu, C., & Wen, D. (2011). Electromechanical Impedance Response of a Cracked Timoshenko Beam. Sensors, 11(7), 7285-7301. https://doi.org/10.3390/s110707285

Article Metrics

Back to TopTop