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

Size: px
Start display at page:

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

Transcription

1 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 an Inna Sharf an Arun K. Misra 3 Mechanical Engineering Department, McGill University, Montreal, QCH3A C3, Canaa inna.sharf@mcgill.ca, 3 arun.misra@mcgill.ca This paper presents strategies for point-to-point reactionless manipulation of a satellite mounte ual-arm robotic system for capture of tumbling orbiting objects, such as out-ofcommission satellites an space ebris. Use of the ual-arm robot coul be more effective than the single arm when there is no provision for a grapple fixture or the object is tumbling. The ual arms can also provie exterous manipulation. As the main objective in capture of orbital objects is to move the en-effector from initial position to the grapple point with esire velocity, the task-level reactionless constraints in terms of en-effector velocities are erive. The trajectory planne using these constraints, however, results in several singular points within the robot s workspace. In orer to overcome this shortcoming, three point-to-point path planning strategies are presente, which improve the reactionless operation of the ual-arm robot. The strategies are illustrate by carrying out simulations for a 6-egree-of-freeom (DOF) ual-arm robotic system mounte on a satellite. corresponing author

2 I. Introuction Autonomous on-orbit services, such as capture of orbiting objects, an refueling, repair an refurbishment of isable satellites, using a robot mounte on a service satellite will be one of the important space operations in the future [-3]. Particularly, active ebris removal [3] has gaine a lot of attention in the recent years ue to the increase in the number of ebris. While performing these on-orbit services, it is esirable that the fuel consumption to overcome any attitue isturbance of the base satellite is negligible, as limite fuel is mainly reserve for orbital transfer maneuvers. The main motivation behin capturing space ebris is to avoi their possible collision with a working satellite in the same orbit. Some research efforts have been irecte towars achieving capture of a satellite with zero attitue isturbances to the base. This is also commonly referre to as reactionless manipulation. Such reactionless manipulation of a robot mounte on a satellite can translate into fuel savings an increase the operating life of the servicing system. In this regar, an optimization technique was propose in [4] for minimization of the base reactions, which i not lea to a satisfactory result on reactionless manipulation. The concept of isturbance map was propose in [5], which provie the minimum attitue isturbance paths leaing to lower fuel consumption. Later, a more effective extene isturbance map was also propose in [6]. These maps inee help in reucing attitue isturbance, but may not lea to nil attitue isturbance. A unique esign of manipulator was provie in [7] that le to reactionless manipulation. However, the major isavantage of such a esign is that it makes the robot architecture very complex. A reaction null space (RNS) base path planning was introuce in [8] for reactionless manipulation of a robot with flexible base. Later, in [9] the RNS base approach was use for control of space robots. It was shown in [8-9] that the joint velocities, obtaine using the RNS formulation, result in zero reaction moments on the satellite-base. Subsequently, the flight valiation with ETS-VII space robot an its extension to a kinematically reunant arm, for the zero reaction maneuvers, were presente in []. It was shown in [] that the constraints for reactionless manipulation, in terms of joint velocities, can be augmente with the task-level constraints in orer to obtain esire motion of the en-effector. The RNS-base concept was also use in [] for capture of a tumbling satellite. Moreover, the RNS-base approach with backwar integration was attempte in [] for trajectory planning of -link robot in the approach phase. Reactionless motion of a space robot when capturing an unknown tumbling target was also presente in [3] base on the momentum conservation equation an the recursive least squares metho.

3 In the above mentione stuies, the focus was mainly on completing the task with a singlearm robot. However, when the orbiting object oes not have provision for grapple or is tumbling, the interception may be very ifficult. In such cases, interception using a ual-arm robotic system can be appealing as this will increase the probability of grasp in comparison to a single-arm robot. Use of a ual-arm robot was escribe in [4], where one arm trace a given path while the secon arm worke both to control the base attitue an to optimize the total operational torque of the system. The path planning of a planar ual-arm free-floating manipulator was also presente in [5]. In both [4] an [5], the secon arm was only involve in attitue control of the base satellite rather than capture of the object. An autonomous approach for berthing of a moving satellite with two robots using a flexible wrist mechanism an impeance control was propose in [6]. Capture of a spinning object by two flexible manipulators using hybri position/force control an vibration suppression control was also presente in [7]. A harware-in-the-loop simulation of the spatial ual-arm space robot was presente in [8]. In [6-8], the robot arms were manipulate without paying any attention to the base attitue. More recently, in [9], the coorinate motion planning of a spatial ual-arm space robot for target capture was presente. In that work also, capture of the target with both arms was carrie out without regar to the base attitue. It is worth noting that in the above cite works, either the ual arms capture the orbiting object regarless of the attitue of the base satellite or the one arm captures the object while the secon arm is specifically use to control the attitue of the satellite. Thus, the reactionless capture using a ual-arm robotic system is least explore in the literature. The reactionless trajectories have to be esigne such that the two arms start from ifferent initial configurations an intercept the esignate points on the object. It is also essential that both en-effectors capture the tumbling object with velocity equal to that of the contact points in orer to avoi any significant impact. Hence, the major challenge in capture of a tumbling object is point-to-point reactionless manipulation of the robot arms with esire en-effector velocities at the capture instant. This is a complex problem as the reactionless constraints are nonholonomic an the resulting path has many singular points within the robot s workspace []. In comparison to the single-arm manipulation, the case of ual-arm manipulation requires that the paths have to be planne such that both arms stay away from singular configurations, which themselves are path epenents. Motivate by these facts, the reactionless point-to-point path planning strategies are propose in this work for ual-arm robotic systems for capture of a tumbling object. In the

