최근 들어 부유식 구조물을 활용한 해양에너지 개발이 활발히 이루어 지고 있다. 부유식 구조물은 수면 위에서 상하동요(heaving), 횡동요(rolling), 선수동요(yawing)등의 움직임을 보이게 되며, 시추 및 생산 설비의 효율을 극대화 하기 위해 초기 설계시 6자유도 운동에 대한 고유진동수 예측이 매우 중요시 되고 있다. 부유체의 6자유도 운동중의 하나인 상하동요에 관하여 이전 연구로는 Lee and Lee (2013)의 주파수 영역 해석을 통해, 그리고 Kim and Lee (2013)는 시간 영역 해석을 통해, 자유수면상에서 상하동요하는 원주의 무차원 고유진동수 Ω₁을 MCFR(Modulus of Complex Frequency Response)의 국부 최대값으로 정의할 수 있음을 보였다.

본 연구에서는 자유수면 상에서 상하동요하는 2차원 주상체의 고유진동수를 주파수 영역 해석을 통해 얻었으며, 원주를 포함한 여러 가지 형상을 주상체의 형상으로는 Lewis 형상을 사용하였다. 자유 상하동요하는 주상체의 운동방정식은 다음과 같이 쓸수 있다.

여기서 m_a는 부가 질량, μ는 조파 감쇠이며, C는 복원력상수이다. 부가 질량과 조파 감쇠를 얻기 위해서는 Tasai (1959)의 Ursell-Tasai 방법을 사용하였다. 차후의 시간 영역 해석을 고려하여 진동수가 매우 큰 영역에서 대해서도 수치 계산을 수행하였으며, 기존의 결과들과 비교하여 매우 큰 진동수들에 대해서도 수치적으로 안정된 결과를 얻었다. 이 결과를 사용하여, 실제 선박에 사용되는 형상들의 범위를 포함하는 폭-흘수 비 H와 단면적 계수 σ의 범위에 대해 고유진동수 Ω_n 을 얻었다.

한편 Lee (2008)는 부가 질량만을 고려한 근사적인 고유진동수 Ω_a를 구하는 방법을 제안한 바 있는데, 이 방법에 의한 근사값도 구하여 위에서 정의된 고유진동수 Ω_n 의 값들과 비교하였다. 또한 통상 많이 사용되는 조화감쇠 동요 근사에 의한 고유진동수 Ω_h 도 구하여 같이 비교하였으며, 일반적으로 Ω_a > Ω_n > Ω_h 의 결과를 얻었다. 계산된 범위 내에서는 H가 감소하거나, σ가 증가하면, 고유진동수 Ω_n 과 이 진동수에서의 MCFR의 값은 증가하는 경향을 보였다. 본 연구에서 얻어진 결과를 토대로 실제 선박에 대한 고유진동수 추정 또한 스트립 방법 등을 사용하여 가능할 것으로 판단된다. 2절, 3절에서는 문제의 정식화와 그 해법에 대해 간단히 설명하고, 4절에서는 계산 결과에 대한 토론, 5절에서는 결론을 맺기로 한다.

연구에서의 좌표계는 Fig. 1 에 나타낸 바와 같으며, 물리적 평면에 대해 ( x, y)와 Lewis 변환을 위한 수학적 평면에 대해( ξ, η )를 사용하였고, 주상체의 반폭은 a, 흘수는 d이다.

본 연구에서는 이상 유체의 비회전성 유동장을 가정하며, 대칭성에 기인하여 제 1사분면만을 고려하면 충분하다. 주상체의 상하동요는 다음 식으로 주어지는데,

여기서 l은 상하동요의 진폭이고, ω는 진동수, ε은 방사파와의 위상차로 미지량이다.

자유수면 상에서 상하동요하는 2차원 주상체에 대한 경계치 문제는, 비압축성 유체에 대한 연속방정식, 선형화된 자유수면 경계조건, 수저조건, 물체 경계조건 및 방사조건으로 구성되며 속도 포텐셜 ∅ ( x, y, t )에 대해 다음과 같다.

여기서 K = ω²/g으로 파수이다.

여기서 ζ = ξ + ie^αe^-iθ로 쓸 수 있으며, 따라서 α = 0는 단위 원의 원주에 상응한다.

마지막으로 주상체의 상하동요에 의해 생성된 파도는 주상체로부터 먼 쪽으로 전파된다는 방사조건을 만족해야 한다.

본 연구에서 사용한 Ursell-Tasai 방법은 매우 잘 알려져 있으므로, 자세한 사항은 Song (2014)으로 미루기로 하고, 다음 절들에서 사용할 물리량들에 대한 혼동을 피하기 위해 정의 위주로 간략히 다루기로 한다.

