Next Article in Journal
Hydrophobin-Based Surface Engineering for Sensitive and Robust Quantification of Yeast Pheromones
Previous Article in Journal
Search Techniques for the Web of Things: A Taxonomy and Survey
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Simulation of Output Response of PVDF Sensor Attached on a Cantilever Beam Subjected to Impact Loading

Department of Civil Engineering, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8552, Japan
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(5), 601; https://doi.org/10.3390/s16050601
Submission received: 10 February 2016 / Revised: 1 April 2016 / Accepted: 20 April 2016 / Published: 27 April 2016
(This article belongs to the Section Physical Sensors)

Abstract

:
Polyvinylidene Flouride (PVDF) is a film-type polymer that has been used as sensors and actuators in various applications due to its mechanical toughness, flexibility, and low density. A PVDF sensor typically covers an area of the host structure over which mechanical stress/strain is averaged and converted to electrical energy. This study investigates the fundamental “stress-averaging” mechanism for dynamic strain sensing in the in-plane mode. A numerical simulation was conducted to simulate the “stress-averaging” mechanism of a PVDF sensor attached on a cantilever beam subjected to an impact loading, taking into account the contribution of piezoelectricity, the cantilever beam’s modal properties, and electronic signal conditioning. Impact tests and FEM analysis were also carried out to verify the numerical simulation results. The results of impact tests indicate the excellent capability of the attached PVDF sensor in capturing the fundamental natural frequencies of the cantilever beam. There is a good agreement between the PVDF sensor’s output voltage predicted by the numerical simulation and that obtained in the impact tests. Parametric studies were conducted to investigate the effects of sensor size and sensor position and it is shown that a larger sensor tends to generate higher output voltage than a smaller one at the same location. However, the effect of sensor location seems to be more significant for larger sensors due to the cancelling problem. Overall, PVDF sensors exhibit excellent sensing capability for in-plane dynamic strain induced by impact loading.

Graphical Abstract

1. Introduction

Polyvinylidene Fluoride (PVDF) is a thin film-type polymer that is mechanically tough, flexible, and low density. The piezoelectric effect in elongated and polarized films of polymers, particularly of PVDF, was discovered by Kawai [1] in 1969. Since then, the fundamental properties of PVDF have been extensively investigated. The semi-crystalline molecular structure of PVDF consists of long chain molecules with a repeating CF2CH2 unit. Application of heat, electrical fields, and pressure can interconvert the four different forms of PVDF’s crystalline domains [2,3,4]. The molecular dipoles in the crystalline parts are oriented by thermal poling or corona poling, thus resulting in a permanent polarization. In the β-phase, PVDF exhibits piezoelectric effect which means mechanical energy can be converted to electrical energy and vice versa. Therefore, PVDF has been frequently used to manufacture sensors and actuators for a number of practical applications such as shock impact and pressure sensors [5,6,7,8], biomedical [9,10,11], acoustic [12,13,14,15], tactile sensors [16,17,18], active vibration control [19], and structural health monitoring of civil and aerospace structures [20,21,22].
Many previous studies have focused on the use of PVDF sensors in the out-of-plane (3-3) mode in which mechanical stress is induced in the thickness (poled) direction. However, the in-plane sensing mode of PVDF sensor has also been frequently employed for dynamic impact sensing. Lee and O’Sullivan [23] developed a uniaxial strain rate gage that measures only strain rate along a specified direction by combining the effective surface electrodes, appropriate skew angle, and the correct polarization profile. Wang and Wang [24] presented a theoretical approach for feasibility analysis of the application of PVDF sensors to cantilever beam modal testing. Sirohi and Chopra [25] investigated the behavior of piezoelectric elements including piezoceramic (PZT) and piezofilm (PVDF) as dynamic strain sensors of which superior performance compared to conventional strain gages in terms of sensitivity and signal-to-noise ratio was demonstrated. Correction factors to account for transverse strain and shear-lag effect due to bonding layer were analytically derived and experimentally validated by the same authors. Ma et al. [26] investigated the effects of a PVDF sensor’s area and the use of a charge amplifier on the measurement capability of a PVDF sensor attached to a cantilever beam that is subjected to impact loading. PVDF sensors proved to be capable of capturing most of the resonant frequencies from transient responses, and their sensitivity was demonstrated to be better than that of conventional strain gages. Kotian et al. [27] presented an analytical investigation on the effects of stress-averaging for both in-plane sinusoidal stress waves and in-plane impact-induced stresses. It was concluded that the error induced by stress averaging becomes more significant as sensor length increases, density of structure’s material increases, and magnitude of input stress increases, although the error induced by stress averaging is minimal for most practical applications. Furthermore, only very high frequencies (in the order of kHz) can cause a significant reduction in a PVDF sensor’s output voltage due to stress-averaging.
In civil engineering structures, dynamic strain is one of the most fundamental measures. To capture dynamic strain, conventional foil strain gages have been frequently used. PVDF sensor with its superior signal-to-noise ratio can be an alternative for strain gages in several practical applications. Since strain gage size is usually minute compared to that of host structures, dynamic strain output can be considered strain at a point on the host structural member. On the other hand, PVDF sensors are much larger in size compared to strain gages. An electrode of a PVDF sensor typically covers a relatively large area of the host structural member over which dynamic strain may vary significantly depending on the relative size of the PVDF sensor compared to that of the host member and strain gradient. To understand the mechanism of strain averaging, both extreme scenarios including minor and large strain variations should be investigated. In case of minor strain variation, a strain averaging mechanism was thoroughly discussed in Kotian et al. [27] using tensile specimens subjected to sinusoidal excitation. Ma et al. [26] indicated that the use of a charge amplifier is indispensable for improvement of a PVDF sensor’s capability in capturing low-frequency vibration modes of a steel cantilever beam subjected to impact loading. However, the conversion mechanism of dynamic strain to output voltage, i.e., strain averaging mechanism, was not explicitly explained in the same study. Moreover, the effect of sensor size was discussed by comparing sensor pairs at a fixed point for the purpose of illustrating charge amplifier’s indispensability for low-frequency measurement.
This study investigates the strain averaging mechanism of a PVDF sensor attached to an area of the host structural member where large strain variation occurs. A numerical simulation based on the governing equations of piezoelectricity, classical beam theory, and electronic signal conditioning to predict the output voltage of a PVDF sensor would be proposed. To verify the numerical simulation results, experimental impact tests and FEM analysis would be conducted using a steel cantilever beam subjected to impact forces since a cantilever beam can be regarded as the most fundamental and flexible structure that has high strain gradient along its length. The proposed numerical simulation would provide an insight into the mechanism of dynamic impact strain sensing by a PVDF sensor attached to a surface with large strain variation. Furthermore, the effects of sensor size and position on the output voltage would be investigated in parametric studies using the proposed numerical simulation.

