Next Article in Journal
Aptamer-Based Single-Step Assay by the Fluorescence Enhancement on Electroless Plated Nano Au Substrate
Previous Article in Journal
Microfluidic-Based Measurement Method of Red Blood Cell Aggregation under Hematocrit Variations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact Mechanical Response of a 2-2 Cement-Based Piezoelectric Sensor Considering the Electrode Layer Effect

1
School of Transportation Science and Engineering, Beihang University, Beijing 100191, China
2
National Institute of Metrology, Beijing 100023, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(9), 2035; https://doi.org/10.3390/s17092035
Submission received: 10 July 2017 / Revised: 20 August 2017 / Accepted: 30 August 2017 / Published: 6 September 2017
(This article belongs to the Section Physical Sensors)

Abstract

:
Cement-based piezoelectric composite, has been widely used as a kind of smart material in structural health monitoring and active vibration control. However, transient dynamic loads such as impact loads may cause serious damage to the composite. Considering the electrode layer effect, this paper aimed to investigate the theoretical response of a 2-2 cement-based piezoelectric composite sensor subjected to an impact load. The vibration behaviors are analyzed by using the mode summation method and the virtual work principle. To simulate the impact load, transient haversine wave loads are assumed in the numerical simulation. Close agreements between theoretical and numerical solutions are found for peak transient haversine wave loads larger than 500 kPa, therefore proving the validity of the theory. Moreover, the influence of the electrode material and geometrical parameters on the dynamic characteristics of this sensor are considered. The present work should be beneficial to the design of this kind of sensor by taking into account the electrode layer effect.

1. Introduction

With the advances in smart or intelligent structure technology, piezoelectric composites are widely used in structural health monitoring [1,2,3]. These smart devices are sensitive to the influences from the external environment and have extraordinary compatibility with the most popular construction materials used in civil engineering such as cement and concrete. The design and fabrication of composite materials open new avenues for optimizing the electrical, magnetic, and mechanical properties of sensors for specific applications [4,5,6].
According to the different connectivity of each phase in the piezoelectric composite, piezoelectric composites can be divided into 10 basic types such as 0-0, 0-1, 0-3, 2-2 and so on. The first number represents the piezoelectric phase and the second number represents the non-piezoelectric phase. A 2-2 cement-based piezoelectric composite is the composite where a two-dimensional piezoelectric plate is embedded in a two-dimensional cement matrix [7]. Li and Zhang fabricated 2-2 cement-based piezoelectric composites using hardened Portland cement and studied the sensing performance of such composites at low frequency [8]. Investigation of the dielectric and acoustic impedance properties of 2-2 barium zirconate titanate-Portland cement composites reveals that the composites have higher piezoelectric voltage coefficients and lower acoustic impedance values than pure ceramic [9]. Cheng et al. studied the effects of composite thickness on the dielectric, piezoelectric and electromechanical properties of the composite [10]. The receiving piezoelectric transducer can be fabricated by decreasing the thickness. Moreover, mathematical models to describe the deformation and electrical behaviors have also been studied by numerous researchers. For example, based on the theory of elasticity, static analyses of piezoelectric curved composites and multi-layered piezoelectric cantilevers were performed by Zhang et al. [11,12,13,14], who further established an analytical model of the dynamic properties of the 2-2 cement-based piezoelectric transducer subjected to a uniformly distributed harmonic load and external harmonic electrical potential [15]. Chen et al. analyzed the free vibration of laminates beams by using the state space method and the differential quadrature method [16]. Wang and Shi presented the dynamic behaviors of piezoelectric composite stack transducers and discussed the influence of the electrode thickness on their dynamic characteristics [17].
Few of the published papers have studied the dynamic performance of piezoelectric composite sensors subjected to impact loads, especially the theoretical aspect. Zhang et al. studied the dynamic properties of piezoelectric structures subjected to simpact load and gave theoretical solutions of the mechanical and electrical behaviors of the piezoelectric structures [18]. Yin et al. presented a mixed finite element formulation for modelling the behavior of a piezoelectric composite and analyzed the response of distributed sensors made of PVDF film when the composite was subjected to low velocity impacts [19]. Ueda studied the transient response of a functionally graded piezoelectric material strip with a vertical crack under the action of normal impacts [20,21]. In many engineering applications, piezoelectric composite may experience transient dynamic loads, which could induce certain hidden damages (matrix cracking, fiber breakage, etc.). It is therefore important to understand the dynamic behaviors of cement-based piezoelectric sensors.
In the present study, we aim to provide a mathematical model to describe the impact mechanical response of a 2-2 cement-based piezoelectric composite subjected to impact load. The fundamental equations and wave equations are summarized in Section 2 based upon piezo-elasticity. In the following section, the vibration behaviors of the composite are obtained by using the mode summation method and the principle of virtual work. In Section 4, the transient harversine wave load is used for the numerical simulation. The numerical results verified the analytical solutions. Moreover, the influence of the material and geometrical parameters on mechanical and electrical behaviors of sensor is discussed. It can be seen that the proposed model would provide guidance for sensor structure design, material selection and impact load design in simulations and experiments.

2. Basic Equations

