KR101512621B1

KR101512621B1 - 응답면 기법을 이용한 해상풍력용 모노파일 신뢰성 해석 방법

Info

Publication number: KR101512621B1
Application number: KR1020130110943A
Authority: KR
Inventors: 윤길림; 김선빈; 이진학
Original assignee: 한국해양과학기술원
Priority date: 2013-09-16
Filing date: 2013-09-16
Publication date: 2015-04-16
Also published as: KR20150032608A

Abstract

응답면 기법을 이용한 해상풍력용 모노파일의 신뢰성 해석 방법을 개시한다. 상기 응답면 기법을 이용한 해상풍력용 모노파일의 신뢰성 해석 방법은 해상풍력터빈 구조물의 모노파일 하부구조가 모래지반에 설치된 경우에 확률론적 해석을 통하여 상기 모노파일의 거동 확률분석 및 이에 따른 민감도를 분석하는 방법으로, 상기 모래지반에 설치된 상기 모노파일의 안전성을 확률적으로 분석하기 위한 고정 변수 및 한계상태별 확률변수를 설정하는 단계(S10); 및 상기 모노파일의 단계 확률변수 범위에서 응답면 기법을 적용하여 한계상태함수를 확률변수에 대한 양함수의 형태로 근사화 하는 단계(S20); 및 단계 근사화된 한계상태함수를 일계신뢰도법(FORM)에 적용하여 파괴확률 (신뢰도지수)를 산출하는 단계(S30)를 포함한다.

Description

응답면 기법을 이용한 해상풍력용 모노파일 신뢰성 해석 방법{Reliability Analysis Method of Monopile for Offshore Wind Turbines Using Response Surface Method}

본 발명은 해상풍력용 모노파일의 신뢰성 해석법을 발명하는 것으로, 보다 상세하게는 해상풍력 기초공법중에서 수심 30m 이하에서 비교적 많이 활용되고 있는 모노파일의 거동해석을 위하여 국외에서 개발된 풍력터빈의 실제 모델을 선정하고 그에 따른 자중 및 운용환경을 고려하여 적용이 가능한 임의의 지반조건하에서 신뢰성 해석을 통하여 모노파일의 안전에 대한 신뢰성을 확률론적으로 분석하기 위한 응답면 기법을 이용한 해상풍력용 모노파일의 신뢰성 해석 방법에 관한 것이다.

탄소배출권 제한과 신재생 에너지 개발이 범 지구적인 화두가 됨에 따라 해상풍력발전단지 건설에 대한 관심과 연구가 확대되고 있다.

특히, 우리나라의 경우 국토면적이 좁은 반면 바람이 풍유하고 삼면이 바다인 관계로 해양공간 활용에 매우 유리한 입지를 가지고 있어 해상풍력단지건설에 대한 관심이 크다.

해상환경은 비교적 고밀도의 바람에너지가 안정적으로 공급되고 소음공해와 전파장애 등 민원문제도 없어 발전시설의 건설과 운영에 우려되는 제한이 매우 적기도 하다. 해상에 풍력발전시설을 건설할 때 나셀, 블레이드 및 타워 등 상부구조물의 자중과 풍하중 및 파랑하중, 터빈가동에 따른 모멘트 등 발생되는 외력에 저항하기 위해서 해저지반에 설치되는 다양한 기초공법이 현재까지 많이 개발되어 왔다. 이러한 기초공법은 주로 경제성과 안전성을 극대화하는 측면에서 수심, 해저지반조건, 타워의 높이 등과 연관되어져 선택되고 설계에 활용되었다.

해상풍력발전 구조물을 시공하여 운용하는 경우에 육상에 설치되는 구조물과 그 환경조건 매우 다르므로 기초공법도 무리말뚝같은 파일기초 이외에 중력식, 자켓식, 모노파일식, 석션파일 및 트라이포드식 등 다양한 형식으로 개발되었다.

그리고 최근에는 7MW급 고효율, 대용량의 풍력터빈 개발이 계속적으로 이루어짐에 따라서 적용되는 기초공법도 지지능력이 보다 큰 형식으로 대형화 추세로 변화하고 있다.

특히, 해상풍력 하부기초는 안전성 및 경제성 측면에서 구조물 설계기술의 수준이 높고 경제적으로 최적화하기 위한 방향으로 변화하고 있는 추세이다. 해상풍력발전 건설의 많은 경험을 가지고 있는 덴마크, 독일 및 영국 등 유럽, 석유시추 등 해상구조물 구축에 오랜 노하우를 가진 미국을 중심으로 신뢰성 기반 한계상태 설계법에 근간한 설계기준이 개발된 것도 그 일환이라 하겠다.

한편, 다양한 기초형식 중 모노파일 기초는 다른 기초형식과 비교할 때 형태가 단순하고 시공성이 탁월하여 해상시공에 유리하며 무엇보다도 경제성이 뛰어난 장점을 갖고 있어 해상풍력 기초구조물로 많이 사용되고 있다.

모노파일의 해석에 있어 해저지반의 강성을 고려하지 않는 지반고정식보다 해저지반의 강성을 적절히 고려할 수 있는 가상고정식과 등분포 스프링식이 해저지반의 특성을 잘 고려하기 때문에 보다 안전한 설계를 할 수 있다. 또한, 풍력터빈 블레이드의 회전운동과 고유진동에 의해 해저지반에 관입된 기초구조물이 반복하중이 작용하게 됨에 따라 하부기초인 말뚝의 지지력에 미치는 지반반력계수가 감소하는 경향이 있다.

일반적으로 해상풍력발전기 기초는 바람과 파랑 및 조류 등의 수평하중이 크게 반복적으로 작용하므로 말뚝 두부에서의 수평변위나 회전도를 허용기준 이하로 설계해야 하며, 따라서 해상풍력발전기 기초거동에 대한 신뢰성 해석시 주요 파괴모드 및 한계상태함수는 말뚝두부의 수평변위와 회전각에 대한 함수로 정의된다.