4 present work, the two arms are move so as to prouce zero net reaction moment in orer to capture the tumbling target. Compare to the single-arm scenario, reactionless manipulation of the ual-arm may bring forth some interesting aspects of path planning. The point-to-point motion is achieve by eriving the task-level reactionless constraints for ual-arm robotic system by projecting the joint-level reactionless constraints [8] using the Generalize Jacobian Matrix []. As the resulting equations are highly constraine an have many singular points within the robot s workspace, three point-to-point path planning strategies are also presente. Extension of reactionless path planning to a ual-arm robotic system an strategies to avoi singular configurations are the main contributions of this work. Numerical illustrations are provie using a 6-DOF ual-arm robotic system. The rest of the paper is organize as follows: Some mathematical preliminaries are provie in Section II. The reactionless manipulation of the ual-arm robotic system is presente in Section III, whereas several point-to-point reactionless path planning strategies are evelope in Section IV. Finally, conclusions are given in Section V. II. Preliminaries an Notation In this section the equations of motion for ual-arm robotic systems an constraints for reactionless manipulation are erive as an extension of those given in [] for a single-arm robotic system. A. Equations of motion The equations of motion for an n-dof robot mounte on a floating-base can be written similarly to those in [] as: T Hb Hbm && xb cb Fb J b T + = + T F Hbm Hm φ&& cm τm Jm h () where H b R 6 6 an H m R n n are the inertia matrices of the base an manipulator, respectively, H bm R 6 n is the coupling inertia matrix, x&& b R 6 is the vector of linear an angular accelerations of the base, φ&& R n is the vector of joint accelerations, c b R 6 an c m R n are the velocity epenent nonlinear terms of the base an manipulator, F b an F h R 6 are the vectors of force an moment exerte on the centroi of the base an en-effector, respectively, τ m R n is the manipulator joint torque, J b R 6 6 an J m R 6 n are the Jacobian

5 matrices for the base an manipulator. For the ual-arm robotic system, H bm, H m, φ&, c m, τ m, J b, J m an F h can be written as follows: H O φ&& c H [ H H ] H φ&& c τm Jb Jm O Fh τm =, Jb =, J m = an Fh = τm Jb O Jm Fh m m m bm = bm bm, m =, =, m =, O H m φ&& m cm (a) where the subscripts an represent arm- an arm-, respectively. The floating-base moule of the Recursive Dynamics Simulator (ReDySim) [], evelope base on (-), is use for the simulations carrie out in this paper. B. Constraints for reactionless manipulation The matrices H b an H bm, of (), may also be expresse as T Hb _ v H b _ c Hbm _ v H b = an Hbm = Hb _ c Hb _ ω Hbm _ ω (3a) Now, the constraints for reactionless manipulation can be written using those in [] as H % φ & = where H % = H H H H (3b) - bm bm bm _ ω b_ v b_ c bm _ v Equation (3b) represents the nonholonomic constraints, where φ& belongs to the subspace of reactionless motion of the robot. The above equation forms the constraint leaing to zero reaction moments but non-zero reaction forces. It is assume that the satellite has reaction jets or thrusters to take care of the base reaction forces. As in the previous literature, hereafter, the term reactionless manipulation will imply manipulation with zero base moments only. For the ual-arm robotic system, (3b) can also be expresse as φ& m H% bm H % bm = or H% bmφ& m+ H% bmφ& m = φ& m (4) It can be seen from (4) that for the ual-arm robot the sum of the coupling angular momenta of both arms, not of the iniviual, has to be zero. Solution of (4) for φ& man φ& m can be obtaine using the pseuo inverse approach [8]. It is note that (3b-4) represent the constraints as a function of joint velocities, an will be referre to as joint-level reactionless constraints.

6 III. Reactionless Manipulation of Dual-Arm Robot In this section, the task-level constraints for reactionless manipulation are obtaine first, followe by a numerical example an iscussion on singularity. A. Task-level constraints for reactionless manipulation Since (4) is represente in terms of joint velocities, it cannot reaily provie a solution for the reactionless motion of the en-effectors in the Cartesian space. In such cases, one may use the constraint augmentation technique propose in [] for a single-arm robotic system. In the present work the point-to-point path planning for the en-effectors of the ual-arm robot is carrie out using task-level reactionless constraints. These constraints are obtaine by projecting the constraints in (4) using the Generalize Jacobian Matrix (GJM) []. For the ual-arm robot, the en-effectors velocities are relate to joint velocities by x& = J φ& + J φ &, where J = J J H H an J = J H H (5) - - h g m g m g m b b bm g b b bm x& = J φ& + J φ &, where J = J J H H an J = J H H - - h g m g m g m b b bm g b b bm (6) where J g an J g are the Generalize Jacobian matrices of the iniviual arms an, respectively, an J g an J g are the coupling Jacobians. The task-level constraints are obtaine by substituting expressions for φ && m an φ&& m from (5-6) into (4) as: x& h H% bh H % bh = or H% bhx& h+ H% bhx& h= x& h (7) In the above, J g J g H% bh H% bh = H% bm H% bm J g an J g = (8) J g J g where J g may be interprete as GJM of the ual-arm system. We assume that J g is invertible; a pseuo inverse can be use otherwise. Equation (7) represents the task-level reactionless constraints an will be use for point-topoint reactionless manipulation of the ual-arm robot. The egree-of-reunancy (DOR), f r, associate with (7) is given by the ifference in the number of columns an rows % %. Solution of (7) belongs to the f r -imensional subspace of R n. In other wors, of H bh H bh f r en-effectors velocities, out of x& h an x& h, can be controlle or prescribe inepenently. Solution of (7) can be obtaine using the pseuo inverse approach as:

