0044-01-0017-0006-2.pdf2.73MB

1. 서 론

 지구 온난화를 방지하고 기후변화에 대응하기 위하여 신재생에너지에 대한 관심이 전세계적으로 급증하고 있다. 다양한 종류의 신재생에너지 중에서도 풍력발전은 최근 상용화가 이루어져 전세계적으로 수많은 풍력발전 단지가 건설되었고 또 건설되고 있다. 풍력발전은 발전기가 설치되는 위치에 따라 육상풍력발전과 해상풍력발전으로 나누어진다. 해상풍력발전은 육상풍력발전에 비하여 용지 확보가 용이하고, 주변 환경에 대한 소음의 영향에서 어느 정도 자유로울 수 있다는 장점을 가지고 있다. 또한 육상에 비하여 해상에서는 지속적으로 고품질의 바람을 얻을 수 있다. 그러므로 현재까지 육상풍력발전 위주로 진행되어 온 풍력발전은 향후 해상풍력발전중심으로 건설되고 운영될 것으로 예상된다. 하지만 해상풍력발전기는 그 건설 비용이 육상풍력발전기보다 높으므로 비용 효율적으로 발전 시설을 건설하고 운영하는 방안을 모색해야 하고 이를 위해서는 다양한 하중 조건하에서 구조물의 거동을 정확하게 예측할 수 있어야 한다[1-4].

 과거 기록에 의하면 지진에 의해서 구조물이 손상되거나 붕괴하면 인명과 재산에 많은 피해가 초래되었다. 만약 해상풍력발전기가 지진으로 인해서 심각한 손상을 입거나 파손된다면 환경에도 큰 피해를 줄 수 있을 뿐 아니라 운영과 에너지 공급에 큰 차질을 빚을 수 있으며 막대한 복구비용을 유발시킬 수 있다. 그러므로 해상풍력발전기는 지진에 대해서 일반 구조물보다는 한층 더 높은 수준의 안전성을 보장할 수 있도록 설계되고 유지되어야 할 필요가 있다.

 해상풍력발전기의 지진 거동을 실무적으로 정확하게 예측하기 위해서는 복잡하지 않으면서도 중요한 역학적 관계를 잘 표현할 수 있는 해석기법이 확보되어야 한다. Lee et al.[5]은 실제 해양구조물의 내진설계와 내진성능평가에 직접 활용할 수 있고, 한층 더 복잡한 문제에도 용이하게 확대하여 적용할 수 있는 유체-구조물-지반 상호작용을 고려한 해상풍력발전기의 지진응답 해석기법을 개발하였다. 이 연구에서는 이미 개발된 해석기법의 역량을 확장하여 구조물의 기하 비선형과 해수의 항력을 고려할 수 있도록 하였다.

 일반적으로 풍력발전기의 tower는 세장비가 큰 구조물이므로, 이러한 구조물을 설계할 때는 기하 비선형 등의 영향을 고려하여야 한다. 이를 고려하기 위해 이 연구에서는 구조물의 변형과 기초의 회전운동에 의해 발생하는 기하 강성을 유도하여 기하 비선형의 효과를 고려하였다.

 해양구조물에 가해지는 동수압 또는 동수압력은 다양한 방법에 의하여 산정될 수 있다. 일반적으로 가장 많이 사용되는 방법은 Morison 방정식[6]을 이용하는 것이다. Morison 방정식은 동수압의 영향을 관성력에 의한 부분과 항력에 의한 부분으로 나누어 고려하는데, 얇은(slender) 구조물에는 잘 적용할 수 있는 것으로 알려져 있다[7]. 하지만, 입사파의 파장이 짧아져서 유체의 회절과 방사 감쇠의 영향이 중요해지게 되면 Morison 방정식을 적용하기는 힘들어지게 된다. 이러한 경우에는 acoustic wave equation에 기반한 회절이론으로부터 해양구조물에 가해지는 동수압을 엄밀히 산정하여야 한다. Lee et al.[5]은 유체의 회절, 방사 감쇠, 유체의 압축성을 엄밀히 고려하여 동수압을 산정하였고, 이를 활용하여 해상풍력발전기의 지진응답 해석기법을 개발하였다. 하지만, 이 방법에 의해서는 유체의 점성 등으로 인해 발생하는 항력의 효과는 고려할 수 없고 항력의 효과가 우세한 얇은 구조물에 적용하는 것은 한계가 있으므로, 이를 고려하기 위해서는 Navier-Stokes 방정식에 기반한 전산유체역학(computational fluid dynamics)이론을 사용하여야 하지만 이는 많은 계산량을 요구하기 때문에 적용에 한계가 있다. 이 연구에서는 이상과 같은 단점을 보완하고자 회절이론과 Morison 방정식을 결합하여 해상풍력발전기 tower에 가해지는 동수압을 산정하고자 하였다[8]. 즉, 유체의 회절, 방사 감쇠, 유체의 압축성을 엄밀히 고려하여 회절이론으로부터 산정한 동수압과 Morison 방정식에서 항력에 의한 동수압력을 결합하여 지진에 의해 발생하는 해수의 동수압력을 산정하였다.

 이상과 같이 이 연구에서는 기존에 개발된 해석기법을 확장하여 구조물의 기하 비선형과 해수의 항력을 고려한 해상풍력발전기의 지진응답 해석기법을 개발하였다. 개발된 해석기법을 사용하여 해상풍력발전기와 같은 복잡한 유체-구조물-지반 상호작용계의 지진응답에 영향을 미치는 인자를 파악하고 그 인자들의 영향의 유의성을 조사해 보고자 한다.