현재까지 말뚝의 수평거동은 수평하중에 대하여 말뚝의 변형과 말뚝주변의 흙의 저항력을 표현하는 p-y curve에 의해 비선형해석법이 가장 합리적이라 알려져 있으며, 이러한 비선형 p-y 거동해석은 보통 유한차분법(FDM)에 의해 해석하기도 한다. 따라서, 하지만 FDM 해석으로 얻을 수 있는 종속변수인 말뚝의 수평변위나 회전각 등을 포함하는 한계상태함수는 음함수(implicit function)의 형태이므로 신뢰성 해석에 어려움이 따른다. 응답면 기법(response surface method; RSM)은 이러한 음함수 형태의 한계상태함수를 양함수(explicit function)의 형태로 근사화할 수 있는 기법이다.

이에 본 발명에서는 해상풍력 기초공법 가운데 30m 이하의 수심에서 비교적 많이 활용되고 있는 모노파일에 대해 한계상태 기반의 확률론적 해석을 용이하도록 하기 위해 응답면 기법을 이용한 해상풍력용 모노파일의 신뢰성 해석 방법을 제공하고자 한다.

1. API (2005). Recommended Practice for Planning, Design and Constructing Fixed Offshore Platforms-Working Stress Design. American Petroleum Institute Publishing Service, Washington D.C., pp. 1-263. 2. Bucher, C.G. and Bourgund, U. (1987). Efficient Use of Response Surface Method, Report No. 9-87, Institute fur Mechanik, Universitat Innsbruck, Innsbruck, Austria. 3. Bush, E. and Manuel, L. (2009). Foundation Models for Offshore Wind Turbines, 47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition, Orlando, Florida, pp. 1~7. 4. Dewaikar, D. and Patil, P. (2006). A New Hyperbolic p-y Curve Model for Laterally Loaded Piles in Soft Clay. Foundation Analysis and Design, pp. 152-158. 5. DNV (2007). Design of Offshore Wind Turbine Structures. Offshore Standard DNV-OS-J101, Det Norske Veritas, Hik, Norway, p. 153. 6. Emmanuel, F. and Guillaume, C. (2008). Reservoir flow uncertainty assessment using response surface constrained by secondary information. Journal of Petroleum Science and Engineering, Vol. 60, pp. 170-182. 7. Faravelli, L. (1989). Response Surface Approach for Reliability Analysis. Journal of Engineering Mechanics, ASCE , Vol. 115, No. 12, pp. 2763-2781. 8. Jonkman, J., Butterfield, S., Larsen, T. J., Passon, P., Camp, T., Nichols, J., Azcona, J. and Martinez, A. (2007). Offshore Code Comparison Collaboration within IEA Wind Annex XXIII: PhaseII Results Regarding Monopile Foundation Modeling. 2007 European Offshore Wind Conference & Exhibition, Berlin, Germany, pp. 1-12. 9. Khuri, A. I. and Cornell, J. A. (1996). Response Surfaces Designs and Analyses, Second Edition, Marcel Dekker. 10. KIOST (2013). HSRBD ver.3.1 Harbour Structure Reliability Based Design (in Korean). 11. Krolis, V. (2007). Foundation Design of Monopile Support Structures. 2007 European Offshore Wind Energy Conference, Nice, France, pp. 1~46. 12. Kuo, Y. S., Achmus, M. and Kao, C. S. (2008). Practical design considerations of monopile foundations with respect to scour. Global WindPower, Beijing, pp. 29-31. 13. Nicholas, Z. (2010). Importance Sampling (Lecture Note). Cornell University. pp. 13-14. 14. Yoon, G., Yoon, Y., Kim, H. and Kim, B. (2010). Partial Safety Factors for Geotechnical Bearing Capacity of Port Structures. Journal of Korean Society of Coastal and Ocean Engineers, Vol. 22, No. 3, pp. 156-162 (in Korean).

본 발명이 해결하고자 하는 과제는 해상풍력 기초공법 가운데 30m이하의 수심에서 많이 활용되고 있는 모노파일의 신뢰성 해석을 위하여 기시공된 풍력터빈의 실제 모델을 선정하고 그에 따른 자중 및 운용환경을 고려하여 적용이 가능한 임의의 지반조건 하에서 신뢰성 해석을 통하여 모노파일의 거동을 확률론적으로 분석가능하도록 응답면 기법을 이용한 해상풍력용 모노파일의 신뢰성 해석 방법을 제공하는 것을 목적으로 한다.

상기 과제를 해결하기 위한 본 발명의 실시 예에 따른 응답면 기법을 이용한 해상풍력용 모노파일의 신뢰성 해석 방법은 해상풍력터빈 하부구조 모노파일의 거동에 관한 확률론적 해석에 응답면 기법(response surface method; RSM)을 연계하여 적용한 해석방법으로, 해상에 설치된 상기 모노파일의 안전성을 확률적으로 분석하기 위한 고정 변수 및 한계상태별 확률변수를 설정하는 단계(S10); 및 상기 모노파일의 확률변수 범위에서 응답면 기법을 적용하여 한계상태함수를 확률변수에 대한 양함수의 형태로 근사화 하는 단계(S20); 및 근사화된 한계상태함수를 일계신뢰도법(FORM)에 적용하여 파괴확률(신뢰도지수)를 산출하는 단계(S30)를 포함하는 것을 특징으로 한다.

상기 고정 변수는 상기 모노파일의 연직하중, 길이, 직경, 수심 작용력, 수평 작용력, 지층조건을 포함하는 것을 특징으로 한다.

상기 확률 변수는 상기 모노파일의 수평력, 모멘트, 탄성계수, 해저 기초지반의 점착력과 내부마찰각을 포함하는 것을 특징으로 한다.

상기 S20에서 근사화된 한계상태함수는, 2차 다항식 형태의 수학식 1을 통해 도출되며, [수학식 1]

여기서,

와

은 확률변수와 그 개수이고,

,

및

