Experimental Robustness Study of a Second-Order Sliding Mode Controller

Size: px
Start display at page:

Download "Experimental Robustness Study of a Second-Order Sliding Mode Controller"

Transcription

1 Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans A.G.e.Jager@wfw.wtb.tue.nl Fax: The esign of a Secon-Orer Sliing Moe Controller is iscusse an guielines are given for tuning. The robustness to unmoele ynamics an parameter-errors is investigate an teste in an experimental case stuy. The experimental results, as far as robustness to unmoele ynamics is concerne, are not better than for a traitional PD-controller. When robustness to parameter-errors is concerne the Secon-Orer Sliing Moe Controller performs slightly better. 1 Introuction Mechanical manipulators are controlle to make their en-effector track a esire trajectory. The control is often base on a mathematical moel that represents the ynamic behavior of the system. In practice, this moel is never an exact representation of reality. There are always phenomena like unmoele ynamics, inaccurate parameters an measurement noise that cause moel imprecisions. These imprecisions may come from actual uncertainty about the system, or from the eliberate choice for a simplifie representation of the system ynamics. The presence of these imprecisions often requires a robust control algorithm. One class of robust controllers is VSS (Variable Structure System), see Utkin (1977). This class of controllers can be use when the moel structure itself is inaccurate, but the inaccuracies are boune with known bouns. VSS controllers are often use as Sliing Moe Controllers. Characteristic components of Sliing Moe Controllers are sliing hyperplanes s(x, t) = 0 in the state space. Recently the traitional sliing control (applying a first-orer sliing conition) has been further enhance, resulting in a Secon-Orer Sliing Moe Control (SOSMC). The SOSMC was presente by Chang (1990) for MIMO systems in Controllability Canonical Form. Elmali an Olgac (1992) extene the SOSMC-technique to general nonlinear MIMO systems, applying I/O-linearization first. In this paper the theory of SOSMC in combination with I/O-linearization is teste in an experimental environment to investigate the robustness to unmoele ynamics (i.e., system of higher orer than the moel) an parametererrors. To tune the SOSMC, a set of guielines is given to achieve best performance. An expression is given for the tracking accuracy an from this expression we recognize a trae-off between tracking accuracy an amping. 2 Preliminaries Consier a nonlinear system that can be escribe (exactly) by a MIMO moel in Controllability Canonical Form, linear in the control u (affine): x (n) = f (x) + B(x)u (1) Publishe in J. Dynamic Systems, Measurement, an Control, Vol. 118, No. 1, March 1996, pp Author to whom all corresponence shoul be sen. 1

2 with x (n) = [y (n 1) 1, y (n 2) 2,...,y (n k) k ] T, n 1 + n 2 + +n k = n with: x R n (state vector) u R m (input vector) y = [y 1,...,y k ] T R k (output vector) Superscripts in parenthesis inicate the orer of time erivatives. Uner certain conitions, a nonlinear system can be transforme into the Controllability Canonical Form by a technique calle I/O-linearization. This technique an the conitions are escribe by Elmali an Olgac (1992), base on Isiori (1989). Assume that to control the system a mathematical moel of the system is available: x (n) = ˆ f (x) + ˆB(x)u (2) For simplicity we assume that the states can be ientifie with those of (1). This is not necessary. Because sliing control requires the uncertainties to be boune with known bouns, a general assumption is: f (x) = f ˆ(x) + f (x) with f (x) α x B(x) = ˆB(x) + B(x) with B(x) β x (3) In properly controlle systems the state vector x will behave boune, so the uncertainty bouns can be etermine. 3 MIMO Sliing Control With Secon-Orer Sliing Conition A secon-orer Sliing Moe Control strategy by Chang (1990) efines a zero z 0, in the error ynamics : ( ) t + z 0 i s i = n i j=1 ( ) t t + λ ji e i τ (4) 0 e i = y i y i for i = 1,...,k This equation is a set of ban-pass filters where the break-frequencies are etermine by the selection of the poles (λ ji ) an zero (z 0i ). An integral term in the equation assures zero steay-state errors. Writing (4) in the unfactore (polynomial) form we get (see Chang (1990)): t ṡ + Z 0 s = e (n 1) + C n 1 e (n 2) + +C 1 e + C 0 e τ (5) where ṡ = 0 represents the sliing hyperplanes. Taking time erivatives of (5) yiels: s + Z 0 ṡ = e (n) + e p (6) with e p = C n 1 e (n 1) + +C 1 ė + C 0 e anthisisusetorelates with the control input u. The error-vector e p can be compute if all states are measure. Substituting (1) in (6) yiels: s + Z 0 ṡ = f + Bu x (n) + e p (7) Stability is guarantee if the control is esigne as: u = ˆB 1 (û k sign(ṡ)) (8) û = ˆ f + x (n) e p + Z 0 ṡ s with = iag(ω 2 n ) 0 2