2. 유체-구조물-지반 상호작용계

 Fig. 1은 일반적인 해상풍력발전 시스템을 보여주고 있다. Fig. 1(a)에 보인 바와 같이 해상풍력발전 시스템은 유연한 구조물, 해수, 해수로 포화된 지반이 서로 상호작용하는 유체-구조물-지반 상호작용계이다. 해수의 동수압은 tower와의 접촉면에 작용하여 tower의 동적 거동에 영향을 주고, 이 영향을 받은 tower의 동적 거동은 다시 해수의 동수압에 영향을 주게 된다. 유연한 지반과 tower의 동적 거동도 foundation을 통해 서로 영향을 주고 받게 된다. 이러한 유체-구조물-지반 상호작용계는 Fig. 1(b)와 같이 3개의 부분구조로 분리하여 모델링 하는 것이 가능하다. 각 부분구조의 모델링 방법은 Lee et al.[5]에 자세히 설명되어 있으므로 여기서는 주요한 특징들만 간략히 설명하고자 한다.

Fig. 1. Offshore wind turbine system

2.1 구조물 모델링

 해상풍력발전기의 지지구조인 tower는 축대칭 캔틸레버 보로 가정하였다. Tower의 길이(높이)는 L이고, 중공원형 단면으로서 바깥 반지름과 두께는 각각 R과 h로 높이에 따라서 일정하며, 보의 단위 길이당 질량 m과 휨강성 EI도 높이에 따라서 일정하다고 전제한다. 상부구조인 Rotor와 power train을 포함한 nacelle은 tower의 정점에 부착된 집중질량으로 모델링 되었는데, 그 질량과 회전관성은 M_t 와 I_t이다. Foundation은 축대칭형상의 강체기초로 가정하고, 기초의 질량과 회전관성은 M_b와 I_b이다.

 구조물에는 해수의 동수압 p(r, z,t)와 유연한 지반과의 상호작용력 와 가 작용한다. Tower의 분포질량은 양단에 균등하게 집중된다고 볼 수 있고 이러한 가정에서 축력이 계산될 수 있다. 따라서 상부의 집중질량과 tower의 분포질량으로 인해 의 축력이 빔의 중간 높이에서 발생한다. Fig. 2에 보인 바와 같이, 구조물의 변위는 강체기초의 수평운동에 의한 변위 Δ(t), 회전운동에 의한 변위 zφ (t), 이에 대한 상대변위 u(z,t)의 합으로 표현할 수 있다. Tower의 상대변위 u(z,t)는 다음과 같이 캔틸레버 보의 고유모드의 선형조합으로 표현할 수 있다[9].

