과거 해군의 전투력은 인간의 노동력이나 화약무기를 바탕으로 구성되어 왔으나 과학기술의 발달에 따라 미사일기반의 무기체계를 갖추면서 전자 장비를 이용한 무기체계로 변화되어왔다[1].

이러한 변화는 최첨단 탐지체계 및 고성능 전자 장비 등의 무기와 센서를 지원할 수 있는 전력 체계의 발달을 야기한다. 미 해군의 차세대 해군 함정 로드맵에 따르면 기계식 추진체계기반에서 복합식 추진체계 이후, 통합전력시스템(IPS, Integrated Power System)을 적용할 것으로 보인다. 이는 함정의 소요전력 증가에 따라 그 필요성이 높아지고 있는 것으로 보인다[2][3].

최근 전력변환사업 및 기술의 발달에 따라 추진전동기(EPM, Electric Propulsion Motor)를 이용한 복합식 추진체계 및 전기식 추진체계 추진방식이 적용되어지고 있다[4][5]. 이러한 추진전동기를 적용한 추진체계를 사용함으로써 기존의 기계식 추진체계에 비하여 가지는 장점으로 수중방사소음의 절감으로 인한 대잠능력, 생존력 강화 및 전력운영의 유연함 등이 있다. 또한 함정의 경우 작전 및 훈련 상황에 따라 25노트 이상의 높은 속력을 요구하지만 고속으로 운항하는 시간 및 빈도가 매우 낮은반면 10~20노트의 저속에서 운항하는 빈도는 매우 높다[6][7]. 내연기관 엔진 및 가스터빈만을 추진장비로 사용하는 단일 기계식 추진체계의 경우 저속에서 효율이 저조하지만 전기추진모터를 적용한 하이브리드 추진체계에서는 모든 속력구간에서 높은 효율로 운항이 가능하다.

그러나 이러한 하이브리드 추진시스템은 기계식 추진시스템에 비해 설계 및 생산이 복잡하고 연동 장비의 수가 증가함에 따라 연동 공정에서 발행할 수 있는 문제점이 많아 추진시스템에 대한 검증이 필수적이다. 특히 통합기관제어장치인 ECS(Engineering Control System)는 함정 운용에 매우 중요한 장비로 높은 신뢰성이 요구된다. 이러한 이유로 현재 하이브리드 추진시스템이 적용된 함정에 적용되는 통합기관제어장치는 L3, SEASTEMA와 같은 외국 전문 회사에 의지하고 있는 실정이다. 이러한 주요 기술의 국외 의존도는 국내 함정기술의 발달을 저하시키고 있으며 나아가 국내 함정 기밀정보정보의 해외유출로 인한 보안 문제도 함께 발생할 수 있다.

본 연구에서는 cRIO 및 LabVIEW를 활용하여 주요 추진장비 Modeling을 통한 CODLAD(COmbined Diessel-eLectric Or Gas turbine) 추진체계의 ECS 검증 및 함정의 건조 이전에 발생할 수 있는 문제점을 사전에 분석 및 검증하기 위한 RT-HILS(Real-Time Hardware In the Loop System) 기반의 시뮬레이터를 개발한다.

함정의 추진체계는 추진기관의 조합에 따라 다양한 형태를 갖춘다. 디젤엔진이나 가스터빈과 같이 기계식 추진 장비 중 하나의 추진 장비를 적용한 추진체계를 단일 기계식 추진체계라 하며, 두 개 이상의 기계식 추진체계 또는 기계식 추진체계와 전기모터 등의 전기식 추진체계를 혼용하여 적용한 것을 복합식 추진체계 그리고 추진 장비로 전기모터를 사용하며 복수의 발전기 또는 배터리 등을 이용한 추진체계를 완전전기추진체계 또는 통합전력시스템이라 한다. 복합식 추진체계는 추진 장비의 구성에 따라 아래 Table 1과 같이 분류할 수 있다.

과거 해상 전투에 있어 높은 기동성은 함정의 생존성을 높여주었으나 현대의 전투에서는 최첨단 탐지체계를 활용하여 전투를 수행하므로 적에게 발각 위험을 낮추는 수중 방사소음 저감이 함정의 속도만큼 중요하게 되었다[2]. 또한, 첨단 장비의 탑재는 함정의 소요 전력량을 높이게 됨에 따라 여유전력을 확보할 필요가 있어 추진전동기를 탑재한 복합식 추진체계의 필요성은 더욱 높아지고 있다.

