We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
|
|
- Marilyn Kelly
- 5 years ago
- Views:
Transcription
1 We ae IntechOpen, the wold s leading publishe of Open Access books Built by scientists, fo scientists 3,500 08,500.7 M Open access books available Intenational authos and editos Downloads Ou authos ae among the 54 Counties deliveed to OP % most cited scientists.% Contibutos fom top 500 univesities Selection of ou books indexed in the Book Citation Index in Web of Science Coe Collection (BKCI) Inteested in publishing with us? Contact book.depatment@intechopen.com Numbes displayed above ae based on latest data collected. Fo moe infomation visit
2 8 Contol of Redundant Robotic Manipulatos with State Constaints Miosław Galicki Open Access Database Intoduction Recently, edundant manipulatos have attacted much attention due to thei potential abilities which ae inteesting fom both a theoetical and pactical point of view. Redundant degees of feedom make it possible to pefom some useful objectives such as collision avoidance in the wok space with both static and moving obstacles, joint limit avoidance, and/o avoidance of singula configuations when the manipulato moves. Most of pactical tasks, fo example, inseting a shaft into co-opeating elements (beaing, sleeve, o atchet-wheel), equie the knowledge of geometic paths (given in the wok space) and a pope toleance of matching that specifies the coesponding accuacy of the path following. In many othe industial tasks such as lase cutting o ac welding, the accuacy of path following is vital, and it is easonable to assume that designes and manufactues will specify pecision using an absolute toleance on tacking eo. he application of edundant manipulatos to such tasks complicates thei pefomance, since these manipulatos in geneal do not povide unique solutions. Consequently, some objective citeia should be specified to solve the obot tasks uniquely. he minimization of pefomance time is mostly consideed in the liteatue. Seveal appoaches may be distinguished in this context. Using the concept of a egula tajectoy and the extended state space, the stuctue of path-constained time-optimal contols has been studied in the woks (Galicki, 998b Galicki, 000) fo kinematically edundant manipulatos. Moeove, the efficient numeical pocedues able to find such contols wee also poposed in the woks (Galicki, 998a Galicki & Pajak, 999). Nevetheless, these algoithms equie full knowledge of manipulato Jacobian matix and obot dynamic equations, too. Although all the afoementioned algoithms poduce optimal solutions, they ae not suitable to eal-time computations due to thei computational complexity. heefoe, it is natual to attempt othe techniques in ode to contol the obot in eal-time. Using on-line tajectoy time scaling, a dynamic and computed toque laws espectively, a nealy time-optimal path tacking contol fo Souce: Industial-Robotics-heoy-Modelling-Contol, ISBN , pp. 964, ARS/plV, Gemany, Decembe 006, Edited by: Sam Cubeo 499
3 500 Industial Robotics: heoy, Modelling and Contol non-edundant obotic manipulatos with patially uncetain dynamics has been pesented in woks (Dahl, 994 Kiefe et al., 997). Howeve, these algoithms equie the solution of invese kinematic poblem along the path. A technique which avoids solving an invese of obot kinematic equations and uses the exact Jacobian matix, has been offeed in (Galicki, 00 Galicki 004) fo detemining a collision-fee tajectoy of edundant manipulatos opeating in both a static envionment (Galicki, 00) and in a dynamic one (Galicki, 004). Recently, a genealized tanspose Jacobian contolle with gavity compensation and a non-linea (satuating) deivative tem has been intoduced in (Galicki, 006a) to geneate obot contols subject to geometic path and actuato constaints. As is known, many obotic contolles have been poposed to solve both a set point contol poblem (a egulation task) (akegaki & Aimoto, 98 Aimoto, 996 Canudas de Wit et al., 996 Sciavicco & Siciliano, 996 Aimoto, 990 Kelly, 999 Galicki, 00 Galicki, 005) and the tajectoy tacking (Slotine & Li, 987 Slotine &Li, 99 Feng & Palaniswami, 99 Beghuis et al., 993 Lewis et al., 993 omei, 000), espectively. Howeve, most of these contolles have assumed full knowledge of manipulato kinematic equations. Recently, seveal appoximate Jacobian setpoint contolles have been poposed (Cheach et al., 999 Yazael & Cheach, 00 Cheach et al., 003) to tackle uncetainties in both obot kinematics and dynamics. he contolles poposed do not equie the exact knowledge of Jacobian matix and dynamic equations. Howeve, the esults in (Cheach et al., 999 Yazael & Cheach, 00 Cheach et al., 003) ae applicable only to a setpoint contol of a obot. his pape, which is based on ou ecent wok (Galicki, 006b), intoduces a new class of adaptive path following contolles not equiing the full knowledge of both kinematic and dynamic equations in the contol laws. Consequently, they ae suitable fo contolling uncetain obotic manipulatos. Motivated in pat by the dissipativity and adaptivity methodology (Slotine & Li, 987), we develop path following contolles whose stuctue is composed of tanspose adaptive Jacobian contolle plus a non-linea tem including an estimated contol. Unde the assumption of the full ank adaptive Jacobian matix, the poposed contol scheme has been deived based on the Lyapunov stability theoy. By using sensoy feedback of the end-effecto position, it is also shown that the end-effecto is able to follow a pescibed geometic path fo obots with both uncetain kinematics and dynamics. Futhemoe, new adaptive laws extending the adaptive algoithm fom (Slotine & Li, 987) to tackle kinemetic uncetainties too, ae poposed. It is to notice that appoximate Jacobian setpoint contolles fom (Cheach et al., 999 Yazael & Cheach, 00 Cheach et al., 003) can not be diectly applicable to ou task. he eason is that appoximate Jacobian matix in (Cheach et al., 999 Yazael & Cheach, 00 Cheach et al., 003) does not include kinematic paametes to be adapted and the eo of appoximation is a'pioi bounded. On the othe
4 Contol of Redundant Robotic Manipulatos with State Constaints 50 hand, the contolle poposed in ou wok adaptively vaies both kinematic paametes of the Jacobian matix and the dynamic ones in such a way as to stably follow by the end-effecto a geometic path. he pape is oganized as follows. Section fomulates the obotic task to be accomplished in tems of a contol poblem. Section 3 descibes how to employ the Lyapunov stability theoy to detemine contols (if they exist). Section 4 povides us with a compute example of geneating the obot contols in a two dimensional task space fo a plana edundant manipulato compising thee evolute kinematic pais. Finally, some conclusions ae dawn.. Fomulation of the adaptive contol poblem he contol scheme designed in the next section is applicable to holonomic mechanical systems compising both non-edundant and edundant manipulatos consideed hee which ae descibed, in geneal, by the following dynamic equations, n expessed in genealized co-odinates (joint co-odinates) q = (q,...,q n ) R ( q) + C( q, ) + G( q) u M = () whee M (q) denotes the n n inetia matix C( q, ) is the n -dimensional vecto epesenting centifugal and Coiolis foces G (q) stands fo the n - dimensional vecto of gavity foces u ( u,..., ) = is the vecto of contols (toques/foces) n denotes the numbe of kinematic pais of the V-th class. In most applications of obotic manipulatos, a desied path fo the end-effecto is specified in task space such as visual space o Catesian space. he aim is to follow by the end-effecto a pescibed geometic path (given in the m - dimensional task space) descibed by the following equations p ( q) Θ( s) = 0 () whee p n m : R R denotes an m -dimensional, non-linea (with espect to vecto q ) mapping constucted fom the kinematic equations of the manipulato = ( p ( q),..., p ( q ) m ( s ) = ( Θ ( s),..., ( ) ) Θ p q) ) ( m s u n Θ stands fo a given geometic path s is the cuent paamete of the path (e.g. its length) s [ 0, smax ] s max is the maximal path length. he mapping Θ is assumed to be bounded togethe dθ with the fist and second deivatives and not degeneated, i.e. > 0. ds
5 50 Industial Robotics: heoy, Modelling and Contol It should be: he kinematic equations of a manipulato ae independent p p ank = m. In geneal (i.e. when ank m, we should equie that q q p p dθ ank = ank in ode to guaantee consistency of the obotic q q ds task (). he poblem is to detemine contol u which geneates manipulato tajectoy q = q(t) and path paameteization s = s(t) satisfying the equation () fo each t [ 0, ], whee denotes an (unknown) time hoizon of task pefomance. It is natual to assume that at the initial moment t = 0, fo which s ( 0) = 0, a given (by definition) initial configuation q ( 0) = q0 satisfies (), i.e. p ( q 0 ) Θ(0) = 0 (3) Final path paameteization fulfils the equality s ( ) smax = 0 (4) Futhemoe, at the initial and the final time moment, the manipulato and path velocities equal zeo, i.e. ( 0) = ( ) = 0 (5) and s ( 0) = s ( ) = 0 (6) As is known, task space velocity p is elated to joint space velocity as follows p = J ( q, Y ) (7) whee p J(q,Y) = q is the mn Jacobian matix Y stands fo an odeed set Y = Y,..., Y k such as link lengths, joint offsets k denotes the numbe of kinematic paametes. Seveal impotant popeties of dynamic equations () may be deived (Spong & Vidyasaga, 989) of kinematic paametes ( )
6 Contol of Redundant Robotic Manipulatos with State Constaints 503. he inetia matix M ( q ) is symmetic and positive definite fo all q R n. Matix M(q) C(q,q) is skew-symmetic so that n M(q) v,q, R v, C(q, v = 0 (8) whee, is the scala poduct of vectos. 3. he dynamic equations () ae linea with espect to an odeed set of physical paametes X = (X,...,X ), i.e. d M ( q) + C( q, ) + G( q) = D( q,,, )X (9) whee whee D( q,,, ) is called the ( d ) n dynamic egesso matix d stands fo the numbe of the physical paametes such as link masses and inetias. Diffeential equation (7) has the following popety. 4. he ight hand side of (7) is linea with espect to e Y. Consequently, equation (7) can be expessed as follows p = J ( q, Y ) = K( q, )Y (0) whee K( q,q) is called the (m k)kinematic egesso matix. In ode to simplify futhe computations s max is assumed to be equal to, i.e. s max =. Let us define eos e and em+ of path following (task eos) as e = e ( e,..., e ) = p( q) Θ ( s) = ( p Θ,..., p Θ ) = m + s m. m m () Many commecial sensos ae available fo measuement of end-effecto position p, such as vision systems, electomagnetic measuement systems, position sensitive detectos o lase tacking systems. Hence, path following eo e in () is also assumed to be available (fo a given s ) fom measuement. Fo evolute kinematic pais, consideed hee, mapping p(.) is bounded. Consequently, we have the following popety.
7 504 Industial Robotics: heoy, Modelling and Contol 5. Boundedness of mappings p( ),Θ( ), implies that task eo e is bounded. Expessions ()-(6) fomulate the obot task as a contol poblem. he fact that thee exist state equality constaints makes the solution of this poblem difficult. he next section will pesent an appoach that endes it possible to solve the contol poblem ()-(6) making use of the Lyapunov stability theoy. 3. Adaptive path contol of the manipulato Ou aim is to contol the manipulato such that the end-effecto fulfils ()-(6). heefoe, we popose adaptive Jacobian path following contolles fo obots with both uncetain kinematics and dynamics. In ou appoach, the exact knowledge of both obot dynamic equations and Jacobian matix is not equied in updating the uncetain paametes. In the pesence of kinematic uncetainty, the paametes of the Jacobian matix ae uncetain and hence equality (0) can be expessed as follows Jˆ = K( q, )Yˆ () whee Jˆ = J(q,Y) ˆ R mn is an adaptive Jacobian matix and Ŷ stands fo the vecto of estimated kinematic paametes. In ode to show the stability of the path following contol system in the pesence of both kinematic and dynamic uncetainties, we define an adaptive joint-space sliding vectoz as z = λjˆ e + (3) whee λ is a positive scala coefficient. he adaptivity of vecto z is undestood in the sense that the paametes of the adaptive Jacobian matix will be updated by a paamete update law, defined late. Diffeentiating equation (3) with espect to time yields z = Jˆ e + λjˆ e + λ (4) Based on equations () and (3)-(4), we can easily expess obot dynamic e- quations by means of vecto z and its time deivative, as ( q, ) z + M ( q) + C( q, ) q + G( q u M ( q) z + C ) = (5)
8 Contol of Redundant Robotic Manipulatos with State Constaints 505 ˆ Whee = ( ) = λj e J e and q ( q q ) = J e, n λ ˆ, n λ =,...,,,...,, Moeove, we know fom popety 3., that the last thee tems of equation (5) become linea with espect to vecto X and hence they can be witten as follows ˆ ( q, ) + G( q) D( q,,, ). M ( q) + C = X (6) Inseting the ight hand side of (6) into (5), obot dynamic equations (5) take the following fom ( q, ) z + D( q,, ) X u M ( q) z + C, = (7) Based on (3), (4) and (7), lets popose the following adaptive Jacobian contolle u = kz k P Jˆ e + D( q,,, )Xˆ (8) whee ˆX stands fo the estimated physical paamete vecto defined below k and k p ae some positive scalas which could be eplaced by diagonal matices of positive constants without affecting the stability esults obtained futhe on (this should lead to impoved pefomance). Hee, scala constants ae chosen fo simplicity of pesentation. In ode to measue eo e, path paameteization s= s(t) is equied which is computed by solving the following scala diffeential equation s = k s s k P dθ( s) dγ e, + γ ( s ) + s ds ds ( ) (9) whee k s denotes a positive coefficient γ is assumed to be a stictly positive ( ) 0 function { } inf γ > of s with bounded fist and second deivatives fo any γ0 s +, whee γ0 γ( 0) { } inf γ = (as will be seen futhe on, γ 0 detemines the uppe bound on the accuacy of the path following and may be specified by the use). he choice of function γ is cucial fo computational effectiveness of contol dγ scheme (8), (9). One possibility is γ = γ(s) with 0 ds >. An altenative choi-
9 506 Industial Robotics: heoy, Modelling and Contol ce could be γ γ( e ) =, whee γ attains its maximum fo e = 0 and smoothly deceases as e inceases. Fo simplicity of futhe consideations, we take the fist fom of γ. Assumption. Function γ is equied not to satisfy diffeential equation dγ γ + ( s ) = 0 ds As will be seen futhe on, Assumption. esults in an asymptotic convegence of s to. Let us note, that the fist two tems fom (8) pesent an adaptive tanspose Jacobian contolle with adaptively vaying kinematic paamete vecto Ŷ. he last tem in dependence (8) is an estimated contol based on equation (6). Es- Ĵ = J q,y ˆ timated kinematic paametes Ŷ of the adaptive Jacobian matix ( ) ae updated accoding to the following law Yˆ = w k k P K ( q, )e (0) and estimated physical paametes ˆX of the dynamic equations ae updated by Xˆ = w d D ( q,,, )z () whee wk, w d, ae, similaly as befoe, positive gains (scalas) which could be eplaced by diagonal matices of positive constants without affecting the stability of contolle (8). Although some kinematic paametes appea in ˆX, we should adapt on them sepaately inŷ to peseve lineaity. Estimated kinematic paametes Ŷ (updated accoding to ule (0)) ae then used to compute adaptive Jacobian Ĵ and its time deivative Ĵ using fo this pupose the ight hand side of (0) which is only a mapping of q,e, and. Having obtained the adaptive Jacobian matix and its time deivative, we detemine quantities and. It is woth noticing, that thei computation does not equie any pseudoinvese of matix Ĵ which esults in numeical stability of contolle (8). Finally, based on q,q,q,, and, we may detemine dy- D q,q,q,q, which is then used to update estimated namic egesso matix ( ) physical paamete vecto ˆX.
10 Contol of Redundant Robotic Manipulatos with State Constaints 507 Let us note, that setpoint contolles poposed in woks (Cheach et al., 999 Yazael & Cheach, 00 Cheach et al., 003), which ae computationally somewhat simple, can not be applicable to ou task. he eason is that the eo of appoximation in (Cheach et al., 999 Yazael & Cheach, 00 Cheach et al., 003) is bounded by a constant and appoximate Jacobian matix does not include paametes to be adapted. Due to adaptation law (0), one can not guaantee in ou task to satisfy assumption egading the appoximation eo made in woks (Cheach et al., 999 Yazael & Cheach, 00 Cheach et al., 003). he closed-loop eo dynamics is obtained by inseting the ight hand side of equation (8) into equation (7) M ( q) z = C ( q, ) z + D( q,,, ) dθ( s) e = J ( q, Y ) e m+ ds dθ( s) e m+ = kse m+ kp e, ds ~ Y = wkkpk ( q, ) e ~ X = w D ( q,,, )z d ~ X kz k P Jˆ e dγ + γ em+ + em+ () ds ~ = ˆ whee Y Y Y X X X deive the following esult. ~ = ˆ. Applying the Lyapunov stability theoy, we heoem. If thee exists a solution to the poblem ()-(6) and adaptive Jacobian matix Ĵ is non-singula along end-effecto path () and function γ fulfils Assumption., then contol scheme (8) geneates manipulato tajectoy whose limit point ( ( ), e m+ ( ), e( ), em+ ( ) ) = (0,0,0,0),i.e. satisfying state constaints ()-(6), is asymptotically stable. Poof. Conside a Lyapunov function candidate V ~ ~,. Mz z + X + Y + e m+ + kpγem + kpe wd w (3) k = + he time deivative of V is given by
11 508 Industial Robotics: heoy, Modelling and Contol V = z, Mz + e M z, z + k P e, + w ~ ~ X, X dθ + e m+ + kp e, e m+ + kpγem+ e m+ + k ds J + w dγ e ds ~ ~ Y, Y + m P m+ em +. d k (4) Substituting Mz,.e,.e m+,y and X fom V fo the ight-hand sides of closedloop eo dynamics () and using the skew-symmetic popety of matix M C [popety. eqn. (8)], we obtain afte simple calculations, that V ˆˆ = k z, z kp JJ e, e kse m +. Since V 0, function V is bounded. heefoe, z,x, and Y ae bounded vectos. his implies that ˆX and Y ae bounded, too. Consequently, it follows fom (3), that is also bounded. Moeove,s and s ae bounded, too. As can be seen, V is negative fo all ( z,e,e m + ) ( z,e,e m + ) = 0 (Kstic et al., 995) that ( z,e,e m + ) z( ) 0,e( ) 0, and e 0,as. and is zeo only when 0, which implies (using LaSalle-Yoshizawa invaiant theoem tends asymptotically to zeo, i.e. m+, as. By diffeentiating e m + in () 3 de with espect to time, it is also easy to see, that m + is bounded function by 3 dt assumptions egading Θ and γ. his means, that e m + is unifomly continuous. Hence, em+ ( ) 0,, as,, too. he convegence of path velocity and acceleation yields the following equation γe dγ ( ) + e + ( ) = 0. ds m+ m (5) dγ m 0 + m+ = 0. On account of Assumption ds, the second equality is not fulfilled. hus, e m + ( ) = 0 (o equivalently s = ). On account of (3), q( ) 0, as, too. Consequently, Consequently, e + ( ) = o γ e ( ) ( ) bounday conditions (4)-(6) ae (asymptotically) fulfilled and limit point q,e,e,e = 0000,,, is asymptotically stable. Finally, it ( ( ) m+ ( ) ( ) m+ ( )) ( ) should be emphasized, that the chosen Lyapunov function does not guaantee convegence of paamete estimations ˆX and Ŷ to thei tue values.
12 Contol of Redundant Robotic Manipulatos with State Constaints 509 X t 0 Y k 0 pγ = t= On account of (3)-(6), we have Vt = 0 = + + w w Fo sufficiently lage w d and w k, the fist two tems in this equality may be kpγ0 omitted. Hence, we obtain V t = 0 SinceV kpγ0 is not positive, function V fulfils the inequality V. Consequently, the following bound on e may easily be obtained, based on (3) and the last dependence e γ 0 (6) d k 0. An impotant emak may be deived fom the poof caied out. Namely, Inequality (6) pesents an uppe bound (path independent) on the accuacy of path following by the end-effecto accoding to the contol law (8). Let us note that estimation of the uppe bound on path following eo (6) is vey consevative. Consequently, contol gains w and w do not equie lage values to achieve a good path following accuacy, as the numeical simulations (given in the next section) show. Moeove, seveal obsevations can be made egading the contol stategy (8). Fist note, that the poposed contol law equies, in fact no infomation concening the obot kinematic and dynamic equations. Second, the choice of contolle paametes k,k p,k s, w d and w k accoding to dependencies (8)-() guaantees asymptotic stability of the closed-loop eo dynamics () duing the manipulato movement. Moeove, the tanspose of Ĵ (instead of a pseudoinvese) in contol scheme (8) does not esult in numeical instabilities due to (possible) kinematic singulaities met on the obot tajectoy. Nevetheless, (8) has been deived unde the assumption of full-ank adaptive Jacobian matix along the path. Futhemoe, contolle (8) does not equie the knowledge of task space velocity. Due to consevative estimation of the path following accuacy, contol algoithm (8) esults in a bette accuacy of the path following as compaed to uppe bound given fom (6), as the numeical computations caied out in the next section show. In ode to pevent contol (toque) oscillations at the vey beginning of time histoies (caused by e.g. the non-zeo initial path following eo) k s fom (9) should be a bounded, quickly deceasing time dependent function as t (see the next section). d k
13 50 Industial Robotics: heoy, Modelling and Contol Due to eal-time natue of obot contolle (8), we shall ty to estimate the numbe of aithmetic opeations equied to implement the algoithm pesented in this section. he dimension of the obot task space is assumed in estimation to be constant. Opeations equied fo computation of sin,cos, and Θ() functions ae not taken into account. Futhemoe, matices Ĵ, K and D ae assumed in estimation to be given. Moeove, estimations ae caied out at any time instant of the obot task accomplishment. It follows fom (3) and (8) that tems kz,kpˆ J e, equie On ( ) opeations. Computation of the ight hand sides (9) and (0) involves O () and On ( ) opeations, espectively assuming that k = O( n). Computational complexity fo the ight hand side of () equals O(n ) by assumption that d = O( n). Computation of estimated contol ˆ D(q,q,q,q )X equies also the same ode of complexity, i.e. O(n ) opeations. Finally, computational complexity of the whole obot contolle (8) is of the ode O(n ). 4. A numeical example he aim of this section is to illustate the pefomance of the poposed adaptive contol algoithm using a dynamic thee-joint diect-dive am (n = 3) of SCARA-type obotic manipulato opeating in a two-dimensional (m = ) task space. Kinematic scheme of this manipulato and the task to be accomplished is shown in Fig.. In the simulations, SI units ae used. he components of dynamic equations of this manipulato ae as follows (Spong & Vidyasaga, 989): M = [ M ij ] i, j 3 = + + ci=cos( q ), ( ) whee M X X4c3 X6c3 i i s ij = sin( q i + q j ) c ijk = cos( q ) i + qj + qk M = X + X 4c + X 5c3 M = X + X 6 c3 M 33 = X 3 M = M M 3 = X 3 M 3 = M 3 M 3 = M 3 si=sin q, cij=cos(q i+q j ) s ijk = sin qi + qj + qk M 3 = X 3 + X 6c3 ( ).
14 Contol of Redundant Robotic Manipulatos with State Constaints 5 p [m] l q l q q 3 l 3 the geometic path p [m] Figue. A kinematic scheme of the manipulato and the task to be accomplished C = [ C ij ] i, j 3 whee C = ( X 4s + X 5s ) ( X 5s3 + X 6s3 ) 3 C = ( X 4 s + X 5 s)( + ) ( X 5 s + X 6 s) 3 C3 = ( X 5s + X 6s)( + + 3) C = ( X 5s + X 5s) + X 6s3 3 C = ( X 5 s + X 6 s) 3 C3 = X 6s3( ) C3 = ( X 4s + X 5s) X 6s3 C = X s ( + ) q C 33 = 0. = G, ( G G ) G, 3 whee G = X 7c + X 8c + X 9c3 Paametes X, i : 9 i G = X c + X 9c3 8 = take the following nominal values: G = 3 X 9c 3.
15 5 Industial Robotics: heoy, Modelling and Contol X =.956, X = , X 3 = 0.05, X 4 = 0.475, X 5 = 0.80, X 6 = 0.5, X 7 = ( ml c + ml + m3l )g = ( m l m l )g X 8 c + 3 X 9 = m3lc3g,,, whee g stands fo the gavity acceleation m,l, i i and l ci, i = : 3denote link mass, length and location of the mass cente which is assumed to be equal to l ci = l i / l = 0.4 l = 0.36 l 3 = 0.3 m = 3.6 m =.6 m 3 =. Jacobian matix J(q,Y) equals J Ys Ys Y3s3 Ys Y3s3 Y = Yc + Yc + Y3c3 Yc + Y3c3 Y3c3 3 s3 Whee Y = l i = : 3 i i and the kinematic egesso matix takes the following fom K sq = cq c s( + ) s3( ) ( + ) c3( + + ) 3 he end-effecto position p = (pp ) (see Fig. ) epesents in the simulations the task space coodinates (m = ). he uppe bound on the accuacy of the path following in all the compute simulations, is assumed to be equal to 7 γ 0 = 006., whee γ(s) = 000, + 000, e, s. Let us intoduce path following eos
16 Contol of Redundant Robotic Manipulatos with State Constaints 53 e e e 3 p Θ = p Θ s to evaluate the pefomance of the obot contolle (8). In ode to examine the effects of both kinematic and dynamic uncetainties, initial values fo vectos ˆX and Ŷ wee set in the simulations as ˆX( 0) = ( ), Ŷ( 0) = ( ). he task of the obot is to tansfe the end-effecto a- long the geometic path (the dotted line in Fig. ), expessed by the following equations Θ ( s) = s Θ ( s) = s. he initial configuation q 0 equals q 0 = ( π / π / π / ) 0t Paametes k 0,.k 3000,.λ k = 6, e + 30 have been whee, s [ 0, ] = p = = and s ( ) chosen expeimentally to achieve pactically easonable time hoizon of task w = diag 000 and pefomance and elatively small contols with k ( ) w diag ( ) d =. he esults of compute simulation ae pesented in Figs x 0 3 e [m] t [sec] Figue. Path following eo e vs. time
17 54 Industial Robotics: heoy, Modelling and Contol 5 x e [m] t [sec] Figue 3. Path following eo e vs. time e t [sec] Figue 4. Path following eo e3 vs. time 30 5 u [Nm] t [sec] Figue 5. Input toque u vs. time
18 Contol of Redundant Robotic Manipulatos with State Constaints u [Nm] t [sec] Figue 6. Fig. 6 Input toque u vs. time 4 3 u 3 [Nm] t [sec] Figue 7. Input toque u3 vs. time Figue 8. ime couse of adaptive estimate Ŷ
19 56 Industial Robotics: heoy, Modelling and Contol Figue 9. ime couse of adaptive estimate Ŷ Figue 0. ime couse of adaptive estimate Ŷ 3 Figue.. ime couse of adaptive estimate ˆX
20 Contol of Redundant Robotic Manipulatos with State Constaints 57 Figue. ime couse of adaptive estimate ˆX Figue 3. ime couse of adaptive estimate ˆX 3 Figue 4. ime couse of adaptive estimate ˆX 4
21 58 Industial Robotics: heoy, Modelling and Contol Figue 5. ime couse of adaptive estimate ˆX 5 Figue 6. ime couse of adaptive estimate ˆX 6 Figue 7. ime couse of adaptive estimate ˆX 7
22 Contol of Redundant Robotic Manipulatos with State Constaints 59 Figue 8. ime couse of adaptive estimate ˆX 8 Figue 9. ime couse of adaptive estimate ˆX p [m] p [m] Figue 0. Manipulato motion along the geometic path
23 50 Industial Robotics: heoy, Modelling and Contol As might be expected, the path following eos fom Figs -3 ae much smalle than those obtained fom the consevative dependence (6). Moeove, as one can obseve fom Figs -7, the time dependent damping function k s deceases (eliminates) eos and toques oscillations at the vey beginning of time histoies. Futhemoe, as seen fom Figs 8-9, estimations Xˆ, Yˆ do not convege to thei eal (nominal) values. 5. Conclusion his study has pesented an adaptive obot contolle fo the path following by the end-effecto. he contol geneation scheme has been deived using the Lyapunov stability theoy. An advantage of the poposed contol law (8) is that it equies, in fact no infomation egading the paametes of the obot dynamic equations. he contol stategy (8) is shown to be asymptotically stable (by fulfilment of pactically easonable assumptions). he poposed obot contolle has been applied to a plana edundant manipulato of thee evolute kinematic pais opeating in a two dimensional task space. Numeical simulations have shown that the esults obtained ae in accodance with the theoetical analysis. he novelty of the stategy poposed lies in its elative simplicity in design, pogam code and eal-time implementation. he appoach pesented hee will also be in futue diectly applicable to coopeating kinematically edundant manipulatos. Acknowledgement. his wok was suppoted by the DFG Ga 65/--,. 6. Refeences Aimoto, S. (990). Design of obot contol systems, Adv. Robot., vol. 4, no., pp Aimoto, S. (996). Contol theoy of non-linea mechanical systems, Oxfod, U. K., Claendon. Beghuis, H. R. Otega, & H. Nijmeije (993). A obust adaptive obot contolle, IEEE ans. on Robotics and Automat., vol. 9, no. 6, pp Canudas, de Wit, C. B. Siciliano & G. Bastin (996). heoy of obot contol, New Yok: Spinge-Velag. Cheah, C. C. S. Kawamoa & S. Aimoto (999). Feedback contol fo obotic manipulatos with an uncetain Jacobian matix, J. Robot. Syst., vol. 6, no., pp
24 Contol of Redundant Robotic Manipulatos with State Constaints 5 Cheach, C. C. M. Hiano S. Kawamua & S. Aimoto (003). Appoximate jacobian contol fo obots with uncetain kinematics and dynamics, IEEE ans. on Robotics and Automation, vol. 9, no. 4, pp Dahl, O. (994). Path-constained obot contol with limited toques - Expeimental evaluation, IEEE J. Robot. Automat., vol. 0, pp Feng, G. & M. Palaniswami (99). Adaptive contol of obot manipulatos in task space, IEEE ans. on Automatic Contol, vol. 38, no., pp Galicki, M. (998a). he planning of obotic optimal motions in the pesence of obstacles, Int. J. Robotics Res., vol. 7, no. 3, pp Galicki M. (998). he stuctue of time optimal contols fo kinematically edundant manipulatos with end-effecto path constaints, Poc. IEEE Conf. Robotics Automat., pp.0-06, Leuven, Belgium. Galicki, M. & I. Pajak (999). Optimal motion of edundant manipulatos with state equality constaints, Poc. IEEE Int. Symp. on Assembly and ask Planning, pp. 8-85, Poto, Potugal. Galicki, M. (000). ime-optimal contols of kinematically edundant manipulatos with geometic constaints, IEEE ans. Robotics Automat., vol. 6, no., pp Galicki, M. (00). Real-time tajectoy geneation fo edundant manipulatos with path constaints, Int. J. Robotics Res., vol. 0, no. 7, pp Galicki, M. (00). Motion contol of obotic manipulatos in task space, Poc. IEEE/RSJ Int. Conf. Intell. Robots and Systems, pp , Lausanne. Galicki, M. (004). Path following by the end-effecto of a edundant manipulato opeating in a dynamic envionment, IEEE ans. Robotics, vol. 0, no 6, pp Galicki, M. (005). Collision-fee contol of obotic manipulatos in the task space, J. Robot., Syst., vol., no 8, pp Galicki, M. (006a). Path-constained contol of a edundant manipulato in a task space, Robotics and Autonomous Systems, vol. 54, pp Galicki, M. (006b). Adaptive path-constained contol of a obotic manipulato in a task space, Robotica (to appea). Kelly R. (999). Regulation of manipulatos in geneic task space: An enegy shaping plus damping injection appoach, IEEE ans. Rob. Automat., vol. 5, no., pp Kiefe, J. A. J. Cahill & M. R. (997). James, Robust and accuate time-optimal path-tacking contol fo obot manipulatos, IEEE ans. Automat. Robot., vol. 3, pp Kstic,M. I. Kanellakopoulos & P. Kokotovic (995). Nonlinea and adaptive contol design, J. Wiley and Sons, New Yok. Lewis, F. L. C.. Abdallach & D. M. Dawson (993). Contol of obot manipulatos, New Yok: Macmillan Publishing Company.
25 5 Industial Robotics: heoy, Modelling and Contol Sciavicco, L. & B. Siciliano (996). Modeling and contol of obot manipulatos, New Yok: McGaw-Hill. Slotine, J. J. E. & W. Li (987). On adaptive contol of obot manipulatos, Int. J. Robotics Res., no. 6, pp Slotine, J. J. E. & W. Li (99). Applied nonlinea contol, Englewood Cliffs, New Jesey: Pentice Hall. Spong, M. & M. Vidyasaga (989). Robot dynamics and contol, Wiley, New Yok. akegaki, M. & S. Aimoto (98). A new feedback method fo dynamic contol of manipulatos, ans. ASME: J. Dyn. Syst., Meas., Cont., vol. 0, pp omei, P. (000). Robust adaptive fiction compensation fo tacking contol of obot manipulatos, IEEE tans. on Automatic Contol, vol. 45, no., pp , 000. Yazael, H. & C. C. Cheach (00). ask-space adaptive contol of obotic manipulatos with uncetainties in gavity egesso matix and kinematics, IEEE ans. on Automatic Contol, vol. 47, no. 9, pp
26 Industial Robotics: heoy, Modelling and Contol Edited by Sam Cubeo ISBN Had cove, 964 pages Publishe Po Liteatu Velag, Gemany / ARS, Austia Published online 0, Decembe, 006 Published in pint edition Decembe, 006 his book coves a wide ange of topics elating to advanced industial obotics, sensos and automation technologies. Although being highly technical and complex in natue, the papes pesented in this book epesent some of the latest cutting edge technologies and advancements in industial obotics technology. his book coves topics such as netwoking, popeties of manipulatos, fowad and invese obot am kinematics, motion path-planning, machine vision and many othe pactical topics too numeous to list hee. he authos and edito of this book wish to inspie people, especially young ones, to get involved with obotic and mechatonic engineeing technology and to develop new and exciting pactical applications, pehaps using the ideas and concepts pesented heein. How to efeence In ode to coectly efeence this scholaly wok, feel fee to copy and paste the following: Mioslaw Galicki (006). Contol of Redundant Robotic Manipulatos with State Constaints, Industial Robotics: heoy, Modelling and Contol, Sam Cubeo (Ed.), ISBN: , Inech, Available fom: botic_manipulatos_with_state_constaints Inech Euope Univesity Campus SeP Ri Slavka Kautzeka 83/A 5000 Rijeka, Coatia Phone: +385 (5) Fax: +385 (5) Inech China Unit 405, Office Block, Hotel Equatoial Shanghai No.65, Yan An Road (West), Shanghai, 00040, China Phone: Fax:
27 006 he Autho(s). Licensee IntechOpen. his chapte is distibuted unde the tems of the Ceative Commons Attibution-NonCommecial-ShaeAlike-3.0 License, which pemits use, distibution and epoduction fo non-commecial puposes, povided the oiginal is popely cited and deivative woks building on this content ae distibuted unde the same license.
A NEW VARIABLE STIFFNESS SPRING USING A PRESTRESSED MECHANISM
Poceedings of the ASME 2010 Intenational Design Engineeing Technical Confeences & Computes and Infomation in Engineeing Confeence IDETC/CIE 2010 August 15-18, 2010, Monteal, Quebec, Canada DETC2010-28496
More informationC e f paamete adaptation f (' x) ' ' d _ d ; ; e _e K p K v u ^M() RBF NN ^h( ) _ obot s _ s n W ' f x x xm xm f x xm d Figue : Block diagam of comput
A Neual-Netwok Compensato with Fuzzy Robustication Tems fo Impoved Design of Adaptive Contol of Robot Manipulatos Y.H. FUNG and S.K. TSO Cente fo Intelligent Design, Automation and Manufactuing City Univesity
More informationImplicit Constraint Enforcement for Rigid Body Dynamic Simulation
Implicit Constaint Enfocement fo Rigid Body Dynamic Simulation Min Hong 1, Samuel Welch, John app, and Min-Hyung Choi 3 1 Division of Compute Science and Engineeing, Soonchunhyang Univesity, 646 Eupnae-i
More informationPearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms
Peason s Chi-Squae Test Modifications fo Compaison of Unweighted and Weighted Histogams and Two Weighted Histogams Univesity of Akueyi, Bogi, v/noduslód, IS-6 Akueyi, Iceland E-mail: nikolai@unak.is Two
More informationResearch Article On Alzer and Qiu s Conjecture for Complete Elliptic Integral and Inverse Hyperbolic Tangent Function
Abstact and Applied Analysis Volume 011, Aticle ID 697547, 7 pages doi:10.1155/011/697547 Reseach Aticle On Alze and Qiu s Conjectue fo Complete Elliptic Integal and Invese Hypebolic Tangent Function Yu-Ming
More informationGradient-based Neural Network for Online Solution of Lyapunov Matrix Equation with Li Activation Function
Intenational Confeence on Infomation echnology and Management Innovation (ICIMI 05) Gadient-based Neual Netwok fo Online Solution of Lyapunov Matix Equation with Li Activation unction Shiheng Wang, Shidong
More informationAalborg Universitet. Load Estimation from Natural input Modal Analysis Aenlle, Manuel López; Brincker, Rune; Canteli, Alfonso Fernández
Aalbog Univesitet Load Estimation fom atual input Modal Analysis Aenlle, Manuel López; Bincke, Rune; Canteli, Alfonso Fenández Published in: Confeence Poceedings Publication date: 005 Document Vesion Publishe's
More informationCASCADE OPTIMIZATION AND CONTROL OF BATCH REACTORS
CASCADE OPIMIZAION AND CONROL OF BACH REACORS Xiangming Hua, Sohab Rohani and Athu Jutan* Depatment of Chemical and Biochemical Engineeing Univesity of Westen Ontaio, London, Canada N6A B9 * ajutan@uwo.ca
More informationTerminal Sliding Control for a SCARA Robot
Intenational Confeence on Contol, Instumentation and Mechatonics Engineeing (CIM '7), Joho Bahu, Joho, Malaysia, May 8-9, 7 Teminal Sliding Contol fo a SCARA Robot E. Muñoz, C. Gaviia, A. Vivas Depatment
More informationtime [s] time [s]
ROBUST ATTITUDE STABILIZATION OF AN UNDERACTUATED AUV K. Y. Pettesen and O. Egeland Depatment of Engineeing Cybenetics Nowegian Univesity of Science and Technology N- Tondheim, Noway Fax: + 9 99 E-mail:
More informationLocalization of Eigenvalues in Small Specified Regions of Complex Plane by State Feedback Matrix
Jounal of Sciences, Islamic Republic of Ian (): - () Univesity of Tehan, ISSN - http://sciencesutaci Localization of Eigenvalues in Small Specified Regions of Complex Plane by State Feedback Matix H Ahsani
More informationEM Boundary Value Problems
EM Bounday Value Poblems 10/ 9 11/ By Ilekta chistidi & Lee, Seung-Hyun A. Geneal Desciption : Maxwell Equations & Loentz Foce We want to find the equations of motion of chaged paticles. The way to do
More informationJ. Electrical Systems 1-3 (2005): Regular paper
K. Saii D. Rahem S. Saii A Miaoui Regula pape Coupled Analytical-Finite Element Methods fo Linea Electomagnetic Actuato Analysis JES Jounal of Electical Systems In this pape, a linea electomagnetic actuato
More informationCOUPLED MODELS OF ROLLING, SLIDING AND WHIRLING FRICTION
ENOC 008 Saint Petesbug Russia June 30-July 4 008 COUPLED MODELS OF ROLLING SLIDING AND WHIRLING FRICTION Alexey Kieenkov Ins ti tu te fo P ob le ms in Me ch an ic s Ru ss ia n Ac ad em y of Sc ie nc es
More informationModel and Controller Order Reduction for Infinite Dimensional Systems
IT J. Eng. Sci., Vol. 4, No.,, -6 Model and Contolle Ode Reduction fo Infinite Dimensional Systems Fatmawati,*, R. Saagih,. Riyanto 3 & Y. Soehayadi Industial and Financial Mathematics Goup email: fatma47@students.itb.ac.id;
More informationRigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018
Rigid Body Dynamics 2 CSE169: Compute Animation nstucto: Steve Rotenbeg UCSD, Winte 2018 Coss Poduct & Hat Opeato Deivative of a Rotating Vecto Let s say that vecto is otating aound the oigin, maintaining
More informationOn the global uniform asymptotic stability of time-varying dynamical systems
Stud. Univ. Babeş-Bolyai Math. 59014), No. 1, 57 67 On the global unifom asymptotic stability of time-vaying dynamical systems Zaineb HajSalem, Mohamed Ali Hammami and Mohamed Mabouk Abstact. The objective
More informationQuaternion Based Inverse Kinematics for Industrial Robot Manipulators with Euler Wrist
Quatenion Based Invese Kinematics fo Industial Robot Manipulatos with Eule Wist Yavuz Aydın Electonics and Compute Education Kocaeli Univesity Umuttepe Kocaeli Tukey yavuz_98@hotmailcom Seda Kucuk Electonics
More informationState tracking control for Takagi-Sugeno models
State tacing contol fo Taagi-Sugeno models Souad Bezzaoucha, Benoît Max,3,DidieMaquin,3 and José Ragot,3 Abstact This wo addesses the model efeence tacing contol poblem It aims to highlight the encouteed
More informationPHYS 110B - HW #7 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased
PHYS 0B - HW #7 Sping 2004, Solutions by David Pace Any efeenced euations ae fom Giffiths Poblem statements ae paaphased. Poblem 0.3 fom Giffiths A point chage,, moves in a loop of adius a. At time t 0
More informationPerturbation to Symmetries and Adiabatic Invariants of Nonholonomic Dynamical System of Relative Motion
Commun. Theo. Phys. Beijing, China) 43 25) pp. 577 581 c Intenational Academic Publishes Vol. 43, No. 4, Apil 15, 25 Petubation to Symmeties and Adiabatic Invaiants of Nonholonomic Dynamical System of
More informationRotor Flux Estimation of Induction Motors Using Sliding-Mode Observer
5th Intenational Confeence on Sustainable Enegy and Envionment Engineeing (ICSEEE 2016) Roto Flux Estimation of Induction Motos Using Sliding-Mode Obseve Yong Feng1,a, Minghao Zhou1,b and Fengling Han2,c
More informationChapter 5 Linear Equations: Basic Theory and Practice
Chapte 5 inea Equations: Basic Theoy and actice In this chapte and the next, we ae inteested in the linea algebaic equation AX = b, (5-1) whee A is an m n matix, X is an n 1 vecto to be solved fo, and
More informationAnalysis of high speed machining center spindle dynamic unit structure performance Yuan guowei
Intenational Confeence on Intelligent Systems Reseach and Mechatonics Engineeing (ISRME 0) Analysis of high speed machining cente spindle dynamic unit stuctue pefomance Yuan guowei Liaoning jidian polytechnic,dan
More informationA scaling-up methodology for co-rotating twin-screw extruders
A scaling-up methodology fo co-otating twin-scew extudes A. Gaspa-Cunha, J. A. Covas Institute fo Polymes and Composites/I3N, Univesity of Minho, Guimaães 4800-058, Potugal Abstact. Scaling-up of co-otating
More information3.1 Random variables
3 Chapte III Random Vaiables 3 Random vaiables A sample space S may be difficult to descibe if the elements of S ae not numbes discuss how we can use a ule by which an element s of S may be associated
More informationASTR415: Problem Set #6
ASTR45: Poblem Set #6 Cuan D. Muhlbege Univesity of Mayland (Dated: May 7, 27) Using existing implementations of the leapfog and Runge-Kutta methods fo solving coupled odinay diffeential equations, seveal
More informationON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION. 1. Introduction. 1 r r. r k for every set E A, E \ {0},
ON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION E. J. IONASCU and A. A. STANCU Abstact. We ae inteested in constucting concete independent events in puely atomic pobability
More information4/18/2005. Statistical Learning Theory
Statistical Leaning Theoy Statistical Leaning Theoy A model of supevised leaning consists of: a Envionment - Supplying a vecto x with a fixed but unknown pdf F x (x b Teache. It povides a desied esponse
More informationChapter 5 Force and Motion
Chapte 5 Foce and Motion In chaptes 2 and 4 we have studied kinematics i.e. descibed the motion of objects using paametes such as the position vecto, velocity and acceleation without any insights as to
More informationMagnetometer Calibration Algorithm Based on Analytic Geometry Transform Yongjian Yang, Xiaolong Xiao1,Wu Liao
nd Intenational Foum on Electical Engineeing and Automation (IFEEA 5 Magnetomete Calibation Algoithm Based on Analytic Geomety ansfom Yongjian Yang, Xiaolong Xiao,u Liao College of Compute Science and
More informationCOMPUTATIONS OF ELECTROMAGNETIC FIELDS RADIATED FROM COMPLEX LIGHTNING CHANNELS
Pogess In Electomagnetics Reseach, PIER 73, 93 105, 2007 COMPUTATIONS OF ELECTROMAGNETIC FIELDS RADIATED FROM COMPLEX LIGHTNING CHANNELS T.-X. Song, Y.-H. Liu, and J.-M. Xiong School of Mechanical Engineeing
More informationChapter 5 Force and Motion
Chapte 5 Foce and Motion In Chaptes 2 and 4 we have studied kinematics, i.e., we descibed the motion of objects using paametes such as the position vecto, velocity, and acceleation without any insights
More informationDesigning a Sine-Coil for Measurement of Plasma Displacements in IR-T1 Tokamak
Designing a Sine-Coil fo Measuement of Plasma Displacements in IR-T Tokamak Pejman Khoshid, M. Razavi, M. Ghoanneviss, M. Molaii, A. TalebiTahe, R. Avin, S. Mohammadi and A. NikMohammadi Dept. of Physics,
More informationAn Exact Solution of Navier Stokes Equation
An Exact Solution of Navie Stokes Equation A. Salih Depatment of Aeospace Engineeing Indian Institute of Space Science and Technology, Thiuvananthapuam, Keala, India. July 20 The pincipal difficulty in
More informationSTABILITY AND PARAMETER SENSITIVITY ANALYSES OF AN INDUCTION MOTOR
HUNGARIAN JOURNAL OF INDUSTRY AND CHEMISTRY VESZPRÉM Vol. 42(2) pp. 109 113 (2014) STABILITY AND PARAMETER SENSITIVITY ANALYSES OF AN INDUCTION MOTOR ATTILA FODOR 1, ROLAND BÁLINT 1, ATTILA MAGYAR 1, AND
More informationTHE INFLUENCE OF THE MAGNETIC NON-LINEARITY ON THE MAGNETOSTATIC SHIELDS DESIGN
THE INFLUENCE OF THE MAGNETIC NON-LINEARITY ON THE MAGNETOSTATIC SHIELDS DESIGN LIVIU NEAMŢ 1, ALINA NEAMŢ, MIRCEA HORGOŞ 1 Key wods: Magnetostatic shields, Magnetic non-lineaity, Finite element method.
More information7.2. Coulomb s Law. The Electric Force
Coulomb s aw Recall that chaged objects attact some objects and epel othes at a distance, without making any contact with those objects Electic foce,, o the foce acting between two chaged objects, is somewhat
More informationDuality between Statical and Kinematical Engineering Systems
Pape 00, Civil-Comp Ltd., Stiling, Scotland Poceedings of the Sixth Intenational Confeence on Computational Stuctues Technology, B.H.V. Topping and Z. Bittna (Editos), Civil-Comp Pess, Stiling, Scotland.
More informationAnalytical Solutions for Confined Aquifers with non constant Pumping using Computer Algebra
Poceedings of the 006 IASME/SEAS Int. Conf. on ate Resouces, Hydaulics & Hydology, Chalkida, Geece, May -3, 006 (pp7-) Analytical Solutions fo Confined Aquifes with non constant Pumping using Compute Algeba
More informationF-IF Logistic Growth Model, Abstract Version
F-IF Logistic Gowth Model, Abstact Vesion Alignments to Content Standads: F-IFB4 Task An impotant example of a model often used in biology o ecology to model population gowth is called the logistic gowth
More informationHua Xu 3 and Hiroaki Mukaidani 33. The University of Tsukuba, Otsuka. Hiroshima City University, 3-4-1, Ozuka-Higashi
he inea Quadatic Dynamic Game fo Discete-ime Descipto Systems Hua Xu 3 and Hioai Muaidani 33 3 Gaduate School of Systems Management he Univesity of suuba, 3-9- Otsua Bunyo-u, oyo -0, Japan xuhua@gssm.otsua.tsuuba.ac.jp
More information1 Spherical multipole moments
Jackson notes 9 Spheical multipole moments Suppose we have a chage distibution ρ (x) wheeallofthechageiscontained within a spheical egion of adius R, as shown in the diagam. Then thee is no chage in the
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II September 15, 2012 Prof. Alan Guth PROBLEM SET 2
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Depatment Physics 8.07: Electomagnetism II Septembe 5, 202 Pof. Alan Guth PROBLEM SET 2 DUE DATE: Monday, Septembe 24, 202. Eithe hand it in at the lectue,
More informationDetermining solar characteristics using planetary data
Detemining sola chaacteistics using planetay data Intoduction The Sun is a G-type main sequence sta at the cente of the Sola System aound which the planets, including ou Eath, obit. In this investigation
More informationAbsorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere
Applied Mathematics, 06, 7, 709-70 Published Online Apil 06 in SciRes. http://www.scip.og/jounal/am http://dx.doi.og/0.46/am.06.77065 Absoption Rate into a Small Sphee fo a Diffusing Paticle Confined in
More informationConservative Averaging Method and its Application for One Heat Conduction Problem
Poceedings of the 4th WSEAS Int. Conf. on HEAT TRANSFER THERMAL ENGINEERING and ENVIRONMENT Elounda Geece August - 6 (pp6-) Consevative Aveaging Method and its Application fo One Heat Conduction Poblem
More informationRelating Branching Program Size and. Formula Size over the Full Binary Basis. FB Informatik, LS II, Univ. Dortmund, Dortmund, Germany
Relating Banching Pogam Size and omula Size ove the ull Binay Basis Matin Saueho y Ingo Wegene y Ralph Wechne z y B Infomatik, LS II, Univ. Dotmund, 44 Dotmund, Gemany z ankfut, Gemany sauehof/wegene@ls.cs.uni-dotmund.de
More informationON THE TWO-BODY PROBLEM IN QUANTUM MECHANICS
ON THE TWO-BODY PROBLEM IN QUANTUM MECHANICS L. MICU Hoia Hulubei National Institute fo Physics and Nuclea Engineeing, P.O. Box MG-6, RO-0775 Buchaest-Maguele, Romania, E-mail: lmicu@theoy.nipne.o (Received
More informationSTATE VARIANCE CONSTRAINED FUZZY CONTROL VIA OBSERVER-BASED FUZZY CONTROLLERS
Jounal of Maine Science and echnology, Vol. 4, No., pp. 49-57 (6) 49 SAE VARIANCE CONSRAINED FUZZY CONROL VIA OBSERVER-BASED FUZZY CONROLLERS Wen-Je Chang*, Yi-Lin Yeh**, and Yu-eh Meng*** Key wods: takagi-sugeno
More information( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.
9. LAGRANGIAN OF THE ELECTROMAGNETIC FIELD In the pevious section the Lagangian and Hamiltonian of an ensemble of point paticles was developed. This appoach is based on a qt. This discete fomulation can
More informationEFFECTS OF FRINGING FIELDS ON SINGLE PARTICLE DYNAMICS. M. Bassetti and C. Biscari INFN-LNF, CP 13, Frascati (RM), Italy
Fascati Physics Seies Vol. X (998), pp. 47-54 4 th Advanced ICFA Beam Dynamics Wokshop, Fascati, Oct. -5, 997 EFFECTS OF FRININ FIELDS ON SINLE PARTICLE DYNAMICS M. Bassetti and C. Biscai INFN-LNF, CP
More informationChapter Introduction to Finite Element Methods
Chapte 1.4 Intoduction to Finite Element Methods Afte eading this chapte, you should e ale to: 1. Undestand the asics of finite element methods using a one-dimensional polem. In the last fifty yeas, the
More information12th WSEAS Int. Conf. on APPLIED MATHEMATICS, Cairo, Egypt, December 29-31,
th WSEAS Int. Conf. on APPLIED MATHEMATICS, Caio, Egypt, Decembe 9-3, 7 5 Magnetostatic Field calculations associated with thick Solenoids in the Pesence of Ion using a Powe Seies expansion and the Complete
More informationPhysics 235 Chapter 5. Chapter 5 Gravitation
Chapte 5 Gavitation In this Chapte we will eview the popeties of the gavitational foce. The gavitational foce has been discussed in geat detail in you intoductoy physics couses, and we will pimaily focus
More informationELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS
THE 9 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS R. Sbulati *, S. R. Atashipou Depatment of Civil, Chemical and Envionmental Engineeing,
More informationAPPLICATION OF MAC IN THE FREQUENCY DOMAIN
PPLICION OF MC IN HE FREQUENCY DOMIN D. Fotsch and D. J. Ewins Dynamics Section, Mechanical Engineeing Depatment Impeial College of Science, echnology and Medicine London SW7 2B, United Kingdom BSRC he
More informationAppendix B The Relativistic Transformation of Forces
Appendix B The Relativistic Tansfomation of oces B. The ou-foce We intoduced the idea of foces in Chapte 3 whee we saw that the change in the fou-momentum pe unit time is given by the expession d d w x
More informationLight Time Delay and Apparent Position
Light Time Delay and ppaent Position nalytical Gaphics, Inc. www.agi.com info@agi.com 610.981.8000 800.220.4785 Contents Intoduction... 3 Computing Light Time Delay... 3 Tansmission fom to... 4 Reception
More informationMODULE 5a and 5b (Stewart, Sections 12.2, 12.3) INTRO: In MATH 1114 vectors were written either as rows (a1, a2,..., an) or as columns a 1 a. ...
MODULE 5a and 5b (Stewat, Sections 2.2, 2.3) INTRO: In MATH 4 vectos wee witten eithe as ows (a, a2,..., an) o as columns a a 2... a n and the set of all such vectos of fixed length n was called the vecto
More informationMultiple Criteria Secretary Problem: A New Approach
J. Stat. Appl. Po. 3, o., 9-38 (04 9 Jounal of Statistics Applications & Pobability An Intenational Jounal http://dx.doi.og/0.785/jsap/0303 Multiple Citeia Secetay Poblem: A ew Appoach Alaka Padhye, and
More informationPROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.
POBLM S # SOLUIONS by obet A. DiStasio J. Q. he Bon-Oppenheime appoximation is the standad way of appoximating the gound state of a molecula system. Wite down the conditions that detemine the tonic and
More informationRight-handed screw dislocation in an isotropic solid
Dislocation Mechanics Elastic Popeties of Isolated Dislocations Ou study of dislocations to this point has focused on thei geomety and thei ole in accommodating plastic defomation though thei motion. We
More informationEncapsulation theory: radial encapsulation. Edmund Kirwan *
Encapsulation theoy: adial encapsulation. Edmund Kiwan * www.edmundkiwan.com Abstact This pape intoduces the concept of adial encapsulation, wheeby dependencies ae constained to act fom subsets towads
More informationHow to Obtain Desirable Transfer Functions in MIMO Systems Under Internal Stability Using Open and Closed Loop Control
How to Obtain Desiable ansfe Functions in MIMO Sstems Unde Intenal Stabilit Using Open and losed Loop ontol echnical Repot of the ISIS Goup at the Univesit of Note Dame ISIS-03-006 June, 03 Panos J. Antsaklis
More informationAST 121S: The origin and evolution of the Universe. Introduction to Mathematical Handout 1
Please ead this fist... AST S: The oigin and evolution of the Univese Intoduction to Mathematical Handout This is an unusually long hand-out and one which uses in places mathematics that you may not be
More informationAn Adaptive Neural-Network Model-Following Speed Control of PMSM Drives for Electric Vehicle Applications
Poceedings of the 9th WSEAS Intenational Confeence on Applied Mathematics, Istanbul, Tuey, May 27-29, 2006 (pp412-417) An Adaptive Neual-Netwo Model-Following Speed Contol of PMSM Dives fo Electic Vehicle
More informationarxiv: v1 [physics.pop-ph] 3 Jun 2013
A note on the electostatic enegy of two point chages axiv:1306.0401v1 [physics.pop-ph] 3 Jun 013 A C Tot Instituto de Física Univesidade Fedeal do io de Janeio Caixa Postal 68.58; CEP 1941-97 io de Janeio,
More informationKinematic Motion Control of Wheeled Mobile Robots Considering Curvature Constraints
8 IEEE Intenational Confeence on Robotics and Automation Pasadena, CA, USA, May 9-3, 8 Kinematic Motion Contol of Wheeled Mobile Robots Consideing Cuvatue Constaints Youngshik Kim, Student Membe, IEEE
More informationworking pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50
woking pages fo Paul Richads class notes; do not copy o ciculate without pemission fom PGR 2004/11/3 10:50 CHAPTER7 Solid angle, 3D integals, Gauss s Theoem, and a Delta Function We define the solid angle,
More informationPredictive Filtering for Nonlinear Systems
Pedictive Filteing fo Nonlinea Systems John L. Cassidis F. Landis Makley Abstact In this pape, a eal-time pedictive filte is deived fo nonlinea systems. he majo advantage of this new filte ove conventional
More informationSIO 229 Gravity and Geomagnetism. Lecture 6. J 2 for Earth. J 2 in the solar system. A first look at the geoid.
SIO 229 Gavity and Geomagnetism Lectue 6. J 2 fo Eath. J 2 in the sola system. A fist look at the geoid. The Thee Big Themes of the Gavity Lectues 1.) An ellipsoidal otating Eath Refeence body (mass +
More informationIntroduction to Nuclear Forces
Intoduction to Nuclea Foces One of the main poblems of nuclea physics is to find out the natue of nuclea foces. Nuclea foces diffe fom all othe known types of foces. They cannot be of electical oigin since
More informationAnalytical time-optimal trajectories for an omni-directional vehicle
Analytical time-optimal tajectoies fo an omni-diectional vehicle Weifu Wang and Devin J. Balkcom Abstact We pesent the fist analytical solution method fo finding a time-optimal tajectoy between any given
More informationTo Feel a Force Chapter 7 Static equilibrium - torque and friction
To eel a oce Chapte 7 Chapte 7: Static fiction, toque and static equilibium A. Review of foce vectos Between the eath and a small mass, gavitational foces of equal magnitude and opposite diection act on
More informationParticle Systems. University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell
Paticle Systems Univesity of Texas at Austin CS384G - Compute Gaphics Fall 2010 Don Fussell Reading Requied: Witkin, Paticle System Dynamics, SIGGRAPH 97 couse notes on Physically Based Modeling. Witkin
More informationOn a quantity that is analogous to potential and a theorem that relates to it
Su une quantité analogue au potential et su un théoème y elatif C R Acad Sci 7 (87) 34-39 On a quantity that is analogous to potential and a theoem that elates to it By R CLAUSIUS Tanslated by D H Delphenich
More informationSolving Problems of Advance of Mercury s Perihelion and Deflection of. Photon Around the Sun with New Newton s Formula of Gravity
Solving Poblems of Advance of Mecuy s Peihelion and Deflection of Photon Aound the Sun with New Newton s Fomula of Gavity Fu Yuhua (CNOOC Reseach Institute, E-mail:fuyh945@sina.com) Abstact: Accoding to
More informationFractional Tikhonov regularization for linear discrete ill-posed problems
BIT manuscipt No. (will be inseted by the edito) Factional Tikhonov egulaization fo linea discete ill-posed poblems Michiel E. Hochstenbach Lotha Reichel Received: date / Accepted: date Abstact Tikhonov
More informationAn Application of Fuzzy Linear System of Equations in Economic Sciences
Austalian Jounal of Basic and Applied Sciences, 5(7): 7-14, 2011 ISSN 1991-8178 An Application of Fuzzy Linea System of Equations in Economic Sciences 1 S.H. Nassei, 2 M. Abdi and 3 B. Khabii 1 Depatment
More informationA new approach in classical electrodynamics to protect principle of causality
A new appoach in classical electodynamics to potect pinciple of causality Biswaanjan Dikshit * Lase and Plasma Technology Division Bhabha Atomic Reseach Cente, Mumbai-400085 INDIA * Coesponding autho E-mail:
More information1D2G - Numerical solution of the neutron diffusion equation
DG - Numeical solution of the neuton diffusion equation Y. Danon Daft: /6/09 Oveview A simple numeical solution of the neuton diffusion equation in one dimension and two enegy goups was implemented. Both
More informationFresnel Diffraction. monchromatic light source
Fesnel Diffaction Equipment Helium-Neon lase (632.8 nm) on 2 axis tanslation stage, Concave lens (focal length 3.80 cm) mounted on slide holde, iis mounted on slide holde, m optical bench, micoscope slide
More information6 Matrix Concentration Bounds
6 Matix Concentation Bounds Concentation bounds ae inequalities that bound pobabilities of deviations by a andom vaiable fom some value, often its mean. Infomally, they show the pobability that a andom
More informationFUNDAMENTAL PROPERTIES OF LINEAR SHIP STEERING DYNAMIC MODELS
Jounal of Maine Science and Technology, Vol. 7, No. 2, pp. 79-88 (1999) 79 FUNDAMENTAL PROPERTIES OF LINEAR SHIP STEERING DYNAMIC MODELS Ching-Yaw Tzeng* and Ju-Fen Chen** Keywods: Nomoto model, Contollability
More informationHigh precision computer simulation of cyclotrons KARAMYSHEVA T., AMIRKHANOV I. MALININ V., POPOV D.
High pecision compute simulation of cyclotons KARAMYSHEVA T., AMIRKHANOV I. MALININ V., POPOV D. Abstact Effective and accuate compute simulations ae highly impotant in acceleatos design and poduction.
More informationPhysics 107 TUTORIAL ASSIGNMENT #8
Physics 07 TUTORIAL ASSIGNMENT #8 Cutnell & Johnson, 7 th edition Chapte 8: Poblems 5,, 3, 39, 76 Chapte 9: Poblems 9, 0, 4, 5, 6 Chapte 8 5 Inteactive Solution 8.5 povides a model fo solving this type
More informationAdaptive Backstepping Output Feedback Control for SISO Nonlinear System Using Fuzzy Neural Networks
Intenational Jounal of Automation and Computing 6(), May 009, 45-53 DOI: 0.007/s633-009-045-0 Adaptive Backstepping Output Feedback Contol fo SISO Nonlinea System Using Fuzzy Neual Netwoks Shao-Cheng Tong
More informationMoment-free numerical approximation of highly oscillatory integrals with stationary points
Moment-fee numeical appoximation of highly oscillatoy integals with stationay points Sheehan Olve Abstact We pesent a method fo the numeical quadatue of highly oscillatoy integals with stationay points.
More informationGeneral Railgun Function
Geneal ailgun Function An electomagnetic ail gun uses a lage Loentz foce to fie a pojectile. The classic configuation uses two conducting ails with amatue that fits between and closes the cicuit between
More informationSteady State and Transient Performance Analysis of Three Phase Induction Machine using MATLAB Simulations
Intenational Jounal of Recent Tends in Engineeing, Vol, No., May 9 Steady State and Tansient Pefomance Analysis of Thee Phase Induction Machine using MATAB Simulations Pof. Himanshu K. Patel Assistant
More informationIntroduction to Mathematical Statistics Robert V. Hogg Joeseph McKean Allen T. Craig Seventh Edition
Intoduction to Mathematical Statistics Robet V. Hogg Joeseph McKean Allen T. Caig Seventh Edition Peason Education Limited Edinbugh Gate Halow Essex CM2 2JE England and Associated Companies thoughout the
More informationHomework # 3 Solution Key
PHYSICS 631: Geneal Relativity Homewok # 3 Solution Key 1. You e on you hono not to do this one by hand. I ealize you can use a compute o simply look it up. Please don t. In a flat space, the metic in
More informationGeometry of the homogeneous and isotropic spaces
Geomety of the homogeneous and isotopic spaces H. Sonoda Septembe 2000; last evised Octobe 2009 Abstact We summaize the aspects of the geomety of the homogeneous and isotopic spaces which ae most elevant
More informationA Power Method for Computing Square Roots of Complex Matrices
JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 13, 39345 1997 ARTICLE NO. AY975517 A Powe Method fo Computing Squae Roots of Complex Matices Mohammed A. Hasan Depatment of Electical Engineeing, Coloado
More informationA Comparison and Contrast of Some Methods for Sample Quartiles
A Compaison and Contast of Some Methods fo Sample Quatiles Anwa H. Joade and aja M. Latif King Fahd Univesity of Petoleum & Mineals ABSTACT A emainde epesentation of the sample size n = 4m ( =, 1, 2, 3)
More informationSurveillance Points in High Dimensional Spaces
Société de Calcul Mathématique SA Tools fo decision help since 995 Suveillance Points in High Dimensional Spaces by Benad Beauzamy Januay 06 Abstact Let us conside any compute softwae, elying upon a lage
More informationMotion in One Dimension
Motion in One Dimension Intoduction: In this lab, you will investigate the motion of a olling cat as it tavels in a staight line. Although this setup may seem ovesimplified, you will soon see that a detailed
More informationChapter 3 Optical Systems with Annular Pupils
Chapte 3 Optical Systems with Annula Pupils 3 INTRODUCTION In this chapte, we discuss the imaging popeties of a system with an annula pupil in a manne simila to those fo a system with a cicula pupil The
More informationA Relativistic Electron in a Coulomb Potential
A Relativistic Electon in a Coulomb Potential Alfed Whitehead Physics 518, Fall 009 The Poblem Solve the Diac Equation fo an electon in a Coulomb potential. Identify the conseved quantum numbes. Specify
More information