Control of industrial robots

Size: px
Start display at page:

Download "Control of industrial robots"

Transcription

1 Control o industrial robots Control with eternal sensors Pro. Paolo Roo (paolo.roo@polimi.it) Politenio di Milano Dipartimento di Elettronia, Inormazione e Bioingegneria

2 Introdution So ar we have assumed that the robot is ontrolled in a ree environment However the robot has to interat with the environment either to manipulate it, or to avoid to ollide with it, or to interat with other devies or humans In order to eetively interat with the environment, the robot must be endowed with eternal (eteroeptive) sensors, whih make the robot ontroller aware o the situation Common sensors used to measure suh interation are ore sensors and vision sensors In the ollowing we will study with some detail the ontrol systems based on ore measurements. We will also give a short outline on vision based ontrol systems Control o industrial robots Control with eternal sensors P. Roo [2]

3 Interation with the environment A robot ommonly interats with the working environment, manipulating objets and perorming operations on suraes (think o grinding, deburring, part assembly). In modern robotis, situations are ommon where the robot physially interats with the operator (think o telemanipulation or rehabilitation). Control o industrial robots Control with eternal sensors P. Roo [3]

4 Passive ontrol o ompliane A irst way to address the interation with the environment is to endow the manipulator with devies that ailitate the eeution o the task in a passive way: RCC: Remote Center o Compliane The RCC is used in assembly tasks (peg-in-a-hole) It is typially plaed between the robot s wrist and the gripper The RCC lets the gripper assembly move in the plane perpendiular to the peg s ais. This allows the peg to rotate With the RCC, the ores generated by any misalignment an be ompensated in a passive way Control o industrial robots Control with eternal sensors P. Roo [4]

5 Ative ontrol o the interation During more general interation tasks, the use o a purely positional ontrol strategy (the same adopted in ree motion) may lead to problems due to positioning errors and unertainties in task planning, related to a inomplete knowledge o the environment. Eamples where an ative ontrol o the interation is desirable inlude: mehanial mahining (deburring, surae inishing, polishing, assembly,...) tele-manipulation ooperation o multi-manipulator systems deterous robot hands physial human-robot interation Eperiments in ore/position ontrol with a rigid surae at POLIMI. For an ative ontrol o interation it is neessary to use devies or ore (and moment) measurement. Control o industrial robots Control with eternal sensors P. Roo [5]

6 Fore sensors Fore/torque sensors are devies that return the measurements o three omponents o ore and three omponents o moment with respet to a loal rame. Control o industrial robots Control with eternal sensors P. Roo [6]

7 Internal behaviour Sensors are based on strain gauges, i.e. devies that an measure a strain. The strain indues a variation in the resistane o a iruit. Strain gauges are suitably mounted, in suh a way as to return inormation rom whih three omponents o ore and moment an be omputed. Control o industrial robots Control with eternal sensors P. Roo [7]

8 Taonomy Fore based interation ontrol strategies an be divided into two ategories: Impedane/admittane ontrol The goal is to assign a presribed dynami relation between interation ores and position errors We want that the manipulator, in ase o unplanned interations, behaves like a generalized mass-spring-damper system Hybrid position and ore ontrol We separate diretions whih are only position ontrolled rom diretions whih are only ore/torque ontrolled In the diretions onstrained by the environment we want that a ore/torque with a speiied value be established Furthermore we will talk about: Epliit ontrol The ore (or the impedane) is diretly assigned, ating on the ontrol variables (joint torques) Impliit ontrol The ore (or the impedane) is indiretly assigned, ating on the set-points o the position ontrol loops Control o industrial robots Control with eternal sensors P. Roo [8]

9 Mehanial impedane The low o power among two physial systems an always be deined as the produt o two onjugate quantities: a (generalized) eort and a (generalized) low In the eletrial ramework, the low is the urrent, while the eort is the voltage. In the mehanial ramework, instead, the low is the (linear or angular) veloity, the eort is the ore (or the moment) The mehanial impedane is then deined as dynamial relation that is established between ore and veloity (or displaement) or a mehanial system. The admittane is just the reiproal o the impedane The impedane ontrol aims at making the manipulator, position ontrolled and in interation with the environment, assume a desired mehanial impedane, like a generalized mass-spring-damper system Control o industrial robots Control with eternal sensors P. Roo [9]

10 Mehanial impedane An eample o impedane/admittane ontrol or physial human-robot interation: Control o industrial robots Control with eternal sensors P. Roo [10]

11 A 1 d.o.. system In order to understand how an impedane ontroller operates, let us irst onsider the simple eample o a body with mass M interating with the environment, whih eerts a ore : u M The system an be modeled with the blok diagram: u + 1 M.. 1 s. 1 s The model o the interation with the environment is not speiied Control o industrial robots Control with eternal sensors P. Roo [11]

12 A 1 d.o.. system Let us onsider an inverse dynamis motion ontrol law: u = My with: y = && + K ~ & + K d D P ~ where: ~ = d The losed loop dynamis are obtained: M ~ && + D ~ & d + Kd ~ = where: D = MK, K = d D d MK P The equation establishes a mehanial impedane relation between the ore ating on the system and the positional error. However the inertial omponent o the impedane (M) annot be arbitrarily assigned. Control o industrial robots Control with eternal sensors P. Roo [12]

13 A 1 d.o.. system Suppose now that the ore is measured and let: u = My + with: y = ( ) M && + D ~ + K ~ 1 & d d d M d.. d d This orresponds to the blok diagram: The losed loop dynamis are:. d d + + D d K d M d + y M + + u 1 M.. 1 s. 1 s M ~ & d + Dd ~ & + Kd ~ = It is an impedane relation where all the terms an be arbitrarily assigned. Control o industrial robots Control with eternal sensors P. Roo [13]

14 Impliit impedane ontrol (admittane ontrol) When a positional ontrol is already applied to the mehanial system, it might be onvenient to lose the impedane ontrol outside the position ontrol. Consider again the ase o single d.o.. system, position ontrolled: d R u M The system an be modeled with the blok diagram: + u + d R(s) G(s) ( ) 1 G s = 2 Ms Control o industrial robots Control with eternal sensors P. Roo [14]

15 Impliit impedane ontrol (admittane ontrol) Consider the ollowing blok diagram, where the ore applied by the environment is measured and iltered through a system with transer untion G d (s), to give a variation o the reerene d : d + + R(s) u + G(s) G d (s) impedane ontroller Elaborating the blok diagram: = G s 1+ R s ( ) ( ) G( s) + ( G ( s) ) d d ( ) G( s) ( ) G( s) R s 1+ R s Also alled admittane ontroller (it aepts a ore and delivers a displaement) Control o industrial robots Control with eternal sensors P. Roo [15]

16 Impliit impedane ontrol (admittane ontrol) In a band o requenies well inside the bandwidth o the position ontrol system, we an assume the ore disturbane well rejeted and the position reerene orretly ollowed. Thus: d Gd ( s) = G ( s) d 1~ I we impose: G d ( s) = M d s D d s + K d M ~ && d + Dd ~ & + Kd ~ = we assign to the manipulator the desired impedane behaviour. The advantage o this sheme is that the position ontrol is kept unhanged: the impedane ontrol modiies the set-points o the position ontrol, whih is usually the only openness allowed by an industrial motion ontrol system. Control o industrial robots Control with eternal sensors P. Roo [16]

17 Statis o the manipulator Beore applying the onepts o impedane ontrol to a omplete manipulator, we need to inlude the eet o a ore applied at the end-eetor in the model. Suppose then that a ore is applied at the end-eetor. What are the joint torques that keep the system in equilibrium???? Control o industrial robots Control with eternal sensors P. Roo [17]