7 + x& = ( E H % H % ) ξ & (9) h bh bh T T T where x& h =[ x& h x& h], H % bh = [ H % bh H % bh ], ( + E H% ) bhh % bh is the null-space projector an ξ & is an arbitrary velocity vector mappe into the null space. Alternatively, one may use coorinate partitioning as: i x h H% bh H% & bh = i x& h () where i x& h an x& h are the inepenent or free an epenent velocities. In (), the number of inepenent velocities is equal to f r, whereas the epenent velocities are obtaine as x& = H% H% x& () i i h bh bh h Equations (-) are preferre over (9), as the former allow one to irectly prescribe f r en-effector velocities inepenently. For example, for the 6-link planar ual-arm robot, shown in Fig., the linear an angular velocities of the en-effectors of arm-, i.e., (v Ex, v Ey, ω E ), an arm-, i.e., (v Ex, v Ey, ω E ), are the six variables in the task-space. The DOR associate with (7) is five, an hence, any five out of the above six can be chosen inepenently. On the other han, solution of () can be obtaine, provie Otherwise, () will lea to a singularity as illustrate in section III-C. H % is invertible. bh.5.5 φ 3 Arm- (,.5) Y (m).5 Satellite (Base) φ φ (.3,.) Orbiting object (,.5) Grapple points φ 4 φ 5 φ 6 (.7, -.9) Arm X (m) Figure. A ual-arm robotic system mounte on a satellite

8 B. A numerical example In orer to illustrate the use of (-), a planar ual-arm robotic system mounte on a satellite, as shown in Fig., is consiere. Each arm is comprise of three rigi links an 3- DOF. The centers-of-mass of the satellite an orbiting object lie at ( m, m), an ( m, m), respectively. The ual arms are initially in a non-symmetric configuration, as in practice it is not possible to achieve perfect symmetry. The points to be grapple on the object are also shown in Fig.. For this particular example, the object is assume to be stationary relative to satellite. The geometric an inertia properties of the ual arms an satellite are given in Table I. TABLE I. MODEL PARAMETERS OF THE DUAL-ARM AND SATELLITE Satellite Arm Arm Orbital () object Mass (Kg) 5 - Length (m) m ia. I zz (Kg.m ) Equations (-) are use to obtain the reactionless motion in the Cartesian space. As the ual-arm robot has DOR equal to five, the five velocities, v Ex, v Ey, v Ex, v Ey an ω E, are i chosen as free or inepenent, i.e., x&, while ω E is chosen as the epenent velocity, i.e., h The inepenent velocities are esigne using a fourth orer interpolating polynomial as follows: 3 4 t t t x& ( t) = a+ b + 4c + 5 () T T T where a= x&, b= e (6x& + 4 x& ), c= 5 e+ (8x& + 7 x& ), = 6 e (3x& + 3 x& ), an I I F I F I F e =[ x x ]/ T. Moreover, (x I an x F ) an ( x& an x& ) are the initial an final positions an F I I velocities, respectively. Equation () ensures zero initial an final acceleration. The initial an final positions are summarize in Table II, while the initial velocities of the en-effectors are assume to be zero. In orer to avoi high impact forces, the final velocities are set to match those of the grapple points, which are zeros in this particular example. With the above, the epenent velocity leaing to reactionless manipulation is obtaine using (). F x& h. TABLE II. INITIAL AND FINAL CONFIGURATIONS OF THE DUAL-ARM p Ex (m) p Ey (m) p Ex(m) p Ey (m) φ E (ra) t = s T = s Note: Initial orientation of the en-effector was taken as taken as.47 ra

9 Simulation of the ual-arm robotic system moving to the grapple points on the orbiting object was carrie out for the time perio of sec. The torque inputs following Proportional an Derivative (PD) control law given by (3) were use for carrying out the ynamic simulations: τ = K ( φ φ ) + K ( φ& φ & ), =, (3) mk p k ak k ak k In (3), K p an K are the iagonal matrices of proportional an erivative gains, an φ& k vector an φ& are the vectors of esire an actual joint rates, respectively, of the k th arm. The ak φ& k is obtaine using (5-6) from alreay esigne en-effectors velocities in (-). The iagonal elements of K p an K are obtaine empirically an are assume here to be 49 an 6, respectively. The inverse ynamics control for contacting robot, as propose in [- 3], can also be use, however, (3) is use here as ual-arm in approach phase is not subjecte to any external constraints. The paths travelle by the en-effectors are shown in Fig. by the otte curves. Figure also shows the final configuration of the ual-arm where the final positions of the eneffectors match with those provie in Table II. The angular velocity of the base, as epicte in Fig. 3, is of O( -6 ), proving reactionless manipulation of the ual-arm. Moreover, the final velocities of both en-effectors are zero as shown in Fig. 4. To provie partial valiation of the results obtaine using the ReDySim, the linear an angular momenta are also plotte in Fig. 5. The change in linear an angular momenta is of O( -7 ), confirming conservation of momenta by the ynamics simulator. x -6 Y (m) ω z - - X (m) Figure. Motion in X-Y plane -4 5 Figure 3. Angular velocity of base

