J. Ocean Eng. Technol. Search

CLOSE


J. Ocean Eng. Technol. > Volume 27(2); 2013 > Article
반응표면법 및 비지배 분류 유전자 알고리즘을 이용한 취배수문의 응력 및 변형 최적화

Abstract

A valve is a marine structure that is subjected to multiple seawater loads. Therefore, it is necessary to define the kind of loads applied to it to confirm whether the structure has sufficient strength. In this research, we aimed to find the optimal solution for the stress and deformation of valves under various loads. We first selected design variables and implement a finite element analysis according to changes in the thickness of each component of a valve based on a central composite design. Next we developed a regression model of the response surface. Using this model, we calculated the optimal objective value based on NSGA-II. Finally, to confirm the correspondence between the optimal objective value and the real FEM value, we compared the optimal result and structural analysis result to verify the performance of NSGA-II.

1. 서 론

항만시설은 외력(파랑, 태풍, 지진에 의한 동수압 등)이 크게 작용하기 때문에 대형 구조물로 이루어진다. 파랑, 태풍 또는 지진등은 실험이 불가능하기 때문에 대부분 전산 프로그램을 이용해서 구조물의 안전성을 평가한다(Hur et al., 2005). 그러나 대형 구조물의 전산구조해석을 수행할 경우 많은 해석 시간과 비용이 소요되기 때문에 최소의 구조해석 작업으로 최적화하는 것이 필요하다.
항만시설의 하나로써 취배수문(Culvert gate)은 갑문(Lockgate)을 이용해서 해측과 내측의 수위를 일정하게 유지하기 위해서 물이 지나가는 통로에 설치되는 장비이다. 일반적인 해양 장비와 같이 취배수문 또한 여러 종류의 하중이 동시에 작용하게 되므로 각각의 개별하중 및 하중 조합은 구조해석 시 중요한 요소가 된다. 해수의 유동운동에 의해 하중이 작용하기 때문에 동적 유동해석은 반드시 필요하다. 그러나 대형 구조물에 대한 유동해석은 상당한 시간과 노력을 필요로 하기 때문에 이러한 문제점을 해결하기 위해서 동적 운동을 정적운동으로 근사화하는 방법들이 제시되고 있다. Naudascher et al.(1986)은 취배수문이 닫혀 있다가 상승할 경우 바닥과 상부의 해수유동의 압력차이에 의해서 발생하는 견인력(Downfull force)을 정적하중의 식으로 근사화 하였고, Westergaaed(1933)는 지진에 의해서 발생하는 지진 동수압을 등가정적(Equivalent-static) 해석법을 이용하여 근사화한 식을 제공하고, 구조해석 시 유한요소해석을 수행하였다.
취배수문과 같은 구조물의 안전성을 평가함에 있어서 응력과 변형은 아주 중요한 요소이기 때문에 본 논문에서는 취배수문에 작용하는 응력과 변형의 근사최적화를 수행함으로써 많은 시간과 비용을 줄이고자 하였으며 근사 최적화를 수행 후 얻은 결과값을 유한요소해석 결과값과 비교하였다.
본 논문에서는 다양한 하중이 작용하는 조건에서 취배수문의 무게증감이 한계 범위를 벗어나지 않도록 응력과 변형을 최적화하고자 한다. 응력과 변형은 목적함수로 정의하고 제한조건은 무게로 정하였다. 먼저 중심합성계획을 통해 2차 회귀모형의 반응표면을 생성하고, 생성된 근사식들은 다중목적 함수 최적화를 수행하기 위해서 다양한 분야에서 사용되어 그 유용성이 증명된 NSGA-II(non- dominated sorting genetic algorithm II)에 적용되었다.

2. 취배수문의 구조 및 하중 종류

2.1 취배수문의 구조