는 응답면의 형태를 나타내기 위한 계수이며,

,

및

는

개의 미지의 계수로서, 이들 계수 값은 실제 한계상태함수

로부터 일련의 표본점을 이용하여 특이값 분해(singular value decomposition)를 통해 결정되는 것을 특징으로 한다.

본 발명의 실시 예에 따른 응답면 기법을 이용한 해상풍력용 모노파일의 신뢰성 해석 방법을 통해 수평변위 및 회전각에 대해 허용치를 증가시키는 경우에 모노파일 기초의 파괴확률은 대수적으로 감소하였다. 이는 말뚝기초 설계규정인 파일직경(D=6m)의 1%를 허용수평변위(60mm)로 적용하는 경우에 파일의 파괴확률은 1.5%이고 허용수평변위를 최소기준인 15mm로 하는 경우는 파괴확률이 60%로서 그 차이가 크므로 해상풍력단지에 적용하는 말뚝의 수평변위 및 회전각에 대한 허용기준의 정립은 현장여건을 반영하여 정하여야 하며 시급한 과제인 것으로 평가된다.

또한, 모노파일의 직경이 커지면 수평변위 및 회전각에 대한 파괴확률은 감소하나 직경이 7.0m 이상의 경우에는 파괴확률이 점차 수렴하여 그 저감효과는 미미하였다. 이는 모래지반에서 모노파일이 일정한 직경이상 확보되면 말뚝두부의 응력이 점차 하부로 전이되면서 나타나는 현상으로 판단된다. 아울러 파일직경 변화에 대한 영향정도는 말뚝두부의 수평변위에 비하여 회전각이 보다 큰 것으로 나타났으며 이를 최적설계에 활용할 수 있다.

또한, 해저면으로부터 기반암층의 심도가 깊어 파일이 모래지반에 근입되는 마찰말뚝인 경우 파일의 길이가 증가함에 따라 수평변위에 대한 파괴확률은 감소하나 그 감소폭은 작아지며, 파일길이가 일정 근입깊이 이상이 되는 조건, 즉 60m 이상이 되는 경우 수평변위의 감소효과는 작게 평가되었다.

또한, 이계신뢰성해석법인 FORM을 이용한 민감도 분석결과, 변동성이 큰 사질토 지반에서는 흙의 내부마찰각이 파일의 거동에 가장 큰 영향을 미치는 것으로 나타났다. 따라서, 불확실성이 클 것으로 예측되는 기초지반에서는 정밀한 지반조사뿐 아니라 토질의 변동성 분석이 경제성과 안정성을 최적화하는 측면에서 매우 중요한 것으로 판단된다.

도 1은 본 발명의 실시 예에 따른 응답면 기법을 이용한 해상풍력용 모노파일의 신뢰성 해석 방법을 설명하기 위한 플로우 챠트이다.
도 2는 하중(X1) 과 저항(X2)에 대한 구조물이 파괴되는 개념을 설명한 그림으로 하중이 저항보다 큰 경우에는 파괴영역이 되고 작은 경우에서는 안전영역이 되며, 같은 경우에는 한계상태로 표현된다. 또한 한계상태는 선형식이나 비선형 곡선이 된다. 도 2에서 비선형곡선에 이르는 하중과 저항의 변수로 이루어 지는 응답면이 한계상태면에서 이루어 지는 경우를 나타내고 있다
도 3은 MCS와 IS의 정밀도를 비교한 그래프이다.
도 4는 본 발명에 적용된 해상풍력용 모노파일의 레이아웃을 나타낸 예시도이다.
도 5는 허용수평변위의 변화에 따른 파괴확률의 변화를 나타낸 그래프이다.
도 6은 허용회전각의 변화에 따른 파괴확률의 변화를 나타낸 그래프이다.
도 7은 모노파일의 직경 변화에 따른 파괴확률의 변화를 나타내 그래프이다.
도 8은 모노파일의 길이 변화에 따른 파괴확률의 변화를 나타낸 그래프이다.
도 9는 내부마찰각의 불확실성에 따른 신뢰도지수의 변화를 나타낸 그래프이다.
도 10은 확률변수(내부마찰각, 외력, 모노파일의 탄성계수)의 민감도를 나타낸 그래프이다.
도 11은 도 10에 도시된 확률변수의 변동계수 변화에 따른 민감도지수의 변화를 나타낸 그래프이다.

이하, 본 발명의 바람직한 실시 예의 상세한 설명은 첨부된 도면들을 참조하여 설명할 것이다. 하기에서 본 발명을 설명함에 있어서, 관련된 공지 기능 또는 구성에 대한 구체적인 설명이 본 발명의 요지를 불필요하게 흐릴 수 있다고 판단되는 경우에는 그 상세한 설명을 생략할 것이다.

본 발명의 개념에 따른 실시 예는 다양한 변경을 가할 수 있고 여러 가지 형태를 가질 수 있으므로 특정 실시 예들을 도면에 예시하고 본 명세서 또는 출원에 상세하게 설명하고자 한다. 그러나, 이는 본 발명의 개념에 따른 실시 예를 특정한 개시 형태에 대해 한정하려는 것이 아니며, 본 발명의 사상 및 기술 범위에 포함되는 모든 변경, 균등물 내지 대체물을 포함하는 것으로 이해되어야 한다.

어떤 구성요소가 다른 구성요소에 "연결되어" 있다거나 "접속되어" 있다고 언급된 때에는, 그 다른 구성요소에 직접적으로 연결되어 있거나 또는 접속되어 있을 수도 있지만, 중간에 다른 구성요소가 존재할 수도 있다고 이해되어야 할 것이다. 반면에, 어떤 구성요소가 다른 구성요소에 "직접 연결되어" 있다거나 "직접 접속되어" 있다고 언급된 때에는, 중간에 다른 구성요소가 존재하지 않는 것으로 이해되어야 할 것이다. 구성요소들 간의 관계를 설명하는 다른 표현들, 즉 "~사이에"와 "바로 ~사이에" 또는 "~에 이웃하는"과 "~에 직접 이웃하는" 등도 마찬가지로 해석되어야 한다.