Lewis 변환은 다음 식으로 주어진다.

여기서 M_f는 배율 계수(magnification factor)이고, a₁ , a₃는 단면형상 계수(sectional form parameter)이다. 따라서 다음 관계를 얻으며,

주상체의 반폭 a와 흘수 d를 다음과 같이 얻는다.

또한 주상체의 단면적 계수(sectional area coefficient) σ와 폭-흘수 비(beam-draft ratio) H는 다음과 같이 주어진다.

식 (13)에서의 A_s는 주상체의 단면적이다.

Ursell (1949)에 따르면 속도 포텐셜은 다음과 같다.

여기서 b는 주상체의 상하동요에 의해 생성 전파되는 물체로부터 상당히 떨어진 곳에서의 진행파, 즉 방사파의 진폭이다. p_2n과 q_2n은 미정으로 물체 경계조건을 만족하도록 결정되어야하며, ∅_c와 ∅_s는 진행파, R_2n은 국부교란 유동을 나타낸다. 또한 식 (15)의 속도 포텐셜에 공액(conjugate)인 유량 함수는 다음과 같다.

여기서 ψ_c와 ψ_s는 각각 ∅_c 와 ∅_s의 공액이다. κ = Ka라고 하면 유량 함수를 사용한 물체 경계조건은 다음과 같다.

여기서 α_L(θ)는 다음과 같은 함수이고,

ci , si 는 각각 코싸인 적분, 싸인 적분이다.

한편 방사파의 진폭과 상하동요하는 2차원 주상체의 진폭의 비, b/l은 다음과 같이 주어지는데,

여기서 κ_L = κ/ ( 1 +a₁ + a₃ )라고 하면, A와 B는 각각 다음과 같다.

식 (17), (18)을 풀어 p_2n과 q_2n을 구하면, 비정상 유동에 대한 베르누이 방정식을 사용하여 물체표면에서의 압력을 구할 수 있고, 압력의 연직방향 성분을 적분하여 부가 질량과 조파 감쇠를 얻을 수 있다. 먼저 물체표면 상에서의 압력은 다음과 같고,

여기서 하첨자 S는 물체표면을 뜻하며, M, N 또한 물체표면에서 계산된 값이다. 부가 질량과 조파 감쇠는 다음과 같이 얻을 수 있으며,

여기서 M₀, N₀는 다음과 같다.

또한 수치 계산의 정확성을 검사하기 위해 다음 식을 이용할 수 있다.

따라서 부가 질량계수와 조파 감쇠계수는 다음과 같이 주어진다.

위의 식에서 a₁ = a₃ = 0인 경우, 원주에 대한 결과를 얻는 것을 확인할 수 있다.

상하동요하는 주상체의 부가 질량과 조파 감쇠를 고려하지 않은 경우에 대한 원고유진동수(pre-natural frequency) ω₀와 그에 상응하는 무차원 파수 κ₀를 다음과 같이 정의한다.

Lee (2008)는 부가 질량만을 고려한 경우의 근사적인 고유진동수 ω_a를 다음 식의 해로 정의하였다.

이하에서는 모든 진동수를 원고유진동수 ω₀ 로 무차원화하고, Ω_a와 같이 나타내기로 한다.

식 (1)로부터 무차원화된 MCFR, R^* ( Ω )를 다음과 같이 얻을 수 있는데,

여기서 이며, 일 때의 무차원 진동수 Ω_n을 무차원 고유진동수로 정의한다.

위에서 정의한 무차원 진동수들을 표로 정리하면 Table 1과 같다.

실제 선박의 횡단면 형상을 고려하여, 폭-흘수 비 H = 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 또 단면적 계수 σ = 0.6, 0.7, 0.8, 0.9에 대해, Lee and Lee (2013)에 의해 개발된 공치법(collocation method)을 사용하여 수치 계산을 수행하였다.

계산 결과의 정확성을 검증하기 위해 식 (28)을 사용하여 오차 E를 다음과 같이 정의하고,

비교적 오차가 큰 경우인 H = 3.5, σ = 0.8 인 형상에 대해 공치법에서의 미정 계수의 수인 N = 10, 30, 50에 대한 오차를 Fig. 2에 보였다. 여기서 N = 50에 대한 수치 결과는 κ <에 대해 1.0% 이하의 오차를 주고 있음을 확인할 수 있으며, 본 연구의 계산에서는 N = 50으로 취한 공치법을 사용하였다.

