Ứng dụng thuật toán lập trình bậc hai tuần tự trong tối ưu hóa mô phỏng chuyển động tàu

CHÚC MỪNG NĂM MỚI 2019 
Tạp chí khoa học Công nghệ Hàng hải Số 57 - 01/2019 45 
ỨNG DỤNG THUẬT TOÁN LẬP TRÌNH BẬC HAI TUẦN TỰ TRONG TỐI ƯU 
HÓA MÔ PHỎNG CHUYỂN ĐỘNG TÀU 
APPLICATION OF SEQUENTIAL QUADRATIC PROGRAMMING ALGORITHM 
FOR SHIP MOTION SIMULATION OPTIMIZATION 
TRẦN KHÁNH TOÀN 
Khoa Công trình, Trường Đại học Hàng hải Việt Nam 
Email liên hệ: toantk.ctt@vimaru.edu.vn 
Tóm tắt 
Trong bài báo này, tác giả giới thiệu qui trình lập trình dựa trên ứng dụng phương pháp tối 
ưu hóa đơn mục tiêu để tối ưu hóa các thành phần chuyển động của tàu từ số liệu thực 
nghiệm của các phép thử tàu điển hình: phép thử quay vòng (Turning Circle), phép thử dích 
dắc (Zigzag). Qui trình này được kết hợp giữa mô hình mô phỏng chuyển động tàu và mô 
hình toán ứng dụng phương pháp tối ưu hóa đơn mục tiêu. Dựa trên thuật toán được đề xuất 
áp dụng là lập trình bậc hai tuần tự đã lập trình mô phỏng và tối ưu hóa để xác định các 
tham số thủy động lực học tối ưu của tàu, từ đó tính toán tối ưu của các thành phần chuyển 
động tàu như vị trí tàu, hướng đi của tàu, vận tốc dọc tàu, vận tốc ngang tàu, 
Từ khóa: Tối ưu hóa đơn mục tiêu, thuật toán SQP, mô phỏng chuyển động tàu, hệ số thủy động 
lực học tàu, phép thử quay vòng, phép thử dích dắc. 
Abtracts 
In this paper, a numerical procedure based on the application of single objective optimization is 
proposed to optimize the ship motion components by using the experimental data of Turning 
Circle and Zigzag tests. This procedure is coupled between the ship motion simulation model 
and single objective optimization mathematical model. Based on the Sequential Quadratic 
Programming algorithm, a numerical simulation and optimization model is carried out to 
determine the optimal ship hydrodynamic coefficients and the optimal of ship motion components 
such as position, heading angle, longitudinal velocity, lateral velocity, etc. 
Keywords: Single objective optimization, SQP algorithms, ship motion simulation, ship hydrodynamic 
coefficient, Turning Circle, Zigzag. 
1. Giới thiệu chung 
Trong lập trình mô phỏng quĩ đạo tàu, thường kết quả mô phỏng có sự sai khác so với kết 
quả đo thực nghiệm, ảnh hưởng tới độ chính xác và độ tin cậy của các hệ thống mô phỏng hàng 
hải. Nguyên nhân của sự sai khác này có thể là do các tác động gây nhiễu của môi trường (sóng, 
gió, dòng chảy, tương tác,..) chưa được tính toán đến trong mô hình toán lập trình mô phỏng, hoặc 
có thể do các hệ số trong hệ phương trình chuyển động tàu chưa đạt đến giá trị tối ưu, Xuất phát 
từ thực tế đó, trên cơ sở áp dụng thuật toán tối ưu hóa đơn mục tiêu, tác giả đã tính toán, lập trình 
mô phỏng chuyển động tàu và xác định các thành phần chuyển động tối ưu của tàu dựa trên nguyên 
tắc cực tiểu hóa giá trị sai lệch về vị trí tàu và các thông số chuyển động tàu giữa dữ liệu mô phỏng 
và dữ liệu thực nghiệm, xác định giá trị tối ưu của các hệ số thủy động lực học trong hệ phương 
trình chuyển động tàu. Kết quả nghiên cứu đạt được dựa trên việc áp dụng các phương pháp nghiên 
cứu lý thuyết, phương pháp lập trình mô phỏng, phương pháp phân tích độ nhạy và phương pháp 
tối ưu hóa. 
2. Mô hình toán chuyển động của tàu 
Mô hình toán chuyển động tàu xét cho tàu chuyển động trên mặt nước với 3 bậc tự do (3 
DOF), gồm 3 thành phần chuyển động: chuyển động thẳng (surge), chuyển động ngang (sway) và 
chuyển động xoay/đảo mũi (yaw) như Hình 1. 
Hình 1. Hệ qui chiếu mô hình hóa chuyển động tàu 
CHÚC MỪNG NĂM MỚI 2019 
46 Tạp chí khoa học Công nghệ Hàng hải Số 57 - 01/2019 
Hệ gồm 3 phương trình chuyển động của tàu ở dạng không thứ nguyên được biểu diễn như sau 
[3, 4]: 
2
" " " "; ; Zv ur gY v ur gY Lk r gLN (1) 
Trong đó: 
u, v, r lần lượt là các thành phần vận tốc chuyển động thẳng, ngang và xoay của tàu; 
�̇�, �̇�, �̇� lần lượt là các thành phần gia tốc chuyển động thẳng, ngang và xoay của tàu; 
" 1 Z
Z
I
k
L m
 là mô men quán tính không thứ nguyên của tàu; 