본 논문에서 개발하고자 하는 추진체계 장비 모델 및 시뮬레이터는 대구급 호휘함인 FFX(Future Frigate eXperimental)함을 CODLOG 추진체계기반의 가상함정을 대상으로 한다. 대상함정은 가스터빈 1기, 추진전동기 2기, 클러치 2기, 가변피치프로펄러 2기로 구성되어 Figure 1과 같이 구성되며 주요 장비 사양은 Table 2와 같다.

Figure 1: 
Configuration of propulsion system (CODLOG)

함정 및 함정 추진체계의 모든 동특성을 모델링하기 위해서는 모든 장비 및 특성에 대한 자료가 확보되어야 한다. 이러한 주요 장비의 자료는 제조사의 보안상 문제로 확보의 어려움이 있다. 따라서 본 논문에서는 시스템 동정(SI, System Identification)을 이용하여 비교적 적은 데이터로 시뮬레이션 모델을 구현하고자 한다. 시스템 동정은 모델 종류에 따라 White box model, Gray box model, Black box model로 나눌 수 있으며, 본 연구에서는 Figure 2와 같이 black box model을 적용하여 구현한다.

Figure 2: 
Black box model (system identification)

시스템 동정의 black box model은 입력(u(k)) 데이터와 출력(y(k)) 데이터의 응답특성으로 시간 영역에 대한 모델을 알아낼 수 있다.

3.1 선속 모델

함정의 속력은 대상 함정인 FFX-Batch III의 ROC(Requirement Of Concept)에 따라 35,000 kW에서 최대 속도 30 knot이상의 속력를 가질 수 있도록 모델을 설계한다. 함속과 출력의 관계는 Froude 상사법칙과 Reynold 상사법칙을 적용하여 식 (1) 및 식 (2)와 같이 나타낼 수 있으며, 함정의 저항성분 중 조파저항과 와류저항 및 공기저항은 무시하였다.

Rf=12×f×S×V-1.825

(1)

IHP∝RfV∝D23V

(2)

여기서 R_f는 마찰저항, f는 실험계수, S는 함정의 침수 표면적, V는 함속, IHP는 지시마력, D는 함정의 배수량을 나타낸다. 즉, 함정의 마찰저항은 선속의 거듭제곱에 비례하며 침수표면적의 경우 배수량의 2/3승에 비례한다.

Figure 3은 대상 모델선박의 속력에 따른 필요 출력을 나타낸다.

Figure 3: 
Ship speed-power curve

3.2 가스터빈 모델

본 논문에서는 Rolls-royce사의 MT30을 대상으로 가스터빈 모델을 시뮬레이션 하였다. MT30 가스터빈은 최대 3,600 rpm에서 약 40 MW의 출력을 내며 군함 및 항공기 등 다양하게 적용된다[8], 또한, MT30 가스터빈의 출력은 본 논문에 적용하고자 하는 함정의 최대속력에 필요한 출력을 갖추므로 적합하다. Table 3는 MT30 가스터빈의 특성을 나타내며, Figure 4은 가스터빈의 출력 관계를 나타낸다.

Table 3: 
MT30 gas turbine model factors and coefficients

Input factors	Value
ECS power demand	0~100 [%]
Internal power demand	0~40 [MW]
Feedback power	0~40 [MW]
Feedback GT rpm	0~3,500 [rpm]
CPP P/D	0~100 [%]
Output factors	Value
External power	0~40 [MW]
Shaft torque	0~500 [kNm]
GT RPM	0~3,500 [rpm]
Gas turbine model coefficients
Ship’s torque coefficient according to pitch
Ship’s torque coefficient according to RPM
Ship’s torque coefficient according to ship speed

Figure 4: 
Gas turbine power diagram

가스터빈 모델은 ECS로부터 내려오는 TIC(Throttle Input Command)입력신호에 따라 가스터빈 출력을 변환시킨다. 이러한 가스터빈 출력은 함정 모델 및 가변피치프로펠러의 P/D(Pitch/Diagram ratio)값에 따라 축토크 및 축회전수를 도출한다. 가스터빈의 출력-축속도-토크의 관계는 식 (3)과같이 나타낼 수 있다.

PW=TNm×w[rad/sec]=TNm×nrpm×2π60

(3)

