Transactions of the Korean Society for Noise and Vibration Engineering
[ Article ]
Transactions of the Korean Society for Noise and Vibration Engineering - Vol. 28, No. 4, pp.410-416
ISSN: 1598-2785 (Print) 2287-5476 (Online)
Print publication date 20 Aug 2018
Received 24 Apr 2018 Revised 18 Jul 2018 Accepted 18 Jul 2018
DOI: https://doi.org/10.5050/KSNVE.2018.28.4.410

두 영역 분할법을 이용한 임의 형상 음향 공동의 고유치 해석

강상욱
Eigenvalue Analysis of Arbitrarily Shaped, Acoustic Cavities Using Two-domain Method
Sang Wook Kang

Correspondence to: Member, Hansung University E-mail: swkang@hansung.ac.kr Recommended by Editor Jong Kwan Ryu


ⓒ The Korean Society for Noise and Vibration Engineering

초록

A two-domain method for eigenvalue analysis of arbitrarily-shaped acoustic cavities with convex or concave shapes is proposed in this paper. This method divides the concave acoustic cavity of interest into two convex regions. The sound pressure (approximate solution) of each convex region is assumed by linearly superposing plane waves generated at edges of the region. A sub-system matrix for each convex region is extracted by applying a provisional boundary condition to the approximate solution. Finally, a system matrix, of which the determinant provides the eigenvalues of the concave cavity, is made by considering the rigid-wall boundary condition at edges and the compatibility condition (the condition of continuity in sound pressure and its slope) at the interface between the two regions. In addition, a practical method for reducing the size of the system matrix is presented to facilitate the calculation of the determinant of the system matrix. Case studies show that the proposed method is valid and accurate when the eigenvalues by the proposed method are compared to those by the exact method or FEM.

Keywords:

Eigenvalues, Two-domain Method, Plane Wave, Concave Acoustic Cavity, Arbitrarily Shaped Acoustic Cavity

키워드:

고유치, 두 영역 분할법, 평면파, 오목 음향 공동, 임의 형상 음향 공동

1. 서론

저자는 임의 형상 음향 공동의 고유치와 고유모드를 정확히 구하기 위한 방법으로 NDIF법(non-dimensional dynamic influence function method)을 2000년 처음으로 개발하였다(1). 그리고 최근에 오목 형상 음향 공동에 대해서 해의 정밀도가 많이 떨어지는 NDIF법의 단점을 극복하기 위해, 분할영역법 기반 NDIF법을 새로이 제안하였다(2).

NDIF법은 유한요소법(3)과 경계요소법(4)에 비해 해의 정밀도 면에서 탁월한 우수성을 가지고 있음이 저자의 선행 연구에서 증명되었다(1,2). 그러나, NDIF법은 유한요소법이나 경계요소법과는 달리 최종 시스템 행렬이 주파수 파라미터에 종속되는 단점을 가진다. 결과적으로 최종 시스템 행렬식은 대수 고유치 문제(algebraic eigenvalue problem)의 형태를 가지지 못하게 되어서, 시스템 행렬로부터 고유치와 고유모드를 추출하는 과정이 매우 번거롭게 된다(1,2).

저자는 이러한 문제점을 해결하기 위해, 최근에 기존 NDIF법 이론을 개선하여 임의 볼록 형상 음향 공동의 고유치와 고유모드를 수월하게 추출할 수 있는 개선된 NDIF법인 MNDIF법을 개발하였다(5). MNDIF법은 최종 시스템 행렬식이 NDIF법과는 차별적으로 대수 고유치 문제의 형태를 가지는 장점을 가진다. 나중에 상기의 방법을 확장하여, 오목 형상 음향 공동에 대한 MNDIF법을 정립하였다(6).

그러나, NDIF법과 MNDIF법은 해석 대상물의 경계를 이산화하기 위해 노드의 개수를 증가시킬 경우, 저차의 고유치가 수렴하지 않고 발산하는 단점을 가진다. 저자는 이러한 단점을 극복하기 위해, 평면파(plane wave) 형태의 시리즈 함수(series function)를 이용하여 임의 형상 평판과 멤브레인의 고유치와 고유모드를 추출하는 방법(7,8)에 대한 연구도 병행하여 왔으며, 이 논문에서는 기존 연구(7,8)를 확장하여 임의 형상 음향 공동에 대한 이론을 정립하고자 한다.

