Control of Generalized Contact Motion and Force in Physical Human-Robot Interaction

Size: px
Start display at page:

Download "Control of Generalized Contact Motion and Force in Physical Human-Robot Interaction"

Transcription

1 15 IEEE International Conference on Robotics an Automation ICRA Washington State Convention Center Seattle, Washington, May -3, 15 Control of Generalize Contact Motion an Force in Physical Human-Robot Interaction Emanuele Magrini Fabrizio Flacco Alessanro De Luca Abstract During human-robot interaction tasks, a human may physically touch a robot an engage in a collaboration phase with exchange of contact forces an/or requiring coorinate motion of a common contact point. Uner the premise of keeping the interaction safe, the robot controller shoul impose a esire motion/force behavior at the contact or explicitly regulate the contact forces. Since intentional contacts may occur anywhere along the robot structure, the ability of controlling generalize contact motion an force becomes an essential robot feature. In our recent work, we have shown how to estimate contact forces without an explicit force sensing evice, relying on resiual signals to etect contact an on the use of an external epth sensor to localize the contact point. Base on this result, we introuce two control schemes that generalize the impeance an irect force control paraigms to a generic contact location on the robot, making use of the estimate contact forces. The issue of human-robot task compatibility is pointe out in case of control of generalize contact forces. Experimental results are presente for a KUKA LWR robot using a Kinect sensor. I. INTRODUCTION In the robotics community at large, there is great excitement about the recent possibility of realizing safe physical collaboration between human users an a new inustrial generation of lightweight, compliant, an frienly robots [1], []. This technical achievement is the outcome of avances in mechanics, actuation, sensing, an control that have emerge in the last few years, from statements of basic requirements see, e.g., [3] [] to novel research results see, e.g., [7], [] for the original collision etection algorithm implemente in the KUKA LWR an iiwa robots. Inee, more of such scientific evelopments an technological transfers are neee in orer to improve further the autonomy, safety, an ease of use of robotic coworkers in inustrial settings an of personal assistants in service applications. Along this path, we have propose a control architecture for safer physical human-robot collaboration, base on a hierarchy of three consistent robot behaviors [9]. In this framework, safety is the innermost layer with the highest priority of execution: it implies the ability of etecting an fast reacting to acciental collisions [1] an istinguishing between these an intentional contacts [11], [1]; coexistence is the robot capability of sharing workspaces, without interfering with humans: it calls for collision avoiance features [13], [1], an robot operative flexibility for The authors are with the Dipartimento i Ingegneria Informatica, Automatica e Gestionale, Sapienza Università i Roma, Via Ariosto 5, 15 Roma, Italy {magrini,fflacco,eluca}@iag.uniroma1.it. This work is supporte by the European Commission, within the FP7 ICT-7513 SAPHARI project moifying on line, abanoning, an resuming the originally planne tasks; finally, collaboration requires the robot to perform complex tasks involving irect physical interaction an coorination with humans. One can consier ifferent possible collaborative tasks where physical human-robot interaction occurs: carrying a heavy payloa, haning over objects, manipulating limp materials, holing firmly a part while the partner is operating on it, an so on. In most cases, regulation of the forces exchange at the contact is neee an/or accurate control of relative motion while keeping in place the original contact. As aitional convenient features, whole-boy manipulation shoul be possible, with single or multiple contacts occurring anywhere along the robot structure, an a minimum amount of extra sensing shoul be requeste, in particular without the nee of a F/T sensor at the wrist or of istribute tactile sensing for the measurement of contact/interaction forces. In [15], a metho has been evelope for estimating the interaction forces in irection an intensity uring ynamic contacts between a robot an a human, with no explicit force/torque sensing. The metho is base on the integrate use of moel-base resiual signals that etect the occurrence of collisions [] an of one or more external RGB-D sensors e.g., Kinects to approximately localize the contact point on the surface of the robot links. This virtual force sensor was shown to perform with reasonably goo estimation accuracy. Its output has been use alreay in a series of humanrobot interaction experiments base on an amittance control scheme efine at the contact point. In this paper, we illustrate further the use of the estimate contact forces for the esign of two human-robot interaction controllers that exten the impeance control metho [1] an the irect force control metho [17], respectively to the case of generic contact locations on the robot. For this reason, we collectively enote this approach as generalize contact motion an force control. While irect force control nees always a measure or an estimate of the contact force, stanar joint or Cartesian impeance control laws require a force measure only if the natural robot inertia at the proper level has to be change to a esire one. When the robot inertia is left unchange, force sensing is no longer require the term compliance control shoul be use then. The same situation hols true also for our generalize contact impeance control laws, simply replacing actual measurements with the use of force estimates see also Fig. 1. When realizing a force control law at the contact level, task compatibility issues may arise, mainly ue to the /15/$ IEEE 9

