해저케이블과 파이프라인은 육지에서 도서지역으로 전력, 식수, 통신 등의 수송을 담당하는 사회기반시설이며 설치 위치(천해 또는 심해)에 따라 발생 가능한 손상을 고려한 보호공법이 요구된다. 특히 천해부의 경우 각종 어로 활동이나 선박의 정박활동이 위해요소로 작용한다. 그리고 포설되는 수심의 변화에 따라 위해요소의 종류와 특성이 상이하기 때문에 상황에 맞는 보호 공법의 적용이 매우 중요하다. 천해부에서는 항행하는 선박에 의한 앵커와 어구에 의한 손상을 방지하기 위하여 매설, 보호 구조물 등 다양한 보호공법이 요구된다(Ahn et al., 2007).

기존의 보호공법에 관한 연구를 살펴보면 해저케이블의 매설에 대해 위해요소의 분류 및 대상 지질에 대한 어구와 앵커의 관통깊이 평가가 진행되었다(Allan, 1998). 또한 관련 설계지침으로 보호공법의 위해성 평가, 파이프라인의 규준과 해저케이블의 설치와 보호에 따른 권고사항들이 제시되고 있다(DNV, 2010; DNV, 2012). 그리고 컴퓨터 기술의 발전과 유한요소해석 프로그램의 개발로 인하여 컴퓨터 시뮬레이션을 이용한 연구가 진행되고 있으며 Woo et al.(2009)은 아치형 보호구조물의 앵커 충돌 시뮬레이션을 수행하였다. 또한 위해요소의 충돌력 산정과 관련하여 전산유체해석을 이용하여 앵커의 항력계수를 계산하여 앵커의 수중낙하속도를 산정하는 연구도 진행되었다(Woo et al., 2013).

상기의 보호공법 중 돌을 쌓는 형태인 rock-berm은 구조물의 상부에 돌을 쌓아서 흐름에 의한 유체력과 충돌체에 의해 발생 가능한 외적하중으로부터 해저케이블, 파이프라인을 보호하는 역할을 수행한다. Rock-berm에 대한 기초 연구들은 berm 유형의 구조물의 파력 안정성 해석에 집중되어 있다(Corkum and Martin, 2004; Tørum et al., 2012). 또한 rock-berm의 기본재료가 되는 암석 재료의 수치해석에 관한 연구들도 폭발, 고속 충돌과 같은 고변형률속도를 대상으로 하는 암석의 손상 및 파괴에 관한 연구들이 대부분이다(Hao and Hao, 2013; Zhu et al., 2007). 반면 자유 낙하하는 앵커의 충돌과 같은 저속 충돌에 관한 연구는 Woo et al.(2009)에 의해 연구가 수행된바 있으나 이는 구조물의 충돌해석이고 rock-berm과 같이 작은 입자로 구성된 구조체를 대상으로 하는 연구는 찾아보기 어렵다. 특히 암석으로 이루어진 rock-berm의 충돌해석에 관한 연구는 종전에 수행된 바가 없다. 기존의 Lagrange 기법의 경우 재료내부의 공극을 묘사할 수 없는 단점이 있다. 그와 달리 SPH(Smooth Particle Hydrodynamics) 기법은 무메쉬(meshless)법으로 입자의 이동을 통하여 재료의 거동을 기술하는 방법이며 파편의 비산과정에 대한 수치모사가 실제적으로 재현 가능한 기법이다(Noh and Park, 2007). 이러한 SPH 기법을 이용하는 경우 rock-berm과 같이 입자로 구성된 내부공극이 있는 재료의 충돌해석이 가능하리라 사료된다.