식 (3)을 기본수식으로 하여 가스터빈 TIC 값에 따른 지연 PID제어기를 설계하여 토크 및 축회전수, 출력을 나타내는 가스터빈 모델을 적용하였다. 토크[T]는 가스터빈에서 발생하는 토크와 프로펠러 피치 각에 의해 발생되는 토크를 고려하여 수학적 모델을 구현하였으며 가스터빈과 축의 기어비는 15.54를 적용하였다. Figure 5는 가스터빈의 블록 다이어그램을 나타낸다.

Figure 5: 
Gas turbine model block diagram

3.3 추진전동기 모델

대상 함정인 FFX Batch-II 함정모델의 경우 2축 2기의 추진전동기를 탑재하는 모델을 구현하여야 하나, 본 논문에서는 시뮬레이션의 단순화를 위해 1대의 추진전동기에 대하여 모델을 구현하되 두 배의 설계출력을 갖는 추진전동기를 모델을 적용하였다. Figure 6과 Figure 7은 추진전동기 모델 추력 및 블록 다이어그램을 나타낸다.

Table 4: 
Propulsion motor specification

Item	Specification
Type	3,592 ton
Rating Power	2,800 ton

Figure 6: 
EPM power diagram

Figure 7: 
EPM model block diagram

3.4 가변피치 프로펠러

프로펠러는 추진기관으로부터 나오는 출력을 축계 및 함체에 전달하여 추력을 얻도록 한다. 프로펠러의 출력(P_P[kW]) 즉, 추력(T_P[N]) 및 토크(Q_P[Nm])는 프로펠러의 회전속도(N_P[rps])와 프로펠러의 직경(D[m]) 및 해수의 밀도(ρ[kg/m³])로 식 (4), 식 (5), 식 (6)과 같이 나타낼 수 있다[9]-[11].

PP=2πQP=2πKQρD5Np2

(4)

Tp=KTρD4Np2

(5)

QP=KQρD5NP2

(6)

이러한 프로펠러의 특성은 아래 매개 변수에 의해 결정된다. 여기서 V_A는 함정의 전진속력이다. 가변피치프로펠러의 경우 차원계수 K_P와 K_Q는 전진계수 J와 프로펠러 블레이드의 방향각 P에 대한 함수로 표현할 수 있다.

Forward coefficient J=VA/(NPD)

(7)

Trust coefficient KT=TP/(ρNP2D4)

(8)

Torque coefficient KQ=QP/ρNP2D5

(9)

Propeller efficiency no=VATP/2πNPQP=J/2πKT/KQ

(10)

KQ=KQJ,P

(11)

KT=KTJ,P

(12)

여기서 P는 기준피치(P₀)에서 시작하는 프로펠러 블레이드의 방향각을 의미하므로 고정피치프로펠러의 경우 P가 고정되어있으므로 J를 다항식으로 표현할 수 있으나, 가변피치프로펠러의 경우 P의 변화에 따라 J와 프로펠러의 회전속도에 따라 지속적으로 변화하게 된다. 따라서 본 논문에서는 가변피치프로펠러의 특성을 Table 5 및 Figure 8과 같이 가정하여 시뮬레이션을 진행한다.

Table 5: 
Controllable pitch propeller specification

Item	Specification
Type	CPP
Number of blade	5
Propeller diameter	4.3 m
Pitch to diameter ration	1.8
Area ration	0.57

Figure 8: 
CPP characteristic curve

Figure 9는 토크와 속도, P/D를 출력하는 가변피치 프로펠러의 블록다이어그램이다.

Figure 9: 
CPP model block diagram

3.5 감속기어 및 클러치

감소기어는 주 추진기관의 동력전달 및 동력전환 등의 기능을 수행한다. 또한 함정의 경우 프로펠러의 회전수가 낮을 때 효율이 높으므로 중속 및 고속기관의 속력를 낮추고 토크를 높이는 역할을 수행한다.

최근 개발되는 감속기어의 경우 저속구간을 제외한 대부분의 구간에서 98%이상의 에너지 전달효율을 가지므로 본 논문에서는 전달효율을 100%로 가정하여 모델을 적용한다. 또한 하이브리드 추진기관이 적용됨에 따라 함정 추진 장비의 안정적인 절환을 위해 SSS-Clutch 동작절차를 모델링하여 적용한다.

Figure 10은 SSS-clutch모델의 운전 모드별 블록다이어그램을 나타낸다.

Figure 10: 
SSS-clutch model block diagram

