An experiment on using the Pioneer3-DX robot in a multichannel measurement system
|
|
- Laura Elfreda Garrison
- 5 years ago
- Views:
Transcription
1 8th WSEAS Internatonal Conference on SYSTEMS THEORY and SCIENTIFIC COMUTATION (ISTASC 8 Rhodes, Greece, August -, 8 An experment on usng the oneer3-dx robot n a multchannel measurement system G. VASILIU *, A. FILIESCU *,AL.STANCU *, S. FILIESCU ** * Department of Automaton and Industral Informatcs ** Student n Electrcal Engneerng Unversty Dunarea de Jos of Galat Domneasca, 7, Galat, 88, ROMANIA Abstract: Increasng the precson of multchannel measurement systems requres a mnmum executon tme to select and analyze many probes, accessble from a depost located nsde the measurement system. Moble robot oneer3 handles all the movements of selected probes between the depost and measurement stand. Sldng-mode control of the robotc manpulator, and fuzzy technques for trajectory trackng, are used for a mnmum manpulaton tme and hgher securty n selecton and placement of probes. The dynamcal model of the moble robot wth one par of actve wheels, tme varyng mass and moment of nerta have been used n sldng-mode control. Two closed-loop, on-lne parameter estmators have been desgned to acheve robustness aganst parameter uncertantes. Two sldng-mode adaptve controllers correspondng to angular and translaton moton have been desgned. Keywords: Multchannel measurement system, oneer3-dx, Real-tme fuzzy control, Wheeled moble robot, Sldng-mode control.. Introducton to autoous, ntellgent, multchannel measurement systems Multchannel measurement systems are an mportant component of ntegrated steel manufacturng process. The nput of these systems are the physcal probes and the output s the nformaton concernng percentual concentraton of varous chemcal component. For a relatvely long perod of tme (days, weeks, such system can operate wthout the nterventon of a human operator. Usually, the overall precson of an ntellgent measurement system s ncreased f the result s corrected after comparng t wth several control probes. These control probes are stored n a structured probe-holder, nsde the measurement system. Fgure shows how the oneer 3-DX robot grasps the selected probe wth hs grpper. Fg.. oneer 3-DX grpper wth a steel probe. Fgure shows the same robot selectng the probe from a set of control probes. After the selecton, the robot carres the probes to the analsys stand and back to the probe holder. The probe holder s structured as a lbrary and every probe s assocated wth an ID code, a spatal poston and a n-dmensonal vector representng the composton. ISSN: ISBN:
2 8th WSEAS Internatonal Conference on SYSTEMS THEORY and SCIENTIFIC COMUTATION (ISTASC 8 Rhodes, Greece, August -, 8 Fg.. oneer 3-DX wheeled moble robot wth DOF robotc manpulator selectng a probe from the control set The avalable space n an autoous measurement system must be effcently used. Fgure 3 llustrates a compact block of probes, organzed as parallel matrx, every par of matrx are avalable n an moble stand D k. The optmal trajectory mnmzes the tme needed for the auxlary measurements of control probes. The precson of measurement for the current probe ncreases when the measurement s wthn a short tme nterval (necessary n order to approxmate the state of the system as stable wth a great number of auxlary probes, havng the composton confrmed by several laboratores. Once the destnaton ponts of robot trajectory s determned, the mnmum tme to transport all the selected probes can be obtaned usng the correspondng values for parameters lke: Dstance between the moble blocks wth probes; Velocty of the robot body; Angular velocty of the robotc arm; All these parameters must have correspondng numercal values adapted to the varable mass of the robot (when the robot grasps a new probe, the total mass ncreases, and when the probes are deposted back n the probe holder the mass decreases (the mass of a steel probe s.. kg. A smplfed dagram of movement problem s presented n fgure. COMACT BLOCK OF ROBES COMACT BLOCK OF ROBES Fg. 3. A compact block of probes organzed n parallel matrx. The ntal space between the moble stands D and D 3, after the dsplacement of stands D 3 D k becomes fnally the space between the moble stands D k and D k+. The unque space between all the stands of the probe holder and the moblty of all stands, represent the practcal soluton to keep a large number of probes n the lmted space of the measurement system. The space between two stands allows the access of the robot oneer 3-DX for handlng the selected control probe. In ths experment, fuzzy logc and sldng mode control s used to solve the problem of optmal trajectory of robot between the blocks wth measurement probes. IONEER3-DX ROBOT BODY ARM AND GRIER Fg.. Smplfed dsplacement of robot between the blocks of the probe holder. ISSN: ISBN:
3 8th WSEAS Internatonal Conference on SYSTEMS THEORY and SCIENTIFIC COMUTATION (ISTASC 8 Rhodes, Greece, August -, 8 The soluton of suspended blocks wth probes allows a substantal decrease (around procents of the necessary space between the moble blocks of probes. (In fgure, the dmensons of the robot body s double when ncludng the dmenson of the arm and grpper Actve zone of grpper ARM AND GRIER ROBOT BODY Fg.. Smplfed dagram of the robot body, arm, and actve zone of the grpper.. Fuzzy control of the robot For every possble process stuaton, the trajectory of the robot between the measurement stand and the locaton of the probe n the probe holder s determned, and can be represented as a lst of spatal coordnates (x,y. The optmum trajectory s assumed to be the one tracked n the least tme nterval. At any gven moment, the poston error of the current pont (x,y to the segment of the generalzed trajectory (Ax+by+C s determned by the dstance: Ax + By + C d ( A + B The outputs must be controlled so that the dstance from the current poston to the nearest pont of generalzed trajectory s mnmum. The error dot correspondng to the error derved from ( s: Δe e' ( Δt provded that Δt s small. Ths leads to three obvous domans of varaton for e(t: negatve (N, zero (Z and postve (. Note that, all the three membershp functons assocated wth e(t are fully determned by the parameter M, and the degree of membershp (DOM to a specfc doman can be easly calculated by lnear nterpolaton. By defnng a smlar set of membershp functons for the error dot (e (tδe/δt, one gets a total of 9 rules (3x3 of the followng type: If the error s postve and the error dot s postve, then v R must be HIGH and v L must be LOW. The entre rule base descrbng the fuzzy controller s presented n Table. Each cell of table contans a logc sentence and should be read as: If e(t s Negatve AND e (t s negatve, THEN V L must be HIGH and V R must be LOW Table. Rule base for the Fuzzy controller error dot e (t error e(t N Z N HL LH LM Z HL MM MH ML MH LH where H,M,L desgnate the sngleton values for HIGH, MEDIUM and LOW fuzzy domans of the outputs v L, v R respectvely. The truth value of the antecedent of the above sentence s mn(n(e,n(e, where N(e, N(e are the degrees of membershp of e(t and e (t to the doman N. The crsp output of the fuzzy controller s a combnaton of all the rules n the rule base as follows: v where: K out K z z S z mn(e, E' (3 ( ISSN: ISBN:
4 8th WSEAS Internatonal Conference on SYSTEMS THEORY and SCIENTIFIC COMUTATION (ISTASC 8 Rhodes, Greece, August -, 8 S s the correspondng sngleton value of the fuzzy output, and K s the total number of rules n the rule base. E, E are the degrees of membershp of e(t and e (t to the doman correspondng to the cell. Also generates the values of the left and rght wheel speed, v L and v R, for the next tme nterval untl the entre parabolc segment s fnshed. 3. Dynamc model and parameter uncertantes The vehcle dynamcs s fully descrbed by a three dmensonal vector of generalzed coordnates q ( t consttuted by the coordnates (( x ( t, y( t of the mdpont between the two drvng wheels, and by the orentaton angle Φ( t. The velocty constrant (non-holoc constrant of vehcle moton s x & sn Φ y& cos Φ. Defne by τ r and τ l the torques provded by DC motors to the rght and left wheel, respectvely. x x, x y, x3 Φ x x&, x y&, x Φ& (7 Defne the control nput correspondng to angular, u A τ r τ l and poston moton, u τ r + τ l, respectvely. x x + Tx x x( k + Tx x3 x3 + Tx x x Tx x + Tπ cos( x u (8 3 x x + Tx x + Tπ sn( x3 u x x + Tα ( ku A ( x x x x sn 3 cos( 3 (9 The control nput for angular moton has two terms: the frst one, denoted compensaton part u comp A (k, has to compensate the rotatonal dynamcs; the second one, denoted sldng mode part,, corresponds to system evoluton nsde of u sm A ( k sldng surface neghborhood. The whole control nput for angular moton s u comp sm u u A A + A ( Fg.. WMR confguraton varables for angular and poston moton. The vehcle s descrbed by the followng dynamcal model(see fgure where m, I, D, r are the robot mass, moment of nerta, dstance between wheels and wheels radus, respectvely τ r + τ l mx && my& Φ & + cos Φ r τ r + τ l my && mx& Φ & + sn Φ ( r D IΦ && ( τ r τ l r The real values of the parameters are tme-varyng wth upper bounded uncertantes α π Let real real max ( t α Δα ( t ; Δα Δα ( max ( t π Δπ ( t ; Δπ Δπ x R be the state vector, whose elements are Expressng the estmated value for angular moton control nput parameter, ˆ α α Δ ˆ α, the next sequence, correspondng to recursve least squares method, (Ljung, 999; Stoca and Ahgren,, can be used to provde an estmaton of the uncertanty scalar term Δ α ( k at the k th step Δ α ( k u A( k LΔα ( + u k k [ ( ] ( A Δα ( k L u ( k ( k Δ Δα Δα A Δα α ( ( k ( k L ( k u ( k ( k Δα Δα Δα A Δα (3 u k, sm u k u k. For ˆ π k π Δ ˆ π k, The control nput for poston moton, ( has only sldng-mode part, ( ( the correspondng parameter, ( ( smlar updatng law s used L Δ π ( k + Δ π ( k u ( k [ u ( k ] ( k Δ π ( k L u ( k ( k ( Δ π Δπ Δπ π Δπ ( ISSN: ISBN:
5 8th WSEAS Internatonal Conference on SYSTEMS THEORY and SCIENTIFIC COMUTATION (ISTASC 8 Rhodes, Greece, August -, 8 ( k u ( k u ( k S TΔ ˆ π Δ ˆ π Δ ˆ π( k + LΔπ + π ( + S where L Δ π, π have the same meanng as prevously and wll be defned later. Δ S. Sldng-mode adaptve control of angular moton The followng stable sldng surface has been chosen, n order to desgn the control nput for angular moton S A A( k + μ A (7 where ( ( ( ( x k δ e k A k x3 k arctg x ( ( k δe k (8 wth: μ (, δ, δ. arameter μ T and the poston errors, e, e, establsh the dynamcs of sldng surface. The nterval set of δ and δ assures the stablty of poston errors. For computng the control nput, the followng attractveness condton, (Furuta, 99; Yu and Xu,, has been used: S A ΔS A ( k + < ΔS A ( k + (9 where ΔS A ( k + S A ( k + S A ( The followng sldng surface s proposed: S ( k [ x ( k ] + [ x ( k ] ( [ x ( k δ e ( k ] + [ x ( k δ e ( k ] ( Startng wth the thrd equaton of model (, usng a trgoetrc equalty and the non-holoc constrant (, the followng equalty holds x x x ( ( x x tgtx k + ( x x x Moreover, ntroducng the expressons of the state varables, from state model (8, and usng the constrant(9, the above equalty becomes tg ( + Tπ u ([ x k ] + [ x ] ( Tx [ x ] [ x ] Tx ( Let defne S ( k cos ( Tx ( k (3 [ ] ([ x ( k ] + [ x ( k ] [ x ( k + δ e ( k ] + [ x ( k + δ e ( k ] ( The sldng moton on the surface concerns the reduced order system of the robotc model, wthout of 3 rd and th equaton. The same attractveness condton, for computng the poston moton control nput has been consdered S ΔS ( k + < ΔS ( k + ( ΔS ( k + S ( k + S (k ( An approxmate sldng mode evoluton on the surface can be assured. Consequently of sldngmode evoluton, the angular state x 3 tends to hold the followng expressons costx k x k δ e k ( 3( ( ( ( [ x δe ( k ] + [ x δe ( k ] ( Tx3 ( x δe ( k [ x δe ( k ] + [ x δe ( k ] sn (7 (8 Usng (, the followng expresson can be obtaned [ x ( k + ] + [ x ( k + ] [ ( Tx ] cos ([ ( ] [ ( ] ( ( ( x k + x k T π Δπ k u k Wth (7 and (, (3 and ( become S T Δ ( k + S ( k cos( Tx π Δπ ( k u S ( k + S ( k S T cos( Tx ( k π Δπ ( k ( ( k ( u ( k (9 (3 (3 Usng (9, (3 and upper bound of poston moton uncertanty, the second degree nequalty can be wrtten ISSN: ISBN:
6 8th WSEAS Internatonal Conference on SYSTEMS THEORY and SCIENTIFIC COMUTATION (ISTASC 8 Rhodes, Greece, August -, 8 If max ( π Δπ T u + S [ S ] < cos( Tx (3 > and S > S, then the sldng u sm control nput for poston moton s u ρ T where ( S S max [ cos( Tx ] ( π Δπ (33 ρ. The control nput stll can be computed usng onlne estmates for Δ π. Theore, the approxmate sldng mode condton s satsfed, S ( k +,.e. T cos Tx k π Δ ˆ π k u k + S k [ ( ( ] ( ( ( ( (3 From above, the control nput can be expressed as u ( k S ( k T cos( Tx ( k ( [ ] π ˆ π ( k Δ (3 When the system evolves n sldng-mode, can express the followngs x k x k δ e k (3 x ( ( ( x e ( k δ (37 Theore, output trackng error dynamcs assocated to the reduced order system can be expressed as: e ( k + e Te ( k ( k + e ( k Te ( k δ (38 e δ (39 For, the above errors are stable. δ, δ T For testng the proposed dscrete-tme sldngmode adaptve controller oneer 3-DX wth on board C and wreless adapter has been used n crcular trajectory trackng. The two motors use 38.3: gear ratos and contan -tck encoders. The followng parameters of model (3 were used: mkg, Dcm, I, kgm, T.3s.The moment of nerta has been computed assumng the mass unformly dstrbuted. A lnear-tme varyng mass addtonally to the nal one has been consdered from kg to kg. y(t 3 poston trajectory - 3 x(t 8 erence robot angular trajectory erence robot 8 t [sec] - control nput - 8 t [sec] trackng error angular poston poston angular - 8 t [sec] Fg.7. WMR closed loop response for crcular erence and ntal condtons x(33; x(33; x3( π /7; x(-.; x(.; x(.. The crcle trajectory trackng, shown n fgures 7, was obtaned for Δα max., Δπ max.33. The followng values have been chosen for the constants: μ., ρ ρ.99, δ δ 3. 33, ( ( A. Concluson Δ α Δπ. Dscrete-tme, sldng-mode adaptve controllers and parameter estmators for trajectory trackng appled to control angular and poston moton of oneer 3-DX one par of actve wheels moble robot, have presented n ths paper. The tmevaryng mass and moment of nerta dynamcal state space model have been undertaken n order to desgn the controllers. Two sldng-mode adaptve controllers have been desgned, for angular and poston moton, respectvely. The robustness s guaranteed by sldng-mode controllers and by onlne parameter estmators. Controllers parameters, on-lne updated, assure an approxmate sldngmode evoluton even f the attractveness condton s not satsfed and contrbute to an ncreased robustness. References [] Ljung, L., (999 System Identfcaton ( nd edton. rentce Upper saddle Rver, NJ. ISSN: ISBN:
7 8th WSEAS Internatonal Conference on SYSTEMS THEORY and SCIENTIFIC COMUTATION (ISTASC 8 Rhodes, Greece, August -, 8 [] Stoca,. and Ahgren., ( Exact ntalzaton of the recursve least squares algorthm. Internatonal Journal of Adaptve Control and Sgnal rocessng, vol., pp [3] Furuta, K., (99, Sldng mode control of a dscrete system. System & Control Letters, vol., pp.-. [] Yu, X., Xu, J.,X., ( Varable Structures Systems: Towards st Century, Sprnger_Verlag Berln Hedelberg. ISSN: ISBN:
On-line Parameter Estimation in Sliding-mode Control of Pioneer 3-DX Wheeled Mobile Robot
roceedngs of the 7th WSES Internatonal Conference on Systems Theory and Scentfc Computaton, thens, Greece, ugust -, 7 7 On-lne arameter Estmaton n Sldng-mode Control of oneer 3-DX Wheeled Moble Robot.
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationAn Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors
An Algorthm to Solve the Inverse Knematcs Problem of a Robotc Manpulator Based on Rotaton Vectors Mohamad Z. Al-az*, Mazn Z. Othman**, and Baker B. Al-Bahr* *AL-Nahran Unversty, Computer Eng. Dep., Baghdad,
More informationWeek 9 Chapter 10 Section 1-5
Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,
More informationModeling of Dynamic Systems
Modelng of Dynamc Systems Ref: Control System Engneerng Norman Nse : Chapters & 3 Chapter objectves : Revew the Laplace transform Learn how to fnd a mathematcal model, called a transfer functon Learn how
More information829. An adaptive method for inertia force identification in cantilever under moving mass
89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationStudy on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI
2017 2nd Internatonal Conference on Electrcal and Electroncs: echnques and Applcatons (EEA 2017) ISBN: 978-1-60595-416-5 Study on Actve Mcro-vbraton Isolaton System wth Lnear Motor Actuator Gong-yu PAN,
More informationDesign and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm
Desgn and Optmzaton of Fuzzy Controller for Inverse Pendulum System Usng Genetc Algorthm H. Mehraban A. Ashoor Unversty of Tehran Unversty of Tehran h.mehraban@ece.ut.ac.r a.ashoor@ece.ut.ac.r Abstract:
More informationTime-Varying Systems and Computations Lecture 6
Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy
More information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationTracking with Kalman Filter
Trackng wth Kalman Flter Scott T. Acton Vrgna Image and Vdeo Analyss (VIVA), Charles L. Brown Department of Electrcal and Computer Engneerng Department of Bomedcal Engneerng Unversty of Vrgna, Charlottesvlle,
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationThe Chaotic Robot Prediction by Neuro Fuzzy Algorithm (2) = θ (3) = ω. Asin. A v. Mana Tarjoman, Shaghayegh Zarei
The Chaotc Robot Predcton by Neuro Fuzzy Algorthm Mana Tarjoman, Shaghayegh Zare Abstract In ths paper an applcaton of the adaptve neurofuzzy nference system has been ntroduced to predct the behavor of
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationSolutions to exam in SF1811 Optimization, Jan 14, 2015
Solutons to exam n SF8 Optmzaton, Jan 4, 25 3 3 O------O -4 \ / \ / The network: \/ where all lnks go from left to rght. /\ / \ / \ 6 O------O -5 2 4.(a) Let x = ( x 3, x 4, x 23, x 24 ) T, where the varable
More information4DVAR, according to the name, is a four-dimensional variational method.
4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The
More informationIterative General Dynamic Model for Serial-Link Manipulators
EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators 1. Introducton Iteratve General Dynamc Model for Seral-Lnk Manpulators In ths set of notes, we are gong to develop a method for computng a general
More informationUncertainty as the Overlap of Alternate Conditional Distributions
Uncertanty as the Overlap of Alternate Condtonal Dstrbutons Olena Babak and Clayton V. Deutsch Centre for Computatonal Geostatstcs Department of Cvl & Envronmental Engneerng Unversty of Alberta An mportant
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 31 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 6. Rdge regresson The OLSE s the best lnear unbased
More information10.34 Numerical Methods Applied to Chemical Engineering Fall Homework #3: Systems of Nonlinear Equations and Optimization
10.34 Numercal Methods Appled to Chemcal Engneerng Fall 2015 Homework #3: Systems of Nonlnear Equatons and Optmzaton Problem 1 (30 ponts). A (homogeneous) azeotrope s a composton of a multcomponent mxture
More information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationAn Interactive Optimisation Tool for Allocation Problems
An Interactve Optmsaton ool for Allocaton Problems Fredr Bonäs, Joam Westerlund and apo Westerlund Process Desgn Laboratory, Faculty of echnology, Åbo Aadem Unversty, uru 20500, Fnland hs paper presents
More informationUncertainty in measurements of power and energy on power networks
Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:
More informationThe Study of Teaching-learning-based Optimization Algorithm
Advanced Scence and Technology Letters Vol. (AST 06), pp.05- http://dx.do.org/0.57/astl.06. The Study of Teachng-learnng-based Optmzaton Algorthm u Sun, Yan fu, Lele Kong, Haolang Q,, Helongang Insttute
More informationChapter 11 Angular Momentum
Chapter 11 Angular Momentum Analyss Model: Nonsolated System (Angular Momentum) Angular Momentum of a Rotatng Rgd Object Analyss Model: Isolated System (Angular Momentum) Angular Momentum of a Partcle
More informationChapter Newton s Method
Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationPart C Dynamics and Statics of Rigid Body. Chapter 5 Rotation of a Rigid Body About a Fixed Axis
Part C Dynamcs and Statcs of Rgd Body Chapter 5 Rotaton of a Rgd Body About a Fxed Axs 5.. Rotatonal Varables 5.. Rotaton wth Constant Angular Acceleraton 5.3. Knetc Energy of Rotaton, Rotatonal Inerta
More informationArmy Ants Tunneling for Classical Simulations
Electronc Supplementary Materal (ESI) for Chemcal Scence. Ths journal s The Royal Socety of Chemstry 2014 electronc supplementary nformaton (ESI) for Chemcal Scence Army Ants Tunnelng for Classcal Smulatons
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationIndeterminate pin-jointed frames (trusses)
Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all
More informationDO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.
EE 539 Homeworks Sprng 08 Updated: Tuesday, Aprl 7, 08 DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. For full credt, show all work. Some problems requre hand calculatons.
More informationHarmonic Detection Algorithm based on DQ Axis with Fourier Analysis for Hybrid Power Filters
Harmonc Detecton Algorthm based on DQ Axs wth Fourer Analyss for Hybrd Power Flters K-L. AREERAK Power Qualty Research Unt, School of Electrcal Engneerng Insttute of Engneerng, Suranaree Unversty of Technology
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationAdditional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty
Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,
More informationIrregular vibrations in multi-mass discrete-continuous systems torsionally deformed
(2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected
More informationPhysics 111: Mechanics Lecture 11
Physcs 111: Mechancs Lecture 11 Bn Chen NJIT Physcs Department Textbook Chapter 10: Dynamcs of Rotatonal Moton q 10.1 Torque q 10. Torque and Angular Acceleraton for a Rgd Body q 10.3 Rgd-Body Rotaton
More informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More informationSimulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests
Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth
More informationMoments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.
Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these
More informationDesign and Analysis of Landing Gear Mechanic Structure for the Mine Rescue Carrier Robot
Sensors & Transducers 214 by IFSA Publshng, S. L. http://www.sensorsportal.com Desgn and Analyss of Landng Gear Mechanc Structure for the Mne Rescue Carrer Robot We Juan, Wu Ja-Long X an Unversty of Scence
More informationχ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
More informationˆ (0.10 m) E ( N m /C ) 36 ˆj ( j C m)
7.. = = 3 = 4 = 5. The electrc feld s constant everywhere between the plates. Ths s ndcated by the electrc feld vectors, whch are all the same length and n the same drecton. 7.5. Model: The dstances to
More informationAdiabatic Sorption of Ammonia-Water System and Depicting in p-t-x Diagram
Adabatc Sorpton of Ammona-Water System and Depctng n p-t-x Dagram J. POSPISIL, Z. SKALA Faculty of Mechancal Engneerng Brno Unversty of Technology Techncka 2, Brno 61669 CZECH REPUBLIC Abstract: - Absorpton
More informationOutline. Communication. Bellman Ford Algorithm. Bellman Ford Example. Bellman Ford Shortest Path [1]
DYNAMIC SHORTEST PATH SEARCH AND SYNCHRONIZED TASK SWITCHING Jay Wagenpfel, Adran Trachte 2 Outlne Shortest Communcaton Path Searchng Bellmann Ford algorthm Algorthm for dynamc case Modfcatons to our algorthm
More informationHigh resolution entropy stable scheme for shallow water equations
Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal
More informationLecture Note 3. Eshelby s Inclusion II
ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte
More information1 Derivation of Point-to-Plane Minimization
1 Dervaton of Pont-to-Plane Mnmzaton Consder the Chen-Medon (pont-to-plane) framework for ICP. Assume we have a collecton of ponts (p, q ) wth normals n. We want to determne the optmal rotaton and translaton
More information11. Dynamics in Rotating Frames of Reference
Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons
More informationStructure and Drive Paul A. Jensen Copyright July 20, 2003
Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationcoordinates. Then, the position vectors are described by
Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,
More informationResearch on the Fuzzy Control for Vehicle Semi-active Suspension. Xiaoming Hu 1, a, Wanli Li 1,b
Advanced Materals Research Onlne: 0-0- ISSN: -9, Vol., pp -9 do:0.0/www.scentfc.net/amr.. 0 Trans Tech Publcatons, Swterland Research on the Fuy Control for Vehcle Sem-actve Suspenson Xaomng Hu, a, Wanl
More informationPivot-Wheel Drive Crab with a Twist! Clem McKown Team November-2009 (eq 1 edited 29-March-2010)
Pvot-Wheel Drve Crab wth a Twst! Clem McKown Team 1640 13-November-2009 (eq 1 edted 29-March-2010) 4-Wheel Independent Pvot-Wheel Drve descrbes a 4wd drve-tran n whch each of the (4) wheels are ndependently
More informationModelling of the precise movement of a ship at slow speed to minimize the trajectory deviation risk
Computatonal Methods and Expermental Measurements XIV 29 Modellng of the precse movement of a shp at slow speed to mnmze the trajectory devaton rsk J. Maleck Polsh Naval Academy, Poland Faculty of Mechancs
More informationMath1110 (Spring 2009) Prelim 3 - Solutions
Math 1110 (Sprng 2009) Solutons to Prelm 3 (04/21/2009) 1 Queston 1. (16 ponts) Short answer. Math1110 (Sprng 2009) Prelm 3 - Solutons x a 1 (a) (4 ponts) Please evaluate lm, where a and b are postve numbers.
More informationSo far: simple (planar) geometries
Physcs 06 ecture 5 Torque and Angular Momentum as Vectors SJ 7thEd.: Chap. to 3 Rotatonal quanttes as vectors Cross product Torque epressed as a vector Angular momentum defned Angular momentum as a vector
More informationNMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582
NMT EE 589 & UNM ME 48/58 ROBOT ENGINEERING Dr. Stephen Bruder NMT EE 589 & UNM ME 48/58 7. Robot Dynamcs 7.5 The Equatons of Moton Gven that we wsh to fnd the path q(t (n jont space) whch mnmzes the energy
More informationMEV442 Introduction to Robotics Module 2. Dr. Santhakumar Mohan Assistant Professor Mechanical Engineering National Institute of Technology Calicut
MEV442 Introducton to Robotcs Module 2 Dr. Santhakumar Mohan Assstant Professor Mechancal Engneerng Natonal Insttute of Technology Calcut Jacobans: Veloctes and statc forces Introducton Notaton for tme-varyng
More information= 1.23 m/s 2 [W] Required: t. Solution:!t = = 17 m/s [W]! m/s [W] (two extra digits carried) = 2.1 m/s [W]
Secton 1.3: Acceleraton Tutoral 1 Practce, page 24 1. Gven: 0 m/s; 15.0 m/s [S]; t 12.5 s Requred: Analyss: a av v t v f v t a v av f v t 15.0 m/s [S] 0 m/s 12.5 s 15.0 m/s [S] 12.5 s 1.20 m/s 2 [S] Statement:
More informationNotes on Analytical Dynamics
Notes on Analytcal Dynamcs Jan Peters & Mchael Mstry October 7, 004 Newtonan Mechancs Basc Asssumptons and Newtons Laws Lonely pontmasses wth postve mass Newtons st: Constant velocty v n an nertal frame
More informationElshaboury SM et al.; Sch. J. Phys. Math. Stat., 2015; Vol-2; Issue-2B (Mar-May); pp
Elshabour SM et al.; Sch. J. Phs. Math. Stat. 5; Vol-; Issue-B (Mar-Ma); pp-69-75 Scholars Journal of Phscs Mathematcs Statstcs Sch. J. Phs. Math. Stat. 5; (B):69-75 Scholars Academc Scentfc Publshers
More informationFormulas for the Determinant
page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use
More informationEN40: Dynamics and Vibrations. Homework 7: Rigid Body Kinematics
N40: ynamcs and Vbratons Homewor 7: Rgd Body Knematcs School of ngneerng Brown Unversty 1. In the fgure below, bar AB rotates counterclocwse at 4 rad/s. What are the angular veloctes of bars BC and C?.
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationConvexity preserving interpolation by splines of arbitrary degree
Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete
More informationHongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)
ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of
More informationLecture 12: Classification
Lecture : Classfcaton g Dscrmnant functons g The optmal Bayes classfer g Quadratc classfers g Eucldean and Mahalanobs metrcs g K Nearest Neghbor Classfers Intellgent Sensor Systems Rcardo Guterrez-Osuna
More informationMATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018
MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng
More information8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS
SECTION 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS 493 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS All the vector spaces you have studed thus far n the text are real vector spaces because the scalars
More informationSecond Order Analysis
Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to
More informationParameter Estimation for Dynamic System using Unscented Kalman filter
Parameter Estmaton for Dynamc System usng Unscented Kalman flter Jhoon Seung 1,a, Amr Atya F. 2,b, Alexander G.Parlos 3,c, and Klto Chong 1,4,d* 1 Dvson of Electroncs Engneerng, Chonbuk Natonal Unversty,
More informationSpin-rotation coupling of the angularly accelerated rigid body
Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s
More informationLagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013
Lagrange Multplers Monday, 5 September 013 Sometmes t s convenent to use redundant coordnates, and to effect the varaton of the acton consstent wth the constrants va the method of Lagrange undetermned
More informationMathematical Preparations
1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the
More informationCollege of Computer & Information Science Fall 2009 Northeastern University 20 October 2009
College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:
More informationLectures on Multivariable Feedback Control
Lectures on Multvarable Feedback Control Al Karmpour Department of Electrcal Engneerng, Faculty of Engneerng, Ferdows Unversty of Mashhad June 200) Chapter 9: Quanttatve feedback theory Lecture Notes of
More informationAmplification and Relaxation of Electron Spin Polarization in Semiconductor Devices
Amplfcaton and Relaxaton of Electron Spn Polarzaton n Semconductor Devces Yury V. Pershn and Vladmr Prvman Center for Quantum Devce Technology, Clarkson Unversty, Potsdam, New York 13699-570, USA Spn Relaxaton
More informationVARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES
VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES BÂRZĂ, Slvu Faculty of Mathematcs-Informatcs Spru Haret Unversty barza_slvu@yahoo.com Abstract Ths paper wants to contnue
More informationModule 3: Element Properties Lecture 1: Natural Coordinates
Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers
More informationStatistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA
14 th Internatonal Users Conference Sesson: ALE-FSI Statstcal Energy Analyss for Hgh Frequency Acoustc Analyss wth Zhe Cu 1, Yun Huang 1, Mhamed Soul 2, Tayeb Zeguar 3 1 Lvermore Software Technology Corporaton
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationRobert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations
Quantum Physcs 量 理 Robert Esberg Second edton CH 09 Multelectron atoms ground states and x-ray exctatons 9-01 By gong through the procedure ndcated n the text, develop the tme-ndependent Schroednger equaton
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationNote 10. Modeling and Simulation of Dynamic Systems
Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada
More informationCS 468 Lecture 16: Isometry Invariance and Spectral Techniques
CS 468 Lecture 16: Isometry Invarance and Spectral Technques Justn Solomon Scrbe: Evan Gawlk Introducton. In geometry processng, t s often desrable to characterze the shape of an object n a manner that
More informationChapter 11: Simple Linear Regression and Correlation
Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests
More informationSOLVING CAPACITATED VEHICLE ROUTING PROBLEMS WITH TIME WINDOWS BY GOAL PROGRAMMING APPROACH
Proceedngs of IICMA 2013 Research Topc, pp. xx-xx. SOLVIG CAPACITATED VEHICLE ROUTIG PROBLEMS WITH TIME WIDOWS BY GOAL PROGRAMMIG APPROACH ATMII DHORURI 1, EMIUGROHO RATA SARI 2, AD DWI LESTARI 3 1Department
More informationFUZZY FINITE ELEMENT METHOD
FUZZY FINITE ELEMENT METHOD RELIABILITY TRUCTURE ANALYI UING PROBABILITY 3.. Maxmum Normal tress Internal force s the shear force, V has a magntude equal to the load P and bendng moment, M. Bendng moments
More information6.3.4 Modified Euler s method of integration
6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from
More informationA Tale of Friction Basic Rollercoaster Physics. Fahrenheit Rollercoaster, Hershey, PA max height = 121 ft max speed = 58 mph
A Tale o Frcton Basc Rollercoaster Physcs Fahrenhet Rollercoaster, Hershey, PA max heght = 11 t max speed = 58 mph PLAY PLAY PLAY PLAY Rotatonal Movement Knematcs Smlar to how lnear velocty s dened, angular
More informationGeneralized Linear Methods
Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set
More information