따라서 본 연구에서는 해양에서 발생 가능한 위해요소로부터 해저케이블을 보호하는 rock-berm 구조의 저속 충돌 유한 요소해석 시뮬레이션을 수행하였다. 이때 암석의 모델링에 Smooth Particle Hydrodynamics (SPH) 기법을 이용하여 충돌해석을 수행하고 이를 기존의 Lagrange 기법을 사용한 해석결과와 비교하여 SPH법의 적용 가능성을 알아보았다. 유한요소 해석법의 차이에 따른 동적 유한 요소 해석을 수행하여 응답을 얻었으며 이를 비교하였다. 또한 SPH의 입자 직경변화에 따른 rock-berm의 변위(혹은 충돌) 응답 특성 변화를 알아보고자 하였다.

충돌 해석에 관한 컴퓨터 시뮬레이션의 경우 충돌체와 피충돌체의 충돌에 관한 여러 가지 설계인자가 포함되어야한다. 본 연구에 사용된 ANSYS AUTODYN은 재료의 상태방정식과 외연적(explicit) 해석법을 이용하여 충돌에 의해 발생하는 물체의 대변위, 대변형 문제를 빠르게 해결할 수 있는 유한요소해석 프로그램이다. 일반적으로 Lagrange, Euler와 SPH 기법을 사용한다. 본 연구에서는 앵커와 rock-berm의 충돌해석을 모델링하기 위해 Lagrange, SPH 기법을 사용하였다. AUTODYN solver상의 계산 사이클은 Fig. 1로 간단히 나타낼 수 있다.

암석과 같은 입자들로 구성된 구조체의 거동을 일반적인 요소망을 이용한 유한요소해석으로 구현하기 어렵다. 이는 입자로 구성된 구조체의 거동 해석시 입자간의 상호작용을 무시할 수 없기 때문이다. 본 연구에서는 입자로 구성된 구조체(rock-berm)의 분할에 Smooth Particle Hydrodynamics(SPH) 기법을 이용하여 충돌해석을 수행하고 이를 기존의 Lagrange 기법을 사용한 해석결과와 비교하여 SPH기법의 적용 가능성을 알아보았다.

Lagrange 기반의 유한요소해석 기법은 구조체를 유한개의요소망(mesh)으로 나누는 방식이며 일반적인 구조해석에 많이 이용된다. Lagrange 기법의 경우 구조체를 유한개의 절점(node)을 가지는 요소망으로 분할한다. Fig. 5(a)는 정육면체 구조물을 Lagrange 기반 기법의 요소 분할형태로 나타낸 그림이다. 특히, 공극이 없고 균일한 구조체의 해석에 주로 이용되는 방법으로 모래, 유체와 같은 재료의 해석에 적합하지 못하다는 단점이 있다.

SPH 기법은 요소망이 필요 없는 방법(meshless)이다. 장점은 요소망을 형성하지 않고 요소를 생성할 수 있으며 Lagrange 기법에서 발생하는 대변형에 따른 요소망의 엉킴, 왜곡과 같은 문제점을 해결할 수 있다(Ansys, Inc., 2009). 또한 물체를 유한개의 구형의 입자(particle)로 모델링하기 때문에 유체와 모래와 같이 작은 입자들로 구성된 재료의 해석에 적합한 장점이 있다. Fig. 5(b)는 정육면체 구조물의 SPH 기법의 유한 요소 모델링의 형상을 나타낸다. SPH 기법은 무요소법(meshless method)이라 불리며 kernel estimation을 이용하여 입자간의 분포 특성을 표현한다. 다음 식은 kernel estimation을 이용한 입자의 밀도계산식을 나타낸다.

여기서, ρ^I는 I번째 입자의 밀도, m^J은 J번째 입자의 질량, W^IJ (x^I - x^J, h)는 입자의 직경 h에 따른 가중함수이다. 그리고 상첨자 x^I는 I번째 입자의 x좌표와 x^J는 J번째 입자의 x좌표를 나타낸다. 즉, SPH법의 입자 밀도는 주변 입자의 질량과 입자간의 거리에 따른 가중함수의 곱으로 나타내어진다. Fig. 6은 밀도계산식과 가중함수의 개념도를 나타낸다(Clegg et al., 1997). Kernel B-spline을 가중함수로 사용하며 거리 q에 따른 가중함수의 기본형은 다음 식과 같다.