Fig. 2. Tower Model

 여기서 ω_i는 캔틸레버 보의 i번째 고유진동수이다

2.2 Tower에 작용하는 동수압

 해수는 비점성의 압축성 이상 유체로 모델링 하였고 그 깊이는 H로 일정하며 반경방향으로 무한하다고 가정하였다. 지진 시에는 3축 방향 운동성분이 모두 작용하나 이 연구에서는 일축으로 작용하는 수평지반운동만 고려한다. 이때 원통형 좌표계에서의 동수압은 원주방향으로 cosθ 의 분포를 가지게 되고, 이러한 분포를 가지는 해수의 동수압의 진폭 p(r, z,ω )는 다음과 같이 표현된다[5].

 여기서 는 제2종 제1차 Hankel 함수이고, k_n은 양의 실수이거나 음의 허수부를 가지는 복소수이어야 한다. 식(2d)의 C_w 는 해수에서의 압축파의 전파속도이다.

2.3 지반-구조물 상호작용력

 구조물이 위치하는 지반은 해수로 포화된 다공성 2상 매질로 가정하고, 구조물과 지반의 상호작용을 모델링하였다. 지반의 일정 깊이부터는 아주 단단한 기반암이 존재한다고 가정하였다. 강체 기반암 조건은 continuedfraction absorbing boundary condition 등을 활용하여 완화할 수 있다[10, 11].

 Fig. 3에 묘사된 바와 같이 불규칙하고 비균질인 근역 지반은 유한요소로 모델링하고, 깊이가 일정하고 수평방향으로 균질인 원역 지반은 전달경계를 사용하여 모델링하면[10-13], 이러한 지반에 놓인 무질량 강체기초의 운동방정식과 상부구조물에 가해지는 상호작용력을 다음과 같이 얻을 수 있다[14].

Fig. 3. Soil-structure interaction

 여기서, , i, j = Δ,φ 는 기초의 동적강성이고, Δ⁰와 φ⁰는 각각 기초의 수평운동과 회전운동에 대한 입력운동을 의미한다. 기초 입력운동의 위첨자 0는 상부구조물이 없고 기초만 있는 경우를 의미하고, 기초의 입력운동은 어떠한 외력도 작용하지 않을 때 입사 지진파에 대한 기초의 응답을 의미한다[14].

