임성현
(Sung-Hyun Lim)
*iD
이수호
(Su-Ho Lee)
*iD
나준영
(Jun-Young Na)
*iD
황재은
(Jae-Eun Hwang)
*iD
유재성
(Jae-Sung Yoo)
*iD
박혜리
(Herie Park)
†iD
-
(Dept. of Electrical Engineering, Dong-A University, Republic of Korea.)
Copyright © The Korean Institute of Electrical Engineers
Key Words
COMSOL Multiphysics, Piezoelectric energy harvester, Thickness mode, Active power
1. 서 론
현재 압전 에너지 하베스팅 기술은 대규모 발전보다는 소형기기의 반영구 충전 및 자가 구동을 위한 방향으로 연구가 집중되고 있다. 그러나, 특정 조건에서의
높은 에너지 밀도와 구조물 통합 가능성 덕분에, 도로, 철도, 교량 등 대형 인프라에 적용하여 분산형 에너지 생산원으로 활용하려는 연구가 지속되고
있다[1].
압전 하베스트 발전은 압전 직접효과를 이용하여 기계적 진동을 전기적인 에너지로 변환하는 방식이므로 우수한 압전 특성을 가지는 압전 세라믹스를 요구하고
있다. 사용되는 재료로는 PMN, PVDF, PZT 등이 있으며, PZT 세라믹스가 다른 재료에 비해 가장 우수한 에너지 변환 효율을 가지고 있기
때문에 산업계에서 PZT계 세라믹스를 주로 사용하고 있다. 소자의 사용처에 따라서 특정 특성을 강화시키기 위해 소자 조성을 다양화하는데, 압전 하베스트
발전의 경우에는 FOM(Figure of Merit) = d33 * g33[2]도 중요하지만, 실제 응용 환경에서 진동원 주파수가 공진 주파수와 일치해야 에너지 수확 효율을 최대로 높일 수 있기에 압전 소자의 공진-반공진 주파수
대역폭 또한 중요하다[3].
본 연구에서는 에너지 변환 효율과 구조적 안정성, 그리고 실제 인프라 적용의 용이성을 고려하여 디스크형 압전 소자를 모델로 선정하였다. 디스크형 구조는
기계적 또는 전기적 입력에 대해 변형과 응력 분포가 균일하게 발생하기 때문에 높은 일관성과 재현성을 확보할 수 있다. 또한 직경/두께 비와 같은 간단한
설계 변경만으로 출력 특성의 조정이 가능하기에 공진 특성 설계와 최적화에 유리하다[4].
본 연구에서는 Pb(Ni1/3Nb2/3)O3-Pb(Zr Ti)O3(PMN-PZT)계 세라믹 조성을 선정하여[5] 디스크형 압전 소자로 직경/두께 비($\phi$/T)에 따라 유한요소 해석 소프트웨어인 COMSOL Multiphysics를 이용해 여러 모델로
출력 특성을 시뮬레이션하고 실제 압전소자의 특성과 비교하고자 한다.
2. 본 론
2.1 두께 진동 모드 압전 소자의 설계
본 연구에서는 유한요소 해석 프로그램인 COMSOL Multiphysics를 사용하여 PMN-PZT계 압전소자의 물성치를 기반으로 디스크형 소자를
모델링하여 주파수 영역에서 임피던스 분석을 실시하였다. 공진($f_r$)-반공진($f_a$) 주파수를 기반으로 전기기계 결합계수($k_{33}$)를
도출하였다[6].
다음 식(1)은 $k_{33}$의 계산에 적용된 식이다[7].
에너지 하베스팅을 위한 두께방향 진동모드를 갖는 압전소자 모델은 그림 1과 같은 형태로 설계하였다. 실제 응용 환경에서 바닥 면을 고정하여 사용한다는 가정하에 하부를 고정 구속조건으로 설정하고 Polarization 설정을
위해 소자의 하부를 접지로 설정하고 상부를 전위 터미널로 두고 초기값을 1 V로 설정했다[8]. 표 1은 모델링한 압전 소자에 적용된 재료 물성을 나타낸다.
그림 1. 두께 진동 모드 압전 소자 사양
Fig. 1. Specification of thickness vibration mode piezoelectric device
표 1. Pb(Ni1/3Nb2/3)O3–Pb(Zr Ti)O3의 물성표
Table 1. Physical properties of Pb(Ni1/3Nb2/3)O3–Pb(Zr Ti)O3
|
Material constants
|
Pb(Ni1/3Nb2/3)O3–Pb(Zr Ti)O3 |
|
Density($kg/m^3$)
|
7800
|
|
$Q_m$
|
70.06
|
|
Dielectric constant
|
1660
|
|
$s_{11}^E (10^{-12} m^2/N)$
|
16.4
|
|
$s_{12}^E (10^{-12} m^2/N)$
|
-5.74
|
|
$s_{13}^E (10^{-12} m^2/N)$
|
-7.22
|
|
$s_{33}^E (10^{-12} m^2/N)$
|
18.8
|
|
$s_{44}^E (10^{-12} m^2/N)$
|
47.5
|
|
$d_{15} (10^{-12} C/N)$
|
658
|
|
$d_{13} (10^{-12} C/N)$
|
-147
|
|
$d_{33} (10^{-12} C/N)$
|
387
|
2.2 시뮬레이션 결과
그림 2는 COMSOL Multiphysics 시뮬레이션을 통하여 구현한 두께 방향 진동모드로 동작하는 압전 소자의 진동 형태이다. 해당 그림은 직경($\phi$)
28mm, 두께(T) 6.5mm의 압전 소자를 외부 압력 주파수 영역에서 해석한 것으로 각각 (a) 84kHz, (b) 85kHz, (c) 86kHz에서의
변위 크기를 색상 분포로 나타낸 것이다. 84kHz, 86kHz에서의 최대변위가 각 1.87×10-5mm, 1.32×10-5mm로 공진주파수인 85kHz에서
최대변위가 4.73×10-5mm로 다른 주파수보다 최대 변위가 크고 가장자리 진동이 뚜렷이 나타나는 것을 볼 수 있다[9].
그림 2. 시뮬레이션을 통한 압전 하베스팅 소자 진동 현상
Fig. 2. Simulated vibration behavior of a piezoelectric harvesting element
압전 하베스팅 소자의 공진 특성을 검토하기 위하여 그림 3과 같이 두께 6.5mm인 소자를 대상으로 시뮬레이션을 수행하고 임피던스 특성을 구하였다. 이는 실제 제작한 6.5mm의 압전 세라믹 하베스팅 소자의
공진을 구하였을 때 잘 일치하고 있어 시뮬레이션의 신뢰성을 기대할 수 있게 하였다. 공진 주파수가 85kHz로 그림 2에서 최대 변위가 나타난 주파수와 일치함을 확인할 수 있으며[10], 이로부터 공진 주파수($f_r$)에서 최대 변위가 발생한다는 것을 알 수 있다.
그림 3. 시뮬레이션을 통한 압전 소자의 임피던스 특성 곡선
Fig. 3. Simulated impedance characteristic curve of a piezoelectric element
그림 4 는 $\phi$/T비의 변화에 따라 공진주파수($f_r$)의 변화를 나타낸 것으로 $\phi$/T비가 증가함에 따라 공진주파수가 다소 감소하고 있음을
알 수 있다. 이는 $\phi$/T비가 클수록 경방향 진동에 의한 간섭현상이 발생되고 있음을 보여준다. 또한 그림 5에 전극 터미널(Forced)의 전압(V)과 전류(A) 값을 COMSOL 자체 기능을 이용하여 찾아내고 realdot함수로 계산한 값을 통하여 얻은
유효전력(P)를 나타내고 있다. 유효전력(P)는 $\phi$/T가 4까지는 증가하다가 그 이상에서는 감소하는 특성을 나타내고 있어 본 연구에서는 표 3과 같이 $\phi$/T를 세부적으로 나누어 시뮬레이션 결과를 토대로 $\phi$/T가 4.3이 되도록 시편을 제작하였다. 그림 6은 전기기계 결합계수($k_{33}$)의 시뮬레이션 결과를 나타낸 것으로 $\phi$/T가 증가함에 따라 낮아지는 경향을 확인할 수 있다. $\phi$/T가
4.0 지점에서 유효전력이 가장 크게 나타나 이 지점을 중심으로 세부적으로 시뮬레이션 한 결과를 표 2로 나타내었다[11].
그림 4. $\phi$/T에 따른 공진 주파수
Fig. 4. Resonant frequency according to $\phi$/T ratio
그림 5. $\phi$/T에 따른 유효전력
Fig. 5. Active power according to $\phi$/T ratio
그림 6. $\phi$/T에 따른 전기기계 결합계수
Fig. 6. Electromechanical coupling factor according to $\phi$/T ratio
표 2는 두께(T)를 6.5mm로 고정하고 직경($\phi$)을 변화시킬 때 COMSOL Multiphysics 시뮬레이션을 통하여 그림 3~6에서 얻어진 압전소자의 공진 주파수와 유효전력 및 를 정리하여 나타낸 것이다.
표 2. $\phi$/T에 따른 공진-반공진 주파수, 유효전력 및 $k_{33}$
Table 2. Resonant and anti-resonant frequency, active power and $k_{33}$ as a function
of $\phi$/T
Diameter
($\phi$)[mm]
|
Thickness
(T)[mm]
|
Ratio
($\phi$/T)
|
Resonant
Frequency
($f_r$)[kHz]
|
Anti-Resonant
Frequency
($f_a$)[kHz]
|
Active
Power
[W]
|
$k_{33}$
|
|
19.5
|
6.5
|
3.0
|
97
|
124
|
0.0115
|
0.66
|
|
26
|
6.5
|
4.0
|
87
|
103
|
0.0160
|
0.57
|
|
28
|
6.5
|
4.3
|
85
|
99
|
0.0184
|
0.55
|
|
32.5
|
6.5
|
5.0
|
80
|
91
|
0.0093
|
0.51
|
그림 7은 표 2에 해당하는 소자들에 대해 시뮬레이션을 통하여 얻은 임피던스 특성 곡선을 보여준다. $\phi$/T의 수치가 증가함에 따라 공진-반공진 주파수의 대역대와
대역폭이 낮아지는 경향이 나타났다.
그림 7. $\phi$/T에 따른 임피던스 특성 곡선
Fig. 7. Impedance characteristic curve according to $\phi$/T ratio
표 3은 두께(T)를 6.5mm로 고정하고 직경($\phi$)을 변화시켜 압전소자의 공진주파수, 유효전력 및 을 $\phi$/T = 4.0과 4.5 구간에서
최적의 비율을 얻고자 $\phi$/T = 0.1 간격으로 시뮬레이션 결과를 세분화하여 나타낸 것이다. $\phi$/T = 4.2인 부분에서 가장 높은
유효전력 수치를 얻을 수 있었다.
표 3. $\phi$/T에 따른 공진-반공진 주파수, 유효전력 및 $k_{33}$
Table 3. Resonant and anti-resonant frequency, active power and $k_{33}$ as a function
of $\phi$/T
Diameter
($\phi$)[mm]
|
Thickness
(T)[mm]
|
Ratio
($\phi$/T)
|
Resonant
Frequency
($f_r$)[kHz]
|
Anti-Resonant
Frequency
($f_a$)[kHz]
|
Active
Power
[W]
|
$k_{33}$
|
|
26.65
|
6.5
|
4.1
|
85.5
|
100.5
|
0.0155
|
0.56
|
|
27.3
|
6.5
|
4.2
|
85.5
|
100
|
0.0196
|
0.55
|
|
28
|
6.5
|
4.3
|
85
|
99
|
0.0184
|
0.55
|
|
28.6
|
6.5
|
4.4
|
84
|
94.5
|
0.0191
|
0.49
|
|
29.25
|
6.5
|
4.5
|
83.5
|
96
|
0.0072
|
0.53
|
2.3 실험 시편 소자 제작
표 2와 표 3에서 두께 진동모드에서의 소자 설계에서 유효 전력 수치가 높은 소자는 $\phi$/T=4.2라는 것을 알 수 있었다. 그러나 $\phi$/T=4.3,
$\phi$/T=4.4에서의 유효전력 수치가 큰 차이가 나지 않았고 소결을 진행하는 압전 소자의 직경을 고려하여 $\phi$/T=4.3인 소자를 시편으로
제작하였다.
그림 8은 제작된 시편의 소결 온도에서 얻은 SEM 사진이다. 평균 입자의 크기가 6.28$\mu$m로 고르게 입자가 성장한 것을 관측할 수 있다.
그림 9는 앞서 채택된 $\phi$/T=4.3의 비율인 압전 소자를 시편으로 제작한 것이다. 해당 시편의 임피던스 측정을 위해 그림 1에서와 동일하게 시편 양면에 전극을 부착하였고, 전극에 전선을 연결하여 다음과 같은 시편 소자를 만들었다.
그림 8. 소결 온도에서의 SEM 사진
Fig. 8. SEM micrograph at the sintering temperature
그림 9. 제작된 시편 소자
Fig. 9. Manufactured specimens
2.4 실험 및 특성 비교 분석
그림 10은 시뮬레이션 결과와 제작한 시편 소자와의 비교를 위해 최적의 모델로 채택된 두께 6.5mm일 때 $\phi$/T=4.3 비율의 압전소자를 제작한
후 임피던스 분석기로 측정한 소자 임피던스를 시뮬레이션의 임피던스 특성과 비교한 그래프이다. 측정 결과와 시뮬레이션 결과가 비교적 잘 일치하고 있음을
알 수 있다.
표 4는 제작된 시편의 물성치를 측정한 결과이다. 이때 측정된 시편의 공진 및 반공진 주파수는 각각 86.445kHz, 98.921kHz로 나타났으며 전기기계
결합계수($k_{33}$)는 0.52로 나타났다. 시뮬레이션 결과와 비교하면 다소 오차가 있는데 이는 시뮬레이션을 수행한 소자의 물성과 실제 제작한
압전소자의 가공 정밀도로 인해 생긴 것으로 추정된다.
그림 10. 시뮬레이션과 제작된 시편의 임피던스 특성 곡선
Fig. 10. Impedance characteristic curves of simulation and manufactured specimens
표 4. 제작된 시편의 물성 측정표($\phi$/T=4.3)
Table 4. Physical properties of measured specimen ($\phi$/T=4.3)
|
|
$f_r$[kHz]
|
$f_a$[kHz]
|
C[nF]
|
$k_{33}$
|
P[mW]
|
Piezoelectric
element
|
86.445
|
98.921
|
1.2661
|
0.52
|
1.4636
|
그림 11과 그림 12는 압전 소자의 충격 응답을 측정하기 위해 구성된 실험 장치의 블록도 및 실제 사진을 나타낸 것이다. 스프링 장력을 조절하여 압전 세라믹에 가해지는
충격량 및 타격 시간을 제어하였으며, 인가 하중은 로드셀을 통해 측정하였다. 충격 인가 후 세라믹 시편 내부에서 발생한 전기적 신호를 오실로스코프를
이용하여 계측하였으며, 측정된 출력 파형의 분석을 통해 최종적으로 발생한 전기 에너지를 정량적으로 도출했다.
그림 11. 충격 응답 측정을 위한 실험 장치 블록도
Fig. 11. Block diagram of the experimental testing setup for measuring the impact
response
그림 12. 레버와 스프링 시스템이 적용된 자체 제작 수직 충격 인가 장치
Fig. 12. Self-built vertical impact application device with lever and spring system
제작된 시편에 자체 제작한 스프링 수직 충격 인가 장치를 통해 40N의 힘을 가했을 경우 최대 전압(V)이 오실로스코프로 측정했을 때 340V로 나타났다.
식 2에 따라 최대 피크 전압($V_{peak}$)은 534V이어야 하지만 이보다 적은 수치로 나타났다. $Q_m$은 기계적 증폭비로 감쇠에 반비례한다[12].
식 3에 따라서 C(capacitance)가 1.2661nF이고 피크전압에 따라 발생된 전류는 I는 4.30474×10-6A이다. 그러므로 실험에서 계산된 유효전력 P는 1.4636116mW로 시뮬레이션 결과에 비해 적은 수치로 나타났다[13].
3. 결 론
본 연구에서는 직경/두께 비에 따른 두께 진동모드 압전 하베스터 소자의 공진-반공진 주파수와 유효전력 및 전기기계 결합계수를 비교하였다.
압전소자의 물성은 PMN-PZT계 세라믹의 값을 적용하였으며 유한요소해석 프로그램인 COMSOL Multiphysics를 이용하여 시뮬레이션으로 최적의
비율을 찾은 후 압전소자를 제작하여 그 특성을 비교하였다.
두께 진동모드를 갖는 두께 6.5mm의 압전소자에서 $\phi$/T비가 증가할수록 공진 주파수가 낮아짐을 알 수 있었다. 또한, 시뮬레이션으로부터
$\phi$/T=4.3의 시편에 대해 높은 유효전력 값을 얻었으며, 해당 비로 제작한 시편의 물성이 시뮬레이션과 실험 결과 간 일치함을 알 수 있었다.
본 연구는 압전소자의 설계 및 제작 시 두께 진동모드에서 직경 및 두께의 조건이 공진-반공진 주파수와 출력 특성에 영향을 미친다는 사실을 제시하고
있으며, 이러한 시뮬레이션 수치가 실험치와 잘 일치함을 확인하였다. 따라서 향후 다양한 환경에서 사용될 압전 하베스트 소자의 출력 특성을 예측 설계하기
위한 소자 출력 설계 예측 모델로 연구에 활용될 수 있을 것으로 기대된다.
Acknowledgements
Following are results of a study on the "Busan Regional Innovation System & Education(RISE)"
Project, supported by the Ministry of Education and Busan Metropolitan City (2025-RISE-02-033-034)
References
K. Y. Kong, D. U. Kim, M. H. Kim, J. P. Won, "The Evaluation of Energy Harvesting
of Concrete Composite," pp. 597-598, 2021.