2 k F Fc k " k F m F c k " "#$%&"#$%&'$ '$$%*+,'--,./$*$'*,.* F c F ˆ c '"$%%&#,'-"#$%&'$ 1-,-3,.$-,./$*$'*'5 1-$*7"&7,',3-$,'-&-3,.$ '*%+,#$&"#$%&'$ '$$%*91.*-*$'*,. " k Fig. 1. A schematic comparison of impeance control schemes for a fully rigi robot arm, when a esire mass-spring-amper behavior is impose at the joint-level, at the Cartesian-level, or at the propose contact-level, in response to a contact force F c. The estimate contact force F b c can be obtaine with the metho in [15]. generality of the possible locations of the contact point. A first reason follows from the arguments alreay presente in [15]. When contacts occur on a link close to the robot base, the estimation metho may have too little information to recover a complete estimate of the applie contact force in particular, the contact Jacobian may not be full rank. The secon reason is even more funamental, an refers to the type of engage human-robot contact an its compatibility with the specification of a esire contact force. A complete specification of three inepenent scalar references for the contact force may be incompatible with the way the human is pushing on the robot. The situation is similar to efining a correct task frame in hybri force/motion control [1], [19, Chap. 9. 7]: the robot cannot control forces along contact irections where the human is not proviing a reaction. While the contact amittance control scheme in [15] specifies joint velocity references in response to estimate contact forces, an relies upon low-level servo loops for their execution, the two control esigns presente in this paper generate torque commans, an assume thus a torquecontrolle robot viz., motor currents can be impose. In aition, an accurate robot ynamic moel will be neee for control implementation. All evelopments in the paper are presente for rigi robot manipulators. Accoringly, no joint torque sensors are use nor neee an the virtual force sensor concept an the associate contact force/impeance controllers can be implemente using only the joint encoers as proprioceptive sensing. On the other han, the presente control experiments have been performe on a KUKA LWR where the use of harmonic rive introuces joint elasticity []. In practice, however, the control implementation for this robot is fully transparent with respect to the presence of compliant joints. The user efine torque-controlle moality of the FastResearchInterface FRI [1] implements the torque controller originally propose in [], an allows to use the available joint torque sensors both for resiual computation an for control. Moreover, for all moel-base computations we have use our own ientifie ynamic moel of the LWR link ynamics [3]. The paper is organize as follows. Section II sets the notation an inclues a short summary on the contact force estimation metho [15] use in the remainer of the paper. Section III recaps in a compact way the stanar impeance an irect force control laws for a single mass subject to an external contact force. The generalization of these control schemes to hanle contact situations occurring at a generic point of a multi-of robot manipulator is etaile in Sec. IV. Section V reports on several experimental results obtaine with the propose control methos uring physical humanrobot interaction with a KUKA LWR robot an a Kinect sensor. II. PRELIMINARIES For a robot manipulator having n serial joints, with associate generalize coorinates q R n, assume that a motion/force interaction task is efine for a contact point x c R 3 whose irect kinematics map is given by x c = fq. Differentiating this twice w.r.t. time yiels ẍ c = J c q q + J c q q, 1 where J c = fq/ q is the 3 n contact Jacobian matrix. Note that if the contact point is locate on link k k = 1,..., n, the last n k columns of J c will be ientically zero. In most cases, the robot is reunant respect to the given interaction task, being k > 3. At a given robot state q, q, all joint accelerations associate to a esire acceleration ẍ c of the contact point can be written as q = J # c ẍ c J c q + P c q, where J # c is a n 3 generalize inverse of the contact Jacobian, q R n is an arbitrary joint acceleration, an P c = I J # c J c is a projector in the null space of J c. The ynamic moel of a rigi manipulator interacting with the human environment at a robot point x c is given by Mq q + Cq, q q + gq = τ + J T c qf c, 3 with positive efinite inertia matrix M, Coriolis an centrifugal velocity terms C q factorize by the Christtoffel symbols, gravitational terms g, control torque τ R n, an joint torque τ c = J T c qf c resulting from the contact interaction force F c R 3. We will use also the compact notation nq. q = Cq, q q + gq. In the framework of operational space control [1], it has been shown that using the inertia-weighte pseuoinverse as the generalize inverse in guarantees consistency of 99

3 the transformation of task forces in the robot ynamics. The inertia-weighte pseuoinverse of contact Jacobian is J # c,m = M 1 J T c J c M 1 J T c 1. In the following, to simplify notation, each pseuoinverse J # c is to be intene as the inertia-weighte one given by. A. Contact force estimation The resiual vector r R n for a robot with ynamics 3 is efine as [7] t rt = K I p τ + C T q, q q gq + r s, 5 where p = Mq q is the generalize momentum of the robot an K I > is a iagonal gain matrix. The ynamic evolution of r has the stable, first-orer filter structure ṙ = K I τ c r. Thus, for sufficiently large gains, we can assume that r τ c = J T c qf c. Equation forms the basis for the estimations of the unknown contact force F c R 3. Using external sensing, the location of the contact point x c can be etermine, an thus the associate contact Jacobian J c can be compute. Depening on which link in the kinematic chain is involve in the contact, may consist of a square, uner-, or overetermine linear system. In any case, the contact force is estimate by pseuoinversion as F c = J T c q #r. 7 Inee, the estimate F c will be limite only to those components of F c that can be etecte by the resiual r, namely those contact forces that o not belong to the null space N J T c q. For further etails an for the analysis of cases when the contact Jacobian is not full rank, see [15]. III. IMPEDANCE AND FORCE CONTROL FOR A SINGLE MASS As an introuction to the control esign for the general multi-of robotic case, consier first a very simple 1-of example. As shown in Fig. 3, a single mass m moves on a frictionless horizontal plane uner the action of a control force f an of a contact force f c. The ynamic equation is mẍ = f + f c. Fig. 3. A mass m subject to a control force f an a contact force f c. A. Impeance control The inverse ynamics control law f = ma f c, 9 transforms system into the ouble integrator ẍ = a. 1 The auxiliary input a has to be esigne so that the controlle mass m, uner the action of the contact force f c, matches the behavior of an impeance moel characterize by a esire apparent mass m >, esire amping k >, an esire stiffness k p >, all with respect to a smooth motion reference x t, or m ẍ ẍ + k ẋ ẋ + k p x x = f c. 11 This is obtaine by using the control force f = m m ẍ + k ẋ ẋ + k p x x+ 1 f c. m m 1 For k < k pm, the secon-orer linear system 11 is characterize by a pair of asymptotically stable complex poles with natural frequency an amping ratio given by Fig.. A view of the result of the on-line estimation of the contact force ω n = kp k, ζ = m. 13 k p m Reucing the esire mass m, for given values of stiffness an amping, will increase both the angular frequency an the amping ratio, an thus improve transients. On the other han, for a given mass, an increase of the stiffness k p shoul be accompanie by an increase of the amping k in orer to prevent more oscillatory transients. If the esire mass equals the natural original mass, i.e., m = m, a measure of the contact force f c is no longer neee in the impeance controller 1 which is then also calle a compliance control law, since the main esign parameter left is the esire stiffness k p. In particular, for 3