3. 유체-구조물-지반 상호작용계의 운동방정식

 자중으로 인한 기하 비선형의 효과와 tower와 해수의 상대속도로 인한 항력의 효과를 고려한 유체-구조물-지반 상호작용계의 운동방정식을 유도하였다. 상부구조에서 발생하는 운동에너지, 변형에너지, 감쇠에 의해 소산되는 에너지, 외력의 일함수(work function)를 Lagrange equation에 적용하여 운동방정식을 구한다.

 기초와 구조물에 발생하는 운동에너지는 다음과 같다.

 Tower에 저장되는 변형에너지는 다음과 같다.

 Tower에 작용하는 동수압, 기초에 작용하는 상호작용력, 자중으로 인해 tower에 발생하는 축력, tower와 해수의 상대속도로 인한 항력의 일함수는 다음과 같다.

 여기서 항력은 Morison 방정식을 이용하여 다음과 같이 선형화하여 표현할 수 있다[7].

 여기서 D = 2R은 tower의 지름, ≈ constant는 해수의 선형화된 항력계수이고, C_D 는 해수의 항력계수이다. 또한, 식(6)에서 L(1-cos[φ(t)])와 dz는 각각 기초의 회전과 tower의 변형으로 인해 tower의 최상단이 수직 방향으로 이동하는 거리이고, Taylor series를 사용하여 L(1− cos[φ (t)]) ≈ [φ(t)]²으로 근사할 수 있다. 이와 같이 근사하면 식(6)은 다음과 같이 표현된다.

 각 일반화된 변위에 대한 Lagrange equation은 다음과 같이 주어진다.

 구조물의 고유모드 중 처음 N개만 고려하고, 식(4), (5), (8)을 식(9)에 대입한 후, 식(3)으로부터 상호작용력 과 의 표현을 활용하면 다음을 얻을 수 있다.

 여기서

 식(10)과 (11)에서 확인할 수 있듯이 식(8)에서 자중의 효과를 고려한 부분은 기하 강성의 형태로 표현되게 된다. 식(2)의 동수압과 식(7)의 항력을 식(10)에 대입하면 최종 운동방정식을 얻을 수 있다.

 여기서

 구조물의 감쇠를 이력감쇠로 가정하고, 식(12)를 행렬의 형태로 표현하면 다음을 얻을 수 있다.

 여기서

 식(15)에서 ξ 는 구조물의 이력감쇠에 대한 등가 점성감쇠비이다. 식 (15)가 해수와 해수로 포화된 지반과의 상호작용, 해수 항력, 구조물 기하비선형의 효과를 고려한 해상풍력발전기 모델의 최종 운동방정식이다.

4. 적용 예제

 개발된 해석기법을 이용하여 Fig. 4에 보인 해상풍력발전기의 지진응답을 계산하였다. 해석 대상은 National Renewable Energy Laboratory(NREL)에서 제시한 5 MW급 해상풍력발전기[15, 16] 를 일부 수정한 형태로서, tower, nacelle 등을 포함한 상단구조물, 해수, 지반의 물성치는 Table 1과 Fig. 4에 주어져 있다. Table 1에 주어진 rotor와 power train을 포함한 nacelle의 회전관성 I_t는 rotor의 회전축과 평행한 축에 대한 회전관성이다. NREL에서 대상으로 한 해상풍력발전기는 monopole 형태의 기초를 가지고 있지만, 이 연구에서는 Fig. 4와 같이 단순화된 gravity-type foundation을 가정하였다. 하지만, 복잡한 형태의 다양한 기초 형식에 대해서도 동일한 방법론을 적용 가능하다[17]. 3차원 Cantilever의 고유모드 중 1 ~ 7차 모드를 사용하여 Tower의 응답을 식(1)과 같이 표현하였고, 해수의 동수압에 대해서는 1 ~ 40차 모드까지 사용하여 식(2)와 같이 표현하였다. 또한, 4절점 유한요소를 사용하여 Fig. 4에 보인 것과 같이 지반의 유한요소망을 구성하였다. 입력지반운동 은 Fig. 5의 El Centro 지진기록(1940년, N-S성분, PGA=0.319g)을 사용하였다.

Fig. 4. 5 MW offshore wind turbine

Table 1. Properties of 5 MW Offshore Wind Turbine

Fig. 5. Input ground motion: El Centro earthquake(1940, N-S component, PGA=0.319g)