본 명세서에서 사용한 용어는 단지 특정한 실시예를 설명하기 위해 사용된 것으로, 본 발명을 한정하려는 의도가 아니다. 단수의 표현은 문맥상 명백하게 다르게 뜻하지 않는 한, 복수의 표현을 포함한다. 본 명세서에서, "포함하다" 또는 "가지다" 등의 용어는 설시된 특징, 숫자, 단계, 동작, 구성요소, 부분품 또는 이들을 조합한 것이 존재함을 지정하려는 것이지, 하나 또는 그 이상의 다른 특징들이나 숫자, 단계, 동작, 구성요소, 부분품 또는 이들을 조합한 것들의 존재 또는 부가 가능성을 미리 배제하지 않는 것으로 이해되어야 한다.

이하, 도면을 참조하여, 본 발명의 실시 예에 따른 응답면 기법을 이용한 해상풍력용 모노파일의 신뢰성 해석 방법을 보다 상세하게 설명하도록 한다.

먼저, 본 발명을 설명함에 앞서, 신뢰성 해석기법의 종류로는 설계자의 편의성을 고려하여 가장 단순화된 방식을 취하는 부분계수법(레벨 1 수준)에서부터 한계상태함수의 평균치 또는 파괴점 부근에서 근사해를 찾는 근사화 방법(레벨 2 수준), 그리고 함수를 근사화하지 않고 직접 적분하거나 수많은 난수를 생성시켜 파괴가 일어나는 경우의 수를 찾는 추출법(레벨 3 수준)에 이르기 까지 매우 다양하다.

이에 본 발명에서는 앞에서 열거한 방법 중 근사화 방법과 추출법을 이용하여 모노파일의 파괴확률을 산정하였다.

또한, 현재까지 해상풍력용 모노파일의 수평거동은 p-y 곡선에 의해 해석함이 가장 합리적이라 알려져 있으며, 이러한 비선형 p-y 거동해석은 유한요소법(FEM)에 의해 이루어진다.

그러나, 유한요소법(FEM)에서 얻을 수 있는 종속변수인 변위나 회전각 등 성능변수의 한계상태함수는 음함수(implicit function)의 형태이므로 직접 한계상태함수로 신뢰성 해석을 수행하기는 복잡하고 많은 시간이 요구된다.

이에, 본 발명에서는 응답면 기법(response surface method; RSM)을 이용하여 음함수 형태의 한계상태함수를 양함수(explicit function)의 형태로 근사화시킨 후, 이로 얻어진 양함수를 통해 모노파일의 파괴확률를 산출 및 산출된 파괴확률의 신뢰성을 검증하고자 한다.

도 1은 본 발명의 실시 예에 따른 응답면 기법을 이용한 해상풍력용 모노파일의 신뢰성 해석 방법을 설명하기 위한 플로우 챠트이다.

도 2는 하중(X1) 과 저항(X2)에 대한 구조물이 파괴되는 개념을 설명한 그림으로 하중이 저항보다 큰 경우에는 파괴영역이 되고 작은 경우에서는 안전영역이 되며, 같은 경우에는 한계상태로 표현된다. 또한 한계상태는 선형식이나 비선형 곡선이 된다. 도 2에서 비선형곡선에 이르는 하중과 저항의 변수로 이루어 지는 응답면이 한계상태면에서 이루어 지는 경우를 나타내고 있다

도 3은 MCS와 IS의 정밀도를 비교한 그래프이다.

도 4는 본 발명에 적용된 해상풍력용 모노파일의 레이아웃을 나타낸 예시도이다.

도 5는 허용수평변위의 변화에 따른 파괴확률의 변화를 나타낸 그래프이다.

도 6은 허용회전각의 변화에 따른 파괴확률의 변화를 나타낸 그래프이다.

도 7은 모노파일의 직경 변화에 따른 파괴확률의 변화를 나타내 그래프이다.

도 8은 모노파일의 길이 변화에 따른 파괴확률의 변화를 나타낸 그래프이다.

도 9는 내부마찰각의 불확실성에 따른 신뢰도지수의 변화를 나타낸 그래프이다.

도 10은 확률변수(내부마찰각, 외력, 모노파일의 탄성계수)의 민감도를 나타낸 그래프이다.

도 11은 도 10에 도시된 확률변수의 변동계수 변화에 따른 민감도지수의 변화를 나타낸 그래프이다.

먼저, 도 1에 도시된 바와 같이, 본 발명의 실시 예에 따른 응답면 기법을 이용한 해상풍력용 모노파일의 신뢰성 해석 방법(S100)은 해상풍력터빈, 하부구조 모노파일의 거동에 관한 확률론적 해석에 응답면 기법(response surface method; RSM)을 연계하여 적용한 해석방법으로, 해상에 설치된 모노파일의 안전성을 확률적으로 분석하기 위한 고정 변수 및 한계상태별 확률변수를 설정하는 단계(S10); 및 상기 모노파일의 확률변수 범위에서 응답면 기법을 적용하여 한계상태함수를 확률변수에 대한 양함수의 형태로 근사화 하는 단계(S20); 및 근사화된 한계상태함수를 일계신뢰도법(FORM)에 적용하여 파괴확률(신뢰도지수)를 산출하는 단계(S30)를 포함하는 것을 특징으로 한다.

상기 고정 변수는 상기 모노파일의 연직하중, 길이, 직경, 수심 작용력, 수평 작용력, 지층조건을 포함하고, 상기 확률 변수는 상기 모노파일의 수평력, 모멘트, 탄성계수, 해저 기초지반의 점착력과 내부마찰각을 포함할 수 있다.

참고로, 현재까지 파일의 수평거동은 p-y 곡선에 의해 해석함이 가장 합리적이라 알려져 있으며, 이러한 비선형 p-y 거동해석은 유한요소법(FEM)에 의해 이루어진다.