18 Systems o ores As a matter o at a system o ores may at at the end-eetor. How an this system be redued? From lassial mehanis One ore applied in a point (the resultant) One torque It ollows that we will always be able to represent a system o ores with two vetors: A ore vetor applied in a point O (the resultant) A moment vetor (the moment o the original system o ores with respet to the ulrum O) Control o industrial robots Control with eternal sensors P. Roo [18]

19 Statis o the manipulator Let us then onsider a manipulator subjet to a set o ores at the end eetor, having a resultant and moment, w.r.t. the origin o the end eetor rame, µ. We then deine vetor: h = µ Let us all τ the vetor o torques (or ores) applied to the joints. The virtual work o the ores/torques at the joints is: T δw τ = τ δq The virtual work o the ores at the end-eetor is: T T δw = δp + µ ωdt = J µ h where δq is the elementary displaement o the joint oordinates. T P T T ( q) δq + J ( q) δq = h J( q) δq where δp is the virtual displaement o the origin o the end-eetor rame, onsistent with the dierential kinematis o the robot, and ω is the angular veloity o the end eetor rame. O Control o industrial robots Control with eternal sensors P. Roo [19]

20 Statis o the manipulator T δw τ = τ δq T δw = J µ h P T T ( q) δq + J ( q) δq = h J( q) δq, O J ( q) = J J P O ( q) ( ) q position Jaobian orientation Jaobian Based on the priniple o virtual work the manipulator is in stati equilibrium i and only i: δw τ = δw h τ = ( q)h J T Notie that the same matri, the Jaobian matri, whih deines the dierential kinematis o the manipulator, is also involved in the statis o the system (kinetostati duality). Control o industrial robots Control with eternal sensors P. Roo [20]

21 Inverse dynamis ontrol The dynami model in interation with the environment is then given by the equation: ( )& T + C( q, q& ) q& + g( q) = τ J ( q)h B q q where h is the vetor o the ontat ores (and moments) applied by the end eetor on the environment. Let us onsider now the adoption o an inverse dynamis ontrol law: τ ( q) y + C( q, q& ) q g( q) = B & + Substituting we obtain: q&& = y B 1 T ( q) J ( q)h nonlinear oupling term Control o industrial robots Control with eternal sensors P. Roo [21]

22 Inverse dynamis ontrol Assume now or y the epression introdued in the inverse dynami ontrol in the operational spae: y By substituting the epression or y we obtain the equation o the dynamis in losed loop: && ~ where = d, while: B A 1 A = J ( ) q& ( q) & + K ~ & + K ~ J& ( q, q& ) & d T 1 ( q) J ( q) B( q) J ( q) D ~ K P ~ 1 K D ~ + + = B A = A A P ( q) h A is the inertia matri in the operational spae and: h A = T T A ( ) h, J = T A ( ) J A A We have obtained an impedane relation whih is oupled and only partially assignable: a ore applied along a Cartesian diretion yields motion in other diretions as well. T J h = J generalized ores perorming work on ẋ Control o industrial robots Control with eternal sensors P. Roo [22] T A h A

23 Impedane ontrol with ore measurement I ore measurements are available, we an set: τ T ( ) y + C( q, q& ) q + g( q) + J ( q)h = B q & and: y = J 1 A ( ) q& h 1 ( q) M M & + D ~ & + K ~ M J& ( q, q& ) d d d d d d A A where M d, D d, K d are diagonal, positive deinite matries. We obtain: M ~ && + D ~ & + K ~ = d d d h A i.e. a ompletely deoupled system. Control o industrial robots Control with eternal sensors P. Roo [23]

24 Impedane ontrol The equation: M ~ && + D ~ & + K ~ = d d d h A deines a mehanial impedane relation between the ores/moments deined at the right hand side o the equations and the position/orientation errors in the operational spae. This impedane has the orm o a mehanial system with mass matri M d, damping D d, and stiness K d. The dynami behaviour along eah diretion o the operational spae an thus be speiied: a ore or a moment along a ertain Cartesian diretion yields a motion in suh diretion, whih is related to the ore and moment by a ully ontrollable dynamis. But what do we eatly mean or impedane in the operational spae? Control o industrial robots Control with eternal sensors P. Roo [24]

25 Translational impedane The irst three degrees o reedom in the operational spae are the position oordinate o the end-eetor. In this ase the deinition o the impedane is rather intuitive: p d : desired position p e : atual position p = p d p e The impedane relation is: M p&& + D p& + K p = p p p We thus dupliate on eah degree o reedom what we have already seen or the single mass. Control o industrial robots Control with eternal sensors P. Roo [25]

26 Rotational impedane The remaining three degrees o reedom in the operational spae onsist in a minimal representation (Euler angles) o the end-eetor orientation. In this ase the deinition o the impedane is less intuitive: φ d : desired orientation φ e : atual orientation φ = φ d φ e The impedane relation is: M && φ + D φ& + K φ = T φ φ φ Τ ( φ )µ e Due to the presene o matri T, this relation is subjet to representation singularities. Better deinition o the rotational impedane have been proposed in the literature. Control o industrial robots Control with eternal sensors P. Roo [26]

27 Natural onstraints and artiiial onstraints An alternative way to perorm an interation task is to assign reerene values o the ores and o the positions, onsistent with the geometry o the environment. Some diretions an be in at subjeted to onstraints on the position, other ones to onstraints on the ores that an be applied. In general, it is not possible to impose both the ore, and the position, along the same diretion. We make the ollowing distintion: Natural onstraints: Artiiial onstraints: they are imposed by the environment along eah degree o reedom o the task; they depend on the geometry o the task. onstraints imposed by the ontrol system, related to the task eeution strategy. For task desription we reer to a onstraint rame, y, z. Control o industrial robots Control with eternal sensors P. Roo [27]

28 Sliding on a lat surae Natural Constraints p& ω ω z y y µ z Artiiial Constraints µ µ z p& p& ω y y z This piture is taken rom the tetbook: B. Siiliano, L. Siavio, L. Villani, G. Oriolo: Robotis: Modelling, Planning and Control, 3rd Ed. Springer, 2009 Control o industrial robots Control with eternal sensors P. Roo [28]

29 Insertion o a peg in a hole Natural Constraints p& p& ω ω y y z Artiiial Constraints µ µ y p& y z This piture is taken rom the tetbook: µ z. Usually we set p z 0 and all the other artiiial onstraints to zero. ω z B. Siiliano, L. Siavio, L. Villani, G. Oriolo: Robotis: Modelling, Planning and Control, 3rd Ed. Springer, 2009 Control o industrial robots Control with eternal sensors P. Roo [29]

30 Rotation o a rank Oriented like the rank Natural Constraints p& p& ω z Artiiial Constraints µ z Idle knob ω y y µ z µ p& ω y y z This piture is taken rom the tetbook: B. Siiliano, L. Siavio, L. Villani, G. Oriolo: Robotis: Modelling, Planning and Control, 3rd Ed. Springer, 2009 Normally artiiial onstraints on ores and moments are assumed to be zero. Control o industrial robots Control with eternal sensors P. Roo [30]

31 Seletion matri We assume that the onstraints an always be epressed along mutually orthogonal diretions o rame C. Thus on eah diretion either a natural onstraint on position, or a natural onstraint on ore, eists. In order to speiy the onstraints we use a seletion matri Σ diagonal and with null or unitary elements: Natural onstraints Artiiial onstraints n ( I ) h h Σv = v Σ = ( I Σ) v = v Σh = h a 0 0 Σ = n a The seletion matri etrats diretions where a natural onstraint in veloity applies and thus an artiiial onstraint in ore an be enored Control o industrial robots Control with eternal sensors P. Roo [31]