400TEU 이상의 큰 용량의 배가 내항으로 들어 올 경우 내항과 외항의 수위차를 일정하게 유지시키기 위해서 워터펌프(Water pump)을 이용해서 수위를 조절한다. 이 때 취배수문은 해수가 드나드는 통로 중간에 갑문과 함께 설치되는 장비로써, 수위가 달라질 경우에는 해수 통로를 막고 있던 취배수문이 리프팅(Lifting)되어 해수가 드나들게 되고, 시간이 지난 후 수위가 일정해지면, 다시 하강하게 되어 통로를 막게 되는 구조이다. 형상은 설치 조건에 따라서 다르지만 설치형태는 대체로 Fig. 1과 같다.
Fig. 1

Layout drawing of valve

HOGHC7_2013_v27n2_27_f001.jpg

2.2 취배수문의 하중조건

취배수문에 작용하는 개별 하중은 정수압, 동수압, 지진동수압, 지진하중 등이 있다. 정수압은 작업 시와 유지보수 시 작용하는 수압이 달라지므로 세분화하였다.
2.2.1 개별하중의 종류
① 작업 시 정수압(Hydrostatic force, operating)
작업 시 취배수문을 중심으로 좌우 다른 수압이 작용하기 때문에, 최종적으로 문비에 작용하는 외력은 좌우 수압의 대수적인 차이만큼 작용하게 되며, 수압은 Fig. 2와 같다.
Fig. 2

Hydrostatic pressure in service

HOGHC7_2013_v27n2_27_f002.jpg
HOGHC7_2013_v27n2_27_e901.jpg
여기서 P는 작업 시 취배수문에 작용하는 최대 수압이고, W는 해수의 단위체적당 중량을 나타내며, h는 문짝의 상·하류에서의 파랑을 고려한 수위에서 정수압이 작용하는 하부수심까지의 거리를 나타낸다.
② 유지보수 시 정수압(Hydrostatic force, maintenance)
수문(Lock-gate)의 유지보수 시 한쪽 면의 해수는 모두 배수 되어야 하므로, 수압은 한쪽 문비(Skin plate)에 대해서만 작용하게 되며, 수압은 Fig. 3 과 같다.
Fig. 3.

Hydrostatic pressure out of service

HOGHC7_2013_v27n2_27_f003.jpg
HOGHC7_2013_v27n2_27_e902.jpg
③ 동수압(Hydraulic downpull force)
동수압이 발생하게 되는 원인은 취배수문을 중심으로 좌우 다른 수압이 작용하고 있을 때, 취배수문이 올라가게 되면 해수의 유동에 의해서 위아래 수압차가 발생하고 순간적으로 아래로 당기는 힘이 발생한다. 동수압 발생에 영향을 미치는 인자는 취배수문 하부의 기하학적 형상과 상부를 지나는 유량에 의해서 결정된다(Naudascher et al., 1986; Weaver and Martin, 1980).
HOGHC7_2013_v27n2_27_e903.jpg
여기서 F는 Downpull force가 되며, αT는 취배수문 상부의 Downpull 계수이고, αB는 취배수문 하부의 Downpull 계수, B는 취배수문의 폭, d는 취배수문의 깊이, ρ 는 밀도 그리고 V는 취배수문의 하부를 지나는 유체의 속도를 나타낸다. 식 (3)의 계수들을 구하기 위해서 Naudascher et al.(1986)의 실험그래프를 사용하고자 한다. Fig. 4는 일반적인 취배수문 구조도를 나타내고 있다. Y는 유체가 흐르는 통로이므로 고정값이 되고 Y0는 취배수문이 상하운동을 하면서 변하는 값이 된다. 유체 유속은 취배수문의 기하학적 요소(d, t, θ)들에 의해서도 영향을 받지만, 여기서는 고려하지 않았다. Fig. 5는 취배수문의 상하운동에 대한 Downfull force계수(αT-αB)변화를 보여주고 있다. 취배수문이 약 40% 정도 개방되었을 때 가장 높은 계수값이 나타내기 때문에 Y/Y0가 0.4일때 Downfull force의 계수 0.4를 이용하여 동수압을 구하였다.
Fig. 4

Diagram of tunnel valve under submerged flow condition (Naudascher et al., 1986)

HOGHC7_2013_v27n2_27_f004.jpg
Fig. 5

