Exact Solutions for Bending Vibration of Beam with Linearly Reduced width Along Its Length

길이를 따라 선형적으로 감소된 폭을 가지는 보의 굽힘 진동에 대한 정확해

Abstract

In this paper a transfer matrix method is developed to solve for bending vibration of beam with linearly reduced width, and subsequently used to determine the exact natural frequencies for such problems. The differential equation, shear force, and bending moment are derived from Hamilton's principle, and the roots of the differential equation are computed using the power series solution of the Frobenius method. The effect of various taper ratio for bending vibration of beam with linearly reduced width is investigated in detail, and to validate the accuracy of the proposed method the results computed are compared with those given from commercial software(ANSYS).

1. 서 론

테이퍼진 보(tapered beam)에 대한 자유진동특성(free vibration characteristics) 해석은 많은 연구자들에 의해서 다양한 방법으로 수행되었고, 공학적 설계에 널리 사용되고 있다(1~11). 테이퍼진 보는 폭(width)과 높이(height)가 동시에 감소된 것과 높이 치수만 감소된 것, 그리고 보의 길이를 따라 폭이 감소된 보로 구분된다. 일반적으로 많이 연구 되는 것은 높이가 감소된 보 요소에 대한 굽힘 진동이다. 테이퍼진 보의 굽힘 진동을 연구하기 위해 대부분의 연구자들이 형상함수(shape function)를 가정하였고, 이를 사용하여 선형적으로 단면 치수가 감소되는 보에 대한 연구가 보편적으로 수행하였다. 또한, 단진보(stepped beam)를 많은 요소로 이상화함으로써 테이퍼진 보를 연구할 수 있다. 하지만 길이를 따라 감소되는 폭을 가진 보 요소에 대한 조사는 극히 제한된 경우에 대해서만 연구되었기 때문에 이 연구에서 다양한 테이퍼 비(taper ratio)에 대하여 자세하게 조사되었다.

이와 관련하여, 동강성법(1~3)을 이용한 유한요소법이 테이퍼진 보에 대한 굽힘 진동 특성을 연구하기 위하여 많은 연구자들에 의하여 연구되었고, Rayleigh-Ritz 법(4,5)과 미분변환법(differential transform method)(6) 그리고 전달행렬법(transfer matrix method)(7) 등이 고려되었다. 그러나 유한요소법이나 다른 근사법과 같은 방법들에 비하여 상대적으로 적은 연구들이 전달행렬법에 대하여 수행되었다. 따라서 이번 연구에서는 테이퍼진 보의 고유진동수를 결정하기 위하여 전달행렬법이 고려되었고, 미분방정식의 근을 결정하기 위하여 멱급수(power series)의 하나인 Frobenius method이 고려되었다(1,8,9).

이번 연구에서의 목적은 전달행렬법을 사용하여 길이를 따라 선형적으로 폭이 감소되는 보의 정확한 고유진동수를 결정하기 위한 것이다. 제안된 방법은 기존의 전달행렬 연구들과는 다르게 장행렬(field matrix)과 점행렬(point matrix)같은 구성요소를 가지지 않지만 기존의 방법과 동일하게 전달행렬의 주요특징들을 가지고 있다. 또한, 한 개의 보 요소로서 무한한 고유진동수를 계산 할 수 있다. 제안된 방법의 타당성을 검증하기 위하여 상업용 소프트웨어(ANSYS)를 사용하여 계산된 결과들과 비교 되었다.

 

2. 이 론

2.1 지배 미분방정식

Fig. 1은 길이를 따라 선형적으로 폭이 감소되는 보를 도시한 것이고, 단면의 치수 보다 충분히 길이가 크다고 가정하였기 때문에 회전관성과 전단변형이 고려되지 않았다. 여기서, XYZ는 전체 좌표계이고 L은 보의 전체길이, c는 테이퍼 비이다.

Fig. 1Notation and coordinate system