10 Linear velocity (m/s).. v Ex v Ey v Ex v Ey 5 5 x -3 ω E ω E 5 (a) Linear velocity (b) Angular velocity Figure 4. Velocity of the en-effectors of the ual-arm Angular velocity (ra/s) 5 L x L y L z P x P y P z x -7 x -7 5 L (Nm.s) P (Nm) (a) Angular momentum (b) Linear momentum Figure 5. Angular an linear momenta C. Illustration of singularity Equation () is highly constraine an the resulting path can have many singular points within the robot s workspace. It is worth noting that these singular points are in general ifferent from those associate with ynamic singularities [4] an occur when the configuration epenent matrix H %, employe in (), becomes non-invertible. Therefore, it bh is important to etect this singularity so that appropriate corrective action can be taken. It is interesting to observe that the point-to-point path planne using (-) can be nonsingular for one set of esire final velocities an singular for another set. This is ue to the nonholonomic nature of the path planner. We illustrate this by attempting to fin a solution for the same example as consiere before, but with the object now spinning at constant angular rate of. ra/s, resulting in the velocities of the two grapple points of -. an. m/s in the X-irection. In this scenario, the path planner an simulation faile at t=6. sec, because of the singularity associate with singularity efine as: H %. Figure 6 plots the commonly use measure of matrix bh

11 T S= et( H% bhh % bh ) (4) while Fig. 7 shows that the torque at joint 3 is increasing rapily as the singularity is nearing. Given that it is not possible to ensure a non-singular path with the planner efine by (), we next propose alternative point-to-point strategies which can avoi such singular configurations while capturing a tumbling object. Singularity inex (S) Joint torques (Nm) 5-5 τ τ τ Time(s) Figure 6. Singularity inex 4 6 Time(s) Figure 7. Joint torques of arm- IV. Point-to-Point Reactionless Path Planning Strategies In this section three strategies are presente to alleviate the singularity problem in achieving reactionless manipulation of the ual-arm robotic system for capture of the tumbling object. A. Reactionless manipulation through choice of inepenent variables For the example consiere in III-C, the number of inepenent variables use in the path planner was chosen equal to DOR, i.e., 5. As objective in point-to-point path planning is mainly to capture the tumbling object at grapple points with esire velocity, one may alternatively choose only the translational task-level velocities v Ex, v Ey, v Ex an v Ey as the inepenent variables, leaving ω E an ω E to be the epenent variables. The compromise on the trajectory of ω E makes () reunant an provies enhance workspace. This is acceptable if the orientation of the en-effectors falls within the allowable limit at the time of capture. Therefore, the epenent velocities are calculate as: + i i + x& = H % H % x& + ( E H % H % ) ξ & (5) h bh bh h bh bh The example consiere in III-C is now simulate by computing the epenent velocities accoring to (5) with ξ& =. The resulting plots are illustrate in Figs. 8 an 9 an show that the en-effectors are able to reach the points on the spinning object with the require

12 velocities. The angular velocity of base is also shown in Fig., which is of O( -6 ) confirming reactionless manipulation of the ual-arm. However, the path travelle by arm- traverses the object as shown in Fig., which in a real setting woul result in collision. In such situations, one may chose alternative value of ξ &. Another possible solution is to move the ual-arm to an intermeiate point before it captures the object, as escribe in the following section..3 v Ex v Ey v Ex v Ey Y (m) Velocity (m/s).. - X (m) Figure 8. Motion in X-Y plane -. 5 Figure 9. En-effector s velocity 6 x -6 ω z 4 Y (m) 5 Figure. Angular velocity of base - - X (m) Figure.Snapshots of arm- B. Reactionless manipulation through intermeiate points In this approach, the intermeiate points between initial an final points are consiere when the en-effectors are not able to reach the grapple points accoring to (5) ue to either singular configuration or collision. Then, we propose that the motion is planne as follows: + i i + x& hj = H% bh H% bhx& hj + ( E H% bh H% bh ) ξ & j, j=,..., p+ (6) where x& hj an i x& hj are the epenent an esire inepenent velocities from (j-) th to j th point, an p is total number of intermeiate points. It may be note that i x& hj are esigne using (), which allows one to prescribe initial an final velocities for each intermeiate segment.

13 Equation (6) enables moving from the initial to intermeiate point an intermeiate to next point in a reactionless manner. The proceure for the propose piecewise path planner is illustrate by simulating the example of IV-A, however, with the en effectors of arm- an - commane to travel via intermeiate points (.,.8) an (.5, -.3), respectively. Moreover, the en-effector velocities of the two arms in the X-irection are graually increase through values of -.8m/s an.8 m/s at the intermeiate points. Figures an 3 show the simulation results, which epict that the en-effectors are able to reach the points on the spinning object with the require velocity. The ual-arm is able to move in a reactionless manner as the change in angular velocity of base is of O( -5 ) as shown in Fig. 4. The upper arm (arm-) is also able to avoi a collision, as epicte in Fig.5. For collision avoiance a metho base on repulsion potential fiel [5] can also be applie, which however is beyon the scope of the paper. v Ex v Ey v Ex v Ey.3 Y (m) Velocity (m/s).. - X (m) Figure. Motion in X-Y plane -. 5 Figure 3. En-effector s velocity x -5 ω z - Y (m) Figure 4. Angular velocity of base - - X (m) Figure 5. Snapshots of arm- It can be observe from Fig. 5 that the joint angle φ (efine in Fig. ) of arm-traverses through a value φ =8, which correspons to a ynamic singularity. Aitional insight is