32 Hybrid ore/position ontrol Hybrid ore/position ontrol is a ontrol strategy that uses the artiiial onstraints to speiy the goals o the ontrol system, in terms o both ore and position. The 6 diretions o the operational spae are divided into position/veloity ontrolled diretions and ore/moment ontrolled diretions: Σ h A R h Ad v d Σ I Σ INTERACTION CONTROL MOTION CONTROL y F y v + + y y u h A R NONLINEAR MANIPULATOR COMPENSATIONS AND ENVIRONMENT v I Σ v R Control o industrial robots Control with eternal sensors P. Roo [32]

33 Soures o inonsisteny Hybrid position/ore ontrol is based on a nominal model o the interation. Inonsisteny may however our in the measurements, due e.g. to: 1. rition at the ontat (a ore is deteted in a nominally ree diretion) 2. ompliane in the robot struture and/or at the ontat (a displaement is deteted in a diretion whih is nominally onstrained in motion) 3. unertainty in the environment geometry at the ontat The irst two soures o inonsisteny are automatially iltered out through the seletion matries. The third soure o inonsisteny an be mitigated by real time estimation proess. Control o industrial robots Control with eternal sensors P. Roo [33]

34 A 1 d.o.. system In order to have a better idea on how the ontrol system operates along the ore ontrolled diretions, let s onsider the simple eample o a body o mass M in interation with the environment: u M We assume a purely elasti model o the environment. Blok diagram: u + G(s) K u G (s) G K ( s) = 2 Ms + K Control o industrial robots Control with eternal sensors P. Roo [34]

35 A 1 d.o.. system u G (s) G K ( s) = 2 Ms + K The system presents a unitary gain and resonant, high requeny (or a sti ontat) dynamis. Notie that = u or K. The best way to ontrol this system is through an integral ontroller: + u d R (s) G (s) R ( s ) = K s i This way the bandwidth o the ore ontrol is diretly assigned with the integrator, leaving outside the bandwidth the high requeny dynamis o the interation with the environment. Control o industrial robots Control with eternal sensors P. Roo [35]

36 Impliit ore ontrol As with impedane ontrol, when a position ontrol is already ative, it might be onvenient to lose the ore ontrol around the position ontrol loop. Consider or simpliity only the ase o a, position ontrolled, single degree o reedom: d R u M The system an be modelled through the ollowing blok diagram: d + R(s) u + G(s) K Control o industrial robots Control with eternal sensors P. Roo [36]

37 Impliit ore ontrol Under the assumption o an ininitely sti ontat (K ) we obtain =0 and u: d u R(s) 1 I a PID position ontroller is used: R ( s) we an adopt or the ore ontrol the sheme: d + R (s) C(s) d The position ontrol loop opens! = 2 KD s + KPs + K s u R(s) 1 I with: ( ) =, R ( s) C s K D s K P s + K I = K p + K s i The advantage is that the position ontrol is let unhanged. Control o industrial robots Control with eternal sensors P. Roo [37]

38 Industrial appliation o ore ontrol ABB RobotWare Mahining Fore Control (mahining o blades): Control o industrial robots Control with eternal sensors P. Roo [38]

39 Visual measurements Artiiial vision devies are useul sensors or robotis beause they mimi the human sense o sight and allow to gather inormation rom the environment without ontat. Nowadays several roboti ontrollers integrate vision systems. The typial use o vision in industrial robotis is to detet an objet in the robot s sene, whose position (and orientation) is then used or online path planning in order to drive the robot to the identiied objet. Online re-planning o the path an also be perormed when the vision system detets some unepeted hange in the path the robot is supposed to ollow (or eample a orner in a ontouring task). Alternatively, visual measurements an be used in a real time eedbak loop in order to improve position ontrol o the end eetor: this is the onept o visual servoing. Control o industrial robots Control with eternal sensors P. Roo [39]

40 An eample o path re-planning Following a orner during a ontouring task: View rom an eternal amera View rom the onboard amera Control o industrial robots Control with eternal sensors P. Roo [40]

41 The hardware amera video data rame grabber industrial PC zoom lens Control o industrial robots Control with eternal sensors P. Roo [41]

42 The amera The amera is a devie that an measure the intensity o light, onentrated by a lens on a plane, the image plane, whih ontains a matri o piels (CCD: Charge Coupled Devie). Piels then transorm the light energy in eletrial energy. The amera perorms a 2D projetion o the sene. This projetion entails a loss o depth inormation: eah point in the image plane orresponds to a ray in the 3D spae. In order to determine the 3D oordinates o a point orresponding to a 2D point in the plane additional inormation is needed: multiple views with a single amera multiple ameras knowledge o harateristi relations between relevant points o the ramed objets Control o industrial robots Control with eternal sensors P. Roo [42]

43 Perspetive projetion A point with oordinates (X, Y, Z) in the amera rame is projeted into the ollowing point in the image plane, epressed in piels: ξ = y = α λ Z X Y λ: oal length α: sale ator in piel/m Other methods eist to represent the projetion on the 2D plane (saled orthographi projetion, aine projetion) Control o industrial robots Control with eternal sensors P. Roo [43]

44 Image eatures In artiiial vision we denote with image eature any harateristis that an be etrated rom an image (e.g. an edge or a orner). We then deine a parameter o an image eature a quantity, epressed by a real numeri value, whih an be omputed rom one or more image eatures. Parameters o an image eature an be gathered in a vetor: ξ = [ξ 1, ξ 2, ξ k ]. Eamples o parameters o image eatures: point oordinates length and orientation o a line onneting two points entroids and higher order moments parameters o an ellipse Control o industrial robots Control with eternal sensors P. Roo [44]

45 Calibration The amera has to be alibrated beore usage in a roboti vision system: Internal alibration: Determination o the intrinsi parameters o the amera (like the oal length λ) as well as o some additional distortion parameters due to lens imperetions and misalignments in the optial system Eternal alibration: Determination o the etrinsi parameters o the amera like the position and the orientation o the amera with respet to a reerene rame Control o industrial robots Control with eternal sensors P. Roo [45]

46 Camera oniguration The irst deision to be made when setting up a vision-based ontrol system is where to plae the amera. The amera an be: mounted in a ied loation in the workspae (eye-to-hand oniguration) so that it an observe the manipulator and any objets to be manipulated attahed to the robot above the wrist (eye-in-hand oniguration) Eye-To-Hand Eye-In-Hand Control o industrial robots Control with eternal sensors P. Roo [46]

47 Eye-to-hand oniguration Advantages the ield o view does not hange as the manipulator moves the geometri relationship between the amera and the workspae is ied and an be alibrated oline Eye-To-Hand Disadvantages as the manipulator moves through the workspae it an olude the amera s ield o view Control o industrial robots Control with eternal sensors P. Roo [47]

48 Eye-in-hand oniguration Advantages the amera an observe the motion o the end eetor at a ied resolution and without olusion as the manipulator moves through the workspae Eye-In-Hand Disadvantages the geometri relationship between the amera and the workspae hanges as the manipulator moves the ield o view an hange drmatially or even small motion o the manipulator Control o industrial robots Control with eternal sensors P. Roo [48]

49 Control arhitetures: lassiiation Vision ontrol systems an be lassiied based on various riteria. A irst lassiiation is based on the ollowing question: Is the ontrol struture hierarhial, with the vision system providing set-points as input to robot s joint-level ontroller, or does the visual ontroller diretly ompute the joint-level inputs? In the irst ase: In the seond ase: dynami look and move diret visual servoing Advantages o the dynami look and move approah: the redued sampling rate o the visual signal does not ompromise the overall perormane o the position ontrol system in several industrial robot ontrollers it is only allowed to operate at the position setpoints level the robot an be seen as an ideal positioner in the Cartesian spae, thus simpliying the design o the vision ontrol system Control o industrial robots Control with eternal sensors P. Roo [49]