속도 포텐셜의 수렴을 확인하기 위해 H = 3.5, σ =0.8인 경우에 대해 p_2n과 q_2n의 변화를 Fig. 3, Fig. 4에 보였다.

여기서 수평축, 이며, 을 기준으로 좌우의 수직축의 크기를 다르게 나타내었다. 식 (17)과 (18)의 급수에서 계수 p₂와 q₂가 지배적임을 확인할 수 있으며, 또한 각각의 계수들이 κ에 따라 진폭이 서서히 감소하는 조화함수의 형태를 보이고 있고, 이는 원주에 대한 결과와 본질적으로 유사함을 확인하였다.

각각의 형상에 대해 b/l은 κ의 함수이며, 국부 최대값을 가지는데, σ가 변화함에 따른 b/l의 최대값을 Fig. 5에 보였다. H가 증가할수록, 또 σ가 감소할수록 각 형상에 대한 b/l의 국부 최대값이 증가함을 알 수 있으며, 이는 폭이 크고, 또 수면 부근에서의 형상 변화가 클수록 b/l의 최대값이 증가하는 것을 뜻한다.

Tasai (1959)는 K₄ ( κ )를 다음과 같이 정의하였고,

K₄ ( κ )에 대한 Tasai의 결과와 비교하여 Fig. 6, 7에 보였다. Tasai의 결과보다 훨씬 넓은 κ의 범위에 대해 충분히 만족할 만한 정확성을 가지는 결과를 얻고 있음을 확인할 수 있다.

각 경우에 대해 부가 질량 계수 C_a ( κ )는 국부 최소값을 가지며, 이들 최소값의 각 H 에 대해 σ에 따른 변화 경향을 최소 값에 상응하는 κ의 값, 즉 κ_min 의 함수로 Fig. 8에 보였다. C_a의 최소값과 그에 상응하는 κ의 값, 즉 κ_min은 H가 증가할수록, 또 H ≥ 2.0인 경우에는 σ가 감소할수록 커지는 것을 확인할 수 있으며, 원주에 대한 C_a의 최소값은 H ≥ 2.0인 경우에 대해서는 모든 형상 중 가장 작은 값을 가짐에 유의한다.

Lee (2008)는 부가 질량만을 고려한 근사적인 고유진동수 ω_a를 다음과 같은 반복자(iterator)를 통해 정의하였다.

위 식을 이용하여 얻은 무차원 고유진동수 Ω_a 를 각 H에 대해 σ의 함수로 Fig. 9에 보였다.

Ω_a 는 H가 감소할수록, 또 H ≥ 2.0인 경우에는 σ가 증가할수록 커지는 것을 확인할 수 있다. Lee and Lee (2013)는 2차원 원주에 대해 Ω_a=0.7919를 얻었고, 이는 Ω_n =0.7970과 비교하여 매우 좋은 근사임을 보였는데, 위의 그림에서 원주에 대한 Ω_a 는 H ≥ 2.0 인 모든 형상 중 최대값을 가짐을 알 수 있다.

각 경우에 대해 부가 질량 계수 C_a가 최소값을 가지는 진동수 Ω_m을 Fig. 9에서와 같은 방식으로 Fig. 10에 보였다. Ω_a와 Ω_m이 전반적으로 매우 유사한 경향을 나타냄을 확인할 수 있다.

Ω_a 와 Ω_m 의 차이를 %로 Fig. 11에 보였다. Ω_a 와 Ω_m 의 값은 H ≥ 1.0인 모든 형상에 대해 2% 이하의 차이를 보이며, Ω_a 는 Ω_m으로 근사하여도 큰 오차를 범하지 않음을 알수 있고, 또 대체적으로 σ가 감소하면 오차가 증가함을 확인할수 있다. Lee, et al. (2013)에 따르면 Ω_a ≃ Ω_m 이기 위한 조건은 인데, 예를 들어 H =5.0, σ = 0.6인 경우, κ_min =2.4, C_{a min} =2.2이고, = 0.12κ_min 이며, 따라서 δ=0.0245를 얻어, 위의 조건이 만족되는 것을 확인할 수 있다.

식 (33)의 MCFR의 국부 최대값에 상응하는 진동수, 즉 본 연구에서 상하동요의 고유진동수로 정의한 Ω_n 을 각 형상에 대해 Fig. 9에서와 같은 방식으로 Fig. 12에 보였다. 전반적으로는 Ω_a 및 Ω_m 과 비슷한 경향, 즉 Ω_n 은 H 가 감소할수록, 또 H ≥ 2.0인 경우에는 σ 가 증가할수록 큰 값을 가지는 것을 볼 수 있으며, Lee and Lee (2013)가 얻은 2차원 원주에 대한 값인 Ω_n =0.7970이 H ≥ 2.0인 모든 형상 중 최대값임을 알 수 있다.