본 논문에서는 함정 추진체계의 검증 및 ECS 검증을 위한 가상의 함정 시뮬레이션을 진행한다. 시뮬레이션은 함정의 추진 장비모델과 추진체계를 가시적으로 보여줄 수 있는 하드웨어 유닛과 가상의 ECS로 Figure 11과 같이 구성한다.

Figure 11: 
Propulsion system simulator layout

시뮬레이션 및 시뮬레잍의 제작은 NI사의 LabVIEW와 Real-Time기반의 OS를 활용한 CRIO를 활용하여 진행하였다. 시뮬레이터 사양은 Table 6과 같다.

시뮬레이션의 제어명령을 받는 각 추진 장비 및 함체 모델은 함정의 축회전수, 토크, 함속 등을 나타내게 되며 신호 및 추력의 전달체계는 Figure 12과 같다.

Figure 12: 
Propulsion

본 논문에서는 개발한 추진체계 장비 모델 및 시뮬레이터의 성능을 시험하기 위해 시뮬레이션을 통한 장비 특성과 시뮬레이터의 특성을 비교 및 분석한다. 실제 함정의 동특성을 분석하기 위해서는 열역학적, 동역학적 모델이 필요한다. 그러나 모델의 구현을 위해 많은 변수들에 대한 정확한 값을 알아야 하며 시뮬레이션에 많은 시간이 필요하여 실시간 시뮬레이션에 적용하기 어렵다는 단점을 가진다. 따라서 본 논문에서는 시스템 동정기법(SI)의 Black box 모델을 활용하여 실제 추진 장비의 입력 및 출력 값을 바탕으로 전달함수의 신뢰성을 높이고 모델을 단순화하여 실시간 운영체제에 적합하도록 개발하였다.

개발한 시뮬레이션은 가스터빈 기본모델, 가스터빈 내부제어기 모델, 추진전동기 기본모델, 추진전동기 내부 제어기 모델, 가변피치프로펠러 모델, 기관통합제어장치(ECS) 모델, 데이터 측정 및 저장 프로그램으로 구성되며, 이를 하드웨어적으로 표현하기 위해 전동기 2대, 전자식 클러치 1대, 베벨기어 2대, 포토센서, 베어링 등으로 Figure 11과 같이 구성하였다. 시뮬레이션 루프 속도는 실제 함정시간의 1/10 scale로 진행하여 사용자가 보다 빨리 시뮬레이션 결과를 확인할 수 있도록 구현하였다.

해외의 경우 함정의 선도함 건조 및 LBTS (Land Based Test Site) 등을 이용하여 함정의 추진체계를 검증하거나 HILS(Hardware In the Loop System)를 이용하여 함정 추진체계를 검증하고 있으나, 국내의 경우 예산 및 기술의 해외 의존으로 인하여 복합식 추진체계의 검증에 어려움을 겪고 있다.

본 논문에서는 복합식 추진 함정의 추진체계검증을 위하여 추진 장비에 대한 물리적 특성 값 확보가 어려움에 따라 시스템 동정을 이용한 추진 장비 모델을 구현하여 추진체계 및 통합기관제어장치의 특성을 검증 할 수 있도록 구현하였다.

Figure 13과 Figure 17는 실제 장비를 바탕으로 SI기법을 활용한 시뮬레이션 결과와 하드웨어 기반의 시뮬레이터 모델의 결과를 나타낸다. 실제 장비의 데이터의 값과 시뮬레이션 결과 값을 비교하여 가스터빈과 추진전동기의 실제 데이터와 추진 장비 모델 간의 오차는 각각 3.04 %와 3.26 %임을 확인하였으며, 하드웨어 기반 시뮬레이터 모델과 시뮬레이션 결과와의 오차는 가스터빈 2.97 %, 추진전동기 3.38 %의 평균 오차율로 실제 장비 데이터와 유사한 결과 값을 가짐을 확인하였다. 본 시뮬레이터는 해군 함정의 설계 단계에서 추진 성능을 예측할 수 있도록 하여 함정 추진체계의 분석 및 최적화와 통합기관제어장치 제어기의 성능 및 신뢰도 검증 등의 핵심기술을 확보할 수 있는 기반기술로 활용할 수 있다. 추후 추진 장비 모델의 정밀도 향상을 위해 추진 장비 모델의 최적화 및 함정 선형 모델연동, 전력 계통 모델과의 연동에 대한 연구가 필요할 것으로 사료된다.