4 regulation tasks with x t = x = constant the control law collapses just to a PD action on the position error e = x x, f = k p x x k ẋ. 1 The system will converge to x = x if an only if there is no contact force. With f c but constant, there will be convergence to the close-loop equilibrium x e satisfying k p x x e + f c = x e = x + f c k p. 15 B. Force control If we esire to control irectly the contact force to a esire possibly, constant value f, it is necessary to buil a force error e f = f f c into the control law. After using eq. 9, efine the auxiliary input a as a = 1 m k f f f c k ẋ 1 with force error gain k f > an velocity amping k >. The associate control force is then f = m m k f f f c k ẋ f c. 17 A contact force measure is always neee in this case, even if we choose m = m. The close-loop system becomes m ẍ + k ẋ = k f f f c, 1 which shows that the force error e f goes to zero if a constant contact force is applie. However, the actual position x = x e reache at the equilibrium is not known a priori, an will epen on the actual history of the external contact force. IV. GENERALIZED CONTACT MOTION AND FORCE CONTROL FOR ROBOTS We exten here the control esign results of Sec. III to a generic robot contact point x c of a robot manipulator, an consier also the use of the contact force estimation F c, as provie by 7, in place of the real measurement, each time this will be neee. A. Impeance control Starting from robot moel 3, an ropping from now on epenencies for the sake of compactness, consier the inverse ynamics control law τ = Ma + n J T c F c, 19 where F c is the estimate contact force. In ieal esign conitions, the above control law realizes a feeback linearization, leaing to a system of ouble integrators q = a. The auxiliary control input a is chosen so as to impose a mechanical impeance moel between the estimate contact forces F c an the motion of the contact point x c, or M ẍ c ẍ + K ẋ c ẋ + K p x c x = F c, being M > the esire inertia matrix, K > the esire amping matrix, K p > the esire stiffness matrix very often, these matrices are taken as iagonal, an x t the smooth esire trajectory of the contact point. Solving for ẍ c from an using, the control input a is obtaine as a=j # c M 1 M ẍ +K ė+k p e M J c q+ F c +P q, 1 with contact position error e = x x c. Thus, the resulting torque control input is τ = MJ # c M 1 M ẍ + K ė + K p e M J c q + F c + MP q + n J T F c c. In, observe that the contribution of F c is given by τ =... + MJ # c M 1 J T c F c. 3 As expecte, when choosing the esire inertia matrix M equal to the natural Cartesian inertia of the robot at the contact point, also the estimation of the contact forces F c will no longer be neee. At a given configuration q, the natural inertia of the robot at a generic contact point is 1 M = J c M 1 J T c. By using, it is easy to see that MJ # c M 1 J T c =. Whenever we esire to impose a ifferent inertial behavior at the contact, e.g., so as to achieve ifferent behaviors in ifferent irections of the operative space, we shall make use of the contact force estimate F c. B. Force control The impeance metho realizes only an inirect control of the forces exchange uring the interaction with the human, ynamic balancing with the esire position x c of the contact point. The evelopment of a irect controller for the contact force is crucial for collaborative tasks that nees accurate execution in uncertain conitions. With the reference to 19, choose in place of the auxiliary comman 1, the following acceleration a = J # c M 1 K f F F c K ẋ c M J c q +P q, 5 where F is the esire force at the contact point, the contact force error is e f = F F c, an K f > is the force error gain matrix. Replacing 5 in 19, the final torque control input is τ = MJ # c M 1 + MP q + n J T c F c. K f e f K ẋ c M J c q Similarly to 1, the close-loop system is escribe by M ẍ c + K ẋ c = K f F F c, 7 which shows again that the force error shoul go to zero whenever a constant contact force is applie. Nonetheless, the issue of task compatibility remains open. In the above expression of the control law for the human-robot contact forces, the specification of the esire force F is apparently an unconstraine one. We shall see in the experiments that 31

