Modeling the Dynamics of an Ice Tank Carriage

Similar documents
d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

RC, RL and RLC circuits

Electrical and current self-induction

6.01: Introduction to EECS I Lecture 8 March 29, 2011

INVERSE RESPONSE COMPENSATION BY ESTIMATING PARAMETERS OF A PROCESS COMPRISING OF TWO FIRST ORDER SYSTEMS

Vehicle Arrival Models : Headway

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Some Basic Information about M-S-D Systems

STATE-SPACE MODELLING. A mass balance across the tank gives:

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

A First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18

Notes on Kalman Filtering

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

SUPPLEMENTARY INFORMATION

Sliding Mode Controller for Unstable Systems

L1, L2, N1 N2. + Vout. C out. Figure 2.1.1: Flyback converter

Computation of the Effect of Space Harmonics on Starting Process of Induction Motors Using TSFEM

Sub Module 2.6. Measurement of transient temperature

1. VELOCITY AND ACCELERATION

2.4 Cuk converter example

Determination of the Sampling Period Required for a Fast Dynamic Response of DC-Motors

Lab 10: RC, RL, and RLC Circuits

04. Kinetics of a second order reaction

Hall effect. Formulae :- 1) Hall coefficient RH = cm / Coulumb. 2) Magnetic induction BY 2

ECE 2100 Circuit Analysis

Principle and Analysis of a Novel Linear Synchronous Motor with Half-Wave Rectified Self Excitation

Two Coupled Oscillators / Normal Modes

i L = VT L (16.34) 918a i D v OUT i L v C V - S 1 FIGURE A switched power supply circuit with diode and a switch.

Embedded Systems and Software. A Simple Introduction to Embedded Control Systems (PID Control)

Chapter 2: Principles of steady-state converter analysis

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

V AK (t) I T (t) I TRM. V AK( full area) (t) t t 1 Axial turn-on. Switching losses for Phase Control and Bi- Directionally Controlled Thyristors

Learning Objectives: Practice designing and simulating digital circuits including flip flops Experience state machine design procedure

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing

R.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder#

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

The equation to any straight line can be expressed in the form:

Numerical Dispersion

Class Meeting # 10: Introduction to the Wave Equation

Dead-time Induced Oscillations in Inverter-fed Induction Motor Drives

Chapter 7 Response of First-order RL and RC Circuits

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

INDEX. Transient analysis 1 Initial Conditions 1

Estimation of Poses with Particle Filters

I. Define the Situation

ES 250 Practice Final Exam

Chapter 4. Truncation Errors

Topic Astable Circuits. Recall that an astable circuit has two unstable states;

Analysis and design of a high-efficiency zero-voltage-switching step-up DC DC converter

V L. DT s D T s t. Figure 1: Buck-boost converter: inductor current i(t) in the continuous conduction mode.

Chapter Floating Point Representation

Module 4: Time Response of discrete time systems Lecture Note 2

NEWTON S SECOND LAW OF MOTION

Time Domain Transfer Function of the Induction Motor

Electrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit

3. Alternating Current

(b) (a) (d) (c) (e) Figure 10-N1. (f) Solution:

Optimal Control of Dc Motor Using Performance Index of Energy

AC Circuits AC Circuit with only R AC circuit with only L AC circuit with only C AC circuit with LRC phasors Resonance Transformers

ECE 2100 Circuit Analysis

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits

Silicon Controlled Rectifiers UNIT-1

dv 7. Voltage-current relationship can be obtained by integrating both sides of i = C :

LAB 5: Computer Simulation of RLC Circuit Response using PSpice

Linear Response Theory: The connection between QFT and experiments

Signal and System (Chapter 3. Continuous-Time Systems)

not to be republished NCERT MATHEMATICAL MODELLING Appendix 2 A.2.1 Introduction A.2.2 Why Mathematical Modelling?

The motions of the celt on a horizontal plane with viscous friction

Lab #2: Kinematics in 1-Dimension

Structural Dynamics and Earthquake Engineering

CHAPTER 12 DIRECT CURRENT CIRCUITS

2.9 Modeling: Electric Circuits

Unsteady Flow Problems

Notes 04 largely plagiarized by %khc

Unified Control Strategy Covering CCM and DCM for a Synchronous Buck Converter

Lecture 20: Riccati Equations and Least Squares Feedback Control

Homework-8(1) P8.3-1, 3, 8, 10, 17, 21, 24, 28,29 P8.4-1, 2, 5