Resultant downpull and discharge coefficients versus gate opening (Weaver and Martin, 1980)

HOGHC7_2013_v27n2_27_f005.jpg
④ 지진동수압(Seismic hydrodynamics)
지진동수압은 지진에 의한 지반운동에 의해서 해수가 동적 움직임을 발생시켜 취배수문에 수압을 발생하므로 동적응답해석(Dynamic response analysis)방법을 수행하여야 한다. 그러나 동적 응답해석 시 상당한 시간이 소요되므로 본 논문에서는 이러한 방법대신 동적운동을 정적운동으로 근사한 등가정적방법(Equivalent static method)을 사용해서 외력을 구한다. 이러한 근사화 방법은 Westergaard(1933)에 의해서 제안이 되었으며 식은 (4)와 같다.
HOGHC7_2013_v27n2_27_e904.jpg
여기서 Pu는 수면으로부터 깊이 y까지의 수압을 나타내며, α, w, H, y는 각각 수평지진가속도, 해수의 단위중량, 해수와 접촉하고 있는 취배수문의 높이, 수면으로부터의 깊이를 나타낸다. 식 (4)는 댐(Dam)과 같이 연성(Flexibility)이 거의 없는 물체(Body)에 적용이 되는 식이므로 철 구조물로 이루어진 취배수문에 바로 적용할 수 없다. 그래서 연성을 고려하여 수정된 식 (5)를 사용하며, 이 식은 Westergaard(1933)에 의해서 제시 되었다.
HOGHC7_2013_v27n2_27_e905.jpg
여기서 H1, H2는 물과 접촉하고 있는 취배수문의 높이를 나타내며, 아래첨자 1, 2는 취배수문의 외항편, 내항편을 각각 나타낸다. T는 진동주기이며, MF(Magnification factor)는 지진에 의한 연성 구조물에 발생하는 가속도비이다. 구조물은 강체(Rigid body)보다 연성일 경우에 가속도에 의한 동적응답은 연성 구조물이 강체 구조물보다 크기 때문에 MF값을 고려하여 전체동수압을 증가시키게 된다. 진동주기(T)는 유한요소해석을 수행하여 얻었으며 해석결과는 1.35sec이다. 이 결과값에 의해 MF 값은 5로 설정되었다.
⑤ 지진하중(Seismic inertial force)
지진 하중은 지진에 의해서 구조물이 흔들리는 현상으로 외력은 식 (6)과 같이 구할 수 있다. 여기서 Ef는 지진하중이고 R, D, g는 각각 지진구역계수, 위험도 계수, 중력가속도를 나타내며 각각 0.11, 1.4, 그리고 9.8m/s2으로 설정하였다.
HOGHC7_2013_v27n2_27_e906.jpg
2.1.2 하중 조합
개별 하중들(①~⑤)을 모두 구한 후 각각의 개별 하중을 합리적으로 조합하여 조합 하중을 만들어야 한다. 본 논문은 하중조합을 사용환경의 하중조합에 따라서 적합하게 구성하였으며, 하중 조합을 Table 1과 같이 구성하였다. 하중 조합은 기본적으로 취배수문의 자중(Self-weight)과 정수압이 항상 작용하고 있기 때문에 모든 경우에 대해서 고려하였다. 정수압이 작용하고 있는 상태에서 취배수문이 개방(Opening)되면 동수압이 작용하기 때문에 Case 3과 같이 조합하였으며, 정수압 조건에서 지진이 발생할 경우 지진 동수압 및 지진 관성력이 발생하기 때문에 하중 조합을 Case 4와 Case 5와 같이 구성하였다.
Table 1

Analysis parameters

HOGHC7_2013_v27n2_27_t001.jpg

3. 유한요소 해석