3 an the gain k is quantifie as (see Chang (1990) an Elmali an Olgac (1992)): α + β ˆB 1 û k > if β ˆB 1 sign(ṡ) < 1 (9) 1 β ˆB 1 sign(ṡ) We get the s-ynamics by substituting the control (8) in (6): s + k sign(ṡ) + s + B ˆB 1 (k sign(ṡ) Z 0 ṡ + s) = f + B ˆB 1 (x (n) ˆ f e p ) (10) As we can see from (8) the control law is iscontinuous across ṡ = 0, which leas to chattering. In general, chattering must be eliminate for the controller to perform properly. This can be achieve by smoothing out the control iscontinuity in a bounary layer neighboring the switching surface (ṡ = 0). Therefore we on t use the signum -function (sign), but we apply the saturation -function (sat) instea for each element ṡ i of ṡ: 1 ṡ i >φ ṡ i sat(ṡ i,φ) = ṡ i φ (11) φ 1 ṡ i < φ Then, the s-ynamics within the bounary layer become: s + k 1 ṡ + s + B ˆB 1 (k 1 ṡ Z 0 ṡ + s) = f + B ˆB 1 (x (n) This equation represents a set of secon-orer low-pass filters. ˆ f e p ) (12) 4 Tuning the SOSMC The lack of tools in nonlinear systems theory now creates a problem; a systematic way of selecting the parameters Z 0, ω n,, C n 1,...,C 0 an k oes not exist. Yet, the influence of the tuning parameters can be evaluate. First, consier the gain k. This controlparameteris completelyetermineby the confinementson the uncertainties an can be calculate with (9). For control esign purposes, the minimum value of k is selecte, since the least control effort is esire. Next, we investigate the influence of the zero z 0 in the error-ynamics. Chang (1990) suggeste that higher values of the zeros provie more amping in the s-ynamics. As we can see from (10/12), z 0 inee contributes in the amping of the s-ynamics, but the exact influence is not clear, since the total amping can increase an ecrease, epening on the sign of the matrix B ˆB 1 Z 0. Simulations, however, inee show more amping when z 0 is increase. This is very likely ue to the error-ynamics (4), in whose response z 0 has a amping contribution, recognize as a reuction of the rise-time. Thir, we investigate the influence of ω n. As we can see from the s-ynamics insie the bounary layer (12), ω n sets the break-frequency. Preferably we choose this break-frequency (etermining the s-ynamics banwith) smaller than the lowest unmoele structural resonant moe. Fourth, to simplify the choice of the poles in the error ynamics (by setting C n 1,..., C 0 ) we can choose the banwith of the error ynamics the same as the banwith (ω n )ofthes-ynamics. This is not necessary. As we will show, a larger error-ynamics-banwith has a positive effect on the tracking error, but may involve a higher control effort. Fifth, we will motivate the choice of, by eriving an expression for the tracking error. Within a finite time ṡ <. Noting that ṡ contains no frequencies higher than ω n (as an approximation) we fin: s <ω n,sothat the maximal tracking error will be (once ṡ < an for the SISO case): ε < φ(ω n + z 0 ) (13) λ 1 λ 2 λ n with a guarantee precision ε (Blom, 1992). has a minimum value, since the implementation has a limite sample frequency. We see from (13) that the zero affects the tracking accuracy in a negative sense also. However, this zero provies amping, which is especially important uring transient. We see that there exists a trae-off between tracking accuracy an amping. 3