여기서, N_δ와 h^δ는 차원상수이고 q = |x^I - x^J| / h를 나타낸다(Omang et al., 2005).

SPH 기법은 rock berm의 모델링에만 사용되었다. SPH 기법을 이용한 해석에서 입자의 직경은 중요한 요소 중의 하나로 이는 단위체적당 입자의 수에 따라 결과의 정확성이 달라지기 때문이다. 본 연구에서는 입자의 직경 변화에 따른 응답변화를 알아보기 위하여 설계상의 암석의 평균직경인 200 mm 이외에 50, 100, 150, 300 mm의 직경을 추가로 고려하였다. 입자의 상호작용(interaction)은 고려하였으나 연결(connection)은 고려하지 않았기 때문에 실제 암석의 거동과 유사할 것이라 사료된다.

유한요소해석 기법의 종류와 입자의 직경에 따른 해석 case는 Table 3에 나타내었다. 해석 case에 따라 Lagrange 기법은 대문자 L로 표기하였으며, SPH 기법은 대문자 S에 입자의 직경을 표기하였다. 예를 들어 L은 Lagrange 기법의 해석 case를 나타내고 S200은 SPH 기법의 입자 직경은 200 mm인 해석 case를 의미한다. SPH 기법의 경우 입자의 충돌위치에 따라 결과가 상이하게 나타날 수 있음으로 입자의 상부중앙에 충돌이 발생하도록 모델링을 수행하였다.

Fig. 7의 (a)와 (b)는 Lagrange 기법(L)과 SPH 기법(S200)의 von-Mises 응력 분포도를 나타낸다. 결과를 비교하였을 때 해석 기법의 차이로 인하여 응력 분포도만으로 그 차이를 알아보기 어렵다. 이는 SPH 기법의 경우 입자의 거동에 의하여 각 입자에 발생하는 응력의 방향과 크기가 벡터의 형태로 나타나기 때문이다. 따라서 모든 해석에서 동일한 위치에 5개의 게이지 포인트를 선정하여 응답 결과를 비교하였다. 게이지 포인트는 rock berm의 하부면을 기준으로 상부로 500 mm의 간격으로 선정하였다(0 m, 0.5 m, 1 m, 1.5 m, 2 m). Table 4와 5에 응답결과를 나타내었다.

Fig. 8과 9는 유한요소해석 기법의 변화에 따른 최대 von-Mises 응력과 최대 수직 변위의 분포를 나타낸다. 실제 암석의 평균직경인 S200과 L의 2가지 케이스를 비교하였다. 모든 결과에서 충돌부에서 제일 가까운 2 m 위치에서 최대 응답(von-Mises응력, 최대변위)이 발생한다. 최대 von-Mises 응력은 L에서 132 MPa이 발생하였다. 2 m의 결과를 제외한 나머지 위치에서는 S200의 응력이 크게 발생한다. 최대 수직 변위의 경우 모든 위치에서 S200의 변위가 크게 나타난다. S200의 변위는 위치가 낮아질수록 감소하는 경향을 보이나, L의 결과에서는 0.5와 1 m의 변위가 1.5 m에 비해 크게 발생한다. 또한 응력의 경우에도 0 m에 비해 0.5 m의 응력이 작게 나오는 현상을 보인다. 이는 해석 기법의 차이에 기인한 것으로 Lagrange 기법을 사용한 L의 경우 Fig. 10 (a)와 같이 충돌에 의한 탄성파가 발생하여 일정한 주기로 진동하는 현상을 보이며 SPH기법인 S200은 Fig. 10 (b)와 같이 발생한 변위가 유지되는 결과를 보인다. 일반적인 암석의 경우 공극에 의해 탄성파가 소산되는데 반해 Lagrange 요소법을 사용하는 경우 내부공극의 부재로 인하여 탄성재료의 거동을 나타낸다.