L là chiều dài tàu, g là gia tốc trọng trường. 
X”, Y”, N” lần lượt là các thành phần lực theo phương dọc, ngang và mô men đảo lái ở dạng 
không thứ nguyên tác dụng lên tàu, được biểu diễn dưới dạng đầy đủ trong các biểu thức (2), (3) và 
(4). Trong đó, 
" " " " " " ", ,..., , ,..., , ,...,u v uv r uvu u c cX X Y Y N N N     là các hệ số thủy động lực học 
(TĐLH) của tàu, đóng vai trò là các biến trong các hàm sai lệch (hàm mục tiêu) thể hiện giá trị trị sai 
lệch về vị trí tàu và các thông số chuyển động tàu giữa số liệu tính toán mô phỏng với số liệu thực 
nghiệm. 
" " " " "
" 2 " "
" " " " 2 2
1 1 1
.
1 1
(1 )
1 1 1 1
u vr vvu u
dc c c c
u vr vvu u
g X X u X u u X vr X v v
L L L
X c c X c c gT t
L L
X u X u u X vr X v
L L L L
 
  
 
   
(2) 
" " " " " "
" " " "
" " "
1 1 1 1
.
1 1 1
1 1
v uv urv v c c
T ur uvc c
v v v c c
g Y Y v Y uv Y v v Y c c Y ur
L L L L
Y c c Y gT Y ur Y uv
L L L
Y u Y v v Y c c
L L

  
   

   
   
(3) 
" " " "
2 2
" " "
2 2
" " " "
2
" "
2
1 1 1
. .
1 1 1
1 1 1
1 1
r uv v r
urc c c c
T ur r uv
vr c c
g L N N r N uv N v r
LL L
N c c N ur N c c
LL L
N gT N ur N r N uv
L L L
N rv N c c
L L
  
  
    
  
  
    
(4) 
Hệ phương trình (1) được giải bởi phương pháp số (phương pháp Euler), sử dụng ngôn ngữ lập 
trình Matlab [5]. Trong nghiên cứu này, tác giả lựa chọn mô hình tàu thực mang tên Esso Bernicia trọng 
tải 193000 DWT, tham khảo dữ liệu thực nghiệm của các phép thử quay vòng tàu (Turning Circle) và 
dích dắc (Zigzag) [2] của Hiệp hội nghiên cứu tàu Anh BSRA (British Ship Research Association) [9]. 
3. Phương pháp tối ưu hóa đơn mục tiêu và thuật toán SQP 
Trong nghiên cứu này, phương pháp tối ưu hóa đơn mục tiêu được áp dụng đối với duy nhất 
một hàm mục tiêu thỏa mãn kèm theo điều kiện ràng buộc bất đối xứng (ở dạng bất đẳng thức) đối 
với các biến tối ưu [6]. Các biến tối ưu (các hệ số TĐLH) trong quá trình tính toán phải đặt trong 
quan hệ toán học với các thành phần lực và mô men thủy động lực học của tàu. Khi đó, bài toán tối 
ưu hóa đơn mục tiêu được định nghĩa như sau: 
CHÚC MỪNG NĂM MỚI 2019 
Tạp chí khoa học Công nghệ Hàng hải Số 57 - 01/2019 47 
Cực tiểu hóa hàm mục tiêu: 
1/2
2
1
N
obj i
i
F f 
 (5) 
