NDIF법(non-dimensional dynamic influence function method 또는 무차원 동영향 함수법)은 저자가 임의 형상 고정단 평판의 고정밀도 자유진동해석을 위해 개발한 고유치 해석 기법이다⁽¹⁾. 저자는 혼합 경계 조건을 가진 임의 형상 평판과 자유단 평판의 고정밀도 고유치 해석을 위한 NDIF법도 성공적으로 개발하였다^(2~4). 이후, 저자는 평판의 고유치 해석에 있어서 NDIF법의 해석 효율성과 정확도를 높이기 위한 다양한 연구를 수행하였다^(5~8).

NDIF법은 해석 대상 평판의 내부를 요소(element)로 분할하는 유한요소법(FEM)과는⁽⁹⁾ 달리 해석 대상 평판의 경계만을 이산화하기 때문에, 이산화 방법의 관점에서는 경계요소법(BEM)과⁽¹⁰⁾ 유사하다고 말할 수 있다. 하지만, 경계요소법과는 달리 NDIF법은 평판의 경계를 요소 없이 노드로만 이산화하기 때문에 노드와 노드 사이에 설정된 보간 함수의 특이치 적분 과정이 필요 없는 특징을 가진다. 이러한 특징 덕분에, NDIF법의 이론이 간결하고 그 해석 결과도 유한요소법 및 경계요소법에 비해 더 정확한 장점을 가진다.

그러나 NDIF법은 평판의 고차 고유치와 고유 모드를 추출하기 위해 노드의 개수를 증가시킬 경우 저주파수 영역에서 시스템 행렬이 발산하여 저차 고유치가 구해지지 않는 단점이 존재한다. 그래서 NDIF법에서는 고차 고유치를 추출할 때, 저차 고유치 추출을 위한 별도의 이산화 과정이 추가로 필요한 비효율성이 존재한다.

이 논문에서는 상기의 비효율성을 제거하기 위한 연구가 수행된다. 임의 형상 멤브레인과 음향 공동에 대해 최근에 수행된 연구에서는^(11,12), 노드의 개수를 증가시켰을 때 저주파수 고유치가 구해지지 않는 문제점이 개선된 NDIF법이 성공적으로 개발되었다. 이번 논문에서는 위의 연구 결과를 평판으로 확장하는 시도가 이루어졌다.

이 논문에서는 평판의 시스템 행렬이 노드의 개수와 상관없이 전 주파수 대역에서 정확한 고유치와 고유모드를 제공하는 NDIF법 이론이 성공적으로 개발되었으며, 엄밀해가 존재하는 고정단 원형 평판에 개발된 이론을 적용하여 그 타당성과 정확성이 검증되었다. 참고로, 이 논문의 2부에서는 다양한 형상의 고정단 평판에 대해 개발된 이론을 적용하는 예제 연구가 수행되었다.

이전의 멤브레인에 관한 연구⁽¹¹⁾에서 시스템 행렬의 랭크가 시스템 행렬의 크기 보다 작을 경우에는 시스템 행렬의 판별식 값이 발산함이 확인되었다. 그래서 이러한 발산 문제를 없애기 위해 이전 연구에서 제안한, 시스템 행렬의 랭크를 고려한 시스템 행렬 판별식 계산 방법⁽¹¹⁾을 고려하면, 고정단 원형 평판의 고유치 추출 식인 식 (26)의 좌변, 즉 시스템 행렬((SM_plate(Λ)))의 판별식은 식 (27)의 방법에 의해 계산 가능하다.

여기서 R(Λ)은 주파수 파라미터 Λ의 함수로 주어지는 시스템 행렬 SM_plate(Λ)의 랭크이며, λ_i(Λ)는 식 (28)과 같은 SM_plate(Λ)에 대한 대수 고유치 문제에서 i번째 고유치를 의미한다.

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

더 나아가, 식 (27)을 이용하여 시스템 행렬의 판별식 곡선을 그릴 경우, 시스템 행렬의 랭크가 불연속적으로 변하는 주파수 파라미터 값들 위치에서 판별식 곡선도 불연속적으로 변하는 현상이 발생한다^(11,12). 이러한 판별식 곡선의 불연속 현상을 제거하기 위해, 랭크의 변화에 의해 추가로 곱해지는 고유치들의 곱들을 식 (27)의 판별식 계산식에서 나누어 주는 방안^(11,12)인 식 (29)가 이용된다.