4 5 Experiments To investigate the SOSMC in an experimental environment, the control law was implemente in the control software of an XY-table. In Fig. 1 a schematic top view representation of the XY-table is given. The en effector is a slie with mass m e, which can move in the XY-plane by three slieways. Two of them slie in x-irection an one in y-irection. The belt wheels of both slieways are riven by servomotors, exerting torques T 1 an T 3. Coulomb friction appears in all slies an is represente by friction-torques W 1, W 2,anW 3 [Nm]. Viscous amping is represente by D 1, D 2,anD 3 [Nms]. Figure 1: Schematic representation XY-table Unfortunately it is not possible to measure the position of the en effector irectly. Only three encoer signals are available: x 1, x 2,any. We therefore restrict ourselves to the control of motor positions: x 1 an y. The belt wheels of the slieways in x-irections can be connecte in two ways: 1. With a rigi bar (k 1 = ), resulting in a (stiff) moel with two egrees of freeom: x 1, y, since the translations x 1 an x 2 are equal. 2. With a torsion spring with stiffness k 1, resulting in a (flexible) moel with three egrees of freeom: x 1, x 2,any. The equations of motion in case of the two egrees of freeom (stiff) moel are represente by: a 1 ẍ 1 + a 3 sign(ẋ 1 ) + a 5 ẋ 1 = T 1 /r x a 2 ÿ + a 4 sign(ẏ) + a 6 ẏ = T 3 /r y (14) with ientifie parameters: a 1 = J 1 /rx 2 + 2m s + m e + m y = 34 [kg] a 2 = J 3 /ry 2 + m e = 2.7 [kg] a 3 = (W 1 + W 2 )/r x = 36 [N] a 4 = W 3 /r y = 9 [N] a 5 = (D 1 + D 2 )/rx 2 = 50 [Ns/m] a 6 = D 3 /ry 2 = 8 [Ns/m] r x = = 0.01 [m] r y = = 0.01 [m] The equations of motion in case of the three egrees of freeom (flexible) moel are much more complex an will therefore not be presente. 4

5 The esire trajectory to be tracke by the en-effector uring all experiments is chosen to be a circle: x = 0.5 r cos(ωt) y = r sin(ωt) with r = 0.25 [m] an ω = π [ra/s]. The setting of the SMC is tune up, to get the best results, i.e., the controlle system banwith is chosen maximal, with a relative amping in the error ynamics of β = The gain k has been etermine by assuming that the available three egrees of freeom moel is an exact representation of reality, so that the uncertainty bouns on the two egrees of freeom moel can be calculate (3). The control parameters are liste below. c 1 = 63.6 [ra/s] c 0 = 2025 [ra/s 2 ] ω n = 45 [ra/s] z 0 = 90 [ra/s] k = 70 [m/s 2 ] = 1.11 [m/s] To assess robustness of the SOSMC to unmoele ynamics, experiments are one with several torsion springs k 1. The results are compare with a traitional PD-controller, whose setting is also tune up, resulting in a controlle system break frequency of 45 [ra/s]. To eliminate the trae-off between amping an tracking accuracy, experiments are one while moifying the zero on-line, from the initial (i.e., high) value uring transient (amping) to half the initial value (i.e., low) afterwars (high tracking accuracy). If switching is only one once, stability is guarantee. The results for the x-irection only are shown in Fig. 2. In this figure the tracking error RMS is plotte against the stiffness of torsion spring k 1 (8 springs were available with stiffnesses from 0.19 to 213 [Nm/ra]). RMS of tracking error [mm] Spring stiffness [Nm/ra] Figure 2: Robustness to unmoele ynamics; soli: SOSMC, ashe: SOSMC with moification z 0, otte: PD control We see that for all stiffnesses k 1 the SOSMC with moification of z 0 realizes a smaller tracking error than an orinary SOSMC, as expecte (13). For relatively stiff torsion springs (right sie) both SOSMC controllers (with an without moification of z 0 ) realize a smaller tracking error RMS than a PD controller. However, the tracking error RMS for the PD-controller is more constant for a wie range of stiffnesses, in contrast with the SOSMC. For one torsion spring (k 1 = 0.5 [Nm/ra]) the tracking error realize with a PD-controller is even smaller than realize with a SOSMC. We therefore conclue that a PD controller is more robust to unmoele ynamics (ue to stiffness k 1 ) than SOSMC. The level of robustness for both SOSMC is approximately the same. A strange phenomenon is that for a ecreasing stiffness the tracking error RMS increases, but the weakest spring k 1 again yiels a small tracking error RMS. This phenomenon can be seen in all experiments an is because we control the motor positions. 5