Matlab and Python programming: how to get started

Involute Gear Tooth Bending Stress Analysis

CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK

Maintenance Models. Prof. Robert C. Leachman IEOR 130, Methods of Manufacturing Improvement Spring, 2011

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model

The problem with linear regulators

Technical Report Doc ID: TR March-2013 (Last revision: 23-February-2016) On formulating quadratic functions in optimization models.

Chapter 10 INDUCTANCE Recommended Problems:

MATHEMATICAL MODELING OF THE TRACTOR-GRADER AGRICULTURAL SYSTEM CINEMATIC DURING LAND IMPROVING WORKS

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL

Problem Set #1. i z. the complex propagation constant. For the characteristic impedance:

Detection of Tire Lateral Force Based on a Resolver Mechanism

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

EE100 Lab 3 Experiment Guide: RC Circuits

Physics for Scientists & Engineers 2

Analyze patterns and relationships. 3. Generate two numerical patterns using AC

Dual Current-Mode Control for Single-Switch Two-Output Switching Power Converters

Section 3.8, Mechanical and Electrical Vibrations

Section 4.4 Logarithmic Properties

UNIVERSITY OF CALIFORNIA AT BERKELEY

EECS 2602 Winter Laboratory 3 Fourier series, Fourier transform and Bode Plots in MATLAB

Transcription:

Modeling he Dynamics of an Ice Tank Carriage The challenge: To model he dynamics of an Ice Tank Carriage and idenify a mechanism o alleviae he backlash inheren in he design of he gearbox. Maplesof, a division of Waerloo Maple Inc., 28 Execuive Summary Inroducion Problem Saemen 1. Model Descripion of Ice Tank Carriage 2. Verifying Model Accuracy 3. Improving Model Dynamics Resuls ** This applicaion was developed using Maple and MapleSim www.maplesof.com/appsbriefs Page 1 of 9

Execuive Summary Tradiional ships due o heir design can no be used o sail hrough ice-covered waers such as hose presen in he Arcic. Maneuvering hrough ice-covered waers requires a special design of ship commonly known as an icebreaker. An icebreaker has hree disinguishing feaures ha enable i o crush he ice beneah i. These are: 1) a srenghened hull, 2) an ice-clearing shape, and 3) a powerful engine. Before an icebreaker is consruced, he design of each of hese hree feaures is rigorously esed and opimized in an ice ank. An ice ank is a refrigeraed waer ank in which a scale model of he ship's hull is pulled by a ganry carriage hrough a layer of formed ice. The ganry carriage is driven by an elecric moor hrough a gearbox ha is elecronically conrolled o obain he required orque and speed performance for he scale model. Tha said, problems inheren in he design of he ice ank could poenially disor he inegriy of any es resuls ha are obained. This was he case, when i was noed ha he speed of he model ship became unsable when he ganry carriage was operaing a low speeds. Invesigaion ino his problem led o he source of he error; i was found ha he force generaed by he ship s hull when breaking ice would creae backlash beween he gear eeh driving he ganry carriage. The conroller, in an aemp o correc for he backlash, would cause he speed of he model ship o overshoo and hen undershoo, gradually increasing in magniude o a poin ha would cause evenual damage o he equipmen. To recify he error, a high-fideliy mahemaical model of he ice ank sysem was creaed using MapleSim and Maple. The model was simulaed using wo differen parameer ses; he firs parameer se simulaed he model wihou perurbaions from breaking he ice, while he second se included he effecs of breaking he ice. This was done o ensure ha he model accuraely capured he vehicle dynamics refleced in he real sysem. To alleviae he backlash in he gearbox, differen soluions were esed on he model. Re-calibraing he PID conroller reduced he backlash presen in he sysem, however i was no able o remove i o a level ha was deemed saisfacory by he eam. As an alernaive, a orque preload componen was added o he oupu shaf of he gearbox as a means o preven he gear eeh from slipping. Using he Opimizaion Templae in MapleSim, i was found ha applying a orque preload of -3.5 Nm was sufficien o sabilize he speed of he ganry carriage and hus he model ship. Inroducion An icebreaker is used o sail hrough ice-covered waers because is design enables i o break and clear ice along is pah. An icebreaker, which ypically moves a speeds ranging from.2 m/s o 4. m/s, is able o crush he ice beneah i by driving is bow up ono he ice and using is massive weigh o break he ice. Before consrucion of an icebreaker begins, i is rigorously esed in an ice ank. An ice ank consiss of a refrigeraed waer basin, which provides a physical modeling environmen o measure he vehicle dynamics of ship's hull as i is pulled hrough a layer of ice. The scale model of he ship's hull is pulled along he ice ank via a suspended ganry mechanism, which is driven by an elecric moor and gearbox assembly. The composiion of he formed ice is modified o reflec he ice/waer condiions of he environmen where he icebreaker will operae. The ice ank presen in our faciliy measures 9m by 12m by 3m. To ge accurae ice ank simulaion resuls, i is of umos imporance ha any problems inheren in he design of he ice ank be correced. Tha said, wihou an accurae model of he ice ank mechanism, disinguishing he source of error beween he ice ank or he ship's design can become very difficul. As a resul, i is www.maplesof.com/appsbriefs Page 2 of 9

