유체-구조 연성(Fluid-structure interaction, FSI)이란 수치적으로 유체의 유동해석과 고체의 구조해석에 있어 상호작용을 고려하는 것이다. 최근 항공(Park et al., 2011), 조선·해양(Koo et al., 2009; Lee and Rlee, 2011; Shin, 2007,: Shin, 2011) 및 의료(Lee et al., 2011b) 등 다양한 분야에서 FSI 해석에 관한 연구가 진행 중이다. 특히, 조선·해양 분야에서는 선박 또는 해양 구조물이 운행 중 갑판침수(Green water), 슬로싱(Sloshing), 또는 슬래밍(Slamming)과 같은 원인에 의해 유체 충격력을 받게 되면, 선체 구조는 변형이나 응력이 발생하게 되는데, 이때 발생한 응력과 변형은 이들의 안전성과 운항성에 영향을 미치게 된다. 그러나 이전의 수치 해석 기법은 유동해석과 구조해석을 각각 독립적으로 수행해 왔으며 상호작용에 대한 연성을 포함하지 않는 경우가 대부분이었다. 그러나 점차 선박의 대형화 및 해양구조물의 개발이 가속화됨에 따라 정도 높은 해석과 유체-구조 상호연성 해석의 필요성이 강조되고 있다.
유체-구조 상호작용 연성해석이란 유동해석과 구조해석의 연동해석을 의미한다. 연동 해석은 먼저 유동 해석을 통해 구조물의 표면에 작용하는 유체 충격력을 계산하고, 구조 해석시 외력에 대한 하중조건으로 부여하여 구조물의 변형과 응력을 계산한다. 계산된 구조물의 변형 및 변형속도는 다음의 유동 해석에서 역시 경계조건이 되어 참조가 되는 일련의 과정을 포함하게 된다. 일반적으로, 탄성체 모델의 지배방정식은 선형 탄성체의 평형 방정식을, 유동 해석에는 Navier-stokes(N-S) 방정식과 연속 방정식을 사용한다. 이러한 해석기법은 최근 상용 해석 프로그램들을 중심으로 다양한 분야에서 연동해석의 지원을 시작하고 있다(Choi et al., 2009; Kennedy and Solbrekken, 2011).
유동 해석의 경우, 지금까지 연산 속도와 계산기 용량의 한계 등으로 효율적인 수치 연산에 유리한 오일러(Eulerian) 관점의 격자기반 해석 방법들이 주로 개발되어 왔다. 하지만, 최근 컴퓨터의 급속한 발전으로 대용량의 수치 연산이 가능해져 이로 인해 기존의 수치연산에서 배제되어 오던 라그란지(Lagrangian) 관점의 해석 기법인 이동 입자 기반의 시뮬레이션 기법들이 제안되고 있다.
Monaghan(1988)이 제안한 SPH(Smooth particle hydro-dynamics)법과 Koshizuka and Oka(1996)에 의해 제안된 MPS(Moving particle simulation)법은 대표적인 입자를 이용한 유동해석 기법이다. 격자법에 비해 입자법은 형상의 재현이 쉽고 격자 생성이 손쉬우며 특히 물질의 분해 및 결합 등과 같은 극도로 비선형성이 강한 물리적 현상을 표현하는데 큰 장점을 갖고 있다. 그러나 입자법의 경우 비물리적인 압력 진동과 막대한 계산 속도는 해결해야 할 단점으로 지적되고 있다.
본 연구에서는 입자법인 MPS법(Kondo et al., 2005)을 이용하여 탄성체 해석 모델을 개발하고 유체와의 연성 해석에 대한 가능성을 확인하였다. 먼저, 탄성체 모델을 외팔보(Cantilever beam)의 처짐 문제에 적용하여 해석해와의 비교를 통해 개발된 모델의 타당성을 검증하였다. 또한, 유동해석에 MPS법(Koshizuka and Oka, 1996)을 개량한 PNU-MPS법(Lee et al., 2011a)을 이용하여 탄성체 모델과의 연성 모델을 개발하고 외팔보의 끝단에서 자유낙하 하는 정사각형의 물덩이와의 유체-구조 상호작용 문제에 적용해 보았다.
선형 탄성체의 평형 방정식은 다음과 같다.
여기서,
단,
선형 탄성체의 평형 방정식의 지배방정식에 관한 이산화에는 MPS법(Kondo et al., 2005)의 연산모델을 사용하기로 한다.
식 (1)의 첫 번째 항은 다음과 같이 표현 할 수 있다.
여기서 주변 입자
단,
여기서,
단,
식 (6)으로부터 구한 (
여기서,
다음으로 식 (1)의 우변 두 번째 항은 전단변형률에 관한 항으로 다음과 같이 정의 할 수 있다.
식 (10)에서 전단 변형률은 변위벡터의 라플라스(Laplacian) 연산 결과와 같아지며, MPS법의 라플라스 모델을 이용하면 다음과 같이 이산화 할 수 있다.
식 (11)의 변위 벡터 차는 다음과 같이 상대 위치 벡터의 차로서 구할 수 있다.
식 (11)와 같이 상대 위치 벡터에 의한 전단변형률 성분은 초기 상대 위치로 되돌아가려는 성질을 가지게 된다. 그러나 탄성체가 강체 회전을 할 경우 강체 회전에 의한 상대위치 변화성분은 제외되어야 한다. 따라서 강체회전에 의한 성분이 제외된 수정된 상대위치 벡터는 아래와 같다.
여기서
최종 변위 벡터의 가속도 성분은 식 (9)와 식 (14)의 합으로 구성되지만, 여기에는 탄성체의 강체 회전성분이 고려되지 않았다. 따라서 추가적으로 각 운동량 보존법칙을 이용하여 각 입자들의 강체회전성분에 대한 연산을 수행할 필요가 있다. 이때 각운동량 보존식은 다음과 같다.
여기서
여기서
따라서 식 (15)는 입자
이렇게 구해진 각속도 변화량과 속도 변화량을 이용하여 아래의 식들과 같이 다음 시간 단계의 위치와 속도, 각속도, 그리고 회전연산자를 구하게 된다.
여기서, 은 회전연산을 나타낸다.
비압축성 유동장의 지배방정식은 다음과 같은 연속방정식과 N-S 방정식이다.
여기서
식 (22)와 (23)은 PNU-MPS법(Lee et al., 2011a)의 입자간 상호작용 모델을 이용하여 시뮬레이션을 수행한다. 유동해석의 이산화방법은 참고문헌(Lee et al., 2011a)에 상세히 나타나 있다.
Fig. 1에는 유체-구조 상호작용에 관한 개략적인 연성해석법을 나타내며, 다음과 같이 정리할 수 있다.
① PNU-MPS을 이용하여 유동해석을 수행하여 각 벽면상에 작용하는 압력을 계산(이때 유동해석의 벽면경계입자들은 강체로 적용하며, 구조해석으로부터 구한 구조입자의 변형과 변형속도를 이용하여 경계조건으로 부여) ② 유동해석을 통해 구한 각 벽면경계입자들에 작용하는 압력에 각 입자의 미소면적을 곱하여 식 (1)의 우변 두 번째 항인 fi = PiΔAi를 산출하여 구조해석의 외력에 대한 하중조건으로 적용(단, 이때 각 입자의 미소면적은 ③ 구조해석을 수행하여 동일한 시각에서 각 구조입자들의 변형 및 변형속도 계산
이로써 한 시간 단계의 계산이 완료되고 일정한 시간까지 반복하여 연동해석을 수행한다. 이때 시간적분 방법은 시차제법(Staggered method)을 사용하여 매 시간 양방향 연동해석을 수행하였다.
2.1 절에서 설명한 탄성체 해석 기법을 이용하여 2차원 외팔보의 처짐 문제에 적용하였다. Fig. 2는 시뮬레이션의 초기 조건을 나타난다. 외팔보의 길이
Fig. 3는 시뮬레이션에서 구한 외팔보의 처짐량을 식 (24)의 해석해(Ugural and Fenster, 2011)와 비교한 결과이다. 이 때 해석해는 세장비가 작은 외팔보 문제에 대한 해석해로 포아송비에 의한 전단 변형이 고려된 식이다. 전체적으로 보의 처짐에 대한 경향이 잘 일치하는 것을 알 수 있다. 이 때, 최대 처짐에 대한 상대 오차는 3.3% 이내이다. 상대 오차의 발생은 식 (24)가 연속체역학에 기반하여 산출된 결과인 반면, 본 연구에서 적용한 시뮬레이션 기법과 같은 이산화 수치해석법 상의 입자구성 과정 중에 발생하는 상대적 수치 오차에 기인한 것으로 볼수 있다. Fig. 4는 외팔보에 작용하는 응력 분포를 나타내며 식(25)의 해석해와 비교하였다. 여기서
본 절에서는 탄성체 해석 기법을 유동해석 기법인 PNU-MPS법과 연동하여 FSI 해석에 관한 시뮬레이션을 수행하였다. 시뮬레이션의 초기조건은 Fig. 5와 같다. 단, 외팔보의 조건은 3.1절과 동일하다. 외팔보의 끝단에는
Fig. 6은
Fig. 7은 유체의 밀도 차에 의한 외팔보 끝단의 처짐 변위에 대한 시간 변화를 나타낸다. 밀도가 커짐에 따라 전체적인 경향은 다소 유사하지만 최대 처짐량이 증가하며 진동하는 것을 알 수 있다. 이때의 최대 처짐량을 Table 1에 정리하여 비교하였다.
[Table 1] Comparison of maximum deflection the end of cantilever with various water density
Comparison of maximum deflection the end of cantilever with various water density
본 시뮬레이션은 Intel(R) Core(TM) i7 CPU 2.8GHz, 4GHz Memory의 PC에서 수행되었고 0.5초까지의 총 CPU 시간은 2,527초가 소요 되었다.
본 연구에서는 입자법을 이용하여 유체-구조 상호연성 해석을 위한 탄성체 해석 모델을 개발하였다. 개발된 탄성체 모델은 MPS법의 입자간 상호작용 모델(Kondo et al., 2005)을 이용하여 이산화 되었다. 먼저, 모델의 타당성을 검증하기 위하여 단순한 외팔보 끝단에 집중하중이 작용할 경우의 처짐 문제에 적용하였다. 시뮬레이션 결과는 해석해와 정성적·정량적으로 잘 일치함을 알 수 있었다. 다음으로, 유동해석에 PNU-MPS법(Lee et al., 2011a)을 이용하여 탄성체 모델과의 연성 모델을 개발하고 다양한 밀도의 유체 덩어리가 자유낙하 하여 외팔보에 충돌하는 문제에 적용하여 보았다. 이를 통해, 입자법을 이용한 유체-구조 상호작용 해석 기법에 대한 가능성을 확인할 수 있었다. 향후 포아송 비를 포함한 3차원 처짐량에 대한 고려 및 이산화 과정의 고도화를 통해 다양한 유체-구조 해석 문제에 적용하여 개발된 탄성체 모델과 유체해석 기법간의 상호작용 연성 기법의 타당성 검증을 수행할 예정이다.
