Development of a Formula for Acoustic Transmission Loss According to Partition Positions in a Simple Expansion Chamber
Development of a Formula for Acoustic Transmission Loss According to Partition Positions in a Simple Expansion Chamber
Transactions of the Korean Society for Noise and Vibration Engineering. 2015. Feb, 25(2): 116-123
• Received : December 12, 2014
• Accepted : January 27, 2015
• Published : February 20, 2015
PDF
e-PUB
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
정 필, 박
Graduate School of Mechanical Engineering, Pusan National University
경 원, 배
Graduate School of Mechanical Engineering, Pusan National University
의 봉, 정
Corresponding Author ; Member, School of Mechanical Engineering, Pusan National University E-mail :wbjeong@pusan.ac.krTel : +82-51-510-2337, Fax : +82-51-517-3805
세 진, 안
Memver, School of Mechanical&Electrical Engineering, Uiduk University

Abstract
This paper proposed a new formula for predicting acoustic power transmission loss of simple expansion chamber with two internal partitions. Seventeen cases of computational results were performed and the database was constructed for the TL according to the various positions of internal partitions. Using this database, firstly, the formula for the peak values of TL at certain frequencies was developed using the least square estimation. Secondly, the formula for the TL curve could be obtained automatically with the input data of the positions of two internal partitions. The formula of TL developed in this paper showed good agreement with computational results. This formula will be helpful for the positioning of internal partitions to improve TL at target frequencies.
Keywords
1. 서 론
자동차, 중장비 및 건설 기계 등에서 발생하는 배기소음을 줄이기 위해서, 배기계에 소음기(silencer)를 설치한다.
소음기의 성능을 평가하는 방법은 대표적으로 삽입손실(insertion loss)과 투과손실(transmission loss)이 있다 (1) . 삽입손실은 음원과 소음기의 특성을 동시에 고려한 평가 방법이고, 투과손실은 소음기만의 고유한 특성을 평가하는 방법이다. 소음기의 투과손실 성능에 관한 연구로 Kim (2) 이 소음기의 입·출구 위치와 단면의 크기가 투과손실에 미치는 영향을 조사하였으며, Park (3) 이 내부 파티션의 위치와 개수가 투과손실에 미치는 영향을 조사하였다. 그리고 Kwon (4) 은 유동속도가 투과손실에 미치는 영향을 조사하였다. 그러나 3차원 음향전산해석을 통한 머플러의 투과손실 예측은 모델마다 새로운 메시(mesh)를 생성하여 해석을 해야 하기 때문에 많은 시간이 필요하므로 간편한 계산식에 의한 효율적이고 빠른 투과손실 예측법이 필요하다.
이 연구에서는 소음기의 특성으로 투과손실을 사용하여 소음기의 성능을 평가하고자 한다. Lee (5) 가 연구한 건설장비용 머플러를 이 논문의 머플러 모델 외형 크기로 정하였으며, 소음기의 전체 길이와 단면적은 일정하게 유지하는 것으로 가정한다. 머플러의 내부는 단순확장관(simple expansion chamber)으로 단순화하고 2개의 내부 파티션이 삽입되어 있다. 그리고 내부 파티션 위치에 따른 투과손실 해석결과 들을 데이터베이스화시켜 내부 파티션 위치에 따른 투과손실을 빠르고 효율적으로 추정할 수 있는 계산식을 도출하는 것이 목적이다.
2. 소음기 성능 평가
- 2.1 음향파워 투과손실
소음기의 성능을 측정하기 위한 구성도는 Fig. 1 과 같다. 입구에 x 12 의 간격으로 2개의 마이크로폰을 설치하고 출구에는 1개의 마이크로폰을 설치한다. 소음기의 출구에 무반사단을 연결한다.
PPT Slide
Lager Image
Model concept
3개의 지점에 설치된 마이크로폰에서 측정된 음압의 파워스펙트럼을 S 11 , S 22 , S 33 이라 하고, 1번과 2번 위치에서 측정된 음압의 상호스펙트럼을 S 12 = C 12 + j D 12 라 한다. 이 값들을 이용해 진행파의 파워스펙트럼 SAA 와 반사파의 파워스펙트럼 SBB 을 다음 식으로 나타낼 수 있다 (6,7) .
PPT Slide
Lager Image
PPT Slide
Lager Image
출구단에 무반사단을 연결하였으므로 3번 마이크로폰에서 측정한 음압의 파워스펙트럼 S 33 은 투과음의 파워스펙트럼이다.
입사파의 음향파워 Wi 와 투과파의 음향파워 Wt 는 다음 식으로 표현된다.
PPT Slide
Lager Image
여기서 pi 는 입사파음압, pt 는 투과파음압이고, Ai Ao 는 입출구 단면적이며 𝜌 c 는 공기의 임피던스이다.
음향파워 투과손실(TL)은 다음과 같이 표현할 수 있다.
PPT Slide
Lager Image
- 2.2 실험과 해석의 검증
이 논문에서 실험 및 해석을 위한 모델의 제원은 다음과 같다. 단순확장관의 길이 L =732 mm, 확장관 지름 D =394 mm, 입·출구단의 지름 d =112.5 mm. 내부에 삽입된 2개의 파티션의 구멍의 지름은 입·출구의 지름과 동일한 112.5 mm이다. Fig. 1 과 동일하게 입구에는 0.2 m의 간격을 가지는 2개의 마이크로폰을 설치하고 출구에는 1개의 마이크로폰을 설치한다. 입구단에 설치된 스피커에서는 백색잡음을 생성하고 출구단에는 무반사단을 적용하였다. Fig. 2 는 머플러의 투과손실 성능을 측정하는 실험 사진이다.
PPT Slide
Lager Image
Experimental setup for the measurement of TL
전산해석은 실험과 동일한 조건으로 진행하였다. 음향해석 소프트웨어인 LMS. Sysnoise를 사용하여 전산 해석을 진행하였으며, Fig. 3 은 model 2:2:4의 모델링을 나타낸다. 모델의 mesh 개수는 약 18만개이다.
PPT Slide
Lager Image
Modeling of the model 2:2:4 for computer analysis
신뢰성을 검증하기 위해 실험과 전산해석을 비교한 결과를 Fig. 4 에 나타내었다. 전산해석을 통한 투과손실 결과가 실험의 경향을 잘 표현하고 있으므로 이 논문에서는 해석한 투과손실 결과로 진행하였다.
PPT Slide
Lager Image
TL comparison with experiment and computer analysis(model 2:2:4)
3. 투과손실 예상 곡선 제작 원리
- 3.1 해석 모델
건설장비용 머플러를 Fig. 5 와 같이 단순화하였다. 모델의 길이 및 지름은 2.2절의 모델과 동일하며, 단순확장관의 입구단 공간의 길이 L 1 =x·L과 가운데 공간의 길이 L 2 =y·L와 출구단 공간의 길이 L 3 =(1-x-y)·L는 길이 L과의 길이비로 정의한다. 그리고 파티션 길이 L 1 , L 2 에 따른 투과손실 특성을 분석하였다.
PPT Slide
Lager Image
Model geometry
이 논문에서는 길이비율 L 1 : L 2 : L 3 로 모델을 구별하였다. 해석 결과를 비교한 결과, 길이 L 1 L 3 이 바뀐 모델의 투과손실 성능이 같다는 것을 확인하였다. 예를 들면 3:3:2 모델과 2:3:3의 모델의 투과손실 성능은 같다. 투과손실 데이터베이스를 구성하기 위해서 17가지의 모델(1:1:6, 1:2:5, 1:3:4, 1:4:3, 1:5:2, 1:6:1, 2:1:5, 2:2:4, 2:3:3, 2:4:2, 3:1:4, 3:2:3, 3:3:2, 3:4:1, 4:1:3, 4:2:2, 4:3:1)에 대하여 전산 해석을 수행하였다.
- 3.2 투과손실 해석
파티션 위치 변경에 따른 17가지 모델 중 몇 개의 모델에 대한 투과손실 해석 결과를 Fig. 6 에 나타내었다. 여기서 평면파가 성립하는 차단주파수, fc 를 구하면 다음과 같다.
PPT Slide
Lager Image
PPT Slide
Lager Image
Numerical results of TL according to partition position
Fig. 6 에서 차단주파수 이내의 주파수 범위에서 TL특성 형태가 일관성이 있음을 알 수 있다. 따라서 이 논문에서는 2번째 최대 TL값을 포함하는 300 Hz까지를 해석주파수 범위로 설정하였다.
L 1 : L 2 : L 3 =1:2:5인 모델에 대한 투과손실 해석 예를 Fig. 7 에 나타낸다. Fig. 7 에서 투과손실의 최대값을 가지는 점을 ( f 1 , ɡ 1 )와 ( f 3 , ɡ 3 )라 한다. 또한, 최소값의 투과손실을 가지는 점을 ( f 2 , ɡ 2 ), ( f 4 , ɡ 4 )라 한다. 파티션 위치를 다양하게 변화시킨 17종류의 모델에 대한 투과손실값 ɡ ij = ɡ i ( x j , 𝓎 j ), j =1~17을 Table 1 에 나타냈으며, 대응하는 주파수 fij = fi ( x j , 𝓎 j )를 Table 2 에 나타냈다. 여기서 i Fig. 7 에 나타난 점의 위치, j 는 17종류의 모델번호를 나타낸다.
PPT Slide
Lager Image
Symbols representing max and min points in TL curve
Value ofTLmaxandTLminaccording to the partition locations
PPT Slide
Lager Image
Value of TLmax and TLmin according to the partition locations
Frequencies at whichTLmaxorTLminoccur according to the partition locations
PPT Slide
Lager Image
Frequencies at which TLmax or TLmin occur according to the partition locations
- 3.3 내부 파티션과 투과손실의 관계
단순확장관의 전체 길이( L )와 단면적은 변하지 않는다고 가정하였으며, 투과손실은 앞부분 공간의 길이의 비( x )와 가운데 공간의 길이의 비( y )의 함수로 가정한다.
내부 파티션의 위치와 투과손실의 관계를 파악하기 위해, 여러 개의 모델 중 j 번째 모델의 길이 비인 변수(𝑥,𝓎)에 의해 투과손실 ɡi (𝑥,𝓎)과 주파수 fi (𝑥,𝓎)가 변하는 특징을 함수화시켜 다음 식으로 표현했다.
PPT Slide
Lager Image
PPT Slide
Lager Image
Table 1 Table 2 에서 제시한 17종류의 모델의 길이 비 (𝑥 j ,𝓎 j ) 와 ɡi (𝑥 j ,𝓎 j ), fi (𝑥 j ,𝓎 j ) , j = 1 ∼ 17를 이용하여 식 (6), (7)을 행렬식으로 표현하면 다음 식 (8)과 같다.
PPT Slide
Lager Image
단,
PPT Slide
Lager Image
여기서 17개의 모델에 대한 파티션 위치 행렬 [ T ]는 [17×10], 미지수 벡터 { Ci }와 { Di }는 [10×1], 벡터 { Gi }와 { Fi }는 ɡij = ɡi (𝑥 j ,𝓎 j ) 와 fij = fi (𝑥 j ,𝓎 j )를 나타내는 [17×1] 벡터이다.
예측 오차의 제곱식은 다음과 같다.
PPT Slide
Lager Image
식 (10)의 오차제곱을 최소화하는 계수 { Ci }와 { Di }는 다음과 같다.
PPT Slide
Lager Image
최소제곱법을 적용하여 구한 식 (11)의 미지수 벡터 { Ci }와 { Di }는 Table 3 Table 4 와 같다. 미지수 벡터 { Ci }, { Di }, i =1~4를 구하였으므로, 변수( x , y )와 ɡi , fi 의 관계를 함수화하여 특징을 파악하였다고 볼 수 있다.
Coefficients of polynomials for determining value of TL given inTable 1
PPT Slide
Lager Image
Coefficients of polynomials for determining value of TL given in Table 1
Coefficients of polynomials for determining frequencies given inTable 2
PPT Slide
Lager Image
Coefficients of polynomials for determining frequencies given in Table 2
임의의 길이 비 x = L 1 /L, y = L 2 /L값과 Table 3 Table 4 에서 구한 계수를 이용하여 식 (6)과 식 (7)에 대입하면, 길이 비 ( x , y )에 따른 투과손실의 최대 및 최소값과 그 때의 주파수 ( ɡi (𝑥,𝓎), fi (𝑥,𝓎), i =1~4)를 예측할 수 있다.
- 3.4 투과손실 곡선 계산과정
앞의 과정을 통해 전체길이( L )에 대한 파티션의 길이 비율인 x y 값을 결정하면 투과손실 곡선의 최대 및 최소값과 그 때의 주파수를 예측할 수 있다. 투과손실 성능 예상 곡선을 4개의 구간으로 나누어 다음 식 (12)와 같이 가정하였다. 구간에 대한 정보는 Table 5 에 나타내었다.
Information of value in frequency range in Eq. (12)
PPT Slide
Lager Image
Information of value in frequency range in Eq. (12)
PPT Slide
Lager Image
단,
PPT Slide
Lager Image
여기서 주파수 f = fp 에서 TL은 최대값 ɡp 를 가지며, 주파수 f = fq 에서 TL은 최소값 ɡq 를 갖는다. 식 (12)에서 좌변과 우변 오차의 log스케일을 주파수 fi 마다 나타내면 다음과 같다.
PPT Slide
Lager Image
또는 간단히 표현하면 다음과 같다.
PPT Slide
Lager Image
예측 오차의 제곱을 최소화하도록 최소제곱법을 적용하면 식 (12)의 차수 n 은 다음과 같이 구해진다.
PPT Slide
Lager Image
이러한 방법으로 17개에 모델의 구간마다의 차수를 구하였으며, Table 6 에 나타내었다. 여기서 { nij }는 주파수 i 구간에서 j 번째 모델에 대한 차수를 의미한다.
Value of the order according to the partition locations
PPT Slide
Lager Image
Value of the order according to the partition locations
내부 파티션의 위치( x,y )에 따른 차수 n 을 계산할 수 있는 보간식은 다음과 같이 표현하였다.
PPT Slide
Lager Image
17종류의 모델에 대하여 식 (16)을 행렬식으로 표현하면 다음과 같다.
PPT Slide
Lager Image
단,
• {Ni}T= [N1iN2i⋯N10i]
• {ni}T= [ni1ni2⋯ni17]
PPT Slide
Lager Image
예측 오차의 제곱을 최소화하도록 최소제곱법을 적용하면 계수벡터 { Ni }는 다음과 같이 표현된다.
PPT Slide
Lager Image
식 (18)에서 구한 미지수 벡터 { Ni }을 Table 7 에 나타내었다. 이상과 같이 파티션 위치변수 ( x , y )가 주어지면 식 (12)에서 TL곡선을 자동으로 예측할 수 있게 되었다.
Coefficients of polynomials for determining value of order given in Eq. (16)
PPT Slide
Lager Image
Coefficients of polynomials for determining value of order given in Eq. (16)
4. 결과 및 비교
- 4.1 GUI 개발
내부에 두 개의 파티션이 있는 단순확장관의 모델들을 해석하여 얻은 결과를 바탕으로 내부 파티션 길이의 비만을 이용하여 투과손실의 최대 및 최소점( f 1 , ɡ 1 ), ( f 2 , ɡ 2 ), ( f 3 , ɡ 3 ), ( f 4 , ɡ 4 )을 예상하는 계산식을 개발하였다. 그리고 원점, 최대 및 최소점 중 이웃한 2점을 잇는 곡선의 방정식도 개발하였다. 이로써 컴퓨터 해석을 하지 않고 내부 파티션 길이 비를 이용하여 투과손실 성능을 예측할 수 있는 계산식을 개발하였다. 이 논문에서 개발한 계산식을 바탕으로 만든 GUI를 이용해 내부 파티션의 길이 비 𝑥 = 2/8, 𝓎 = 1/8 인 모델 2:1:5의 투과손실을 예측한 결과를 Fig. 8 에 제시하였다.
PPT Slide
Lager Image
GUI for the estimation of TL curve(2:1:5 model)
- 4.2 투과손실 해석과 계산식의 비교
모델 2:1:5를 컴퓨터 해석으로 구한 투과손실과 내부 파티션 길이 비를 이용해 계산식으로부터 구한 투과손실을 Fig. 9 에 비교하여 나타내었다. 모델 2:1:5는 계산식을 개발할 때 기초자료로 사용되었던 모델이다.
PPT Slide
Lager Image
Comparison of a typical TL curve used for the database(model 2:1:5)
계산식을 개발할 때 사용하지 않은 모델 1:1:1에 대해서, 컴퓨터 해석으로 구한 투과손실과 내부 파티션 길이 비를 이용하여 계산식으로부터 구한 투과손실의 비교를 Fig. 10 에 나타내었다.
PPT Slide
Lager Image
Validation of TL curve by the proposed formula for arbitrary position of partition(model 1:1:1)
Fig. 9 Fig. 10 의 결과를 비교하면 특징을 잘 따라가며 거의 일치하는 것을 볼 수 있다. 따라서 해석을 하지 않고 모델의 내부 파티션의 위치에 따른 길이 비만을 이용해 투과손실 성능을 추정할 수 있는 새로운 계산식이 개발되었음을 확인할 수 있다.
5. 결 론
이 논문에서는 전산해석에 의한 17가지 경우에 대한 투과손실특성을 기초로 최소제곱 추정법을 이용하여, 단순확장관 내부에 있는 2개의 파티션 위치변화에 따른 투과손실 측정곡선을 자동으로 예측할 수 있는 계산식을 개발하였다. 개발된 계산식을 이용하면 전산해석을 수행하지 않아도 투과손실을 정확히 예측할 수 있음을 확인하였으며, 향후 머플러 설계에 필요한 시간과 비용을 줄이는데 도움이 될 것으로 기대된다.
A part of this paper was presented at the KSNVE 2014 Annual Autumn Conference
BIO
JeongPil Park received B.S. degrees from in Pusan National University in 2014 and he is currently majoring in noise and vibration in Pusan National University. He especially is interested in Coil Spring.
KyeongWon Bae received B.S. degrees from in Pusan National University in 2014 and he is currently majoring in noise and vibration in Pusan National University. He especially is interested in Modal Expansion Method.
Weuibong Jeong received B.S. and M.S. degrees from Seoul National University in 1978 and from KAIST in 1980, respectively. He then received his Ph.D. degree from Tokyo Institute of Technology in 1990. Dr. Jeong is currently a professor at the department of mechanical engineering at Pusan National University in Busan, Korea. His research interests are in the area of the measurement and signal processing, finite/boundary element analysis of noise and vibration, fluid-structure interactions and acoustic-structure interactions.
SeJin Ahn is received B.S., M.S. and Ph.D. degrees from Pusan National University in 1994, 1996, and 2003, respectively. Dr. Ahn is currently a professor at the department of energy & electricity in Uiduk University. His research interest is in human vibration.
References
Munjal M. L. 1987 Acoustics of Ducts and Mufflers Wiley-Interscience New York
Kim B. J. , Jeong W. B. , Hwang S. M. 1999 Acoustic Analysis of Circular Simple Expansion Chamber with Arbitrary Location and Cross-section Area of Inlet/Outlet Transactions of the Korean Society for Noise and Vibration Engineering 9 (5) 914 - 921
Park J. P. , Bae K. W. , Jeong W. B. , Ahn S. J. 2014 Analysis of Transmission Loss according to Internal Space Partition of Simple Expansion Chamber Proceedings of the KSNVE Annual Spring Conference 520 - 521
Kwon J. , Jeong W. B. , Hong J. S. 2012 Effects of Fluid Velocity on Acoustic Transmission Loss of Simple Expansion Chamber Transactions of the Korean Society for Noise and Vibration Engineering 22 (10) 994 - 1002    DOI : 10.5050/KSNVE.2012.22.10.994
Lee S. B. , Lee W. T. , Nam K. H. 2009 Development of a Muffler for the Super Heavymachinery in the 70 Ton Class Journal of the Korea Academia-Industrial Cooperation Society 10 (10) 2557 - 2564    DOI : 10.5762/KAIS.2009.10.10.2557
Tao Z. , Seybert A. F. 2003 A Review of Current Techniques for Measuring Muffler Transmission Loss
Seybert A. F. 1988 Two-sensor Methods for the Measurement of Sound Intensity and Acoustic Properties in Ducts The Journal of the Acoustical Society of America 83 (6) 2233 - 2239    DOI : 10.1121/1.396352