Table 1의 5개의 Case에 대해서 유한요소해석을 수행하여 응력값과 변형값을 도출한 후 각 Case 중 가장 큰 값이 나타나는 경우에 대해서 최적화를 진행하였다. 유한요소해석에 사용된 모델은 Fig. 6과 같고, 5개의 Case에 대한 해석결과값(응력과 변형값)은 Table 2와 같다. Table 2에서 Case 4가 응력과 변형값이 가장 크게 발생하였으며 그 결과를 Fig. 7에서 보여주고 있다. Fig. 7은 정수압+지진동수압(Z dir.)의 하중조합을 이용한 결과를 나타내고 있으며, 유한요소 해석 조건은 프아송비는 0.3, 탄성계수는 2.0904e6kg/cm2, 요소(Element)의 종류는 Shell63, 절점(Node)와 요소 수는 각각 8,605개 8,960개로 하였다. 해석결과는 최대응력 2,940kg/cm2, 최대변형량(Z dir.) 0.584cm이며 현재 취배수문에 사용되는 재질은 SM490YA로써 항복응력은 3,600kg/cm2이다. 따라서 Case 4의 최대응력이 항복응력을 초과하지 않았기 때문에 취배수문은 탄성변형영역에 존재한다.
Fig. 6

FEM model using ANSYS (Inman, 2007)

HOGHC7_2013_v27n2_27_f006.jpg
Table 2

FEM analysis results

HOGHC7_2013_v27n2_27_t002.jpg
Fig. 7

FEM (finite elements method) analysis result for Case 3 (using ANSYS)

HOGHC7_2013_v27n2_27_f007.jpg

4. 최적 설계

4.1 설계변수 정의

취배수문에는 여러 종류의 플레이트(Plate)가 사용 되지만, 외력에 대한 응력 및 변형이 발생하는 부재는 Fig. 8과 같이 Skin-plate, Reinforce-plate, Side-plate 세가지로 분류할 수 있다. 여기에서,Sskin-plate는 취배수문을 중심으로 좌우 수압을 유지하도록 해수 통로를 막는 역할을 하며, Reinforce-plate는 수압을 받고 있는 Skin-plate을 지지하게 된다. 그리고 Side-plate는 휠(Wheel)이 부착되어 있는 부재로써 취배수문이 위아래로 이동할 경우 양측면에 Wheel이 부착된 수직통로를 따라서 이동한다.
Fig. 8

Design variables

HOGHC7_2013_v27n2_27_f008.jpg

4.2 최적설계 정식화

취배수문의 목적함수 및 제한조건함수는 식 (7)과 같이 정의하였다. 아래의 조건식 중 변형량 허용치는 한국수자원공사의 권고치(Korea water resources corporation, 2003)를 참고로 하였다. 그리고 제한조건 최대 무게는 사용 모터의 리프팅(Lifting) 최대 용량인 16ton으로 하였다.
HOGHC7_2013_v27n2_27_e907.jpg

4.3 최적화 과정

설계변수들과 목적함수, 제한조건함수의 근사식을 생성하기 위해서 먼저, 중심합성계획법을 통해서 설정된 설계변수(두께)들로 유한요소해석을 실행해서 응력과 변형값을 구하게 된다. 반응표면법에 의해 2차 회귀방정식을 구성하고, 이 2차 회귀방정식을 비지배 분류 유전자 알고리즘(NSGA-II)을 사용하여 목적함수의 최적해를 찾는다. 이 최적해를 유한요소 해석값과 비교하여 정확성을 검증하였다.
4.3.1 중심합성계획법
중심합성계획법은 적은 횟수의 실험으로 곡선을 추정하기 위해서 중심점과 축점을 2k요인실험에 추가시킨 구조로써, 이 방법의 장점은 2k요인 실험계획에 의하여 실험을 한 후에 반응표면이 1차 회귀모형으로 적당하지 않을 경우에 2k요인에 중심점과 축점을 증가시키면 쉽게 2차 회귀모형을 생성할 수 있는 장점이 있다. 총 시험 횟수를 식 (8)에 의해 구할 수 있으며 k는 설계변수이고 n0은 중심점의 수이다.
HOGHC7_2013_v27n2_27_e908.jpg
본 논문에서는 취배수문의 두께 3곳이 설계변수이므로 총 실험횟수는 n=15, 알파값은 1.68로 하였다.
4.3.2 반응표면법
반응표면법은 여러개의 설계 변수가 복합적으로 작용함으로써 어떤 반응변수에 영향을 주고 있을 때 이러한 반응의 변화가 이루는 반응표면에 대한 통계적인 분석방법이다. 본 논문에서는 2차 회귀모형을 사용하였다.
HOGHC7_2013_v27n2_27_e909.jpg
만약, 설계변수(k)가 3개인 경우 추정량에 대한 방정식은 식 (10)과 같다.
HOGHC7_2013_v27n2_27_e910.jpg
Table 3Table 4을 기반으로 해서 2차 회귀방정식을 만들면 식 (11), (12), (13)과 같다.
Table 3.