14 gleane by consiering, singularity inices for H % an J bh g (the GJM of arm-) plotte in Figs. 6 an 7, respectively. It is seen that the peaks in Fig. 6 correspon to singular points in Fig. 7. Such situations are not esirable as they can result in high joint torque requirements. Even though for the problem uner stuy, the joint torques shown in Fig. 8 remaine in a reasonable range, high joint torques were observe when ifferent intermeiate points were chosen. The concept of hybri manipulation is propose next, which is foun to be more avantageous when such ynamic singularities are unavoiable. τ τ τ 3 Singularity inex (S) Singularity inex (S) Joint torques (Nm) Time(s) Time(s) Time(s) Figure 6. Singularity inex S( H% bh ) Figure 7. Singularity inex S(J g ) Figure 8.Joint torques for arm- C. Hybri reactionless manipulation In this approach, the path planning is carrie out using both joint- an task-level constraints for reactionless manipulation, erive in (4) an (7), respectively. The joint-level constraints are mainly use to cross the singular configurations, whereas the task-level constraints are use to reach the esire points. Equation (5) is use for the task-level path planning, an the joint-level path planning is carrie out similarly to (5) as In (7), + i i + φ& = H % H % φ& + ( E H % H % ) ξ & (7) h bm bm h bm bm H % bm is not function of the GJM, an hence, (7) is unaffecte by ynamic singularity. The concept of hybri reactionless manipulation is illustrate using the example solve in IV-B, however, with the motion of the ual arm robot ivie into three stages. In stage, the en-effectors travel form initial points to the intermeiate points, illustrate in Table III, using the task-level planner of Eq. (5) The intermeiate points are selecte closer to the object such that that φ or 8. Note that the angle φ was -.47 ra in the beginning an ra at the en of stage. In stage, the angle φ is brought from to.533 ra using joint-level planner of Eq. (7). In the thir an final stage, the eneffectors travel from current to final position using task-level planner.

15 TABLE III. INITIAL AND FINAL CONFIGURATOIN FOR THE VARIOUS STAGES IN HYBRID MANIPULATION p Ex (m) p Ey (m) p Ex(m) p Ey (m) Stage t= t= φ (ra) φ 4 (ra) φ 4 (ra) φ 5 (ra) φ 6 (ra) Stage t= t= p Ex (m) p Ey (m) p Ex(m) p Ey (m) Stage 3 t= T=.5.5 Note:. Joint angles correspon to the en-effector s position at the en of stage. En-effector s position correspon to the joint angles at the en of stage Figure 9 shows that the en-effectors are able to reach the points on the spinning object. The en-effectors are also able to match the velocities of grapple points as epicte in Fig.. Figure epicts that the angular velocity of the base is of O( -5 ), which ensure reactionless manipulation of the ual-arm. The singularity inex for stages -3 is plotte in Fig., which illustrate that H % is non-singular in stages an 3, whereas bh This shows the effectiveness of the hybri reactionless manipulation. H % is non-singular in stage. bm.4 v Ex v Ey v Ex v Ey Y (m) Velocity (m/s). - X (m) Figure 9. Motion in X-Y plane 5 Figure. En-effector s velocity ω z (ra/s) x Figure. Angular velocity of base Singularity inex (S) 8 6 Stage Stage Stage3 S( H % ) bh S( H % ) bm 5 Time(s) Figure. Singularity inex