5 a careful choice of F shoul be mae so as to avoi inconsistent an rifting behavior of the robot, with serious consequences for safety. " V. E XPERIMENTAL RESULTS $" Fig.. The set-up an the Cartesian reference frame for the contact impeance control experiments. The color of the X, Y, an Z axes are the same use for the associate components in the following plots. 15 Resiual [Nm] 1 Similarly for the force control law it is 1 τ = MJ# K f ef K x c M J c q c M b c. M P K N q + C q J Tc F 5 r1 15 r r3 r r5 r r Fbc,x b c. M P K N q + C q J Tc F 5 1 Fbc,y Fbc,z.. 9. Error [m] c #" Dynamic human-robot interaction experiments using the propose controllers have been performe on a KUKA LWR. The workspace is monitore by a Microsoft Kinect epth sensor, positione at a istance of 1.7 m behin the robot an at a height of 1.1 m w.r.t. the robot base frame. The Kinect provies epth images at 3 Hz rate. All algorithms are execute on a qua-core CPU. The complete process of contact force estimation an motion an/or force control run on the robot at 5 ms cycle time. With reference to the impeance control law, efine x = xc tc R3 as the initial an constant position of the contact point when the interaction with the human begins at t = tc. The esire velocity an acceleration are =. Due to the reunancy of this set then to zero, x = x robot with respect to many tasks, a null-space acceleration = K N q, with K N >, in vector has been chosen as q orer to boun/amp out self-motion movement of the arm. Finally, recall that the KUKA LWR has a built-in gravity compensation. Thus, the actual impeance control law use to comman the robot was bc τ = M J # M 1 K p e K x c M J c q + F All experiments presente in this paper are recore in the accompanying vieo clip. The plots can be even better appreciate when looking in parallel at the vieo.. ex ey ez A. Impeance control with natural contact inertia In the first experiment, the human pushes the robot at ifferent contact points an on ifferent links, as shown in the snapshot of Fig.. The control law is given by, with the esire inertia matrix at the contact chosen equal to. The other impeance parameters were set to K = I 3, K p = 5 I 3, an K N = I 7, where I k enotes the k k ientity matrix. Figure 5 shows the behavior of the resiual vector, of the contact force estimate, an of the Cartesian position error of the contact point with respect to its initial position x = xc tc. B. Impeance control with moifie contact inertia In this experiment, the impeance scheme was use again uring the human-robot interaction. However, the esire inertia matrix at the contact point was set to M = 3, 1 Fig. 5. Contact impeance control with esire inertia matrix at the contact equal to the natural one. Resiual vector components [top]. Estimate contact force components [center]. Contact position error components [bottom]. i.e., ifferent from the natural one which was also couple. Non-uniform values were chosen on the iagonal, so as to obtain ifferent behaviors along the three Cartesian axes see also Fig.. To enhance this effect, the human pushes the robot always on link, although in ifferent irections. The other impeance parameters use were the same as in Sec. Vb c are neee A. In this case, the estimate contact forces F by the control law in orer to obtain the esire mass-springamper system. Figure shows the same quantities of the previous impeance control experiment. When the contact is initially etecte, the robot will move the contact point with a ynamics that epens on the irection of the external force. When the han is remove, the robot brings back the contact point to the original initial position. As shown in Fig. thir 3

6 plot, the error returns closer to zero on the y axis, where M,y = 3. By increasing the esire inertia, the error natural frequency an amping ecrease, as inicate by 13, an transients will worsen see, e.g., the x axis at t = 15 s. " #" $" 3 Resiual [Nm] 1 1 r 1 r r 3 r r 5 r r Fig. 7. The set-up an the Cartesian reference frame for the contact force control experiments. On the left, the human is pushing along the Y irection where a non-zero esire force is specifie F y, = 15 N. On the right, a rift effect occurs ue the absence of a human reaction in the Y irection a force is applie approximately along the X irection only. Error [m] F c,x Fc,y Fc,z Fig.. Contact impeance control with moifie inertia matrix at the contact. Resiual vector components [top]. Estimate contact force components [center]. Contact position error components [bottom]. C. Contact force control with possible rift effects The purpose of this experiment was to show that it is highly not recommene to try to regulate human-robot contact forces at arbitrary values an in multiple fixe irections of the Cartesian space. When the human oes not push or resist along a irection where a non-zero force reference has been specifie, the robot typically rifts in space in the attempt to regulate the incompatible esire force. This coul be angerous because an unexpecte movement occur, through which the robot may collie with the human. With this in min, the parameters were set to K = 1.5 I 3, K f =.1 I 3, K N = 15 I 7, an the esire contact force was chosen as F = 15 T, i..e., only in the Y irection. Figure shows the behavior of the resiual vector r, of the contact force estimate F c, an of the Cartesian velocity ẋ c of the contact point. When the human pushes the robot on link along the Y axis, the robot reacts properly an regulate the force components to the esire level. However, when the human pushes the robot along the X irection, a rift occurs along the Y irection, as shown in Fig. e.g., at t = 1. s. D. Contact force control with task compatibility The contact force control scheme 9 was use in this interaction control experiment, where the human pushes the e x e y e z Resiual [Nm] Velocity [m/s] r 1 r r 3 r r 5 r r F c,x Fc,y Fc,z F,y ẋ c,x ẋ c,y ẋ c,z Fig.. Contact force control with rift motion ue to task incompatibility. Components of the resiual vector [top]. Moulus of the estimate contact force [center]. Components of the Cartesian contact velocity ẋ c [bottom]. When the human pushes along the Y irection, regulation occurs see the first 1 s. The human pushing along the X irection where F x, = will be accommoate by a robot retraction, whereas the absence of a human reaction in the Y irection where F y, = 15 N oes not allow force regulation an generates large velocity rifts in that same irection see the last s in the thir plot. robot successively on link 3, link 5, an then link see the associate behaviors of the resiuals in Fig. 9, [top]. To avoi incompatibility with the human, only the norm of the esire contact force F has been, here as a constant F = 15 N. However, F will change irection accoring to the human-robot contact type, an being always aligne with the estimate contact force vector which may now vary without any restriction. We will have thus, F,x = 15 F c,x F c, F,y = 15 F c,y F c, F,z = 15 F c,z F c. 33