4.1 유체-구조물 상호작용의 영향

 유체-구조물 상호작용이 해상풍력발전기의 거동에 미치는 영향을 조사하기 위해, 해수가 없는 경우(case ‘dry’)와 해수로 인해 유체-구조물 상호작용이 발생하는 경우(case ‘FSI’)에 강체 지반에 설치된 해상풍력발전기의 지진응답을 계산하였다. 해수의 영향을 고려할 때 식(7)의 선형화된 항력계수 C′_D 는 반복계산을 수행하여 그 값을 0.0191로 설정하였다. 앞서 언급한 방법으로 기하 비선형의 효과를 고려하였다.

 지진응답 해석을 수행하여 tower 최상단의 변위, tower 밑면에서 휨모멘트의 전달함수와 시간이력을 Fig. 6과 7에 비교하였다. 해석 결과로부터 유체-구조물 상호작용으로 인해 구조물의 응답은 상당한 영향을 받을 수 있음을 확인할 수 있다. 이 적용예제에서는 밑면에서의 휨모멘트가 크게 증가하였는데 이는 다음과 같이 설명할 수 있다. 유체-구조물 상호작용의 영향은 식(13)의 부가질량의 형태로 나타나는데, 그 상대적인 영향은 질량 M_ij + (i,j = 1,…, N)에 대한 부가질량 ( i, j = 1,…,N )의 비로 표현된다. 단, 이때 식(13)의 부가질량은 복소수이므로 그 실수부만을 고려하여야 한다. Fig. 8은 이 부가질량의 상대비 /(M_ii + )를 나타내는데, 여기서 확인할 수 있듯이 유체-구조물 상호작용의 상대적인 영향은 저차 진동모드보다는 고차 진동모드에서 더 크게 나타난다. 단, 여기서 주의할 점은 Fig. 8과 같은 부가질량의 상대비 변화가 모든 유체-구조물 상호작용계에서 발생하는 것은 아니므로, 일반적인 현상으로 간주하여서는 안 된다. Fig. 8에서 관찰된 사항은 Fig. 6과 7의 전달함수에서도 확인할 수 있는데, 유체-구조물 상호작용으로 인해 1차 고유모드의 변화는 거의 없지만 고차 고유모드는 그렇지 않다. 그런데, 입력지반운동의 변위는 가속도보다 고진동수 성분의 영향이 감소하게 된다. Fig. 6은 지반 변위에 대한 구조물 최상단 변위의 전달함수이고, Fig. 7은 지반 가속도에 대한 구조물 밑면에서의 휨모멘트의 전달함수이므로, 결국 구조물 최상단 변위보다는 구조물 밑면에서의 휨모멘트에 대하여 유체-구조물 상호작용의 영향이 증가하게 되는 것이다.

Fig. 6. Displacements at the top of the tower: effects of the fluidstructure interaction

Fig. 7. Base overturning moments of the tower: effects of the fluidstructure interaction

Fig. 8. Ratio of added mass

 이상과 같이 유체-구조물 상호작용으로 인해 해상풍력발전기의 지진응답이 증가하게 되고, 특히 밑면에서의 휨모멘트가 크게 증가하게 된다. 이는 시스템의 요구 설계력이 증가된다는 것을 의미하고 안전한 구조물의 설계를 위해서는 유체-구조물 상호작용을 고려해서 그 영향을 검증하여야 함을 입증하고 있다.

4.2 지반-구조물 상호작용의 영향

 지반-구조물 상호작용이 해상풍력발전기의 거동에 미치는 영향을 조사하기 위해, 유연한 지반이 아닌 강체 지반에 구조물이 설치된 경우(case ‘FSI’)와 유연한 지반에 구조물이 설치된 경우 (case ‘FSSI’)에 해상풍력 발전기의 지진응답을 계산하였다. 유연한 지반의 전단파 속도는 200 m/sec로 가정하였고, 두 경우 모두 유체-구조물 상호작용의 영향을 고려하였다. 선형화된 항력계수 C′_D 는 그 값을 각각 0.0191과 0.0351로 설정하였다. 기하 비선형의 효과를 고려하였다.

 지진응답 해석을 수행하여 tower 최상단의 변위, tower 밑면에서 휨모멘트의 전달함수와 시간이력을 Fig. 9와 10에 비교하였다. Fig. 9와 10의 전달함수를 살펴보면 시스템의 1차 고유진동수는 지반-구조물 상호작용의 영향에 의해서 크게 영향을 받지 않았음을 확인할 수 있다. 이 예제의 구조물은 1차 고유진동수가 0.113Hz인 아주 유연한 구조물이지만, 이 구조물이 설치된 지반의 고유진동수는 0.5Hz이다. 이와 같은 상호작용계의 고유진동수와 같은 특성은 지반-구조물 상호작용의 효과로 인해 크게 영향을 받지 않는다. 하지만, 기반암에 작용한 입력지반운동은 유연한 지반을 전파하면서 크게 증폭되게 되고 이로 인해 상호작용계의 지진응답은 크게 영향을 받게 된다. Fig. 9와 10의 전달함수에서 0.5Hz 성분이 크게 증폭된 것은 이러한 지반 증폭을 반증하는 것이다. 그리고, 이로 인해 해상풍력발전기의 지진응답은 크게 증가하게 되는데, 이는 Fig. 9와 10의 시간이력에서도 확인할 수 있다.