J. E. Kim, "Analysis of Vibration - powered Piezoelectric Energy Harvesters by Using
Equivalent Circuit Models," Transactions of the Korean Society for Noise and Vibration
Engineering, vol. 20, no. 4, pp. 397-404, 2010.

D. Lee, Y. Yee, H. Song, "Broadband Piezoelectric Energy Harvesting Technology," CERAMIST,
vol. 22, no. 1, pp. 56-69, 2019.

Xu. Chaoqun, Li. Yuanbo, Yang. Tongqing, "Optimization of Non-Uniform Deformation
on Piezoelectric Circular Diaphragm Energy Harvester with a Ring-Shaped Ceramic Disk,"
Micromachines, vol. 11, no. 11, pp. 963, 2020.

G. Kim, J. Yoo, W. S. A, L. S. HO, P. HERIE, I. I. HO, C. W. Oh, "Dielectric and Piezoelectric
Characteristics of Pb(Ni1/3Nb2/3)O3–Pb(Zr Ti)O3 System Ceramics for the Application
of Energy Harvesting Device," Journal of the Korean Institute of Electrical and Electronic
Material Engineers, vol. 38, no. 5, pp. 580-585, 2025.

S. Lim, J. Na, J. Hwang, M. Lee, S. Lee, H. Park, "Analysis of Output Characteristics
of a Disk-Type Piezoelectric Energy Harvester," pp. 1880-1881, 2025.