이 논문에서는 볼록 형상뿐만 아니라 임의 오목 형상 음향 공동에 대해서도 적용 가능한, 평면파를 이용한 영역 분할법(sub-domain method) 기반 음향 공동 고유치 해석 기법을 제안하고자 한다.


2. 이론 정식화

2.1 지배방정식과 경계조건

음향 공동의 지배방정식은 식 (1)과 같이 헤름홀츠 방정식(Helmholtz equation)의 형태를 가진다(9).

2pr+k2pr=0(1) 

여기서 r은 음향 공동 내부의 한 점에 대한 위치 벡터를 의미하고, p(r)은 음향 공동 내부의 한 점에서의 음압(sound pressure)을 나타내며, k는 파수(wavenumber)이다. 그리고 강체벽(rigid wall) 음향 공동의 경계 r에서의 경계 조건은 다음과 같이 주어진다.

prΓn=0(2) 

여기서 rΓ는 음향 공동 경계상의 한 점에 대한 위치 벡터이며, n은 경계에서의 법선 방향을 의미한다.

2.2 분할된 영역에 대한 시스템 행렬식

Fig. 1에서와 같이 해석 대상 음향 공동(실선으로 표시된 다각형)을 두 개의 볼록 다각형 영역 DI와 DII로 분할한다. 영역 DINa개의 경계 모서리 Γ1(I), Γ2(I), …, ΓNa(I)로 둘러싸여 있다고 가정하며, ΓNa(I)은 분할된 두 영역 사이의 접경(interface)인 Γc에 해당한다. 그리고 이들 Na개의 경계 모서리 각각에 대해 지역(local) 좌표계 (x1(I), y1(I)), (x2(I), y2(I)), …, (xNa(I), yNa(I))를 설정한다. 마찬가지로, 영역 DIINb개의 경계 모서리 Γ1(II), Γ2(II), …, ΓNb(II)로 둘러싸여 있다고 가정하며, ΓNb(II)는 분할된 두 영역 사이의 접경인 Γc에 해당한다. 그리고 이들 Nb개의 경계 모서리 각각에 대해 국부 직각 좌표계 (x1(II), y1(II)), (x2(II), y2(II)), …, (xNb(II), yNb(II))를 설정한다.

Fig. 1

Arbitrarily shaped, concave acoustic cavity divided into two regions DI and DII

이 논문에서는 분할된 영역 각각에 대한 시스템 행렬식을 먼저 추출할 것이다. 이를 위해 영역 DI에서의 음압 p(I)식 (3)과 같이, DI의 경계 모서리들에서 발생한 평면파(plane wave)의 선형 중첩으로 가정한다.

pI=i=1Nam=0NsAmicosmπxiILi(I)×expjk2-mπLiI2yi(I)(3) 

여기서 Am(i)은 평면파의 진폭 미지 계수, Li(I)i번째 경계 모서리의 길이, (x1(I), y1(I))는 i번째 경계 모서리에 설정된 지역 좌표계, Ns는 평면파에 사용된 코사인 시리즈 함수의 개수, j는 순허수를 나타낸다.

다음으로 영역 DI를 둘러싸고 있는 경계 모서리들에서의 음압의 기울기는 코사인 시리즈 함수의 선형 중첩으로 가정한다.

pIyrI(xrIyrI=0)=n=1NsUnrcosnπxrILrIr=1,2,···,Na(4) 

여기서 Un(r)은 미지 계수, (xr(I), yr(I))은 r번째 경계 모서리에 설정된 지역 좌표계를 나타낸다.

이제 경계조건 식 (4)를 가정된 음압 식 (3)에 적용한 후, 경계 모서리를 따라 적분을 수행하면 식 (5)를 얻을 수 있다.

i=1Na(m=0NsAm(i)SMIq,mi,r)=Uqrq=0,1,,Ns; r=1,2,,Na(5) 

여기서

SMIq,mi,r=dLrI0LrIcosmπxrILrI×expjk2-mπLiI2yiIsinqπxrILrIyrI=0dxrIford=1,q=0d=2,q=1,2,,Ns(6) 

나중에 두 영역의 접경 Γc에서 음압과 음압 기울기 연속 조건을 적용하기 위하여, 식 (5)는 다음과 같이 경계 모서리 Γ1(I), Γ2(I), …, ΓNa-1(I)와 접경 Γc(=ΓNa(I))에 대한 식으로 분리된다.

i=1Na-1m=0NsAmiSMIq,mi,r+m=0NsAmNaSMIq,mi,r=Uqrq=0,1,,Ns;r=1,2,,Na-1(7) 
i=1Na-1(m=0NsAmiSMIq,mi,Na)+m=0NsAmNaSMIq,mNa,Na=UqNaq=0,1,,Ns(8) 