The 2-2 cement-based piezoelectric sensor consists of a piezoelectric layer P # 3 , elastic electrode layers E # 2 , E # 4 and cement layers C # 1 , C # 5 , which are arranged in an alternating manner as shown in Figure 1.
The thickness of layer i are determined by l i + 1 l i ( i   =   1 ,   2 ,   3 ,   4 ,   5 ). The composite is excited longitudinally by the impact load δ ( t ) at its free end. The piezoelectric layer is polarized along its thickness direction and the electrode surface is perpendicular to the thickness direction of the piezoelectric layer. In order to analyze the transient phase of the composite vibration, the following assumptions are adopted: (1) all of the constitutive materials are linear and the deformations are small; (2) one-dimensional assumption. When the composite is deformed, the cross-section maintain a plane and the axial stress is uniformly distributed on the cross-section; (3) damping is not considered. For the piezoelectric layer P # 3 , the body force and body charge are neglected, the governing equations can be expressed as follows:
{ 2 w p 3 t 2 = 1 ρ p σ p 3 z σ p 3 = C 33 p ε p 3 e 33 E z ε p 3 = w p 3 z
{ D z = e 33 ε p 3 + ε 33 S E z E z = ϕ z D z z = 0
where σ p 3 , ε p 3 and w p 3 are the stress, strain, and displacement of the piezoelectric layer P # 3 along z-direction, respectively. D z and φ represent the electric displacement and electric potential in the along z-direction. C 33 p ,    e 33 and ε 33 S are the elastic stiffness coefficient, piezoelectric coefficient and permittivity coefficient of the piezoelectric material, respectively. ρ p is the density of the piezoelectric material. Combining Equations (1) and (2), the following equation are obtained:
{ ρ p 2 w p 3 t 2 = C 33 p 2 w p 3 z 2 + e 33 2 ϕ z 2 ε 33 S 2 ϕ z 2 = e 33 2 w p 3 z 2
Combining Equation (3), one gets the following wave equation:
2 w p 3 t 2 C a 2 2 w p 3 z 2 = 0
where C a   = E 0 / ρ p   , E 0   =   C 33 p   +   e 33 2 / ε 33 S . C a represents the propagation velocity of the vibration wave in the piezoelectric material.
For the cement layer C # i   ( i   =   1 ,   5 ) , the body charge are neglected, and the basic equations can be expressed as follow:
{ 2 w c i t 2 = 1 ρ c σ c i z σ c i = C 33 c ε c i ε c i = w c i z
where σ c i , ε c i and w c i are the stress, strain, and displacement of the cement layer C # i along z-direction, respectively. C 33 c and ρ c are the elastic stiffness coefficient and density of the cement material. Combining Equation (5), one gets the following wave equation:
2 w c i t 2 C b 2 2 w c i z 2 = 0
where C b   =   C 33 c ρ c and represents the propagation velocity of the vibration wave in the cement material.
For the elastic electrode layer E # i   ( i   =   2 ,   4 ) , the body charge are neglected, and the basic equations can be expressed as follow:
{ 2 w E i t 2 = 1 ρ E σ E i z σ E i = C 33 E ε E i ε E i = w E i z
where σ E i , ε E i and w E i are the stress, strain, and displacement of the elastic electrode layer E # i along z-direction, respectively. C 33 E and ρ E are the elastic stiffness coefficient and density of the elastic electrode material. Similarly, wave equation can be easily obtained:
2 w E i t 2 C c 2 2 w E i z 2 = 0
where C c   =   C 33 E ρ E and represents the propagation velocity of the vibration wave in the elastic electrode material.

3. Accurate Vibration Analysis of 2-2 Cement Based Piezoelectric Composite

Figure 1 shows the 2-2 cement-based piezoelectric composite sensor subjected to the impact load. In order to obtain the accurate vibration solution, we use the mode summation method and the principle of the virtual work in this section. It is assumed that the solutions for the displacements are:
{ w p i ( z , t ) = Z p i ( z ) T n ( t ) ;    i = 3 w c i ( z , t ) = Z c i ( z ) T n ( t ) ;    i = 1 , 5 w E i ( z , t ) = Z E i ( z ) T n ( t ) ;    i = 2 , 4
where Z c i ( z ) , Z E i ( z ) and Z p i ( z ) are normal modes of the cement, piezoelectric and elastic electrode layers for the longitudinal vibration. T n is the function of t which can be determined by the initial conditions.
Substituting the above equations into Equations (4), (6) and (8), we obtain:
{ Z p i ( z ) + λ i Z p i ( z ) = 0 ; l i 1 z l i ,    i = 3 Z c i ( z ) + λ i Z c i ( z ) = 0 ; l i 1 z l i ,    i = 1 , 5 Z E i ( z ) + λ i Z E i ( z ) = 0 ; l i 1 z l i ,    i = 2 , 4
{ T n ( t ) + λ i C a 2 T n ( t ) = 0 ; t 0 ,    i = 3 T n ( t ) + λ i C b 2 T n ( t ) = 0 ; t 0 ,    i = 1 , 5 T n ( t ) + λ i C c 2 T n ( t ) = 0 ; t 0 ,    i = 2 , 4
where λ i are eigenvalues.
By solving Equation (10), the solutions of the normal modes are shown as follow:
{ Z p i ( z ) = a i n cos λ i n z + b i n sin λ i n z ; l i 1 z l i ,    i = 3 Z c i ( z ) = a i n cos λ i n z + b i n sin λ i n z ; l i 1 z l i ,    i = 1 , 5 Z E i ( z ) = a i n cos λ i n z + b i n sin λ i n z ; l i 1 z l i ,    i = 2 , 4
where a i n and b i n will be determined by the boundary conditions.
In order to solve a i n and b i n , the boundary and connecting conditions are written as (including the electrical boundary conditions and the initial conditions):
{ (13a) w c 1 ( 0 , t ) = 0 ; w c 5 ( l 5 , t ) z = 0 ; t 0 (13b) w c 1 ( l 1 , t ) = w E 2 ( l 1 , t ) ;    w E 2 ( l 2 , t ) = w p 3 ( l 2 , t ) ; (13c) w p 3 ( l 3 , t ) = w E 4 ( l 3 , t ) ;    w E 4 ( l 4 , t ) = w c 5 ( l 4 , t ) ; t 0 (13d) C 33 c w c 1 ( l 1 , t ) z = C 33 E w E 2 ( l 1 , t ) z ;    C 33 E w E 2 ( l 2 , t ) z = E 0 w p 3 ( l 2 , t ) z ; (13e) E 0 w p 3 ( l 3 , t ) z = C 33 E w E 4 ( l 3 , t ) z ;    C 33 E w E 4 ( l 4 , t ) z = C 33 c w c 5 ( l 5 , t ) z ;    t 0 (13f) w p i ( z , 0 ) = w p i ( z , 0 ) t = 0 ;    i = 3 (13g) w c i ( z , 0 ) = w c i ( z , 0 ) t = 0 ;    i = 1 , 5 (13h) w E i ( z , 0 ) = w E i ( z , 0 ) t = 0 ;    i = 2 , 4 (13i) ϕ i | z = l 2 = 0 ; D 0
Substitution of Equation (9) into Equation (13) leads to:
{ (14a) Z c 1 ( l 1 ) = Z E 2 ( l 2 ) ;    C 33 c Z c 1 ( l 1 ) = C 33 E Z E 2 ( l 1 ) (14b) Z E 2 ( l 2 ) = Z p 3 ( l 2 ) ;    C 33 E Z E 2 ( l 2 ) = E 0 Z p 3 ( l 2 ) (14c) Z p 3 ( l 3 ) = Z E 4 ( l 3 ) ;    E 0 Z p 3 ( l 3 ) = C 33 E Z E 4 ( l 3 ) (14d) Z E 4 ( l 4 ) = Z c 5 ( l 4 ) ;    C 33 E Z E 4 ( l 4 ) = C 33 c Z c 5 ( l 4 ) (14e) Z c 1 ( 0 ) = 0 ;    Z c 5 ( l 5 ) = 0
By combining Equations (11), (13b,c) and (14a–c), one obtains the following relation:
λ 1 C b = λ 2 C c = λ 3 C a = λ 4 C c = λ 5 C b
Substitution of Equation (12) into Equation (14) leads to the following equations:
{ (16a) a 1 n = 0 (16b) a i n cos λ i n l i + b i n sin λ i n l i a i + 1 , n cos λ i + 1 , n l i b i + 1 , n sin λ i + 1 n l i = 0 (16c) C i λ i n ( a i n sin λ i n l i + b i n cos λ i n l i ) C i + 1 λ i + 1 , n ( a i + 1 , n sin λ i + 1 , n l i + b i + 1 , n cos λ i + 1 , n l i ) = 0 (16d) a 5 n λ 5 n sin λ 5 n l 5 + b 5 n λ 5 n cos λ 5 n l 5 = 0
here   i = 2, 3, 4, and correspondingly C 1   =   C 5   =   C 33 c ,   C 2   =   C 4   =   C 33 E ,   C 3   =   E 0 .
By solving the linear Equation (16b,c) with two unknowns, the expressions of a i n and b i n are obtained:
{ a i n = a i 1 , n ( cos λ i n l i 1 cos λ i 1 , n l i 1 + C i 1 λ i 1 , n C i λ i n sin λ i n l i 1 sin λ i n l i ) +                   b i 1 , n ( sin λ i 1 , n l i 1 cos λ i n l i 1 C i 1 λ i 1 , n C i λ i n cos λ i 1 , n l i 1 sin λ i n l i 1 ) b i n = a i 1 , n ( cos λ i 1 , n l i 1 sin λ i n l i 1 C i 1 λ i 1 , n C i λ i n sin λ i 1 , n l i 1 cos λ i n l i 1 ) +               b i 1 , n ( sin λ i n l i 1 sin λ i 1 , n l i 1 + C i 1 λ i 1 , n C i λ i n cos λ i 1 , n l i 1 cos λ i n l i 1 )
here i =   2, 3, 4 and we assume that b 1 n   =   1 .
It should be noted that when i =   4, there are three linear Equation (16b–d) with two unknowns:
{ a 5 n cos λ 5 n l 4 + b 5 n sin λ 5 n l 4 = a 4 n cos λ 4 n l 4 + b 4 n sin λ 4 n l 4 a 5 n sin λ 5 n l 4 b 5 n cos λ 5 n l 4 = C 33 E λ 4 n C 33 C λ 5 n ( a 4 n sin λ 4 n l 4 b 4 n cos λ 4 n l 4 ) a 5 n λ 5 n sin λ 5 n l 5 + b 5 n λ 5 n cos λ 5 n l 5 = 0
Nonhomogeneous equations have nontrivial solutions, which requires:
| cos λ 5 n l 4 sin λ 5 n l 4 a 4 n cos λ 4 n l 4 + b 4 n sin λ 4 n l 4 sin λ 5 n l 4 cos λ 5 n l 4 C 33 E λ 4 n C 33 c λ 5 n ( a 4 n sin λ 4 n l 4 b 4 n cos λ 4 n l 4 ) sin λ 5 n l 5 cos λ 5 n l 5 0 | = 0
By solving the above determinant, we obtain the following equation:
a 4 n ( C 33 E λ 4 n C 33 c λ 5 n sin λ 4 n l 4 cos λ 5 n h 5 + cos λ 4 n l 4 sin λ 5 n h 5 ) + b 4 n ( C 33 E λ 4 n C 33 c λ 5 n cos λ 4 n l 4 cos λ 5 n h 5 + sin λ 4 n l 4 sin λ 5 n h 5 ) = 0
where h 5   =   l 5 l 4 .
Based on Equation (15), we define:
(21a) λ 1 n C b = λ 2 n C c = λ 3 n C a = λ 4 n C c = λ 5 n C b = λ n ¯ (21b) h 1 = l 1 , h 2 = l 2 l 1 , h 3 = l 3 l 2 , h 4 = l 4 l 3 , h 5 = l 5 l 4 (21c) t 1 = l 1 C b t 2 = h 2 C c t 3 = h 3 C a t 4 = h 4 C c , t 5 = h 5 C b , T 0 = t 1 + t 2 + t 3 + t 4 + t 5    λ n ¯ T 0 = λ n ¯ ¯
By combining Equation (21), one obtains:
λ 1 n = λ 5 n = λ n ¯ ¯ A , λ 2 n = λ 4 n = λ n ¯ ¯ B , λ 3 n = λ n ¯ ¯ C
where A   =   T 0 C b ,   B   =   T 0 C c , C   =   T 0 C a .
By substituting Equation (22) into Equation (20), the value of λ n ¯ ¯ is given by:
a 4 n ( A · C 33 E B · C 33 c sin λ n ¯ ¯ B l 4 cos λ n ¯ ¯ A h 5 + cos λ n ¯ ¯ B l 4 sin λ n ¯ ¯ A h 5 ) + b 4 n ( A · C 33 E B · C 33 c cos λ n ¯ ¯ B l 4 cos λ n ¯ ¯ A h 5 + sin λ n ¯ ¯ B l 4 sin λ n ¯ ¯ A h 5 ) = 0
After we obtain the value of λ n ¯ ¯ by Equation (23), λ 1 n , λ 2 n , λ 3 n ,   λ 4 n , and λ 5 n can be obtained by Equation (21a).
To solve the expressions of T n ( t ) in Equation (9), we consider the initial condition of the sensor. By substituting Equation (9) into Equation (13f–h), we have:
T n ( t ) = T ˙ n ( t ) = 0
At t   =   0 , there are three forces, namely the inertia forces, the elasticity force due to the deformation in each element of the composite, and the impact force δ ( t ) loaded at the free end. By taking any displacement δ u satisfying the boundary and connecting conditions as a virtual displacement, according to Equation (12), δ u can be given as follows:
{ (25a) δ u p i = Z p i ( z ) = a i n cos λ i n z + b i n sin λ i n z ; l i 1 z l i ,    i = 3 (25b) δ u c i = Z c i ( z ) = a i n cos λ i n z + b i n sin λ i n z ; l i 1 z l i ,    i = 1 , 5 (25c) δ u E i = Z E i ( z ) = a i n cos λ i n z + b i n sin λ i n z ; l i 1 z l i ,    i = 2 , 4
The virtual work δ W i of the inertia force on the virtual displacement is expressed as follows:
δ W i = 0 l 1 ( ρ c A d x ) · 2 w c 1 t 2 · δ u c 1 + l 2 l 2 ( ρ E A d x ) · 2 w E 2 t 2 · δ u E 2 + l 2 l 3 ( ρ p A d x ) · 2 w p 3 t 2 · δ u p 3 + l 3 l 4 ( ρ E A d x ) · 2 w E 4 t 2 · δ u E 4 + l 4 l 5 ( ρ c A d x ) · 2 w c 5 t 2 · δ u c 5
By substituting Equations (9) and (25) into Equation (26), δ W i can be obtained as:
δ W i = A T ¨ n ( t ) · ( ρ c B 1 + ρ E B 2 + ρ p B 3 + ρ E B 4 + ρ c B 5 )
where B i   =   l i 1 l i ( a i n cos λ i n z   +   b i n sin λ i n z ) 2 d z , i =   1, 2, 3, 4, 5.
The virtual work δ W E of the elastic force on the virtual displacement is expressed as follows:
δ W E = 0 l 1 ( C 33 c A d x ) · 2 w c 1 z 2 · δ u c 1 + l 2 l 2 ( C 33 E A d x ) · 2 w E 2 z 2 · δ u E 2 + l 2 l 3 ( E 0 A d x ) · 2 w p 3 z 2 · δ u p 3 + l 3 l 4 ( C 33 E A d x ) · 2 w E 4 z 2 · δ u E 4 + l 4 l 5 ( C 33 c A d x ) · 2 w c 5 z 2 · δ u c 5
By substituting Equations (9) and (25) into Equation (28), δ W E can be obtained as:
δ W E = A T n ( t ) · ( λ 1 n C 33 c B 1 + λ 2 n C 33 E B 2 + λ 3 n E 0 B 3 + λ 4 n C 33 E B 4 + λ 5 n C 33 c B 5 )
In order to obtain the virtual work δ W δ ( t ) of the impact load   δ ( t ) , by substituting l 5 into Equation (25b) to obtain the virtual displacement at the free end, δ W δ ( t )   is expressed as:
δ W δ ( t ) = δ ( t ) A ( a 5 n cos λ 5 n l 5 + b 5 n sin λ 5 n l 5 )
Summation of δ W i , δ W E and δ W δ ( t ) gives the total virtual work, equating it to zero:
δ W i + δ W E + δ W δ ( t ) = 0
Substituting Equations (27), (29) and (30) into Equation (31), one obtains:
T ¨ n ( t ) + λ n ¯ T n ( t ) = δ ( t ) M n · ( a 5 n cos λ 5 n l 5 + b 5 n sin λ 5 n l 5 )
where M n   =   ρ c B 1   +   ρ E B 2   +   ρ p B 3   +   ρ E B 4   +   ρ c B 5 and the Dirac function δ ( t ) can be expressed as:
δ ( t ) = { , t = 0 0 , t 0
Writing the solution of the Equation (32) in the form of the Duhamel’s integral, one obtains:
T n ( t ) = ( a 5 n cos λ 5 n l 5 + b 5 n sin λ 5 n l 5 ) M n λ n ¯ · 0 t δ ( τ ) sin λ n ¯ ( t τ ) d τ = ( a 5 n cos λ 5 n l 5 + b 5 n sin λ 5 n l 5 ) M n λ n ¯ · sin λ n ¯ t
Therefore, the accurate vibration analysis of the 2-2 cement-based piezoelectric composite sensor excited by the impact load can be obtained as:
Displacement functions:
{ w c i ( z , t ) = n = 1 D n λ n ¯ · sin λ n ¯ t ( a i n cos λ i n z + b i n sin λ i n z ) ; t 0 , l i 1 z l i ,    i = 1 , 5 w E i ( z , t ) = n = 1 D n λ n ¯ · sin λ n ¯ t ( a i n cos λ i n z + b i n sin λ i n z ) ; t 0 , l i 1 z l i ,    i = 2 , 4 w p i ( z , t ) = n = 1 D n λ n ¯ · sin λ n ¯ t ( a i n cos λ i n z + b i n sin λ i n z ) ; t 0 , l i 1 z l i ,    i = 3
where D n   =   ( a 5 n cos λ 5 n l 5   +   b 5 n sin λ 5 n l 5 ) / M n ;
Stress functions:
{ σ c i ( z , t ) = n = 1 C 33 c D n λ i n λ n ¯ · sin λ n ¯ t ( b i n cos λ i n z a i n sin λ i n z ) ; t 0 , l i 1 z l i ,    i = 1 , 5 σ E i ( z , t ) = n = 1 C 33 E D n λ i n λ n ¯ · sin λ n ¯ t ( b i n cos λ i n z a i n sin λ i n z ) ; t 0 , l i 1 z l i ,    i = 2 , 4 σ p i ( z , t ) = n = 1 E 0 D n λ i n λ n ¯ · sin λ n ¯ t ( b i n cos λ i n z a i n sin λ i n z ) ; t 0 , l i 1 z l i ,    i = 3
Strain functions:
{ ε c i ( z , t ) = n = 1 D n λ i n λ n ¯ · sin λ n ¯ t ( b i n cos λ i n z a i n sin λ i n z ) ; t 0 , l i 1 z l i ,    i = 1 , 5 ε E i ( z , t ) = n = 1 D n λ i n λ n ¯ · sin λ n ¯ t ( b i n cos λ i n z a i n sin λ i n z ) ; t 0 , l i 1 z l i ,    i = 2 , 4 ε p i ( z , t ) = n = 1 D n λ i n λ n ¯ · sin λ n ¯ t ( b i n cos λ i n z a i n sin λ i n z ) ; t 0 , l i 1 z l i ,    i = 3
Velocity functions:
{ v c i ( z , t ) = n = 1 D n · cos λ n ¯ t ( a i n cos λ i n z + b i n sin λ i n z ) ; t 0 , l i 1 z l i ,    i = 1 , 5 v E i ( z , t ) = n = 1 D n · cos λ n ¯ t ( a i n cos λ i n z + b i n sin λ i n z ) ; t 0 , l i 1 z l i ,    i = 2 , 4 v p i ( z , t ) = n = 1 D n · cos λ n ¯ t ( a i n cos λ i n z + b i n sin λ i n z ) ; t 0 , l i 1 z l i ,    i = 3
Acceleration functions:
{ a c i ( z , t ) = n = 1 D n λ n ¯ · sin λ n ¯ t ( a i n cos λ i n z + b i n sin λ i n z ) ; t 0 , l i 1 z l i ,    i = 1 , 5 a E i ( z , t ) = n = 1 D n λ n ¯ · sin λ n ¯ t ( a i n cos λ i n z + b i n sin λ i n z ) ; t 0 , l i 1 z l i ,    i = 2 , 4 a p i ( z , t ) = n = 1 D n λ n ¯ · sin λ n ¯ t ( a i n cos λ i n z + b i n sin λ i n z ) ; t 0 , l i 1 z l i ,    i = 3
The electric potential of the piezoelectric layer:
ϕ ( z , t ) = n = 1 e 33 D n ε 33 S λ n ¯ · sin λ n ¯ t [ a 3 n ( cos λ 3 n z cos λ 3 n l 2 ) + b 3 n ( sin λ 3 n z sin λ 3 n l 2 ) ] ; t 0 , l 2 z l 3
Electric field intensity of piezoelectric layer:
E ( z , t ) = n = 1 e 33 D n λ 3 n ε 33 S λ n ¯ · sin λ n ¯ t ( b 3 n cos λ 3 n z a 3 n sin λ 3 n z ) ; t 0 , l 2 z l 3
Till now, all the mechanical and electrical solutions have been obtained by using the mode summation method and the principle of virtual work.

4. Comparisons and Discussions

In this section, by comparing the theoretical and the numerical solutions, the modal interception of the mode summation method, and the critical impact load value in the numerical simulation are discussed. We also analyze the influence of the material and geometrical parameters on the mechanical and electrical behaviors of the sensor.

4.1. Comparions between the Theoretical Solutions and the Numerical Solutions

The thickness of the sensor is taken as 0.02 m. The material properties used in Li’s experiments [22] are adopted and are listed in Table 1. ANSYS is used for the numerical simulation, where an analytical model of the size 0.001   m   ×   0.001   m   ×   0.02   m is considered. The piezoelectric layer is defined as a Solid5 unit, cement layers and electrode layers are both defined as Solid45 units. The interlayer contact is glue. The unit partition of the model is divided into five segments along the x-axis and y-axis, and 100 segments along the z-axis by using the free meshing method. The upper and lower surfaces of the piezoelectric layer in the z-axis direction are subjected to the piezoelectric coupling. The electric potential of the lower surface of the piezoelectric layer is set to zero. The model is loaded and solved after the symmetrical boundary conditions are set on the four sides of the model. The impact load Q ( t )   used in the numerical simulation has the form of a transient haversine wave and is shown in Figure 2. The relation + Q ( t ) d t   =   1 holds. The theoretical solutions are worked out by using common programming language with only about 6s. The calculation time of ANSYS simulation is 13 min.
Figure 3a,b show the influences of the number of modes n involved in the summation on the displacement. It is found that the participation of the high-order modes would cause slight increase in the displacement amplitude and longer vibration duration. There are not much dramatic changes in the modal curves when high-order modes join in. The theoretical solutions of the electric potential and stress are plotted in Figure 3c,d, where n   =   3 .
In order to obtain the critical impact load value in the numerical simulation, comparisons between theoretical solutions and the numerical solutions with different peak value of the impact load Q ( t )   are shown in Figure 4. The peak value of Q ( t ) takes 300 kPa, 400 kPa, 500 kPa and 600 kPa respectively. It is observed that the numerical solution is closer to the theoretical solution as the peak value of Q ( t ) becomes larger. In particular, the electric potential amplitude in the simulation with 300 kPa peak value is approximately 50 V less than the theoretical amplitude, and the stress amplitude with 400 kPa peak value is less about 400 kPa than the theoretical amplitude (see Figure 4c,d). Therefore, it is appropriate to use the impact load Q ( t ) with 500 kPa or larger peak value in the numerical simulation.

4.2. Material and Geometrical Parameters of Composite Properties

It is assumed that three types of elastic electrode layers are made of H62-brass, aluminum and (gold-tin, 80% wt % Au-20 wt % Sn), respectively. The material constants of H62-brass [17], aluminum [23] and (gold-tin, 80% wt % Au-20 wt % Sn) [24] are listed in Table 2.
Figure 5 demonstrates the displacement w c ( l 5 , t 0 ) , electric potential ϕ ( l 3 , t 0 )   and stress σ c ( l 3 , t 0 ) versus the thickness ratio h 3 / h 2 . Here t 0   =   0.70   ×   10 5   s, the first displacement amplitude appeared at t 0 (shown in Figure 2a). h 3 and h 2 are the thickness of the piezoelectric layer and electrode layer, respectively. Meanwhile, the influence of the elastic electrode material on the displacement, electric potential and stress are also shown. By keeping the total thickness l 5 and the thickness of piezoelectric layer h 3   constant, it is observed that when h 3 / h 2   =   4 , the internal stresses and electric potential are minimal. The tip displacement is larger with a thinner electrode layer. This is an indication that the mechanical and electrical characteristics can be enhanced by tailoring the geometry of the sensor. The impact of aluminum on the mechanical and electrical behaviors are relatively larger.
Figure 6 shows the displacement w c ( l 5 , t 0 ) , electric potential ϕ ( l 3 , t 0 ) and stress σ c ( l 1 , t 0 ) as functions of C 33 p respectively. The tip displacement and electric potential both decrease, while the internal stress increases, as C 33 p grows. It can be concluded that a sensor with smaller value of C 33 p can provide larger displacement and electric potential, therefore further causing smaller internal stress.
Figure 7 illustrates the displacement w c ( l 5 , t 0 ) , electric potential ϕ ( l 3 , t 0 ) and stress σ c ( l 1 , t 0 ) as functions of e 33 respectively. It can be easily seen that the influence of e 33 on the electric potential is more obvious. The larger e 33 is, the larger the electric potential generated is.
Figure 8 plots the displacement w c ( l 5 , t 0 ) , electric potential ϕ ( l 3 , t 0 ) and stress σ c ( l 1 , t 0 ) as functions of ε 33 S / ε , respectively. It can be described that the electric potential decreases as ε 33 S / ε increases and the trend flattens when ε 33 S / ε is larger than 50. Figure 6, Figure 7 and Figure 8 also reveal the influence of the thickness ratio h 3 / l 5 on the mechanical and electrical behaviors. Composites with thicker piezoelectric layers can generate larger electric potential as well as larger internal stresses and smaller tip displacements. A composite with a thinner piezoelectric layer has better mechanical behaviors. h 3 / l 5 of 0.50 is a relatively good geometrical parameters for the sensor.
The influence of the piezoelectric material on the displacement w ( z , t 0 ) electric potential ϕ ( z , t 0 ) stress σ ( z , t 0 ) are shown in Figure 9. PZT-5H, PZT-4 and PVDF [25] are discussed as common piezoelectric materials, and a comparison with Li’s experiments [22] is presented. The material constants of PZT-5H, PZT-4 and PVDF are listed in Table 3.
For PVDF, obvious differences in mechanical and electrical behaviors are observed as compared to the case for PZT, as shown in Figure 9. z is the direction along the thickness direction of the composite. The tip displacement and electric potential are large enough for the sensor under impact load. Meanwhile, the resulting stress is relatively small compared with PZT.
Figure 10 demonstrates the influence of material parameters on the fundamental frequency. It can be seen that the fundamental frequency ω increases as C 33 p , e 33 and ε 33 S / ε grow, decreases as ρ p increases. Figure 10 also shows that the thicker the piezoelectric layer, the greater the effect of the piezoelectric coefficients on the frequency (larger rate of change).

5. Conclusions

Based on the theory of piezo-elasticity, an accurate mechanical and electrical analysis of the 2-2 cement-based piezoelectric sensor are presented in this paper. Theoretical solutions are obtained with the mode summation method and the principle virtual work. Through comparisons with numerical solutions, the following conclusions can be drawn:
(1) For theoretical solutions, the vibration modal curves of the sensor subjected to the impact load have no obvious change after the addition of high-order modes. It’s sufficient to analyze and summate the first three modes. Numerical simulations have good agreement with the theoretical solutions when the peak value of impact load Q ( t ) is larger than 500 kPa.
(2) By keeping the total thickness of the sensor and the thickness of piezoelectric layer h 3 constant, the sensor shows good mechanical properties with a thickness ratio h 3 / h 2   =   4 and good electrical property with a thickness ratio h 3 / h 2   =   2 . Aluminum as the elastic electrode material has a relatively large impact on the tip displacement, electric potential and internal stress.
(3) Through adjusting the thickness of layers and material parameters, the displacement, electric potential and stress of the sensor could be optimized. The coefficient C 33 p has obvious influence on w c ( l 5 , t 0 ) , σ c ( l 1 , t 0 ) and ϕ ( l 3 , t 0 ) ; while e 33 and   ε 33 S / ε have larger influence on ϕ ( l 3 , t 0 ) . For the sensor, smaller C 33 p and ε 33 S / ε , larger e 33 would provide better mechanical and electrical behaviors. PVDF as piezoelectric material can provide stronger electric power as well as causing smaller internal stress.
(4) The frequency of the composite could also be controlled by choosing different materials and tailoring the geometry of the composite. The thicker the piezoelectric layer, the greater the effect of the piezoelectric coefficients on the frequency (lager changing rate).
By analyzing the dynamic characteristics of the sensor, the present work would provide certain guidance for the sensor structure design, material selection and impact load design, both in simulations and experiments.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (51278517). Support is also supplied by the Fundamental Research Funds for the Central Universities (KG12002801).

Author Contributions

Taotao Zhang and Keping Zhang designed the model and worked on theoretical analyses. Yangchao Liao conceived of, designed and performed the simulation. Keping Zhang wrote the paper with the help of Wende Liu and Yangchao Liao.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ayres, J.W.; Lalande, F.; Chaudhry, Z.; Rogers, C.A. Qualitative impedance-based health monitoring of civil infrastructures. Smart Mater. Struct. 1998, 7, 599. [Google Scholar] [CrossRef]
  2. Qin, Y.; Mao, W.; Huang, Z.; Guo, D. Study on damping property of epoxy-matrix with PbZr0.47Ti0.53O3 piezoelectric film effect. In Proceedings of the 2nd International Conference on Smart Materials and Nanotechnology in Engineering, Weihai, China, 8–11 July 2009; p. 749356. [Google Scholar]
  3. Staszewski, W.; Boller, C.; Tomlinson, G.R. Health Monitoring of Aerospace Structures: Smart Sensor Technologies and Signal Processing; John Wiley & Sons: Hoboken, NJ, USA, 2004. [Google Scholar]
  4. Duan, W.H.; Wang, Q.; Quek, S.T. Applications of piezoelectric materials in structural health monitoring and repair: Selected research examples. Materials 2010, 3, 5169–5194. [Google Scholar] [CrossRef]
  5. Dong, B.; Liu, Y.; Qin, L.; Wang, Y.; Fang, Y.; Xing, F.; Chen, X. In-situ structural health monitoring of a reinforced concrete frame embedded with cement-based piezoelectric smart composites. Res. Nondestruct. Eval. 2016, 27, 216–229. [Google Scholar] [CrossRef]
  6. Yang, Y.; Miao, A. Effect of external vibration on PZT impedance signature. Sensors 2008, 8, 6846–6859. [Google Scholar] [CrossRef] [PubMed]
  7. Chee, C. A Review on the Modeling of Piezoelectric Sensors and Actuators Incorporated in Intelligent Structures. J. Intell. Mater. Syst. Struct. 1998, 9, 3–19. [Google Scholar] [CrossRef]
  8. Li, Z.; Zhang, D.; Wu, K. Cement matrix 2–2 piezoelectric composite—Part 1. Sensory effect. Mater. Struct. 2001, 34, 506–512. [Google Scholar]
  9. Potong, R.; Rianyoi, R.; Ngamjarurojana, A.; Chaipanich, A. Fabrication and performance investigation of 2-2 connectivity lead-free barium zirconate titanate—Portland cement composites. Ceram. Int. 2014, 40, 8723–8728. [Google Scholar] [CrossRef]
  10. Cheng, X.; Xu, D.Y.; Guo, L.L.; Huang, S.F. Influence of thickness parameters on 2-2 cement based piezoelectric composite. Adv. Mater. Res. 2009, 79–82, 31–34. [Google Scholar] [CrossRef]
  11. Zhang, T.T.; Shi, Z.F. Bending behavior of 2-2 multi-layered piezoelectric curved actuators. Smart Mater. Struct. 2007, 16, 634–641. [Google Scholar] [CrossRef]
  12. Zhang, T.; Shi, Z. Two-dimensional exact analysis for piezoelectric curved actuators. J. Micromech. Microeng. 2006, 16, 640. [Google Scholar] [CrossRef]
  13. Shi, Z.; Zhang, T. Static analyses for 2-2 multi-layered piezoelectric curved composites. Int. J. Eng. Sci. 2007, 45, 509–524. [Google Scholar] [CrossRef]
  14. Zhang, T.; Shi, Z. Exact solutions of the piezoelectric transducer under multi loads. Smart Mater. Struct. 2011, 8, 413–431. [Google Scholar] [CrossRef]
  15. Zhang, T.; Shi, Z. Exact analysis of the dynamic properties of a 2-2 cement based piezoelectric transducer. Smart Mater. Struct. 2011, 20, 085017. [Google Scholar] [CrossRef]
  16. Chen, W.; Lv, C.; Bian, Z. Elasticity solution for free vibration of laminated beams. Compos. Struct. 2003, 62, 75–82. [Google Scholar] [CrossRef]
  17. Wang, J.; Shi, Z.; Han, Z. Analytical solution of piezoelectric composite stack transducers. J. Intell. Mater. Syst. Struct. 2013, 24, 1626–1636. [Google Scholar] [CrossRef]
  18. Zhang, T.; Ma, K. Analysis of Dynamic Properties of Piezoelectric Structure under Impact Load. Micromachines 2015, 6, 1577–1587. [Google Scholar] [CrossRef]
  19. Yin, L.; Shen, Y. Strain sensing of composite plates subjected to low velocity impact with distributed piezoelectric sensors: A mixed finite element approach. J. Sound Vib. 1997, 199, 17–31. [Google Scholar] [CrossRef]
  20. Ueda, S. Impact response of a functionally graded piezoelectric plate with a vertical crack. Theor. Appl. Fract. Mech. 2005, 44, 329–342. [Google Scholar] [CrossRef]
  21. Ueda, S. Electromechanical impact of an impermeable parallel crack in a functionally graded piezoelectric strip. Eur. J. Mech. A Solids 2007, 26, 123–136. [Google Scholar] [CrossRef]
  22. Li, Z.; Zhang, D.; Wu, K. Cement-Based 0-3 Piezoelectric Composites. J. Am. Ceram. Soc. 2002, 85, 305–313. [Google Scholar] [CrossRef]
  23. Chung, S.; Chou, C. Electromechanical analysis of an asymmetric piezoelectric/elastic laminate structure: Theory and experiment. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 1999, 46, 441–451. [Google Scholar] [CrossRef] [PubMed]
  24. Rassaian, M.; Beranek, M.W. Quantitative characterization of 96.5Sn3.5Ag and 80Au20Sn optical fiber solder bond joints on silicon micro-optical bench substrates. IEEE Trans. Adv. Packag. 1999, 22, 86–93. [Google Scholar] [CrossRef]
  25. Measurement Specialties, Inc. Piezo Film Sensors Technical Manual; Measurement Specialties, Inc.: Hampton, VA, USA, 1999. [Google Scholar]
Figure 1. Schematic 2-2 cement-based piezoelectric sensor under impact load.
Figure 1. Schematic 2-2 cement-based piezoelectric sensor under impact load.
Sensors 17 02035 g001
Figure 2. Schematics of the transient haversine wave load Q ( t ) .
Figure 2. Schematics of the transient haversine wave load Q ( t ) .
Sensors 17 02035 g002
Figure 3. Theoretical solution of displacement, electric potential and stress: (a) Displacement w c ( l 5 , t ) ; (b) Displacement w p ( l 3 , t ) ; (c) Electric potential ϕ ( l 3 , t ) ; (d) Stress σ p ( l 3 , t ) .
Figure 3. Theoretical solution of displacement, electric potential and stress: (a) Displacement w c ( l 5 , t ) ; (b) Displacement w p ( l 3 , t ) ; (c) Electric potential ϕ ( l 3 , t ) ; (d) Stress σ p ( l 3 , t ) .
Sensors 17 02035 g003
Figure 4. Comparisons of theoretical solutions and the numerical solutions of displacement, electric potential and stress: (a) Displacement w c ( l 5 , t ) ; (b) Displacement w p ( l 3 , t ) ; (c) Electric potential ϕ ( l 3 , t ) ; (d) Stress σ p ( l 3 , t ) .
Figure 4. Comparisons of theoretical solutions and the numerical solutions of displacement, electric potential and stress: (a) Displacement w c ( l 5 , t ) ; (b) Displacement w p ( l 3 , t ) ; (c) Electric potential ϕ ( l 3 , t ) ; (d) Stress σ p ( l 3 , t ) .
Sensors 17 02035 g004aSensors 17 02035 g004b
Figure 5. Displacement w c ( l 5 , t 0 ) , electric potential ϕ ( l 3 , t 0 ) and stress σ c ( l 1 , t 0 ) changing with versus the thickness ratio h 3 / h 2 : (a) Displacement w c ( l 5 , t 0 ) ; (b) Electric potential ϕ ( l 3 , t 0 ) ; (c) Stress σ c ( l 1 , t 0 ) .
Figure 5. Displacement w c ( l 5 , t 0 ) , electric potential ϕ ( l 3 , t 0 ) and stress σ c ( l 1 , t 0 ) changing with versus the thickness ratio h 3 / h 2 : (a) Displacement w c ( l 5 , t 0 ) ; (b) Electric potential ϕ ( l 3 , t 0 ) ; (c) Stress σ c ( l 1 , t 0 ) .
Sensors 17 02035 g005
Figure 6. Displacement w c ( l 5 , t 0 ) , electric potential ϕ ( l 3 , t 0 ) and stress σ c ( l 1 , t 0 ) changing with C 33 p : (a) Displacement w c ( l 5 , t 0 ) ; (b) Electric potential ϕ ( l 3 , t 0 ) ; (c) Stress σ c ( l 1 , t 0 ) .
Figure 6. Displacement w c ( l 5 , t 0 ) , electric potential ϕ ( l 3 , t 0 ) and stress σ c ( l 1 , t 0 ) changing with C 33 p : (a) Displacement w c ( l 5 , t 0 ) ; (b) Electric potential ϕ ( l 3 , t 0 ) ; (c) Stress σ c ( l 1 , t 0 ) .
Sensors 17 02035 g006aSensors 17 02035 g006b
Figure 7. Displacement w c ( l 5 , t 0 ) , electric potential ϕ ( l 3 , t 0 ) and stress σ c ( l 1 , t 0 ) changing with e 33 : (a) Displacement w c ( l 5 , t 0 ) ; (b) Electric potential ϕ ( l 3 , t 0 ) ; (c) Stress σ c ( l 1 , t 0 ) .
Figure 7. Displacement w c ( l 5 , t 0 ) , electric potential ϕ ( l 3 , t 0 ) and stress σ c ( l 1 , t 0 ) changing with e 33 : (a) Displacement w c ( l 5 , t 0 ) ; (b) Electric potential ϕ ( l 3 , t 0 ) ; (c) Stress σ c ( l 1 , t 0 ) .
Sensors 17 02035 g007
Figure 8. Displacement w c ( l 5 , t 0 ) , electric potential ϕ ( l 3 , t 0 ) and stress σ c ( l 1 , t 0 ) changing with   ε 33 S / ε : (a) Displacement w c ( l 5 , t 0 ) ; (b) Electric potential ϕ ( l 3 , t 0 ) ; (c) Stress σ c ( l 1 , t 0 ) .
Figure 8. Displacement w c ( l 5 , t 0 ) , electric potential ϕ ( l 3 , t 0 ) and stress σ c ( l 1 , t 0 ) changing with   ε 33 S / ε : (a) Displacement w c ( l 5 , t 0 ) ; (b) Electric potential ϕ ( l 3 , t 0 ) ; (c) Stress σ c ( l 1 , t 0 ) .
Sensors 17 02035 g008
Figure 9. Influences of piezoelectric materials on the dynamic behaviors of the sensor. (a) Displacement w ( z , t 0 ) ; (b) Electric potential ϕ ( z , t 0 ) ; (c) Stress σ ( z , t 0 ) .
Figure 9. Influences of piezoelectric materials on the dynamic behaviors of the sensor. (a) Displacement w ( z , t 0 ) ; (b) Electric potential ϕ ( z , t 0 ) ; (c) Stress σ ( z , t 0 ) .
Sensors 17 02035 g009
Figure 10. Influence of material parameters on ω : (a) C 33 p ; (b) e 33 ; (c) ε 33 S / ε 0 ; (d)   ρ p .
Figure 10. Influence of material parameters on ω : (a) C 33 p ; (b) e 33 ; (c) ε 33 S / ε 0 ; (d)   ρ p .
Sensors 17 02035 g010aSensors 17 02035 g010b
Table 1. Properties of piezoelectric ceramics, elastic electrode and cement.
Table 1. Properties of piezoelectric ceramics, elastic electrode and cement.
MaterialThicknessDensityElastic Stiffness CoefficientPoisson’s RatioPiezoelectric CoefficientPermittivity Coefficient
Ordinary Portland Cement0.004 m 2500   kg / m 3 2.5 × 10 10   Pa 0.2//
Piezoelectric Ceramics0.010 m 5700   kg / m 3 6.0 × 10 10   Pa / 0.75   C / m 2 52.5 ε 0   1
H62-Brass0.001 m 8430   kg / m 3 10.0 × 10 10   Pa 0.34//
1 ε 0   =   8.85   ×   10   12   F / m is the vacuum dielectric constant.
Table 2. Properties of elastic electrode.
Table 2. Properties of elastic electrode.
MaterialDensityElastic Stiffness Coefficient
H62-Brass 8430   kg / m 3 10 × 10 10   Pa
Aluminum 2536   kg / m 3 7.0 × 10 10   Pa
Gold-tin 80% wt % Au 20 wt % Sn 16,900   kg / m 3 13.73 × 10 10   Pa
Table 3. Properties of piezoelectric materials.
Table 3. Properties of piezoelectric materials.
MaterialDensityElastic Stiffness CoefficientPiezoelectric CoefficientPermittivity Coefficient
PZT-5H 7500   kg / m 3 11.7 × 10 10   Pa 23.3   C / m 2 1470 ε 0   1
PZT-4 7500   kg / m 3 11.5 × 10 10   Pa 15.1   C / m 2 635 ε 0   1
PVDF 1780   kg / m 3 0.25 × 10 10   Pa 0.16   C / m 2 13 ε 0   1
1 ε 0   =   8.85   ×   10   12   F / m is the vacuum dielectric constant.

Share and Cite

MDPI and ACS Style

Zhang, T.; Zhang, K.; Liu, W.; Liao, Y. Impact Mechanical Response of a 2-2 Cement-Based Piezoelectric Sensor Considering the Electrode Layer Effect. Sensors 2017, 17, 2035. https://doi.org/10.3390/s17092035

AMA Style

Zhang T, Zhang K, Liu W, Liao Y. Impact Mechanical Response of a 2-2 Cement-Based Piezoelectric Sensor Considering the Electrode Layer Effect. Sensors. 2017; 17(9):2035. https://doi.org/10.3390/s17092035

Chicago/Turabian Style

Zhang, Taotao, Keping Zhang, Wende Liu, and Yangchao Liao. 2017. "Impact Mechanical Response of a 2-2 Cement-Based Piezoelectric Sensor Considering the Electrode Layer Effect" Sensors 17, no. 9: 2035. https://doi.org/10.3390/s17092035

APA Style

Zhang, T., Zhang, K., Liu, W., & Liao, Y. (2017). Impact Mechanical Response of a 2-2 Cement-Based Piezoelectric Sensor Considering the Electrode Layer Effect. Sensors, 17(9), 2035. https://doi.org/10.3390/s17092035

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