여기서, 종속변수인 변위나 회전각 등 성능변수의 한계상태함수는 음함수(implicit function)의 형태이므로 직접 한계상태함수로 신뢰성 해석을 수행하기는 복잡하고 많은 시간이 요구되는 문제점이 있어, 본 발명에서는 응답면 기법을 파괴확률을 산출하는 데 사용한다.

이에, 응답면 기법(response surface method; RSM)을 이용하면 음함수 형태의 한계상태함수를 양함수(explicit function)의 형태로 근사화가 가능하고 여기서 얻어진 양함수를 이용하면 기존의 신뢰성 해석방법으로 가능하다.

한계상태함수를 근사화할 때 근사화된 한계상태함수

즉, 응답면 형태의 결정은 실제 한계상태함수

의 형태와 그 비선형성을 감안해야 한다. 음함수 형태에서는 대개

를 알 수 없으므로 현재까지 응답면은 광범위한 문제들에 걸쳐 적용할 수 있도록 일반적인 형태로 개발된 바 있다. 가장 일반적인 형태는 다음과 같은 2차 다항식이다.

[수학식 1]

여기서,

와

은 확률변수와 그 개수이고,

,

및

는 응답면의 형태를 나타내기 위한 계수이다.

만약,

가

보다 훨씬 높은 차수를 갖게 된다면 비정상적인 식이 도출될 수 있을 뿐 아니라 표본점(sample point) 영역 밖에서

와

사이에 큰 차이가 나타날 수 있으므로 대개 높은 차수의 다항식은 사용하지 않는다.

비록 2차 다항식을 사용함이 물리적 타당성은 부족하더라도 그것이 하나의 설계점(design point)을 쉽게 확인할 수 있도록 하는 단순하고 얻기 쉬운 매끄러운 면(surface)을 제공해 준다.

이는 특히 근사화된 한계상태면 상에 설계점을 위치시키는 응답면 근사화의 반복적인 설계와 근사화된 한계상태면 상에서 FORM과 SORM을 적용하는데 유용하다. 식에서

,

및

는

개의 미지의 계수이다. 이들 계수 값은 실제 한계상태함수

로 부터 일련의 표본점을 이용하여 특이값 분해(singular value decomposition)를 통해 결정할 수 있다.

표본점의 수는 적어도 계수의 수와 같거나 커야한다. 다양한 샘플링 방법이 있지만 통상적으로 에 나타난바와 같이

과

의

개 조합에서

를 평가한다. 여기서,

와

는

의 평균과 표준편차이고

는 임의의 계수이다.

실제 한계상태의 비선형성을 보다 정확하게 근사화하기 위하여 다음과 같이 2차 다항식에 혼합 항이 포함되기도 한다.

[수학식 2]

이렇게 근사화된 2차방정식 형태의 응답면 함수를 이용하면 기존의 일계신뢰도법(FORM)을 이용하여 신뢰성 해석이 가능하게 된다.

다음으로, 응답면기법을 이용한 신뢰성해석 결과는 Level 3의 중요도 추출법(Important sampling)을 이용한 결과와 상호 비교하여 결과의 타당성 여부를 검증하였다.

참고로, Level 3의 대표적인 추출법인 몬테카를로 시뮬레이션(Monte Carlo simulation; MCS)은 확률변수 발생빈도가 상대적으로 낮아질 경우 즉, 파괴확률이 작은 경우 정확한 파괴확률 추정을 위해 상당한 모의횟수가 필요하게 되며 그에 따라 많은 시간이 필요하게 된다.

이에 반해, 본 발명에서 사용되는 중요도 추출법(importance sampling; IS)은 확률변수의 발생이 파괴면 주위에서 많이 일어나도록 추출(sampling)에 사용하는 확률밀도함수를 조정함으로써 모의횟수 및 시간을 단축시킬 수 있는 방법이다.

보다 구체적으로, 주어진 확률분포

를 갖는 일련의 샘플

를 가정하면, 수학식 3에서

의 기대치의 형식은 수학식 4와 같이 그 샘플들로부터 산정된

의 평균으로 근사화할 수 있다. 여기서는 표기상 편의를 위하여 기대치를 대표하여

를 사용한다.

[수학식 3]

[수학식 4]

주어진

에서

의 값을 평가는 할 수 있지만, 산정에 필요한 분포로부터 샘플을 추출하기는 쉽지 않다. 대신에 샘플

를 추출할 일명 표본분포(sampling distribution)라는 또 다른 분포

를 도입한다.

[수학식 5]

중요도 추출법의 기본 개념은

와는 다르지만 유사한 분포, 소위

로부터 샘플을 추출한 다음 잘못된 분포에서의 샘플링으로 도입된 편의(bias)를 보정하기 위하여 얻어진 식을 수정하는 것이다.

한 점에서

를 평가할 수 있다고 가정하였으므로 수학식 5를 통하여 주어진

에 대하여 편의를 수정하거나 중요도 가중치

를 정확히 결정할 수 있다.

실제로

나

는 종종 비정규화된(unnormalized) 즉,

및

형태가 되며,

여기서,

[수학식 6]

이고, 수학식 6에서는 정규화 계수(normalization factor)를 산정하기 위하여

로부터 얻은 샘플을 이용한다.

만약, 중요도 가중치를

로 재정의 한다면 두 가지 형태의 근사화로 나타낼 수 있고, 여기서 수학식 8은 단순히 수학식 5에서 유도된 근사치이다.

[수학식 7]

[수학식 8]

수학식 7은

의 확률분포가 두 번 근사화된 근사식임에 주목할 필요가 있다(주요 항에서 한번, 그리고 수학식 6에서 중요도 가중치를 정규화하기 위하여 한번).

이러한 근사화는

나

가 정규화 된다는 보장이 없을 때 중요하다. 근사수하식 8의 두 번째 항은 두 분포가 정규화됨을 전제하고 통상 기대치에서

인 경우를 제외하고는 본래의 값과 같지 않게 된다.