여기서 접경을 제외한 나머지 경계 모서리는 강체벽 경계라는 사실을 고려하면, 식 (7)의 우변은 0이 된다(즉, Uq(r)=0, r=1,2,…,Na-1). 마지막으로 식 (7)식 (8)은 다음과 같이 시스템 행렬식의 형태로 각각 표현될 수 있다.

SMI(bb)AIb+SMIbcAIc=0(9) 
SMI(cb) AI(b)+SMI(cc)AI(c)=UI(c)(10) 

여기서 위첨자 b와 c는 경계 모서리(boundary edge)와 접경(interface)을 각각 의미한다.

이제 영역 DII에 대한 시스템 행렬식을 구하기 위하여, 식 (3)식 (4)와 같은 방식으로 영역 DII의 내부 음압과 경계 모서리에서의 음압 기울기를 다음과 같이 각각 가정한다.

pII=i=1Nbm=0NsBmicosmπxiIILiII×expjk2-mπLiII2yiII(11) 
p(II)yr(I)(xr(II),yr(II)=0)=n=1NsVn(r)cosnπxr(II)Lr(II)r=1,2,,Nb(12) 

여기서 Bm(i)Vn(r)은 미지 계수를 나타내며, (xi(II), yi(II))는 영역 DIIi번째 경계 모서리에 설정된 지역 좌표계를 나타낸다.

다음으로, 영역 DI에 대해 수행된 식 (3~10)의 과정을 식 (11, 12)에 적용을 하면, 다음과 같은 영역 DII에 대한 시스템 행렬식을 얻을 수 있다.

SMII(bb)BII(b)+SMII(bc)BII(c)=0(13) 
SMII(cb)BII(b)+SMII(cc)BII(c)=VII(c)(14) 

2.3 경계조건과 연속조건 적용

접경 Γc에서의 음압과 음압 기울기 연속 조건을 고려하는 방법에 의해, 시스템 행렬식 식 (9, 10, 13, 14)를 하나의 시스템 행렬식으로 통합하고자 한다. 접경 Γc에서 음압 기울기 연속 조건을 고려하면, 식 (10)의 우변과 식 (14)의 우변이 같아야 한다. 즉,

UIc=-VIIc(15) 

여기서 우변에 -1이 곱해진 이유는 접경 ΓNa(I)ΓNa(II)에 각각 설정된 지역 좌표계의 yNa(I)yNb(II) 방향이 반대이기 때문이다. 식 (15)식 (10, 14)에 적용하면 다음의 식을 얻을 수 있다.

SMI(cb)AI(b)+SMI(cc)AI(c)+SMII(cb)BII(b)+SMII(cc)BII(c)=0(16) 

추가적으로 접경 Γc에서 음압 연속 조건을 고려하기 위하여, 영역 DI의 관점에서 접경 ΓNa(I)(=Γc)에서의 음압을 다음과 같이 가정한다.

p(I)xNa(I),yNa(I)=0=n=0NsU^nNacosnπxNa(I)LNa(I)(17) 

여기서 U^nNa는 미지 계수이다. 다음으로 영역 DI의 음압 식 (3)식 (17)에 대입한 후, 접경을 따라 적분을 수행하면, 다음의 식을 구할 수 있다.

i=1Nam=0NsAm(i)PMIq,mi,Na)=U^qNa,q=0,1,,Ns(18) 

여기서

PMIq,mi,Na=dLrI0LrIxrIcosmπxiILiIexpjk2-mπLiI2yiIsinqπxrILrIyrI=0dxrI(19) 

마지막으로 식 (19)는 다음과 같이 행렬식의 형태로 간단히 표현된다.

PMI(cb)AI(b)+PMI(cc)AI(c)=U̯I(c)(20) 

마찬가지로, 영역 DII의 관점에서 접경 ΓNb(II)(=Γc)에서의 음압은 다음과 같이 가정된다.

p(II)(xNb(II),yNb(II)=0)=n=0NsV^nNbcosnnpixNb(II)LNb(II)(21) 

여기서 V^nNa는 미지 계수이다. 다음으로 영역 DII의 음압 식 (11)식 (21)에 대입한 후, 접경을 따라 적분을 수행하면, 다음의 식을 구할 수 있다.