2. Theoretical Background

2.1. Governing Equation for Piezoelectricity

Consider the constitutive governing equation that is reduced from the tensor expression for the piezoelectric effect induced by one-dimensional mechanical deformation [28]:
D 3 = d 31 T 1 + ε 33 T E 3  
where D3 is the electrical displacement component; d31 is the piezoelectric constant; T1 is the axial stress component, ε 33 T is the permittivity component at constant stress, and E3 is the electrical field component. Subscripts 1 and 3 indicate the longitudinal direction and the poling direction, respectively (Figure 1). In case the external electrical field is absent (E3 = 0), the relation is further reduced to:
D 3 = d 31 T 1 = d 31 E p S 1
where S 1 is the bending strain and Ep is the Young’s modulus of the PVDF layer.
The bending strain is expressed by the bending moment/curvature differential equation assuming small deflection and rotation [29]:
S 1 ( x ,  t ) = y p 0 2 y ( x ,  t ) x 2
where yp0 is the distance from the neutral axis of the cross-section to the center of the PVDF layer; y(x, t) is the transverse displacement of the cantilever beam, which can be represented by a convergent series of the eigenfunctions as
y ( x ,  t ) = r = 1 ϕ r ( x ) η r ( t )
where ϕ r ( x ) is the mass-normalized eigenfunction; η r ( t ) is the modal coordinate of the cantilever beam for the r-th vibration mode. If the cantilever beam is assumed to be proportionally damped, the eigenfunctions denoted by ϕ r ( x ) are the mass-normalized eigenfunctions of the undamped free vibration [30,31].
ϕ r ( x ) = 1 mL [ cosh λ r L x cos λ r L x σ r ( sinh λ r L x sin λ r L x ) ]
where m is the mass per unit length of the cantilever beam; λr is the dimensionless frequency number for each mode obtained from the following characteristic equation:
1 + cos λ cosh λ = 0
σ r = s i n h λ r s i n λ r c o s h λ r + c o s λ r
The mass-normalized eigenfunctions satisfy the following orthogonality conditions:
x = 0 L m ϕ s ( x ) ϕ r ( x ) dx = δ rs ; x = 0 L EI ϕ s ( x ) d 4 ϕ r ( x ) dx 4 dx = ω r 2 δ rs
where E is the Young’s modulus of the cantilever beam; δrs is the Kronecker delta, δrs = 1 for s = r and δrs = 0 for s ≠ r; and ωr is the undamped natural frequency of the r-th mode:
ω r = λ r 2 E I m L 4

2.2. Governing Mechanical Equation

