고정단 평판의 전 주파수 대역 고유치 추출이 가능한 NDIF법 개발: 제2부(예제 연구 및 고찰)

2.1 기존 NDIF법

Fig. 1에서 점선으로 표시한 모양이 해석 대상 고정단 평판의 경계에 해당된다. 고정단 평판의 경계는 Fig. 1과 같이 N개의 노드로 이산화된다. 평판의 진동 지배 방정식에 고정단 경계 조건을 적용하면, 고정단 평판의 시스템 행렬 $$ {\mathbf{S}\mathbf{M}}_{\text{plate }}^N\left(\mathrm{\Lambda }\right)$$을 추출할 수 있다. 여기서 위 첨자 N은 정사각 시스템 행렬의 크기, Λ는 주파수 파라미터를 의미한다⁽¹⁾.

고정단 평판의 고유치는 시스템 행렬의 판별식이 0이 되게 하는 주파수 파라미터에 해당되므로, 식 (1)로부터 구할 수 있다.

Fig. 1

Arbitrarily shaped clamped plate whose boundary Γ is discretized with boundary nodes P1, P2,⋯, PN

식 (1)을 이용할 경우, 노드 수가 증가하면 저주파수 고유치가 구해지지 않는 문제점이 존재한다⁽⁴⁾.

2.2 논문 1부에서 제안된 개선된 NDIF법

논문 1부에서 새로이 제안된 방법을 이용하면, 크기가 2N인 시스템 행렬 $$ {\mathbf{SM}}_{\text{plate }}^{2N}\left(\mathrm{\Lambda }\right)$$을 추출할 수 있다. 고정단 평판의 고유치는 시스템 행렬의 판별식의 해, 즉 식 (2)의 해에 해당된다⁽⁴⁾.

식 (2)의 판별식 값은 일반적으로 식 (3)에 의해 계산된다.

여기서 λ_i(Λ)는 $$ {\mathbf{SM}}_{\text{plate }}^{2N}\left(\mathrm{\Lambda }\right)$$에 대한 대수 고유치 문제(algebraic eigenvalue problem) 식인 식 (4)에서 i 번째 고유치를 의미한다.

여기서 v_i는 i번째 고유벡터를 의미한다.

한편, 식 (3)을 이용하더라도 노드 수가 증가하면 저주파수 영역에서 고유치가 구해지지 않음이 확인되었다⁽⁴⁾. 그래서 식 (3) 대신에 식 (5)가 논문 1부에서 제안되었다.

여기서 R(Λ)은 주파수 파라미터 Λ의 함수로 주어지는 시스템 행렬 $$ {\mathbf{SM}}_{\text{plate }}^{2N}\left(\mathrm{\Lambda }\right)$$의 랭크 함수이다.

식 (5)를 사용할 경우 저주파수 영역에서 고유치가 정확히 구해짐이 확인되었다⁽⁴⁾. 그러나 시스템 행렬의 랭크가 불연속적으로 변하는 지점들이 많아 질 경우, 고유치 식별이 어려운 문제점이 발생하였다. 이 문제점을 개선하기 위해, 랭크의 변화에 의해 추가로 곱해지는 고유치들의 곱들을 식 (5)의 판별식에서 나누어 주는 방안인 식 (6)이 제안되었다⁽⁴⁾.

3.1 직사각형 고정단 평판

Fig. 2는 엄밀해가 존재하지 않는 고정단 직사각형 평판(가로 1.2 m, 세로 0.9 m)을 보여준다. 먼저 Fig. 2(a)와 같이 평판의 경계를 20개의 노드로 이산화하고, 코너 노드에서의 법선 벡터는 양쪽 변의 두 법선 벡터의 평균으로 구해진다. 20개 노드로 이산화 된 평판에 대해 기존 NDIF법의 식 (1)을 이용하여 판별식 곡선을 그리면 Fig. 3의 실선과 같다. 이 실선 곡선에서 6개 고유치들(p₁ ~ p₆)이 추출되었음을 확인할 수 있다. 이들 고유치들은 Table 1의 두번째 열에 요약되었으며, FEM(ANSYS)에 의해 구한 고유치들과 0.11 % 이내의 적은 오차를 가짐을 확인할 수 있다. 참고로 m₁ ~ m₉는 허위 고유치(spurious eigenvalue)를 의미하며, 본 평판과 같은 형상을 가진 멤브레인 고유치와 일치한다⁽¹⁾.