6 To assess robustness to parameter variations, experiments are one with aitional mass attache to the eneffector (m e ). The results are shown in Fig. 3. In this figure the tracking error RMS is plotte against the aitional mass. RMS of tracking error [mm] Aitional mass [kg] Figure 3: Robustness to parameter variations; soli: SOSMC, ashe: SOSMC with moification z 0, otte: PD control For all mass-variations the SOSMC with moification of z 0 realizes the smallest, an the PD-controller the highest tracking error. The tracking error, realize with both SOSMC is more constant for variations in the aitional mass, in contrast with the PD-controller. We conclue that both SOSMC are more robust to mass-variations than a PD-controller. Again, the robustness of both SOSMC is approximately the same. 6 Conclusions The main conclusion of this investigation into robustness of the SOSMC is that as far as robustness to unmoele ynamics is concerne, there is no avantage in using a SOSMC technique in favor of a traitional PD-controller. As far as robustness to parameter errors is concerne, the SOSMC performs slightly better than a PD controller. However, for a large range of operating conitions the SOSMC has far better tracking properties, ue to its moel base structure, an shoul be given preference. References Blom, A., 1993, Robustness of a Secon Orer Sliing Moe Controller, Masters thesis, Einhoven University of Technology, Department of Mechanical Engineering, WFW Report Chang, L.-W., 1990, A MIMO Sliing Control with a Secon-Orer Sliing Conition, ASME WAM, paper no. 90 WA/DSC-5, Dallas, Texas. Elmali, H an Olgac, N., 1992, Robust Output Tracking Control of Nonlinear MIMO Systems via Sliing Moe Technique, Automatica, Vol 28, No. 1, pp Isiori, A., 1989, Nonlinear control systems, an introuction, Berlin: Springer-Verlag. Utkin, V. I., 1997, Variable Structure Systems with Sliing Moes, IEEE Transactions on Automatic Control, Vol. AC-22, Apr., pp

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign

More information

Some Remarks on the Boundedness and Convergence Properties of Smooth Sliding Mode Controllers