결론적으로 파괴확률은 수학식 5와 같이 파괴가 발생하는 구간 내에서 확률밀도함수를 적분함으로써 얻을 수 있다. 중요도 추출법에서는

를 도입함으로써 파괴확률을 구하게 된다.

는 파괴면상의 값(설계점)을 평균으로 하는 가상의 확률밀도함수이며, 이 값은 FORM 등의 근사해법을 이용하여 얻을 수 있다(도 3 참조). 도 3은 MCS와 IS의 정밀도를 비교한 예로서 가로축은 모의횟수, 세로축은 확률을 나타낸다.

상대적으로 모의횟수(N)에 관계없이 IS의 오차가 적게 나타나며, 추출수(N)가 많을수록 정도가 향상됨을 알 수 있다.

이하에서는 해상풍력 모노파일의 신뢰성 해석을 실제 시공된 사례분석을 통해 설명하도록 한다.

1.해석조건

해상조건에서 풍력터빈 구조물 기초의 안정성을 확률론적으로 분석하기 위하여 최근 Jonkman et al.(2009)이 해상시스템에 적용할 목적으로 개발한 5MW 규모의 풍력터빈 제원을 이용하였다. 개발된 풍력터빈은 나셀, 블레이드, 타워 및 기초직경의 제원이 제시되어 있다.

[표 1]

표 1은 신뢰성 분석에 사용되는 모노파일의 환경조건을 나타낸 표이다.

따라서, 표 1과 도 4에 나타난 바와 같이 기초설계에 필요한 연직하중은 터빈구조물의 제원으로부터 산출하였고 수심, 수평 작용력과 지층조건은 임의로 가정하되 수평 작용력은 Bush et al.(2009)이 적용한 수평하중을 인용하였다.

특히, 기초지반은 조밀한 모래가 해저면으로부터 35m 이상의 심도에 존재하고 수심은 15m가 유지되는 것으로 가정하였다.

확률론적 해석에서 작용력(자중, 외력 등)과 저항요소(재료특성 등)의 불확실성 정량화는 매우 중요하다. 특히, 그 중에서도 확률분포와 변동성은 해석결과에 주요하게 영향을 미치는 인자이다. 확률론적 해석에 적용한 확률변수들과 그 특성을 표 2에 나타내었다.

본 실험에서는 여러 가지 설계변수 가운데 하중으로서 수평력과 모멘트, 저항요소로서 파일의 탄성계수와 지반의 내부마찰각을 확률변수로 취급하였다.

모든 확률변수의 분포는 정규(normal)분포를 따르고 표 2에 나타난 바와 같이 5∼15%의 불확실성을 가지는 것으로 가정하였다. 표 2는 임의의 변수들의 특성을 나타낸 표이다.

[표 2]

2. 신뢰성 해석과 매개변수 분석

모노파일의 두부 수평변위에 대한 파괴확률 또는 신뢰도지수를 계산하기 위하여 응답면 기법(RSM)과 중요도 추출법(IS)을 이용하였다. IS는 응답면 주변에서 10,000개의 난수를 생성하여 한계상태함수를 계산한 후 파괴확률을 산정하였고 계산된 한계상태함수 값들의 변동성이 0.01% 이내가 될 때까지 반복하여 정밀도를 확보하였다.

신뢰성 해석은 한국해양과학기술원(KIOST)에서 개발한 신뢰성 설계 프로그램 HSRBD ver3.1(2013)을 이용하였다. HSRBD(Harbor Structure Reliability Based Design)은 p-y 모델을 이용하여 말뚝의 신뢰성 설계가 가능한 프로그램으로서, 말뚝 이외에도 방파제와 안벽, 호안 등 항만구조물을 허용응력 및 신뢰성 기법을 이용하여 해석할 수 있는 모듈을 갖추고 있다.

신뢰성 해석기능은 부분안전계수법(레벨 1)으로부터 FORM, RSM의 레벨 2 및 MCS 및 IS를 포함한 레벨 3까지 다양하게 해석이 가능하다.

여기서, 본 발명에서 제공된 모노파일 구조물은 모래지반에 설치되는 파일기초로서 가상 고정식 내지는 등분포 스프링식 해석모델을 이용하므로 p-y 곡선을 이용하여 변위와 회전각을 계산하게 된다.

따라서, 여기서는 사질토 지반에서 p-y 모델로 많이 사용되는 API(2005) 방법과 쌍곡선 모델(Dewaikar et al., 2006)을 적용하였다. 또한, 허용응력 설계법에서 파일 두부의 수평변위와 회전각에 대한 안정여부는 각각 허용변위 및 허용 회전각 기준과 계산된 수평변위 및 회전각을 비교하여 판단한다.

따라서, 한계상태함수는 다음과 같이 정의할 수 있다. 여기서,

와

는 각각 외력에 의해 발생되는 수평변위와 회전각이고

와

는 허용 수평변위 및 회전각으로서 과거 연구사례(DNV, 2007; Kuo et al., 2008)를 참고하여

,

로 설정하였다.

[수학식 9]

[수학식 10]

응답면 기법을 적용하기 위하여 수학식 9와 수학식 10의 한계상태함수를 응답면 함수로 변환하면 다음과 같다.

여기서,

는 모래의 내부마찰각,

는 파일의 탄성계수이고

은 구조해석을 통하여 회귀분석으로 구해지는 계수이다.

[수학식 11]

[수학식 12]

응답면 기법(RSM)과 중요도 추출법(IS)을 이용한 해석결과를 표 3에 비교하였다.

두부 수평변위(

)가 회전각(

) 파괴모드에 비하여 신뢰도지수가 작게 산정되어 지배 파괴모드인 것으로 나타났다.

API(2005) 기준에 제시된 p-y 함수를 이용한 결과와 비교할 때 쌍곡선(hyperbolic) 함수를 이용한 경우 신뢰도지수는 14∼18% 포인트 가량 낮게 산정되었고 이는 구성식의 차이에서 기인한다.