Fig. 9. Displacements at the top of the tower: effects of the soilstructure interaction

Fig. 10. Base overturning moments of the tower: effects of the soilstructure interaction

 이상의 해석결과로부터 유체-구조물-지반 상호작용은 해상풍력발전기의 지진응답을 크게 증가시킬 수 있음을 확인할 수 있다. 즉, 해상에 건설되는 해상풍력발전 시스템의 정확한 동적 거동을 예측하기 위해서는 유체-구조물-지반 상호작용의 영향을 엄밀히 고려하여야 한다.

4.3 지반 강성의 영향

 앞에서 살펴보았듯이 지반-구조물 상호작용은 해상풍력발전기의 거동에 많은 영향을 미친다. 여기서는 지반의 강성에 따라 지반-구조물 상호작용의 효과가 어느 정도인지 살펴보기 위해, 지반의 전단파 속도를 200 m/sec, 500 m/sec, 1000 m/sec로 변화시켜 가며 그 영향을 살펴보았다. 이 경우에도 모두 유체-구조물 상호작용과 기하 비선형의 효과를 고려하였고, 선형화된 항력계수 C′_D 는 그 값을 각각 0.0351, 0.0344, 0.0228로 설정하였다.

 Fig. 11과 12는 각각 tower 최상단의 변위, tower 밑면에서 휨모멘트의 전달함수와 시간이력을 보여주고 있다. 유연한 지반의 고유 진동수는 전단파 속도에 따라 달라지므로, 기반암에 작용한 입력지반운동이 유연한 지반을 전파하면서 증폭되는 특성도 달라지게 된다. Fig. 11과 12의 전달함수는 이러한 변화를 보여주고 있다. 특히, 지반의 전단파 속도가 200 m/sec인 경우에는 지반의 증폭으로 인해 지진응답이 크게 증가하지만, 지반의 전단파속도가 커질수록 이러한 증폭 현상은 작아짐을 확인할 수 있다.

Fig. 11. Displacements at the top of the tower: effects of the soil stiffness

Fig. 12. Base overturning moments of the tower: effects of the soil stiffness

 원자력 시설의 경우에는 지반의 전단파 속도가 2440 m/sec 이하이면 지반-구조물 상호작용을 고려하여야 한다[18]. 하지만, 해상풍력발전기의 경우에는 구조물의 고유 주기가 상당히 길므로, 지반-구조물 상호작용을 고려하여야 하는 지반의 최대 전단파 속도가 이 보다는 더 작을 것이다. 이 기준에 대한 연구는 향후 좀 더 다양한 구조 형식을 고려하여 충분한 사례 분석을 통하여 결정하여야 할 것이다.

4.4 기하 비선형의 영향

 기하 비선형의 효과를 살펴보기 위해, 식(15)에서 기하 강성을 포함한 경우와 포함하지 않은 경우 해상풍력발전기의 지진응답이 어떻게 변화하는지 살펴보았다. 지반의 전단파 속도는 200m/sec로 가정하여 유체-구조물-지반 상호작용의 효과를 고려하였고, 선형화된 항력계수 C′_D  는 기하 비선형을 고려한 경우와 고려하지 않은 경우 각각 0.0351과 0.0379로 설정하였다.

 Fig. 13과 14는 각각 tower 최상단의 변위, tower 밑면에서 휨모멘트의 전달함수와 시간이력을 보여주고 있다. 기하 비선형을 고려하면 구조물의 고유 진동수가 고려하지 않은 경우보다 작아지므로, 구조물의 응답은 다소 영향을 받게 된다. 하지만, 그 영향은 그리 심각하지 않음을 Fig. 13과 14에서 확인할 수 있다.