common pracice in indusry o develop a mahemaical model of he ice ank which, when simulaed accuraely predics he performance seen in he ice ank basin. Problem Saemen This documen presens he seps involved in developing an accurae model of he ice ank carriage. In addiion, i illusraes he approach ha was aken o eliminae chaer in he gearbox, due o backlash beween he gear eeh. 1. Model Descripion of Ice Tank Carriage 1.1 Free-Body Diagram The free-body diagram shown in Figure 1 is a simplified represenaion of a model ship in an ice ank. As shown, he force of he waer F waer and he force from he impac of breaking he ice F ice are opposing he force driving he ganry carriage F drive. Alhough i is no depiced, i is imporan o noe ha he moor/drive conains a gearbox. The gearbox ranslaes he moor orque ino a linear force ha pulls he ganry carriage. In he acual physical sysem, he rail is locaed along he op of he ice ank basin (ha is, he inerface beween he ice and he model ship). In his diagram, he rail is locaed below he Moor/Drive componen for he purposes of illusraing he ineracion of he forces in he sysem. Figure 1: Free-Body diagram of Ice Tank Mechanism www.maplesof.com/appsbriefs Page 3 of 9

Physical Model Represenaion The free body diagram in Figure 1 can be easily convered ino physical model represenaion using MapleSim. The MapleSim diagram depiced in Figure 2 is a physical model represenaion of he sysem. The orque sep componen in his represenaion was added o eliminae backlash in he gear during ice breaking. The reason for adding he orque sep componen is discussed in he nex secion. As can be seen in he diagram, he physical model is made up of a plan and conroller submodels. The conroller submodel adjuss he volage level driving he moor in response o changes (or errors) in he speed of he model ship. These changes are relayed back o he PID conroller hrough a speed sensor ha is conneced beween he ganry carriage and model ship submodels. www.maplesof.com/appsbriefs Page 4 of 9

The plan model consiss of 5 submodels. These are: 1) DC Moor, 2) Gear, 3) Ganry Carriage, 4) Model Ship, and 5) Breaking Ice. The DC Moor submodel was modeled using a simplified model of a DC Moor as an equivalen elecric circui The Gear submodel was modeled wih he buil-in siffness, damping and dead-space necessary o represen he backlash in he gear eeh The Ganry submodel was represened as an ideal roaional o ranslaion gear The Model Ship was represened as a sliding mass of approximaely 8, kg, wih a resising force o he damping effec of waer The Breaking Ice submodel was included in he model as a means o measure he sysem response due o breaking ice. This submodel was modeled using a simple single saw-ooh cycle: a linearly increasing force ha releases o zero when he ice breaks. Each of hese submodels, and he componens ha were used o model hem can be seen in he.msim file. www.maplesof.com/appsbriefs Page 5 of 9

2. Verifying Model Accuracy The model in Figure 2 (wihou he preload orque componen) was simulaed using wo differen parameer ses. This was done o es he validiy and accuracy of he model agains real-world resuls. The firs parameer se simulaed he ice ank model wihou any perurbaions due o breaking ice, whereas he second parameer se simulaed he ice ank model aking ino accoun he effecs of breaking ice. The resuls of hese wo simulaions are shown in he nex wo secions, respecively. Sysem response wihou ice breaking perurbaions As can be seen from he four plos below, when he sysem was no perurbed (ha is, no breaking ice), he speed of he ganry carriage followed a smooh rajecory. In addiion, he error beween he sep inpu and he velociy feedback signal quickly wen o zero. In fac, he only reason ha here is an error is because i akes he sysem approximaely 35s o reach seady sae. 12 8 6 4 2 Probe1.GearTorque 2 4 6 8 Probe3.SignalError.2.15.1.5 Probe2.GanrySpeed.2.15.1.5 2 4 6 8 'Moor::Probe4'.MoorVolage 2 15 5 2 4 6 8 2 4 6 8 Noe: These plos were obained by simulaing he model wih he parameer se defined in "IniialParamSe_NoPerubaions.params" file. www.maplesof.com/appsbriefs Page 6 of 9

