항만시설은 외력(파랑, 태풍, 지진에 의한 동수압 등)이 크게 작용하기 때문에 대형 구조물로 이루어진다. 파랑, 태풍 또는 지진등은 실험이 불가능하기 때문에 대부분 전산 프로그램을 이용해서 구조물의 안전성을 평가한다(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)에 적용되었다.
400TEU 이상의 큰 용량의 배가 내항으로 들어 올 경우 내항과 외항의 수위차를 일정하게 유지시키기 위해서 워터펌프(Water pump)을 이용해서 수위를 조절한다. 이 때 취배수문은 해수가 드나드는 통로 중간에 갑문과 함께 설치되는 장비로써, 수위가 달라질 경우에는 해수 통로를 막고 있던 취배수문이 리프팅(Lifting)되어 해수가 드나들게 되고, 시간이 지난 후 수위가 일정해지면, 다시 하강하게 되어 통로를 막게 되는 구조이다. 형상은 설치 조건에 따라서 다르지만 설치형태는 대체로 Fig. 1과 같다.
취배수문에 작용하는 개별 하중은 정수압, 동수압, 지진동수압, 지진하중 등이 있다. 정수압은 작업 시와 유지보수 시 작용하는 수압이 달라지므로 세분화하였다.
2.2.1 개별하중의 종류
① 작업 시 정수압(Hydrostatic force, operating)
작업 시 취배수문을 중심으로 좌우 다른 수압이 작용하기 때문에, 최종적으로 문비에 작용하는 외력은 좌우 수압의 대수적인 차이만큼 작용하게 되며, 수압은 Fig. 2와 같다.
여기서
② 유지보수 시 정수압(Hydrostatic force, maintenance)
수문(Lock-gate)의 유지보수 시 한쪽 면의 해수는 모두 배수 되어야 하므로, 수압은 한쪽 문비(Skin plate)에 대해서만 작용하게 되며, 수압은 Fig. 3 과 같다.
③ 동수압(Hydraulic downpull force)
동수압이 발생하게 되는 원인은 취배수문을 중심으로 좌우 다른 수압이 작용하고 있을 때, 취배수문이 올라가게 되면 해수의 유동에 의해서 위아래 수압차가 발생하고 순간적으로 아래로 당기는 힘이 발생한다. 동수압 발생에 영향을 미치는 인자는 취배수문 하부의 기하학적 형상과 상부를 지나는 유량에 의해서 결정된다(Naudascher et al., 1986; Weaver and Martin, 1980).
여기서 F는 Downpull force가 되며,
[Fig. 5] Resultant downpull and discharge coefficients versus gate opening (Weaver and Martin, 1980)
④ 지진동수압(Seismic hydrodynamics)
지진동수압은 지진에 의한 지반운동에 의해서 해수가 동적 움직임을 발생시켜 취배수문에 수압을 발생하므로 동적응답해석(Dynamic response analysis)방법을 수행하여야 한다. 그러나 동적 응답해석 시 상당한 시간이 소요되므로 본 논문에서는 이러한 방법대신 동적운동을 정적운동으로 근사한 등가정적방법(Equivalent static method)을 사용해서 외력을 구한다. 이러한 근사화 방법은 Westergaard(1933)에 의해서 제안이 되었으며 식은 (4)와 같다.
여기서
여기서
⑤ 지진하중(Seismic inertial force)
지진 하중은 지진에 의해서 구조물이 흔들리는 현상으로 외력은 식 (6)과 같이 구할 수 있다. 여기서
2.1.2 하중 조합
개별 하중들(①~⑤)을 모두 구한 후 각각의 개별 하중을 합리적으로 조합하여 조합 하중을 만들어야 한다. 본 논문은 하중조합을 사용환경의 하중조합에 따라서 적합하게 구성하였으며, 하중 조합을 Table 1과 같이 구성하였다. 하중 조합은 기본적으로 취배수문의 자중(Self-weight)과 정수압이 항상 작용하고 있기 때문에 모든 경우에 대해서 고려하였다. 정수압이 작용하고 있는 상태에서 취배수문이 개방(Opening)되면 동수압이 작용하기 때문에 Case 3과 같이 조합하였으며, 정수압 조건에서 지진이 발생할 경우 지진 동수압 및 지진 관성력이 발생하기 때문에 하중 조합을 Case 4와 Case 5와 같이 구성하였다.
Analysis parameters
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의 최대응력이 항복응력을 초과하지 않았기 때문에 취배수문은 탄성변형영역에 존재한다.
[Table 2] FEM analysis results
FEM analysis results
취배수문에는 여러 종류의 플레이트(Plate)가 사용 되지만, 외력에 대한 응력 및 변형이 발생하는 부재는 Fig. 8과 같이 Skin-plate, Reinforce-plate, Side-plate 세가지로 분류할 수 있다. 여기에서,Sskin-plate는 취배수문을 중심으로 좌우 수압을 유지하도록 해수 통로를 막는 역할을 하며, Reinforce-plate는 수압을 받고 있는 Skin-plate을 지지하게 된다. 그리고 Side-plate는 휠(Wheel)이 부착되어 있는 부재로써 취배수문이 위아래로 이동할 경우 양측면에 Wheel이 부착된 수직통로를 따라서 이동한다.
취배수문의 목적함수 및 제한조건함수는 식 (7)과 같이 정의하였다. 아래의 조건식 중 변형량 허용치는 한국수자원공사의 권고치(Korea water resources corporation, 2003)를 참고로 하였다. 그리고 제한조건 최대 무게는 사용 모터의 리프팅(Lifting) 최대 용량인 16ton으로 하였다.
설계변수들과 목적함수, 제한조건함수의 근사식을 생성하기 위해서 먼저, 중심합성계획법을 통해서 설정된 설계변수(두께)들로 유한요소해석을 실행해서 응력과 변형값을 구하게 된다. 반응표면법에 의해 2차 회귀방정식을 구성하고, 이 2차 회귀방정식을 비지배 분류 유전자 알고리즘(NSGA-II)을 사용하여 목적함수의 최적해를 찾는다. 이 최적해를 유한요소 해석값과 비교하여 정확성을 검증하였다.
4.3.1 중심합성계획법
중심합성계획법은 적은 횟수의 실험으로 곡선을 추정하기 위해서 중심점과 축점을 2
본 논문에서는 취배수문의 두께 3곳이 설계변수이므로 총 실험횟수는
4.3.2 반응표면법
반응표면법은 여러개의 설계 변수가 복합적으로 작용함으로써 어떤 반응변수에 영향을 주고 있을 때 이러한 반응의 변화가 이루는 반응표면에 대한 통계적인 분석방법이다. 본 논문에서는 2차 회귀모형을 사용하였다.
만약, 설계변수(
Table 3과 Table 4을 기반으로 해서 2차 회귀방정식을 만들면 식 (11), (12), (13)과 같다.
[Table 3.] Central composite design
Central composite design
FEM results
Fitness
4.3.3 적합도 평가
중심합성계획법(Na et al., 2005; Park et al., 2005)을 이용하여 반응표면 모델을 구성하였으며, 반응표면 모델의 통계적인 신뢰성을 검증하기 위해서 결정계수
4.3.4 최적화 결과
NSGA-II (Deb et al., 2002)는 빠른 엘리티스트 비지배 분류 유전자 알고리즘(Fast elitist non-dominated sorting genetic algorithm)으로도 불리며, 기존에 사용되던 비지배 분류 알고리즘(Horn et al., 1994)의 단점(계산의 복잡성, 엘리트즘의 부족, 공유 파라미터 지정의필요성)을 개선한 알고리즘으로 본 논문에서는 세대수 300, 인구수 150으로 설정하였다. 반응표면법에 의해 생성된 목적함수 응력과 변형, 그리고 제한조건함수 무게의 근사식은 식 (13)과 같다. 다중목적함수 응력(
[Table 6] Compared with NSGA-II and FEM results
Compared with NSGA-II and FEM results
기존의 취배수문의 기본 두께인 경우와 최적해의 중간값의 두께인 경우를 비교하였으며 Table 7과 같다. 두께의 변화는
[Table 7] Compared with standard and optimization value
Compared with standard and optimization value
본 연구에서는 취배수문의 응력과 변형을 최적화하기 위해서 응력과 변형을 다중목적함수로 설정하고 무게를 제한조건함수로 설정하여 근사최적설계를 하였다. 중심합성계획법을 이용하여 반응표면을 생성하였으며 목적함수, 제한조건함수 모두 99%이상의 높은 적합도를 나타내었다. 반응표면법에 의해 구한 2차 근사모델과 NSGA-II를 이용하여 다중목적함수를 최적화하여 Pareto 해를 얻었으며, 이 Pareto 해 중에서 중간에 해당하는 해를 FEM 해석 결과와 비교하였다. 응력과 변형은 각각 98.5%, 97.9% 정도의 좋은 일치를 보여 주었으며 제한조건을 만족하면서 취배수문의 응력과 변형량을 각각 2.2%와 24.8%를 감소시켰다. 본 연구는 취배수문과 같은 대형 구조물의 다중목적 최적설계에 효과적으로 적용될 수 있을 것이다. 해양 구조물은 한번 설치 된 후 교체시까지 약 20~30년을 사용하므로, 본 연구를 통한 응력 및 변형량 감소는 구조물의 사용 수명을 연장에 영향을 줄 수 있을 것으로 기대된다.