Trong đó: 1 2, ,...,
T
N là véc tơ chứa các biến tối ưu (các hệ số TĐLH), N là số lượng 
các biến, fi(α) là các hàm sai lệch của từng thành phần chuyển động tàu giữa giá trị mô phỏng và 
giá trị thực nghiệm. Bài toán tối ưu hóa (5) được giải bằng phương pháp số trên cơ sở áp dụng thuật 
toán tối ưu lập trình bậc hai tuần tự (SQP - Sequential Quadratic Programming). Trong nghiên cứu 
này, tác giả lựa chọn thuật toán SQP vì trong kỹ thuật tối ưu hóa nói chung, thuật toán SQP sử dụng 
phương pháp lập trình bậc hai tuần tự là thuật toán tương đối thông dụng và khả dụng đối với các 
bài toán khoa học kỹ thuật có hàm mục tiêu và điều kiện ràng buộc biểu diễn dưới dạng các hàm 
phi tuyến tính [7]. 
Ngoài ra, vì số lượng các biến lớn (35 biến tối ưu tương ứng với 35 hệ số TĐLH), nên trước khi 
lập trình giải bài toán tối ưu hóa, tác giả áp dụng thêm kỹ thuật phân tích độ nhạy (Sensitive Analysis) 
[8] để phân tích mức độ ảnh hưởng của từng biến tối ưu tới sự biến thiên của hàm mục tiêu. Phân tích 
này đóng vai trò là phép lọc để lựa chọn ra một nhóm các biến có độ nhạy lớn nhất, các biến này có 
ảnh hưởng lớn nhất tới sự biến thiên của các hàm mục tiêu. Phép lọc này có tác dụng rút gọn số lượng 
các biến bằng việc loại bỏ các biến có độ nhạy thấp, chỉ giữ lại các biến có độ nhạy cao, giảm thời gian 
tính toán của máy tính, tập trung sự phân tích ảnh hưởng của các hệ số TĐLH quan trọng ảnh hưởng 
tới sự sai lệch các thông số chuyển động tàu trong lập trình mô phỏng. 
Trên cơ sở đó, tác giả xây dựng biểu thức biểu diễn hàm mục tiêu cho 2 phép thử tàu điển 
hình như sau [13, 14, 15, 16]: 
(1) Phép thử quay vòng: 
𝐹𝑜𝑏𝑗 = (∑ 𝑓𝑖(𝛼)
2
𝑛
𝑖=1
)
1/2
= (∑ ∆𝑆𝑖
2
𝑛
𝑖=1
)
1/2
 (4) 
Trong đó: 𝑓𝑖(𝛼) = ∆𝑆𝑖 là hàm sai lệch khoảng cách giữa quĩ đạo mô phỏng và quĩ đạo thực 
nghiệm tại điểm đo thực nghiệm thứ (i); 𝛼 = [𝛼1, 𝛼2, 𝛼3,  , 𝛼𝑀]
𝑇: là véc tơ của M hệ số thủy động lực 
học của tàu (M = 35) ảnh hưởng đến giá trị độ sai lệch quĩ đạo; N là số điểm đo thực nghiệm. 
(2) Phép thử dích dắc: 
𝐹𝑜𝑏𝑗 = (∑ 𝑓𝑖(𝛼)
2
𝑚
𝑖=1
)
1/2
= (∑ ∆𝜓𝑖
2
𝑚
𝑖=1
)
1/2
 (5) 