16 V. Conclusions The point-to-point reactionless path planning for interception of the orbiting object using a ual-arm robotic system has been presente. It is shown that the reactionless paths contain singular points in the transformation between inepenent an epenent task velocities. To remey this problem, three strategies, namely, through a choice of fewer than maximum inepenent variables, reactionless manipulation through intermeiate points an hybri reactionless manipulation, have been presente. By changing a subset of inepenent variables one can exploit the reunancy in the ual-arm robot to enhance the reactionless manipulation. On the other han, reactionless manipulation through intermeiate points allows one to obtain an alternative non-singular or collision-free path, whereas the hybri reactionless manipulation helps in overcoming the ynamic singularity, which cannot be completely avoie using the former two approaches. The above three path planning strategies are simple yet effective for capturing the tumbling target using a ual-arm system in a reactionless manner. Even though the above methos allow to move the ual-arm in a reactionless manner, the following issues nee to be aresse: ) evelopment of a systematic approach for selection of the number of intermeiate points an their locations, ) calculation of optimum path out of the multiple reactionless paths resulting out of the point-to-point strategies, e.g., the path giving favourable orientation of the en-effector or minimum energy consumption, 3) use of inverse ynamics control for trajectory planning. These will be carrie out as future work. References [] F. Sellmaier, T. Boge, J. Spurmann, S. Gully, T. Rupp an F. Huber, On-Orbit Servicing Missions: Challenges an Solutions for Spacecraft Operations, AIAA SpaceOps Conference,, pp [] On-Orbit Satellite Servicing Stuy, NASA Project Report, ( [3] J. Liou, An active ebris removal parametric stuy for LEO environment remeiation, Avances in Space Research, 47(), ,. [4] C. L. Chung, S. Desa, an C. W. esilva, Base reaction optimization of reunant manipulators for space applications, The Robotic Institute CMU-RI-TR, pp. 88-7, 988.

17 [5] S. Dubowsky, an M. A. Torres, Path planning for space manipulators to minimize spacecraft attitue isturbances, IEEE International Conference on Robotics an Automation, 99, pp [6] M. A. Torres an S. Dubowsky, Minimizing spacecraft attitue isturbances in space manipulator systems, Journal of Guiance, Control, an Dynamics, vol. 5, no. 4, pp. -7, 99. [7] E. Papaopoulos, an A. Abu-Abe, Design an motion planning for a zero-reaction manipulator, IEEE International Conference on Robotics an Automation, 994,pp [8] D. N. Nenchev, K. Yoshia, an M. Uchiyama, Reaction null-space base control of flexible structure mounte manipulating systems, IEEE CDC, 996, Kobe, Japan, pp [9] D. N. Nenchev an K. Yoshia, Impact analysis an post-impact motion control issues of a free-floating space robot subject to a force impulse, IEEE Transactions on Robotics an Automation, vol. 5, no. 3, pp , 999. [] K. Yoshia, K. Hashizume an S. Abiko, Zero reaction maneuver: Flight valiation with ETS-VII space robot an extension to kinematically reunant arm, IEEE International Conference on Robotics an Automation,, pp [] K. Yoshia, D. Dimitrov an H. Nakanishi, On the capture of tumbling satellite by a space robot, IEEE International Conference on Intelligent Robots an Systems, 6, pp [] P. Piersigilli, I. Sharfan A. K. Misra, Reactionless capture of a satellite by a two egree-of-freeom manipulator, Acta Astronautica, vol. 66, no. (-), pp. 83-9,. [3] T. Nguyen-Huynh an I. Sharf, Aaptive reactionless motion for space manipulator when capturing an unknown tumbling target, IEEE International Conference on Robotics an Automation,, pp [4] K. Yoshia, R. Kurazumean Y. Umetani, Dual arm coorination in space free-flying robot, IEEE International Conference on Robotics an Automation, 99, pp [5] S. K. Agrawal an S. Shirumalla, Planning motions of a ual-arm free-floating manipulator keeping the base inertially fixe, Mechanism an Machine Theory, vol. 3, no., pp. 59-7, 995.

18 [6] A. Ejiri, I. Watanabe, K. Okabayashi, M. Hashima, M. Tatewaki, T. Aoki, an T. Maruyama, Satellite berthing experiment with a two-arme space robot, IEEE International Conference on Robotics an Automation,994,pp [7] T. Miyabe, A. Konno, an M. Uchiyama, Automate object capturing with a two-arm flexible manipulator, IEEE International Conference on Robotics an Automation, 3, pp [8] R. Takahashi, H. Ise, A. Konno, M. Uchiyama an D. Sato, Hybri simulation of a ual-arm space robot colliing with a floating object, IEEE International Conference on Robotics an Automation, 8, pp. -6. [9] W. Xu, Y. Liu, an Y.Xu, The coorinate motion planning of a ual-arm space robot for target capturing, Robotica, vol. 3, no. 5, pp ,. [] Y. Umetani an K. Yoshia, Resolve motion rate control of space manipulators with generalize jacobian matrix, IEEE Transactions on Robotics an Automation, vol. 5, no. 3, pp , 989. [] S. V. Shah, P. V. Nanhial an S. K. Saha, Recursive Dynamics Simulator (ReDySim)- A Multiboy Dynamics Solver, Theoretical an Applie Mechanics Letters, (6),pp. 63:-6,. [] L. Righetti, M. Mistry, J. Buchli an S. Schaal, Inverse ynamics control of floatingbase robots with external constraints: an unifie view, IEEE International Conference on Robotics an Automation, pp. 85-9,. [3] E. Papaopoulos an S. Dubowsky, Dynamic singularities in free-floating space manipulators, Journal of Dynamic Systems, Measurement an Control, vol. 5, no., pp. 44-5, 993. [4] N. Mansar, A Deicate Solver for Fast Operational-Space Inverse Dynamics, IEEE International Conference on Robotics an Automation, pp ,. [5] A. Dietrich, T. Wimböck, H. Täubig, A. Albu-Schäffer, an G. Hirzinger, Extensions to Reactive Self-Collision Avoiance for Torque an Position Controlle Humanois, IEEE International Conference on Robotics an Automation, pp ,.

Energy Optimum Reactionless Path Planning for Capture of Tumbling Orbiting Objects using a Dual-Arm Robot

Energy Optimum Reactionless Path Planning for Capture of Tumbling Orbiting Objects using a Dual-Arm Robot Energy Optimum Reactionless Path Planning for Capture of Tumbling Orbiting Objects using a Dual-Arm Robot S. V. Shah and A. Gattupalli Robotics Research Center International Institute of Information Technology

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

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

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

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

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

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 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

An Analysis of the Operational Space Control of Robots

An Analysis of the Operational Space Control of Robots EEE nternational Conference on Robotics an Automation Anchorage Convention District May 3-8,, Anchorage, Alaska, USA An Analysis of the Operational Space Control of Robots Ngoc Dung Vuong, Marcelo H Ang

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

Zero Reaction Maneuver: Flight Validation with ETS-VII Space Robot and Extension to Kinematically Redundant Arm

Zero Reaction Maneuver: Flight Validation with ETS-VII Space Robot and Extension to Kinematically Redundant Arm Proceedings of the 2001 IEEE International Conference on Robotics & Automation Seoul, Korea May 21-26, 2001 Zero Reaction Maneuver: Flight Validation with ETS-VII Space Robot and Extension to Kinematically

More information

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

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

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

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

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

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

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

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

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

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

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

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

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 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

Vibration Suppression Control of a Space Robot with Flexible Appendage based on Simple Dynamic Model*

Vibration Suppression Control of a Space Robot with Flexible Appendage based on Simple Dynamic Model* 213 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS) November 3-7, 213 Tokyo, Japan Vibration Suppression Control of a Space Robot with Flexible Appendage based on Simple Dynamic

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

