Axial vibration analysis is very important for a deep-seabed mining system. In this study, an axial vibration analysis was carried out to estimate the natural frequencies and tensions of the umbilical cable using experimental data obtained from the first pre-pilot mining test. The axial vibrations of the umbilical cable with a pilot mining robot at the bottom end were analytically determined. The range of the added mass coefficients of the pilot mining robot is estimated by comparing the experimental and analytical data. The natural frequencies and maximum tensions are calculated using four estimated added mass coefficients.
한국해양과학기술원 부설 선박해양플랜트연구소에서는 2007년 심해저 망간단괴 채광기술 실용화의 첫 단계로서 자주식 시험집광기를 개발하여, 2009년 우리나라 동해 후포항 부근의 수심 약 100m 해역에서 첫 번째 집광 실증시험(Hong et al., 2009), 2010년 동일 해역에서 경로주행제어시험(Yeu et al., 2010)을 각각 수행하였다.
2차례의 근해역 실증 시험을 통해 획득된 시험집광기 기술을 바탕으로, 2011년부터 5,000m급 파일럿 규모(상업 채광용광의 1/5)의 채광로봇(미내로, MineRo)을 설계하여, 2012년 9월 제작을 완료하였다. 제작된 파일럿 채광로봇(Fig. 1)은 동일한 구조를 가진 2개의 단위모듈로 구성되는 것이 특징이다. 파일럿 채광로봇의 예비성능 검증을 위해, 선박해양플랜트연구소가 보유하고 있는 심해저 집광실험동에서 단위모듈별 전력 · 신호 통전 시험(Electric-electronic system tset), 유압구동시험(Hydraulic system test), 주행시험(Driving test), 파쇄시험(Crushing test), 채집시험(Collecting test) 등을 각각 수행하였다. 2012년 11월, 동해안 수심 130m 해역에서 1차 예비 파일럿 집광성능시험(Pre-pilot mining test, PPMT)을 수행하였다(Hong et al., 2013). 1차 성능시험은 파일럿 채광로봇의 기계적 성능과 모조단괴(Artifical nodule)를 이용한 채집성능을 시험하는 것에 초점을 두었다
예비 파일럿 시험에서의 채광로봇은 와이어로 보호된 엄빌리컬 케이블에 연결되어 진 · 회수가 된다. 엄빌리컬 케이블은 선박 갑판에 설치된 윈치에 연결되어 선박의 상하운동에 의해 축방향 운동을 하게 된다. 이러한 동적 운동은 엄빌리컬 케이블에 작용하는 장력의 변화를 발생시킬 수 있다. 특히, 외력에 의한 가진 주파수가 엄빌리컬 케이블의 고유주파수에 근접할 경우, 엄빌리컬 케이블에 작용하는 장력은 급격히 상승하게 된다. 따라서 파일럿 채광로봇의 진 · 회수 작업에서의 안정성 검토를 위해 케이블의 동적해석은 매우 중요하다.
물체가 매달린 해양 케이블의 동적해석에 관한 많은 연구가 수행되어왔다. Berteaux(1975)는 유체 점성력을 무시하고 케이블 끝단에 물체가 매달린 경우 케이블 상부와 끝단에서 발생되는 장력을 계산하였다. Jung et al.(2008)은 케이블의 장력을 수치해석 모델을 이용하여 계산하였다. Namba et al.(2009)은 라이저의 안쪽과 바깥쪽의 유체감쇠력의 상호작용을 고려하여 라이저의 장력을 계산하였으며, Boe and Nestegard(2010)은 케이블 하단에 물체가 매달려 진수되는 경우의 축방향 하중을 계산하였다. Park et al.(2010)은 유체 점성력과 구조감쇠력을 고려하여 해양 파이프의 축방향 진동해석을 수행하였다. Nam et al.(2012)와 Nam et al.(2013)은 상하동요 저감장치(Heave compensator)를 고려한 해상 크레인의 시간영역 연성 운동해석을 수행하였다. 하지만, 위의 많은 연구에도 불구하고 물체의 부가질량계수 추정에 관한 연구는 매우 제한적으로 수행되었다.
본 논문에서는 보다 정확한 축진동 해석을 위해 파일럿 채광로봇 부가질량 계수 추정에 관한 연구를 수행하였다. 파일럿 채광로봇에서 계측된 축방향 가속도 데이터를 분석하여, 파일럿 채광로봇의 부가질량 계수의 범위를 추정하였다. 추정된 부가질량 계수를 적용하여, 1차 예비 채광성능시험에서의 엄빌리컬 케이블-파일럿 채광로봇의 축방향 고유주기와 최대 인장력을 계산하였다. 끝으로 계측된 데이터를 분석하여 채광시스템의 안정성을 검토하였다.
1차 예비 집광성능시험에서 파일럿 채광로봇의 진수 및 회수 과정에서의 축방향 가속도 데이터는 파일럿 채광로봇 몸체에 부착된 RLG(Ring laser gyro) 센서(MiniRLG2, CDL)를 통해 계측된다. RLG 센서의 사양을 Table 1에 정리하였다.
[Table 1] Specification of RLG
Specification of RLG
Fig. 2와 Fig. 3은 2012년 11월 6일, 로봇의 진수 과정에서 계측된 수심 데이터와 수직 방향의 가속도 데이터이다. Fig. 3에서 1번은 로봇이 해수면으로부터 진수가 시작되는 부분이며, 2번과 3번은 윈치에 의해 케이블이 점차 풀려져나가면서 로봇이 점차 깊이 진수되는 부분이다. 4번은 로봇이 해저면에 안착되어 진수가 끝난 부분이다. 2번 영역에서 계측된 값이 다른 영역의 값보다 큰 변화폭을 갖는 이유는 진수와 동시에 트랙이 구동되어 로봇에 추가적 진동을 발생시킨 것으로 분석된다.
5번의 사각 점선으로 표시된 영역은 진수 과정에서 로봇이 일시적으로 정지된 지점으로, 금번 시험에서의 부가질량 계수 및 축방향 고유주기를 계산할 수 있는 가장 적합한 영역이다. 여기서 계측된 가속도를 Fig. 4에 자세히 나타내었다.
Fig. 5와 같이 길이가
여기서
여기서
식 (3)을 식 (1)에 대입하면 아래와 같다.
식 (4)를 변수분리법으로 풀이하면
여기서
케이블 하단(
여기서,
여기서
식 (8) ~ (10)을 식 (6)에 대입하면 식 (11)을 얻을 수 있다.
식 (11)을
식 (12)를 식 (5)에 대입하면,
케이블 상단(
여기서,
식 (14)를 고려하여 식 (13)을 풀이하면 케이블의 변위 함수를 다음과 같이 얻을 수 있다.
케이블 하단에서의 변위 함수와 가속도 함수는 다음과 같다.
또한, 케이블 상단에서의 변위 함수와 가속도 함수는 다음과 같다.
케이블 하단과 상단에 작용하는 인장하중은 아래의 식으로 각각 정의될 수 있다.
여기서,
2.3 1차 예비 집광성능 엄빌리컬 케이블 축진동 해석
1차 예비 집광성능에 사용된 엄빌리컬 케이블의 물성치와 형상정보를 Table 2에 나타내었다. 파일럿 채광로봇의 수중무게는 9381 kg이며, 크기는 폭 4m, 길이 6.5m, 높이 5m 이다. 케이블의 부가질량 계수(
[Table 2] Properties and geometric information of the umbilical cable
Properties and geometric information of the umbilical cable
파일럿 채광로봇의 부가질량 계수(
Fig. 7 ~ Fig. 9에서는 계측된 가속도와 계산된 가속도의 크기가 유사한 반면, Fig. 10 ~ Fig. 12에서는 계산된 가속도의 크기가 계측된 가속도 보다 매우 크게 나타나고 있다. 이러한 데이터 비교를 통해
[Table 3] Resonance wave periods of axial vibration with different added mass coefficient
Resonance wave periods of axial vibration with different added mass coefficient
1차 예비 집광성능시험의 엄빌리컬 케이블의 축진동 공진주기는 파일럿 채광로봇의 부가질량 계수 변화에 따라 1.1초 ~ 3.1초로 계산되었다. 실해역 실험 시, 파주기 6.3초에 비해 케이블의 축진동 공진주기가 매우 작으므로 안정된 영역에서 1차 예비 집광성능실험이 수행되었음을 알 수 있다. 파주기가 6.3초인 경우 엄빌리컬 케이블 상부에 작용하는 최대 인장하중은 부가질량 계수 0, 0.5, 1, 1.5에 따라 113.16kN, 143.44kN, 178.88kN, 220.93kN으로 각각 계산되었다.
본 논문에서는 1차 예비 집광성능시험에서 계측된 엄빌리컬 케이블의 축방향 가속도 응답을 분석하였으며, 고유주기와 최대인장하중을 계산하였다. 엄빌리컬 케이블의 축진동 분석을 통해 다음과 같은 결론을 얻을 수 있었다.
(1) 계측된 가속도 데이터를 이용하여 파일럿 채광로봇의 부가질량 계수를 0.5 ~ 1.5 범위로 산정하였다.
(2) 파일럿 채광로봇의 부가질량 계수가 커짐에 따라 축방향 공진주기는 1.1초 ~ 3.1초로 변화되었다.
(3) 실해역 시험 시 발생된 해양파의 주기 6.3초는 공진영역에서 매우 떨어진 영역이므로, 1차 예비 집광성능시험에서 파일럿 채광로봇의 진/회수 시 케이블의 축진동 공진이 발생할 가능성이 매우 낮았음을 알 수 있었다.
(4) 엄빌리컬 케이블의 상부에 작용하는 축방향 인장하중은 최대 220.93kN으로 실험에 사용된 A-프레임의 안전운용하중 255kN의 범위로 계산되었다.
향후 보다 정확한 파일럿 채광로봇의 부가질량 계수 추정을 위하여 케이블 상부 변위의 실측이 필요하며, 이를 고려한 케이블 해석 기법의 연구가 수행되어야 한다.