응답면 기법(RSM)과 중요도 추출법(IS)를 이용한 신뢰도지수는 오차가 약 0.3% 이내로 거의 일치하는 결과를 보였다. 이는 중요도 추출법(IS)을 이용한 추출 시 이미 응답면 기법(RSM)에서 근사화된 파괴점 주변에서 값들을 추출하기 때문인 것으로 판단된다.

[표 3]

수학식 9과 수학식 10의 한계상태함수에 나타났듯이 허용 수평변위와 허용 회전각이 각각의 파괴모드에서 저항에 해당하므로 허용 변위 및 회전각 기준의 설정에 따라 파괴확률은 달라지게 된다. 따라서, 이들을 변화시키며 해석함으로써 계산된 파괴확률의 변화를 관찰하였다.

도 5에 나타난 바와 같이 허용변위가 15mm에서 100mm까지 변화할 때 파괴확률은 거의 대수적(logarithmic)으로 감소하는 것을 알 수 있다. 국내의 말뚝설계 시 이용되는 규정 가운데 파일직경의 1%에 해당하는 60mm에서 파괴확률은 약 1.5%인 반면, 최소 기준변위인 15mm로 설정한 경우의 파괴확률은 약 60%로서 매우 큰 것을 알 수 있다.

도 6에 도시된 바와 같이, 허용 회전각을 0.2∼0.6°까지의 범위로 설정할 때의 파괴확률로서 역시 대수적으로 감소하나 그 범위는

%로 수평변위 파괴에 비하여 매우 작은 범위에서 변화하였다.

모노파일의 직경은 수평지지력뿐 아니라 변위 감소효과와 경제성에도 민감한 요소이므로 최적설계에서 파일의 직경은 주요한 고려사항이다. 따라서, 모노파일의 직경을 변화시키며 해석하여 파괴확률 변화를 관찰하였다.

도 7에 도시된 바와 같이, 직경이 증가할수록 파괴확률이 감소하며 그 감소폭은 커져 지배 파괴모드인 수평변위의 경우 7.0m부터는 직경을 증가시켜도 최대 파괴확률의 차이는 0.27% 포인트 이하로 나타남으로써 변위 감소효과가 적어지는 것으로 나타났다.

이는 모래지반에서 모노파일이 일정한 직경이상 확보되면 두부의 응력이 점차 하부로 전이되면서 나타나는 현상으로 판단된다. 아울러, 회전각의 경우는 직경증가에 따라 파괴확률이 급격하게 작아짐으로써 상대적으로 직경변화에 따라 민감함을 알 수 있다.

모노파일의 길이는 기반암층의 심도, 지층분포와 수심 등에 영향을 받으며 경제성에 큰 영향을 미친다. 따라서, 파일의 길이에 따라 파괴확률의 변화를 관찰하였다.

도 8에 도시된 바와 같이, 파일의 길이가 증가함에 따라 파괴확률은 감소하고 그 감소폭은 점차 작아진다. 두부 수평변위와 회전각 파괴 모두 파일길이 60m부터 파괴확률의 감소폭은 매우 작아져 그 차이는 약 0.07% 이하로 나타난다. 즉, 파일길이 60m부터는 길이가 증가하여도 변위 감소효과는 미미함을 알 수 있다.

토질정수의 불확실성은 비교적 크나 모래의 내부마찰각에 대한 변동성은 일반적으로 10% 내외인 것으로 알려져 있다. 내부마찰각의 변동성에 따른 민감도를 파악하기 위하여 변동계수(COV)를 감소시키며 해석하여 그 결과를 도 8에 나타내었다.

다른 해석조건은 동일한 상태에서 COV를 당초 15%에서 10%로 감소하여 해석한 결과, 신뢰도지수(β)는 약 0.35% 포인트 증가하는데 그쳤다. 반면 10% 이하의 COV 범위에서는 매우 다른 양상을 나타내었다. 즉, 내부마찰각의 COV 8%에서 신뢰도지수 62.6% 포인트 증가, COV 6%에서는 72.5% 포인트 증가하는 등 내부마찰각의 COV 10% 이내에서는 신뢰도지수 상승 폭이 매우 큰 것으로 나타났다.

이는 모래지반에 모노파일이 설치되는 경우 내부마찰각의 변동성 수준이 파괴확률과 공사비에 큰 영향요인임을 나타내는 것이다.

3. 확률변수의 민감도 분석

2차식 형태로 근사화된 응답면 함수를 가지고 FORM에 의해 해석하는 과정에서 신뢰도지수와 민감도지수(sensitivity factor)를 얻을 수 있다. 수학식 13의 신뢰도지수(reliability index,

)와 각 확률변수의 민감도지수(

)의 관계로부터 두 변수의 반복계산을 통하여 일정한 값에 수렴하는

와

를 결정한다(Yoon et al., 2010).

식에서

는 통계적으로 독립인 확률변수,

는 표준정규분포로 선형 변환된 확률변수,

는 한계상태함수를 의미한다. 수학식 13은 확률변수 사이에 서로 상관성이 없을 때의 관계이며,

는 표준화된 공간상의 각 확률변수 축에서 신뢰도지수의 방향여현(direction cosine)을 의미한다.

[수학식 13]

또한, 모든 확률변수의

는 수학식 14의 관계를 만족하도록 결정해야 한다.

[수학식 14]

민감도지수는 각 확률변수가 파괴확률에 미치는 기여도를 의미하며, 일반적으로 외력의 변수에 대하여 양(+)의 값을, 저항의 변수에 대하여 음(-)의 값을 나타낸다.

도 10에 도시된 바와 같이, 본 발명에서 확률변수로 취급한 세 가지 변수 가운데 흙의 내부 마찰각에 대한 민감도가 지배적인 것으로 나타났으며 다음으로 외력(수평력 및 모멘트)과 파일의 탄성계수 순이었다.

이는 흙의 내부마찰각이 파괴확률에 미치는 영향이 가장 크므로 설계시 지반조사 및 시험의 정도를 높여 불확실성을 최소화 하여야 경제성을 제고할 수 있음을 의미하기도 한다.