Central composite design

HOGHC7_2013_v27n2_27_t003.jpg
Table 4.

FEM results

HOGHC7_2013_v27n2_27_t004.jpg
HOGHC7_2013_v27n2_27_e911.jpg
HOGHC7_2013_v27n2_27_e912.jpg
HOGHC7_2013_v27n2_27_e913.jpg
Table 5

Fitness

HOGHC7_2013_v27n2_27_t005.jpg
4.3.3 적합도 평가
중심합성계획법(Na et al., 2005; Park et al., 2005)을 이용하여 반응표면 모델을 구성하였으며, 반응표면 모델의 통계적인 신뢰성을 검증하기 위해서 결정계수 R2을 주로 사용한다. 그러나 결정 계수는 그 중요성의 유무와 상관없이 변수가 증가할수록 높아지는 단점을 가지고 있으므로 이를 보완하기 위해 조정결정계수 R2adj을 사용한다. 조정결정계수는 불필요한 변수가 추가되었을 때 그 값이 감소함으로써 필요한 변수들만을 선택하는데 유용하게 사용된다. 결정계수는 식 (14)와 같고 n은 실험횟수이고 p는 반응표면 모델 계수의 개수이다.
HOGHC7_2013_v27n2_27_e914.jpg
4.3.4 최적화 결과
NSGA-II (Deb et al., 2002)는 빠른 엘리티스트 비지배 분류 유전자 알고리즘(Fast elitist non-dominated sorting genetic algorithm)으로도 불리며, 기존에 사용되던 비지배 분류 알고리즘(Horn et al., 1994)의 단점(계산의 복잡성, 엘리트즘의 부족, 공유 파라미터 지정의필요성)을 개선한 알고리즘으로 본 논문에서는 세대수 300, 인구수 150으로 설정하였다. 반응표면법에 의해 생성된 목적함수 응력과 변형, 그리고 제한조건함수 무게의 근사식은 식 (13)과 같다. 다중목적함수 응력(f1)과 변형(f2)의 최적화를 위해 NSGA-II를 적용하였으며, 최적화 결과로 얻은 pareto 해를 Fig. 9에서 보여주고 있다. 이 해의 정확성을 검증하기 위해서 양극점과 중간점을 선택해서 정확성을 비교한 결과는 Table 6과 같이 나타났다.
Fig. 9

Result of NSGA-II using Matlab (Attaway, 2009)

HOGHC7_2013_v27n2_27_f009.jpg
Table 6

Compared with NSGA-II and FEM results

HOGHC7_2013_v27n2_27_t006.jpg
기존의 취배수문의 기본 두께인 경우와 최적해의 중간값의 두께인 경우를 비교하였으며 Table 7과 같다. 두께의 변화는 x1의 경우 2.2cm의 두께에서 2.12cm 로 감소, x2의 경우는 1.50cm에서 1.78cm로 증가, 그리고 x3의 경우 2.40cm에서 2.00cm로 변경되었다. 부재 두께에 따른 응력과 변형량을 비교해 보면 기존의 기본 두께에서는 2,940kg/cm2와 0.584cm로 나타났으며, 최적화를 통해서 나온 값은 각각 2,875kg/cm2와 0.439cm로 나타났다. 그리고 취배수문의 전체 무게 변화량은 기본 두께에서 14,499kg, 최적화가 이루어진 상태에서 15,304kg으로 나타났으며, 제한조건을 만족하였다.
Table 7

Compared with standard and optimization value