50 Control arhitetures: lassiiation A seond lassiiation is based on the ollowing question: Is the error signal deined in 3D (task spae) oordinates or diretly in terms o image eatures? In the irst ase : In the seond ase : position based ontrol image based ontrol Position based ontrol: vision data are used to build a partial 3D representation o the world pose estimation algorithms are omputationally intensive (a real-time implementation is required) and sensitive to errors in amera alibration Image based ontrol: uses the image data diretly to ontrol the robot motion an error untion is deined in terms o quantities that an be diretly measured in an image, and a ontrol law is onstruted that maps this error diretly to robot motion Control o industrial robots Control with eternal sensors P. Roo [50]

51 Typial ontrol arhitetures Dynami position-based look-and-move d + Cartesian ontrol law Joint ontrollers Atuators Camera Pose estimation Image eature etration Video Dynami image-based look-and-move + Control law in the image eature spae Joint ontrollers Atuators Camera Image eature etration Video Control o industrial robots Control with eternal sensors P. Roo [51]

52 Typial ontrol arhitetures Position-based visual servoing d + Cartesian ontrol law Atuators Camera Pose estimation Image eature etration Video Image-based visual servoing + Control law in the image eature spae Atuators Camera Image eature etration Video Control o industrial robots Control with eternal sensors P. Roo [52]

Control of industrial robots. Control of the interaction

Control of industrial robots. Control of the interaction Control o industrial robots Control o the interation Pro. Paolo Roo (paolo.roo@polimi.it) Politenio di Milano Dipartimento di Elettronia, Inormazione e Bioingegneria Introdution So ar we have assumed that

More information

v = fy c u = fx c z c The Pinhole Camera Model Camera Projection Models

v = fy c u = fx c z c The Pinhole Camera Model Camera Projection Models The Pinhole Camera Model Camera Projetion Models We will introdue dierent amera projetion models that relate the loation o an image point to the oordinates o the orresponding 3D points. The projetion models

More information

f 2 f n where m is the total mass of the object. Expression (6a) is plotted in Figure 8 for several values of damping ( ).

f 2 f n where m is the total mass of the object. Expression (6a) is plotted in Figure 8 for several values of damping ( ). F o F o / k A = = 6 k 1 + 1 + n r n n n RESONANCE It is seen in Figure 7 that displaement and stress levels tend to build up greatly when the oring requeny oinides with the natural requeny, the buildup

More information

Finite Formulation of Electromagnetic Field

Finite Formulation of Electromagnetic Field Finite Formulation o Eletromagneti Field Enzo TONTI Dept.Civil Engin., Univ. o Trieste, Piazzale Europa 1, 34127 Trieste, Italia. e-mail: tonti@univ.trieste.it Otober 16, 2000 Abstrat This paper shows

More information

Control Theory association of mathematics and engineering

Control Theory association of mathematics and engineering Control Theory assoiation of mathematis and engineering Wojieh Mitkowski Krzysztof Oprzedkiewiz Department of Automatis AGH Univ. of Siene & Tehnology, Craow, Poland, Abstrat In this paper a methodology

More information

The gravitational phenomena without the curved spacetime

The gravitational phenomena without the curved spacetime The gravitational phenomena without the urved spaetime Mirosław J. Kubiak Abstrat: In this paper was presented a desription of the gravitational phenomena in the new medium, different than the urved spaetime,

More information

Force/Vision Based Active Damping Control of Contact Transition in Dynamic Environments

Force/Vision Based Active Damping Control of Contact Transition in Dynamic Environments Fore/Vision Based Ative Damping Control of Contat Transition in Dynami Environments Tomas Olsson, Rolf Johansson, Anders Robertsson Department of Automati Control, Lund Inst. of Teh., Lund University,

More information

Practice Exam 2 Solutions

Practice Exam 2 Solutions MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department o Physis Physis 801T Fall Term 004 Problem 1: stati equilibrium Pratie Exam Solutions You are able to hold out your arm in an outstrethed horizontal position

More information

Development of moving sound source localization system

Development of moving sound source localization system Eletroni Journal «Tehnial Aoustis» http://www.ejta.org 2006, 11 Doh-Hyoung Kim 1, Youngjin Park 2 Korea Advaned Institute o Siene and Tehnology, ME3076, 373-1, Guseong-Dong, Yuseong-Gu, Daejon, 305-701,

More information

Doppler effect of the rupture process of the great M W 7.9 Wenchuan earthquake

Doppler effect of the rupture process of the great M W 7.9 Wenchuan earthquake Earthq Si (1)3: 535 539 535 539 Doi: 1.17/s11589-1-75-4 Doppler eet o the rupture proess o the great M W 7.9 Wenhuan earthquake Ge Jin 1,, Youai Tang 1 Shiyong Zhou 1 and Yongshun John Chen 1 1 Institute

More information

Advanced Computational Fluid Dynamics AA215A Lecture 4

Advanced Computational Fluid Dynamics AA215A Lecture 4 Advaned Computational Fluid Dynamis AA5A Leture 4 Antony Jameson Winter Quarter,, Stanford, CA Abstrat Leture 4 overs analysis of the equations of gas dynamis Contents Analysis of the equations of gas

More information

Aharonov-Bohm effect. Dan Solomon.

Aharonov-Bohm effect. Dan Solomon. Aharonov-Bohm effet. Dan Solomon. In the figure the magneti field is onfined to a solenoid of radius r 0 and is direted in the z- diretion, out of the paper. The solenoid is surrounded by a barrier that

More information

Hankel Optimal Model Order Reduction 1

Hankel Optimal Model Order Reduction 1 Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both

More information

Acoustic Attenuation Performance of Helicoidal Resonator Due to Distance Change from Different Cross-sectional Elements of Cylindrical Ducts

Acoustic Attenuation Performance of Helicoidal Resonator Due to Distance Change from Different Cross-sectional Elements of Cylindrical Ducts Exerpt rom the Proeedings o the COMSOL Conerene 1 Paris Aousti Attenuation Perormane o Helioidal Resonator Due to Distane Change rom Dierent Cross-setional Elements o Cylindrial Duts Wojieh ŁAPKA* Division

More information

DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS

DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS CHAPTER 4 DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS 4.1 INTRODUCTION Around the world, environmental and ost onsiousness are foring utilities to install

More information

A novel Infrared Thermography (IRT) based experimental technique for distributed temperature measurements in hot gas flows

A novel Infrared Thermography (IRT) based experimental technique for distributed temperature measurements in hot gas flows QIRT 1 1 http://dx.doi.org/1.1611/qirt.1.45 th International Conerene on Quantitative InraRed Thermography July 7-3, 1, Québe (Canada A novel Inrared Thermography (IRT based experimental tehnique or distributed

More information

Relativity in Classical Physics

Relativity in Classical Physics Relativity in Classial Physis Main Points Introdution Galilean (Newtonian) Relativity Relativity & Eletromagnetism Mihelson-Morley Experiment Introdution The theory of relativity deals with the study of

More information

Contact State Estimation using Multiple Model Estimation and Hidden Markov Models

Contact State Estimation using Multiple Model Estimation and Hidden Markov Models Thomas J. Debus Pierre E. Dupont Department of Aerospae and Mehanial Engineering Boston University Boston, MA 02445, USA Robert D. Howe Division of Engineering and Applied Sienes Harvard University Cambridge,

More information

PREDICTION OF CONCRETE COMPRESSIVE STRENGTH