Trong đó: 𝑓𝑖(𝛼) = ∆𝜓𝑖 là hàm sai lệch góc xoay giữa giá trị mô phỏng và giá trị thực nghiệm 
tại điểm đo thực nghiệm thứ (i); 𝛼 = [𝛼1, 𝛼2, 𝛼3,  , 𝛼𝑀]
𝑇: là véc tơ của M hệ số thủy động lực học của 
tàu ảnh hưởng đến giá trị độ sai lệch góc xoay; m là số điểm đo thực nghiệm. 
4. Kết quả lập trình mô phỏng và tối ưu hóa 
4.1. Dữ liệu đầu vào 
a) Thông số mô hình tàu thực: 
Các thông số chính của mô hình tàu thực Esso Bernicia 193000 DWT thể hiện trong Bảng 1: 
Bảng 1. Thông số mô hình tàu thực 
STT Thông số tàu Giá trị STT Thông số tàu Giá trị 
1 Chiều dài (Lpp) 304,8 m 6 Tỉ số B/T 2,56 
2 Chiều rộng (B) 47,17 m 7 Hệ số béo (CB) 0,83 
3 Mớn nước (T) 18,46 m 8 Vận tốc tàu thiết kế (U0) 16 knots 
4 Lượng giãn nước (∇) 220.000 m3 9 Vòng quay chân vịt thiết kế (n0) 80 rpm 
5 Tỉ số Lpp/B 6,46 
b) Dữ liệu đầu vào mô phỏng của các phép thử: 
Dữ liệu đầu vào mô phỏng của 2 phép thử Turning Circle và Zigzag thể hiện trong Bảng 2: 
CHÚC MỪNG NĂM MỚI 2019 
48 Tạp chí khoa học Công nghệ Hàng hải Số 57 - 01/2019 
Bảng 2. Dữ liệu đầu vào mô phỏng của 2 phép thử Turning Circle và Zigzag 
STT Thông số đầu vào 
Turning 
Circle 
Zigzag 
1 (x0, y0): vị trí tàu ban đầu (0,0) m (0,0) m 
2 ψ0: góc xoay ban đầu 00 00 
3 U0: Vận tốc tàu ban đầu 5,3 m/s 7,5 m/s 
4 0: góc bánh lái ban đầu 00 00 
5 max : vận tốc góc lớn nhất của bánh lái 2,70/s 2,70/s 
6 n0: vận tốc vòng quay chân vịt ban đầu 57 rpm 80 rpm 
7 nc: lệnh máy (vận tốc vòng quay chân vịt cần đạt tới) 57 rpm 80 rpm 
8 c: lệnh lái (góc quay bánh lái cần đạt tới) -350 deg ±200 
4.2. Kết quả mô phỏng và tối ưu hóa 
a) Phép thử Turning Circle: 
Hình 2. Quĩ đạo ban đầu Hình 3. Quĩ đạo tối ưu Hình 4. Bước lặp tính toán tối ưu 
Kết quả tối ưu hóa quĩ đạo: Quĩ đạo mô phỏng ban đầu thể hiện trong Hình 2, với độ lệch trung 
bình giữa quĩ đạo mô phỏng và quĩ đạo thực nghiệm là 68 m. Quĩ đạo sau khi tối ưu hóa thể hiện trong 
Hình 3, với độ lệch trung bình quĩ đạo đạt được là 5,8 m. Hình 4 biểu diễn sự suy giảm của giá trị hàm 
mục tiêu (qui đổi về giá trị đơn vị) với giá trị ban đầu là 1 và giá trị đạt được sau 30 vòng lặp là 0,084. 
Các thông số tối ưu hóa được tổng hợp trong Bảng 4. 
Bảng 4. Tổng hợp kết quả tối ưu hóa cho phép thử Turning Circle 
Thông số tính toán tối ưu Giá trị 
Sai số giới hạn hàm mục tiêu (giá trị đơn vị) 0,0001 
Sai số giới hạn các biến tối ưu 0,0001 
Số bước lặp tính toán 30 
Giá trị cực tiểu đạt được của hàm mục tiêu 
(giá trị đơn vị) 
0,084 
Độ lệch trung bình quĩ đạo trước tối ưu (m) 68 
Độ lệch trung bình quĩ đạo sau tối ưu (m) 5,8 
b) Phép thử Zigzag: 
Kết quả tối ưu hóa góc xoay: góc xoay mô phỏng ban đầu thể hiện trong Hình 5, với độ lệch 
trung bình giữa góc xoay mô phỏng và góc xoay thực nghiệm là 17,30. Góc xoay sau khi tối ưu hóa 
thể hiện trong Hình 6, với độ lệch trung bình góc xoay đạt được là 6,60, dữ liệu mô phỏng tiệm cận 
dần dữ liệu thực nghiệm từ chu kỳ bẻ lái thứ 2. Hình 7 biểu diễn sự suy giảm của giá trị hàm mục 
tiêu (qui đổi về giá trị đơn vị) với giá trị ban đầu là 1 và giá trị đạt được sau 18 vòng lặp là 0,365. 
Các thông số tối ưu hóa được tổng hợp trong Bảng 5. 
Hình 5. Góc xoay ban đầu Hình 6. Góc xoay tối ưu Hình 7. Bước lặp tính toán tối ưu 
CHÚC MỪNG NĂM MỚI 2019 
Tạp chí khoa học Công nghệ Hàng hải Số 57 - 01/2019 49 
Bảng 5. Tổng hợp kết quả tối ưu hóa cho phép thử Zigzag 
Thông số tính toán tối ưu Giá trị 
Sai số giới hạn hàm mục tiêu (giá trị đơn vị) 0,0001 
Sai số giới hạn các biến tối ưu 0,0001 
Số bước lặp tính toán 18 
Giá trị cực tiểu đạt được của hàm mục tiêu (giá trị đơn vị) 0,365 
Độ lệch trung bình của góc xoay trước tối ưu (độ) 17,3 
Độ lệch trung bình của góc xoay sau tối ưu (m) 6,6 
c) Kết quả lập trình tính toán các hệ số TĐLH tối ưu: 
Bảng 6. Tổng hợp kết quả lập trình tính toán các hệ số TĐLH tối ưu 
STT 
Hệ số 
TĐLH 
Giá trị 
ban đầu 
Giá trị tối 
ưu 
(Turning 
Circle) 
Giá trị tối 
ưu 
(Zigzag) 
Giá trị 
tối ưu 
tổng 
hợp 
STT 
Hệ số 
TĐLH 
Giá trị 
ban 
đầu 
Giá trị tối 
ưu 
(Turning 
Circle) 
Giá trị tối 
ưu 
(Zigzag) 
Giá trị 
tối ưu 
tổng 
hợp 
1 "uX -0,0500 - - -0,050 19 
"
vrN  -0,1200 - - -0,120 
2 
"
vrX 1,0200 - - 
1,0200 20 
"
c c
Y  
0,2080 0,1761 - 0,1761 
3 "vY -0,0200 - - -0,020 21 
"
uvY  0,0000 - - 0,0000 
4 
"
c c
Y   -2,16 - - -2,16 22 
"
uvN  -0,2410 -0,2823 -0,2910 
-0,287 
5 "TY 0,0400 - 0,0300 0,030 23 
"
c c
X  0,1520 0,1684 - 0,1684 
6 "TN -0,0200 -0,0240 -0,0160 -0,020 24 
"
c c
N  -0,0980 -0,0805 -0,0800 -0,080 
7 "rN -0,0728 - -0,0878 -0,087 25 
"
vvX  0,0125 - - 0,0125 
8 
"
v v
Y
 -2,4000 - - -2,400 26 