The governing equation of motion can be written as [31]:
EI 4 y ( x ,  t ) x 4 + c s I 5 y ( x ,  t ) x 4 t + c a y ( x ,  t ) t + m 2 y ( x ,  t ) t 2 = p ( t )
where y(x, t) is the transverse displacement of the cantilever beam; cs is the equivalent coefficient of strain rate damping; and I is the equivalent area moment of inertia of the cross-section. In terms of damping, csI represents the equivalent damping term of the cross section due to structural viscoelasticity while ca is the viscous air damping coefficient. Both of the above-mentioned damping mechanisms satisfy the proportional damping criterion [30,31]; p(t) is the time history of external excitation. For an impact, p(t) can be represented as:
p ( t ) = F δ ( x x F ) δ ( t τ )
where F is the magnitude of the impact force; xF is the location of impact force; and δ is the direct delta function.
Substitution of y(x, t) in Equation (4) into the governing equation of motion (Equation (10)) and using the orthogonality conditions (Equation (8)), the modal response of the cantilever beam can be obtained from the following electromechanically coupled ordinary differential equation:
d 2 η r ( t ) dt 2 + 2 ζ r ω r d η r ( t ) dt + ω r 2 η r ( t ) = N r ( t )
where ζ r is the mechanical damping ratio including both effects of strain rate damping and viscous air damping [32],
ζ r = c s I ω r 2 EI + c a 2 m ω r
The assumption of proportional damping was discussed in [32]. In this assumption, once the proportionality constants cs and ca are identified using the modal properties (i.e., natural frequencies and damping ratios) of two vibration modes, the other mode’s damping ratio is not arbitrary but mathematically derived from Equation (13). Identified damping ratios of the vibration modes of interest might also be used directly without obtaining the cs and ca values since the resulting electromechanical expressions only need the ζ r values.
The modal mechanical forcing function, Nr(t) , can be expressed as,
N r ( t ) = x = 0 L ϕ r ( x ) p ( t ) dx
The solution of Equation (12) can be expressed using the unit impulse response function in the form of Duhamel integration as
η r ( t ) = 1 ω rd τ = 0 t N r ( τ ) e - ζ r ω r ( t - τ ) sin ω rd ( t τ ) d τ
where ωrd is the damped natural frequency of the r-th mode,
ω rd = ω r 1 ζ r 2

2.3. Governing Electrical Equation

There are two equally valid equivalent electrical models of the piezofilm element—one is a voltage source in series with a capacitance that is equal to the capacitance of the sensor, the other a charge generator in parallel with a capacitance [33]. The latter is used in this study.
Capacitance of a PVDF sensor is expressed as:
C p = ε A t
where ε is the permittivity, which can also be expressed in the form of ε = ε r ε 0 where ε r is the relative permittivity (about 12 for PVDF) and ε 0 is the permittivity of free space (constant, 8.854 × 10−12 F/m); A is the active area of the film’s electrodes; and t is the film thickness.
The signal conditioning circuit is shown in Figure 2. The advantages of connecting a PVDF sensor to a charge amplifier have been emphasized in [25]. First, the charge generated by the sensor is transferred onto the feedback capacitance CF. The voltage output of a charge amplifier is proportional to the input strain, irrespective of the sensor’s capacitance. The gain is controlled by the feedback capacitance of the charge amplifier. Second, the value of time constant, defined as RFCF, can be selected in order to navigate the loading effect in order to obtain the desired dynamic frequency range. Ma et al. [26] also indicated that a charge amplifier is indispensable for a PVDF sensor, especially a small-sized one, to improve the low-frequency responses of a PVDF sensor.
The governing equation of a charge amplifier’s output voltage can be found in [19]:
V 0 = j ω A V s C p [ j ω ( ( A + 1 ) C F + C p + C c ) + ( 1 R a + ( A + 1 ) 1 R F ) ]
where ω is the angular frequency (rad/s); Vs is the voltage generated by the PVDF sensor; Cp is the equivalent capacitance of the PVDF sensor; Ra is the output impedance of the PVDF sensor; Cc is the equivalent capacitance of the electric wire; A is the gain of the charge amplifier; and CF and RF are the feedback capacitance and impedance of the charge amplifier, respectively.
When the magnitude of gain is large enough, the output voltage equation can be reduced to:
V 0 = j ω V s C p j ω C F + 1 R F
In the high frequency region, the output voltage can be further reduced to:
V 0 = V s C p C F = q C F
Meanwhile, in the low-frequency region, the amplitude of output voltage is expressed as:
| V 0 | = ω q 1 R F 2 + ω 2 C F 2
where q is the electrical charge accumulated on the PVDF sensor’s electrodes, which is given as
q = A d 31 E p S 1 n d A = b p x p 1 x p 2 d 31 E p ( y p 0 2 y ( x ,   t ) x 2 ) dx = y p 0 d 31 E p b p η r ( t ) r = 1 d ϕ r ( x ) dx | x p 1 x p 2
The cut-off frequency is written as:
f c = 1 2 π R F C F