PREDICTION OF CONCRETE COMPRESSIVE STRENGTH PREDICTION OF CONCRETE COMPRESSIVE STRENGTH Dunja Mikuli (1), Ivan Gabrijel (1) and Bojan Milovanovi (1) (1) Faulty o Civil Engineering, University o Zagreb, Croatia Abstrat A ompressive strength o onrete

More information

Motor Sizing Application Note

Motor Sizing Application Note PAE-TILOGY Linear Motors 70 Mill orest d. Webster, TX 77598 (8) 6-7750 ax (8) 6-7760 www.trilogysystems.om E-mail emn_support_trilogy@parker.om Motor Sizing Appliation Note By Jak Marsh Introdution Linear

More information

Measuring & Inducing Neural Activity Using Extracellular Fields I: Inverse systems approach

Measuring & Inducing Neural Activity Using Extracellular Fields I: Inverse systems approach Measuring & Induing Neural Ativity Using Extraellular Fields I: Inverse systems approah Keith Dillon Department of Eletrial and Computer Engineering University of California San Diego 9500 Gilman Dr. La

More information

An Adaptive Optimization Approach to Active Cancellation of Repeated Transient Vibration Disturbances

An Adaptive Optimization Approach to Active Cancellation of Repeated Transient Vibration Disturbances An aptive Optimization Approah to Ative Canellation of Repeated Transient Vibration Disturbanes David L. Bowen RH Lyon Corp / Aenteh, 33 Moulton St., Cambridge, MA 138, U.S.A., owen@lyonorp.om J. Gregory

More information

Relativistic Dynamics

Relativistic Dynamics Chapter 7 Relativisti Dynamis 7.1 General Priniples of Dynamis 7.2 Relativisti Ation As stated in Setion A.2, all of dynamis is derived from the priniple of least ation. Thus it is our hore to find a suitable

More information

Dynamic Programming and Multi Objective Linear Programming approaches

Dynamic Programming and Multi Objective Linear Programming approaches 0 NSP Vol. () (0), 9 Dynami Programming and Multi Objetive Linear Programming approahes P. K. De and Amita Bhinher Department o Mathematis, National Institute o Tehnology SILCHAR - 788 00 (Assam) India

More information

Normative and descriptive approaches to multiattribute decision making