해밀턴의 원리(Hamilton's Principle)로부터 테이퍼진 보에 대한 미분방정식 및 다양한 힘을 유도 할수 있다. 해밀턴의 원리의 필수요소인 변형에너지(U)와 운동에너지(T)는 다음과 같다.

여기서, EI(x)와 mI(x)는 테이퍼 비로 인해 보 요소의 길이를 따라 변화되는 굽힘 강성과 단위 길이당 질량의 변화를 표현 한 것으로 다음의 식과 같이 쓸 수 있다(1,8,10).

여기서, EI0와 m0는 c=0일 때의 굽힘 강성과 단위 길이당 질량이다.

해밀턴의 원리는 다음과 같이 표현된다.

식 (1)과 (2)를 식 (5)에 대입하고 부분적분을 통하여 지배 미분방정식, 전단력(shear force)과 굽힘 모멘트(bending moment)를 다음과 같이 결정 할 수 있다.

여기서, V는 전단력이고 M은 굽힘 모멘트이다.

식 (6)이 각주파수(ω)를 가지는 조화진동이라고 가정한다면 구조물의 진폭은 다음과 같은 형태로 표현된다.

식 (9)를 미분하여 식 (6)에 대입하고 변수를 분리한다면 미분방정식은 다음과 같이 주어진다.

여기서, 이고 이다.

2.2 미분방정식의 일반해

식 (10)의 정확한 근들은 멱급수의 하나인 Frobenius Method을 사용하여 결정할 수 있으며 일반해는 다음과 같이 쓸 수 있다.

여기서, ai+1은 멱급수의 계수이다.

식 (11)을 미분하여 식 (10)에 대입하고 정리하면 지수 k와 계수 ai+1값들을 결정 할 수 있고, 다음과 같이 멱급수의 결정방정식(indicial equation)이 얻어진다.

또한, 순환관계(recurrence relationship)에 의해 일반화된 계수들에 대한 식이 다음과 같이 결정된다.

그리고 계수 a1, a2, a3, a4, 그리고 a5의 값들은 다음과 같이 주어진다.

따라서, 식 (12)에서 주어진 지수 k의 4개에 값에 대해 미분방정식의 일반해가 다음과 같이 주어진다.

식 (19)의 함수 f(ζ ,k)는 다음과 같다.

2.3 전달행렬

식 (20)을 미분함으로서 변형곡선의 기울기(𝛷)는 다음과 같이 된다.

식 (7)과 식 (8)에 식 (11)을 미분하여 대입하면 전단력과 굽힘 모멘트는 다음과 같이 표현된다.

고정단과 자유단에서의 보 요소의 길이 x는 0과 L이 이므로 x=0을 식 (19), (21)~(23)에 대입하고 행렬 형태로 표현하였을 경우 다음과 같은 형태를 가진다.

식 (24)는 식 (25)와 같이 단수하게 표현할 수 있고,

이 식으로부터 상수 Aj는 다음과 같이 결정된다.

동일한 방법으로 x = L을 식 (19), (21)~(23)에 대입하고 행렬형태로 표현한다면 다음과 같이 된다.

식 (27)을 다음과 같이 단순하게 표현하고

식 (26)를 대입하면 폭이 선형적으로 감소된 보의 굽힘 진동에 대한 전체 전달행렬이 얻어진다.

여기서, [Tij] = [Hij] [Cij]-1 이다.

식 (29)의 전체 전달행렬에 다양한 경계조건을 적용함으로서 시스템의 고유진동수를 계산 할 수 있으나, 하나의 예로써 Cantilever 경계조건이 고려되었다. Cantilever 경계조건은 고정단에서 변위들이 0이 되고

자유단에서 힘이 0이 된다.

 

3. 결과 및 고찰

3.1 적용 예

제안된 방법의 정확성을 논리적으로 검증하기 위하여 Table 1에서 주어진 물성치가 사용되었고, 제안된 방법으로 계산된 고유진동수들이 ANSYS를 사용하여 얻어진 결과들과 비교되었다. Table 1에서 b0와 h0는 테이퍼 비가 0일 때, 다시 말해, 감소되지 않은 사각 단면의 폭과 높이 이다. 또한, 회전관성과 전단변형의 영향을 최소화하기 위해 폭과 높이에 비해 상당히 긴 구조물이 선택되었다.

Table 1Material properties

3.2 결과 및 고찰

검토된 테이퍼 비는 0부터 0.9까지 0.1간격으로 증가 되었고, Table 2에 테이퍼 비가 0일 때 계산된 고유진동수들이 나타나있고, ANSYS로부터 주어진 것들과의 비교 결과 잘 일치하는 것을 알 수 있다. 또한, Tables 3~5에 테이퍼 비가 0.1부터 0.9까지 계산된 결과들이 비교 되어있다.

Table 2The first three natural frequencies when c = 0

Table 3Effects of tapering on the first three natural frequencies when c=0.1, 0.2 and 0.3

Table 4Effects of tapering on the first three natural frequencies when c=0.4, 0.5 and 0.6

Table 5Effects of tapering on the first three natural frequencies when c=0.7, 0.8 and 0.9

Table에서 비교된 결과들을 보아서 알 수 있듯이 제안된 방법으로부터 주어진 결과와 ANSYS를 사용하여 계산된 결과가 아주 잘 일치하는 것을 알 수 있다. 하지만, 이 결과를 얻기 위하여 제안된 방법은 1개의 요소를 사용하였고, ANSYS는 평균 7370개의 요소가 사용되었다. 또한, 결과들 사이의 차이(Diff. (%))는 다음의 식으로 계산되었다.

선택된 예제를 통하여 검토된 결과들로부터 제안된 방법의 정확도가 논리적으로 증명되었고, 다양한 테이퍼 비에 대한 영향이 검토 되어 Fig. 2에 그 결과들이 보여진다. 일반적으로 높이가 테이퍼진 보는 1차 고유진동수가 증가하고, 2차 및 3차 고유진동수는 감소하는 것으로 문헌(8)에서 보고되었으나 폭이 테이퍼진 보는 모든 고유진동수가 테이퍼의 영향으로 증가 되는 것으로 확인되었고, 1차 고유진동수가 테이퍼 비의 영향이 가장 큰 것으로 확인되었다.

Fig. 2Effects of tapering on the natural frequencies of beam

 

4. 결 론

이 길이를 따라 선형적으로 감소되는 폭을 가진 보에 대해 정확한 고유진동수를 제안된 전달행렬법을 이용하여 해석하였으며 다음과 같은 결론을 얻었다.

(1) 질량행렬과 강성행렬을 가지는 기존의 전달행렬법과는 다르게 구성요소로 이들 행렬을 가지지 않는 새로운 전달행렬법을 개발하였다.

(2) 테이퍼 보에 대해 정확한 형상함수를 가정함으로서 멱급수를 통하여 미분방정식의 근을 결정하였고, 제안된 방법에 적용함으로서 1개의 요소를 사용하여 폭이 테이퍼진 보에 대해 정확한 고유진동수들을 결정할 수 있는 전달행렬법이 개발되었다.

(3) 길이를 따라 선형적으로 감소되는 폭을 가진 보에 대한 테이퍼의 영향을 조사하기 위해 테이퍼 비를 0에서 0.9까지 0.1단위로 증가시켜 그 영향을 검토하였다. 또한, 검토된 3개의 고유진동수가 모두 테이퍼의 영향으로 증가되는 것으로 확인되었고, 1차 고유진동수는 다른 고유진동수들에 비하여 그 영향이 더 큰 것으로 확인되었다.