Scribd
Upload
182 views5 pages

Model 01

1) The document describes mathematical models for marine vehicle motion, including equations of motion for surge, sway and yaw for a Mariner Class vessel. 2) Key parameters and dimensions are provided for the Mariner Class vessel, along with non-dimensional coefficients for the hydrodynamic forces and moments. 3) A MATLAB function is presented that calculates the time derivative of the vessel state vector based on the vessel model and input rudder angle.

Uploaded by

hoangductuan
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
182 views5 pages

Model 01

1) The document describes mathematical models for marine vehicle motion, including equations of motion for surge, sway and yaw for a Mariner Class vessel. 2) Key parameters and dimensions are provided for the Mariner Class vessel, along with non-dimensional coefficients for the hydrodynamic forces and moments. 3) A MATLAB function is presented that calculates the time derivative of the vessel state vector based on the vessel model and input rudder angle.

Uploaded by

hoangductuan
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 5

Extracted from Fossen (1994), Guidance and Control of Ocean Vehicles, Wiley and Sons Inc.

Mathematical Models of Marine Vehicles In order to verify a good control design it is useful to simulate the control law against a realistic model of the vessel. The following motion parameters will be used to describe the different mathematical models:

u = u 0 + u v = v 0 + v

p = p 0 + p r = r0 + r

= 0 + = 0 + (A3.5.1)

For instance, the definition implies that u is a small perturbation from a nominal (constant) surge velocity u0 while u denotes the total surge velocity. The total speed of the vessel is defined according to:
U = u 2 + v2 =

(u 0 + u )2 + v 2

(A3.5.2)

1. Mariner Class Vessel

The hydro- and aerodynamics laboratory in Lyngby, Denmark, has performed both planar motion mechanism (PMM) tests and full-scale steering and manoeuvring predictions for a Mariner Class Vessel. The main data and dimension of the Mariner Class Vessel are (Fossen, 1994):
Mariner Class Vessel Length overall (Loa) Length between perpendiculars (Lpp) Maximum beam (B) Design draft (T) Design displacement ( ) Design speed (u0) 171.80 (m) 160.93 (m) 23.17 (m) 8.23 (m) 18541.00 (m3) 15.00 (knots)

For this vessel the dynamic equations of motion in surge, sway and yaw are: 0 0 m X u u X 0 m Yv mx Yr v = Y G x N v z N r r N 0 m G I

(A3.5.3)

where the nonlinear forces and moment X , Y and N are defined as (Prime System I with Lpp and U as normalisation variables, see attached document):

X = Xu u + Xuu u 2 + Xuuu u 3 + Xvv v 2 + Xrr r 2 + Xrv r v + X 2 + Xu u 2 + Xv v + Xuv u v


Y = Yv v + Yrr + Yvvv v 3 + Yvvr v 2 r + Yvu vu + Yru r u

Extracted from Fossen (1994), Guidance and Control of Ocean Vehicles, Wiley and Sons Inc.

+ Y + Y 3 + Yu v + Yuu u 2 + Yv v 2 + Yvv v 2 + Y0 + Y0u u + Y0uu u 2

N = Nv v + Nr r + Nvvv v 3 + N vvr v 2 r + Nvu vu + Nrru r u

+ N + N 3 + N u v + N uu u 2 + N v v 2 + N vv v 2 + N0 + N0u u + N0uu u 2

The non-dimensional coefficients in the model are:


m = 798.10 5 , Iz = 39.2 10 5 , x = 0.023 G
Table 1 Non-dimensional hydrodynamic coefficients for the Mariner Class Vessel

X-equation u = 840 10 5 X X u = 184 10 5 X uu = 110 10 5

Y-equation = 1546 10 5 Yv Yr = 9 10 5 Yv = 1160 10 5 Yr = 499 10 5 Yvvv = 8078 10 5

N-equation v = 23 10 5 N N r = 83 10 5 N v = 264 10 5 N r = 166 10 5 N vvv = 1636 10 5

X uuu = 215 10 5
X vv = 899 10 5 X rr = 18 10 5 X = 95 10 5

Yvvr = 15356 10 5 Yvu = 1160 10 5


Yru = 499 10 5

Nvvr = 5483 10 5 Nvu = 264 10 5


N ru = 166 10 5

X u = 190 10 5 X rv = 798 10 5
X v = 93 10 5

Y = 278 10 5 Y = 90 10 5
Yu = 556 10 5

N = 139 10 5 N = 45 10 5
N u = 278 10 5

X uv = 93 10 5

Yuu = 278 10 5
Yv = 4 10 5

N uu = 139 10 5
N v = 13 10 5

Yvv = 1190 10 5
Y0 = 4 10 5 Y0u = 8 10 5
Y 0uu = 4 10 5

N vv = 489 10 5 N0 = 3 10 5 N0u = 6 10 5
N 0uu = 3 10 5

Extracted from Fossen (1994), Guidance and Control of Ocean Vehicles, Wiley and Sons Inc.

MATLAB M-file for Mariner Class Vessel (Fossen, 2004)


