Parametric Time-domain models based on frequency-domain data
|
|
- Alvin Stone
- 6 years ago
- Views:
Transcription
1 Parametric Time-domain models based on frequency-domain data (Module 7) Dr Tristan Perez Centre for Complex Dynamic Systems and Control (CDSC) Prof. Thor I Fossen Department of Engineering Cybernetics 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 1
2 Time-domain modelling approaches Two approaches can be distinguished for timedomain modelling: Full time-domain hydrodynamic codes, Time-domain models based on frequencydomain data. Here, we will focus on the second approach. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia
3 Frequency-domain Eq. of Motion In the hydrodynamic literature it is common to find the following model: [ M + A( ω)] & ξ ( t) + B( ω) ξ& ( t) + Gξ( t) = τ ( t) RB This is not a true time-domain model (Cummins, 196) exc This is valid to describe the steady-state response to sinusoidal excitations i.e., Frequency Response: ~ 1 ( ) ω [ MRB + A( ω)] + jω B( ) + G τexc ξ = ω ~ ~ denotes complex variable Force to Motion RAO 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 3
4 Cummins s equation Cummins (196), took a different modelling approach and consider the radiation problem ab initio in the time domain : The added mass matrix is constant frequency and speed independent. The convolution term accounts for fluid-memory effects. The kernel of the convolution is a matrix of retardation functions or impulse responses. This is a true linear time-domain model. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 4
5 Cummins s equation with forward speed The convolution terms depend on the forward speed. With forward speed appears a constant damping term. The restoring forces are affected by hydrodynamic pressure Lift, changes in trim. (Usually ignored for Fn<.3) 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 5
6 Ogilvie s relations If Cummins s Equation is valid for any input, it must then be valid for sinusoids in particular (Ogilvie, 1964). Ogilvie (1964) transformed the Cummins Equation to the frequency domain, and found that A( ω) = A 1 ω From the Riemann-Lesbesgue Lemma: B( ω) = B( U ) + ω B( U ) = K( t)sin( ωt) dt K( t)cos( ωt) dt A = lim A( ω) : = A( ) lim B( ω) : = ω B( ) 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 6
7 Non-parametric Representations Time-domain K( t) = [ B( ω) B( )]cos( ωt) dω π Frequency-domain K( jω) = [ B( ω) B( )] + jω [ A( ω) A( )] 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 7
8 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 8 Parametric Representations Because the convolution is a dynamic linear operation, it can be represented by a linear ordinary differential equation state-space model: = ) ( ) ( τ τ τ d t r ξ K μ & C x μ B ξ A x x c r c c = + = & & = x n x x M 1 x
9 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 9 Parametric Representations From the state-space representation, it follow = = ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( s Q s P s Q s P s Q s P s Q s P s s c c c L M O M L B A I C K c c c t t B A C K ) exp( ) ( = Impulse Response Frequency-domain model (Transfer Function matrix) ) ( ) ( q s q s b s b s b s Q s P n n n m m m m ij ij = L L Rational TF
10 Convolution replacement The convolution (non-parametric model) in the Cummins equation can be time and memory consuming for simulation. For analysis and design of a control system, the convolutions are not very well suited. The parametric state-space representation of appropriate order (n-order) eliminates the above problems. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 1
11 Convolution replacement τ w + μ r [ M + A( )] 1 K ( t τ ) ξ& ( τ ) dτ ξ & ξ G μ = K( t τ ) ξ& ( τ ) dτ r x& = A x + B ξ& μ r = c C x c c 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 11
12 Properties of the convolution terms Property limk( jω) = B( ) ω Implication on parametric models K(s) is zero at s= for U=. lim K ( jω) = ω + K( t = ) = [ B( ω) B( )] dω TFs strictly proper TFs relative degree 1 limk ( t) = t Re{ K ii ( jω)} Note: Bold symbols denote matrices. TF BIBO stable K(s) is passive => diagonal terms are positive real; off diagonal terms stable. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 1
13 Low-frequency limit (U=) K( jω) = [ B( ω) B( )] + jω [ A( ω) A( )] In the limit at low freq, B(ω) is zero, since there cannot be waves (Faltinsen, 199); thus the real part is zero for U= and -B( ) for U>. The imaginary part tends to zero as the following difference is finite: Note that regularity conditions for the exchange of limit and integration are satisfied. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 13
14 High-frequency limit (U=) K( jω) = [ B( ω) B( )] + jω [ A( ω) A( )] In the limit at high frequency, the real part is zero. The imaginary part also tends to zero by Ogilvie s relation and Riemann-Lebesgue Lemma: 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 14
15 Initial and final time (U=) Initial time: lim K( t) = lim [ B( ω) B( )] cos( ω t) dω = [ B( ω) B( )] dω + + t t π π Regularity conditions for the exchange of limit and integration are satisfied. The last relation follows from energy considerations (Faltinsen, 199): Final time: limk ( t) = lim [ B( ω) B( )] cos( ω t) dω = t t π Which follows by Ogilvie s relation and Riemann-Lebesgue Lemma. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 15
16 Passivity For U= and no current, the damping matrix is symmetric and positive semi-definite: From this follows the positive realness of the convolution terms and thus the passivity; that is these terms cannot generate energy. From energy considerations, it also follows that the diagonal terms of K(s) are passive. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 16
17 Parametric model identification The convolution replacement can be posed in different ways, which in theory should provide the same answer: B (ω) K(t) Time-domain identification Aˆ C ˆ c c Bˆ Dˆ c c A( ω), B( ω) Frequency-domain identification K ( jω) K ˆ ( s ) Model conversion Aˆ C ˆ c c Bˆ Dˆ c c In practice one method can be more favourable than the other. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 17
18 Parametric model identification Different proposals have appeared in the literature: Time-domain identification: LS-fitting of the impulse response (Yu & Falnes, 1998) Realization theory (Kristiansen & Egeland, 3) Frequency-domain identification: LS-fitting of the frequency response K(jω) (Jeffreys, 1984),(Damaren ). LS-fitting of added mass and damping (Soding 198), (Xia et. al 1998), (Sutulo & Guedes-Soares 6). 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 18
19 Time-domain identification From K(t) to state-space models. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 19
20 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia Numerical computations of K(t) A key issue for time-domain identification is to start with a good impulse response computed from the damping. Numerical codes can only provide accurate computations of added mass and damping up to a certain frequency, say Ω. This introduces an error in the computation of the retardation functions: ω ω ω π ω ω ω π d t U d t U U t ) )]cos( ( ) ( [ ) )]cos( ( ) ( [ ), ( Ω Ω + = B B B B K ω ω ω π d t U U t ) )]cos( ( ) ( [ ), ( Ω B B K ω ω ω π d t U U t ) )]cos( ( ) ( [ ), ( Ω = B B Error
21 High-frequency values of A(ω) and B(ω) In the limit at high frequency the following tendencies are observed for the 3D damping and added mass: B αik βik ( ω) as ω Aik ( ω) Aik ( ) as ω ω ω ik As commented by Damaren (), this seems at odds with what is generally stated in the hydrodynamic literature! Note that there are no expansions involved to obtain these results, the only assumption is the linearity which results in a rational representation and the relative degree 1, which results from the integration of damping over the frequencies. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 1
22 Example Containership (Taghipour et al., 7a) The panel sizing was done to be able to compute frequencies up to.5 rad/s. Rule of thumb: characteristic panel length < 1/8 min wave length (Faltinsen, 1993). 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia
23 Numerical computations Extending the damping with tail prop to 1/ω B33 [Kg/s] B53 [Kgm/s rad] 6 x 17 4 B33 B33 ext w Freq. [rad/s] 6 x 18 4 B53 B53 ext w - B53 ext w - +w Freq. [rad/s] B35 [Kgm/s rad] B55 [Kg m /s] 6 x 18 4 B35 B35 ext w - B35 ext w - +w Freq. [rad/s] 3 x B55 B55 ext w Freq. [rad/s] In this example, the tail α/ω is not a very good for B35 and B53 (.5rad/s is too low), whereas it is ok for B33 and B55. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 3
24 Numerical computations Sometimes it is necessary to extend the damping with asymptotic values: x 1 11 B55 B55 ext w x 111 Ω K( t) Bext ( ω)cos( ωt) dω π K55 K55ext B55 [Kg m /s] K Freq. [rad/s] time [s] When doing time-domain identification, it is important to start from a good approx of the retardation function! Note the differences at t= +, and the errors at different time instants. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 4
25 Impulse response curve fitting Given the SISO SS realization of order n: The parameters can be obtained from The application of this method to marine structures was proposed by Yu and Falnes (1998). 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 5
26 Impulse response curve fitting It is hard to guess the order of the system by looking at the impulse response alone one should start with lower order and increase it to improve the fit. The LS-problem is non-linear in the parameters. This problem can be solved with Gaussian-Newton methods. The Gaussian-Newton methods are known to work well if the parameters initial guess are close to the optimal parameters. Since the initial value of the parameters is difficult to obtain from the impulse response, and these depends on the particular realization being chosen, the method is not very practical. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 6
27 Realization theory A key result of realization theory is the following factorization (Ho and Kalman, 1966): Hankel matrix of the impulse response values (constant along the anti-diagonals) Extended controllability matrix Extended observability matrix 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 7
28 Realization theory Kung s Algorithm (Kung, 1978): Singular value decomposition The number of significant singular values give the order of the system: 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 8
29 Realization theory Kung s Algorithm (Kung, 1978): The application of this method to marine structures was proposed by Kristiansen and Egeland (3). 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 9
30 Realization theory The problem is solved in discrete time Kung s algorithm obtains the model based on a SVD-factorization of the Hankel matrix of samples of the impulse response. If the impulse response is not accurate, it may in very large order systems The conversion from discrete to continuous often gives a matrix Dc in the state-space realization, which is inconsistent with the dynamics of the problem for the retardation function (relative degree 1). The MATLAB command impss implements Kung s algorithm, and chooses the order by neglecting singular values less than 1% of the largest one. impss requires using model order reduction afterwards (Kristiansen et al 5). The resulting models may not be passive. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 3
31 Example Container ship Singular values of the Hankel matrix of the samples of the impulse response. Normalised Singular Values K33 K33ext Normalised Singular Values K35 K35ext These suggest (based on the blue plots) Normalised Singular Values K53 K53ext Normalised Singular Values K55 K55ext Order K33(s) = 3 or 4 Oder K35(s) = 5 or 6 Order K55(s) = 3 or /9/7 One-day Tutorial, CAMS'7, Bol, Croatia 31
32 Example Container ship Impulse response fitting for K33(t) with a system of order 4. Identification method: impss + balmr (model order reduction). K33(t) K33(jw) Retardation function 3.5 x Data Approximation order Time [s] K33(jw) db Phase K33(jw) [deg] Convolution Frequency Response Freq. [rad/s] K33(jw) K33(jw) order Freq. [rad/s] The identified model is not passive and does not have a zero at s=. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 3
33 Example Container ship Impulse response fitting for K33(t) with a system of order 3. Identification method: impss + balmr (model order reduction). Retardation function 3.5 x K33(t) Data Approximation order Time [s] K33(jw) db Phase K33(jw) [deg] Convolution Frequency Response K33(jw) Freq. [rad/s] K33(jw) K33(jw) order Freq. [rad/s] The identified model is passive, but still does not have a zero at s=. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 33
34 Example Container ship Impulse response fitting for K35(t) with a system of order 5. 1 x 18 K35(t) Data Approximation order K35(jw) Convolution Frequency Response 8 K35(jw) db K35(jw) K35(jw) order 5 Retardation function Freq. [rad/s] Phase K35(jw) [deg] Time [s] Freq. [rad/s] The identified model is passive, but still does not have a zero at s=. Note that the off-diagonal terms do not necessarily have to passive for K(s) to be passive (Unneland & Perez, 7). 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 34
35 Comments about Realization Theory Depending on the hydrodynamic data, it may be necessary to extend the damping at high freq. to have a good estimate the of the impulse response function before doing the idetification. Looking at the impulse response fitting alone is not a good criteria most properties are evident from the freq. response. Impss may require using model order reduction afterwards. The models almost never satisfy the low frequency asymptotic values (zero at s=). High-order models may not be passive; this can solved trying different orders or using a model order reduction method that enforces passivity. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 35
36 Frequency-domain Identification From K(jω) to K(s) 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 36
37 Frequency-domain identification We can fit a SISO TF to each entry of K(jω): Kˆ ( s) ij = P Q ij ij ( s) ( s) Where, P ij (s) deg Q ( s) = deg P ( s) + 1 has a zero at s= for U= or constant for U> Relative degree = 1; i.e., Stable ij ij 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 37
38 Relative degree condition From the finite initial time of the impulse response: lim K ( ) = lim [ ( ) ( )] cos( ) = [ ( ) ( )] + ij t + Bij ω Bij ω t dω Bij ω Bij dω t t π π From the Initial-value Theorem of the Laplace Transform: lim t + K ij ( t) = lim sk s ij ( s) = lim s sp ( s) Q ij ij ( s) = b m s s m+ 1 n. Hence for this to be finite n=m+1; relative degree = 1. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 38
39 Minimum order transfer function Because of the restriction of relative degree 1, the minimum order TF that can represent a convolution term is min K ij ( s) = s s b1s + a s + a 1 b1s + b + a s + a 1 U U = >,, Therefore, we can start with a system of order n=, and then increase the order until we improve the fitting at an appropriate level. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 39
40 Regression in the frequency domain In this method, a rational transfer function is fitted to the frequency response data: The application of this method to marine structures was proposed by Jeffreys (1984) and Damaren (). 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 4
41 Quasi-linear regression Levi (1959) proposed the following linearization: This can be obtained If we chose the weights in the nonlinear problem as which is affine in the parameters and reduces to a linear LS problem. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 41
42 Iterative Quasi-linear regression The quasi-linear regressor tend to have a poor fit a low freq. This can be avoided by solving the linear LS problem iteratively, starting with the quasi-linear regressor and using the parameters obtained to compute a weighing: k =,3,K After a few iterations, nonlinear problem is recovered. and the 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 4
43 Containership example Frequency-domain identification of K33(jw) order 3 Identification method: iterative quasi-linear regression (invfreqs.m) 16 Convolution Frequency Response 15 K33(jw) Freq. [rad/s] K33(jw) K33(jw) order 3 The model is passive and satisfy the asymptotic values. 1 Phase K33(jw) [deg] Freq. [rad/s] 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 43
44 Containership example Frequency-domain identification of K35(jw) order 5 and K55(jw) order 3 Identification method: iterative quasi-linear regression (invfreqs.m) 18 Convolution Frequency Response 3 Convolution Frequenc y Response 17 K35(jw) K35(jw) K35(jw) order 5 K55(jw) 1 19 K55(jw) K55(jw) order Freq. [rad/s] Freq. [rad/s] 1 1 Phase K35(jw) [deg] 5-5 Phase K55(jw) [deg] Freq. [rad/s] Freq. [rad/s] 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 44
45 Reconstructing B(ω) B33 [Kg/s] 6 x 17 4 B33 B33 FD ident, order 3 B35 [Kg/s] 6 x 18 4 B35 B35 FD ident, order Freq. [rad/s] 4 6 Freq. [rad/s] B53 [Kg/s] 6 x 18 4 B53 B53 FD ident, order 5 B55 [Kg/s] 3 x B55 B55 FD ident, order Freq. [rad/s] Freq. [rad/s] 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 45
46 Reconstructing A(ω).6 x A33 A33 FD ident, order 3 A33inf 3 x A35 A35 FD ident, order 5 A35inf A33 [Kg] A35 [Kgm] Freq. [rad/s] Freq. [rad/s] 3 x 19 9 x 111 A53 [Kgm] A53 A53 TD ident, order 5 A53inf A55 [Kg m ] A55 A55 TD ident, order 3 A55inf Freq. [rad/s] Freq. [rad/s] 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 46
47 Comments about Freq.-dom. regression FD-identification avoids having to compute the impulse response from the damping The identification method is simplest: a series of linear LS problems easy to programme. The zero at s= and the relative degree can be enforced in the structure of the model, so the asymptotic values are always ensured. The resulting models may be unstable: this is fixed by reflecting the unstable poles about the imaginary axis. The resulting models may not be passive: this can be solved using weights in the LS problem. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 47
48 Simulink Model implementation After obtaining a state-space representation or the transfers functions, we can assemble a complete model: 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 48
49 References Perez, T. and T. I. Fossen (6) Time-domain Models of Marine Surface Vessels Based on Seakeeping Computations. 7th IFAC Conference on Manoeuvring and Control of Marine Vessels MCMC, Portugal, September. Perez, T. (7) Identification and Validation of a Time-domain Hydrodynamic Model of a Prototype Hull for a Wave Energy Converter (WEC). Technical Report. Centre for Complex Dynamic Systems and Control. Australia. Cummins, W., 196. The impulse response function and ship motions. Schiffstechnik 9 (1661), Ogilvie, T., Recent progress towards the understanding and prediction of ship motions. In: 6th Symposium on Naval Hydrodynamics. Kristiansen, E., Egeland, O., 3. Frequency-dependent added mass in models for controller design for wave motion damping. In: Proceedings of 6th Conference on Maneoeuvering and Control ofmarine Craft, Girona, Spain. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 49
50 References Damaren, C.J. (). Time-domain floating body dynamics by rational approximations of the radiation impedance and diffraction mapping. Ocean Engineering 7, Yu, Z., Falnes, J., State-space modelling of dynamic systems in ocean engineering. Journal of Hydrodynamics B(1), Sutulo, S., Guedes-Soares, C., 5. An implementation of the method of auxiliary state variables for solving seakeeping problems. International shipbuilding progress (Int. shipbuild. prog.) 5, Taghipour, R., Perez, T., Moan, T., 7a. Hybrid Frequency-Time Domain Models for Dynamic Response Analysis of Marine Structures. (To appear in Ocean Engieering) Taghipour, R., Perez, T., Moan, T., 7b. Time domain hydroelastic analysis of a flexible marine structure using state-space models. In: 6th International Conference on Offshore Mechanics and Arctic Engineering-OMAE 7, San Diego, CA, USA. 19/9/7 One-day Tutorial, CAMS'7, Bol, Croatia 5
A DERIVATION OF HIGH-FREQUENCY ASYMPTOTIC VALUES OF 3D ADDED MASS AND DAMPING BASED ON PROPERTIES OF THE CUMMINS EQUATION
Journal of Maritime Research, Vol. V. No., pp. 65-78, 8 Copyright 8. SEECMAR Printed in Santander (Spain). All rights reserved ISSN: 697-484 A DERIVATION OF HIGH-FREQUENCY ASYMPTOTIC VALUES OF 3D ADDED
More informationModelling and Identification of Radiation-force Models for Ships and Marine Structures
Modelling and Identification of Radiation-force Models for Ships and Marine Structures A/Prof Tristan Perez Leader Autonomous Systems & Robotics School of Engineering, the University of Newcastle AUSTRALIA
More informationTime- vs. Frequency-domain Identification of Parametric Radiation Force Models for Marine Structures at Zero Speed
Modeling, Identification and Control, Vol. 9, No., 8, pp. 9 Time- vs. Frequency-domain Identification of Parametric Radiation Force Models for Marine Structures at Zero Speed Tristan Pérez Thor I. Fossen
More informationSeakeeping Models in the Frequency Domain
Seakeeping Models in the Frequency Domain (Module 6) Dr Tristan Perez Centre for Complex Dynamic Systems and Control (CDSC) Prof. Thor I Fossen Department of Engineering Cybernetics 18/09/2007 One-day
More informationIdentification of Dynamic Models of Marine Structures from Frequency-domain Data Enforcing Model Structure and Parameter Constraints
ARC Centre of Excellence for Complex Dynamic Systems and Control, pp. 1 8 Identification of Dynamic Models of Marine Structures from Frequency-domain Data Enforcing Model Structure and Parameter Constraints
More informationA Discussion About Seakeeping and Manoeuvring Models For Surface Vessels
A Discussion About Seakeeping and Manoeuvring Models For Surface Vessels Tristan Perez, Thor I. Fossen and Asgeir Sørensen Technical Report (MSS-TR-001) Marine System Simulator (MSS) Group (http://www.cesos.ntnu.no/mss/)
More information2 nd SDWED Symposium Advances in Modelling of Wave Energy Devices - WP4
2 nd SDWED Symposium Advances in Modelling of Wave Energy Devices - WP4 Andrew Zurkinden Wave Energy Research Group, Civil Engineering Department, Aalborg University, Denmark April 26, 212 Structural Design
More informationTHE EFFECT OF MEMORY IN PASSIVE NONLINEAR OBSERVER DESIGN FOR A DP SYSTEM
DYNAMIC POSIIONING CONFERENCE October 1-13, 1 DESIGN SESSION HE EFFEC OF MEMORY IN PASSIVE NONLINEAR OBSERVER DESIGN FOR A DP SYSEM By A. Hajivand & S. H. Mousavizadegan (AU, ehran, Iran) ABSRAC he behavior
More informationLecture 1: Introduction to System Modeling and Control. Introduction Basic Definitions Different Model Types System Identification
Lecture 1: Introduction to System Modeling and Control Introduction Basic Definitions Different Model Types System Identification What is Mathematical Model? A set of mathematical equations (e.g., differential
More informationThe Laplace Transform
The Laplace Transform Generalizing the Fourier Transform The CTFT expresses a time-domain signal as a linear combination of complex sinusoids of the form e jωt. In the generalization of the CTFT to the
More informationPredictive Control of Gyroscopic-Force Actuators for Mechanical Vibration Damping
ARC Centre of Excellence for Complex Dynamic Systems and Control, pp 1 15 Predictive Control of Gyroscopic-Force Actuators for Mechanical Vibration Damping Tristan Perez 1, 2 Joris B Termaat 3 1 School
More informationA system that is both linear and time-invariant is called linear time-invariant (LTI).
The Cooper Union Department of Electrical Engineering ECE111 Signal Processing & Systems Analysis Lecture Notes: Time, Frequency & Transform Domains February 28, 2012 Signals & Systems Signals are mapped
More informationQualification of tabulated scattering parameters
Qualification of tabulated scattering parameters Stefano Grivet Talocia Politecnico di Torino, Italy IdemWorks s.r.l. stefano.grivet@polito.it 4 th IEEE Workshop on Signal Propagation on Interconnects
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 12: I/O Stability Readings: DDV, Chapters 15, 16 Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology March 14, 2011 E. Frazzoli
More informationKari Unneland. Identification and Order Reduction of Radiation Force Models of Marine Structures. Theses at NTNU, 2007:188 Kari Unneland
Kari Unneland Theses at NTNU, 27:188 Kari Unneland Identification and Order Reduction of Radiation Force Models of Marine Structures ISBN 978-82-471-414-3 (printed ver.) ISBN 978-82-471-4154- (electronic
More informationTime Response Analysis (Part II)
Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary
More informationDynamic Response. Assoc. Prof. Enver Tatlicioglu. Department of Electrical & Electronics Engineering Izmir Institute of Technology.
Dynamic Response Assoc. Prof. Enver Tatlicioglu Department of Electrical & Electronics Engineering Izmir Institute of Technology Chapter 3 Assoc. Prof. Enver Tatlicioglu (EEE@IYTE) EE362 Feedback Control
More informationMAE 143B - Homework 9
MAE 143B - Homework 9 7.1 a) We have stable first-order poles at p 1 = 1 and p 2 = 1. For small values of ω, we recover the DC gain K = lim ω G(jω) = 1 1 = 2dB. Having this finite limit, our straight-line
More informationRaktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Dynamic Response
.. AERO 422: Active Controls for Aerospace Vehicles Dynamic Response Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. . Previous Class...........
More informationAMME3500: System Dynamics & Control
Stefan B. Williams May, 211 AMME35: System Dynamics & Control Assignment 4 Note: This assignment contributes 15% towards your final mark. This assignment is due at 4pm on Monday, May 3 th during Week 13
More informationEfficient Prediction of Dynamic Response for Flexible and Multi- Body Marine Structures
Reza Taghipour Efficient Prediction of Dynamic Response for Flexible and Multi- Body Marine Structures Thesis for the degree of philosophiae doctor Trondheim, October 28 Norwegian University of Science
More informationThe Laplace Transform
The Laplace Transform Syllabus ECE 316, Spring 2015 Final Grades Homework (6 problems per week): 25% Exams (midterm and final): 50% (25:25) Random Quiz: 25% Textbook M. Roberts, Signals and Systems, 2nd
More informationThe Laplace Transform
The Laplace Transform Introduction There are two common approaches to the developing and understanding the Laplace transform It can be viewed as a generalization of the CTFT to include some signals with
More informationProblem Set 5 Solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 5 Solutions The problem set deals with Hankel
More informationPrüfung Regelungstechnik I (Control Systems I) Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 29. 8. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 13: Stability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 13: Stability p.1/20 Outline Input-Output
More informationSIGNALS AND SYSTEMS LABORATORY 4: Polynomials, Laplace Transforms and Analog Filters in MATLAB
INTRODUCTION SIGNALS AND SYSTEMS LABORATORY 4: Polynomials, Laplace Transforms and Analog Filters in MATLAB Laplace transform pairs are very useful tools for solving ordinary differential equations. Most
More informationTrajectory Tracking of a Near-Surface Torpedo using Numerical Methods
ISSN (Print) : 2347-671 An ISO 3297: 27 Certified Organization Vol.4, Special Issue 12, September 215 Trajectory Tracking of a Near-Surface Torpedo using Numerical Methods Anties K. Martin, Anubhav C.A.,
More informationSubspace Identification
Chapter 10 Subspace Identification Given observations of m 1 input signals, and p 1 signals resulting from those when fed into a dynamical system under study, can we estimate the internal dynamics regulating
More informationChap 4. State-Space Solutions and
Chap 4. State-Space Solutions and Realizations Outlines 1. Introduction 2. Solution of LTI State Equation 3. Equivalent State Equations 4. Realizations 5. Solution of Linear Time-Varying (LTV) Equations
More informationIntroduction & Laplace Transforms Lectures 1 & 2
Introduction & Lectures 1 & 2, Professor Department of Electrical and Computer Engineering Colorado State University Fall 2016 Control System Definition of a Control System Group of components that collectively
More information( ) ( = ) = ( ) ( ) ( )
( ) Vρ C st s T t 0 wc Ti s T s Q s (8) K T ( s) Q ( s) + Ti ( s) (0) τs+ τs+ V ρ K and τ wc w T (s)g (s)q (s) + G (s)t(s) i G and G are transfer functions and independent of the inputs, Q and T i. Note
More informationA simplified method for calculating propeller thrust decrease for a ship sailing on a given shipping lane
POLISH MARITIME RESEARCH 2(82) 2014 Vol 21; pp. 27-33 10.2478/pomr-2014-0015 A simplified method for calculating propeller thrust decrease for a ship sailing on a given shipping lane Katarzyna Zelazny,
More informationDESIGN OF A HYBRID POWER/TORQUE THRUSTER CONTROLLER WITH LOSS ESTIMATION. Øyvind N. Smogeli, Asgeir J. Sørensen and Thor I. Fossen
DESIGN OF A HYBRID POWER/TORQUE THRUSTER CONTROLLER WITH LOSS ESTIMATION Øyvind N. Smogeli, Asgeir J. Sørensen and Thor I. Fossen Department of Marine Technology Norwegian University of Science and Technology
More informationGATE EE Topic wise Questions SIGNALS & SYSTEMS
www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)
More informationOutline. Classical Control. Lecture 1
Outline Outline Outline 1 Introduction 2 Prerequisites Block diagram for system modeling Modeling Mechanical Electrical Outline Introduction Background Basic Systems Models/Transfers functions 1 Introduction
More informationÜbersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 3.. 24 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -
More informationExercise 1 (A Non-minimum Phase System)
Prof. Dr. E. Frazzoli 5-59- Control Systems I (HS 25) Solution Exercise Set Loop Shaping Noele Norris, 9th December 26 Exercise (A Non-minimum Phase System) To increase the rise time of the system, we
More informationDynamic circuits: Frequency domain analysis
Electronic Circuits 1 Dynamic circuits: Contents Free oscillation and natural frequency Transfer functions Frequency response Bode plots 1 System behaviour: overview 2 System behaviour : review solution
More informationEvaluation of Hydrodynamic Performance of a Damaged Ship in Waves
Evaluation of Hydrodynamic Performance of a Damaged Ship in Waves Sa Young Hong, Seok-Kyu Cho, Byoung Wan Kim, Gyeong Jung Lee, Ki-Sup Kim Maritime and Ocean Engineering Research Institute/KORDI, Daejeon,
More informationFrequency domain analysis
Automatic Control 2 Frequency domain analysis Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011
More informationSubject: Optimal Control Assignment-1 (Related to Lecture notes 1-10)
Subject: Optimal Control Assignment- (Related to Lecture notes -). Design a oil mug, shown in fig., to hold as much oil possible. The height and radius of the mug should not be more than 6cm. The mug must
More informationControl System Design
ELEC ENG 4CL4: Control System Design Notes for Lecture #4 Monday, January 13, 2003 Dr. Ian C. Bruce Room: CRL-229 Phone ext.: 26984 Email: ibruce@mail.ece.mcmaster.ca Impulse and Step Responses of Continuous-Time
More informationIntroduction. Performance and Robustness (Chapter 1) Advanced Control Systems Spring / 31
Introduction Classical Control Robust Control u(t) y(t) G u(t) G + y(t) G : nominal model G = G + : plant uncertainty Uncertainty sources : Structured : parametric uncertainty, multimodel uncertainty Unstructured
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 2.6 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS
More informationProblem Set 3: Solution Due on Mon. 7 th Oct. in class. Fall 2013
EE 56: Digital Control Systems Problem Set 3: Solution Due on Mon 7 th Oct in class Fall 23 Problem For the causal LTI system described by the difference equation y k + 2 y k = x k, () (a) By first finding
More informationSimplified Numerical Model for the Study of Wave Energy Convertes (WECs)
Simplified Numerical Model for the Study of Wave Energy Convertes (WECs) J. A. Armesto, R. Guanche, A. Iturrioz and A. D. de Andres Ocean Energy and Engineering Group IH Cantabria, Universidad de Cantabria
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014
Prof. Dr. Eleni Chatzi System Stability - 26 March, 24 Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can
More information1 POTENTIAL FLOW THEORY Formulation of the seakeeping problem
1 POTENTIAL FLOW THEORY Formulation of the seakeeping problem Objective of the Chapter: Formulation of the potential flow around the hull of a ship advancing and oscillationg in waves Results of the Chapter:
More informationExercise 1 (A Non-minimum Phase System)
Prof. Dr. E. Frazzoli 5-59- Control Systems I (Autumn 27) Solution Exercise Set 2 Loop Shaping clruch@ethz.ch, 8th December 27 Exercise (A Non-minimum Phase System) To decrease the rise time of the system,
More informationAn Internal Stability Example
An Internal Stability Example Roy Smith 26 April 2015 To illustrate the concept of internal stability we will look at an example where there are several pole-zero cancellations between the controller and
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering Dynamics and Control II Fall 2007
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering.4 Dynamics and Control II Fall 7 Problem Set #9 Solution Posted: Sunday, Dec., 7. The.4 Tower system. The system parameters are
More informationFEEDBACK CONTROL SYSTEMS
FEEDBAC CONTROL SYSTEMS. Control System Design. Open and Closed-Loop Control Systems 3. Why Closed-Loop Control? 4. Case Study --- Speed Control of a DC Motor 5. Steady-State Errors in Unity Feedback Control
More informationScientific Computing
Scientific Computing Direct solution methods Martin van Gijzen Delft University of Technology October 3, 2018 1 Program October 3 Matrix norms LU decomposition Basic algorithm Cost Stability Pivoting Pivoting
More informationSimplified formulas of heave added mass coefficients at high frequency for various two-dimensional bodies in a finite water depth
csnak, 2015 Int. J. Nav. Archit. Ocean Eng. (2015) 7:115~127 http://dx.doi.org/10.1515/ijnaoe-2015-0009 pissn: 2092-6782, eissn: 2092-6790 Simplified formulas of heave added mass coefficients at high frequency
More informationLecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.
ISS0031 Modeling and Identification Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. Aleksei Tepljakov, Ph.D. September 30, 2015 Linear Dynamic Systems Definition
More informationLaplace Transform Part 1: Introduction (I&N Chap 13)
Laplace Transform Part 1: Introduction (I&N Chap 13) Definition of the L.T. L.T. of Singularity Functions L.T. Pairs Properties of the L.T. Inverse L.T. Convolution IVT(initial value theorem) & FVT (final
More informationComparison of Present Wave Induced Load Criteria with Loads Induced by an Abnormal Wave
Rogue Waves 2004 1 Comparison of Present Wave Induced Load Criteria with Loads Induced by an Abnormal Wave C. Guedes Soares, N. Fonseca, R. Pascoal Unit of Marine Engineering and Technology, Technical
More informationLecture Note #6 (Chap.10)
System Modeling and Identification Lecture Note #6 (Chap.) CBE 7 Korea University Prof. Dae Ryook Yang Chap. Model Approximation Model approximation Simplification, approximation and order reduction of
More informationProceedings of the ASME 27th International Conference on Offshore Mechanics and Arctic Engineering OMAE2008 June 15-20, 2008, Estoril, Portugal
Proceedings of the ASME 7th International Conference on Offshore Mechanics and Arctic Engineering OMAE8 June 15-, 8, Estoril, Portugal OMAE8-57935 DISTRIBUTION OF MAXIMA OF NON-LINEAR ROLLING IN CASE OF
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 2 Laplace Transform I 1/52
1/52 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 2 Laplace Transform I Linear Time Invariant Systems A general LTI system may be described by the linear constant coefficient differential equation: a n d n
More informationLINEAR RESPONSE THEORY
MIT Department of Chemistry 5.74, Spring 5: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff p. 8 LINEAR RESPONSE THEORY We have statistically described the time-dependent behavior
More informationDesign and Tuning of Fractional-order PID Controllers for Time-delayed Processes
Design and Tuning of Fractional-order PID Controllers for Time-delayed Processes Emmanuel Edet Technology and Innovation Centre University of Strathclyde 99 George Street Glasgow, United Kingdom emmanuel.edet@strath.ac.uk
More informationTime domain assessment of nonlinear coupled ship motions and sloshing in free surface tanks
Time domain assessment of nonlinear coupled ship motions and sloshing in free surface tanks 1 outline 1.Motivation & state-of-the-art 2.Simulation approach 1.SHIXDOF: nonlinear ship motion TD 6DOF 2.AQUAgpusph:
More informationAvailable online at ScienceDirect. Energy Procedia 94 (2016 )
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 94 (16 ) 9 35 13th Deep Sea Offshore Wind R&D Conference, EERA DeepWind 16, - January 16, Trondheim, Norway Parametric Wave Excitation
More informationMEM 355 Performance Enhancement of Dynamical Systems
MEM 355 Performance Enhancement of Dynamical Systems Frequency Domain Design Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 5/8/25 Outline Closed Loop Transfer Functions
More informationRichiami di Controlli Automatici
Richiami di Controlli Automatici Gianmaria De Tommasi 1 1 Università degli Studi di Napoli Federico II detommas@unina.it Ottobre 2012 Corsi AnsaldoBreda G. De Tommasi (UNINA) Richiami di Controlli Automatici
More information2.010 Fall 2000 Solution of Homework Assignment 1
2. Fall 2 Solution of Homework Assignment. Compact Disk Player. This is essentially a reprise of Problems and 2 from the Fall 999 2.3 Homework Assignment 7. t is included here to encourage you to review
More informationDesign of a Heading Autopilot for Mariner Class Ship with Wave Filtering Based on Passive Observer
Design of a Heading Autopilot for Mariner Class Ship with Wave Filtering Based on Passive Observer 1 Mridul Pande, K K Mangrulkar 1, Aerospace Engg Dept DIAT (DU), Pune Email: 1 mridul_pande000@yahoo.com
More informationLinear State Feedback Controller Design
Assignment For EE5101 - Linear Systems Sem I AY2010/2011 Linear State Feedback Controller Design Phang Swee King A0033585A Email: king@nus.edu.sg NGS/ECE Dept. Faculty of Engineering National University
More informationOutline. Classical Control. Lecture 5
Outline Outline Outline 1 What is 2 Outline What is Why use? Sketching a 1 What is Why use? Sketching a 2 Gain Controller Lead Compensation Lag Compensation What is Properties of a General System Why use?
More informationHankel Optimal Model Reduction 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Hankel Optimal Model Reduction 1 This lecture covers both the theory and
More informationEE102 Homework 2, 3, and 4 Solutions
EE12 Prof. S. Boyd EE12 Homework 2, 3, and 4 Solutions 7. Some convolution systems. Consider a convolution system, y(t) = + u(t τ)h(τ) dτ, where h is a function called the kernel or impulse response of
More informationf (t) K(t, u) d t. f (t) K 1 (t, u) d u. Integral Transform Inverse Fourier Transform
Integral Transforms Massoud Malek An integral transform maps an equation from its original domain into another domain, where it might be manipulated and solved much more easily than in the original domain.
More informationELG 3150 Introduction to Control Systems. TA: Fouad Khalil, P.Eng., Ph.D. Student
ELG 350 Introduction to Control Systems TA: Fouad Khalil, P.Eng., Ph.D. Student fkhalil@site.uottawa.ca My agenda for this tutorial session I will introduce the Laplace Transforms as a useful tool for
More informationEE/ME/AE324: Dynamical Systems. Chapter 7: Transform Solutions of Linear Models
EE/ME/AE324: Dynamical Systems Chapter 7: Transform Solutions of Linear Models The Laplace Transform Converts systems or signals from the real time domain, e.g., functions of the real variable t, to the
More informationModel Reduction for Unstable Systems
Model Reduction for Unstable Systems Klajdi Sinani Virginia Tech klajdi@vt.edu Advisor: Serkan Gugercin October 22, 2015 (VT) SIAM October 22, 2015 1 / 26 Overview 1 Introduction 2 Interpolatory Model
More informationCHAPTER 1 Basic Concepts of Control System. CHAPTER 6 Hydraulic Control System
CHAPTER 1 Basic Concepts of Control System 1. What is open loop control systems and closed loop control systems? Compare open loop control system with closed loop control system. Write down major advantages
More informationControl Systems I. Lecture 5: Transfer Functions. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 5: Transfer Functions Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 20, 2017 E. Frazzoli (ETH) Lecture 5: Control Systems I 20/10/2017
More informationPD, PI, PID Compensation. M. Sami Fadali Professor of Electrical Engineering University of Nevada
PD, PI, PID Compensation M. Sami Fadali Professor of Electrical Engineering University of Nevada 1 Outline PD compensation. PI compensation. PID compensation. 2 PD Control L= loop gain s cl = desired closed-loop
More informationHydrodynamics: Setting the Scene*
Hydrodynamics: Setting the Scene* *a selective view Prof. Rodney Eatock Taylor and Prof. Paul Taylor University of Oxford The Lloyd s Register Educational Trust (LRET) Marine & Offshore Research Workshop
More informationReview of Linear Time-Invariant Network Analysis
D1 APPENDIX D Review of Linear Time-Invariant Network Analysis Consider a network with input x(t) and output y(t) as shown in Figure D-1. If an input x 1 (t) produces an output y 1 (t), and an input x
More informationz Transform System Analysis
z Transform System Analysis Block Diagrams and Transfer Functions Just as with continuous-time systems, discrete-time systems are conveniently described by block diagrams and transfer functions can be
More informationsc Control Systems Design Q.1, Sem.1, Ac. Yr. 2010/11
sc46 - Control Systems Design Q Sem Ac Yr / Mock Exam originally given November 5 9 Notes: Please be reminded that only an A4 paper with formulas may be used during the exam no other material is to be
More informationCITY UNIVERSITY LONDON. MSc in Information Engineering DIGITAL SIGNAL PROCESSING EPM746
No: CITY UNIVERSITY LONDON MSc in Information Engineering DIGITAL SIGNAL PROCESSING EPM746 Date: 19 May 2004 Time: 09:00-11:00 Attempt Three out of FIVE questions, at least One question from PART B PART
More information3 Rigid Spacecraft Attitude Control
1 3 Rigid Spacecraft Attitude Control Consider a rigid spacecraft with body-fixed frame F b with origin O at the mass centre. Let ω denote the angular velocity of F b with respect to an inertial frame
More informationTransform Solutions to LTI Systems Part 3
Transform Solutions to LTI Systems Part 3 Example of second order system solution: Same example with increased damping: k=5 N/m, b=6 Ns/m, F=2 N, m=1 Kg Given x(0) = 0, x (0) = 0, find x(t). The revised
More informationChapter 7. Digital Control Systems
Chapter 7 Digital Control Systems 1 1 Introduction In this chapter, we introduce analysis and design of stability, steady-state error, and transient response for computer-controlled systems. Transfer functions,
More informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Norms for Signals and Systems
. AERO 632: Design of Advance Flight Control System Norms for Signals and. Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Norms for Signals ...
More information10 Transfer Matrix Models
MIT EECS 6.241 (FALL 26) LECTURE NOTES BY A. MEGRETSKI 1 Transfer Matrix Models So far, transfer matrices were introduced for finite order state space LTI models, in which case they serve as an important
More informationSTABILITY. Have looked at modeling dynamic systems using differential equations. and used the Laplace transform to help find step and impulse
SIGNALS AND SYSTEMS: PAPER 3C1 HANDOUT 4. Dr David Corrigan 1. Electronic and Electrical Engineering Dept. corrigad@tcd.ie www.sigmedia.tv STABILITY Have looked at modeling dynamic systems using differential
More informationSubspace-based Identification
of Infinite-dimensional Multivariable Systems from Frequency-response Data Department of Electrical and Electronics Engineering Anadolu University, Eskişehir, Turkey October 12, 2008 Outline 1 2 3 4 Noise-free
More informationIntro to Frequency Domain Design
Intro to Frequency Domain Design MEM 355 Performance Enhancement of Dynamical Systems Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Closed Loop Transfer Functions
More informationNonlinear Tracking Control of Underactuated Surface Vessel
American Control Conference June -. Portland OR USA FrB. Nonlinear Tracking Control of Underactuated Surface Vessel Wenjie Dong and Yi Guo Abstract We consider in this paper the tracking control problem
More informationHomework 7 - Solutions
Homework 7 - Solutions Note: This homework is worth a total of 48 points. 1. Compensators (9 points) For a unity feedback system given below, with G(s) = K s(s + 5)(s + 11) do the following: (c) Find the
More informationMin-Guk Seo, Dong-Min Park, Jae-Hoon Lee, Kyong-Hwan Kim, Yonghwan Kim
International Research Exchange Meeting of Ship and Ocean Engineering in Osaka, December 1-, Osaka, Japan Comparative Study on Added Resistance Computation Min-Guk Seo, Dong-Min Park, Jae-Hoon Lee, Kyong-Hwan
More informationBIBO STABILITY AND ASYMPTOTIC STABILITY
BIBO STABILITY AND ASYMPTOTIC STABILITY FRANCESCO NORI Abstract. In this report with discuss the concepts of bounded-input boundedoutput stability (BIBO) and of Lyapunov stability. Examples are given to
More informationCollocated versus non-collocated control [H04Q7]
Collocated versus non-collocated control [H04Q7] Jan Swevers September 2008 0-0 Contents Some concepts of structural dynamics Collocated versus non-collocated control Summary This lecture is based on parts
More informationSEAKEEPING AND MANEUVERING Prof. Dr. S. Beji 2
SEAKEEPING AND MANEUVERING Prof. Dr. S. Beji 2 Ship Motions Ship motions in a seaway are very complicated but can be broken down into 6-degrees of freedom motions relative to 3 mutually perpendicular axes
More informationANNEX A: ANALYSIS METHODOLOGIES
ANNEX A: ANALYSIS METHODOLOGIES A.1 Introduction Before discussing supplemental damping devices, this annex provides a brief review of the seismic analysis methods used in the optimization algorithms considered
More information