Advanced
Dynamic Response Analysis of Top-tensioned Riser Under Sheared Current Load
Dynamic Response Analysis of Top-tensioned Riser Under Sheared Current Load
Journal of Ocean Engineering and Technology. 2013. Aug, 27(4): 83-89
Copyright © 2013, Korean Society of Ocean Engineers
  • Received : August 07, 2013
  • Accepted : August 12, 2013
  • Published : August 31, 2013
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
국현 김

Abstract
A numerical scheme based on a mode superposition method is presented for the dynamic response analysis of a top-tensioned riser (TTR) under sheared current loads. The natural frequencies and mode shapes of the TTR have been calculated analytically for a beam with a slowly varying tension and pinned-pinned boundary conditions at the top and bottom ends. The lift coefficients and corresponding amplitudes used to estimate the vortex-induced modal force and damping for each mode were predicted via iterative calculations based on the input and output power balancing concept. Here, the power-in regions were controlled by the normal distribution function, for which the center was coincident with the lock -in location by local vortex-shedding, and the range was defined by the constant standard deviation for the reduced velocity by the local current speed. Finally, dynamic responses such as root-mean-squared displacement and stress were calculated using the mode superposition technique. In order to verify the presented scheme, a numerical calculation was performed for a TTR under an arbitrary linearly sheared current and linearly varying tension. A comparison with the results of the existing software showed that the presented scheme could give reliable and feasible solutions. Case studies were performed to investigate the effects of various current loads and tensions.
Keywords
1. 서 론
상부장력 라이저(Top-tensioned riser)는 해양에서의 석유 및 가스 생산을 위한 해양플랜트(Offshore plant)에 있어서 원유의 이송을 위해 해저시스템과 장력 다리 플랫폼(Tension leg platform) 등의 부유식 플랜트를 연결하는 핵심 구성품 중 하나이다. 상부장력 라이저는 수명기간 중 구조적 안전성 확보를 위해 해류, 파랑등 해양환경에 의한 외부하중을 잘 견디도록 설계되어야 한다. 상부장력 라이저가 설치되는 해역에 해류가 발생하는 경우, 라이저 주위에는 유체박리에 의한 와류 방출(Vortex shedding) 현상에 의해 라이저 표면에 주기적인 하중이 발생하며, 이러한 현상을 와류유기 진동(Vortex-induced vibration)이라 한다. 이 때, 와류방출 주파수가 라이저의 고유주파수(Natural frequency)에 가까워질수록 큰 진폭의 라이저 동적응답(Dynamic response)이 나타나는 Lock-in 현상이 나타나며( Park et al., 1999 ), 상부장력 라이저에 주기적인 응력을 발생시켜 장기적으로는 피로손상의 원인이 된다. 따라서 전통적으로 와류유기 진동은 상부장력 라이저 설계에 있어서 그 피로수명을 예측하기 위한 주요 하중인자로 고려되어 왔다. 특히, 최근 해양자원 개발영역이 천해에서 심해로 이동하면서 상부장력 라이저의 길이는 증가하고 있으며, 그에 따른 와류유기 진동 현상도 보다 광범위하게 나타나고 있다. 이에 와류유기 진동을 고려한 상부장력 라이저의 정도 높은 동적응답 해석기법에 대한 연구 필요성이 더욱 증대되고 있다.
와류유기 진동을 고려한 상부장력 라이저의 동적응답(Dynamic response)에 대한 해석모델로는 전산유체역학 모델(Computational fluid dynamics model)과 경험 모델(Empirical model)이 있다. 전산유체역학 모델로는 라이저를 2차원 스트립들(2 dimensional strips)의 조합으로 간주하고 각각의 스트립에 작용하는 유체동력학적 힘과 감쇠를 구하는 방법이 주로 적용되고 있으며, 보다 실제에 가까운 현상을 시뮬레이션할 수 있는 장점이 있으나 과도한 연산부담으로 인해 실용적으로는 적용에 한계가 있다. 반면, 경험 모델은 대부분 모드중첩법(Mode superposition method)을 기반으로 하고 있으며, 각각의 모드에 작용하는 유체동력학적 힘과 감쇠를 실험을 통해 기 확보된 데이터베이스를 활용하여 산정하여 반영하는 방법이다. 전산유체역학 모델에 비해 연산부담이 매우 적고, 모형 또는 실물시험 등과의 보정을 통해 지속적인 정도 향상이 가능하므로 실무적 관점에서 매우 유용한 방법이라 할 수 있다. 최근 들어 경험 모델을 기반으로 한 다양한 해석기법들이 제안되었으며, 이는 VIVA( Triantafyllou, 2003 ), VICoMo ( Moe et al., 2001 ), VIVANA( Larsen et al., 2005 ), SHEAR7 ( Vandiver and Li, 2005 ) 등의 프로그램으로 구현되어 실무에 활용되고 있다.
한편, 경험 모델을 이용한 와류유기 진동응답에 있어서 파워입력 범위 산정은 양력계수, 감쇠계수 추정 등과 같은 입력 요소들 못지 않게 매우 중요한 요소이다. 파워 입력 범위를 과도하게 넓게 산정할 경우 동적 응답이 실제에 비해 과도하게 추정될 수 있고 반대로 그 범위를 너무 좁게 산정할 경우 과소하게 추정될 수 있다. VIVANA, SHEAR7 등의 프로그램에서는 각각의 모드별 파워 입력 범위를 산정한 후 모드가 중첩되어 가진되지 않도록 범위를 조정하여 해석을 수행한다( Larsen et al., 2005 ; Vandiver and Li, 2005 ). 이 때, 조정된 범위 내에서는 파워가 균일(Uniform)하게 입력되는 것으로 간주하여 해석을 수행하는데, 최근 실해역 시험결과와의 비교연구에서 실제 응답보다 과도하게 추정된다는 보고가 있었으며, 이에 대한 대안으로 모드별 기여도를 고려해 확률적으로 입력파워를 분배하는 소위 시간분배(Time sharing technique) 기법이 제안되고 현재 SHEAR7 4.6 버전에 반영된 바 있다(Vandiver et al., 2011).
본 연구에서는 전단류 하중(Sheared current load)을 받는 상부장력 라이저의 와류유기진동에 의한 동적 응답 해석을 위한 기초연구로써 기존 프로그램에서 채택하고 있는 해석기법을 참고하여 모드중첩법 기반의 수치해석기법을 정립한다. 제안된 수치해석기법은 고유진동해석, 각 모드에 작용하는 양력 및 감쇠 추정, 모드 중첩에 의한 동적 응답 해석 등의 절차를 따른다. 이때, 와류방출에 의한 파워 입력 범위 산정을 위해 정규분포함수(Normal distribution function)를 적용한다. 수치해석기법의 타당성을 검토하기 위해 가상의 상부장력 라이저 모델에 대한 수치해석을 수행하고 그 결과를 기존 수치해석 프로그램에 의한 결과와 비교하여 적절한 해석인자를 도출한다. 이를 이용해 전단류와 장력의 변화에 따른 동적 응답 해석을 추가적으로 수행하고 그 결과 고찰을 통해 실무 적용 가능성을 검토한다.
2. 해석이론
- 2.1 지배방정식
길이가 지름에 비해 매우 큰 경우(약 200배 이상), 상부장력 라이저는 유연한 케이블로 간주할 수 있으며, 와류 방출에 의해 외력이 작용하는 경우 전체 동적 응답에 대한 지배방정식은 식 (1)과 같다( Vandiver and Li, 2005 ).
PPT Slide
Lager Image
여기서, m ( s )는 접수에 의한 부가질량을 포함하는 단위 길이당 등가질량, R ( s ) 는 유체동력학적 감쇠(Hydrodynamic damping)와 구조감쇠(Structural damping)를 모두 포함하는 단위 길이당 감쇠, T ( s ) 는 유효장력(Effective tension), P ( s , t ) 는 외력을 의미하며, 이 때 s 는 위치좌표, t 는 시간을 나타낸다.
전체 시스템의 동적 응답을 식 (2)와 같이 고유 모드의 중첩으로 가정하고 식 (1)에 대입하여 정리하면, 식 (3)과 같은 모드좌표계에서의 운동방정식을 얻을 수 있다.
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서, Y i ( s ) 는 i 번째 고유모드, q i ( t ) 는 모드좌표계를 의미하며, Mi, Ri, Ki , Pi ( t ) 는 식 (4)에서 식 (7)까지 정의된 i 번째 모달 질량(Modal mass), 모달 감쇠(Modal damping), 모달 강성(Modal stiffness), 모달 힘(Modal force)이다.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
- 2.2 입력-출력파워 평형
해류의 와류 방출에 의해 양력이 외력으로 작용하는 경우, 외력은 식 (8)으로부터 계산된다. 또한, 국부에 작용하는 힘과 i 번째 모달 속도(Modal velocity)가 동위상(In-phase)인 것으로 간주하면 라이저에 입력되는 파워
PPT Slide
Lager Image
는 식 (9)가 되며, 한 주기동안의 시간 평균 입력파워 를 구하면 식 (10)이 된다.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서, ρf 는 유체 밀도, L 은 라이저 길이, D 는 라이저 지름, V ( s ) 는 해류 속도 분포, CL(s;ωi) i 번째 고유모드에 대한 양력계수(Lift coefficient), Ai i 번째 모드변위의 진폭(Mode displacement amplitude)을 의미한다. 또한, 𝛷 ( s ) 는 양력에 의한 가진 범위(excitation region)와 기여도를 정의하기 위한 함수이며, 본 연구에서는 평균중심에서 그 크기가 1이고 표준편차(standard deviation: STD)를 갖는 정규분포함수를 이용하였다.
이와 마찬가지 방법으로 i 번째 모드에 대한 시간평균 출력파워 <
PPT Slide
Lager Image
> 를 유도하면 식 (11)이 된다.
PPT Slide
Lager Image
여기서, 𝛷′( s ) (=1−𝛷( s ) )는 감쇠범위(Damping region)를 정의하는 함수이며, Rs ( s ) 는 구조감쇠계수이다. 또한, Rh ( s ) 는 유체동력학적 감쇠계수로써, i 번째 모드변위의 진폭 Ai 에 의해 결정된다.
각각의 모드에 대해 입력 파워와 출력 파워의 평형(Input- output power balancing,
PPT Slide
Lager Image
)을 고려하면 식 (12)를 얻을 수 있다.
PPT Slide
Lager Image
식 (12)에서 좌변과 우변의 분모 항은 i 번째 모드변위의 진폭 Ai 에 의존적이며, 우변의 분자 항은 양력계수 CL(s;ωi) 에 의존적이다. 따라서 이 식을 문제에 적합한 Ai CL(s;ωi) 를 찾는 반복 식으로 활용 가능하다. 이 때 양력계수와 유체동력학적 감쇠계수는 실험 결과로부터 도출된 경험 모델을 적용한다.
- 2.3 양력계수 모델
양력계수는 무차원 응답진폭( A/D , Non-dimensional response amplitude)에 대한 함수로 표현된다. Fig. 1 은 본 연구에서 적용한 포물선 기반의 양력계수 모델을 나타낸 것으로써( Larsen et al. 2005 ), 양력계수가 최대가 되는 점( A/D max , CL max )을 꼭지점 으로 하고 좌측으로는 (0, CL ,0 )인 점을 지나는 포물선과 우측으로 는 ( A/D 0 , 0)인 점을 지나는 두 개의 포물선의 조합으로 구성되며, 양력계수가 무한정 작아지는 것을 방지하기 위한 최소값 CL ,min 이 추가로 정의된다.
PPT Slide
Lager Image
Lift coefficient model
- 2.4 유체동력학적 감쇠모델
유체동력학적 감쇠계수는 경험 모델을 적용하되 감소속도(Reduced velocity, UR ( s ) = 2π V ( s )/( i ))가 5인 해류 속도를 기준으로 저속구간과 고속구간으로 나누어 적용한다. 식 (13)과 식 (14)는 저속구간과 고속구간인 경우에 대한 유체동력학적 감쇠 산정 식을 각각 나타낸다( Vikestad et al. 2000 ).
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서, Re ω (=ω D 2 / ν )는 진동 레이놀즈수(Vibration Reynolds number), ν 는 유체점성(Fluid viscosity)이며, Csw , Crl , Crh 은 경험계수이다.
- 2.4 모드중첩법에 의한 동적 응답 산정
식 (3)에서 감쇠항을 무시하고 조화가진을 가정한 고유치 문제의 해를 구하면 고유주파수와 고유모드를 구할 수 있으며, 식 (12)의 반복계산에 따라 구한 양력계수와 감쇠계수를 활용하면 모드중첩법에 의한 동적 응답을 산정할 수 있다.
상부장력 라이저에 대한 RMS 변위(Root-mean-squared displacement,
PPT Slide
Lager Image
( s )) 와 RMS 응력(Root-mean-squared stress,
PPT Slide
Lager Image
( s ))는 식 (15) 및 식 (16)을 이용해 산정할 수 있다( Vandiver and Li, 2005 ).
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서, E 는 탄성계수(Young's modulus),
PPT Slide
Lager Image
j 번째 모드의 곡률(Curvature), Ds 는 강력 지름(Strength diameter)이다. 또한,
PPT Slide
Lager Image
Hji(s) 는 식 (17)과 식 (18)에 정의한 바와 같다.
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서, sgn [•]는 내용(Kernel)이 음인 경우 −1, 0인 경우 0, 양인 경우 +1인 부호함수(Sign function)이다. 또한, 𝜁j j 번째 모드에 대한 모달 감쇠비(Modal damping ratio, 𝜁j = 𝜁j,s + Rj,h /(2 ωjMj )), 𝜁j,s 는 구조감쇠비(Structural damping ratio)이다.
3. 수치해석
정립된 수치해석기법의 타당성을 검토하기 위해 가상의 상부장력 라이저 모델에 대한 수치해석을 수행하고 그 결과를 기존 프로그램인 SHEAR7(Vandiver et al., 2011)에 의한 결과와 비교하였다. Table 1 은 수치해석 모델을 요약한 것이다. 해석모델은 양단이 단순지지(Pinned-pinned)되어 있고 강(Steel)으로 되어 있는 라이저이며, 그 길이는 1,000m, 지름과 강력 지름은 공히 0.6m로 설정하였다. 해류는 하단(Bottom) 0.2 m/s 에서 상단(Top) 1.2 m/s 까지 선형적으로 점차 증가하는 전단류를 가정하였으며, 유효 장력 또한 하단기준 5.00e6 N 에서 상단 5.05 N 까지 선형적으로 점차 증가하는 것으로 가정하였다. 참고로, 양력계수 및 감쇠계수 추정하기 위해 Table 2 에 정리된 계수들을 사용하였다.
Virtual top-tensioned riser model
PPT Slide
Lager Image
Virtual top-tensioned riser model
Lift coefficient and damping model for numerical calculation
PPT Slide
Lager Image
Lift coefficient and damping model for numerical calculation
한편 해석모델에 대한 n 번째 고유주파수 ωn 는 식 (19), 고유모드 형상 Yn(s) 와 곡률 Yn″(s) 은 식 (20) 및 식 (21)을 이용하여 산정하였다( Vandiver and Li, 2005 ).
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
Table 3 은 10차까지의 고유주파수 산정결과를 SHEAR7의 결과와 함께 정리한 것으로써, SHEAR7의 결과가 다소 높게 나타나며 고차모드로 갈수록 차이가 커지나 전반적으로 2.3% 차이 이내의 타당한 결과를 주고 있다. 이러한 차이는 본 연구에서는 라이저를 굽힘 강성이 없는 유연한 케이블로 간주한 반면, SHEAR7에서의 결과는 굽힘강성을 갖는 보(Beam)로 가정하였기 때문인 것으로 판단되며, 고유진동해석 시 정도 높은 해석기법을 활용할 경우 충분히 개선될 수 있을 것으로 사료된다. Fig. 2 는 차수별 고유모드 형상과 곡률을 도시한 것이다.
Natural frequency calculation results
PPT Slide
Lager Image
Natural frequency calculation results
PPT Slide
Lager Image
Shape of (a) natural mode and (b) curvature
라이저의 모드별 가진범위(Power-in region)에 따른 해석결과의 변화를 살펴보기 위해 STD가 0.2, 0.3, 0.4, 0.5, 0.6인 정규분포 함수로 가진범위를 정의한 경우에 대한 수치해석을 수행하였다. Fig. 3 은 STD가 0.2와 0.4인 경우의 가진범위를 대표적으로 도시한 것이며, STD가 커질수록 가진범위가 넓어짐을 확인할 수 있다. Fig. 4 는 RMS 변위와 RMS 응력을 STD 별로 도시한 것이며, 본 연구에서 정립한 해석기법의 타당성과 적절한 STD값을 도출을 위해 SHEAR7에 의한 해석결과와 함께 나타내었다. 이로 부터 STD가 클수록 동적 응답이 과도하게 산정되고 있으며, 특히 일부 프로그램에서 적용하고 있는 균일 분포 함수(Uniform distribution function)를 가정하는 경우 상당한 오차가 발생함을 알 수 있다. 전반적으로 하단(Bottom)부근에서는 STD가 0.6인 경우, 상부(Top)부근에서는 STD가 0.4인 경우, SHEAR7의 결과와의 부합성이 높게 나타나며, 평균적으로는 STD가 0.4인 경우의 부합성이 크게 나타나고 있다. 이는 본 연구에서 제안한 정규분포함수를 이용해 파워입력 범위를 정의하되 적절한 STD을 선정할 경우 신뢰성 높은 해석이 가능함을 의미한다. 참고로, 본 해석에서는 전단류에 의해 가진 가능성이 있는 주파수를 고려할 때 최고차 모드는 8번으로 파악되었으나 실제 계산에 있어서는 고차 모드영향을 고려해 1.5배에 해당하는 13번까지의 모드를 포함하여 해석을 수행하였다.
PPT Slide
Lager Image
Power-in region by mode: (a) STD=0.2 and (b) STD=0.4
PPT Slide
Lager Image
Root-mean-squared (a) displacement and (b) stress distribution
한편, 본 연구에서 정립한 해석기법을 이용해 유효장력(T-1, T-2), 전단류(C-1, C-2) 변화에 따른 CASE 별 수치해석을 수행하고 그 결과를 살펴보았다. 이 때, 입력파워 범위 정의를 위한 정규분포함수의 STD는 0.4를 적용하였다. Table 4 는 CASE별 해석에 적용된 유효 장력과 전단류를 정리한 것이다.
Tension and current variation for case study
PPT Slide
Lager Image
Tension and current variation for case study
Fig. 5 는 전단류 변화(C-1, C-2)에 따른 라이저 동적응답해석 결과를 Fig. 4 의 STD=0.4(Reference)에 대한 것과 비교하여 나타낸 것이다. 변위( Fig. 5(a) )의 경우, 전단류가 강해질수록 커지고 약할수록 작아지는 특성이 나타나며, 최대값의 위치도 함께 변화하고 있다. 응력( Fig. 5(b) )의 경우, 유사한 양상을 보이지만 응력의 변화폭이 매우 크게 나타나고 있다. 전단류가 강해질수록 고차모드에 의한 기여도가 강하게 나타나고 있음을 확인할 수 있으며, 이는 변위 변동폭이 유사한 수준이더라도 고차모드가 지배적인 경우 응력의 변동폭이 크게 나타날 수 있으므로 응력관점에서 저차모드 보다는 고차모드에 대한 제어가 우선되어야 함을 의미한다.
PPT Slide
Lager Image
Root-mean-squared (a) displacement and (b) stress distribution w.r.t current speed variation (C-1, C-2)
Fig. 6 은 유효장력 변화(T-1, T-2)에 따른 라이저 고유주파수를 기존결과( Fig. 4 의 STD=0.4(Reference))와 함께 막대그래프로 나타낸 것이다. 기존 결과에 비해 큰 고유주파수를 보이며, 이는 유효 장력의 증가에 따른 영향을 판단된다. Fig. 7 은 유효장력 변화에 따른 동적응답해석결과를 Fig. 4 의 STD=0.4(Reference)에 대한 것과 비교하여 나타낸 것이다. 기존 결과에 비해 유효장력이 증가함에 따라 변위와 응력이 감소하고 있음을 알 수 있다. 이는 유효장력의 증가가 고유주파수를 증가시키므로 동일한 전단류 하중 하에서 동적응답에 기여하는 고유모드의 수가 감소했기 때문인 것으로 판단된다. 또한, 유효장력 증가에 따라 변위보다는 응력의 변화폭이 크게 나타나고 있다. 이는 유효장력이 증가하면서 응력에 크게 기여하는 고차모드의 영향이 사라지기 때문인 것으로 판단된다.
PPT Slide
Lager Image
Natural frequency variation w.r.t. tension (T-1, T-2)
PPT Slide
Lager Image
Root-mean-squared (a) displacement and (b) stress distribution w.r.t tension variation (T-1, T-2)
5. 결 론
본 연구에서는 전단류 하중(Sheared current load)을 받는 상부장력 라이저의 와류유기진동에 의한 동적 응답 해석을 위한 기초연구로써 기존의 프로그램에서 해석기법으로 채택하고 있는 모드중첩법 기반의 수치해석기법을 정립하였다. 이 때, 와류 방출에 의한 파워 입력 범위 산정을 위해 정규분포함수를 도입하고 가상의 상부장력 라이저에 대한 수치해석 결과와 기존 프로그램에 의한 결과를 비교하여 그 타당성을 검토하였다. 또한, 정립된 해석기법을 활용해 전단류와 유효 장력의 변화에 따른 CASE 연구를 수행하고 이들이 라이저의 동적 거동 특성에 미치는 영향을 검토하였으며, 다음과 같은 현상을 관찰할 수 있었다.
(1) 전단류가 강해질수록 라이저의 동적응답은 커지며, 특히 고차모드의 영향으로 변위보다는 응력의 변화폭이 크게 나타났다. 이는 변위 변화폭이 크지 않더라도 저차모드보다는 고차모드가 응력에 크게 기여함을 의미한다.
(2) 유효 장력이 증가하면, 라이저 고유주파수는 증가하고 응력은 감소한다. 이는 라이저의 강성이 증가하면서 고차모드의 영향이 감소하기 때문인 것으로 판단된다.
이상의 결과로부터 제안된 해석기법이 실무에 효과적으로 활용될 수 있을 것으로 기대되며, 향후, 보다 높은 정도 확보를 위해 모형 또는 실 규모 실험자료와의 비교를 통한 타당성 연구가 추가적으로 필요할 것으로 사료된다.
References
Larsen C.M. , Vikestad K. , Yttervik R. , Passano E. , Baarholm G.S. 2005 VIVANA - Theory Manual Version 3.4 Marinetek
Moe G. , Arntsen Ø. , Hoen C. 2001 VIV Analysis of Risers by Complex Modes Proceedings of the 11th International Offshore and Polar Engineering Conference 3 426 - 430
Park H.I. , Cho H.J. , Jung D.H. 1999 Response Analysis of Deep Ocean Risers to Vortex-Shedding by Numerical Analysis Journal of Ocean Engineering and Technology 13 (3B) 65 - 72
Triantafyllou M.S. 2003 VIVA Extended User's Manual Massachusetts Institute of Technology
Vandiver J.K. , Li L. 2005 SHEAR7 V4.4 Program Theoretical Manual Massachusetts Institute of Technology
Vikestad K. , Larsen C.M. , Vandiver J.K. 2000 Damping of Vortex-Induced Vibrations, Paper No. OTC 11998 Offshore Technology Conference Houston, Texas