식 (29)를 이용하여 30개의 노드로 이산화한 고정단 원형 평판의 시스템 행렬(SM_plate)에 대한 판별식 곡선을 그린 결과는 Fig. 3의 실선과 일치하는 것으로 확인되었다. 그 이유는 식 (24)에 의해 주어지는 SM_plate의 랭크가 전체 주파수 구간에서 30으로 일정하기 때문인 것으로 분석되었다.

위에서 살펴본 바와 같이, 식 (29)를 이용하는 최근의 연구 방법⁽¹¹⁾으로는 평판 시스템 행렬의 저주파수 발산 문제를 해결할 수 없음을 알 수 있다. Fig. 3의 실선으로 주어진 판별식 곡선의 저주파수 영역에서의 발산 문제를 해결하기 위해선, 저주파수 영역에서 시스템 행렬의 랭크의 크기가 시스템 행렬의 크기보다 작아야 한다. 그런데 앞에서 설명한 바와 같이, 현재의 경우는 시스템 행렬의 랭크의 크기가 시스템 행렬의 크기와 일치하기 때문에, 식 (29)가 그 효용성을 발휘하지 못하는 것으로 판단된다.

SM_plate의 랭크가 30으로 일정한 원인은, 식 (24)에서 SM_plate을 계산할 때 로컬 시스템 행렬의 역행렬이 사용되었기 때문으로 판단된다. 그래서 논문에서는 시스템 행렬을 다른 방법으로 구할 방안이 제안된다. 식 (12)와 식 (13)에 의해 주어진 로컬 시스템 행렬식은 다음과 같이 역행렬을 사용한 행렬 축약 과정 없이 식 (30)과 같은 하나의 시스템 행렬식으로 나타낸다.

식 (31, 32)는 여기서,

식 (30)에서 위첨자 2N은 시스템 행렬의 크기를 의미한다. 참고로, 식 (24)와 같이 역행렬을 이용해서 구해진 시스템 행렬은 축약이 되어서 그 크기가 N이다. 식 (30)의 시스템 행렬을 식 (29)에 대입하면 식 (33)과 같다.

식 (33)을 이용하여 판별식 곡선을 구하면 Fig. 4와 같다. 이 곡선에 의해 구해진 고유치 값들은 Fig. 4에 표시되어 있으며, 이들 값들은 Table 1의 엄밀해와 정확히 일치함이 확인된다. 이러한 사실은, 이 논문에서 제안된 방법이 노드의 개수가 증가하더라도 전 주파수 대역에서 고정단 평판의 고유치들을 정확하게 추출할 수 있음을 보여준다. 참고로 Fig. 5를 보면, 시스템 행렬의 랭크가 저주파수 영역에서 시스템 행렬의 크기보다 작은 정상적인 현상을 확인할 수 있다.

Fig. 4 
Determinant curve for the clamped circular plate with 30 boundary nodes, obtained by the proposed method using Eq. (33)

Fig. 5 
Rank of SMplate 2N for the clamped circular plate discretized using 30 boundary nodes

제안된 방법에 의해 구해진 고유모드 형상은 Fig. 6에 제시되었으며 이들 형상은 엄밀해 고유모드 형상과 정확히 일치함이 확인되었다.

Fig. 6 
Mode shapes of the clamped circular plate obtained by the proposed method using Eq. (33) and Eq. (7)

기존의 NDIF법에서 고정단 평판의 고유치를 추출할 때, 노드의 수를 증가시키면 저주파수 영역에서 시스템 행렬이 발산하여 저차 고유치가 구해지지 않는 문제점이 존재한다. 이 논문에서는 상기의 문제점을 개선하기 위하여 시스템 행렬의 판별식을 계산할 때 시스템 행렬의 랭크를 고려하는 방안이 새로이 제안되었다. 제안된 방법의 타당성과 정확성은 엄밀해가 존재하는 고정단 원형 평판에 대한 예제 검증을 통해 확인되었다. 논문의 2부에서는 다양한 형상을 가진 임의 형상 고정단 평판에 대한 예제 연구들이 수행되었다.