Fig. 13. Displacements at the top of the tower: effects of the geometric nonlinearity

Fig. 14. Base overturning moments of the tower: effects of the geometric nonlinearity

4.5 해수 압축성의 영향

  일반적으로 해상풍력발전기 지진응답 해석 시 해수의 영향은 Morison 방정식을 사용하여 고려한다[7]. 이 때 유체의 부가질량 계수 A C 는 1.0을 사용하는데, 이 값은 무한한 길이를 가지는 원통형 구조물에 가해지는 비압축성 이상 유체의 동수압으로부터 얻어진 값이다[19]. Morison 방정식을 사용하게 되면 식(15)의 운동방정식에서 부가질량의 표현만이 다음과 같이 바뀌게 된다.

 하지만, 이 연구에서는 압축성 이상 유체의 파동방정식으로부터 해수의 경계조건과 방사조건을 만족하는 동수압을 식(2)와 같이 구하여 사용한다. 그러므로, 해수 압축성의 효과가 해상풍력발전기의 지진응답에 어떠한 영향을 미치는지 살펴보고자 한다. 식(2)의 엄밀한 표현을 사용한 경우(case ‘compressible’)와 Morison 방정식을 사용한 경우(case ‘incompressible’)의 지진응답을 각각 구하였다. 지반의 전단파 속도는 200 m/sec로 가정하여 유체-구조물-지반 상호작용의 효과를 고려하였고, 기하 비선형 효과를 고려하였다. 선형화된 항력계수 C′_D는 엄밀한 동수압을 사용한 경우와 Morison 방정식을 사용한 경우 각각 0.0351과 0.0347로 설정하였다.

 Fig. 15와 16은 각각 tower 최상단의 변위, tower 밑면에서 휨모멘트의 전달함수와 시간이력을 보여주고 있다. 이 예제에서는 유체의 압축성을 고려한 경우와 그렇지 않은 경우의 응답 차이가 아무 미미하다. 이 두가지 경우의 차이를 명확하게 조사하기 위해, 유체로 인한 부가질량, 즉 식(13)과 (17)의 부가질량을 계산하여 Fig. 17에 비교하였다. 단, 식(13)의 부가질량은 복소수이므로 그 실수부만을 비교하였다. Fig. 17에서 확인할 수 있듯이 유체 압축성 고려 여부에 따라 유체의 부가질량은 영향을 받는다. 하지만, 그 차이는 전체 질량 M_ij + (i,j = 1,…, N,Δ,φ)와 비교하였을 때 그리 크지 않다.

Fig. 15. Displacements at the top of the tower: effects of the fluid compressibility

Fig. 16. Base overturning moments of the tower: effects of the fluid compressibility

Fig. 17. Added mass from the sea water

 식(13)의 부가질량은 복소수인데, 실수부가 부가질량의 의미를 가진다면 허수부는 유체의 방사 감쇠(radiation damping)와 연관된다. 즉, 해수의 무한 영역으로 발생하는 에너지 방사로 인한 감쇠를 의미한다. Morison 방정식, 즉, 식(17)의 부가질량을 사용하여서는 이러한 방사 감쇠를 고려할수 없다. Fig. 18은 유체의 압축성으로 인한 방사감쇠 Im[-ω](i,j = 1,…, N,Δ,φ)를 항력에 의한 감쇠 (i,j = 1,…, N,Δ,φ) 및 전체 시스템의 감쇠와 비교한 그림이다. 전체 시스템의 감쇠는 구조물의 재료 감쇠 2ξK_i,j / ω(i, j = 1, …, N)와 지반의 감쇠/ω(i, j = Δ,φ )를 포함한다. Fig. 18에서 유체의 방사감쇠는 cut-off frequency C_w / 4H= 7.195Hz 이후에 주요하게 발현되는 것을 확인할 수 있다.