Fig. 2

Rectangular clamped plates discretized with 20 and 28 nodes, respectively

Fig. 3

Determinant curves for the rectangular clamped plates by the original NDIF method using Eq. (1) (dotted line: 20 nodes, solid line: 28 nodes)

Table 1

Eigenvalues of the rectangular clamped plate by the original NDIF method, the proposed method, and FEM(ANSYS) (parenthesized values denote errors[%] with respect to FEM results)

Fig. 2(b)와 같이 노드의 개수를 28개로 늘인 후 기존 NDIF법의 식 (1)을 이용하여 판별식 곡선을 그리면 Fig. 3의 점선과 같다. 이 점선 곡선을 살펴보면, 저주파수 영역에서 판별식 값이 발산하여 단지 4개의 고유치(p₃ ~ p₆)만이 추출되는 문제점을 확인할 수 있으며, 이 결과는 Table 1의 3열에 요약되었다. 한편, 28개의 노드를 사용한 경우에 대한 시스템 행렬의 랭크 함수 R(Λ)를 구해 보면 Fig. 4와 같다. 시스템 행렬의 풀 랭크(full rank) 값은 56인데, 랭크 함수 곡선의 저주파수 영역을 살펴 보면 랭크 값이 56보다 훨씬 작은 값을 가짐을 알 수 있다. 이러한 이유로 인해 Fig. 3의 점선에서 저주파수 고유치가 추출되지 않는다는 것을 논문 1부에서 이미 설명하였다.

Fig. 4

Rank function R(Λ) of the rectangular clamped plate discretized with 28 nodes

상기의 문제점을 개선하기 위해 논문 1부에서 제안한 식 (5)를 이용하여 판별식 곡선을 그리면 Fig. 5의 점선과 같으며, 식 (5)의 계산에 필요한 시스템 행렬의 랭크 R(Λ)은 앞의 Fig. 4에 의해 주어진다. Fig. 5의 점선 곡선을 살펴보면, 시스템 행렬의 랭크 변화에 의해 발생한 불연속 지점들과 평판의 고유치에 해당하는 골(crest)들이 서로 겹쳐 있어서, 고유치 식별이 어려운 상황임을 알 수 있다. 그래서, 랭크 변화에 의한 불연속 지점들을 제거하기 위해 제안된 식 (6)을 이용하여 판별식 곡선을 그리면 Fig. 5의 실선과 같이 된다. 이 실선 곡선을 살펴보면, 평판의 6개 고유치들(p₁ ~ p₆)이 정확히 추출되었음을 확인할 수 있다. 이들 고유치들은 Table 1의 네번째 열에 요약되었으며, 다섯번째 열의 FEM(ANSYS) 해석 결과와 비교해볼 때 거의 정확히 일치함도 확인 가능하다. FEM 해석은 961개의 많은 노드를 사용한 반면에, 제안된 방법은 단지 28개의 노드를 사용했다는 점을 고려하면, 제안된 방법은 FEM 보다 수치 계산량에 있어서 더 효율적인 방법이라 말할 수 있다.

Fig. 5

Determinant curves for the rectangular clamped plate with 28 nodes obtained by the proposed method (dotted and solid lines use Eq. (5) and Eq. (6), respectively)

Fig. 6은 논문 1부의 식 (7)과 식 (33)을 사용하여 제안된 방법에 의해 추출된 고정단 사각 평판의 고유 모드를 보여준다. 이들 모드들은 Fig. 2(b)와 같이 단지 28개의 노드만이 사용되었지만 FEM(ANSYS)으로 구한 모드들과 노달 라인과 형상이 정확히 일치함이 확인되었다. 참고로 지면 관계 상 FEM(ANSYS)으로 구한 모드들은 제시하지 않았다.

Fig. 6

Mode shapes of the rectangular clamped plate obtained by the proposed method using Fig. 2(b)

3.2 임의 형상 고정단 평판

