Advanced
Wave Transmission Analysis of Semi-infinite Mindlin Plates Coupled at an Arbitrary Angle
Wave Transmission Analysis of Semi-infinite Mindlin Plates Coupled at an Arbitrary Angle
Transactions of the Korean Society for Noise and Vibration Engineering. 2014. Dec, 24(12): 999-1006
Copyright © 2014, The Korean Society for Noise and Vibration Engineering
  • Received : October 28, 2014
  • Accepted : December 03, 2014
  • Published : December 20, 2014
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
영호 박
Corresponding Author ; Member, Department of Naval Architecture and Marine Engineering, Changwon National University E-mail :parkyh@changwon.ac.krTel : +82-55-213-3684, Fax : +82-55-213-3689

Abstract
Mindlin plate theory includes the shear deformation and rotatory inertia effects which cannot be negligible as exciting frequency increases. The statistical methods such as energy flow analysis(EFA) and statistical energy analysis(SEA) are very useful for estimation of structure-borne sound of various built-up structures. For the reliable vibrational analysis of built-up structures at high frequencies, the energy transfer relationship between out-of-plane waves and in-plane waves exist in Mindlin plates coupled at arbitrary angles must be derived. In this paper, the new wave transmission analysis is successfully performed for various energy analyses of Mindlin plates coupled at arbitrary angles.
Keywords
1. 서 론
산업이 발달하고 삶의 질에 대한 관심이 증대함에 따라 고주파수 대역에서의 소음진동현상에 중요성이 증대되고 있다. 하지만 고주파수 대역에서는 시스템 자체 물리적 파라미터(physical parameters)의 불확실성과 모드 밀도(modal density)가 높아지기 때문에 기존 저주파수 대역에서의 유한요소법(FEM)이나 경계요소법(BEM)은 예측 결과의 타당성 측면에서나 해석 비용 측면에서 적절하지 않은 것으로 알려져 있다. 이런 이유로 고주파수 대역에서는 에너지에 기반한 통계적에너지해석법(SEA)이나 에너지흐름해석법(EFA)과 같은 통계적 방법(statistical approach)이 대안으로 여겨지고 있다 (1~3) .
에너지 변수 기반의 통계적에너지해석법이나 에너지흐름해석법은 복합구조물(built-up structures)의 진동에너지를 예측하기 위해 불연속 연성구조물 간의 에너지전달관계(energy transfer relationship) 확립이 요구된다. 일반적으로 에너지 기반 통계적 기법에서 약한 연성(weakly coupling) 연성구조물 간의 주파수 평균된 분산장(diffuse field)에서의 에너지전달관계는 반무한장(semi-infinite wavefield)에서의 파워투과반사계수(power transmission and reflection coefficient)로 표현될 수 있다 (3) .
주파수가 높아짐에 따라 대표적인 구조요소인 보(beam)와 평판(plate)의 면외운동(out-of-plane motion)에서 회전관성효과(rotatory inertia effect)와 전단변형효과(shear distortion effect)가 매우 중요해진다 (4~9) . 따라서 고주파수 대역에서 고체음(structure borne noise) 전달 예측에 효과적인 통계적에너지해석법이나 에너지흐름해석법은 보다 신뢰성 있는 소음진동응답의 예측을 위해 이와 같은 효과를 반영할 수 있는 Timoshenko 보와 Mindlin 판에 대한 에너지 모델 적용이 반드시 요구된다. 박영호와 홍석윤은 (6~8) Timoshenko 보와 Mindlin 판의 에너지흐름모델을 개발했고 Timoshenko 연성 구조물의 에너지흐름해석을 위해 파동전달해석을 수행하였다. 또한 연성된 Mindlin 판의 에너지전달관계를 얻기 위해 박영호는 (9) 동일 평면상에 존재하는 Mindlin 판의 파동전달해석(wave transmission analysis)을 수행하였다.
이 논문에서는 고주파수 대역에서 일반적인 형상의 구조물의 신뢰성 있는 소음진동해석을 위해 필수적인 임의의 각으로 연성된 Mindlin 판의 파동전달 해석을 수행하였다. 이를 통해 에너지 기반 통계적 해석기법에서 복합구조물의 보다 신뢰성 있는 소음 진동응답 예측이 가능할 것으로 판단된다.
2. Mindlin 판 이론
- 2.1 Mindlin 판의 면외 운동방정식
Mindlin 판의 면외 운동방정식(equation of out-of-plane motion)은 다음과 같이 표현된다.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서 w 은 횡방향 변위를 나타내며, αx αy 는 굽힘에 의한 각 축방향의 회전각, h 는 판의 두께, I = h 3 /12는 단위 폭 당 관성모멘트, D = Eh 3 /[12(1− v 2 )]은 굽힘 강성률, ν 는 푸아송비, E 는 탄성계수, κ = 20(1+ v )/(24+25 v + v 2 )는 전단보정계수, G = E /[2(1+ v )]는 전단계수를 나타낸다.
- 2.2 Mindlin 판의 면외운동
Mindlin 판의 면외 운동방정식은 변위 포텐셜 함수(displacement potential function) ϕ ψ 을 정의하여 연성된 변위의 파동성분을 분리하여 표현할 수 있다 (5) .
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서 회전각과 변위 포텐셜과의 관계는 다음과 같다.
PPT Slide
Lager Image
변위 포텐셜 함수와 수직 변위로 표현된 운동방정식 (4~6)을 이용하면 3가지 종류의 파동해가 얻어진다. 굽힘우세횡파(BDFW: bending dominant flexural wave)는 Kirchhoff 판의 굽힘파와 같이 전 주파수영역에서 원거리 파동해(farfield wave solution)의 특성을 갖고 있으며 면외전단파(OPSW: out-of-plane shear wave)와 전단우세횡파(SDFW: shear dominant flexural wave)는 임계주파수
PPT Slide
Lager Image
위에서만 원거리 해가 존재한다 (8) .
각 파동의 파수(wavenumber)는 다음과 같이 표현될 수 있다.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서 k 1 k 2 , k 3 은 각각 면외전단파, 굽힘우세횡파, 전단우세횡파의 파수를 나타낸다.
- 2.3 평판의 면내 운동방정식과 면내운동
임의의 각으로 연성된 평판의 진동예측을 위해서는 평판 내의 면내운동에 대한 고려가 필요하다. 2차원 평판의 면내 운동방정식은 다음과 같이 서로 직교방향의 면내운동 u v 의 연성식으로 표현된다.
PPT Slide
Lager Image
PPT Slide
Lager Image
면내운동은 Mindlin 판의 면외운동과 유사하게 면내 변위
PPT Slide
Lager Image
를 종파의 특성을 가진 스칼라 변위 포텔셜함수( μ )와 면내전단파의 특성을 갖은 벡터 변위 포텐셜함수
PPT Slide
Lager Image
를 이용하면 면내 운동방정식 (12,13)은 종파와 면내전단파의 파동방정식으로 분리될 수 있다 (9) .
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서
PPT Slide
Lager Image
이고
PPT Slide
Lager Image
이다.
3. 임의의 각으로 연성된 두 개의 Mindlin 판의 파동전달해석
- 3.1 Mindlin 판의 진동 파워
Mindlin 판의 면외진동(out-of-plane vibration)과 면내진동(in-plane vibration)에서 x 방향의 시간평균된파워는 다음과 같은 식으로 나타낼 수 있다.
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서 K = Eh /(1 − v 2 )이고 *는 켤레복소수를 의미한다.
- 3.2 Mindlin 판의 면외/면내 파동
임의의 각으로 연성된 Mindlin 판에서는 모든 평판에서 면외 파동과 면내 파동이 동시에 존재한다. 먼저 입사 평판의 면외 파동을 살펴보면 다음과 같다. Fig. 1 의 입사 평판(incident plate) i 에서 임계주파수보다 높은 주파수 영역일 경우
PPT Slide
Lager Image
의 크기를 가진 원거리 면외전단파(far-field OPSW)가 두 판의 경계로 θi 의 각도로 입사할 때 입사 평판 i 에 존재하는 면외전단파는 다음과 같이 표현할 수 있다.
PPT Slide
Lager Image
Two semi-infinite Mindlin plates coupled at an arbitrary angle
PPT Slide
Lager Image
여기서 θi = θr , k 1x,i = k 1,i cos θ i , k 1y,i = k 1,i sin θi 이다.
마찬가지로 입사 평판 i 에서
PPT Slide
Lager Image
의 크기로 입사하고 D 의 크기로 반사하는 원거리 굽힘우세횡파는 다음과 같다.
PPT Slide
Lager Image
여기서 k 2x,i = k 2,i cos θ i , k 2y,i = k 2,i sin θi 이다.
또한 입사 평판 i 에서, 임계주파수보다 높은 주파수 영역일 경우
PPT Slide
Lager Image
의 크기로 입사하고 F 의 크기로 반사하는 원거리 전단우세횡파는 다음과 같다.
PPT Slide
Lager Image
여기서 k 3x,i = k 3,i cos θ i , k 3y,i = k 3,i sin θi 이다.
입사 평판 i 에서 면내 파동(in-plane wave)은 면내 종파(in-plane longitudinal wave)와 면내 전단파(in-plane shear wave)로 구분되며 평판의 면내 변위 u , v 를 각 파동으로 표현하면 다음과 같다.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서
PPT Slide
Lager Image
,
PPT Slide
Lager Image
는 각각 입사 면내 종파와 면내 전단파의 크기이며 H , J 는 면내 종파와 면내 전단파의 반사파의 크기이다. 또한 cos θli = k lx,i / k l,i , cos θsi = k sx,i / k s,i , sin θli = k ly,i / k l,i , sin θsi = k sy,i / k s,i 이다.
입사 평판 i 에서 입사한 파동의 투과로 인해 투과 평판 t 에서 존재하는 면외, 면내 투과파는 다음과 같이 표현된다.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서 경계에서 y 방향 파수 매칭(matching)과 전반사 조건을 고려하면
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
이다.
Mindlin 판이 두 개가 임의의 각으로 연성이 되면 식 (19)~(32)의 식과 같이 각 평판에서 5개의 미지수가 존재한다. 따라서 연성된 두 개의 Mindlin판에서의 파동전달해석을 위해서는 총 10개의 경계조건이 필요하다. 두 평판에서의 경계에서 변위 연속조건이 5개가 존재하며, 힘-모멘트 연속조건이 5개 존재한다. 두 판의 경계( x = xt = 0)에서의 경계조건을 식으로 표현하면 다음과 같다.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서
PPT Slide
Lager Image
이다.
위와 같은 10개의 경계조건을 적용하면 파동장 식 (19)~(32)의 10개의 미지수를 얻을 수 있고 각 파동의 시간평균된 x 방향의 파워는 다음과 같이 표현된다.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
- 3.3 Mindlin 판의 분산장에서의 파워투과 반사계수
분산장에 대한 파워투과반사계수는 입사각에 따른 각 파동의 파워를 나타내는 식 (46)~(60)을 입사각에 대한 평균을 취하면 다음과 같은 분산장에 대한 파워투과반사계수가 유도된다 (7,9) .
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서 T mn12 , R mn11 은 각각 “ m ”파형의 파동이 입사할 경우 “ n ”파형의 파동으로 투과, 반사되는 파워 비율을 의미한다. 유도된 동일 평면상에서 연성된 Mindlin 판의 수치해석을 위해 Fig. 1 과 같은 반무한(semi-infinite) 평판구조물에 철(steel) 재질의 각각 0.02 m, 0.04 m 두께의 입사평판, 투과평판이 채택되었다. 각 평판의 임계주파수는 8.1 kHz, 4.1 kHz이다. 두 평판의 연성각도는 45°로 가정하였다.
임의의 각으로 연성된 Mindlin 판에 대한 파동전달해석결과가 Figs. 2 ~ 6 에 나타나 있다. 임의의 각으로 연성된 Mindlin 판에 존재하는 각 5종류의 파동이 입사할 경우에 주파수에 따른 파워투과반사계수를 도시하였다. 일반적으로 파워손실이 없는 보존 조인트(conservative joint)인 경우 에너지 보존으로 인해 파워투과반사계수의 합은 1이다. Figs. 2 ~ 6 에서 모든 파워투과반사계수들의 합은 모든 주파수 대역에서 정확히 1이 됨을 확인할 수 있고, 이를 통해 앞에서 전개한 이론식들의 타당성을 검증할 수 있다. 특히 Figs. 2 , 4 에서는 첫 번째 평판의 임계주파수(8.1 kHz) 아래 영역에서 모든 파동의 파워투과판사계수가 0이다. 이는 임계주파수 아래 영역에서는 입사파인 OPSW, SDFW가 진행파가 아니기 때문에 원거리 영역의 입사 파워가 정의될 수 없어 나타나는 현상이다. 이와 반대로 BDFW와 면내 종파, 면내 전단파가 입사하는 경우에는 전 주파수 영역에서 0이 아닌 파워투과반사계수가 존재한다. 물론 두 번째 평판의 임계주파수(4.1 kHz) 아래에서는 투과판의 OPSW와 SDFW와 연관되는 파워투과계수는 0이 됨을 확인할 수 있다.
PPT Slide
Lager Image
Power transmission and reflection coefficients of coupled coplanar Mindlin plates in incident OPSW
PPT Slide
Lager Image
Power transmission and reflection coefficients of coupled coplanar Mindlin plates in incident BDFW
PPT Slide
Lager Image
Power transmission and reflection coefficients of coupled coplanar Mindlin plates in incident SDFW
PPT Slide
Lager Image
Power transmission and reflection coefficients of coupled coplanar Mindlin plates in incident in-plane longitudinal wave
PPT Slide
Lager Image
Power transmission and reflection coefficients of coupled coplanar Mindlin plates in incident in-plane shear wave
4. 결 론
에너지흐름해석법(EFA)이나 통계적에너지해석법(SEA)과 같은 에너지 기반 통계적 접근법은 모드밀도가 높은 고주파수 대역 진동해석에 유용한 방법이다. 따라서 1, 2차원 구조 부재인 보와 평판으로 이루어진 복합구조물의 신뢰성 있는 고주파수 대역 진동에너지를 예측하기 위해서는, 고주파수 대역 진동응답에서 무시할 수 없는 회전관성효과와 전단변형 효과를 반영할 수 있는 Timoshenko 보와 Mindlin 판에 대한 에너지 기반 통계적 해석 기법의 개발이 필수적이다.
이 연구는 임의의 형상의 복합 Mindlin 평판구조물의 에너지 기반 진동 해석을 위해 반드시 필요한 Mindlin 판에서의 면외파와 면내파 사이의 에너지 전달 관계를 유도함으로써 복합 평판구조물의 고주파수 대역 진동 해석의 신뢰성을 높이는 데 기여할 것으로 예상된다.
기 호 설 명
ωc : Mindlin 판의 임계주파수 k1 : 면외전단파(OPSW)의 파수 k2 : 굽힘우세횡파(BDFW)의 파수 k3 : 전단우세횡파(SDFW)의 파수 kl : 면내 종파(in-plane longitudinal wave)의 파수 ks : 면내 전단파(in-plane shear wave)의 파수 ϕ,ψ : 면외운동 변위 포텐셜 함수 μ : 면내운동 스칼라 변위 포텐셜 함수 : 면내운동 벡터 변위 포텐셜 함수
Acknowledgements
이 논문은 2014~2015년도 창원대학교 연구비에 의하여 연구되었음.
BIO
Young-Ho Park was born in Seoul, Korea, in 1975. He studied Naval Architecture and Ocean Engineering at Seoul National University for a bachelor's degree in 1999. He received a Ph.D. in engineering of Seoul National University in 2006. He has 5 years' experience in automotive NVH engineering by working in Hyundai Motor Company. In 2011, he moved to Changwon National University, Changwon-si, Korea, where he became an assistant professor and continued research on the vibro-acoustics of transport machines in Department of Naval Architecture and Marine Engineering.
References
Lyon R. H. , Dejong R. G. 1995 Theory and Application of Statistical Energy Analysis Second ed. Butterworth-Heinemann London
Wohlever J. C. , Bernhard R. J. 1992 Mechanical Energy Flow Models of Rods and Beams Journal of Sound and Vibration 153 (1) 1 - 19    DOI : 10.1016/0022-460X(92)90623-6
Cho P. E. 1993 Energy Flow Analysis of Coupled Structures, Ph.D. Thesis Purdue University
Liew K. M. , Xiang Y. , Kitipornchai S. 1995 Research on Thick Plate Vibration: a Literature Survey Journal of Sound and Vibration 180 (1) 163 - 176    DOI : 10.1006/jsvi.1995.0072
Nelson H. M. 1978 High Frequency Flexural Vibration of Thick Rectangular Bars and Plates Journal of Sound and Vibration 60 (1) 101 - 118    DOI : 10.1016/0022-460X(78)90404-2
Park Y. H. , Hong S. Y. 2006 Vibrational Energy Flow Analysis of Corrected Flexural Waves in Timoshenko Beam - Part I: Theory of an Energetic Model Shock and Vibration 13 (3) 137 - 165    DOI : 10.1155/2006/308715
Park Y. H. , Hong S. Y. 2006 Vibrational Energy Flow Analysis of Corrected Flexural Waves in Timoshenko Beam - Part II: Application to Coupled Timoshenko Beams Shock and Vibration 13 (3) 167 - 196    DOI : 10.1155/2006/562762
Park Y. H. , Hong S. Y. 2008 Vibrational Power Flow Models for Transversely Vibrating Finite Mindlin Plate Journal of Sound and Vibration 317 (3-5) 800 - 840    DOI : 10.1016/j.jsv.2008.03.049
Park Y.H. 2013 Wave Transmission Analysis of Co-planar Coupled Semi-infinite Mindlin Plate Transactions of the Korean Society for Noise and Vibration Engineering 23 (6) 574 - 580    DOI : 10.5050/KSNVE.2013.23.6.574