Fig. 18. Radiation damping of the sea water

 Fig. 17과 18로부터 가진 진동수가 유체의 cut-off frequency보다 작은 경우에는 유체의 압축성이 시스템의 응답에 큰 영향을 미치지 않고, 큰 경우에는 방사 감쇠로 인한 추가 감쇠가 발생하여 시스템의 응답이 더 감소할 것임을 예상할 수 있다. Fig. 15와 16의 전달함수에서 확인할 수 있듯이 가진 진동수가 유체의 cut-off frequency보다 큰 영역에서의 시스템 응답이 차지하는 비중은 작다. 그러므로, 주어진 예제에서 유체의 압축성 고려 여부가 전체 시스템의 응답에 미치는 영향은 아주 작게 되는 것이다. 일반적으로 해상풍력발전기 tower는 아주 유연한 구조물이고 해수의 깊이가 그리 깊지않으므로, 계산의 효율성을 위해 간단한 표현의 Morison 방정식을 사용하여도 될 것이다. 하지만, 구조물이 충분한 강성을 가지고 유체의 깊이가 충분히 깊을 경우에는 유체의 압축성은 구조물의 거동에 상당한 영향을 끼칠수 있으므로[20], 이러한 경우에는 유체 압축성의 효과가 어느 정도인지 반드시 살펴보아야 할 것이다.

 한편, Fig. 18에서 해수의 항력에 의한 감쇠는 전체 시스템의 감쇠와 비교하였을 때 무시할 만한 수준임을 확인할 수 있다. 이는, 해상풍력발전기의 지진응답 해석 시 해수의 항력은 무시할 수도 있음을 의미한다.

5. 결 론

 이 연구에서는 구조물 기하 비선형의 효과와 해수 항력의 효과를 고려한 해상풍력발전기의 지진응답 해석기법을 개발하였다. 구조물의 변형과 기초의 회전운동에 의해 발생하는 기하 강성을 유도하여 기하 비선형의 효과를 고려하였고, Morison 방정식의 항력 부분을 선형화하여 해수의 항력의 효과를 고려하였다. 또한, 상부구조의 유연성, 상부구조와 해수와의 동적상호작용, 기초에서 유연한 지반과의 상호작용으로 인한 유체-구조물-상호작용을 엄밀히 고려하였다.

 5 MW급 해상풍력발전기의 지진응답 해석을 통하여 유체-구조물 상호작용, 지반-구조물 상호작용, 지반 강성 등이 해상풍력발전기의 지진응답에 상당한 영향을 미칠 수 있지만, 기하 비선형과 유체 압축성의 영향은 그리 심각하지 않음을 확인하였다.

 이 연구에서는 기본적인 구조 형식과 기초 형식만을 고려하였지만, 이 연구를 기반으로 하여 좀 더 다양하고 복잡한 구조물과 하중에 적용이 가능하도록 발전시킬 수 있을 것이다. 향후에는 좀 더 다양한 구조 형식과 기초형식이 해상풍력발전기의 지진응답에 어떠한 영향을 미치는지 조사하여야할 것이다. 또한, 해상풍력발전 시스템에는 지진하중 외에도 파랑하중, 조류하중, 풍하중 등의 환경하중이 지속적으로 가해지고 있다. 특히, 파랑하중은 지진하중과 달리 장주기 성분이 우세하다. 그러므로 이러한 다양한 하중 조건하에서 해상풍력발전기의 동적 응답이 어떠한 영향을 받는지 살펴보아야 할 것이다.

감사의 글

 본 연구는 국토해양부 산하 한국해양과학기술진흥원의 첨단항만건설기술개발사업 연구비 지원 (과제명: 해상풍력 지지구조 설계기준 및 콘크리트 지지구조물 기술 개발/20120093)을 받아 수행되었습니다.