7 The other force control parameters use were the same as in Sec. V-C. Figure 9 shows the same quantities of the previous force control experiment. Note that the controller starts operating after the first etection of a contact. When the human pushes the robot, it gets compliant in orer to regulate the contact force at the esire value. When the human tries to recee, the robot pushes against the human han so as to maintain the esire F c = F. The fact that the human may abanon the contact an then resume it in a ifferent point an along a ifferent irection is no longer a problem for this correctly efine contact force controller. On the other han, this force control law is able to regulate contact forces only in static conitions i.e., when ẋ c =. While the human is moving the robot uring a contact, some error will result as shown in Fig. 9. Resiual [Nm] Velocity [m/s] r 1 r r 3 r r 5 r r F c F ẋ c,x ẋ c,y ẋ c,z Fig. 9. Contact force control is here esigne for automatic task compatibility. Components of the resiual vector [top]. Moulus of the estimate contact force [center]. Components of the Cartesian contact velocity ẋ c [bottom]. VI. CONCLUSION We have expane the portfolio of controllers for physical human-robot interaction tasks, introucing the concept of motion an/or force control at any contact point on the robot. Generalizations of classical impeance an irect force control schemes have been introuce an implemente without the explicit nee of a force sensor, relying on a fast estimation of contact forces. Besie eveloping an refining further such basic robotic skills for safe collaboration, our next step will be their application to concrete an specific tasks e.g,, performing human-robot collaborative assembly within the inustrial use cases of the SAPHARI project. We also wish to exploit the robot reunancy for controlling motion in the plane perpenicular to the estimate contact force i.e., achieving hybri force-motion control. REFERENCES [1] KUKA presents the worl s first lightweight inustrial robot at Hanover Fair, April 1, 13. [Online]. Available: LBR iiwa worl premiere.htm [] ABB unveils collaborative robot: YuMi, September 9, 1. [Online]. Available: [3] O. Khatib, K. Yokoi, O. Brock, K.-S. Chang, an A. Casal, Robots in human environments: Basic autonomous capabilities, Int. J. of Robotics Research, vol. 1, no. 7, pp. 9, [] J. Heinzmann an A. Zelinsky, Quantitative safety guarantees for physical human-robot interaction, Int. J. of Robotics Research, vol., no. 7/, pp. 79 5, 3. [5] A. Bicchi an G. Tonietti, Fast an soft arm tactics: Dealing with the safety-performance trae-off in robot arms esign an control, IEEE Robotics an Automation Mag., vol. 11, no., pp. 33,. [] A. De Santis, B. Siciliano, A. De Luca, an A. Bicchi, An atlas of physical human-robot interaction, Mechanism an Machine Theory, vol. 3, no. 3, pp. 53 7,. [7] A. De Luca an R. Mattone, Sensorless robot collision etection an hybri force/motion control, in Proc. IEEE Int. Conf. on Robotics an Automation, 5, pp [] A. De Luca, A. Albu-Schäffer, S. Haain, an G. Hirzinger, Collision etection an safe reaction with the DLR-III lightweight robot arm, in Proc. IEEE/RSJ Int. Conf. on Intelligent Robots an Systems,, pp [9] A. De Luca an F. Flacco, Integrate control for phri: Collision avoiance, etection, reaction an collaboration, in Proc. IEEE Int. Conf. on Biome. Robotics an Biomechatronics, 1, pp. 95. [1] S. Haain, A. Albu-Schäffer, A. De Luca, an G. Hirzinger, Collision etection an reaction: A contribution to safe physical humanrobot interaction, in Proc. IEEE/RSJ Int. Conf. on Intelligent Robots an Systems,, pp [11] M. Eren an T. Tomiyama, Human-intent etection an physically interactive control of a robot without force sensors, IEEE Trans. on Robotics, vol., no., pp. 37 3, 1. [1] M. Geravan, F. Flacco, an A. De Luca, Human-robot physical interaction an collaboration using an inustrial robot with a close control architecture, in Proc. IEEE Int. Conf. on Robotics an Automation, 13, pp. 7. [13] D. Ebert an D. Henrich, Safe human-robot-cooperation: Imagebase collision etection for inustrial robots, in Proc. IEEE/RSJ Int. Conf. on Intelligent Robots an Systems,, pp. 39. [1] F. Flacco, T. Kröger, A. De Luca, an O. Khatib, A epth space approach to human-robot collision avoiance, in Proc. IEEE Int. Conf. on Robotics an Automation, 1, pp [15] E. Magrini, F. Flacco, an A. De Luca, Estimation of contact forces using a virtual force sensor, in Proc. IEEE/RSJ Int. Conf. on Intelligent Robots an Systems, 1, pp [1] N. Hogan, Impeance control: An approach to manipulation: Part I - Theory, Part II - Implementation, Part III - Applications, ASME J. of Dynamic Systems, Measurement, an Control, vol. 17, pp. 1, 195. [17] B. Siciliano an L. Villani, Robot Force Control. Kluwer, [1] O. Khatib, A unifie approach for motion an force control of robot manipulators: The operational space formulation, IEEE Trans. on Robotics an Automation, vol. 3, no. 1, pp. 3 53, 197. [19] B. Siciliano, L. Sciavicco, L. Villani, an G. Oriolo, Robotics: Moeling, Planning an Control, 3r e. Lonon: Springer,. [] A. De Luca an W. Book, Robots with flexible elements, in Springer Hanbook of Robotics, B. Siciliano an O. Khatib, Es. Springer,, pp [1] KUKA.FastResearchInterface 1., KUKA System Technology KST, Augsburg, D, 11, Version. [] A. Albu-Schäffer, C. Ott, an G. Hirzinger, A unifie passivitybase control framework for position, torque an impeance control of flexible joint robots, Int. J. of Robotics Research, vol., pp. 3 39, 7. [3] C. Gaz, F. Flacco, an A. De Luca, Ientifying the ynamic moel use by the KUKA LWR: A reverse engineering approach, in Proc. IEEE Int. Conf. on Robotics an Automation, 1, pp

Hybrid Force/Velocity Control for Physical Human-Robot Collaboration Tasks

Hybrid Force/Velocity Control for Physical Human-Robot Collaboration Tasks Hybri Force/Velocity Control for Physical Human-Robot Collaboration Tasks Emanuele Magrini Alessanro De Luca Abstract During human-robot collaboration tasks, we may physically touch the robot at a generic