3. Numerical Simulation of PVDF Sensor’s Response to Impact Force

Numerical simulation was conducted using MATLAB. The sampling rate was chosen as 10 kHz and one second of the PVDF sensor’s output voltage response was simulated. The basic input parameters for simulation are shown in Table 1. Figure 3 illustrates the simulation flow. Three vibration modes were included for impact at locations 1 and 2 and five modes for location 3. Natural frequencies of the cantilever beam were determined using Equations (9) and (16) for flexural modes. Damping ratios obtained by the half-power method (Table 3) were employed to create the unit impulse-response function. For mode shape simulation, the length increment of 1 mm was chosen. Impact time-histories at three impact locations illustrated in Figure 4 were employed as input impact loading in the simulation process (Figure 5). The unit impulse-response function, modal response, and mode shape slope difference are determined for each vibration mode. The electrical charge is then summed from the contribution of all included vibration modes. Predicted output voltage is finally determined from the electrical charge determined by Equation (20).

4. Impact Tests

To obtain the actual impact histories that would be input for numerical simulation and to verify the numerical simulation results, impact tests were conducted by applying impact loadings on a steel cantilever beam on which a PVDF sensor was attached. The piezofilm lab amplifier [34] operating in charge mode was used for signal conditioning. Feedback capacitance of 100 nF was selected. The lower cut-off frequency and higher cut-off frequency were selected as 0.1 Hz and 100 kHz, respectively. The NI USB-4431 enabling four 24-bit simultaneous analog inputs was used for data acquisition. 10 kS/s was chosen for sampling rate and 1 kS/s for data windowing for real-time visualization of the PVDF sensor’s output response during the impact tests. One DT4-028 K/L PVDF sensor [35] was attached to the steel cantilever beam by a commercial double-coated adhesive tape. One edge of the PVDF sensor distances 30 mm from the fixed end of the cantilever beam (Figure 4). The impact hammer PCB 086C03 manufactured by PCB Piezotronics with a mounted plastic tip was used to apply impact force on the cantilever beam at the three impact locations illustrated in Figure 4. During the impact tests, impact signal was transmitted to the ICP sensor signal conditioner provided by PCB Piezotronics. A unit gain was set on the ICP signal conditioner. To verify the numerical simulation results, strain gages were also attached to the opposite side of the cantilever beam to obtain “true” strain (Figure 4). Typical impact force time-histories are shown in Figure 5. Impact force histories at locations 1 and 2 contain two peaks which are both due to the physical contact between the impact hammer and the cantilever beam. Particularly, the first peak represents the first time when the cantilever beam is collided by the hammer tip. The energy absorbed by the cantilever beam is then transformed to displacement of the cantilever beam’s free end, which mainly excites the first vibration mode. Right after the first collision, the hammer tip continues to move in the same direction as that of the cantilever beam’s free end due to the effect of inertia induced by the hammer’s weight and velocity. Approximately one to two milliseconds after the first collision, the hammer hits the tip of cantilever beam the second time, which causes the second impact. However, for impact at location 3 which is close to the fixed end, impact energy absorbed by the cantilever beam is transformed mainly into excitation of higher modes rather than displacement at the point of collision. Therefore, only one peak appears in the time history of impact at location 3. However, time history of impact rather than impact force magnitude (i.e., peak) was employed in the numerical simulation, which means impact force energy is taken into account in the prediction of the PVDF sensor’s output voltage in this study.

5. Verification of Numerical Simulation Results by FEM