"
c c
Y   -2,1600 - - -2,160 
9 
"
v r
N
-0,3000 - - -0,300 27 
"
c c
N   0,6880 - - 0,6880 
10 
"
v v
X
 0,3000 - - -1,2050 28 
"
c c
Y   -0,1910 - - -0,191 
11 "uvY -1,2050 - - -0,4510 29 
"
c c
N   0,3440 - - 0,3440 
12 
"
uvN -0,4510 - - 
-0,0500 30 "urY 0,2480 - - 0,2480 
13 "uX  
-0,0500 - - -0,3780 31 
"
urN -0,2070 -0,2105 - -0,210 
14 
"
vY  -0,3780 - - 0,1509 32 
"
u u
X
 -0,0377 - -0,0457 
-0,046 
15 
"
urY  0,1820 0,1598 0,1420 -0,04 33 
"
rN  -0,0045 - -0,0054 -0,005 
16 
"
urN  -0,0470 -0,0416 -0,0380 0,0378 34 
"
u u
X  -0,0061 -0,0073 -0,0073 -0,007 
17 
"
vrX  0,3780 - - -1,500 35 
"
c c
X  
-0,0930 -0,1000 - -0,100 
18 
"
v v
Y  
-1,5000 - - 
Sau khi áp dụng phương pháp tối ưu hóa đơn mục tiêu và lập trình theo thuật toán SQP cho 
2 phép thử Turning Circle và Zigzag, bằng phương pháp phân tích độ nhạy để xác định nhóm các 
hệ số TĐLH có độ nhạy lớn nhất (ảnh hưởng lớn nhất tới sự thay đổi giá trị của hàm mục tiêu trong 
quá trình tối ưu hóa) cho từng phép thử, đã tính toán được các giá trị hệ số TĐLH tối ưu tổng hợp 
chung (giá trị trung bình) từ cả 2 phép thử (Bảng 6). Giá trị tối ưu của các hệ số TĐLH này được sử 
dụng cho các nghiên cứu lập trình mô phỏng theo các ứng dụng cụ thể sau này, như mô phỏng điều 
động, mô phỏng tương tác, mô phỏng cảnh báo đâm va, [10, 11, 16, 17]. 
5. Kết luận 
Xuất phát từ thực tế quá trình xây dựng thuật toán và lập trình mô phỏng của một số hệ thống 
mô phỏng hàng hải trong nước, kết quả dự đoán quĩ đạo tàu thường sai khác nhiều so với kết quả 
thực nghiệm và người lập trình mô phỏng không thể giảm được triệt để sự sai khác này do không 
can thiệp được tới giá trị các hệ số thủy động lực học tàu hoặc không tính toán được đầy đủ các lực 
gây nhiễu của môi trường, tác giả đã nghiên cứu xây dựng mô hình toán và xây dựng chương trình 
mô phỏng chuyển động tàu, trên cơ sở áp dụng phương pháp tối ưu hóa đơn mục tiêu cụ thể là 
thuật toán tối ưu hóa SQP để tính toán tối ưu hóa các hệ số thủy động lực học của tàu và xác định 
quĩ đạo chuyển động tối ưu của tàu. Kết quả tối ưu hóa được kiểm nghiệm dựa trên so sánh số liệu 
CHÚC MỪNG NĂM MỚI 2019 
50 Tạp chí khoa học Công nghệ Hàng hải Số 57 - 01/2019 
mô phỏng với số liệu thực nghiệm thu thập được theo 2 phép thử thực nghiệm quay vòng (Turning 
Circle) và dích dắc (Zigzag) theo tiêu chuẩn của Tổ chức Hàng hải quốc tế IMO. Kết quả nghiên cứu 
có thể được ứng dụng để nâng cấp mô hình toán của các hệ thống mô phỏng đối trong việc mở 
rộng các chủng loại tàu mới không có trong thư viện cơ sở dữ liệu mô phỏng, ngoài ra kết quả nghiên 
cứu có giá trị tham khảo tốt đối với các nhà khoa học, các chuyên gia, giảng viên, học viên cao học, 
nghiên cứu sinh, trong giảng dạy và nghiên cứu khoa học trong lĩnh vực mô phỏng hàng hải. 
TÀI LIỆU THAM KHẢO 
[1] Trần Khánh Toàn, Thuyết minh đề tài cấp Bộ GTVT “Nghiên cứu ứng dụng các thuật toán tối ưu 
hóa để nâng cao độ chính xác dự đoán quĩ đạo tàu của các hệ thống mô phỏng hàng hải”, mã số 
DT154021, nghiệm thu năm 2017. 
[2] International Maritime Organization, Standards for ship manoeuvrability, Resolution MSC 
137(76), 2002. 
[3] Fossen T.I., Guidance and Control of Ocean Vehicles, John Wiley &Sons, 448 pages, 1994. 
[4] Sutulo S., Moreira L. and Soares C.G., Mathematical models for ship path prediction in 
manoeuvring simulation systems, Ocean Engineering, 2002. 
[5] Trần Khánh Toàn, Hoàng Xuân Danh, Lập trình mô phỏng điều động tàu biển cho phép thử 
Turning Circle và Zig-zag theo tiêu chuẩn IMO, Tạp chí Khoa học Công nghệ Hàng hải (ISSN 
1859-316X), số 37 tháng 02/2014. 
[6] Antoniou A. and Lu W.S., Practical Optimization : Algorithms and Enginnering Applications, 
Springer, 202 pages, 2007. 
[7] Zhang J. and Zhang X., A robust SQP method for optimization with inequality constraints, Journal 
of Computational Mathematics, Vol. 21, No. 2, pp. 247-256, 2003. 
[8] Andrea Saltelli, Sensitivity Analysis for Importance Assessment, Risk Analysis, Vol. 22, No. 3, 2002. 
[9] Clarke D., Patterson D.R., Vfooderson R.K., Manoeuvring trials with the 193000 dwt tanker "Esso 
Bernicia"., Paper : Spring Meeting of the Royal Inst, of Naval Architects, No. 10, 1972. 
[10] Trần Khánh Toàn, Mô phỏng tương tác thủy động lực giữa tàu với bờ và đáy luồng nhằm nâng 
cao hiệu quả công tác bảo đảm an toàn giao thông đường thủy nội địa. Tạp chí Giao thông vận 
tải, Số đặc biệt (ISSN 2354-0818), tr.174-176, tháng 11/2015. 
[11] Trần Khánh Toàn, Nghiên cứu tương tác giữa tàu với tàu phục vụ công tác bảo đảm an toàn hàng 
hải trong luồng hẹp. Tạp chí Khoa học Công nghệ Hàng hải, số 49, tháng 01/2017. 
[12] Tran Khanh Toan, Ouahsine A., Naceur H., Hissel F. and Pourplanche A., A fast simulation 
and identification of hydrodynamic parameters for a freely maneuvering ship vessels , 
International Conference on Multiphysics - MULTIPHYSICS 2009, 09-11 December 2009, 
Lille, France, pp.71, 2009. 
[13] Tran Khanh Toan, Ouahsine A., Naceur H., Hissel F. and Pourplanche A., Coefficient Identification for 
Ship Manoeuvring Simulation Model based on Optimization Techniques, International Conference on 
Computational Methods for Coupled Problems in Science and Engineering IV - COUPLED 
PROBLEMS 2011, 20-22 June 2011, Kos, Greece, pp.1261-1272, 2011. 
[14] Tran Khanh Toan, Ouahsine A., Naceur H., Hissel F. and Pourplanche A., Coefficients 
Identification for Ship Manoeuvring Simulation based on Optimization Techniques, International 
Conference on Computational Methods in Marine Engineering IV - MARINE 2011, 28-30 
September 2011, Lisbon, Portugal, pp.369-380, 2011. 
[15] Tran Khanh Toan., Ouahsine A., Naceur H. and El Wasifi K., Assessment of ship manoeuvrability 
by using a coupling between a nonlinear transient manoeuvring model and mathematical 
programming techniques, ELSEVIER/ScienceDirect/Journal of Hydrodynamics/DOI: 
10.1016/S1001-6058(13)60426-6, pp.788-804, ISSN: 1001-6058, 2013. 
[16] Tran Khanh Toan, P.Du, A.Ouahsine, P.Sergent, Simulation of ship manoeuvring in confined 
waterway using a nonlinear model based on optimization techniques, ELSEVIER - OCEAN 
ENGINEERING 142, pp.194-203. ISSN 0029-8018, 2017. 
[17] Tran Khanh Toan, P.Du, A.Ouahsine, P. Sergent, Simulation of the overtaking manoeuver 
between two ships using the non-linear manoeuvring model, SPRINGER - JOURNAL OF 
HYDRODYNAMIC (JHD), pp.791-802, Vol.30, Issue 5, , ISSN: 1001-6058 (Print) 1878-0342 
(Online), Oct. 2018. 
Ngày nhận bài: 22/11/2018 
Ngày nhận bản sửa: 07/01/2019 
Ngày duyệt đăng: 12/01/2019