Sysem response wih ice breaking perurbaions As expeced, when he sysem was ordered o break ice, he speed of he ganry carriage began o oscillae. The sysem response o breaking ice can be seen in he four plos below. Annoaions on he plo describe he sysem response a differen ime seps. Probe1.GearTorque 12 Torque due o resisance of ice 8 6 4 2 K2 Sudden release when ice breaks 2 4 6 8 Probe3.SignalError.2.2.15.1.5 Probe2.GanrySpeed Induced oscillaions in speed 2 4 6 8 'Moor::Probe4'.MoorVolage 2.15.1.5 Increasing error signal 2 4 6 8 15 5 2 4 6 8 Noe: These plos were obained by simulaing he model wih he parameer se defined in "IniialParamSe_Perubaions.params" file. www.maplesof.com/appsbriefs Page 7 of 9

3. Improving Model Dynamics Many differen approaches were considered o remove he induced speed oscillaions in he ganry carriage. One approach ha was aken involved increasing he proporional gain of he PID conroller. Unforunaely, alhough increasing he gain reduced he chaer beween he gear eeh, i also amplified he volage supplied o he DC moor o a level ha was significanly larger han he raing of he moor. A common pracice used in indusry o remove gear backlash is o apply a preload orque o he gearbox. This forces he gear eeh o remain in conac a all imes. We were able o model he effecs of a preload orque by adding an "applied orque" block o he physical model. Deermining he amoun of "applied orque" needed o sabilize he speed of he ganry carriage was accomplished using he Opimizaion emplae (please refer o he Opimizaion.mw file locaed in he Documen Folder for more informaion). A suiable speed response was obained by using a orque load of -3.5Nm. The sysem response o a preload orque is shown in he following plos. Clearly, we can see ha he gear orque no longer shows signs of chaer, and ha he speed of he ganry carriage sabilizes o a desired speed. 12 8 6 4 2 Probe2.GearTorque Same orque disurbance due o breaking of ice 2 4 6 8.2.15.1.5 Probe1.GanrySpeed Speed disurbance due o breaking ice quickly sabilizes o desired speed 2 4 6 8 Probe3.SignalError.2 'Moor::Probe4'.MoorVolage 2.15.1.5 Error due o breaking of ice sabilizes o desired speed 15 5 2 4 6 8 2 4 6 8 Noe: These plos were obained by simulaing he model wih he parameer se defined in "FinalParamSe.params" file. www.maplesof.com/appsbriefs Page 8 of 9

Resuls This documen oulines he seps aken o develop a high-fideliy physical model of he ice ank carriage ha accuraely depics he behavior during esing. The model was validaed by simulaing i using wo differen parameer ses. The firs parameer se simulaed he model wihou perurbaions from breaking he ice, while he second se included he effecs of breaking he ice. As expeced, when he sysem was se o break ice, he speed of he ganry carriage, and hus he model carriage, began o oscillae. Invesigaion ino his problem led o he source of he error. I was found ha he force generaed by he ship s hull when breaking ice would creae backlash beween he gear eeh driving he ganry carriage. The conroller, in an aemp o correc for he backlash, would cause he speed of he model ship o overshoo and hen undershoo. The addiion of a orque preload componen was considered o alleviae he backlash presen in he gearbox. Modeling he orque preload componen was achieved by using a orque sep componen. Using he Opimizaion Templae in MapleSim, i was found ha applying a orque preload of -3.5Nm was sufficien o sabilize he speed of he ganry carriage and, as a resul, he model ship. Legal Noice: Maplesof, a division of Waerloo Maple Inc. 28. Maplesof, Maple, and MapleSim are rademarks of Waerloo Maple Inc. All oher rademarks are he propery of heir respecive owners. This applicaion may conain errors and Maplesof is no liable for any damages resuling from he use of his maerial. www.maplesof.com/appsbriefs Page 9 of 9