CONTENTS Accuracy of a direct estimation method for equivalent material properties of 1-3 piezocomposites ················· Eunghwy Noh, Donghyeon Kim, Hyeongmin Mun, Woosuk Chang, Hongwoo Yoon, Seonghun Pyo, ···············································································Kyungseop Kim, and Yo-Han Cho377 A study on the arrangement of actuators and speaker zones of the panel speaker ·············································· Jung-Han Woo, Seong-Hyun Lee, Yun-Ho Seo, Pyung-Sik Ma, and Dongjoon Kim388 A study on training DenseNet-Recurrent Neural Network for sound event detection ······················································································································ Hyeonjin Cha and Sangwook Park395 ■ Special Issue on Ultrasound-based Multimodal Imaging and Treatment Shear-wave elasticity imaging with axial sub-Nyquist sampling ···························································································································· Woojin Oh and Heechul Yoon403 Cell-cultivable ultrasonic transducer integrated on glass-coverslip ··················································································································· Keunhyung Lee and Jinhyoung Park412 Adaptive quantization for effective data-rate reduction in ultrafast ultrasound imaging ······················································································································· Doyoung Jang and Heechul Yoon422 Ultrasound-optical imaging-based multimodal imaging technology for biomedical applications ······················ Moon Hwan Lee, HeeYeon Park, Kyungsu Lee, Sewoong Kim, Jihun Kim, and Jae Youn Hwang429 Optically transparent ultrasound transducers for combined ultrasound and photoacoustic imaging: A review ······················································································································ Shunghun Park and Jin Ho Chang441 Single beam acoustic tweezers for biomedical applications ···················································································································································· Hae Gyun Lim452 A comparative study on keypoint detection for developmental dysplasia of hip diagnosis using deep learning models in X-ray and ultrasound images ··························· Sung-Hyun Kim, Kyungsu Lee, Si-Wook Lee, Jin Ho Chang, Jae Youn Hwang, and Jihun Kim460 Acoustic outputs from clinical extracorporeal shock wave lithotripsy devices ··· Jong Min Kim, Oh Bin Kwon, Jin Sik Cho, Sung Joung Jeon, Ki Il Nam, Sung Yong Cho, and Min Joo Choi469 ▪Society News and Information ······················································································································ i 본 사업은 기획재정부의 복권기금 및 과학기술정보통신부의 과학기술진흥 기금으로 추진되어 사회적 가치 실현과 국가 과학기술 발전에 기여합니다. THE ACOUSTICAL SOCIETY OF KOREA Vol.42, No.5September 2023I. 서 론 압전재료는 전기 신호를 음향 신호로 변환하거나 반대로 음향 신호를 전기 신호로 변환하는 센서 구 동 재료로서 주로 수중 소나, 의료용 초음파 및 비파 괴 검사 등에 사용된다. [1,2] 그 중에서도 1-3형 압전세 라믹-폴리머 복합재료(이하 압전복합재료)는 압전 세라믹 기둥이 폴리머 기질 내에 배열되고 상하면 에 전극이 입혀진 구조이다. 압전세라믹 재료 자체 는 기본적으로 두께 방향과 측면 방향의 압전상수 가 서로 반대의 부호를 가지므로 저주파의 음압이 안가될 때 서로 상쇄되어 수신감도가 떨어지는 단 1-3형 압전복합재료 등가물성 직접 추출 기법의 정확도 분석 Accuracy of a direct estimation method for equivalent material properties of 1-3 piezocomposites 노응휘, 1† 김동현, 1 문형민, 1 장우석, 1 윤홍우, 1 표성훈, 2 김경섭, 2 조요한 2 (Eunghwy Noh, 1† Donghyeon Kim, 1 Hyeongmin Mun, 1 Woosuk Chang, 1 Hongwoo Yoon, 1 Seonghun Pyo, 2 Kyungseop Kim, 2 and Yo-Han Cho 2 ) 1 LIG넥스원, 2 국방과학연구소 (Received June 22, 2023; revised July 24, 2023; accepted August 3, 2023) 초 록: 본 논문에서는 1-3형 압전복합재료를 단일상의 균질 압전매질로 모델링하기 위해 필요한 등가의 물성을 유한 요소 해석으로 직접 추출하는 기법의 정확도를 다룬다. 직접 추출 기법은 압전재료의 전기적, 기계적 거동과 상호 간 커플링을 기술하는 구성방정식을 기반으로 물성 행렬의 개별 성분들을 직접적으로 산출하는 방법이다. 직접 추출에 사용되는 두 가지 구성방정식 조합 간의 정확도를 비교하기 위하여 단일 1-3형 압전복합재료 하이드로폰을 대상으로 등가물성과 수신감도를 산출하고 전체 영역에 대한 유한요소 해석 결과와 비교한다. 물성 추출의 정확도는 압전복합 재료를 구성하는 폴리머의 탄성 특성과 밀접한 관련이 있음을 확인하고, 오차 원인을 분석하여 정확한 등가물성 산출 을 위한 가이드라인을 제시한다. 압전복합재료 하이드로폰과 주변의 음향구조물을 포함하는 스테이브 규모에 대해서 도 등가 모델링을 적용하여 추출된 물성의 정확도와 연산량 감소를 확인한다. 핵심용어: 1-3형 압전복합재료, 등가물성, 직접 추출 기법, 유한요소 해석 ABSTRACT: This paper presents accuracy of a method that directly estimates equivalent properties of a 1-3 piezocomposite for modeling it into the single phase homogeneous piezomaterial. This direct estimation method finds individual components of a material property matrix based on the piezoelectric constitutive equations, which represent mechanical and electrical behaviors and their couplings. Equivalent properties on a single 1-3 piezocomposite hydrophone are derived, and their accuracy depending on pairing of the constitutive equations is investigated by comparing them with finite element analysis for the whole domain. The accuracy is related to elastic characteristics of a matrix polymer, and the error is analyzed so that some guidelines for correct estimation are suggested. Fidelity of estimated properties and equivalent modeling is shown in a stave scale including hydrophones and surrounding acoustic structures as well, and reduced computational cost is verified. Keywords: 1-3 piezocomposite, Equivalent material property, Direct estimation method, Finite element analysis PACS numbers: 43.38.Ar, 43.38.Fx 한국음향학회지 제42권 제5호 pp. 377～387 (2023) The Journal of the Acoustical Society of Korea Vol.42, No.5 (2023) https://doi.org/10.7776/ASK.2023.42.5.377 pISSN : 1225-4428 eISSN : 2287-3775 †Corresponding author: Eunghwy Noh (eunghwy.noh@lignex1.com) Maritime Sensor System, LIGNex1 & 333, Pangyo-ro, Bundang-gu, Seongnam-si, Gyeonggi-do 13488, Republic of Korea (Tel: 82-31-5178-4615, Fax: 82-31-5179-7100) Copyrightⓒ2023 The Acoustical Society of Korea. This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited. 377노응휘, 김동현, 문형민, 장우석, 윤홍우, 표성훈, 김경섭, 조요한 한국음향학회지 제 42 권 제 5 호 (2023) 378 점이 있다. 반면 압전복합재료는 두께 방향 압력은 세라믹 기둥의 두께 방향 압전 특성으로 최대한 수 신하되, 측면 방향 압력은 세라믹 대비 연한 폴리머 재료의 탄성 특성으로 저감시킨다. 압전복합재료는 세라믹에 비해 두께 대비 감도가 높으며, 음향 임피 던스는 낮아 변환 효율이 좋은 음향센서 소재이다. 압전복합재료는 비등방성이자 비균질 재료이기 때문에 이를 단일상의 균질한 압전재료 모델로 대치 하기 위한 등가물성 추출 연구들이 수행되어 왔다. 초기에는 주로 이론적인 방법들이 사용되었다. Newnham et al. [3] 은 압전복합재료의 구성 요소들 간 구조적, 물리적 연결성에 따른 주요 등가물성을 산 출하였다. Smith와 Auld [4] 는 1-3형 압전복합재료의 횡방향 스케일이 파장 대비 매우 커서 두께 거동이 지배적인 경우, Smith [5] 는 재료의 스케일이 작아 정 수압 거동이 지배적인 경우에 대하여 등가물성을 추 출하였다. Avellaneda와 Swart [6] 는 정수압 거동에 대 해 복합재료 내부 미세 구조의 특성을 반영할 수 있 는 계수를 추가 적용한 유효 매질 모델을 제시하였 다. 이러한 이론적인 등가 모델링 기법들은 적절한 가정이 필요하며 그에 따라 주요한 일부 물성들만 추출이 가능하다. 또한 세라믹 부피비가 낮거나 폴 리머의 포아송 비가 낮은 경우에는 오차가 발생할 수 있으며, [7] 추가적인 구조물의 영향 등을 반영하기 어렵다. 다른 물성 추출 기법으로 압전복합재료의 공진 모드별 특성을 실험 혹은 유한요소 해석으로 구한 후 해당 공진 모드에 기여하는 주요 물성들을 도출하는 공진법을 적용한 사례들이 있다. [8-10] 공진 법은 해당 모드 주파수 부근에서 정확도가 높지만 그 외의 대역에서는 정확도가 낮을 수 있으므로, 공 진을 피해 저주파에서 운용되어 정수압 거동을 보이 는 소형 압전복합재료에 적용하기는 적합하지 않을 수 있다. 유한요소 해석을 이용한 직접 추출 기법은 압전재 료 구성방정식(piezoelectric constitutive equations)을 기반으로 등가물성 행렬의 개별 성분들을 직접적으 로 계산하는 방법이다. 직접 추출 기법은 형상이 바 뀔 때마다 해석을 새로 수행해야 한다는 단점이 있 으나, 물성 행렬의 전체 성분을 완성할 수 있어 재료 단위 이상의 스케일에 대한 등가 모델링 및 성능 예 측에 용이하다. 또한 수학적으로는 별도의 변형 없 이 구성방정식 자체를 따라 진행되기 때문에 비교적 직관적이며, 추가적인 가정이 요구되지 않으므로 다 양하고 복잡한 형상에 대해 강건하게 적용이 가능하 다. 기존에 유한요소 해석을 이용해 등가물성을 직 접 추출한 사례들이 있으나, [7,11-13] 정수압 압전상수 등의 일부 물성만을 계산하였다. 이는 대부분 기존 의 이론적인 물성추출 기법을 검증하거나 [7] 설계 모 델 혹은 시험결과의 신뢰성을 확보하는 목적으로 사 용되었지만, [11-13] 인가 조건이나 폴리머 특성에 따른 기법 자체의 정확도에 대한 분석적인 논의는 아직까 지 이루어지지 않은 것으로 보인다. 이에 본 논문에 서는 유한요소해석을 이용한 등가물성 직접 추출 기 법의 정확도를 중점적으로 분석해보고자 한다. 직접 추출 기법의 기반이 되는 구성방정식의 조합에 따라 두 가지 방식으로 압전복합재료의 등가물성을 추출 한 후, 단일 하이드로폰 소자의 자유음장 수신감도 측면에서 정확도를 비교한다. 또한 압전복합재료 기 질의 폴리머 재료 특성에 따라 추출된 물성과 정확 도가 어떻게 달라지는지 확인하고 원인을 분석한다. 마지막으로, 추출된 물성을 이용해 음향구조물을 포 함하는 스테이브 규모에서 등가 모델링의 정확도와 연산량을 검증해보고자 한다. II. 경질 폴리머 압전복합재료의 등가물성 추출 및 정확도 분석 압전복합재료를 대치하고자 하는 단일상 압전재 료의 구성방정식은 기계적 변수인 응력 T 와 변형률 S , 전기적 변수인 전기장 E 와 전기변위 D 사이의 관 계식이다. 이는 독립변수와 종속변수 조합에 따라 네 가지 방식으로 표현할 수 있는데, [14] 본 논문에서 는 유한요소 해석에서 주로 사용되는 두 가지 형식 에 대해서만 다루고자 한다. 첫 번째 형식은 응력 T 와 전기장 E 를 우변의 독립변수로, 변형률 S 와 전기 변위 D 를 좌변의 종속변수로 표현하는 방법으로서 Eq. (1)과 같이 행렬식으로 나타낼 수 있는데, 이를 종 속변수명을 기준으로 strain-charge form 혹은 S-D form 이라 칭한다.1-3형 압전복합재료 등가물성 직접 추출 기법의 정확도 분석 The Journal of the Acoustical Society of Korea Vol.42, No.5 (2023) 379 S s E T d t E D dT T E (1) 여기서 s 는 기계적 컴플라이언스, d 는 S-D form 압전 상수, 는 유전상수이다. 대문자로 쓰인 위첨자는 해 당 변수의 값이 일정하다는 것을, 소문자로 쓰인 위 첨자 t는 전치행렬(transpose)을 의미한다. 두 번째 형 식은 stress-charge form 혹은 T-D form이다. 변형률 S 와 전기장 E 는 독립변수로, 응력 T 와 전기변위 D 를 종속변수로 하여 T c E S e t E D eS S E (2) 와 같이 기술되며, c 와 e 는 각각 탄성계수와 T-D form 압전상수이다. Eq. (1)과 마찬가지로 위첨자 t는 전치행렬을, 위첨자 E 와 S는 각각 전기장과 변형률 이 일정함을 의미한다. 압전재료는 방향별 특성이 다르기 때문에 각 변수들은 행렬의 형태로 표현되므 로 Eqs. (1)과 (2)에서 이를 볼드체로 표기하였다. 등가물성 직접 추출 기법은 이러한 압전재료 구성 방정식을 기반으로 수행된다. 먼저 S-D form을 이용 한 물성 추출은 Eq. (1)의 독립변수인 T 또는 E 에 단 위벡터 형태의 행렬을 입력했을 때 출력되는 종속변 수 S 와 D 로부터 물성 행렬 s E, d , T의 각 성분들을 구하는 방식이다. 행렬 T 와 E 를 구성하는 총 9개의 성분 중에서 한 가지만 1이고 나머지는 모두 0이 되 는 조건에서 압전복합재료의 유한요소 해석을 수행 한 뒤, S 와 D 의 각 성분들을 구하면 물성을 추출할 수 있다. 예를 들어 s E 행렬의 (1,1) 성분인 의 값은, 1방향 축하중을 제외한 나머지 5개 응력과 3개 전기 장 값이 모두 0이며 1방향 단위 응력( )만이 작용하는 조건에서 계산된 1방향 변형률 과 같다. 마찬가지로 T-D form 기반의 물성 추출 방식에서는 독립변수인 S 또는 E 내에 하나의 성분만이 1이 되 는 조건에서 해석을 통해 T 와 D 의 값을 산출한다. 이러한 과정을 각 성분별로 반복함으로써 물성 행렬 을 채워나갈 수 있다. 직접 추출 기법을 적용하여 등가물성을 추출하고 자 하는 1-3형 압전복합재료 하이드로폰 소자의 형 상은 Fig. 1과 같다. 하이드로폰은 공진을 피해 저주 파 대역에서 구동되는 것을 상정하였다. 최대 운용 주파수 에서의 파장을 이라 할 때, 소자의 가로 와 세로 길이는 약 0.17 , 높이는 약 0.07 정도로 작은 크기이며 약 3.1의 종횡비율을 갖는다. 내부의 원형 세라믹 기둥 직경은 약 0.02 이며, 5×5 형태로 등간격 배치되었고 중심간 거리는 약 0.03 이다. 세 라믹의 재질은 PZT-5 계열이고, 주변의 기질부에는 에폭시 계열의 경질 폴리머(II장) 또는 우레탄 계열 의 연질 폴리머(III장)를 적용하였으며 폴리머의 부 피비는 약 70 %가 되도록 형상을 선정하였다. 두께 방향 효과를 극대화하기 위해 상하에 고강성의 평판 이 부착되었으며 그 두께는 최대 운용주파수에서의 파장 대비 약 0.007로 매우 얇다. 전극은 내부 복합재 료 구조와 평판 사이에 위치한다고 가정하였으며 별 도의 형상 모델링은 하지 않았다. 각 재료의 물성은 Table 1에 정리하였으며, 상하 평판이 압전복합재료 특성에 큰 영향을 미치기 때문에 평판을 포함하는 소자 전체 형상에 대하여 등가물성을 추출하였다. 1-3형 압전복합재료의 구조를 고려하여 육방정의 구조를 갖는 C 6V (혹은 6 mm) 결정계를 기준으로 등 가물성을 도출하였다. 결정 구조의 대칭성을 반영한 S-D form 구성방정식을 행렬식으로 표현하면 아래 Fig. 1. (a) External feature and (b) internal configu- ration of a 1-3 piezocomposite hydrophone. Table 1. Properties of materials used in a piezo- composite hydrophone. Hard polymer matrix Soft polymer matrix Plate Density (kg/m 3 )16007601500 Young’s Modulus (GPa) 60.1360 Poisson’s ratio0.390.490.3 Dielectric constant44-노응휘, 김동현, 문형민, 장우석, 윤홍우, 표성훈, 김경섭, 조요한 한국음향학회지 제 42 권 제 5 호 (2023) 380 Eq. (3)과 같다. ,(3) 여기서 로 계산이 가능한 성분을 제외 하면 9×9의 등가물성 행렬은 총 10개의 독립적인 성 분으로 완전히 정의될 수 있다. 상기된 Eq. (3)에서 독립변수 T 와 E 로 이루어진 우변의 벡터 행렬이 입력변수, 즉 유한요소 해석의 경계조건이 되며, 종속변수 S 와 D 로 구성된 좌변의 벡터 행렬이 출력변수, 즉 해석으로 계산되는 결과 값이 된다. 세부적인 과정은 다음과 같다. ① 이 1, 나머지 8개 성분은 0 인 조건에서 유한요 소 해석을 수행한다. 이것을 Eq. (3)의 관점에서 입 력변수 행렬에 대입해보면 우변은 물성 행렬의 첫째 열 성분들만 남게 되며, 그 값은 좌변 행렬 성 분들의 값과 같아짐을 알 수 있다. ② ①의 해석 결과로부터 , , , 의 값을 산출 하면 각각 , , , 값이 된다. 이 때, 각각에 대해 유한요소해석을 따로 수행할 필요는 없다. ③ 이 1이고 나머지 성분은 0인 조건에서 유한요 소 해석을 수행한다. Eq. (3)의 관점에서는 우변 물 성 행렬의 세 번째 열 성분만 남게 되는 것이다. ④ ③의 결과로부터 , , 의 값을 산출하면 각각 , , 값이 된다. ⑤ 위 절차를 입력조건이 , , 인 경우 각각에 대해 반복하면서, Eq. (3)의 물성 행렬 성분들과 대응되는 출력변수들을 산출한다. 이러한 과정으로 10개의 성분 모두를 추출하기 위 해 필요한 유한요소 해석의 경우의 수는 5가지이다. 각각의 해석 케이스별로 인가되는 입력변수와 출력 되는 변수들에 따라 추출되는 물성들을 Table 2에 정 리하였다. 2방향 변수들은 1방향과 동일하므로 해석 케이스에서 제외하였고, , , , 와 같이 여 러 케이스에서 중복 추출되는 성분들은 모두 평균값 을 적용하였다. T-D form 기반의 물성추출 또한 이와 유사한 과정으로 진행되었으며, 본 논문의 모든 유 한요소 해석은 상용 소프트웨어인 COMSOL Multi- physics를 통해 수행되었다. 먼저 경질 폴리머가 적용된 압전복합재료에 대한 S-D form 기반 등가물성 추출의 한 예시로서, Table 2 의 1번 케이스에 해당하는 입력조건 및 대표적인 해석 결과를 Fig. 2에 나타냈다. 입력조건이므로 Fig. 2(a)와 같이 보라색으로 표현된 1방향의 양면에 1 Pa의 인장응력 경계조건을 설정하였다. 을 제외 한 나머지 응력은 모두 0이 되어야 하므로, 2방향과 3 Table 2. Estimated components of equivalent pro- perties depending on input and output variables. Case Input variable Output variables #1 -- #2 -- #3 --- - #4 --- - #5 - -- Fig. 2. (Color available online) An example of esti- mating properties of a piezocomposite with hard polymer based on the S-D form of constitutive equations.1-3형 압전복합재료 등가물성 직접 추출 기법의 정확도 분석 The Journal of the Acoustical Society of Korea Vol.42, No.5 (2023) 381 방향의 나머지 4개면에는 자유 경계조건을 주었다. 전기적으로는 , , 가 모두 0이 되어야 하므로 모든 외부 면에 접지 조건을 주었다. 이러한 조건에 서 유한요소 해석을 수행하였을 때 변위 벡터의 크 기, 총 변위 값은 Fig. 2(b)와 같으며, 변형된 형태를 함께 도시하였다. 1방향으로 인장이 발생하였고 특 히 폴리머 가장자리에서 다른 방향으로 구속이 없기 때문에 변형이 크게 나타났다. Fig. 2(c)와 (d)에서는 1 방향과 3방향 표면에서 변형률 , 값의 분포를 변형된 형태와 함께 도시하였다. 여기서 변형률은 각 표면에서 해당 방향으로의 변위를 초기 길이로 나눠줌으로써 계산되었다. 포아송 효과로 인해 1방 향의 인장이 2방향과 3방향의 압축 변형을 유발하였 으며, 상대적으로 평판은 영률이 높고 폴리머는 낮 기 때문에 Fig. 2(c)에서 폴리머 부분의 2방향 압축 변 형이 크게 나타났다. 등가물성 와 은 일종의 대 푯값이므로 각각 1방향과 3방향 표면의 변형률 분포 에 대한 면평균 값으로 산출하였다. Fig. 2에서 표현 된 것 이외에 다른 성분들에 대해서도 모두 면평균 값을 통해 등가물성이 산출되었다. s E , d 와 같이 출 력변수가 변형률 S 인 경우에는 해당 표면에서의 변 위를 초기 길이로 나눈 변형률 표면 분포를 구한 뒤 그 평균값을 계산하였고, d 나 와 같이 출력변수가 전기변위 D 인 경우에는 표면 전하밀도를 적분하여 면적으로 나눠준 값으로 등가물성을 산출하였다. 동일한 하이드로폰에 대하여, T-D form 기반 등가 물성 추출 중에서 이 입력조건일 때의 해석 결과 를 Fig. 3에 나타냈다. Fig. 3(a)와 같이 1방향의 양면 에 단위 변형률을 인가하였고, 나머지 표면에서는 1 방향으로 구속 없이 움직이되 표면에 수직한 방향의 변위는 0이 되는 롤러 경계조건을 설정함으로써 변 형률이 0이 되도록 하였다. 총 변위와 변형된 형태는 Fig. 3(b)와 같이 계산되었으며, 1방향과 3방향 표면 에서의 1방향 응력 의 분포는 각각 Fig. 3(c), (d)와 같이 나타났다. T-D form 기반의 등가물성 추출 또한 S-D form의 경우와 마찬가지로 각 입력조건에 따라 총 5가지 경우에 대하여 해석이 수행되었으며, c E나 e 와 같이 응력을 출력함으로써 얻어지는 등가물성 성분들은 해당 표면에서의 응력 면평균 값으로 산출 되었고, 중복되는 성분일 경우 두 값의 평균으로 계 산하였다. S-D form와 T-D form 두 가지 추출 방식의 정확도 는 단일 하이드로폰 소자의 자유음장 수신감도를 통 해 확인하였다. Fig. 1의 하이드로폰 형상과 동일한 외곽 치수를 갖는 단일상 압전재료에 각 방법으로 추 출한 등가물성을 입력하여 해석된 수신감도를 전체 영역에 대한 유한요소 모델(full Finite Element Model, full FEM)의 결과와 함께 Fig. 4에 비교하였다. Fig. 4(a)와 (b)는 각각 full FEM과 등가모델에 대해 자유 음장 수신감도를 계산하기 위한 해석 모델이다. Fig. Fig. 3. (Color available online) An example of esti- mating properties of a piezocomposite with hard polymer based on the T-D form of constitutive equations. Fig. 4. (Color available online) Models and their results of RVS, and first mode shapes of a piezo- composite hydrophone with hard polymer calculated by a full FEM and equivalent models using directly estimated properties.노응휘, 김동현, 문형민, 장우석, 윤홍우, 표성훈, 김경섭, 조요한 한국음향학회지 제 42 권 제 5 호 (2023) 382 4(a)의 Full FEM은 Fig. 1의 압전복합재료 하이드로폰 형상을 그대로 반영하였으며, Fig. 4(b)의 등가모델 은 동일한 외곽 치수를 가지며 앞서 추출된 등가물 성을 갖는 단일상의 압전체 형상을 반영하였다. 1Pa 의 평면파가 입사할 때 전압을 추출하여 자유음장 수신감도를 계산하였으며, 등가물성은 S-D form과 T-D form 각각을 이용해 추출된 두 가지 종류의 물성 을 적용하였다. 각 모델별로 계산된 수신감도 결과 를 Fig. 4(c)에 비교하였다. 가로축 주파수는 최대 운 용주파수 를 기준으로 정규화를 하였고, 세로축 수신감도는 full FEM에서 계산한 최저 주파수에서 의 값을 기준으로 하여 상대적인 값으로 나타내었 다. 검정색 점선은 등가 모델링 없이 전체 영역을 그 대로 해석했기 때문에 정확도 판단의 기준으로 삼았 고, S-D form 기반 추출 물성을 적용한 수신감도 결과 는 빨간색 실선으로, T-D form 기반 추출 물성을 적용 한 수신감도 결과는 파란색 실선으로 나타내었다. 두 가지 등가물성 추출 방식 모두 최대주파수의 3배 범위까지 full FEM 결과와 잘 일치하고 있으며, 오차 는 2 dB 이내로 나타났다. 1차 공진 모드는 모두 부근에서 발생했으며, 해당 모드 형상 또한 Fig. 4의 오른쪽 그림과 같이 가장자리부 굽힘의 형태로 동일 하게 모사되고 있다. 이로부터 경질 폴리머 재료를 적용한 압전복합재료 하이드로폰에 대해서는 S-D form, T-D form 두 가지 방법 모두가 등가물성을 비교 적 정확하게 추출한 것으로 판단할 수 있다. III. 연질 폴리머 압전복합재료의 등가물성 추출 및 정확도 분석 압전복합재료 기질에 적용된 폴리머의 종류에 따 라 하이드로폰의 수신 특성은 달라질 수 있다. 우레 탄 계열의 연질 폴리머는 II장의 에폭시 계열 경질 폴 리머에 비해 탄성계수가 낮고 프아송 비는 높은 특 징을 갖는다. II장과 동일한 형상의 하이드로폰 소자 에 대하여 Table 1에 기술된 물성을 갖는 연질 폴리머 를 적용하였을 때 등가물성을 계산하였고, 전반적인 추출 과정은 II장에 기술한 것과 동일하다. S-D form 기반의 추출 기법을 적용하였을 때, 입력조건에 대한 1방향과 3방향의 변형량 분포를 각각 Fig. 5(a), (b)에 나타내었다. 동일한 조건에서 경질 폴리머의 결과(Fig. 2)와 비교해보면, 평판 및 세라믹 대비 연질 폴리머의 낮은 탄성계수로 인해 의 분포가 외곽부 에 집중되고, 높은 포아송 비로 인해 측면인 2방향으 로의 변위가 비교적 크게 나타남을 확인할 수 있다. 다음으로 T-D form 기반 추출 방식에서 입력조 건에 대한 추출 결과를 Fig. 6에 나타냈다. 1방향 표면 에서의 1방향 응력 의 분포인 Fig. 6(a)를 보면, 연 질 폴리머의 낮은 탄성계수로 인해 동일한 조건의 경질 폴리머 결과[Fig. 3(c)]에 비하여 상하 평판 대비 폴리머에 작용하는 응력의 크기가 훨씬 작게 나타난 다. Fig. 6(b)의 3방향 응력 의 분포에서도 위치별 로 편차가 경질 폴리머[Fig. 3(d)]에 비해 증가한 것을 알 수 있다. 각 기법의 정확도를 판단하기에 앞서 폴리머 종류 에 따른 압전복합재료 등가물성과 특성을 비교하기 위해 S-D form과 T-D form 기반으로 추출된 최종 등 가물성 행렬을 각각 Tables 3과 4에 폴리머 종류별로 비교하였다. 먼저 폴리머 재료의 기계적 특성 차이 에 기인하여 등가물성이 상이하게 나타난 것을 확인 할 수 있다. 전반적으로 연질 폴리머를 적용하였을 Fig. 5. (Color available online) An example of esti- mating properties of a piezocomposite with soft polymer based on the S-D form of constitutive equations. Fig. 6. (Color available online) An example of esti- mating properties of a piezocomposite with soft polymer based on the T-D form of constitutive equations.1-3형 압전복합재료 등가물성 직접 추출 기법의 정확도 분석 The Journal of the Acoustical Society of Korea Vol.42, No.5 (2023) 383 때 강성 성분들은 낮게, 압전상수 성분들의 절대값 은 더 높게 도출되었다. 우선 경질 폴리머보다 연질 폴리머에서 물성행렬 s E의 대각 성분( , )이 더욱 높게 나타났는데, 이 는 측압이 없는 조건 하에서 각 축방향의 기계적 컴 플라이언스이므로 폴리머의 탄성계수 특성이 그대 로 반영된 것이라고 볼 수 있다. 이러한 양상은 3방향 보다 측면방향(1, 2방향)에서 더욱 두드러지는데, 그 원인으로는 3방향에서는 폴리머 대비 훨씬 단단한 세라믹의 특성이 지배적으로 나타나는 반면에 측면 방향에서는 폴리머의 영향성이 더욱 커지기 때문이 라고 할 수 있다. 한편 비대각 성분들을 보면, 값 대비 값의 크기가 경질 폴리머보다 연질 폴리머 에서 매우 크게 산출되었다. 이는 1방향 응력 작용 시 1방향 변형량 대비 2방향 변형량의 비율이므로 폴리 머의 포아송 비 특성을 반영한다. 즉 연질 폴리머의 포아송 비가 높기 때문에 발생하는 현상이라고 볼 수 있다. 3방향 거동이 포함되는 비대각 성분 의 경우에는 이러한 차이가 감소하는데, 이는 대각 성 분에서와 마찬가지로 단단한 세라믹의 특성이 3방 향에서 더 지배적으로 나타나기 때문이다. 압전상수 측면에서는 경질 폴리머에 비해 연질의 탄성계수가 낮아 값이 높게 나타났으나, 포아송 비가 크기 때 문에 의 절대값이 같이 높아짐을 확인할 수 있다. 연질 폴리머 적용 압전복합재료에 대한 등가물성 추출 정확도를 판단하기 위해 추출된 물성으로 계산 한 자유음장 수신감도를 Fig. 7에 도시하였다. 수신 감도는 Fig. 4(a), (b)와 동일한 유한요소 모델을 통해 계산되었다. 검정색 점선은 full FEM, 빨간색 실선은 S-D form, 파란색 실선은 T-D form 기반으로 추출된 수신감도를 의미하며, 주파수와 상대감도의 기준은 Fig. 4에서와 동일한 값을 사용하였다. 연질 폴리머 는 경질폴리머와는 달리 S-D form 기반 등가물성으 로 계산된 수신감도만이 full FEM 결과와 약 2 dB 오 차 이내로 유사하며 그래프 우측에 도시된 것과 같 이 1차 공진 모드의 주파수 및 모드 형상까지 비교적 잘 모사되고 있다. 반면 T-D form 기반 추출 물성은 전 주파수 대역에서 차이가 많이 발생하고 있으며 Table 3. Equivalent material properties of a 1-3 piezocomposite estimated using the S-D form depending on polymer materials. s E (pm 2 /N), d (pC/N), T Hard polymer Soft polymer Table 4. Equivalent material properties of a 1-3 piezocomposite estimated using the T-D form depending on polymer materials. c E (GPa), e (C/m 2 ), S Hard polymer Soft polymer Fig. 7. (Color available online) RVS and first mode shapes of a piezocomposite hydrophone with soft polymer calculated by a full FEM and equivalent models using directly estimated properties.Next >