각 확률변수의 변동계수를 변화시키며 민감도지수(

)의 변화를 관찰한 결과 도 11에 도시된 바와 같이, 내부마찰각에 대한 민감도지수는 1에 가까운 값으로 나타나 큰 변화가 없는 것으로 나타났다. 나머지 확률변수의 경우는 변동성이 작아질수록 저항관련 확률변수는 0으로부터 1로, 하중관련 확률변수는 0으로부터 -1로 근접해 간다. 이는 확률변수의 불확실성이 클수록 파괴확률에 미치는 영향이 큼을 의미한다.

따라서, 본 발명의 실시 예에 따르면 수평변위 및 회전각에 대해 허용치를 증가시키는 경우에 모노파일 기초의 파괴확률은 대수적으로 감소하였다.

국내 말뚝설계 규정인 파일직경(D=6m)의 1%를 허용수평변위(60mm)로 적용하는 경우 파일의 파괴확률은 1.5%이고 최소기준인 15mm에 대해서는 60%로서 그 차이가 크므로 해상풍력단지에 적용하기 위한 파일의 수평변위 및 회전각에 대한 허용기준의 정립이 시급한 과제인 것으로 평가되었다.

또한, 모노파일의 직경이 커지면 수평변위 및 회전각에 대한 파괴확률은 감소하나 직경이 7.0m이상의 경우 점차 수렴하여 파괴확률의 저감효과는 미미하였다. 이는 모래지반에서 모노파일이 일정한 직경이상 확보되면 두부의 응력이 점차 하부로 전이되면서 나타나는 현상으로 판단된다. 아울러 파일직경 변화에 대한 영향은 파일 두부의 수평변위에 비하여 회전각이 보다 큰 것으로 나타났으며 이를 최적설계에 활용할 수 있다.

또한, 해저면으로부터 기반암층의 심도가 깊어 파일이 모래지반에 근입되는 마찰말뚝인 경우 파일길이가 증가함에 따라 수평변위에 대한 파괴확률은 감소하나 그 감소폭은 작아지며, 파일길이가 일정 근입깊이 이상이 되는 조건, 즉 60m 이상이 되는 경우 수평변위의 감소효과는 작게 평가되었다.

또한, FORM을 이용한 민감도 분석결과, 변동성이 큰 사질토 지반의 본 사례에서는 흙의 내부마찰각이 파일의 거동에 가장 큰 영향을 미치는 것으로 나타났다. 따라서, 불확실성이 클 것으로 예측되는 기초지반에서는 정밀한 지반조사뿐 아니라 변동성 분석이 경제성과 안정성을 최적화하는 측면에서 매우 중요한 것으로 판단된다.

이상 본 발명이 양호한 실시 예와 관련하여 설명되었으나, 본 발명의 기술 분야에 속하는 자들은 본 발명의 본질적인 특성에서 벗어나지 않는 범위 내에 다양한 변경 및 수정을 용이하게 실시할 수 있을 것이다. 그러므로 개시된 실시 예는 한정적인 관점이 아니라 설명적인 관점에서 고려되어야 하고, 본 발명의 진정한 범위는 전술한 설명이 아니라 특허청구범위에 나타나 있으며, 그와 동등한 범위 내에 있는 모든 차이점은 본 발명에 포함된 것으로 해석되어야 할 것이다.

S10: 변수 설정 단계
S20: 응답면기법을 이용하여 한계상태함수를 양함수 형태로 근사화하는 단계
S30: 신뢰성 해석을 통한 파괴확률 산출 단계

Claims

해상풍력터빈 하부구조 모노파일이 해저 기초지반에 설치된 경우에 확률론적 해석을 통하여 상기 모노파일의 거동 확률분석 및 이에 따른 민감도를 분석하는 응답면 기법을 이용한 컴퓨팅 장치에 의한 해상풍력용 모노파일의 신뢰성 해석 방법에 있어서,
해상에 설치된 상기 모노파일의 안전성을 확률적으로 분석하기 위한 고정 변수 및 한계상태별 확률변수를 상기 컴퓨팅 장치가 설정하는 변수 설정 단계(S10);
상기 모노파일의 확률변수 범위에서 응답면 기법을 적용하여 상기 컴퓨팅 장치가 한계상태함수를 확률변수에 대한 양함수의 형태로 근사화 하는 단계(S20); 및
근사화된 한계상태함수를 일계신뢰도법 (FORM)에 적용하여 상기 컴퓨팅 장치가 파괴확률(신뢰도지수)를 산출하는 단계(S30)를 포함하는 것을 특징으로 하는 응답면 기법을 이용한 해상풍력용 모노파일의 신뢰성 해석 방법.
제1항에 있어서,
상기 고정 변수는,
상기 모노파일의 연직하중, 길이, 직경, 수심 작용력, 수평 작용력, 지층조건을 포함하는 것을 특징으로 하는 응답면 기법을 이용한 해상풍력용 모노파일의 신뢰성 해석 방법.
제1항에 있어서,
상기 확률 변수는,
상기 모노파일의 수평력, 모멘트, 탄성계수, 해저 기초지반의 점착력과 내부마찰각을 포함하는 것을 특징으로 하는 응답면 기법을 이용한 해상풍력용 모노파일의 신뢰성 해석 방법.
제1항에 있어서,
상기 근사화된 한계상태함수는,
2차 다항식 형태의 수학식 1을 통해 도출되며,
[수학식 1]

여기서,
와
은 확률변수와 그 개수이고,
,
및
는 응답면의 형태를 나타내기 위한 계수이며,
,
및
는
개의 미지의 계수로서, 이들 계수 값은 실제 한계상태함수
로부터 일련의 표본점을 이용하여 특이값 분해(singular value decomposition)를 통해 결정되는 것을 특징으로 하는 응답면 기법을 이용한 해상풍력용 모노파일의 신뢰성 해석 방법.