The Exact Form and General Integrating Factors

The Exact Form and General Integrating Factors 7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily

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

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June

More information

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

Control of Generalized Contact Motion and Force in Physical Human-Robot Interaction 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

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

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

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

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

More information

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

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

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

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

On the Enumeration of Double-Base Chains with Applications to Elliptic Curve Cryptography

On the Enumeration of Double-Base Chains with Applications to Elliptic Curve Cryptography On the Enumeration of Double-Base Chains with Applications to Elliptic Curve Cryptography Christophe Doche Department of Computing Macquarie University, Australia christophe.oche@mq.eu.au. Abstract. The

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

4. Important theorems in quantum mechanics

4. Important theorems in quantum mechanics TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this

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

An Explicit Formulation of Singularity-Free Dynamic Equations of Mechanical Systems in Lagrangian Form Part Two: Multibody Systems

An Explicit Formulation of Singularity-Free Dynamic Equations of Mechanical Systems in Lagrangian Form Part Two: Multibody Systems Moeng, Ientification an Control, Vol 33, No 2, 2012, pp 61 68 An Expcit Formulation of Singularity-Free Dynamic Equations of Mechanical Systems in Lagrangian Form Part Two: Multiboy Systems Pål Johan From

More information

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

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

More information

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

A NONLINEAR SOURCE SEPARATION APPROACH FOR THE NICOLSKY-EISENMAN MODEL

A NONLINEAR SOURCE SEPARATION APPROACH FOR THE NICOLSKY-EISENMAN MODEL 6th European Signal Processing Conference EUSIPCO 28, Lausanne, Switzerlan, August 25-29, 28, copyright by EURASIP A NONLINEAR SOURCE SEPARATION APPROACH FOR THE NICOLSKY-EISENMAN MODEL Leonaro Tomazeli

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

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

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

IN the recent past, the use of vertical take-off and landing

IN the recent past, the use of vertical take-off and landing IEEE TRANSACTIONS ON ROBOTICS, VOL. 27, NO. 1, FEBRUARY 2011 129 Aaptive Position Tracking of VTOL UAVs Anrew Roberts, Stuent Member, IEEE, an Abelhami Tayebi, Senior Member, IEEE Abstract An aaptive position-tracking

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

Angles-Only Orbit Determination Copyright 2006 Michel Santos Page 1

Angles-Only Orbit Determination Copyright 2006 Michel Santos Page 1 Angles-Only Orbit Determination Copyright 6 Michel Santos Page 1 Abstract This ocument presents a re-erivation of the Gauss an Laplace Angles-Only Methos for Initial Orbit Determination. It keeps close

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

AN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A

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

More information

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

Least-Squares Regression on Sparse Spaces

Least-Squares Regression on Sparse Spaces Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction

More information

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com

More information

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

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

More information

Mechatronics by Analogy and Application to Legged Locomotion

Mechatronics by Analogy and Application to Legged Locomotion Mechatronics by Analogy an Application to Legge Locomotion by Victor Ragusila A thesis submitte in conformity with the requirements for the egree of Doctor of Philosophy Grauate Department of Aerospace

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

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

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

More information

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

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

Adaptive Optimal Path Following for High Wind Flights

Adaptive Optimal Path Following for High Wind Flights Milano (Italy) August - September, 11 Aaptive Optimal Path Following for High Win Flights Ashwini Ratnoo P.B. Sujit Mangal Kothari Postoctoral Fellow, Department of Aerospace Engineering, Technion-Israel

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

739. Design of adaptive sliding mode control for spherical robot based on MR fluid actuator

739. Design of adaptive sliding mode control for spherical robot based on MR fluid actuator 739. Design of aaptive sliing moe control for spherical robot base on MR flui actuator M. Yue,, B. Y. Liu School of Automotive Engineering, Dalian University of echnology 604, Dalian, Liaoning province,

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

Schrödinger s equation.

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

More information

u t v t v t c a u t b a v t u t v t b a

u t v t v t c a u t b a v t u t v t b a Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying

More information

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const.

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const. G4003 Avance Mechanics 1 The Noether theorem We alreay saw that if q is a cyclic variable, the associate conjugate momentum is conserve, q = 0 p q = const. (1) This is the simplest incarnation of Noether

More information

Cascaded redundancy reduction

Cascaded redundancy reduction Network: Comput. Neural Syst. 9 (1998) 73 84. Printe in the UK PII: S0954-898X(98)88342-5 Cascae reunancy reuction Virginia R e Sa an Geoffrey E Hinton Department of Computer Science, University of Toronto,

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

EVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION OF UNIVARIATE TAYLOR SERIES

EVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION OF UNIVARIATE TAYLOR SERIES MATHEMATICS OF COMPUTATION Volume 69, Number 231, Pages 1117 1130 S 0025-5718(00)01120-0 Article electronically publishe on February 17, 2000 EVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION

More information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische

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

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

Attitude Control System Design of UAV Guo Li1, a, Xiaoliang Lv2, b, Yongqing Zeng3, c

Attitude Control System Design of UAV Guo Li1, a, Xiaoliang Lv2, b, Yongqing Zeng3, c 4th National Conference on Electrical, Electronics an Computer Engineering (NCEECE 205) Attitue Control ystem Design of UAV Guo Li, a, Xiaoliang Lv2, b, Yongqing eng3, c Automation chool, University of

More information

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

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

More information

A General Formulation of Under-Actuated Manipulator Systems

A General Formulation of Under-Actuated Manipulator Systems A General Formulation of Under-Actuated Manipulator Systems Kazuya Yoshida and Dragomir N. Nenchev Dept. of Aeronautics and Space Engineering, Tohoku University Aobayama campus, Sendai 98-77 (Japan) yoshida@astro.mech.tohoku.ac.jp

More information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

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

More information

Hybrid Fusion for Biometrics: Combining Score-level and Decision-level Fusion

Hybrid Fusion for Biometrics: Combining Score-level and Decision-level Fusion Hybri Fusion for Biometrics: Combining Score-level an Decision-level Fusion Qian Tao Raymon Velhuis Signals an Systems Group, University of Twente Postbus 217, 7500AE Enschee, the Netherlans {q.tao,r.n.j.velhuis}@ewi.utwente.nl

More information

On the design of reactionless 3-DOF planar parallel mechanisms

On the design of reactionless 3-DOF planar parallel mechanisms Mechanism and Machine Theory 41 (2006) 70 82 Mechanism and Machine Theory www.elsevier.com/locate/mechmt On the design of reactionless 3-DOF planar parallel mechanisms Abbas Fattah *, Sunil K. Agrawal

More information

State-Space Model for a Multi-Machine System

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

More information

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

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

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

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

More information

Deriving ARX Models for Synchronous Generators

Deriving ARX Models for Synchronous Generators Deriving AR Moels for Synchronous Generators Yangkun u, Stuent Member, IEEE, Zhixin Miao, Senior Member, IEEE, Lingling Fan, Senior Member, IEEE Abstract Parameter ientification of a synchronous generator

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

A New Backstepping Sliding Mode Guidance Law Considering Control Loop Dynamics

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

More information

Pitch Attitude Control of a Booster Rocket

Pitch Attitude Control of a Booster Rocket Research Inventy: International Journal Of Engineering An Science Issn: 2278-4721, Vol.2, Issue 7 (March 2013), Pp 8-12 Www.Researchinventy.Com Pitch Attitue Control of a Booster Rocket Mr. Debabrata Roy

More information

Power Generation and Distribution via Distributed Coordination Control

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

More information

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

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

More information

Pure Further Mathematics 1. Revision Notes

Pure Further Mathematics 1. Revision Notes Pure Further Mathematics Revision Notes June 20 2 FP JUNE 20 SDB Further Pure Complex Numbers... 3 Definitions an arithmetical operations... 3 Complex conjugate... 3 Properties... 3 Complex number plane,

More information

RELATIVE POSITION ESTIMATION FOR INTERVENTION-CAPABLE AUVS BY FUSING VISION AND INERTIAL MEASUREMENTS

RELATIVE POSITION ESTIMATION FOR INTERVENTION-CAPABLE AUVS BY FUSING VISION AND INERTIAL MEASUREMENTS RELATIVE POSITION ESTIMATION FOR INTERVENTION-CAPABLE AUVS BY FUSING VISION AND INERTIAL MEASUREMENTS Anreas Huster an Stephen M. Rock Aerospace Robotics Lab, Stanfor University {huster,rock}@arl.stanfor.eu

More information

. Using a multinomial model gives us the following equation for P d. , with respect to same length term sequences.

. Using a multinomial model gives us the following equation for P d. , with respect to same length term sequences. S 63 Lecture 8 2/2/26 Lecturer Lillian Lee Scribes Peter Babinski, Davi Lin Basic Language Moeling Approach I. Special ase of LM-base Approach a. Recap of Formulas an Terms b. Fixing θ? c. About that Multinomial

More information

A Second Time Dimension, Hidden in Plain Sight

A Second Time Dimension, Hidden in Plain Sight A Secon Time Dimension, Hien in Plain Sight Brett A Collins. In this paper I postulate the existence of a secon time imension, making five imensions, three space imensions an two time imensions. I will

More information

Prep 1. Oregon State University PH 213 Spring Term Suggested finish date: Monday, April 9

Prep 1. Oregon State University PH 213 Spring Term Suggested finish date: Monday, April 9 Oregon State University PH 213 Spring Term 2018 Prep 1 Suggeste finish ate: Monay, April 9 The formats (type, length, scope) of these Prep problems have been purposely create to closely parallel those

More information

7.1 Support Vector Machine

7.1 Support Vector Machine 67577 Intro. to Machine Learning Fall semester, 006/7 Lecture 7: Support Vector Machines an Kernel Functions II Lecturer: Amnon Shashua Scribe: Amnon Shashua 7. Support Vector Machine We return now to

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