In order to verify the simulation results, FE model of the cantilever beam subjected to impact loading was analyzed using the commercial software ABAQUS. The purposes of FEM analysis were: (1) To acquire the knowledge of the cantilever beam’s vibration modes, as well as mode shapes and natural frequencies that would be used for the numerical simulation; (2) to obtain dynamic strain for verification of the numerical simulation results; and (3) to convert output impact histories of the ICP signal conditioner into force unit. Table 2 indicates that the natural frequencies of the cantilever beam obtained by FEM analysis, classical beam theory, and modal impact test exhibit a good agreement. Most of the flexural modes can be accurately identified by the PVDF sensor, which indicates the capability of PVDF sensors in identifying natural frequencies of the host structure’s fundamental vibration modes. The first torsional mode at 341.41 Hz did not appear in the frequency response functions (Figure 6), which may indicate the dominance of flexural modes for all three impact locations. If torsional modes are excited due to an impact location that deviates from the central line of the cantilever beam, PVDF sensors can also capture such torsional modes. The superiority of PVDF sensors compared to conventional foil strain gages in capturing torsional modes was indicated in Ma et al. [26]. For the second purpose, FEM modal dynamic analysis was conducted using the actual recorded impact that had been previously generated by the impact hammer during the impact tests. Direct damping ratios (Table 3) determined by the half-power point method [36] using the frequencies response functions shown in Figure 6 were used for the first five flexural modes of the cantilever beam in the FEM analysis. To obtain the theoretically simulated strain shown in Figure 7, the transverse displacement in Equation (4) was determined from the mass-normalized eigenfunction written in Equation (5) and the modal coordinate (i.e., the solution of Equation (12)) expressed in Equation (15). Dynamic strain was then calculated in Equation (3) by taking the second derivative of the known transverse displacement. There appears to be a good agreement between the strain values obtained by FEM analysis, those measured in the impact tests and those predicted by the numerical simulation for impacts at locations 1 and 2 (Figure 7a,b). For impact at location 3, only the strain values determined by FEM analysis and those predicted by simulation were compared since the strain amplitude is only as low as the noise level in strain gages. A good agreement was also observed for impact location 3 (Figure 7c).

6. Results of Numerical Simulation

There appears to be a good agreement between the actual output voltage of the PVDF sensor obtained in the impact test and that predicted by the numerical simulation. Even for impact at location 3 for which the amplitude of output signal is much lower than that for locations 1 and 2 (Figure 8a,b), a reasonable prediction by simulation was also obtained (Figure 8c,d).
Parametric studies were conducted using sensor size and sensor position as input parameters for the numerical simulation described in the previous section in order to investigate the effect of sensor size and sensor position on the output voltage of a PVDF sensor. The input impact forces, damping ratios and all remaining input parameters were chosen similar to those used for the prediction of output voltage shown in Figure 8. Five different sensor lengths including 20, 40, 80, 160, and 240 mm were employed in the simulation. The width of simulated sensors is 20 mm. For the investigation on the effect of sensor position, a midpoint position increment of 10 mm was used. Figure 9 shows the simulated output voltage for different sensor sizes and midpoint positions. For impact at location 1, there appears to be one position at about 260 mm from the fixed end where the output voltage reaches minimum regardless of sensor size. PVDF sensors located closer to the fixed end tend to generate higher peak voltage, which may be attributed to the fact that in the first mode of vibration bending moment increases towards the fixed end. The fact that the minimum voltage occurs at about 260 mm might be due to the problem of cancelling of mode shape slope calculated at the two opposite edges of PVDF sensor (Equation (22)). For impact at locations 2 and 3, the minimum output voltage occurs at multiple positions, which may be due to the cancellation problem at multiple nodal points as a result of higher number of excited vibration modes. The first five flexural mode shapes as well as the corresponding nodal points for each mode shape were illustrated in Figure 10. For all three impact locations, a larger sensor tends to generate higher output voltage and vice versa. For impact location 3, a sensor longer than 240 mm appears to result in no increase in output voltage. Furthermore, larger sensors appear to result in larger variation in the output voltage, which means that the effect of sensor location tends to be more significant for sensor with larger sizes.

7. Conclusions

This study investigates the “stress-averaging” mechanism of the in-plane sensing mode of the PVDF sensor for dynamic strain sensing. A numerical simulation was conducted to predict the output voltage response of a PVDF sensor attached to a steel cantilever beam subjected to impact loading based on the fundamental knowledge of piezoelectricity, classical beam theory, and signal conditioning. FEM analysis and impact tests were also conducted to verify the simulation results. The results of impact test indicate the excellent capability of PVDF sensors in capturing the fundamental natural frequencies of the cantilever beam. The PVDF sensor’s output voltage could be reasonably predicted by the numerical simulation. Parametric studies on the effects of sensor size and sensor position indicate that a larger sensor tends to generate higher output voltage and vice versa. Furthermore, the effect of sensor position seems to be more significant for larger sensors. However, when a large number of modes are excited, e.g., impact at location 3, increasing sensor size may not always result in increased output voltage. The sensor size beyond which increasing sensor size does not result in increased output voltage may depend on the number of dominant vibration modes which in turn depends on impact location. The results indicate that, in order to maximize the signal-to-noise ratio for actual measurements in small scale structures (i.e., the size of an attached PVDF sensor is relatively large with respect to the size of the host structural member), vibration mode shapes of the structural member should be considered to choose the optimal sensor location that helps avoid the cancelling problem. Overall, PVDF sensors exhibit an excellent in-plane sensing capability for dynamic strain induced by impact loading.

Author Contributions