function [xdot,U] = mariner(x,ui,U0) % [xdot,U] = mariner(x,ui) returns the speed U in m/s (optionally) and the % time derivative of the state vector: x = [ u v r x y psi delta n ]' for % the Mariner class vessel L = 160.93 m, where % % u = pertubed surge velocity about Uo (m/s) % v = pertubed sway velocity about zero (m/s) % r = pertubed yaw velocity about zero (rad/s) % x = position in x-direction (m) % y = position in y-direction (m) % psi = pertubed yaw angle about zero (rad) % delta = actual rudder angle (rad) % % The inputs are : % % ui = commanded rudder angle (rad) % U0 = nominal speed (optionally). Default value is U0 = 7.7175 m/s = 15 knots. % % Reference: M.S. Chislett and J. Stroem-Tejsen (1965). Planar Motion Mechanism Tests % and Full-Scale Steering and Maneuvering Predictions for a Mariner Class Vessel, % Technical Report Hy-5, Hydro- and Aerodynamics Laboratory, Lyngby, Denmark. % % Author: Trygve Lauvdal % Date: 12th May 1994 % Revisions: 19th July 2001 (Thor I. Fossen): added input/ouput U0 and U, changed order of x-vector % 20th July 2001 (Thor I. Fossen): replaced inertia matrix with correct values % 11th July 2003 (Thor I. Fossen): max rudder is changed from 30 deg to 40 % deg to satisfy IMO regulations for 35 deg rudder execute % Check of input and state dimensions if (length(x) ~= 7),error('x-vector must have dimension 7 !'); end if (length(ui) ~= 1),error('ui must be a scalar input!'); end if nargin==2, U0 = 7.7175; end % Normalization variables L = 160.93; U = sqrt((U0 + x(1))^2 + x(2)^2); % Non-dimensional states and inputs delta_c = -ui; % delta_c = -ui such that positive delta_c -> positive r u v r psi delta = = = = = x(1)/U; x(2)/U; x(3)*L/U; x(6); x(7);

Extracted from Fossen (1994), Guidance and Control of Ocean Vehicles, Wiley and Sons Inc.

% Parameters, hydrodynamic derivatives and main dimensions delta_max = 40; % max rudder angle (deg) Ddelta_max = 5; % max rudder derivative (deg/s) m = 798e-5; Iz = 39.2e-5; xG = -0.023; Xudot Xu Xuu Xuuu Xvv Xrr Xdd Xudd Xrv Xvd Xuvd = = = = = = = = = = = -42e-5; -184e-5; -110e-5; -215e-5; -899e-5; 18e-5; -95e-5; -190e-5; 798e-5; 93e-5; 93e-5; Yvdot Yrdot Yv Yr Yvvv Yvvr Yvu Yru Yd Yddd Yud Yuud Yvdd Yvvd Y0 Y0u Y0uu = -748e-5; =-9.354e-5; = -1160e-5; = -499e-5; = -8078e-5; = 15356e-5; = -1160e-5; = -499e-5; = 278e-5; = -90e-5; = 556e-5; = 278e-5; = -4e-5; = 1190e-5; = -4e-5; = -8e-5; = -4e-5; Nvdot Nrdot Nv Nr Nvvv Nvvr Nvu Nru Nd Nddd Nud Nuud Nvdd Nvvd N0 N0u N0uu = 4.646e-5; = -43.8e-5; = -264e-5; = -166e-5; = 1636e-5; = -5483e-5; = -264e-5; = -166e-5; = -139e-5; = 45e-5; = -278e-5; = -139e-5; = 13e-5; = -489e-5; = 3e-5; = 6e-5; = 3e-5;

% Masses and moments of inertia m11 = m-Xudot; m22 = m-Yvdot; m23 = m*xG-Yrdot; m32 = m*xG-Nvdot; m33 = Iz-Nrdot; % Rudder saturation and dynamics if abs(delta_c) >= delta_max*pi/180, delta_c = sign(delta_c)*delta_max*pi/180; end delta_dot = delta_c - delta; if abs(delta_dot) >= Ddelta_max*pi/180, delta_dot = sign(delta_dot)*Ddelta_max*pi/180; end % Forces and moments X = Xu*u + Xuu*u^2 + Xuuu*u^3 + Xvv*v^2 + Xrr*r^2 + Xrv*r*v + Xdd*delta^2 +... Xudd*u*delta^2 + Xvd*v*delta + Xuvd*u*v*delta; Y = Yv*v + Yr*r + Yvvv*v^3 + Yvvr*v^2*r + Yvu*v*u + Yru*r*u + Yd*delta + ... Yddd*delta^3 + Yud*u*delta + Yuud*u^2*delta + Yvdd*v*delta^2 + ... Yvvd*v^2*delta + (Y0 + Y0u*u + Y0uu*u^2); N = Nv*v + Nr*r + Nvvv*v^3 + Nvvr*v^2*r + Nvu*v*u + Nru*r*u + Nd*delta + ... Nddd*delta^3 + Nud*u*delta + Nuud*u^2*delta + Nvdd*v*delta^2 + ... Nvvd*v^2*delta + (N0 + N0u*u + N0uu*u^2);

Extracted from Fossen (1994), Guidance and Control of Ocean Vehicles, Wiley and Sons Inc.
% Dimensional state derivative detM22 = m22*m33-m23*m32; xdot = [ X*(U^2/L)/m11 -(-m33*Y+m23*N)*(U^2/L)/detM22 (-m32*Y+m22*N)*(U^2/L^2)/detM22 (cos(psi)*(U0/U+u)-sin(psi)*v)*U (sin(psi)*(U0/U+u)+cos(psi)*v)*U r*(U/L) delta_dot

];

You might also like