More information

Numerical Integrator. Graphics

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

More information

On the Feedback Linearization of Robots with Variable Joint Stiffness

On the Feedback Linearization of Robots with Variable Joint Stiffness 8 IEEE International Conference on Robotics an Automation Pasaena, CA, USA, May 19-3, 8 On the Feeback Linearization of Robots with Variable Joint Stiffness G. Palli an C. Melchiorri Dipartimento i Elettronica,

More information

Combining Real and Virtual Sensors for Measuring Interaction Forces and Moments Acting on a Robot

Combining Real and Virtual Sensors for Measuring Interaction Forces and Moments Acting on a Robot 216 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS) Daejeon Convention Center October 9-14, 216, Daejeon, Korea Combining Real and Virtual Sensors for Measuring Interaction Forces

More information

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

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

More information

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

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

More information

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers

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

More information

Toward Torque Control of a KUKA LBR IIWA for Physical Human-Robot Interaction

Toward Torque Control of a KUKA LBR IIWA for Physical Human-Robot Interaction Toward Torque Control of a UA LBR IIWA for Physical Human-Robot Interaction Vinay Chawda and Günter Niemeyer Abstract In this paper we examine joint torque tracking as well as estimation of external torques

More information

Experimental Robustness Study of a Second-Order Sliding Mode Controller

Experimental Robustness Study of a Second-Order Sliding Mode Controller Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans

More information

Tutorial Test 5 2D welding robot

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

More information

Laplacian Cooperative Attitude Control of Multiple Rigid Bodies

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

More information

Chapter 2 Lagrangian Modeling

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

More information

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

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

More information

Short Intro to Coordinate Transformation

Short Intro to Coordinate Transformation Short Intro to Coorinate Transformation 1 A Vector A vector can basically be seen as an arrow in space pointing in a specific irection with a specific length. The following problem arises: How o we represent

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

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

More information

TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS

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

More information

PD Controller for Car-Following Models Based on Real Data

PD Controller for Car-Following Models Based on Real Data PD Controller for Car-Following Moels Base on Real Data Xiaopeng Fang, Hung A. Pham an Minoru Kobayashi Department of Mechanical Engineering Iowa State University, Ames, IA 5 Hona R&D The car following

More information

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

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

More information

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

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

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

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

More information

APPPHYS 217 Thursday 8 April 2010

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

More information

IROS 16 Workshop: The Mechatronics behind Force/Torque Controlled Robot Actuation Secrets & Challenges

IROS 16 Workshop: The Mechatronics behind Force/Torque Controlled Robot Actuation Secrets & Challenges Arne Wahrburg (*), 2016-10-14 Cartesian Contact Force and Torque Estimation for Redundant Manipulators IROS 16 Workshop: The Mechatronics behind Force/Torque Controlled Robot Actuation Secrets & Challenges

More information

Identifying the Dynamic Model Used by the KUKA LWR: A Reverse Engineering Approach