C. Joo, T. Kim, G. Baek, S. Lee, H. Ahn, "A Study on the Underwater Environment Load
Modeling and Impedance Characteristics Analysis of Langevin Type Piezoelectric Transducer,"
The transactions of The Korean Institute of Electrical Engineers, vol. 71, no. 12,
pp. 1812-1819, 2022.

M. S. Kwak, G. Hwang, "Structural Analysis Simulation of Cantilever Shaped Piezoelectric
Energy Harvester Using COMSOL Multiphysics," Journal of the Korean Institute of Electrical
and Electronic Material Engineers, vol. 34, no. 6, pp. 416-425, 2021.

Y. H. Jeong, "Trajectory of Resonant Displacement of Thickness Vibration Mode Piezoelectric
Devices According to Diameter/Thickness Ratio," Journal of the Korean Institute of
Electrical and Electronic Material Engineers, vol. 25, no. 2, pp. 105-109, 2012.

J. Yoo, J. Hong, Y. Jeong, M. H. Park, "Resonant Displacement and Piezoelectric Properties
of Thickness Shear Mode Piezoelectric Devices According to Length/Thickness Ratio,"
Journal of Electrical and Electronic Materials, vol. 24, no. 6, pp. 463-467, 2011.

J Ajitsaria, "Modeling and Analysis of a Bimorph Piezoelectric Cantilever Beam for
Voltage Generation," Smart Mater. Struct., vol. 16, pp. 447-454, 2007.