한편, 식 (33)의 R^*가 국부 최대값을 가지는 Ω_n에서의 부가 질량과 조파 감쇠를 사용하면, 식 (1)에서 보인 운동 방정식의 계수가 상수인 미분방정식으로 변화하며, 이와 같은 근사는 조화 감쇠 동요로 볼 수 있고, 이에 대한 MCFR 는 다음과 같다.

여기서 C_an = C_a ( Ω_n ) , C_μn = C_μ ( Ω_n )이고, 수치 계산 결과를 사용하여 국부 최대값을 가지는 진동수 Ω_h 를 구할 수 있으며, 이 결과를 Fig. 13에 보였다.

일반적으로 Ω_h 는 Ω_n보다 작은 값을 가지며, 각각의 단면형상에 대한 는 R^* (Ω) 와는 최대값의 위치, 최대값의 크기 등이 모두 다른 특징을 가진다.

Ω_a와 Ω_n의 차이, Ω_h와 Ω_n의 차이를 %로 각각 Fig.14, 15에 보였다. 일반적으로 Ω_a > Ω_n > Ω_h 의 결과를 보여 주고 있다. Ω_h가 증가할수록 또 σ가 감소할수록 그 차이는 증가 한다. 이는 조파 감쇠가 Ω_h가 증가할수록 또 σ가 감소할수록 증 가하는 경향과 일치하며, 조파 감쇠가 커지면 고유진동수를 결정 함에 있어 조파 감쇠를 무시할 수 없음을 알 수 있다.

주상체의 중간 형태인 H =3.5, σ =0.8인 형상에 대한 결과를 기준으로 잡아 H 와 σ가 각각 변화할 때, Ω_n의 값과 그에 상응하는 MCFR R^* 의 값이 어떻게 변화하는지를 Fig. 16, 17에 각각 보였다. H가 감소함에 따라 Ω_n 과 그에 상응하는 R^* 의 값이 증가함을 알 수 있으며, 이는 C_amin 이 H가 감소함에 따라 감소하는 경향과 부합하는 결과이다. 또한 σ 가 증가함에 따라 Ω_n과 그에 상응하는 R^*의 값이 증가함을 알 수 있으며, 이는 C_amin 이 σ가 증가함에 따라 감소하는 경향과 부합하는 결과이다(Fig. 8 참조).

본 연구에서는 상하동요하는 2차원 주상체의 고유진동수를 예측하기 위해서 Ursell (1949)과 Tasai (1959)가 2차원 주상체에 적용한 방법을 기본으로, Lewis 변환을 사용해 보다 광범위한 형태를 가지는 주상체의 부가 질량 계수와 조파 감쇠 계수를 공치법에 의한 수치 계산으로 얻었다. 이렇게 얻은 결과를 이용하고 MCFR을 사용하여 2차원 주상체의 고유진동수를 정의함으로써 다음과 같은 결론을 얻었다.

H ≥ 2.0 인 모든 형상에 대해 부가 질량과 조파 감쇠를 모두 고려하여 정의한 고유진동수 Ω_n 은 원고유진동수의 80%보다 작으며, 그 최대값은 원주의 고유진동수로 간주할 수 있음을 확인하였다.

차후의 시간 영역 해석을 위해 매우 큰 κ에 대해서도 조파 감쇠 계수와 부가 질량 계수를 안정되게 계산할 수 있는 수치 코드를 개발하였다.

부가 질량만을 고려한 고유진동수 Ω_a 는 C_amin 에 상응하는 진동수 Ω_m 과 2% 이하의 차이를 가진다.

MCFR의 국부 최대값에 상응하는 고유진동수 Ω_n 은 계산된 범위에서는 대체적으로 Ω_a 보다는 작고 조화감쇠 동요 근사에 의한 Ω_h 보다는 큰 값을 가진다.

주상체의 H가 감소할수록, σ가 증가할수록 고유진동수 Ω_n과 그에 상응하는 MCFR의 크기 R^* 의 값은 증가하며, 이는 C_amin 의 H, σ에 대한 경향과 부합하는 것으로 보인다.

부유체의 운동에 대해 strip 방법 또는 그와 유사한 방법을 사용하는 경우 본 연구에서 사용된 방법은 3차원 물체에 대해서도 큰 어려움 없이 확장될 수 있을 것으로 생각되며, 초기 설계 시부유체의 상하동요 고유진동수 예측에 큰 도움이 될 것으로 생각된다.