HOGHC7_2013_v27n2_27_t007.jpg

5. 결 론

본 연구에서는 취배수문의 응력과 변형을 최적화하기 위해서 응력과 변형을 다중목적함수로 설정하고 무게를 제한조건함수로 설정하여 근사최적설계를 하였다. 중심합성계획법을 이용하여 반응표면을 생성하였으며 목적함수, 제한조건함수 모두 99%이상의 높은 적합도를 나타내었다. 반응표면법에 의해 구한 2차 근사모델과 NSGA-II를 이용하여 다중목적함수를 최적화하여 Pareto 해를 얻었으며, 이 Pareto 해 중에서 중간에 해당하는 해를 FEM 해석 결과와 비교하였다. 응력과 변형은 각각 98.5%, 97.9% 정도의 좋은 일치를 보여 주었으며 제한조건을 만족하면서 취배수문의 응력과 변형량을 각각 2.2%와 24.8%를 감소시켰다. 본 연구는 취배수문과 같은 대형 구조물의 다중목적 최적설계에 효과적으로 적용될 수 있을 것이다. 해양 구조물은 한번 설치 된 후 교체시까지 약 20~30년을 사용하므로, 본 연구를 통한 응력 및 변형량 감소는 구조물의 사용 수명을 연장에 영향을 줄 수 있을 것으로 기대된다.

References

Matlab : A Practical Introduction to Programing and Problem Elsevier Science & Technology Books, 2009.

Deb, K., Pratap, A., Agarwal, S., Meyarivan, T.. (A Fast and Elitist Multiobjective Genetic Algorithm: NSGA-II, IEEE Transaction On Evolutionary Computation, 2002). 6(2):182-197 10.1109/4235.996017.
crossref
Horn, J., Nafploitis, N., Goldberg, D.. (A Niched Pareto Genetic Algorithm for Multiobjective Optimization Proceeding of the First IEEE Conference on Evolutionary Computation 1994). 82-87.

Hur, D.S., Kim, C.H., Yeom, G.S., Kim, D.S.. (Nonlinear Dynamic Responses among Wave, Submerged Breakwater and Seabed, Journal of Ocean Engineering and Technology, 2005). 19(6):35-43.

Inman, D.J.. (Finite Element Analysis Theory and Application with ANSYS : Thoery and Application with ANSYS 3th ed. Pearson Higher Education, 2007.

Design Standard of Steel Facilities Korea Water Resources Corporation, 2003.

Na, S.Y., Song, J.H., Huh, H., Kim, H.S.. (Design Optimization of Auto-body Member for Crashworthiness Enhancement with the Response Surface Method, Transactions of KSAE, 2005). 13(4):81-89.

Naudascher, E., Rao, P.V., Richter, A., Vargas, P.,Wonik, P.. (Prediction and Control of Downpull on Tunnel Gates, Proc. of the ASCE Journal of Hydr. Engineering, 1986). 112(5):392-416 10.1061/(ASCE)0733-9429(1986)112:5(392).
crossref
Park, J.S., Lee, D.J., Im, B.. (The Study for Construction of the Improved Optimization Algorithm, Proc. of KSAF, 2005). 13(3):22-33.

Weaver, D.S., Martin, W.W.. (Hydraulic model study for the design of the Wreck Cove control gates, Proc. of the Canadian Journal of Civil Engineering, 1980). 7(2):304-314 10.1139/l80-038.
crossref
Westergaard, H.M.. (Water Pressure on Dam During Earthquake, Transactions of ASCE, 1933). 98(2):418-433.



ABOUT
ARTICLE CATEGORY

Browse all articles >

BROWSE ARTICLES
AUTHOR INFORMATION
Editorial Office
President O/T room No.1302, 13, Jungang-daero 180beon-gil, Dong-gu, Busan 48821, Korea
Tel: +82-51-759-0656    Fax: +82-51-759-0656    E-mail: ksoehj@ksoe.or.kr                

Copyright © 2019 by The Korean Society of Ocean Engineers. All rights reserved.

Developed in M2community

Close layer
prev next