i=1Nb(m=0NsBm(i)PMIIq,m(i,Nb))=V̯q(Nb),q=0,1,···,Ns(22) 

여기서

PMIIq,mi,Nb=dLrII0LrIIxrIIcosmπxiIILiII×exp(jk2-mπLiII2yiIIsinqπxrIILrIIyrII=0dxrII (23) 

마지막으로 식 (22)는 다음과 같이 행렬식의 형태로 간단히 표현된다.

PMII(cb)AII(b)+PMII(cc)AII(c)=V^II(c)(24) 

접경에서 음압은 연속이므로, 식 (20)식 (24)에서 U^Ic=V^IIc이 성립한다. 그래서 다음의 식을 식 (20)식 (24)으로부터 얻을 수 있다.

PMI(cb)AI(b)+PMI(cc)AI(c)-PMII(cb)BII(b)-PMII(cc)BI(c)=0(25) 

2.4 시스템 행렬의 크기 축소와 고유치 추출

앞에서 구한 행렬식 식 (9, 13, 16, 25)를 하나의 행렬식으로 머지(merge)하면 다음과 같이 하나의 시스템 행렬식을 얻을 수 있다.

SMC=0(26) 

여기서 시스템 행렬 SM와 벡터 C는 다음과 같이 주어진다.

SM=SMI(bb)SMI(bc)0000SMII(bb)SMI(bc)SMI(cb)SMI(cc)SMII(cb)SMI(cc)PMI(cb)PMI(cc)-PMII(cb)-PMII(cc)(27) 
C=AI(b)AI(c)BII(b)BII(c)(28) 

음향 공동의 고유치는 SM의 판별식(determinant)이 영이 되는 조건인 식 (29)로부터 구해질 수 있다.

detSM=0(29) 

한편, 식 (29)에서 시스템 행렬 SM은 크기가 (Na+Nb)Ns인 정사각 행렬이다. 만약 해석 대상 음향 공동의 모서리의 수 또는 사용된 시리즈 함수의 개수가 많아지게 되면, 시스템 행렬의 크기가 커지게 되어 식 (29)를 계산하는 과정에서 수치 계산량이 증가하는 문제점이 발생한다. 이 논문에서는 이러한 문제점을 해결하기 위해, 시스템 행렬의 크기를 축소하는 방안을 제시한다.

먼저 식 (9, 13)은 다음과 같이 변형한다.

AI(b)=-SMI(bb)-1SMI(bc)AI(c)(30) 
BII(b)=-SMII(bb)-1SMII(bc)BII(c)(31) 

다음으로, 식 (30, 31)식 (16, 25)에 각각 대입하면 다음의 식들을 얻을 수 있다.

(SMI(cc)-SMI(cb)SMI(bb)-1SMI(bc))AI(c)+(SMII(cb)SMII(bb)-1SMII(bc)-SMII(cc))BII(c)=0(32) 
(PMI(cc)-PMI(cb)SMI(bb)-1SMI(bc))AI(c)+(PMII(cb)SMII(bb)-1SMII(bc)-PMII(cc))BII(c)=0(33) 

마지막으로, 식 (32, 33)은 다음과 같이 하나의 시스템 행렬식으로 표현될 수 있다.

SMrCr=0(34) 

여기서

SMr=SMI(cc)-SMI(cb)×SMI(bb)-1SMI(bc)SMII(cb)SMII(bb)-1×SMII(bc)-SMII(cc)PMI(cc)-PMI(cc)×SMI(bb)-1SMI(bc)PMII(cb)SMI(bb)-1×SMII(bc)-PMII(cc)(35) 
Cr={AI(c)BII(c)(36) 

축소된(reduced) 정사각 시스템 행렬 SMr의 크기는 2Ns이며, SM에 비해 크기가 2/(Na+Nb)배 만큼 감소했음을 알 수 있다. 마지막으로 고유치는 식 (37)에 의해 구해진다. 참고로 수치 해석 프로그램인 Mathematia를 사용하여 식 (37)로부터 고유치를 추출하였다.

detSMr=0(37) 

3. 검증 예제

3.1 직사각형 음향 공동

이 논문에서 제안된 방법의 타당성 및 정확성을 검증하기 위해 첫 번째 예제로 엄밀해가 존재하는 직사각형 2차원 음향 공동이 고려되었다. Fig. 2와 같이 가로 길이가 1.2 m, 세로 길이가 0.9 m인 직사각형 음향 공동을 의도적으로 두 개의 영역 DI과 DII로 나눈 다음, 이 논문에서 제안된 방법이 적용되었다.

Fig. 2

Rectangular acoustic cavity divided into two regions

Table 1은 이 논문에서 제안된 방법(proposed method), 엄밀해(exact solution)(10), FEM(ANSYS)에 의해 구해진 고유치 결과들을 보여준다. 표에서 알 수 있듯이, 이 논문에서 제안된 방법은 단지 Ns=3일 때 엄밀해에 완전히 수렴된 정확한 결과를 제공함을 확인할 수 있으며, 2500개의 많은 노드를 사용한 FEM보다 더 정확한 결과를 제공함을 확인할 수 있다.

Eigenvalues of the rectangular cavity by the proposed method, the exact method, and FEM (parenthesized values denote errors (%) with respect to the values by the exact method)

3.2 임의 형상 오목 음향 공동

이 예제에서는 Fig. 3과 같이 앞의 예제에서 고려된 직사각형 음향 공동의 위쪽 모서리가 0.2 m 들어간 오목 음향 공동이 고려되었다. 그리고 이 논문에서 제안된 방법을 적용하기 위하여 그림과 같이 두 개의 볼록 영역 DI과 DII로 분할하였다.

Fig. 3

Arbitrarily shaped, concave acoustic cavity divided into two regions

Eigenvalues of the arbitrarily shaped, concave cavity by the proposed method and FEM (parenthesized values denote errors (%) with respect to the values by FEM with 2806 nodes)

Table 2는 이 논문에서 제안된 방법(proposed method)과 FEM(ANSYS)에 의해 구해진 고유치 결과들을 보여준다. Table 2에서 확인할 수 있듯이, 제안된 방법은 단지 Ns=5인 경우에 고유치가 수렴되었음을 알 수 있다. 그리고 이 결과는 2806개의 많은 노드를 사용한 FEM 결과와 비교해볼 때 0.11 % 오차 이내의 정확한 값임을 확인할 수 있다. 참고로, 이 예제는 엄밀해가 존재하지 않는 경우이기에, 많은 노드를 사용한 FEM 해석 결과를 오차를 계산할 때의 기준값으로 사용하였다.

한편, Table 2에서 제안된 방법에 의해 추출된 네 번째 고유치 4.527은 FEM에서는 추출되지 않는 허위 고유치(spurious eigenvalue)이다. 이 허위 고유치는 Fig. 2의 영역 DI(또는 DII)에 해당하는 강체벽 음향 공동의 세 번째 고유치인 것으로 확인되었다. 이러한 허위 고유치가 발생한 원인과 제거 방안에 대한 연구가 현재 진행 중에 있으며, 향후 논문에서 관련 연구 결과를 제시할 예정이다.


4. 결론

이 논문에서는 임의 형상 오목 음향 공동의 고유치 해석을 위해, 평면파를 이용한 영역 분할법 기반 고유치 해석 기법이 제안되었다. 두 가지 예제를 통하여 이 논문에서 제안된 방법의 타당성과 정확성이 검증되었다. 차후의 연구에서는 음향 공동의 고유 모드 추출 방법, 허위 고유치 제거 방안이 추가로 연구될 예정이다.

비록 이 논문에서는 두 가지 예제를 이용하여, 제안된 이론의 타당성을 검증하였지만, 향후 연구에서는 보다 복잡한 형상을 가진 음향 공동에 대한 예제가 다루어질 것이다. 이를 위해 두 개 영역보다 더 많은 영역으로 영역을 분할하는 방법에 대한 이론 정식화가 이루어질 것이다. 한편 강체벽 경계와 개방 경계로 구성된 혼합 경계 음향 공동에 대한 연구는 곧 발표될 예정이며, 3차원 음향 공동으로의 확장 연구도 진행중에 있다.

Acknowledgments

이 연구는 한성대학교 교내학술연구비 지원 과제임.

References

  • Kang, S. W., and Lee, J. M., (2000), Eigenmode Analysis of Arbitrarily Shaped Two-dimensional Cavities by the Method of Point-matching, Journal of the Acoustical Society of America, 107(3), p1153-1160. [https://doi.org/10.1121/1.428456]
  • Kang, S. W., and Yon, J. I., (2012), Development of a Modified NDIF Method for Extracting Highly Accurate Eigenvalues of Arbitrarily Shaped Acoustic Cavities, Transactions of the Korean Society for Noise and Vibration Engineering, 22(8), p742-747. [https://doi.org/10.5050/KSNVE.2012.22.8.742]
  • Bathe, K., (1982), Finite Element Procedures in Engineering Analysis, Prentice-Hall, New Jersey.
  • Brebbia, C. A., Telles, J. C. F., and Wrobel, L. C., (1984), Boundary Element Techniques, Springer-Verlag, New York. [https://doi.org/10.1007/978-3-642-48860-3]
  • Kang, S. W., and Yon, J. I., (2011), Development of an Improved NDIF Method for Efficiently Extracting Eigenvalues and Eigenmodes of Arbitrarily Shaped Acoustic Cavities, Transactions of the Korean Society for Noise and Vibration Engineering, 21(10), p960-966. [https://doi.org/10.5050/KSNVE.2011.21.10.960]
  • Kang, S. W., (2013), New Formulation of MNDIF Method for Accurate Eigenvalue Analysis of Concave Acoustic Cavities, Transactions of the Korean Society for Noise and Vibration Engineering, 23(11), p1003-1011. [https://doi.org/10.5050/KSNVE.2013.23.11.1003]
  • Kang, S. W., (2007), Free Vibration Analysis of Clamped Plates with Arbitrary Shapes Using Series Functions, Transactions of the Korean Society for Noise and Vibration Engineering, 13(10), p531-538.
  • Kang, S. W., and Yon, J. I., (2009), Eigenvalue Analysis of Arbitrarily Shaped, Concave Membranes with a Deep Groove Using a Sub-domain Methods, Transactions of the Korean Society for Noise and Vibration Engineering, 19(10), p1069-1074. [https://doi.org/10.5050/KSNVN.2009.19.10.1069]
  • Meirovitch, L., (1967), Analytical methods in Vibrations, Macmillan Publishing, New York.
  • Blevins, R. D., (1979), Formulas for Natural Frequency and Mode Shape, Litton Educational Publishing, New York.
저자소개

Sang-Wook Kang received a B.S., M.S. and Ph.D. degree in Dept. of Mechanical Design Engineering from Seoul National University in 1992, 1994 and 1999, respectively. He is currently working in Hansung University as a professor of the mechanical system engineering. His research interests include the development of NDIF method for free vibration analysis of membrane, acoustic cavities, and plates with general shapes.

Fig. 1

Fig. 1
Arbitrarily shaped, concave acoustic cavity divided into two regions DI and DII

Fig. 2

Fig. 2
Rectangular acoustic cavity divided into two regions

Fig. 3

Fig. 3
Arbitrarily shaped, concave acoustic cavity divided into two regions

Table 1

Eigenvalues of the rectangular cavity by the proposed method, the exact method, and FEM (parenthesized values denote errors (%) with respect to the values by the exact method)

Proposed method Exact method(10) FEM
(2500 nodes)
Ns = 2 Ns = 3 Ns = 4
1 2.618
(0.00)
2.618
(0.00)
2.618
(0.00)
2.618
2.618
(0.00)
2 3.491
(0.00)
3.491
(0.00)
3.491
(0.00)
3.491
3.492
(0.03)
3 4.362
(0.02)
4.363
(0.00)
4.363
(0.00)
4.363
4.364
(0.02)
4 5.239
(0.06)
5.236
(0.00)
5.236
(0.00)
5.236
5.238
(0.04)
5 6.293
(0.00)
6.293
(0.00)
6.293
(0.00)
6.293
6.295
(0.03)
6 6.981
(0.00)
6.981
(0.00)
6.981
(0.00)
6.981
6.989
(0.11)

Table 2

Eigenvalues of the arbitrarily shaped, concave cavity by the proposed method and FEM (parenthesized values denote errors (%) with respect to the values by FEM with 2806 nodes)

Proposed method FEM(ANSYS)
Ns = 4 Ns = 5 Ns = 6 2806
nodes
1299
nodes
573
nodes
1 2.454
(0.08)
2.455
(0.04)
2.455
(0.04)
2.456
2.457
2.457
2 3.856
(0.00)
3.855
(0.03)
3.855
(0.03)
3.856
3.857
3.860
3 4.416
(0.05)
4.417
(0.02)
4.417
(0.02)
4.418
4.419
4.421
4 4.527 4.527 4.527 None None None
5 5.262
(0.04)
5.262
(0.04)
5.262
(0.04)
5.264
5.267
5.273
6 6.452
(0.03)
6.452
(0.03)
6.452
(0.03)
6.454
6.458
6.466
7 7.547
(0.12)
7.548
(0.11)
7.548
(0.11)
7.556
7.564
7.582