Both authors conceived and designed the experiments; the first author performed the experiment, analyzed the data, proposed the numerical simulation, and wrote the paper. The second author dedicated his valuable suggestions and comments on the interpretation of results and on the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kawai, H. The piezoelectricity of poly (vinylidene fluoride). Jpn. J. Appl. Phys. 1969, 8, 975–976. [Google Scholar] [CrossRef]
  2. Lando, J.B.; Olf, H.G.; Peterlin, A. Nuclear magnetic resonance and x-ray determination of the structure of poly (vinylidene fluoride). J. Polym. Sci. 1966, 4, 941–951. [Google Scholar] [CrossRef]
  3. Sessler, G.M. Piezoelectricity in polyvinylidenefluoride. J. Acoust. Soc. Am. 1981, 70, 1596–1608. [Google Scholar] [CrossRef]
  4. Darasawa, N.; Goddard, W.A. Force fields, structures, and properties of poly (vinylidene fluoride) crystals. Macromolecules 1992, 25, 7268–7281. [Google Scholar]
  5. Lee, L.M.; Graham, R.A.; Bauer, F.; Reed, R.P. Standardized Bauer PVDF piezoelectric polymer shock gauge. J. De Physique 1988, 9, 651–657. [Google Scholar] [CrossRef]
  6. Dutta, P.K.; Kalafut, J. Evaluation of PVDF Piezopolymer for Use as a Shock Gauge; Special Report; U.S. Army Corps of Engineers: Hanover, NH, USA, 1990; pp. 90–23. [Google Scholar]
  7. Bauer, F. Ferroelectric PVDF polymer for high pressure and shock compression sensors. In Proceedings of the 11th International Symposium on Electrets, Melbourne, Australia, 1–3 October 2002; pp. 219–222.
  8. Shirinov, A.V.; Schomburg, W.K. Pressure sensor from a PVDF film. Sens. Actuators A 2008, 142, 48–55. [Google Scholar] [CrossRef]
  9. Wang, F.; Tanaka, M.; Chonan, S. Development of a wearable mental stress evaluation system using PVDF sensor. J. Adv. Sci. 2006, 18, 170–173. [Google Scholar] [CrossRef]
  10. Kim, K.J.; Chang, Y.M.; Yoon, S.; Kim, H.J. A Novel Piezoelectric PVDF Film-Based Physiological Sensing Belt for a Complementary Respiration and Heartbeat Monitoring System. Integr. Ferroelectr. 2009, 107, 53–68. [Google Scholar] [CrossRef]
  11. Chiu, Y.; Lin, W.; Wang, H.; Huang, S.; Wu, H. Development of a piezoelectric polyvinylidene fluoride (PVDF) polymer-based sensor patch for simultaneous heartbeat and respiration monitoring. Sens. Actuators A 2013, 189, 328–334. [Google Scholar] [CrossRef]
  12. Woodward, B.; Chandra, R.C. Underwater acoustic measurements on polyvinylidene fluoride transducers. Electrocomponent Sci. Technol. 1978, 5, 149–157. [Google Scholar] [CrossRef]
  13. Xu, J. Microphone Based on Polyvinylidene Fluoride (PVDF) Micro-Pillars and Patterned Electrodes. Ph.D. Thesis, Ohio State University, Columbus, OH, USA, 2010. [Google Scholar]
  14. Headings, L.M.; Kotian, K.; Dapino, M.J. Speed of sound measurement in solids using polyvinylidene fluoride (PVDF) sensors. In Proceedings of the ASME 2013 Conference on Smart Materials, Adaptive Structures and Intelligent Systems, SMASIS2013, Snowbird, UT, USA, 16–18 September 2013.
  15. Xu, J.; Headings, L.M.; Dapino, M.J. High sensitivity polyvinylidene fluoride microphone based on area ratio amplification and minimal capacitance. IEEE Sens. J. 2015, 15, 2839–2847. [Google Scholar] [CrossRef]
  16. Lee, M.H.; Nicholls, H.R. Tactile sensing for mechatronics—A state of the art survey. Mechatronics 1999, 9, 1–31. [Google Scholar] [CrossRef]
  17. Dargahi, J. A piezoelectric tactile sensor with three sensing elements for robotic, endoscopic and prosthetic applications. Sens. Actuators A 2000, 80, 23–30. [Google Scholar] [CrossRef]
  18. Seminara, L.; Capurro, M.; Cirill, P.; Cannata, G.; Valle, M. Electromechanical characterization of piezoelectric PVDF polymer films for tactile sensors in robotics applications. Sens. Actuators A 2011, 169, 49–58. [Google Scholar] [CrossRef]
  19. Ma, C.; Chuang, K.; Pan, S. Polyvinylidene Fluoride Film Sensors in Collocated Feedback Structural Control: Application for Suppressing Impact-Induced Disturbances. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2011, 12, 2539–2554. [Google Scholar] [CrossRef] [PubMed]
  20. Matsumoto, E.; Biwa, S.; Katsumi, K.; Omoto, Y.; Iguchi, K.; Shibata, T. Surface strain sensing with polymer piezoelectric film. NDT&E Int. 2004, 37, 57–64. [Google Scholar]
  21. Kurata, M.; Li, X.; Fujita, K.; He, L.; Yamaguchi, M. PVDF Piezo Film as Dynamic Strain Sensor for Local Damage Detection of Steel Frame Buildings. In Proceedings of the SPIE Sensors and Smart Structures Technologies for Civil, Mechanical, and Aerospace Systems, San Diego, CA, USA, 10 March 2013.
  22. Bloomfield, P.E. Piezoelectric Polymer Transducers for Detection of Structural Defects in Aircraft; Final Report; Naval Air Development Center: Warminster, PA, USA, 1977. [Google Scholar]
  23. Lee, C.K.; O’Sullivan, T.C. Piezoelectric strain rate gages. J. Acoust. Soc. Am. 1991, 2, 945–953. [Google Scholar] [CrossRef]
  24. Wang, B.; Wang, C. Feasibility analysis of using piezoceramic transducers for cantilever beam modal testing. Smart Mater. Struct. 1997, 6, 106–116. [Google Scholar] [CrossRef]
  25. Sirohi, J.; Chopra, I. Fundamental Understanding of Piezoelectric Strain Sensors. J. Intell. Mater. Syst. Struct 2000, 11, 246–257. [Google Scholar] [CrossRef]
  26. Ma, C.; Huang, Y.; Pan, S. Investigation of the Transient Behavior of a Cantilever Beam Using PVDF Sensors. Sensors 2012, 12, 2088–2117. [Google Scholar] [CrossRef] [PubMed]
  27. Kotian, K.; Headings, L.M.; Dapino, M.J. Stress Averaging in PVDF Sensors for in-Plane Sinusoidal and Impact-Induced Stresses. IEEE Sens. J. 2013, 11, 4444–4451. [Google Scholar] [CrossRef]
  28. Standards Committee of IEEE Ultrasonics. Ferroelectrics, and Frequency Control Society; IEEE Standard on Piezoelectricity: New York, NY, USA, 1987. [Google Scholar]
  29. Beer, F.P.; Johnston, E.R.; Dewolf, J.T. Mechanics of materials, 5th ed.; McGraw-Hill Inc.: New York, NY, USA, 2006. [Google Scholar]
  30. Caughey, T.K.; O’Kelly, M.E.J. Classical Normal Modes in Damped Linear Dynamic Systems. ASME J. Appl. Mech. 1965, 32, 583–588. [Google Scholar] [CrossRef]
  31. Erturk, A.; Inman, D.J. A Distributed Parameter Electromechanical Model for Cantilevered Piezoelectric Energy Harvesters. J. Vib. Acoust. 2008, 130, 1–15. [Google Scholar] [CrossRef]
  32. Erturk, A. Electromechanical Modeling of Piezoelectric Energy Harvesters. Ph.D. Thesis, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA, 2008. [Google Scholar]
  33. Measurement Specialties Inc. Piezo Film Sensors Technical Manual. Available online: http://www.meas-spec.com/downloads/Piezo_Technical_Manual.pdf (accessed on 7 July 2015).
  34. Measurement Specialties Inc. Piezo Film Lab Amplifier. Available online: http://www.meas-spec.com/downloads/Piezo_Film_Lab_Amplifier.pdf (accessed on 25 March 2015).
  35. Measurement Specialties Inc. DT Series Elements with Lead Attachment. Available online: http://www.meas-spec.com/downloads/DT_Series_LeadAttach.pdf (accessed on 5 March 2015).
  36. Butterworth, J.; Lee, J.H.; Davidson, B. Experimental determination of modal damping from full-scale testing. In Proceedings of the 13th World Conference on Earthquake Engineering, Vancouver, BC, Canada, 1–6 August 2004.