Normative and descriptive approaches to multiattribute decision making De. 009, Volume 8, No. (Serial No.78) China-USA Business Review, ISSN 57-54, USA Normative and desriptive approahes to multiattribute deision making Milan Terek (Department of Statistis, University of

More information

Robust Flight Control Design for a Turn Coordination System with Parameter Uncertainties

Robust Flight Control Design for a Turn Coordination System with Parameter Uncertainties Amerian Journal of Applied Sienes 4 (7): 496-501, 007 ISSN 1546-939 007 Siene Publiations Robust Flight ontrol Design for a urn oordination System with Parameter Unertainties 1 Ari Legowo and Hiroshi Okubo

More information

University of Cambridge Engineering Part IIB Module 4F12: Computer Vision Handout 3: Projection

University of Cambridge Engineering Part IIB Module 4F12: Computer Vision Handout 3: Projection University of Cambridge Engineering Part IIB Module 4F2: Computer Vision Handout 3: Projetion Roberto Cipolla Otober 2009 Projetion Orthographi projetion Reall that omputer vision is about disovering from

More information

Orthogonal Complement Based Divide-and-Conquer Algorithm (O-DCA) for Constrained Multibody Systems

Orthogonal Complement Based Divide-and-Conquer Algorithm (O-DCA) for Constrained Multibody Systems Orthogonal Complement Based Divide-and-Conquer Algorithm (O-DCA) for Constrained Multibody Systems Rudranarayan M. Mukherjee, Kurt S. Anderson Computational Dynamis Laboratory Department of Mehanial Aerospae

More information

ETNA Kent State University

ETNA Kent State University % { Eletroni Transations on Numerial Analysis. Volume 21, pp. 20-27, 2005. Copyright 2005,. ISSN 1068-9613. QR FACTORIZATIONS USING A RESTRICTED SET OF ROTATIONS DIANNE P. O LEARY AND STEPHEN S. BULLOCK

More information

Part G-4: Sample Exams

Part G-4: Sample Exams Part G-4: Sample Exams 1 Cairo University M.S.: Eletronis Cooling Faulty of Engineering Final Exam (Sample 1) Mehanial Power Engineering Dept. Time allowed 2 Hours Solve as muh as you an. 1. A heat sink

More information

GEOMETRY FOR 3D COMPUTER VISION

GEOMETRY FOR 3D COMPUTER VISION GEOMETRY FOR 3D COMPUTER VISION 1/29 Válav Hlaváč Czeh Tehnial University, Faulty of Eletrial Engineering Department of Cybernetis, Center for Mahine Pereption 121 35 Praha 2, Karlovo nám. 13, Czeh Republi

More information

Bayesian Optimization Under Uncertainty

Bayesian Optimization Under Uncertainty Bayesian Optimization Under Unertainty Justin J. Beland University o Toronto justin.beland@mail.utoronto.a Prasanth B. Nair University o Toronto pbn@utias.utoronto.a Abstrat We onsider the problem o robust

More information

The Concept of Mass as Interfering Photons, and the Originating Mechanism of Gravitation D.T. Froedge

The Concept of Mass as Interfering Photons, and the Originating Mechanism of Gravitation D.T. Froedge The Conept of Mass as Interfering Photons, and the Originating Mehanism of Gravitation D.T. Froedge V04 Formerly Auburn University Phys-dtfroedge@glasgow-ky.om Abstrat For most purposes in physis the onept

More information

Department of Mechanical Engineering

Department of Mechanical Engineering Department o Mehanial Engineering AMEE41 / ATO4 Aerodynamis Instrutor: Marios M. Fyrillas Email: eng.m@it.a.y Homework Assignment #4 QESTION 1 Consider the boundary layer low on a lat plate o width b (shown

More information

Kinematics and inverse dynamics analysis for a general 3-PRS spatial parallel mechanism Yangmin Li and Qingsong Xu

Kinematics and inverse dynamics analysis for a general 3-PRS spatial parallel mechanism Yangmin Li and Qingsong Xu Robotia (2005) volume 23, pp. 219 229. 2005 Cambridge University Press doi:10.1017/s0263574704000797 Printed in the United Kingdom Kinematis and inverse dynamis analysis for a general 3-PRS spatial parallel

More information

SEISMIC ANALYSIS OF SPHERICAL TANKS INCLUDING FLUID-STRUCTURE-SOIL INTERACTION

SEISMIC ANALYSIS OF SPHERICAL TANKS INCLUDING FLUID-STRUCTURE-SOIL INTERACTION 3 th World Conferene on Earthquake Engineering Vanouver, B.C., Canada August -6, 2004 aper o. 84 SEISMIC AALYSIS OF SHERICAL TAKS ICLUDIG FLUID-STRUCTURE-SOIL ITERACTIO T.L. Karavasilis, D.C. Rizos 2,

More information

Effect of Droplet Distortion on the Drag Coefficient in Accelerated Flows

Effect of Droplet Distortion on the Drag Coefficient in Accelerated Flows ILASS Amerias, 19 th Annual Conerene on Liquid Atomization and Spray Systems, Toronto, Canada, May 2006 Eet o Droplet Distortion on the Drag Coeiient in Aelerated Flows Shaoping Quan, S. Gopalakrishnan

More information

PIEZOELECTRIC materials are widely used for structural

PIEZOELECTRIC materials are widely used for structural Parametri Analysis o the Vibration Control o Sandwih Beams Through Shear-Based Piezoeletri Atuation M. A. TRINDADE,* A. BENJEDDOU AND R. OHAYON Strutural Mehanis and Coupled Systems Laboratory, Conservatoire

More information

A NEW FLEXIBLE BODY DYNAMIC FORMULATION FOR BEAM STRUCTURES UNDERGOING LARGE OVERALL MOTION IIE THREE-DIMENSIONAL CASE. W. J.

A NEW FLEXIBLE BODY DYNAMIC FORMULATION FOR BEAM STRUCTURES UNDERGOING LARGE OVERALL MOTION IIE THREE-DIMENSIONAL CASE. W. J. A NEW FLEXIBLE BODY DYNAMIC FORMULATION FOR BEAM STRUCTURES UNDERGOING LARGE OVERALL MOTION IIE THREE-DIMENSIONAL CASE W. J. Haering* Senior Projet Engineer General Motors Corporation Warren, Mihigan R.

More information

13 Fitzroy Street London W1T 4BQ Telephone: +44 (0) Facsimile: +44 (0)

13 Fitzroy Street London W1T 4BQ Telephone: +44 (0) Facsimile: +44 (0) Oasys AdSe Theory 13 Fitzroy Street London W1T 4BQ Telephone: +44 (0) 0 7755 330 Fasimile: +44 (0) 0 7755 370 Central Square Forth Street Newastle Upon Tyne NE1 3PL Telephone: +44 (0) 191 38 7559 Fasimile:

More information

Where as discussed previously we interpret solutions to this partial differential equation in the weak sense: b

Where as discussed previously we interpret solutions to this partial differential equation in the weak sense: b Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential

More information

Physical Laws, Absolutes, Relative Absolutes and Relativistic Time Phenomena

Physical Laws, Absolutes, Relative Absolutes and Relativistic Time Phenomena Page 1 of 10 Physial Laws, Absolutes, Relative Absolutes and Relativisti Time Phenomena Antonio Ruggeri modexp@iafria.om Sine in the field of knowledge we deal with absolutes, there are absolute laws that

More information

SURFACE WAVES OF NON-RAYLEIGH TYPE

SURFACE WAVES OF NON-RAYLEIGH TYPE SURFACE WAVES OF NON-RAYLEIGH TYPE by SERGEY V. KUZNETSOV Institute for Problems in Mehanis Prosp. Vernadskogo, 0, Mosow, 75 Russia e-mail: sv@kuznetsov.msk.ru Abstrat. Existene of surfae waves of non-rayleigh

More information

9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry

9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry 9 Geophysis and Radio-Astronomy: VLBI VeryLongBaseInterferometry VLBI is an interferometry tehnique used in radio astronomy, in whih two or more signals, oming from the same astronomial objet, are reeived

More information

The simulation analysis of the bridge rectifier continuous operation in AC circuit

The simulation analysis of the bridge rectifier continuous operation in AC circuit Computer Appliations in Eletrial Engineering Vol. 4 6 DOI 8/j.8-448.6. The simulation analysis of the bridge retifier ontinuous operation in AC iruit Mirosław Wiślik, Paweł Strząbała Kiele University of

More information

The Unified Geometrical Theory of Fields and Particles

The Unified Geometrical Theory of Fields and Particles Applied Mathematis, 014, 5, 347-351 Published Online February 014 (http://www.sirp.org/journal/am) http://dx.doi.org/10.436/am.014.53036 The Unified Geometrial Theory of Fields and Partiles Amagh Nduka

More information

Comparison of Alternative Equivalent Circuits of Induction Motor with Real Machine Data

Comparison of Alternative Equivalent Circuits of Induction Motor with Real Machine Data Comparison of Alternative Equivalent Ciruits of Indution Motor with Real Mahine Data J. radna, J. auer, S. Fligl and V. Hlinovsky Abstrat The algorithms based on separated ontrol of the motor flux and

More information

Grasp Planning: How to Choose a Suitable Task Wrench Space

Grasp Planning: How to Choose a Suitable Task Wrench Space Grasp Planning: How to Choose a Suitable Task Wrenh Spae Ch. Borst, M. Fisher and G. Hirzinger German Aerospae Center - DLR Institute for Robotis and Mehatronis 8223 Wessling, Germany Email: [Christoph.Borst,

More information

Conformal Mapping among Orthogonal, Symmetric, and Skew-Symmetric Matrices

Conformal Mapping among Orthogonal, Symmetric, and Skew-Symmetric Matrices AAS 03-190 Conformal Mapping among Orthogonal, Symmetri, and Skew-Symmetri Matries Daniele Mortari Department of Aerospae Engineering, Texas A&M University, College Station, TX 77843-3141 Abstrat This

More information

arxiv:nucl-th/ v1 27 Jul 1999

arxiv:nucl-th/ v1 27 Jul 1999 Eetive Widths and Eetive Number o Phonons o Multiphonon Giant Resonanes L.F. Canto, B.V. Carlson, M.S. Hussein 3 and A.F.R. de Toledo Piza 3 Instituto de Físia, Universidade do Rio de Janeiro, CP 6858,

More information

TENSOR FORM OF SPECIAL RELATIVITY

TENSOR FORM OF SPECIAL RELATIVITY TENSOR FORM OF SPECIAL RELATIVITY We begin by realling that the fundamental priniple of Speial Relativity is that all physial laws must look the same to all inertial observers. This is easiest done by

More information

arxiv:physics/ v4 [physics.gen-ph] 9 Oct 2006

arxiv:physics/ v4 [physics.gen-ph] 9 Oct 2006 The simplest derivation of the Lorentz transformation J.-M. Lévy Laboratoire de Physique Nuléaire et de Hautes Energies, CNRS - IN2P3 - Universités Paris VI et Paris VII, Paris. Email: jmlevy@in2p3.fr

More information

ELECTROMAGNETIC WAVES WITH NONLINEAR DISPERSION LAW. P. М. Меdnis

ELECTROMAGNETIC WAVES WITH NONLINEAR DISPERSION LAW. P. М. Меdnis ELECTROMAGNETIC WAVES WITH NONLINEAR DISPERSION LAW P. М. Меdnis Novosibirs State Pedagogial University, Chair of the General and Theoretial Physis, Russia, 636, Novosibirs,Viljujsy, 8 e-mail: pmednis@inbox.ru

More information

A model for measurement of the states in a coupled-dot qubit

A model for measurement of the states in a coupled-dot qubit A model for measurement of the states in a oupled-dot qubit H B Sun and H M Wiseman Centre for Quantum Computer Tehnology Centre for Quantum Dynamis Griffith University Brisbane 4 QLD Australia E-mail:

More information

Directional Coupler. 4-port Network

Directional Coupler. 4-port Network Diretional Coupler 4-port Network 3 4 A diretional oupler is a 4-port network exhibiting: All ports mathed on the referene load (i.e. S =S =S 33 =S 44 =0) Two pair of ports unoupled (i.e. the orresponding

More information

IMPACT MODELLING OF THE COEFFICIENT OF RESTITUTION OF POTATOES BASED ON THE KELVIN- VOIGHT PAIR

IMPACT MODELLING OF THE COEFFICIENT OF RESTITUTION OF POTATOES BASED ON THE KELVIN- VOIGHT PAIR Bulletin of the Transilvania University of Braşov Series II: Forestry Wood Industry Agriultural Food Engineering Vol. 9 (58) No. - 06 IMPACT MODELLING OF THE COEFFICIENT OF RESTITUTION OF POTATOES BASED

More information

An Integrated Architecture of Adaptive Neural Network Control for Dynamic Systems

An Integrated Architecture of Adaptive Neural Network Control for Dynamic Systems An Integrated Arhiteture of Adaptive Neural Network Control for Dynami Systems Robert L. Tokar 2 Brian D.MVey2 'Center for Nonlinear Studies, 2Applied Theoretial Physis Division Los Alamos National Laboratory,

More information

Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite.

Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite. Leture Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the funtion V ( x ) to be positive definite. ost often, our interest will be to show that x( t) as t. For that we will need

More information

Speed-feedback Direct-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion

Speed-feedback Direct-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion Speed-feedbak Diret-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion Y. Yamamoto, T. Nakamura 2, Y. Takada, T. Koseki, Y. Aoyama 3, and Y. Iwaji 3

More information

3 Tidal systems modelling: ASMITA model

3 Tidal systems modelling: ASMITA model 3 Tidal systems modelling: ASMITA model 3.1 Introdution For many pratial appliations, simulation and predition of oastal behaviour (morphologial development of shorefae, beahes and dunes) at a ertain level

More information

Determination of the reaction order

Determination of the reaction order 5/7/07 A quote of the wee (or amel of the wee): Apply yourself. Get all the eduation you an, but then... do something. Don't just stand there, mae it happen. Lee Iaoa Physial Chemistry GTM/5 reation order

More information

Wave Propagation through Random Media

Wave Propagation through Random Media Chapter 3. Wave Propagation through Random Media 3. Charateristis of Wave Behavior Sound propagation through random media is the entral part of this investigation. This hapter presents a frame of referene

More information

Four-dimensional equation of motion for viscous compressible substance with regard to the acceleration field, pressure field and dissipation field

Four-dimensional equation of motion for viscous compressible substance with regard to the acceleration field, pressure field and dissipation field Four-dimensional equation of motion for visous ompressible substane with regard to the aeleration field, pressure field and dissipation field Sergey G. Fedosin PO box 6488, Sviazeva str. -79, Perm, Russia

More information

Speed Regulation of a Small BLDC Motor using Genetic-Based Proportional Control

Speed Regulation of a Small BLDC Motor using Genetic-Based Proportional Control World Aademy of Siene, Engineering and Tehnology 47 8 Speed Regulation of a Small BLDC Motor using Geneti-Based Proportional Control S. Poonsawat, and T. Kulworawanihpong Abstrat This paper presents the

More information

INFORMATION CONCERNING MATERIALS TO BE USED IN THE DESIGN

INFORMATION CONCERNING MATERIALS TO BE USED IN THE DESIGN TITLE 5 DESIGN CHAPTER 8 INFORMATION CONCERNING MATERIALS TO BE USED IN THE DESIGN Artile 38. Charateristis o steel or reinorements 38.1 General The harateristis o the steel used or the design desribed

More information

Experimental Investigation and FE Analysis of Fiber Woven Layered Composites under Dynamic Loading

Experimental Investigation and FE Analysis of Fiber Woven Layered Composites under Dynamic Loading 2th International LS-DYNA Users Conferene Constitutive Modeling(2) xperimental Investigation and F Analysis of Fiber Woven Layered Composites under Dynami Loading Pavel A. Mossakovsky, Fedor K. Antonov,

More information

Modes are solutions, of Maxwell s equation applied to a specific device.

Modes are solutions, of Maxwell s equation applied to a specific device. Mirowave Integrated Ciruits Prof. Jayanta Mukherjee Department of Eletrial Engineering Indian Institute of Tehnology, Bombay Mod 01, Le 06 Mirowave omponents Welome to another module of this NPTEL mok

More information

Influence of transverse cracks on the onset of delamination: application to L-angle specimens. F. Laurin*, A. Mavel, P. Nuñez, E.

Influence of transverse cracks on the onset of delamination: application to L-angle specimens. F. Laurin*, A. Mavel, P. Nuñez, E. Influene of transverse raks on the onset of delamination: appliation to L-angle speimens F. Laurin*, A. Mavel, P. Nuñez, E. Auguste Composite strutures subjeted to 3D loading Wings Strutures under 3D loadings

More information

Finite Element Modelling of Vibro-Acoustic Systems for Active Noise Reduction

Finite Element Modelling of Vibro-Acoustic Systems for Active Noise Reduction ECHNISCHE MECHANIK, Band 25, Heft 3-4, (25), 241-247 Manuskripteingang: 11. Dezember 25 Finite Element Modelling of Vibro-Aousti Systems for Ative Noise Redution J. Lefèvre, U. Gabbert Over the past years

More information

A Big Bang-Big Crunch Optimization Based Approach for Interval Type-2 Fuzzy PID Controller Design

A Big Bang-Big Crunch Optimization Based Approach for Interval Type-2 Fuzzy PID Controller Design A Big Bang-Big Crunh Optimization Based Approah or Interval Type-2 Fuzzy PID Controller Design Tuan Kumbasar Control Engineering Department Istanbul Tehnial University Istanbul, Turkey kumbasart@itu.edu.tr

More information

DO PHYSICS ONLINE. SPECIAL RELATIVITY Frames of Reference

DO PHYSICS ONLINE. SPECIAL RELATIVITY Frames of Reference DO PHYSICS ONLINE SPACE SPECIAL RELATIVITY Frames of Referene Spae travel Apollo 11 spaeraft: Earth Moon v ~ 40x10 3 km.h -1 Voyager spaeraft: v ~ 60x10 3 km.h -1 (no sling shot effet) Ulysses spaeraft:

More information

Properties of Quarks

Properties of Quarks PHY04 Partile Physis 9 Dr C N Booth Properties of Quarks In the earlier part of this ourse, we have disussed three families of leptons but prinipally onentrated on one doublet of quarks, the u and d. We

More information

The homopolar generator: an analytical example

The homopolar generator: an analytical example The homopolar generator: an analytial example Hendrik van Hees August 7, 214 1 Introdution It is surprising that the homopolar generator, invented in one of Faraday s ingenious experiments in 1831, still

More information

PHY 396 T: SUSY Solutions for problem set #12.

PHY 396 T: SUSY Solutions for problem set #12. PHY 396 T: SUSY Solutions or problem set #. Problem a: In priniple the non-perturbative superpotential o the theory may depend on the dual quark and antiquark ields q and q as well as the singlets Φ but

More information

Beams on Elastic Foundation

Beams on Elastic Foundation Professor Terje Haukaas University of British Columbia, Vanouver www.inrisk.ub.a Beams on Elasti Foundation Beams on elasti foundation, suh as that in Figure 1, appear in building foundations, floating

More information

PhysicsAndMathsTutor.com 1

PhysicsAndMathsTutor.com 1 PhysisAndMathsTutor.om. (a (i beam splitter [or semi-silvered mirror] (ii a ompensator [or a glass blok] allows for the thikness of the (semi-silvered mirror to obtain equal optial path lengths in the

More information

Stabilization of the Precision Positioning Stage Working in the Vacuum Environment by Using the Disturbance Observer

Stabilization of the Precision Positioning Stage Working in the Vacuum Environment by Using the Disturbance Observer Proeedings of the 4th IIAE International Conferene on Industrial Appliation Engineering 216 Stabilization of the Preision Positioning Stage Working in the Vauum Environment by Using the Disturbane Observer

More information

Particle-wave symmetry in Quantum Mechanics And Special Relativity Theory

Particle-wave symmetry in Quantum Mechanics And Special Relativity Theory Partile-wave symmetry in Quantum Mehanis And Speial Relativity Theory Author one: XiaoLin Li,Chongqing,China,hidebrain@hotmail.om Corresponding author: XiaoLin Li, Chongqing,China,hidebrain@hotmail.om

More information

STRUCTURAL AND BEHAVIORAL OPTIMIZATION OF THE NONLINEAR HILL MODEL

STRUCTURAL AND BEHAVIORAL OPTIMIZATION OF THE NONLINEAR HILL MODEL THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Series A, OF THE ROMANIAN ACADEMY Volume 12, Number 3/2011, pp. 213 220 STRUCTURAL AND BEHAVIORAL OPTIMIZATION OF THE NONLINEAR HILL MODEL Tudor

More information

Inference-based Ambiguity Management in Decentralized Decision-Making: Decentralized Diagnosis of Discrete Event Systems

Inference-based Ambiguity Management in Decentralized Decision-Making: Decentralized Diagnosis of Discrete Event Systems Inerene-based Ambiguity Management in Deentralized Deision-Making: Deentralized Diagnosis o Disrete Event Systems Ratnesh Kumar and Shigemasa Takai Abstrat The task o deentralized deision-making involves

More information

Vibration Control of Smart Structure Using Sliding Mode Control with Observer

Vibration Control of Smart Structure Using Sliding Mode Control with Observer JOURNAL OF COMPUERS, VOL. 7, NO., FEBRUARY Vibration Control of Smart Struture Using Sliding Mode Control with Observer Junfeng Hu Shool of Mehanial & Eletrial Engineering, Jiangxi University of Siene

More information

A note on a variational formulation of electrodynamics

A note on a variational formulation of electrodynamics Proeedings of the XV International Workshop on Geometry and Physis Puerto de la Cruz, Tenerife, Canary Islands, Spain September 11 16, 006 Publ. de la RSME, Vol. 11 (007), 314 31 A note on a variational

More information

A multiscale description of failure in granular materials

A multiscale description of failure in granular materials A multisale desription of failure in granular materials Nejib Hadda, François Niot, Lu Sibille, Farhang Radjai, Antoinette Tordesillas et al. Citation: AIP Conf. Pro. 154, 585 (013); doi: 10.1063/1.4811999

More information

An Evolutionary Algorithm for Constrained Optimization

An Evolutionary Algorithm for Constrained Optimization An Evolutionary Algorim or Constrained Optimization Tapabrata Ray, Tai Kang and Seow Kian Chye Center or Advaned Numerial Engineering Simulations Shool o ehanial and Prodution Engineering Nanyang Tehnologial

More information

Active Magnetic Bearings for Frictionless Rotating Machineries

Active Magnetic Bearings for Frictionless Rotating Machineries Ative Magneti Bearings for Fritionless Rotating Mahineries Joga Dharma Setiawan Abstrat Ative magneti bearing (AMB systems an support a rotor without physial ontat and enable users to preisely ontrol rotor

More information

Experiment 3: Basic Electronic Circuits II (tbc 1/7/2007)

Experiment 3: Basic Electronic Circuits II (tbc 1/7/2007) Experiment 3: Basi Eletroni iruits II (tb /7/007) Objetive: a) To study the first-order dynamis of a apaitive iruits with the appliation of Kirhoff s law, Ohm s law and apaitane formula. b) To learn how

More information

Vector Field Theory (E&M)

Vector Field Theory (E&M) Physis 4 Leture 2 Vetor Field Theory (E&M) Leture 2 Physis 4 Classial Mehanis II Otober 22nd, 2007 We now move from first-order salar field Lagrange densities to the equivalent form for a vetor field.

More information

Theoretical and Experimental Research on the Yoke of the Axial Piston Pumps

Theoretical and Experimental Research on the Yoke of the Axial Piston Pumps Theoretial and Experimental Researh on the Yoke of the Axial Piston Pumps Vasiliu Niolae, P.E., Ph.D., Professor University "Politehnia" of Buharest, Fluid Power Department Rosu Cristian Adrian, M.E.,

More information

Sensitivity Analysis in Markov Networks

Sensitivity Analysis in Markov Networks Sensitivity Analysis in Markov Networks Hei Chan and Adnan Darwihe Computer Siene Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwihe}@s.ula.edu Abstrat This paper explores

More information

Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion

Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six

More information

QUANTUM MECHANICS II PHYS 517. Solutions to Problem Set # 1

QUANTUM MECHANICS II PHYS 517. Solutions to Problem Set # 1 QUANTUM MECHANICS II PHYS 57 Solutions to Problem Set #. The hamiltonian for a lassial harmoni osillator an be written in many different forms, suh as use ω = k/m H = p m + kx H = P + Q hω a. Find a anonial

More information

Development of Accessible Hinge Details for Box-Girder Bridges

Development of Accessible Hinge Details for Box-Girder Bridges Development of Aessible Hinge Details for Box-Girder Bridges Kristopher Darnell University of Texas at Austin University of California Berkeley Dr. Khalid Mosalam paper 35 1 Seattle, Washington August

More information

A Spatiotemporal Approach to Passive Sound Source Localization

A Spatiotemporal Approach to Passive Sound Source Localization A Spatiotemporal Approah Passive Sound Soure Loalization Pasi Pertilä, Mikko Parviainen, Teemu Korhonen and Ari Visa Institute of Signal Proessing Tampere University of Tehnology, P.O.Box 553, FIN-330,

More information

Transient wave propagation analysis of a pantograph- catenary system

Transient wave propagation analysis of a pantograph- catenary system Journal of Physis: Conferene Series PAPER OPEN ACCESS Transient wave propagation analysis of a pantograph- atenary system To ite this artile: Kyohei Nagao and Arata Masuda 216 J. Phys.: Conf. Ser. 744

More information

General Closed-form Analytical Expressions of Air-gap Inductances for Surfacemounted Permanent Magnet and Induction Machines

General Closed-form Analytical Expressions of Air-gap Inductances for Surfacemounted Permanent Magnet and Induction Machines General Closed-form Analytial Expressions of Air-gap Indutanes for Surfaemounted Permanent Magnet and Indution Mahines Ronghai Qu, Member, IEEE Eletroni & Photoni Systems Tehnologies General Eletri Company

More information

WAVE-PARTICLE DUALITY: LIGHT

WAVE-PARTICLE DUALITY: LIGHT MISN-0-246 WAVE-PARTICLE DUALITY: LIGHT by E. H. Carlson WAVE-PARTICLE DUALITY: LIGHT PM r PM 1. The Problem Posed by Light a. Overview................................................ 1 b. Classial Partiles.......................................

More information

CHAPTER 26 The Special Theory of Relativity

CHAPTER 26 The Special Theory of Relativity CHAPTER 6 The Speial Theory of Relativity Units Galilean-Newtonian Relativity Postulates of the Speial Theory of Relativity Simultaneity Time Dilation and the Twin Paradox Length Contration Four-Dimensional

More information

MultiPhysics Analysis of Trapped Field in Multi-Layer YBCO Plates

MultiPhysics Analysis of Trapped Field in Multi-Layer YBCO Plates Exerpt from the Proeedings of the COMSOL Conferene 9 Boston MultiPhysis Analysis of Trapped Field in Multi-Layer YBCO Plates Philippe. Masson Advaned Magnet Lab *7 Main Street, Bldg. #4, Palm Bay, Fl-95,

More information

CSC2515 Winter 2015 Introduc3on to Machine Learning. Lecture 5: Clustering, mixture models, and EM

CSC2515 Winter 2015 Introduc3on to Machine Learning. Lecture 5: Clustering, mixture models, and EM CSC2515 Winter 2015 Introdu3on to Mahine Learning Leture 5: Clustering, mixture models, and EM All leture slides will be available as.pdf on the ourse website: http://www.s.toronto.edu/~urtasun/ourses/csc2515/

More information

Combined Electric and Magnetic Dipoles for Mesoband Radiation, Part 2

Combined Electric and Magnetic Dipoles for Mesoband Radiation, Part 2 Sensor and Simulation Notes Note 53 3 May 8 Combined Eletri and Magneti Dipoles for Mesoband Radiation, Part Carl E. Baum University of New Mexio Department of Eletrial and Computer Engineering Albuquerque

More information