K. Y. Kong, "Energy harvesting performance of composites using piezoelectric elements,"
Master's Thesis, Graduate School, Konkuk University, 2021.

Weijian Zhou, Biao Wang, C.W. Lim, Zhengbao Yang, "A distributed-parameter electromechanical
coupling model for a segmented arc-shaped piezoelectric energy harvester," Mechanical
Systems and Signal Processing, vol. 146, pp. 107005, 2021.

저자소개
Sung-Hyun Lim received the B.S. degree in IT Convergence & Application Engineering
from Pukyong National University, Busan, Korea, in 2021. He is currently working toward
the M.S. degree in the school of electrical engineering at Dong-A University, Busan,
Korea.
E-mail: dlatjdgus955@gmail.com
Su-Ho Lee received the B.S., M.S., and Ph.D. degrees in electrical engineering from
Dong-A University, Busan, South Korea, in 1989, 1991, and 1996. He is currently a
Professor in the Department of Electrical Engineering at Dong-A University, Korea.
E-mail: leesuho@dau.ac.kr
Jun-Young Na received the B.S. degree in Electrical Engineering from Dong-A University,
Busan, Korea, in 2025. Currently, he is working toward the M.S. degree in the school
of Electrical Engineering at Dong-A University, Busan, Korea.
E-mail: wnsdud0496@donga.ac.kr
Hwang Jae-Eun received the B.S. degree in Department of Energy Science from Kyungsung
University, Busan, Korea, in 2018. He received the M.S. degree in Electrical Engineering
from Dong-A University, Busan, Korea, in 2022. Currently, he is working toward the
Ph.D. degree in the school of Electrical Engineering at Dong-A University, Busan,
Korea.
E-mail: 2273084@donga.ac.kr
Jae-Sung Yoo received the B.S. degree in Production and Processing Engineering from
Pukyong National University, Busan, Korea, in 2002. And he received the M.S. degree
in Department of Electrical, Electronic and Computer Engineering from Dong-A University,
Busan, Korea, in 2016. Currently, he is working toward the Ph.D. degree in the school
of Electrical Engineering at Dong-A University, Busan, Korea.
E-mail: pknujs@gmail.com
She received the Ph.D. degree from Cergy Paris Universit, France, and Yeungnam University,
Korea, in 2013. She is currently an Associate Professor in the Department of Electrical
Engineering at Dong-A University, Korea. Her research interests include electrical
insulation, complex materials, and energy management.
E-mail : bakery@donga.ac.kr