Figure 1. Schematic configuration of PVDF sensor.
Figure 1. Schematic configuration of PVDF sensor.
Sensors 16 00601 g001
Figure 2. Electrically equivalent model of a PVDF sensor connected to a charge amplifier.
Figure 2. Electrically equivalent model of a PVDF sensor connected to a charge amplifier.
Sensors 16 00601 g002
Figure 3. Flow of numerical simulation.
Figure 3. Flow of numerical simulation.
Sensors 16 00601 g003
Figure 4. Configuration of cantilever beam with attached sensor and impact locations (Unit: mm).
Figure 4. Configuration of cantilever beam with attached sensor and impact locations (Unit: mm).
Sensors 16 00601 g004
Figure 5. Time-histories of impact loadings at three locations.
Figure 5. Time-histories of impact loadings at three locations.
Sensors 16 00601 g005
Figure 6. Frequency response functions for three impact locations.
Figure 6. Frequency response functions for three impact locations.
Sensors 16 00601 g006
Figure 7. Strain measured at 70 mm from the fixed end for three impact locations. (a,b) Impact at location 1; (c,d) Impact at location 2; (e,f) Impact at location 3.
Figure 7. Strain measured at 70 mm from the fixed end for three impact locations. (a,b) Impact at location 1; (c,d) Impact at location 2; (e,f) Impact at location 3.
Sensors 16 00601 g007
Figure 8. Simulation results for three impact locations. (a) Impact at location 1; (b) Impact at location 2; (c,d) Impact at location 3. The sensor size is 155 × 19 × 0.028 (mm). One sensor edge is placed 30 mm from the fixed end of the cantilever beam so that the lead wire goes towards the fixed end.
Figure 8. Simulation results for three impact locations. (a) Impact at location 1; (b) Impact at location 2; (c,d) Impact at location 3. The sensor size is 155 × 19 × 0.028 (mm). One sensor edge is placed 30 mm from the fixed end of the cantilever beam so that the lead wire goes towards the fixed end.
Sensors 16 00601 g008
Figure 9. Simulated output voltage of PVDF sensor with different sensor lengths and positions. (a) Impact at location 1; (b) Impact at location 2; (c) Impact at location 3.
Figure 9. Simulated output voltage of PVDF sensor with different sensor lengths and positions. (a) Impact at location 1; (b) Impact at location 2; (c) Impact at location 3.
Sensors 16 00601 g009
Figure 10. Mode shapes of the first five flexural vibration modes.
Figure 10. Mode shapes of the first five flexural vibration modes.
Sensors 16 00601 g010
Table 1. Input parameters for simulation.
Table 1. Input parameters for simulation.
ItemSymbolParameterValue
SensorLpLength155 mm
bpWidth19 mm
tpThickness0.028 mm
d31Piezoelectric constant23 × 10−12 (C/N)
EpYoung’s modulus4 GPa
CpCapacitance11 nF
Charge AmplifierfcCut-off frequency0.1 Hz
CFFeedback capacitance100 nF
Cantilever beamρDensity7860 kg/m3
EYoung’s modulus200 GPa
Table 2. Natural frequencies of the cantilever beam (Hz).
Table 2. Natural frequencies of the cantilever beam (Hz).
Mode No.Mode TypeFEMTheoryModal TestMaximum Difference (%)
Location 1Location 2Location 3
11st Flexural12.13512.07112.2912.312.322.02
22nd Flexural76.02975.65177.0677.0877.151.94
31st Lateral120.44n/a119.5n/an/a0.79
43rd Flexural212.89211.85215.2215.4215.51.69
51st Torsional341.41n/an/an/an/an/a
64th Flexural417.28415.14422.4422.5422.91.83
75th Flexural690.06686.19698.3698.3698.91.81
82nd Lateral739.77n/an/an/an/an/a
92nd Torsional1026.7n/an/an/an/an/a
106th Flexural1031.31025.11044104410441.81
Table 3. Idenfified damping ratio by half-power method.
Table 3. Idenfified damping ratio by half-power method.
1st Mode2nd Mode3rd Mode4th Mode5th Mode
Impact at 10.0111590.0025660.005434--
Impact at 20.0093910.0026470.006121--
Impact at 30.0135420.0023460.0059200.0013620.002831

Share and Cite

MDPI and ACS Style

Dung, C.V.; Sasaki, E. Numerical Simulation of Output Response of PVDF Sensor Attached on a Cantilever Beam Subjected to Impact Loading. Sensors 2016, 16, 601. https://doi.org/10.3390/s16050601

AMA Style

Dung CV, Sasaki E. Numerical Simulation of Output Response of PVDF Sensor Attached on a Cantilever Beam Subjected to Impact Loading. Sensors. 2016; 16(5):601. https://doi.org/10.3390/s16050601

Chicago/Turabian Style

Dung, Cao Vu, and Eiichi Sasaki. 2016. "Numerical Simulation of Output Response of PVDF Sensor Attached on a Cantilever Beam Subjected to Impact Loading" Sensors 16, no. 5: 601. https://doi.org/10.3390/s16050601

APA Style

Dung, C. V., & Sasaki, E. (2016). Numerical Simulation of Output Response of PVDF Sensor Attached on a Cantilever Beam Subjected to Impact Loading. Sensors, 16(5), 601. https://doi.org/10.3390/s16050601

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