Fig. 7은 반원과 2개의 선분으로 구성된 임의 형상 고정단 평판을 보여준다. 먼저 Fig. 7(a)와 같이 평판의 경계를 20개의 노드로 이산화 후, 기존 NDIF법의 식 (1)을 이용하여 판별식 곡선을 그리면 Fig. 8의 실선과 같다. 이 실선 곡선에서 6개 고유치들(p₁ ~ p₆)이 추출되었음을 확인할 수 있다. 참고로 피크 m₁ ~ m₆은 멤브레인 고유치(허위 고유치)에 해당된다. 위에서 추출한 평판 고유치 p₁ ~ p₆은 Table 2의 두번째 열에 요약되었으며, FEM(ANSYS)에 의해 구한 고유치들과 0.17 % 이내의 적은 오차로 가짐을 확인할 수 있다. 참고로 m₁ ~ m₉는 허위 고유치를 의미하며 멤브레인 고유치와 일치한다⁽¹⁾.

Fig. 7

Arbitrarily shaped clamped plates discretized with 20 and 28 nodes, respectively

Fig. 8

Determinant curves for the arbitrarily shaped clamped plates by the original NDIF method using Eq. (1) (dotted line: 20 nodes, solid line: 28 nodes)

Table 2

Eigenvalues of the arbitrarily shaped clamped plate by the original NDIF method, the proposed method, and FEM(ANSYS) (parenthesized values denote errors[%] with respect to FEM results)

Fig. 7(b)와 같이 노드의 개수를 28개로 늘인 후 기존 NDIF법의 식 (1)을 이용하여 판별식 곡선을 그리면 Fig. 8의 점선 곡선과 같다. 이 판별식 곡선의 저주파수 영역에서 판별식 값이 발산하여 단지 여섯번째 고유치(p₆)만이 추출되는 문제점을 확인할 수 있으며, 이 결과는 Table 2의 세번째 열에 요약되었다.

상기의 문제점을 개선하기 위해 논문 1부에서 제안한 식 (5)를 이용하여 판별식 곡선을 그리면 Fig. 9의 점선과 같다. 그리고 식 (5)를 계산하기 위해 사용된 시스템 행렬의 랭크 R(Λ)은 Fig. 10과 같이 구해졌다. Fig. 9 점선 곡선을 살펴보면, 이전 예제와 마찬가지로 행렬의 랭크 변화에 의해 발생한 불연속 지점들과 평판의 고유치에 해당하는 골(crest)들이 서로 겹쳐 있어서, 고유치 식별이 어려운 상황임을 알 수 있다. 그래서, 랭크 변화에 의한 불연속 지점들을 제거하기 위해 제안된 식 (6)을 이용하여 판별식 곡선을 그리면 Fig. 9 실선과 같이 된다. 이 실선 곡선을 살펴보면, 평판의 6개 고유치들(p₁ ~ p₆)이 정확히 추출되었음을 확인할 수 있다. 이들 고유치들은 Table 2 네번째 열에 요약되었으며, 네번째 열의 FEM(ANSYS) 해석 결과와 비교해볼 때 거의 정확히 일치함도 확인 가능하다. FEM 해석은 961개의 많은 노드를 사용한 반면에 제안된 방법은 단지 28개의 노드를 사용했다는 측면에서, 제안된 방법은 FEM 보다 수치 계산량에 있어서 더 효율적인 방법이라 말할 수 있다.

Fig. 9

Determinant curves for the arbitrarily shaped clamped plate with 28 nodes obtained by the proposed method (dotted and solid lines use Eq. (5) and Eq. (6), respectively)

Fig. 10

Rank function R(Λ) of the arbitrarily shaped clamped plate discretized with 28 nodes

Fig. 11은 제안된 방법(논문 1부의 식 (7)과 식 (33)을 사용)에 의해 추출된 고정단 임의 형상 평판의 고유 모드를 보여준다. 이들 모드들은 FEM(ANSYS)으로 구한 모드들과 노달 라인과 형상이 정확히 일치함이 확인되었다.

Fig. 11

Mode shapes of the arbitrarily shaped clamped plate obtained by the proposed method using Fig. 7(b)