Fig. 11과 12는 SPH 기법의 입자 크기 변화에 따른 최대 von-Mises응력 분포를 나타낸다. 최대 응력값은 S50의 2 m 위치에서 275 MPa가 발생하였다. Fig. 12와 같이 측정위치가 하부로 갈수록 최대응력은 작아지며 입자 크기가 작아질수록 2 m 위치의 최대 von-Mises 응력이 커지는 결과를 보인다. 반대로 0~1.5 m 위치의 최대 von-Mises응력은 입자의 크기가 커질수록 증가하는 양상을 보인다. 이는 SPH 입자의 크기 증가에 따라 gauge point 사이에 존재하는 입자의 수가 감소하여 입자간의 직접적인 응력의 전달이 발생하여 응력이 증가하는 것으로 보인다.

Fig. 13과 14는 입자의 크기 변화에 따른 최대 수직 변위의 분포를 나타낸다. S200에서 1.8 mm의 최대 수직 변위가 발생한다. Fig. 13과 같이 모든 결과에서 최대 변위는 2 m의 위치에서 발생하고 높이가 낮아질수록 변위가 감소하는 경향을 보인다. 2 m 위치의 최대변위는 입자의 크기에 따른 경향성을 찾기 어렵다. 그러나 1.5 m 이하의 변위결과에서 200 mm의 결과를 제외하고 입자의 크기 증가에 따라 최대수직변위도 증가하는 경향을 보인다. 이는 입자의 크기 증가에 인한 것으로 SPH 입자의 크기 증가에 따라 입자간의 변위가 발생하는 것으로 보인다. 상기의 결과로부터 수직 변위가 입자 직경의 변화에 영향을 받는 것을 알 수 있다.

본 연구는 입자로 구성된 구조물의 충돌해석에 관한 연구로 SPH 기법의 충돌해석에 대한 적용성을 알아보기 위하여 rock-berm 구조물의 앵커 충돌 시뮬레이션을 수행하였으며 충돌 응답 결과를 분석하였다. 암석의 모델링에는 Lagrange 기법과 SPH 기법이 사용되었다. 또한 SPH입자의 직경변화(50, 100, 200, 300 mm)에 따른 rock-berm의 응답특성 변화를 분석하였다. 구조응답은 rock-berm에 발생하는 최대 von-Mises 응력과 최대 수직변위로 선정하였다. 이상의 해석으로 얻어진 결론은 다음과 같다.

Rock-berm의 충돌응답(von-Mises 응력, 수직 변위)은 유한요소해석 기법의 변화에 크게 차이를 보였다. Lagrange 기법의 경우 충돌부에 변위가 집중되며 응력은 작게 발생하였으나 응력과 변위가 하부로 갈수록 감소하는 양상은 보이지 않았다. Lagrange 기법의 특성에 의해 rock-berm 하부에 부적절한 응력과 변위의 증가가 발생하였다. 충돌에 의한 탄성파가 발생하여 전파되는 현상이 발생하였으며 Lagrange 기법의 경우 불연속 재료의 충돌응답을 효과적으로 표현하지 못하는 것으로 판단된다.

SPH 기법의 경우 충돌부에 응력이 집중되고 rock-berm의 하부로 갈수록 순차적으로 응력과 변위가 감소하는 결과를 나타내었다. SPH입자의 직경 변화에 대해 충돌응답(von-Mises 응력, 수직변위)은 민감하게 반응한다. SPH입자 직경의 변화에 따라 응답의 분포와 크기가 달라진다.

본 연구로부터 입자로 구성된 재료의 충돌해석에서 SPH 기법이 Lagrange 기법에 비해 사실적인 모사가 가능하다는 것을 알 수 있었다. 그러나 실질적으로 rock-berm에 사용되는 암석의 크기가 다양하므로 이들을 SPH입자로 모사할 경우 입자의 직경을 달리해야하는 단점이 있다. 따라서 향후 암석의 크기를 대표하는 평균직경입자의 개념이 제시될 필요가 있으며 실제 실험을 통한 수치해석의 검증 및 해석결과의 수렴정도의 판단이 요구된다.