Some Remarks on the Boundedness and Convergence Properties of Smooth Sliding Mode Controllers International Journal of Automation an Computing 6(2, May 2009, 154-158 DOI: 10.1007/s11633-009-0154-z Some Remarks on the Bouneness an Convergence Properties of Smooth Sliing Moe Controllers Wallace Moreira

More information

Nested Saturation with Guaranteed Real Poles 1

Nested Saturation with Guaranteed Real Poles 1 Neste Saturation with Guarantee Real Poles Eric N Johnson 2 an Suresh K Kannan 3 School of Aerospace Engineering Georgia Institute of Technology, Atlanta, GA 3332 Abstract The global stabilization of asymptotically

More information

On Using Unstable Electrohydraulic Valves for Control

On Using Unstable Electrohydraulic Valves for Control Kailash Krishnaswamy Perry Y. Li Department of Mechanical Engineering, University of Minnesota, 111 Church St. SE, Minneapolis, MN 55455 e-mail: kk,pli @me.umn.eu On Using Unstable Electrohyraulic Valves

More information

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables*

Total Energy Shaping of a Class of Underactuated Port-Hamiltonian Systems using a New Set of Closed-Loop Potential Shape Variables* 51st IEEE Conference on Decision an Control December 1-13 212. Maui Hawaii USA Total Energy Shaping of a Class of Uneractuate Port-Hamiltonian Systems using a New Set of Close-Loop Potential Shape Variables*

More information

Dynamic Load Carrying Capacity of Spatial Cable Suspended Robot: Sliding Mode Control Approach

Dynamic Load Carrying Capacity of Spatial Cable Suspended Robot: Sliding Mode Control Approach Int J Avance Design an Manufacturing echnology, Vol. 5/ No. 3/ June - 212 73 Dynamic Loa Carrying Capacity of Spatial Cable Suspene Robot: Sliing Moe Control Approach M. H. Korayem Department of Mechanical

More information

and from it produce the action integral whose variation we set to zero:

and from it produce the action integral whose variation we set to zero: Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine

More information

Tutorial Test 5 2D welding robot

Tutorial Test 5 2D welding robot Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.

More information

AN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A

AN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A AN INTRODUCTION TO AIRCRAFT WIN FLUTTER Revision A By Tom Irvine Email: tomirvine@aol.com January 8, 000 Introuction Certain aircraft wings have experience violent oscillations uring high spee flight.

More information

THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS

THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS

More information

Sliding mode approach to congestion control in connection-oriented communication networks

Sliding mode approach to congestion control in connection-oriented communication networks JOURNAL OF APPLIED COMPUTER SCIENCE Vol. xx. No xx (200x), pp. xx-xx Sliing moe approach to congestion control in connection-oriente communication networks Anrzej Bartoszewicz, Justyna Żuk Technical University

More information

Role of parameters in the stochastic dynamics of a stick-slip oscillator

Role of parameters in the stochastic dynamics of a stick-slip oscillator Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society

More information

Investigation of local load effect on damping characteristics of synchronous generator using transfer-function block-diagram model

Investigation of local load effect on damping characteristics of synchronous generator using transfer-function block-diagram model ORIGINAL ARTICLE Investigation of local loa effect on amping characteristics of synchronous generator using transfer-function block-iagram moel Pichai Aree Abstract of synchronous generator using transfer-function

More information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Switching Time Optimization in Discretized Hybrid Dynamical Systems Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set

More information

Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity

Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity Preprints of the 9th Worl Congress The International Feeration of Automatic Control Cape Town, South Africa. August -9, Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity Zhengqiang

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

APPPHYS 217 Thursday 8 April 2010

APPPHYS 217 Thursday 8 April 2010 APPPHYS 7 Thursay 8 April A&M example 6: The ouble integrator Consier the motion of a point particle in D with the applie force as a control input This is simply Newton s equation F ma with F u : t q q

More information

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain

Nonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain Nonlinear Aaptive Ship Course Tracking Control Base on Backstepping an Nussbaum Gain Jialu Du, Chen Guo Abstract A nonlinear aaptive controller combining aaptive Backstepping algorithm with Nussbaum gain

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

Electric Power Systems Research

Electric Power Systems Research Electric Power Systems Research 84 (22 35 43 Contents lists available at SciVerse ScienceDirect Electric Power Systems Research jou rn al h om epage: www.elsevier.com/locate/epsr Observer-base nonlinear

More information

Laplacian Cooperative Attitude Control of Multiple Rigid Bodies

Laplacian Cooperative Attitude Control of Multiple Rigid Bodies Laplacian Cooperative Attitue Control of Multiple Rigi Boies Dimos V. Dimarogonas, Panagiotis Tsiotras an Kostas J. Kyriakopoulos Abstract Motivate by the fact that linear controllers can stabilize the

More information

Optimum design of tuned mass damper systems for seismic structures

Optimum design of tuned mass damper systems for seismic structures Earthquake Resistant Engineering Structures VII 175 Optimum esign of tune mass amper systems for seismic structures I. Abulsalam, M. Al-Janabi & M. G. Al-Taweel Department of Civil Engineering, Faculty

More information

A new approach to explicit MPC using self-optimizing control

A new approach to explicit MPC using self-optimizing control 28 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June -3, 28 WeA3.2 A new approach to explicit MPC using self-optimizing control Henrik Manum, Sriharakumar Narasimhan an Sigur

More information

Simple Electromagnetic Motor Model for Torsional Analysis of Variable Speed Drives with an Induction Motor

Simple Electromagnetic Motor Model for Torsional Analysis of Variable Speed Drives with an Induction Motor DOI: 10.24352/UB.OVGU-2017-110 TECHNISCHE MECHANIK, 37, 2-5, (2017), 347-357 submitte: June 15, 2017 Simple Electromagnetic Motor Moel for Torsional Analysis of Variable Spee Drives with an Inuction Motor

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

G j dq i + G j. q i. = a jt. and

G j dq i + G j. q i. = a jt. and Lagrange Multipliers Wenesay, 8 September 011 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

Robust Adaptive Control for a Class of Systems with Deadzone Nonlinearity

Robust Adaptive Control for a Class of Systems with Deadzone Nonlinearity Intelligent Control an Automation, 5, 6, -9 Publishe Online February 5 in SciRes. http://www.scirp.org/journal/ica http://x.oi.org/.436/ica.5.6 Robust Aaptive Control for a Class of Systems with Deazone

More information

PID Adaptive Control Design Based on Singular Perturbation Technique: A Flight Control Example

PID Adaptive Control Design Based on Singular Perturbation Technique: A Flight Control Example PID Aaptive Control Design Base on Singular Perturbation Technique: A Flight Control Example Valery D. Yurkevich, Novosibirsk State Technical University, 20 K. Marx av., Novosibirsk, 630092, Russia (e-mail:

More information

The effect of dissipation on solutions of the complex KdV equation

The effect of dissipation on solutions of the complex KdV equation Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,

More information

Unit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule

Unit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule Unit # - Families of Functions, Taylor Polynomials, l Hopital s Rule Some problems an solutions selecte or aapte from Hughes-Hallett Calculus. Critical Points. Consier the function f) = 54 +. b) a) Fin

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Consider for simplicity a 3rd-order IIR filter with a transfer function. where

Consider for simplicity a 3rd-order IIR filter with a transfer function. where Basic IIR Digital Filter The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient ifference equation

More information

TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM

TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM 265 Asian Journal of Control, Vol. 4, No. 3, pp. 265-273, September 22 TRACKING CONTROL OF MULTIPLE MOBILE ROBOTS: A CASE STUDY OF INTER-ROBOT COLLISION-FREE PROBLEM Jurachart Jongusuk an Tsutomu Mita

More information

Basic IIR Digital Filter Structures

Basic IIR Digital Filter Structures Basic IIR Digital Filter Structures The causal IIR igital filters we are concerne with in this course are characterie by a real rational transfer function of or, equivalently by a constant coefficient

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

State observers and recursive filters in classical feedback control theory

State observers and recursive filters in classical feedback control theory State observers an recursive filters in classical feeback control theory State-feeback control example: secon-orer system Consier the riven secon-orer system q q q u x q x q x x x x Here u coul represent

More information

Systems & Control Letters

Systems & Control Letters Systems & ontrol Letters ( ) ontents lists available at ScienceDirect Systems & ontrol Letters journal homepage: www.elsevier.com/locate/sysconle A converse to the eterministic separation principle Jochen

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

Dynamics of the synchronous machine

Dynamics of the synchronous machine ELEC0047 - Power system ynamics, control an stability Dynamics of the synchronous machine Thierry Van Cutsem t.vancutsem@ulg.ac.be www.montefiore.ulg.ac.be/~vct These slies follow those presente in course

More information

Robust Tracking Control of Robot Manipulator Using Dissipativity Theory

Robust Tracking Control of Robot Manipulator Using Dissipativity Theory Moern Applie Science July 008 Robust racking Control of Robot Manipulator Using Dissipativity heory Hongrui Wang Key Lab of Inustrial Computer Control Engineering of Hebei Province Yanshan University Qinhuangao

More information

arxiv: v1 [physics.flu-dyn] 8 May 2014

arxiv: v1 [physics.flu-dyn] 8 May 2014 Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa

More information

A PID-Sliding Mode Control Design for a Coupled Tank

A PID-Sliding Mode Control Design for a Coupled Tank International Conference CRATT 0, Raes, Tunisia 0 A PID-Sliing Moe Control Design for a Couple Tank Ahme RHIF, Zohra Karous, Naceur BenHaj Braiek Avance System Laboratory, Polytechnic School of Tunisia

More information

Balancing Expected and Worst-Case Utility in Contracting Models with Asymmetric Information and Pooling

Balancing Expected and Worst-Case Utility in Contracting Models with Asymmetric Information and Pooling Balancing Expecte an Worst-Case Utility in Contracting Moels with Asymmetric Information an Pooling R.B.O. erkkamp & W. van en Heuvel & A.P.M. Wagelmans Econometric Institute Report EI2018-01 9th January

More information

Lower Bounds for the Smoothed Number of Pareto optimal Solutions

Lower Bounds for the Smoothed Number of Pareto optimal Solutions Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.

More information

Strength Analysis of CFRP Composite Material Considering Multiple Fracture Modes

Strength Analysis of CFRP Composite Material Considering Multiple Fracture Modes 5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation

More information

Numerical Integrator. Graphics

Numerical Integrator. Graphics 1 Introuction CS229 Dynamics Hanout The question of the week is how owe write a ynamic simulator for particles, rigi boies, or an articulate character such as a human figure?" In their SIGGRPH course notes,

More information

TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS

TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS Francesco Bullo Richar M. Murray Control an Dynamical Systems California Institute of Technology Pasaena, CA 91125 Fax : + 1-818-796-8914 email

More information

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi

BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica

More information

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,

More information

Final Exam Study Guide and Practice Problems Solutions

Final Exam Study Guide and Practice Problems Solutions Final Exam Stuy Guie an Practice Problems Solutions Note: These problems are just some of the types of problems that might appear on the exam. However, to fully prepare for the exam, in aition to making

More information

Open Access An Exponential Reaching Law Sliding Mode Observer for PMSM in Rotating Frame

Open Access An Exponential Reaching Law Sliding Mode Observer for PMSM in Rotating Frame Sen Orers for Reprints to reprints@benthamscience.ae The Open Automation an Control Systems Journal, 25, 7, 599-66 599 Open Access An Exponential Reaching Law Sliing Moe Observer for PMSM in Rotating Frame

More information

State-Space Model for a Multi-Machine System

State-Space Model for a Multi-Machine System State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

Polynomial Inclusion Functions

Polynomial Inclusion Functions Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl

More information

A Path Planning Method Using Cubic Spiral with Curvature Constraint

A Path Planning Method Using Cubic Spiral with Curvature Constraint A Path Planning Metho Using Cubic Spiral with Curvature Constraint Tzu-Chen Liang an Jing-Sin Liu Institute of Information Science 0, Acaemia Sinica, Nankang, Taipei 5, Taiwan, R.O.C., Email: hartree@iis.sinica.eu.tw

More information

A Review of Feedforward Control Approaches in Nanopositioning for High-Speed SPM

A Review of Feedforward Control Approaches in Nanopositioning for High-Speed SPM Garrett M. Clayton Department of Mechanical Engineering, Villanova University, Villanova, PA 1985 e-mail: garrett.clayton@villanova.eu Szuchi Tien Department of Mechanical Engineering, National Cheng Kung

More information

An algebraic expression of stable inversion for nonminimum phase systems and its applications

An algebraic expression of stable inversion for nonminimum phase systems and its applications Proceeings of the 17th Worl Congress The International Feeration of Automatic Control An algebraic expression of stable inversion for nonminimum phase systems an its applications Takuya Sogo Chubu University,

More information

A New Backstepping Sliding Mode Guidance Law Considering Control Loop Dynamics

A New Backstepping Sliding Mode Guidance Law Considering Control Loop Dynamics pp. 9-6 A New Backstepping liing Moe Guiance Law Consiering Control Loop Dynamics V. Behnamgol *, A. Vali an A. Mohammai 3, an 3. Department of Control Engineering, Malek Ashtar University of Technology

More information

The Ehrenfest Theorems

The Ehrenfest Theorems The Ehrenfest Theorems Robert Gilmore Classical Preliminaries A classical system with n egrees of freeom is escribe by n secon orer orinary ifferential equations on the configuration space (n inepenent

More information

An inductance lookup table application for analysis of reluctance stepper motor model

An inductance lookup table application for analysis of reluctance stepper motor model ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR

More information

16.30/31, Fall 2010 Recitation # 1

16.30/31, Fall 2010 Recitation # 1 6./, Fall Recitation # September, In this recitation we consiere the following problem. Given a plant with open-loop transfer function.569s +.5 G p (s) = s +.7s +.97, esign a feeback control system such

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

Power Generation and Distribution via Distributed Coordination Control

Power Generation and Distribution via Distributed Coordination Control Power Generation an Distribution via Distribute Coorination Control Byeong-Yeon Kim, Kwang-Kyo Oh, an Hyo-Sung Ahn arxiv:407.4870v [math.oc] 8 Jul 204 Abstract This paper presents power coorination, power

More information

Interpolated Rigid-Body Motions and Robotics

Interpolated Rigid-Body Motions and Robotics Interpolate Rigi-Boy Motions an Robotics J.M. Selig Faculty of Business, Computing an Info. Management. Lonon South Bank University, Lonon SE AA, U.K. seligjm@lsbu.ac.uk Yaunquing Wu Dept. Mechanical Engineering.

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

Lecture 6: Control of Three-Phase Inverters

Lecture 6: Control of Three-Phase Inverters Yoash Levron The Anrew an Erna Viterbi Faculty of Electrical Engineering, Technion Israel Institute of Technology, Haifa 323, Israel yoashl@ee.technion.ac.il Juri Belikov Department of Computer Systems,

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek

More information

New Simple Controller Tuning Rules for Integrating and Stable or Unstable First Order plus Dead-Time Processes

New Simple Controller Tuning Rules for Integrating and Stable or Unstable First Order plus Dead-Time Processes Proceeings of the 3th WSEAS nternational Conference on SYSTEMS New Simple Controller Tuning Rules for ntegrating an Stable or Unstable First Orer plus Dea-Time Processes.G.ARVANTS Department of Natural

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy IPA Derivatives for Make-to-Stock Prouction-Inventory Systems With Backorers Uner the (Rr) Policy Yihong Fan a Benamin Melame b Yao Zhao c Yorai Wari Abstract This paper aresses Infinitesimal Perturbation

More information

LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR

LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR 5 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR Tatsunori Kaneko* an Yasuto Asano* * Department of Mechanical Engineering,

More information

6. Friction and viscosity in gasses

6. Friction and viscosity in gasses IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner

More information

A Comparison between a Conventional Power System Stabilizer (PSS) and Novel PSS Based on Feedback Linearization Technique

A Comparison between a Conventional Power System Stabilizer (PSS) and Novel PSS Based on Feedback Linearization Technique J. Basic. Appl. Sci. Res., ()9-99,, TextRoa Publication ISSN 9-434 Journal of Basic an Applie Scientific Research www.textroa.com A Comparison between a Conventional Power System Stabilizer (PSS) an Novel

More information

EE 370L Controls Laboratory. Laboratory Exercise #7 Root Locus. Department of Electrical and Computer Engineering University of Nevada, at Las Vegas

EE 370L Controls Laboratory. Laboratory Exercise #7 Root Locus. Department of Electrical and Computer Engineering University of Nevada, at Las Vegas EE 370L Controls Laboratory Laboratory Exercise #7 Root Locus Department of Electrical an Computer Engineering University of Nevaa, at Las Vegas 1. Learning Objectives To emonstrate the concept of error

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information

Chapter 11: Feedback and PID Control Theory

Chapter 11: Feedback and PID Control Theory Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by

More information

Chapter 11: Feedback and PID Control Theory

Chapter 11: Feedback and PID Control Theory Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by

More information

Domain Decomposition and Model Reduction of Systems with Local Nonlinearities

Domain Decomposition and Model Reduction of Systems with Local Nonlinearities Domain Decomposition an Moel Reuction of Systems with Local Nonlinearities Kai Sun, Rolan Glowinsi, Matthias Heinenschloss, an Danny C. Sorensen Abstract The goal of this paper is to combine balance truncation

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

DYNAMIC PERFORMANCE OF RELUCTANCE SYNCHRONOUS MACHINES

DYNAMIC PERFORMANCE OF RELUCTANCE SYNCHRONOUS MACHINES Annals of the University of Craiova, Electrical Engineering series, No 33, 9; ISSN 184-485 7 TH INTERNATIONAL CONFERENCE ON ELECTROMECHANICAL AN POWER SYSTEMS October 8-9, 9 - Iaşi, Romania YNAMIC PERFORMANCE

More information

Adaptive Back-Stepping Control of Automotive Electronic Control Throttle

Adaptive Back-Stepping Control of Automotive Electronic Control Throttle Journal of Software Engineering an Applications, 07, 0, 4-55 http://www.scirp.org/journal/jsea ISSN Online: 945-34 ISSN Print: 945-36 Aaptive Back-Stepping Control of Automotive Electronic Control Throttle

More information

Placement and tuning of resonance dampers on footbridges

Placement and tuning of resonance dampers on footbridges Downloae from orbit.tu.k on: Jan 17, 19 Placement an tuning of resonance ampers on footbriges Krenk, Steen; Brønen, Aners; Kristensen, Aners Publishe in: Footbrige 5 Publication ate: 5 Document Version

More information

Optimal operating strategies for semi-batch reactor used for chromium sludge regeneration process

Optimal operating strategies for semi-batch reactor used for chromium sludge regeneration process Latest Trens in Circuits, Automatic Control an Signal Processing Optimal operating strategies for semi-batch reactor use for chromium sluge regeneration process NOOSAD DAID, MACKŮ LUBOMÍR Tomas Bata University

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

Design of an Industrial Distributed Controller Near Spatial Domain Boundaries

Design of an Industrial Distributed Controller Near Spatial Domain Boundaries Design of an Inustrial Distribute Controller Near Spatial Domain Bounaries Stevo Mijanovic, Greg E. Stewart, Guy A. Dumont, an Michael S. Davies ABSTRACT This paper proposes moifications to an inustrial

More information

Reactive Power Compensation in Mechanical Systems

Reactive Power Compensation in Mechanical Systems Reactive Power Compensation in Mechanical Systems Carlos Rengifo, Bassel Kaar, Yannick Aoustin, Christine Chevallereau To cite this version: Carlos Rengifo, Bassel Kaar, Yannick Aoustin, Christine Chevallereau.

More information

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics

More information

CONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS

CONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS BIT 0006-3835/00/4004-0001 $15.00 200?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN

More information

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate

More information

Hyperbolic Systems of Equations Posed on Erroneous Curved Domains

Hyperbolic Systems of Equations Posed on Erroneous Curved Domains Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (

More information

Vehicle Stability Improvement Based on Electronic Differential Using Sliding Mode Control

Vehicle Stability Improvement Based on Electronic Differential Using Sliding Mode Control 7th WSEAS International Conference on Electric Power Systems, High Voltages, Electric Machines, Venice, Italy, November 1-3, 007 331 Vehicle Stability Improvement Base on Electronic Differential Using

More information

Adjoint Transient Sensitivity Analysis in Circuit Simulation

Adjoint Transient Sensitivity Analysis in Circuit Simulation Ajoint Transient Sensitivity Analysis in Circuit Simulation Z. Ilievski 1, H. Xu 1, A. Verhoeven 1, E.J.W. ter Maten 1,2, W.H.A. Schilers 1,2 an R.M.M. Mattheij 1 1 Technische Universiteit Einhoven; e-mail:

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information