Identifying the Dynamic Model Used by the KUKA LWR: A Reverse Engineering Approach 4 IEEE International Conference on Robotics & Automation (ICRA Hong Kong Convention and Exhibition Center May 3 - June 7, 4 Hong Kong, China Identifying the Dynamic Model Used by the KUKA LWR: A Reverse

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

VIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK

VIRTUAL STRUCTURE BASED SPACECRAFT FORMATION CONTROL WITH FORMATION FEEDBACK AIAA Guiance, Navigation, an Control Conference an Exhibit 5-8 August, Monterey, California AIAA -9 VIRTUAL STRUCTURE BASED SPACECRAT ORMATION CONTROL WITH ORMATION EEDBACK Wei Ren Ranal W. Bear Department

More information

A Physically-Based Fault Detection and Isolation Method and Its Uses in Robot Manipulators

A Physically-Based Fault Detection and Isolation Method and Its Uses in Robot Manipulators des FA 4.13 Steuerung und Regelung von Robotern A Physically-Based Fault Detection and Isolation Method and Its Uses in Robot Manipulators Alessandro De Luca Dipartimento di Informatica e Sistemistica

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

Interpolated Rigid-Body Motions and Robotics

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

More information

Lecture 6: Control of Three-Phase Inverters

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

More information

Lecture «Robot Dynamics»: Dynamics and Control

Lecture «Robot Dynamics»: Dynamics and Control Lecture «Robot Dynamics»: Dynamics and Control 151-0851-00 V lecture: CAB G11 Tuesday 10:15 12:00, every week exercise: HG E1.2 Wednesday 8:15 10:00, according to schedule (about every 2nd week) Marco

More information

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics This moule is part of the Memobust Hanbook on Methoology of Moern Business Statistics 26 March 2014 Metho: Balance Sampling for Multi-Way Stratification Contents General section... 3 1. Summary... 3 2.

More information

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

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

More information

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

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

More information

Relative Position Sensing by Fusing Monocular Vision and Inertial Rate Sensors

Relative Position Sensing by Fusing Monocular Vision and Inertial Rate Sensors Proceeings of ICAR 2003 The 11th International Conference on Avance Robotics Coimbra, Portugal, June 30 - July 3, 2003 Relative Position Sensing by Fusing Monocular Vision an Inertial Rate Sensors Anreas

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

An analysis of the operational space control of robots*

An analysis of the operational space control of robots* SIMech technical reports (SR_V_N3 MCH) Volume Number 3 Jul-Sep An analysis of the operational space control of robots* N. D. Vuong, H. A. Marcelo,. M. Lim, an S. Y. Lim Abstract heoretically, the operational

More information

Distributed Force/Position Consensus Tracking of Networked Robotic Manipulators

Distributed Force/Position Consensus Tracking of Networked Robotic Manipulators 180 IEEE/CAA JOURNAL OF AUTOMATICA SINICA, VOL. 1, NO. 2, APRIL 2014 Distribute Force/Position Consensus Tracking of Networke Robotic Manipulators Lijiao Wang Bin Meng Abstract In this paper, we aress

More information

The Principle of Least Action

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

More information

ECE 422 Power System Operations & Planning 7 Transient Stability

ECE 422 Power System Operations & Planning 7 Transient Stability ECE 4 Power System Operations & Planning 7 Transient Stability Spring 5 Instructor: Kai Sun References Saaat s Chapter.5 ~. EPRI Tutorial s Chapter 7 Kunur s Chapter 3 Transient Stability The ability of

More information

q 1 F m d p q 2 Figure 1: An automated crane with the relevant kinematic and dynamic definitions.

q 1 F m d p q 2 Figure 1: An automated crane with the relevant kinematic and dynamic definitions. Robotics II March 7, 018 Exercise 1 An automated crane can be seen as a mechanical system with two degrees of freedom that moves along a horizontal rail subject to the actuation force F, and that transports

More information

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

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

More information

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

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

More information

Robust Tracking Control of Robot Manipulator Using Dissipativity Theory

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

More information

Nested Saturation with Guaranteed Real Poles 1

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

More information

Calculus of Variations

Calculus of Variations Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

Application of Nonlinear Control to a Collision Avoidance System

Application of Nonlinear Control to a Collision Avoidance System Application of Nonlinear Control to a Collision Avoiance System Pete Seiler Mechanical Eng. Dept. U. of California-Berkeley Berkeley, CA 947-74, USA Phone: -5-64-6933 Fax: -5-64-663 pseiler@euler.me.berkeley.edu

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Minimum-time constrained velocity planning

Minimum-time constrained velocity planning 7th Meiterranean Conference on Control & Automation Makeonia Palace, Thessaloniki, Greece June 4-6, 9 Minimum-time constraine velocity planning Gabriele Lini, Luca Consolini, Aurelio Piazzi Università

More information

6. Friction and viscosity in gasses

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

More information

Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory

Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory Design A Robust Power System Stabilizer on SMIB Using Lyapunov Theory Yin Li, Stuent Member, IEEE, Lingling Fan, Senior Member, IEEE Abstract This paper proposes a robust power system stabilizer (PSS)

More information

Human-friendly Motion Control of a Wheeled Inverted Pendulum by Reduced-order Disturbance Observer

Human-friendly Motion Control of a Wheeled Inverted Pendulum by Reduced-order Disturbance Observer 2008 IEEE International Conference on Robotics an Automation Pasaena, CA, USA, May 19-23, 2008 Human-frienly Motion Control of a Wheele Inverte Penulum by Reuce-orer Disturbance Observer DONGIL CHOI an

More information

Text S1: Simulation models and detailed method for early warning signal calculation

Text S1: Simulation models and detailed method for early warning signal calculation 1 Text S1: Simulation moels an etaile metho for early warning signal calculation Steven J. Lae, Thilo Gross Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Str. 38, 01187 Dresen, Germany

More information

Optimal LQR Control of Structures using Linear Modal Model

Optimal LQR Control of Structures using Linear Modal Model Optimal LQR Control of Structures using Linear Moal Moel I. Halperin,2, G. Agranovich an Y. Ribakov 2 Department of Electrical an Electronics Engineering 2 Department of Civil Engineering Faculty of Engineering,

More information

Exponential Tracking Control of Nonlinear Systems with Actuator Nonlinearity

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

More information

Chapter 6: Energy-Momentum Tensors

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

More information

Optimization of a point-mass walking model using direct collocation and sequential quadratic programming

Optimization of a point-mass walking model using direct collocation and sequential quadratic programming Optimization of a point-mass walking moel using irect collocation an sequential quaratic programming Chris Dembia June 5, 5 Telescoping actuator y Stance leg Point-mass boy m (x,y) Swing leg x Leg uring

More information

Calculus of Variations

Calculus of Variations 16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t

More information

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

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

More information

DYNAMIC PERFORMANCE OF RELUCTANCE SYNCHRONOUS MACHINES

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

More information

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann Gravitation as the result of the reintegration of migrate electrons an positrons to their atomic nuclei. Osvalo Domann oomann@yahoo.com (This paper is an extract of [6] liste in section Bibliography.)

More information

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

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

More information

The skeleton algorithm for self-collision avoidance of a humanoid manipulator

The skeleton algorithm for self-collision avoidance of a humanoid manipulator The skeleton algorithm for self-collision avoidance of a humanoid manipulator Agostino De Santis, Alin Albu Schäffer, Christian Ott, Bruno Siciliano, and Gerd Hirzinger Abstract For use in unstructured

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

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

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

More information

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1 Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems

More information

ELEC3114 Control Systems 1

ELEC3114 Control Systems 1 ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.

More information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method 1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse

More information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

More information

A Path Planning Method Using Cubic Spiral with Curvature Constraint

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

More information

Contact Distinction in Human-Robot Cooperation with Admittance Control

Contact Distinction in Human-Robot Cooperation with Admittance Control Contact Distinction in Human-Robot Cooperation with Admittance Control Alexandros Kouris, Fotios Dimeas and Nikos Aspragathos Robotics Group, Dept. of Mechanical Engineering & Aeronautics University of

More information

Lecture «Robot Dynamics»: Dynamics 2

Lecture «Robot Dynamics»: Dynamics 2 Lecture «Robot Dynamics»: Dynamics 2 151-0851-00 V lecture: CAB G11 Tuesday 10:15 12:00, every week exercise: HG E1.2 Wednesday 8:15 10:00, according to schedule (about every 2nd week) office hour: LEE

More information

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

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

More information

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the

More information

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

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

More information

Quantum mechanical approaches to the virial

Quantum mechanical approaches to the virial Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from

More information

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS SEQUENTIAL CIRCUITS CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents

More information

Both the ASME B and the draft VDI/VDE 2617 have strengths and

Both the ASME B and the draft VDI/VDE 2617 have strengths and Choosing Test Positions for Laser Tracker Evaluation an Future Stanars Development ala Muralikrishnan 1, Daniel Sawyer 1, Christopher lackburn 1, Steven Phillips 1, Craig Shakarji 1, E Morse 2, an Robert

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

Systems & Control Letters

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

More information

Some vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10

Some vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10 Some vector algebra an the generalize chain rule Ross Bannister Data Assimilation Research Centre University of Reaing UK Last upate 10/06/10 1. Introuction an notation As we shall see in these notes the

More information

Reactionless Path Planning Strategies for Capture of Tumbling Objects in Space Using a Dual-Arm Robotic System

Reactionless Path Planning Strategies for Capture of Tumbling Objects in Space Using a Dual-Arm Robotic System Reactionless Path Planning Strategies for Capture of Tumbling Objects in Space Using a Dual-Arm Robotic System Suril V. Shah Robotics Research Center, IIIT Hyeraba, Anhra Praesh, 53, Inia surilshah@iiit.ac.in

More information

Conservation Laws. Chapter Conservation of Energy

Conservation Laws. Chapter Conservation of Energy 20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action

More information

Practical implementation of Differential Flatness concept for quadrotor trajectory control

Practical implementation of Differential Flatness concept for quadrotor trajectory control Practical implementation of Differential Flatness concept for quarotor trajectory control Abhishek Manjunath 1 an Parwiner Singh Mehrok 2 Abstract This report ocuments how the concept of Differential Flatness

More information

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy, NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which

More information

for Articulated Robot Arms and Its Applications

for Articulated Robot Arms and Its Applications 141 Proceedings of the International Conference on Information and Automation, December 15-18, 25, Colombo, Sri Lanka. 1 Forcefree Control with Independent Compensation for Articulated Robot Arms and Its

More information

Manual Guidance of Humanoid Robots without Force Sensors: Preliminary Experiments with NAO

Manual Guidance of Humanoid Robots without Force Sensors: Preliminary Experiments with NAO Manual Guidance of Humanoid Robots without Force Sensors: Preliminary Experiments with NAO Marco Bellaccini Leonardo Lanari Antonio Paolillo Marilena Vendittelli Abstract In this paper we propose a method

More information

Invariant Extended Kalman Filter: Theory and application to a velocity-aided estimation problem

Invariant Extended Kalman Filter: Theory and application to a velocity-aided estimation problem Invariant Extene Kalman Filter: Theory an application to a velocity-aie estimation problem S. Bonnabel (Mines ParisTech) Joint work with P. Martin (Mines ParisTech) E. Salaun (Georgia Institute of Technology)

More information

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

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

More information

A Review of Multiple Try MCMC algorithms for Signal Processing

A Review of Multiple Try MCMC algorithms for Signal Processing A Review of Multiple Try MCMC algorithms for Signal Processing Luca Martino Image Processing Lab., Universitat e València (Spain) Universia Carlos III e Mari, Leganes (Spain) Abstract Many applications

More information

Predictive Control of a Laboratory Time Delay Process Experiment

Predictive Control of a Laboratory Time Delay Process Experiment Print ISSN:3 6; Online ISSN: 367-5357 DOI:0478/itc-03-0005 Preictive Control of a aboratory ime Delay Process Experiment S Enev Key Wors: Moel preictive control; time elay process; experimental results

More information

Torque OBJECTIVE INTRODUCTION APPARATUS THEORY

Torque OBJECTIVE INTRODUCTION APPARATUS THEORY Torque OBJECTIVE To verify the rotational an translational conitions for equilibrium. To etermine the center of ravity of a rii boy (meter stick). To apply the torque concept to the etermination of an

More information

Dynamics of the synchronous machine

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

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

Thermal runaway during blocking

Thermal runaway during blocking Thermal runaway uring blocking CES_stable CES ICES_stable ICES k 6.5 ma 13 6. 12 5.5 11 5. 1 4.5 9 4. 8 3.5 7 3. 6 2.5 5 2. 4 1.5 3 1. 2.5 1. 6 12 18 24 3 36 s Thermal runaway uring blocking Application

More information

Slide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13)

Slide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13) Slie10 Haykin Chapter 14: Neuroynamics (3r E. Chapter 13) CPSC 636-600 Instructor: Yoonsuck Choe Spring 2012 Neural Networks with Temporal Behavior Inclusion of feeback gives temporal characteristics to

More information

We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors

We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors We are IntechOpen, the worl s leaing publisher of Open Access books Built by scientists, for scientists 4,100 116,000 120M Open access books available International authors an eitors Downloas Our authors

More information

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Bair, faculty.uml.eu/cbair University of Massachusetts Lowell 1. Pre-Einstein Relativity - Einstein i not invent the concept of relativity,

More information

Dead Zone Model Based Adaptive Backstepping Control for a Class of Uncertain Saturated Systems

Dead Zone Model Based Adaptive Backstepping Control for a Class of Uncertain Saturated Systems Milano (Italy) August - September, 11 Dea Zone Moel Base Aaptive Backstepping Control for a Class of Uncertain Saturate Systems Seyye Hossein Mousavi Alireza Khayatian School of Electrical an Computer

More information