A Simplified Dynamics Block Diagram for a Four-Axis Stabilized Platform

Size: px
Start display at page:

Download "A Simplified Dynamics Block Diagram for a Four-Axis Stabilized Platform"

Transcription

1 A Simplified Dynamic Block Diagram for a FourAxi Stabilized Platform Author: Hendrik Daniël Mouton a Univerity of Cape Town, Rondeboch, Cape Town, South Africa, 770 Abtract: t i relatively traightforward to derive comprehenive dynamic block diagram for two and threeaxi balanced tabilized platform. But trying to do it for four and higher order ytem i getting rapidly uch an amount of work, and the reulting block diagram o extenive, that enible implification from the tart if feaible could be very beneficial. Thoe term that are going to be inignificantly mall relative to the other term hould be left out right from the tart when deriving the equation. n thi paper a method i hown to implify the dynamic model of tabilized platform, with the fouraxi platform a the main example. The relevant equation are derived for the fouraxi example, the block diagram i compiled from the equation after adding imple control loop around the dynamic model to get to the final expanded block diagram, and then ome imulation reult are hown. The method i validated by howing that the imulation reult for a threeaxi ytem are exactly the ame for the comprehenive and implified block diagram, uing the ame method. t i propoed that it i therefore reaonable to expect imilar finding for the fouraxi cae. a Aociate Profeor H D Mouton, Mechanical Engineering, UCT, hennie.mouton@uct.ac.za

2 ntroduction: Variou ort of camera could typically be mounted on tabilized platform to give much improved quality image. t would provide tabilization againt rotational movement of the bae a for example naturally experienced during linear movement of the bae (land craft or aircraft or overhead porting camera mount). Reference [] and [] give very good overall view of major iue in tabilized tracking deign. Deigning uch platform require proper dynamic modelling of the mechanical configuration. t i not difficult to find the typical relevant equation to repreent the dynamic reference [3] and [4] are typical example, but converting them to enible block diagram to ait in imulation and to more eaily ee what i happening i not that readily available. So thi paper contain block diagram, being derived from baic dynamic equation uch a given in reference [5], but the main purpoe of the paper i to demontrate a method to develop implified but till valid block diagram. A threeaxi example: To demontrate the poible validity of the implification method to be given in thi paper, it i firt demontrated with a threeaxi tabilized platform. The derivation of the equation for the threeaxi cae i not going to be given here only the definition, the reaoning for ignoring certain term, a block diagram including control loop and then imulation reult. The block diagram will clearly indicate the comprehenive and implified model. The reult will how that for a uitable tet cenario the implified dynamic model give the ame reult a the comprehenive model. f the latter reult are indeed very imilar, it will be aumed to be valid to apply the ame method of implification to four and moreaxi tabilized platform. Of coure different configuration are poible for the threeaxi cae. But the demontration i going to be done on a pecific choen example. The aumption i that other configuration will give imilar reult. An outer to inner yawpitchroll configuration wa choen. t i aumed that the control loop in the end will provide good tabilization of the platform in roll, yaw and pitch a made poible by the three axe. Nomenclature for the threeaxi cae: From outer to inner, thee are the Euler angle and ubcript to be ued:

3 B will be ued a ubcript for the bae. ψ i the yaw angle between the bae and the outer yaw gimbal yaw will be ued a ubcript. θ i the pitch angle between the yaw gimbal and the pitch gimbal pit will be ued a ubcript. φ i the roll angle between the pitch gimbal and the platform to be tabilized rol will be ued a the ubcript. ω denote inertial angular rate. x, y and z correpond baically to the axe about which roll, pitch and yaw rotation happen repectively. The following picture hould clarify the configuration and the ymbol. X pit, X rol Y yaw, Y pit Z B, Z yaw Figure. Threeaxi configuration Small/large ignal ditinction for three axe: Becaue of the aumption above, the following ditinction between potentially large and mall ignal can be made. Table. Small/large ignal for threeaxi configuration ω xrol ω yrol ω zrol ω xpit ω ypit ω zpit ω xyaw ω yyaw ω zyaw Small Large

4 Becaue it i aumed that the tabilization i going to be good, diturbance torque can be ignored that contain multiplication with any of the ignal in thi table that are claified a Small. Block diagram and imulation reult for three axe: Thi lead to comprehenive and implified block diagram. The implified block diagram i the ame a the comprehenive one but with all block coloured with grey omitted ee the firt block diagram to follow. By looking at the increae in diturbance torque when modelling the dynamic of the three axe from the inner axi to the outer axi a repreented in the comprehenive block diagram, an idea can be formed how expanded comprehenive four and fiveaxi platform model would be. Rather imple tabilization loop are added in the block diagram to follow, with ω g ζ g ω g ω g repreenting feedback from gyro mounted on the platform, meauring ω xrol, ω yrol and ω zrol. A proportional controller and one lead compenator cloe each loop. The command input to the three loop are zero becaue the intention of the loop i to control ω xrol, ω yrol and ω zrol at zero depite the bae motion ω xb, ω yb and ω zb applied to the ytem. They are applied one after the other but overlapping a hown in the firt graph below. Their amplitude are.0 rad/ and their frequencie linearly increae from.0 Hz to 50.0 Hz. n the later graph they will be hown on top of one another imply becaue many other ignal mut be hown on the ame graph.

5 Figure. Applied bae motion The amplitude of the meaurement ω xrol, ω yrol and ω zrol are indication of the quality of tabilization the maller thee ignal the better the tabilization. ω zrol i ignificantly larger than ω xrol and ω yrol a can be expected from uch a threeaxi ytem. t will later be een that the fouraxi ytem provide much better tabilization with baically the ame control loop, but that apect i not the focu of thi paper. The initial yaw, pitch and roll angle, ψ, θ and φ, were et at 45 to prevent ome diturbance torque to be very mall becaue of mall angle. f ome diturbance torque were mall becaue of mall angle it would have defied the validation of the implification propoed here. Remember: n the block diagram below the implified block diagram i the ame a the comprehenive one but with all block coloured with grey omitted.

6 gg g g xrol K T rol rol [ m] 00[ rad/ ] g 0.7 g 0.0 K rol T rol 0.T rol T xrol xrol zrol yrol xrol zrol in co in co yrol ypit xpit in co in co K T K T pit pit 3.0 [ m] yaw yaw [ m] zyaw in co xyaw xb zpit xpit in co yb yyaw xyaw zyaw Typitrol Tzyawpitrol zb T pit 0.T pit yzpitrol zpit zpit xyaw xyaw zxyawpitrol T yaw 0.T yaw K pit xpit zpit ( zrol )co xrol xpit zrol yyaw yyaw yzyawpitrol yyaw xyaw Kyaw 0.0 xrol ( yrol xrol )in ( xpit )co 0. 0 xpit ypit yrol ypit yrol co zrol ( zpit ypit)in zpit in zrol ypitrol ypit yrol yrol xrol zrol yrol zrol ( xrol )coco yrol ( zrol )inco xrol ( yrol )in zrol gg g g yrol yyaw xpit zpit co zpitrol in xpitrol xrol yrol zrol xpit ypit zpit xyaw yyaw zyaw 3.0e 4 [ kg. 5.0e 4 [ kg. 5.0e 4 [ kg..0e 4 [ kg..0e 4 [ kg. 3.0e 4 [ kg. 4.0e 4 [ kg. 3.3e 4 [ kg. 3.0e 4 [ kg. xpitrol ypitrol zpitrol yzpitrol zyawpitrol zxyawpitrol yzyawpitro l xpit ypit zpit yrol zyaw xrol zrol zpitrol yrol yzpitrol co co zrol in zpitrol co co xpitrol zrol yrol in in in xpitrol in yrol zrol co co zrol in co yrol zyaw yzyawpitrol xrol in zyawpitrol zb zrol gg g g Figure 3. Comprehenive and implified block diagram of threeaxi yawpitchroll

7 The reult for the comprehenive and implified model are hown below, demontrating that the propoed implification are valid ince the reult are identical. Figure 4. Comprehenive threeaxi graph reult

8 Figure 5. Simplified threeaxi graph reult

9 The choen fouraxi example: Like for all axe, different configuration are alo poible for the fouraxi cae. The principle to be dicued here i applicable to all configuration. An outer to inner yawpitchyawroll configuration wa choen. t i aumed that the control loop in the end will provide good tabilization of the platform in roll, yaw and pitch a made poible by the three inner axe. The outer yaw axi i only for rough alignment typically not controlled by an inertial rate loop like the inner three axe, but by a relative rate loop enuring that the inner yaw gimbal angle remain relatively mall. Nomenclature for the fouraxi cae: From outer to inner, thee are the Euler angle and ubcript to be ued: B will be ued a ubcript for the bae. β i the yaw angle between the bae and the outer yaw gimbal Yaw will be ued a ubcript. θ i the pitch angle between the outer yaw gimbal and the pitch gimbal pit will be ued a ubcript. ψ i the yaw angle between the pitch gimbal and the inner yaw gimbal yaw will be ued a ubcript. φ i the roll angle between the inner yaw gimbal and the platform to be tabilized rol will be ued a the ubcript. ω denote inertial angular rate. x, y and z correpond baically to the axe about which roll, pitch and yaw happen repectively. The following picture hould clarify the configuration and the ymbol.

10 Y Yaw, Y pit Z pit, Z yaw X yaw, X rol Z B, Z Yaw Figure 6. Fouraxi configuration Small/large ignal ditinction for four axe: Becaue of the aumption above, the following ditinction between potentially large and mall ignal can be made. Table. Small/large ignal for fouraxi configuration ω xrol ω yrol ω zrol ω xyaw ω yyaw ω zyaw ω xpit ω ypit ω zpit ω xyaw ω yyaw ω zyaw Small Large Becaue it i aumed that the tabilization i going to be good, diturbance torque can be ignored that contain the term that were claified a Small in the table above. Derivation of equation for the four axe: Apply T = d dt a (W) ω a W with T the net torque and W = ω to the platform (rol), the inner yaw gimbal (yaw), the pitch gimbal (pit) and the outer yaw gimbal (Yaw). d dt a refer to the derivative w.r.t. a noninertial axe ytem. ω a i the inertial angular rate vector of the noninertial axe ytem.

11 ω i the inertial angular rate vector of the body under conideration. i the Moment of nertia matrix of the body under conideration. To the platform ( rol ): xrol 0 0 W = [ 0 yrol 0 ] [ 0 0 zrol ω xrol ω yrol ω xrol ω yrol ] ω a = [ ] ω zrol ω zrol t i aumed that the ymmetry i uch that the inertia matrix can be repreented in the diagonal form. T xrol T rol = [ T yrol ] = [ T zrol xrol ω xrol yrol ω yrol zrol ω zrol ω yrol ω zrol ( zrol yrol ) ] [ ω zrol ω xrol ( xrol zrol ) ω xrol ω yrol ( yrol xrol ) ] T xrol xrol ω xrol () T yrol yrol ω yrol T zrol zrol ω zrol To the inner yaw gimbal ( yaw ): xyaw 0 0 ω xyaw ω xyaw W = [ 0 yyaw 0 ] [ ω yyaw ] ω a = [ ω yyaw ] 0 0 zyaw ω zyaw ω zyaw t i aumed that the ymmetry i uch that the inertia matrix can be repreented in the diagonal form. T xyaw xyaw ω xyaw ω yyaw ω zyaw ( zyaw yyaw ) T yaw = [ T yyaw ] = [ yyaw ω yyaw] [ ω zyaw ω xyaw ( xyaw zyaw )] T zyaw zyaw ω zyaw ω xyaw ω yyaw ( yyaw xyaw ) T xyaw xyaw ω xyaw T yyaw yyaw ω yyaw T zyaw zyaw ω zyaw T yaw i the net torque applied to the inner yaw gimbal. Call the torque applied to thi gimbal from the outide T yawrol

12 T yaw = T yawrol T rol T yawrol = T yaw T rol Now take the roll angle φ between the yaw and rol axe into account: Y yaw Y rol Z rol Z yaw Figure 7. Roll angle tranformation T xyawrol T xyaw T xrol T yawrol = [ T yyawrol ] = [ T yyaw T yrol co φ T zrol in φ] T zyawrol T zyaw T zrol co φ T yrol in φ T xyawrol xyaw ω xyaw xrol ω xrol T yyawrol yyaw ω yyaw yrol ω yrol co φ zrol ω zrolin φ T zyawrol zyaw ω zyaw zrol ω zrol co φ yrol ω yrolin φ The following are the equation for the angular rate and acceleration: ω xrol = ω xyaw φ ω yrol = ω yyaw co φ ω zyaw in φ () ω zrol = ω zyaw co φ ω yyaw in φ ω xrol = ω xyaw φ ω yrol = ω yyaw co φ ω zyaw in φ φ ω zrol ω zrol = ω zyaw co φ ω yyaw in φ φ ω yrol Subtitute ω xrol, ω yrol and ω zrol in T xyawrol, T yyawrol and T zyawrol. T xyawrol ω xyaw( xyaw xrol ) φ xrol T yyawrol ω yyaw( yyaw yrol co Φ zrol in φ) ω zyaw( yrol zrol ) in φ co φ φ ω zrol yrol co φ φ ω yrol zrol in φ

13 T zyawrol ω zyaw( zyaw zrol co Φ yrol in φ) ω yyaw( yrol zrol ) in φ co φ φ ω yrol zrol co φ φ ω zrol yrol in φ Define: xyawrol = xyaw xrol (3) yyawrol = yyaw yrol co φ zrol in φ (4) zyawrol = zyaw zrol co φ yrol in φ (5) yzyawrol = yrol zrol in φ (6) T xyawrol ω xyaw xyawrol φ xrol T yyawrol ω yyaw yyawrol ω zyaw yzyawrol φ ω zrol yrol co φ φ ω yrol zrol in φ T zyawrol ω zyaw zyawrol ω yyaw yzyawrol φ ω yrol zrol co φ φ ω zrol yrol in φ (7) To the pitch gimbal ( pit ): xpit 0 0 ω xpit ω xpit W = [ 0 ypit 0 ] [ ω ypit ] ω a = [ ω ypit ] 0 0 zpit ω zpit ω zpit t i aumed that the ymmetry i uch that the inertia matrix can be repreented in the diagonal form. T xpit xpit ω xpit ω ypit ω zpit ( zpit ypit ) T pit = [ T ypit ] = [ ypit ω ypit] [ ω zpit ω xpit ( xpit zpit )] T zpit zpit ω zpit ω xpit ω ypit ( ypit xpit ) T xpit xpit ω xpit T ypit = ypit ω ypit ω zpit ω xpit ( xpit zpit ) T zpit zpit ω zpit T pit i the net torque applied to the pitch gimbal. Call the torque applied to thi gimbal from the outide T pityawrol. T pit = T pityawrol T yawrol T pityawrol = T pit T yawrol

14 Now take the yaw angle ψ between the pit and yaw axe into account: X pit X yaw Y yaw Y pit Figure 8. Yaw angle tranformation T xpityawrol T xpit T xyawrol co ψ T yyawrol in ψ T pityawrol = [ T ypityawrol ] = [ T yrol T yyawrol co ψ T xyawrol in ψ] T zpityawrol T zpit T zyawrol T xpityawrol xpit ω xpit (ω xyaw xyawrol φ xrol ) co ψ (ω yyaw yyawrol ω zyaw yzyawrol φ ω zrol yrol co φ φ ω yrol zrol in φ) in ψ xpit ω xpit (ω xyaw xyawrol φ xrol ) co ψ (ω yyaw yyawrol ω zyaw yzyawrol ) in ψ with φ ω yrol and φ ω zrol ignored. T ypityawrol ypit ω ypit ω zpit ω xpit ( xpit zpit ) (ω yyaw yyawrol ω zyaw yzyawrol φ ω zrol yrol co φ φ ω yrol zrol in φ) co ψ (ω xyaw xyawrol φ xrol ) in ψ ypit ω ypit ω zpit ω xpit ( xpit zpit ) (ω yyaw yyawrol ω zyaw yzyawrol ) co ψ (ω xyaw xyawrol φ xrol ) in ψ with φ ω yrol and φ ω zrol ignored. T zpityawrol zpit ω zpit ω zyaw zyawrol ω yyaw yzyawrol φ ω yrol zrol co φ φ ω zrol yrol in φ zpit ω zpit ω zyaw zyawrol ω yyaw yzyawrol with φ ω yrol and φ ω zrol ignored. The following are the equation for the angular rate and acceleration: ω xyaw = ω xpit co ψ ω ypit in ψ

15 ω yyaw = ω ypit co ψ ω xpit in ψ (8) ω zyaw = ω zpit ψ ω xyaw = ω xpit co ψ ω ypit in ψ ψ ω yyaw ω yyaw = ω ypit co ψ ω xpit in ψ ψ ω xyaw ω zyaw = ω zpit ψ Subtitute ω xyaw, ω yyaw and ω zyaw in T xpityawrol, T ypityawrol and T zpityawrol. T xpityawrol xpit ω xpit (ω xpit co ψ ω ypit in ψ ψ ω yyaw ) xyawrol co ψ (ω ypit co ψ ω xpit in ψ ψ ω xyaw ) yyawrol in ψ (ω zpit ψ ) yzyawrol in ψ φ xrol co ψ = ω xpit( xpit xyawrol co ψ yyawrol in ψ) ω ypit( xyawrol yyawrol ) in ψ co ψ ω zpit yzyawrol in ψ ψ ω yyaw xyawrol co ψ ψ ω xyaw yyawrol in ψ ψ yzyawrol in ψ φ xrol co ψ T ypityawrol ypit ω ypit ω zpit ω xpit ( xpit zpit ) (ω ypit co ψ ω xpit in ψ ψ ω xyaw ) yyawrol co ψ (ω xpit co ψ ω ypit in ψ ψ ω yyaw ) xyawrol in ψ (ω zpit ψ ) yzyawrol co ψ φ xrol in ψ = ω ypit( ypit yyawrol co ψ xyawrol in ψ) ω zpit yzyawrol co ψ ω xpit( xyawrol yyawrol ) in ψ co ψ ω zpit ω xpit ( xpit zpit ) ψ ω xyaw yyawrol co ψ ψ ω yyaw xyawrol in ψ ψ yzyawrol co ψ φ xrol in ψ T zpityawrol zpit ω zpit (ω zpit ψ ) zyawrol (ω ypit co ψ ω xpit in ψ ψ ω xyaw ) yzyawrol = ω zpit( zpit zyawrol ) ω xpit yzyawrol in ψ ω ypit yzyawrol co ψ ψ ω xyaw yzyawrol ψ zyawrol Define: xpityawrol = xpit xyawrol co ψ yyawrol in ψ (9) ypityawrol = ypit yyawrol co ψ xyawrol in ψ (0) zpityawrol = zpit zyawrol () xypityawrol = xyawrol yyawrol in ψ () yzpityawrol = yzyawrol co ψ (3) zxpityawrol = yzyawrol in ψ (4)

16 T xpityawrol ω xpit xpityawrol ω ypit xypityawrol ω zpit zxpityawrol ψ ω yyaw xyawrol co ψ ψ ω xyaw yyawrol in ψ ψ yzyawrol in ψ φ xrol co ψ T ypityawrol ω ypit ypityawrol ω zpit yzpityawrol ω xpit xypityawrol ω zpit ω xpit ( xpit zpit ) ψ ω xyaw yyawrol co ψ ψ ω yyaw xyawrol in ψ ψ yzyawrol co ψ φ xrol in ψ (5) T zpityawrol ω zpit zpityawrol ω xpit zxpityawrol ω ypit yzpityawrol ψ ω xyaw yzyawrol ψ zyawrol To the outer yaw gimbal (Yaw): xyaw 0 0 W = [ 0 yyaw 0 ] [ 0 0 zyaw ω xyaw ω yyaw ω xyaw ω yyaw ] ω a = [ ] ω zyaw ω zyaw t i aumed that the ymmetry i uch that the inertia matrix can be repreented in the diagonal form. T xyaw T Yaw = [ T yyaw ] = [ T zyaw xyaw ω xyaw yyaw ω yyaw zyaw ω zyaw ω yyaw ω zyaw ( zyaw yyaw ) ] [ ω zyaw ω xyaw ( xyaw zyaw ) ω xyaw ω yyaw ( yyaw xyaw ) ] T Yaw i the net torque applied to the outer yaw gimbal. Call the torque applied to thi gimbal from the outide T Yawpityawrol. T Yaw = T Yawpityawrol T pityawrol T Yawpityawrol = T Yaw T pityawrol Now take the pitch angle θ between the Yaw and pit axe into account: X pit X Yaw Z pit Z Yaw Figure 9. Pitch angle tranformation

17 T xyawpityawrol T xyaw T xpityawrol co θ T zpityawrol in θ T Yawpityawrol = [ T yyawpityawrol ] = [ T yyaw T ypityawrol ] T zyawpityawrol T zyaw T zpityawrol co θ T xpityawrol in θ T zyawpityawrol zyaw ω zyaw ω xyaw ω yyaw ( yyaw xyaw ) (ω zpit zpityawrol ω xpit zxpityawrol ω ypit yzpityawrol ψ ω xyaw yzyawrol ψ zyawrol ) co θ (ω xpit xpityawrol ω ypit xypityawrol ω zpit zxpityawrol ψ ω yyaw xyawrol co ψ ψ ω xyaw yyawrol in ψ ψ yzyawrol in ψ φ xrol co ψ) in θ zyaw ω zyaw ω xyaw ω yyaw ( yyaw xyaw ) (ω zpit zpityawrol ω xpit zxpityawrol ω ypit yzpityawrol ψ ω xyaw yzyawrol ψ zyawrol ) co θ (ω xpit xpityawrol ω ypit xypityawrol ω zpit zxpityawrol ψ ω xyaw yyawrol in ψ ψ yzyawrol in ψ φ xrol co ψ) in θ with ψ ω yyaw ignored. The following are the equation for the angular rate and acceleration: ω xpit = ω xyaw co θ ω zyaw in θ ω ypit = ω yyaw θ (6) ω zpit = ω zyaw co θ ω xyaw in θ ω xpit = ω xyaw co θ ω zyaw in θ θ ω zpit ω ypit = ω zyaw θ ω zpit = ω zyaw co θ ω xyaw in θ θ ω xpit Subtitute ω xpit, ω ypit and ω zpit in T zyawpityawrol. T zyawpityawrol zyaw ω zyaw ω xyaw ω yyaw ( yyaw xyaw ) (ω zyaw co θ ω xyaw in θ θ ω xpit ) zpityawrol co θ (ω zyaw co θ ω xyaw in θ θ ω xpit ) zxpityawrol in θ (ω xyaw co θ ω zyaw in θ θ ω zpit ) zxpityawrol co θ (ω xyaw co θ ω zyaw in θ θ ω zpit ) xpityawrol in θ (ω yyaw θ ) yzpityawrol co θ (ω yyaw θ ) xypityawrol in θ ψ ω xyaw ( yzyawrol co θ yyawrol in ψ in θ) ψ ( zyawrol co θ yzyawrol in ψ in θ) φ xrol co ψ in θ = ω zyaw( zyaw zpityawrol co θ zxpityawrol in θ xpityawrol in θ) ω xyaw[ zxpityawrol (co θ in θ) ( xpityawrol zpityawrol ) in θ co θ] ω yyaw( yzpityawrol co θ xypityawrol in θ)

18 ω xyaw ω yyaw ( yyaw xyaw ) θ ω xpit ( zpityawrol co θ zxpityawrol in θ) θ ω zpit ( zxpityawrol co θ xpityawrol in θ) ψ ω xyaw ( yzyawrol co θ yyawrol in ψ in θ) θ ( yzpityawrol co θ xypityawrol in θ) ψ ( zyawrol co θ yzyawrol in ψ in θ) φ xrol co ψ in θ Define: zyawpityawrol = zyaw zpityawrol co θ zxpityawrol in θ xpityawrol in θ (7) zxyawpityawrol = zxpityawrol co θ xpityawrol zpityawrol in θ (8) yzyawpityawrol = yzpityawrol co θ xypityawrol in θ (9) T zyawpityawrol ω zyaw zyawpityawrol ω xyaw zxyawpityawrol ω yyaw yzyawpityawrol ω xyaw ω yyaw ( yyaw xyaw ) θ ω xpit ( zpityawrol co θ zxpityawrol in θ) θ ω zpit ( zxpityawrol co θ xpityawrol in θ) ψ ω xyaw ( yzyawrol co θ yyawrol in ψ in θ) θ ( yzpityawrol co θ xypityawrol in θ) ψ ( zyawrol co θ yzyawrol in ψ in θ) φ xrol co ψ in θ (0) Now take the yaw angle β between the B and Yaw axe into account: X B X Yaw Y Yaw Y B Figure 0. Outer yaw angle tranformation ω xyaw = ω xb co β ω yb in β ω yyaw = ω yb co β ω xb in β () ω zyaw = ω zb β

19 Block diagram and imulation reult for four axe: From equation () to () and with the addition of friction, the following block diagram can be compiled. Rather imple tabilization loop are added in the block diagram, with ω g ζ g ω g ω g repreenting feedback from gyro mounted on the platform, meauring ω xrol, ω zrol and ω yrol. A proportional controller and one lead compenator cloe each loop. The outer yaw gimbal ha a tachometer, repreented by ω t ζ t ω t ω t, feeding back the relative angular rate β and alo a proportional controller and one lead compenator. t i driven by the inner yaw gimbal angle ψ through a gain K ctrl to effectively keep ψ at a mall angle therefore aiding good tabilization. The ize of the meaurement ω xrol, ω yrol and ω zrol i an indication of the quality of tabilization the maller thee ignal the better the tabilization. They can be een in the graph below to be much maller than the previou threeaxi cae a expected. The initial Yaw, pitch and roll angle β, θ and φ were et at 45 to prevent ome diturbance torque to be very mall becaue of mall angle.

20 gg g g xrol K T rol rol [ m] K ctrl.0[( rad/ )/ rad] 00[ rad/ ] g 0.7 g K rol T rol 0.T rol K yaw Tyaw 0.T yaw T xrol T ypityawrol gg g T zyawrol g xrol xrol zrol in co in co yrol zyaw yyaw zyaw zyawrol zpit yyaw yyaw yzyawrol xyaw yyaw in co in co xyaw ypit xyaw yyaw K T K T yaw yaw pit pit xpit [ m] [ m] xpit zpit xyaw yyaw K T in co in co xyaw Yaw Yaw [ m] zyaw zyaw zxyawpityawrol yzyawpityawrol yyaw xyaw K ctrl yyaw in co in co T zyawpityawrol 00[ rad/ ] t 0.7 t T yaw 0.T yaw K yaw xb yb zb T pit 0.T pit K pit yzpityawrol zpit xpit xypityawro l xpit zpit xpit zpit co in zpityawrol zxpityawro l zxpityawrol co in xpityawrol 0.0 xpit zpit co yyawrol yzpityawro l in xrol xyaw xyaw yzyawrol yzpityawrol zyawrol co yyawrol in in co in xypityawrol co yzyawrolin in co in xrol ypityawrol ypit zyaw zyawpityawrol zb gg g g yrol yyaw t tt t xrol yrol zrol xyaw yyaw zyaw 3.0e 4 [ kg. 5.0e 4 [ kg. 5.0e 4 [ kg..0e 4 [ kg. 3.0e 4 [ kg..0e 4 [ kg. xpit ypit zpit xyaw yyaw zyaw 4.0e 4 [ kgm. ] 3.0e 4 [ kgm. ] 3.3e 4 [ kgm. ] 5.0e 4 [ kgm. ] 4.4e 4 [ kgm. ] 4.0e 4 [ kgm. ] xyawrol yyawrol zyawrol yzyawrol xyaw yyaw zyaw yrol xrol yrol zrol zrol co co in zrol yrol in in xpityawrol ypityawrol zpityawrol xypityawrol yzpityawrol zxpityawrol xpit ypit zpit xyawrol yzyawrol yzyawrol xyawrol yyawrol zyawrol co in co co yyawrol in yyawrol xyawrol in in zyawpityawrol zxyawpitya wrol yzyawpitya wrol zyaw zxpityawro l yzpityawro l zpityawrol co co co xpityawrol xypityawro l zxpityawro l in in zpityawrol in xpityawrol in Figure. Simplified block diagram of fouraxi yawpitchyawroll

21 The imulation reult i the following. Figure. Simplified fouraxi graph reult A could be expected, ω xrol, ω yrol and ω zrol, which are indication of the quality of tabilization, are much maller than for the threeaxi cae. Concluion: Baed on the demontration with the threeaxi example where the comprehenive and implified model were compared and found to give identical reult, and the reaoning for ignoring certain term when deriving the implified equation and model, it i concluded that the propoed implification method i valid alo for moreaxi, like the fouraxi configuration. The aving in the amount of work to be done to derive model for four and higher number of axe platform by uing thi method i ubtantial a can be appreciated a bit by comparing the comprehenive and implified threeaxi block diagram. Actually the amount of work to be aved will increae dramatically a the number of axe i increaed.

22 Summary: n thi paper a method wa hown to implify the dynamic model of tabilized platform, with the fouraxi platform a the example vehicle. The relevant implified equation for a pecific fouraxi configuration were derived, the reulting block diagram wa given, and then ome imulation reult were hown after adding imple control loop around the dynamic model. The validity of the implification method wa demontrated by howing that identical imulation reult were achieved for a comprehenive and implified threeaxi model example. Reference: [] Michael K. Maten nertially Stabilized Platform for Optical maging Sytem Tracking Dynamic Target with Mobile Senor EEE CONTROL SYSTEMS MAGAZNE, February 008, Digital Object dentifier 0.09/MCS [] J. M. Hilkert nertially Stabilized Platform Technology Concept and Principle EEE CONTROL SYSTEMS MAGAZNE, February 008, Digital Object dentifier 0.09/MCS [3] Edmund Gorczycki Dynamic Conideration Relating to the Behavior of nertial SpaceStabilized Platform Journal of the Aeronautical Science (ntitute of the Aeronautical Science), 957, Vol.4: 3038, 0.54/ [4] Melvin Feinberg Rigid Body Dynamic Stable Platform olation Sytem Guidance and Control, 964: 4437, 0.54/ [5] Murray R Spiegel

23 Schaum Outline of Theory and Problem of Theoretical Mechanic, pp. 538 Schaum Publihing Company, 967

Gain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays

Gain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays Gain and Phae Margin Baed Delay Dependent Stability Analyi of Two- Area LFC Sytem with Communication Delay Şahin Sönmez and Saffet Ayaun Department of Electrical Engineering, Niğde Ömer Halidemir Univerity,

More information

Social Studies 201 Notes for November 14, 2003

Social Studies 201 Notes for November 14, 2003 1 Social Studie 201 Note for November 14, 2003 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the

More information

EE Control Systems LECTURE 14

EE Control Systems LECTURE 14 Updated: Tueday, March 3, 999 EE 434 - Control Sytem LECTURE 4 Copyright FL Lewi 999 All right reerved ROOT LOCUS DESIGN TECHNIQUE Suppoe the cloed-loop tranfer function depend on a deign parameter k We

More information

Advanced D-Partitioning Analysis and its Comparison with the Kharitonov s Theorem Assessment

Advanced D-Partitioning Analysis and its Comparison with the Kharitonov s Theorem Assessment Journal of Multidiciplinary Engineering Science and Technology (JMEST) ISSN: 59- Vol. Iue, January - 5 Advanced D-Partitioning Analyi and it Comparion with the haritonov Theorem Aement amen M. Yanev Profeor,

More information

Social Studies 201 Notes for March 18, 2005

Social Studies 201 Notes for March 18, 2005 1 Social Studie 201 Note for March 18, 2005 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the

More information

CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS

CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.1 INTRODUCTION 8.2 REDUCED ORDER MODEL DESIGN FOR LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.3

More information

GNSS Solutions: What is the carrier phase measurement? How is it generated in GNSS receivers? Simply put, the carrier phase

GNSS Solutions: What is the carrier phase measurement? How is it generated in GNSS receivers? Simply put, the carrier phase GNSS Solution: Carrier phae and it meaurement for GNSS GNSS Solution i a regular column featuring quetion and anwer about technical apect of GNSS. Reader are invited to end their quetion to the columnit,

More information

A Simplified Methodology for the Synthesis of Adaptive Flight Control Systems

A Simplified Methodology for the Synthesis of Adaptive Flight Control Systems A Simplified Methodology for the Synthei of Adaptive Flight Control Sytem J.ROUSHANIAN, F.NADJAFI Department of Mechanical Engineering KNT Univerity of Technology 3Mirdamad St. Tehran IRAN Abtract- A implified

More information

Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog

Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog Chapter Sampling and Quantization.1 Analog and Digital Signal In order to invetigate ampling and quantization, the difference between analog and digital ignal mut be undertood. Analog ignal conit of continuou

More information

Lecture 10 Filtering: Applied Concepts

Lecture 10 Filtering: Applied Concepts Lecture Filtering: Applied Concept In the previou two lecture, you have learned about finite-impule-repone (FIR) and infinite-impule-repone (IIR) filter. In thee lecture, we introduced the concept of filtering

More information

POWER SYSTEM SMALL SIGNAL STABILITY ANALYSIS BASED ON TEST SIGNAL

POWER SYSTEM SMALL SIGNAL STABILITY ANALYSIS BASED ON TEST SIGNAL POWE YEM MALL INAL ABILIY ANALYI BAE ON E INAL Zheng Xu, Wei hao, Changchun Zhou Zheang Univerity, Hangzhou, 37 PChina Email: hvdc@ceezueducn Abtract - In thi paper, a method baed on ome tet ignal (et

More information

S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS

S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS by Michelle Gretzinger, Daniel Zyngier and Thoma Marlin INTRODUCTION One of the challenge to the engineer learning proce control i relating theoretical

More information

The Hassenpflug Matrix Tensor Notation

The Hassenpflug Matrix Tensor Notation The Haenpflug Matrix Tenor Notation D.N.J. El Dept of Mech Mechatron Eng Univ of Stellenboch, South Africa e-mail: dnjel@un.ac.za 2009/09/01 Abtract Thi i a ample document to illutrate the typeetting of

More information

Finite Element Analysis of a Fiber Bragg Grating Accelerometer for Performance Optimization

Finite Element Analysis of a Fiber Bragg Grating Accelerometer for Performance Optimization Finite Element Analyi of a Fiber Bragg Grating Accelerometer for Performance Optimization N. Baumallick*, P. Biwa, K. Dagupta and S. Bandyopadhyay Fiber Optic Laboratory, Central Gla and Ceramic Reearch

More information

Design By Emulation (Indirect Method)

Design By Emulation (Indirect Method) Deign By Emulation (Indirect Method he baic trategy here i, that Given a continuou tranfer function, it i required to find the bet dicrete equivalent uch that the ignal produced by paing an input ignal

More information

SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD

SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD S.P. Teeuwen, I. Erlich U. Bachmann Univerity of Duiburg, Germany Department of Electrical Power Sytem

More information

Digital Control System

Digital Control System Digital Control Sytem - A D D A Micro ADC DAC Proceor Correction Element Proce Clock Meaurement A: Analog D: Digital Continuou Controller and Digital Control Rt - c Plant yt Continuou Controller Digital

More information

Homework 12 Solution - AME30315, Spring 2013

Homework 12 Solution - AME30315, Spring 2013 Homework 2 Solution - AME335, Spring 23 Problem :[2 pt] The Aerotech AGS 5 i a linear motor driven XY poitioning ytem (ee attached product heet). A friend of mine, through careful experimentation, identified

More information

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is EE 4G Note: Chapter 6 Intructor: Cheung More about ZSR and ZIR. Finding unknown initial condition: Given the following circuit with unknown initial capacitor voltage v0: F v0/ / Input xt 0Ω Output yt -

More information

Control Systems Analysis and Design by the Root-Locus Method

Control Systems Analysis and Design by the Root-Locus Method 6 Control Sytem Analyi and Deign by the Root-Locu Method 6 1 INTRODUCTION The baic characteritic of the tranient repone of a cloed-loop ytem i cloely related to the location of the cloed-loop pole. If

More information

7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281

7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281 72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition

More information

Bernoulli s equation may be developed as a special form of the momentum or energy equation.

Bernoulli s equation may be developed as a special form of the momentum or energy equation. BERNOULLI S EQUATION Bernoulli equation may be developed a a pecial form of the momentum or energy equation. Here, we will develop it a pecial cae of momentum equation. Conider a teady incompreible flow

More information

Analysis of Prevention of Induction Motors Stalling by Capacitor Switching

Analysis of Prevention of Induction Motors Stalling by Capacitor Switching 16th NTIONL POWER SYSTEMS CONFERENCE, 15th-17th DECEMER, 2010 260 nalyi of Prevention of Induction Motor Stalling by Capacitor Switching S.Maheh and P.S Nagendra rao Department of Electrical Engineering

More information

An estimation approach for autotuning of event-based PI control systems

An estimation approach for autotuning of event-based PI control systems Acta de la XXXIX Jornada de Automática, Badajoz, 5-7 de Septiembre de 08 An etimation approach for autotuning of event-baed PI control ytem Joé Sánchez Moreno, María Guinaldo Loada, Sebatián Dormido Departamento

More information

Solutions. Digital Control Systems ( ) 120 minutes examination time + 15 minutes reading time at the beginning of the exam

Solutions. Digital Control Systems ( ) 120 minutes examination time + 15 minutes reading time at the beginning of the exam BSc - Sample Examination Digital Control Sytem (5-588-) Prof. L. Guzzella Solution Exam Duration: Number of Quetion: Rating: Permitted aid: minute examination time + 5 minute reading time at the beginning

More information

A Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems

A Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems A Contraint Propagation Algorithm for Determining the Stability Margin of Linear Parameter Circuit and Sytem Lubomir Kolev and Simona Filipova-Petrakieva Abtract The paper addree the tability margin aement

More information

arxiv: v2 [nucl-th] 3 May 2018

arxiv: v2 [nucl-th] 3 May 2018 DAMTP-207-44 An Alpha Particle Model for Carbon-2 J. I. Rawlinon arxiv:72.05658v2 [nucl-th] 3 May 208 Department of Applied Mathematic and Theoretical Phyic, Univerity of Cambridge, Wilberforce Road, Cambridge

More information

Moment of Inertia of an Equilateral Triangle with Pivot at one Vertex

Moment of Inertia of an Equilateral Triangle with Pivot at one Vertex oment of nertia of an Equilateral Triangle with Pivot at one Vertex There are two wa (at leat) to derive the expreion f an equilateral triangle that i rotated about one vertex, and ll how ou both here.

More information

Correction for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002

Correction for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002 Correction for Simple Sytem Example and Note on Laplace Tranform / Deviation Variable ECHE 55 Fall 22 Conider a tank draining from an initial height of h o at time t =. With no flow into the tank (F in

More information

Chapter 13. Root Locus Introduction

Chapter 13. Root Locus Introduction Chapter 13 Root Locu 13.1 Introduction In the previou chapter we had a glimpe of controller deign iue through ome imple example. Obviouly when we have higher order ytem, uch imple deign technique will

More information

A Bluffer s Guide to... Sphericity

A Bluffer s Guide to... Sphericity A Bluffer Guide to Sphericity Andy Field Univerity of Suex The ue of repeated meaure, where the ame ubject are teted under a number of condition, ha numerou practical and tatitical benefit. For one thing

More information

Efficient Methods of Doppler Processing for Coexisting Land and Weather Clutter

Efficient Methods of Doppler Processing for Coexisting Land and Weather Clutter Efficient Method of Doppler Proceing for Coexiting Land and Weather Clutter Ça gatay Candan and A Özgür Yılmaz Middle Eat Technical Univerity METU) Ankara, Turkey ccandan@metuedutr, aoyilmaz@metuedutr

More information

Neural Network Linearization of Pressure Force Sensor Transfer Characteristic

Neural Network Linearization of Pressure Force Sensor Transfer Characteristic Acta Polytechnica Hungarica Vol., No., 006 Neural Network Linearization of Preure Force Senor Tranfer Characteritic Jozef Vojtko, Irena Kováčová, Ladilav Madaráz, Dobrolav Kováč Faculty of Electrical Engineering

More information

ME 375 EXAM #1 Tuesday February 21, 2006

ME 375 EXAM #1 Tuesday February 21, 2006 ME 375 EXAM #1 Tueday February 1, 006 Diviion Adam 11:30 / Savran :30 (circle one) Name Intruction (1) Thi i a cloed book examination, but you are allowed one 8.5x11 crib heet. () You have one hour to

More information

What lies between Δx E, which represents the steam valve, and ΔP M, which is the mechanical power into the synchronous machine?

What lies between Δx E, which represents the steam valve, and ΔP M, which is the mechanical power into the synchronous machine? A 2.0 Introduction In the lat et of note, we developed a model of the peed governing mechanim, which i given below: xˆ K ( Pˆ ˆ) E () In thee note, we want to extend thi model o that it relate the actual

More information

A Study on Simulating Convolutional Codes and Turbo Codes

A Study on Simulating Convolutional Codes and Turbo Codes A Study on Simulating Convolutional Code and Turbo Code Final Report By Daniel Chang July 27, 2001 Advior: Dr. P. Kinman Executive Summary Thi project include the deign of imulation of everal convolutional

More information

Coupling of Three-Phase Sequence Circuits Due to Line and Load Asymmetries

Coupling of Three-Phase Sequence Circuits Due to Line and Load Asymmetries Coupling of Three-Phae Sequence Circuit Due to Line and Load Aymmetrie DEGO BELLAN Department of Electronic nformation and Bioengineering Politecnico di Milano Piazza Leonardo da inci 01 Milano TALY diego.ellan@polimi.it

More information

CONTROL SYSTEMS, ROBOTICS AND AUTOMATION Vol. VIII Decoupling Control - M. Fikar

CONTROL SYSTEMS, ROBOTICS AND AUTOMATION Vol. VIII Decoupling Control - M. Fikar DECOUPLING CONTROL M. Fikar Department of Proce Control, Faculty of Chemical and Food Technology, Slovak Univerity of Technology in Bratilava, Radlinkého 9, SK-812 37 Bratilava, Slovakia Keyword: Decoupling:

More information

into a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get

into a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get Lecture 25 Introduction to Some Matlab c2d Code in Relation to Sampled Sytem here are many way to convert a continuou time function, { h( t) ; t [0, )} into a dicrete time function { h ( k) ; k {0,,, }}

More information

Tuning of High-Power Antenna Resonances by Appropriately Reactive Sources

Tuning of High-Power Antenna Resonances by Appropriately Reactive Sources Senor and Simulation Note Note 50 Augut 005 Tuning of High-Power Antenna Reonance by Appropriately Reactive Source Carl E. Baum Univerity of New Mexico Department of Electrical and Computer Engineering

More information

EE 4443/5329. LAB 3: Control of Industrial Systems. Simulation and Hardware Control (PID Design) The Inverted Pendulum. (ECP Systems-Model: 505)

EE 4443/5329. LAB 3: Control of Industrial Systems. Simulation and Hardware Control (PID Design) The Inverted Pendulum. (ECP Systems-Model: 505) EE 4443/5329 LAB 3: Control of Indutrial Sytem Simulation and Hardware Control (PID Deign) The Inverted Pendulum (ECP Sytem-Model: 505) Compiled by: Nitin Swamy Email: nwamy@lakehore.uta.edu Email: okuljaca@lakehore.uta.edu

More information

March 18, 2014 Academic Year 2013/14

March 18, 2014 Academic Year 2013/14 POLITONG - SHANGHAI BASIC AUTOMATIC CONTROL Exam grade March 8, 4 Academic Year 3/4 NAME (Pinyin/Italian)... STUDENT ID Ue only thee page (including the back) for anwer. Do not ue additional heet. Ue of

More information

The Influence of the Load Condition upon the Radial Distribution of Electromagnetic Vibration and Noise in a Three-Phase Squirrel-Cage Induction Motor

The Influence of the Load Condition upon the Radial Distribution of Electromagnetic Vibration and Noise in a Three-Phase Squirrel-Cage Induction Motor The Influence of the Load Condition upon the Radial Ditribution of Electromagnetic Vibration and Noie in a Three-Phae Squirrel-Cage Induction Motor Yuta Sato 1, Iao Hirotuka 1, Kazuo Tuboi 1, Maanori Nakamura

More information

5.5 Application of Frequency Response: Signal Filters

5.5 Application of Frequency Response: Signal Filters 44 Dynamic Sytem Second order lowpa filter having tranfer function H()=H ()H () u H () H () y Firt order lowpa filter Figure 5.5: Contruction of a econd order low-pa filter by combining two firt order

More information

DIFFERENTIAL EQUATIONS

DIFFERENTIAL EQUATIONS DIFFERENTIAL EQUATIONS Laplace Tranform Paul Dawkin Table of Content Preface... Laplace Tranform... Introduction... The Definition... 5 Laplace Tranform... 9 Invere Laplace Tranform... Step Function...4

More information

Department of Mechanical Engineering Massachusetts Institute of Technology Modeling, Dynamics and Control III Spring 2002

Department of Mechanical Engineering Massachusetts Institute of Technology Modeling, Dynamics and Control III Spring 2002 Department of Mechanical Engineering Maachuett Intitute of Technology 2.010 Modeling, Dynamic and Control III Spring 2002 SOLUTIONS: Problem Set # 10 Problem 1 Etimating tranfer function from Bode Plot.

More information

Massachusetts Institute of Technology Dynamics and Control II

Massachusetts Institute of Technology Dynamics and Control II I E Maachuett Intitute of Technology Department of Mechanical Engineering 2.004 Dynamic and Control II Laboratory Seion 5: Elimination of Steady-State Error Uing Integral Control Action 1 Laboratory Objective:

More information

Control of Delayed Integrating Processes Using Two Feedback Controllers R MS Approach

Control of Delayed Integrating Processes Using Two Feedback Controllers R MS Approach Proceeding of the 7th WSEAS International Conference on SYSTEM SCIENCE and SIMULATION in ENGINEERING (ICOSSSE '8) Control of Delayed Integrating Procee Uing Two Feedback Controller R MS Approach LIBOR

More information

ME 375 FINAL EXAM SOLUTIONS Friday December 17, 2004

ME 375 FINAL EXAM SOLUTIONS Friday December 17, 2004 ME 375 FINAL EXAM SOLUTIONS Friday December 7, 004 Diviion Adam 0:30 / Yao :30 (circle one) Name Intruction () Thi i a cloed book eamination, but you are allowed three 8.5 crib heet. () You have two hour

More information

Optimal Coordination of Samples in Business Surveys

Optimal Coordination of Samples in Business Surveys Paper preented at the ICES-III, June 8-, 007, Montreal, Quebec, Canada Optimal Coordination of Sample in Buine Survey enka Mach, Ioana Şchiopu-Kratina, Philip T Rei, Jean-Marc Fillion Statitic Canada New

More information

A PLC BASED MIMO PID CONTROLLER FOR MULTIVARIABLE INDUSTRIAL PROCESSES

A PLC BASED MIMO PID CONTROLLER FOR MULTIVARIABLE INDUSTRIAL PROCESSES ABCM Sympoium Serie in Mechatronic - Vol. 3 - pp.87-96 Copyright c 8 by ABCM A PLC BASE MIMO PI CONOLLE FO MULIVAIABLE INUSIAL POCESSES Joé Maria Galvez, jmgalvez@ufmg.br epartment of Mechanical Engineering

More information

Finding the location of switched capacitor banks in distribution systems based on wavelet transform

Finding the location of switched capacitor banks in distribution systems based on wavelet transform UPEC00 3t Aug - 3rd Sept 00 Finding the location of witched capacitor bank in ditribution ytem baed on wavelet tranform Bahram nohad Shahid Chamran Univerity in Ahvaz bahramnohad@yahoo.com Mehrdad keramatzadeh

More information

SIMPLIFIED MODEL FOR EPICYCLIC GEAR INERTIAL CHARACTERISTICS

SIMPLIFIED MODEL FOR EPICYCLIC GEAR INERTIAL CHARACTERISTICS UNIVERSITY OF PITESTI SCIENTIFIC BULLETIN FACULTY OF ECHANICS AND TECHNOLOGY AUTOOTIVE erie, year XVII, no. ( 3 ) SIPLIFIED ODEL FOR EPICYCLIC GEAR INERTIAL CHARACTERISTICS Ciobotaru, Ticuşor *, Feraru,

More information

Publication V by authors

Publication V by authors Publication Kontantin S. Kotov and Jorma J. Kyyrä. 008. nertion lo and network parameter in the analyi of power filter. n: Proceeding of the 008 Nordic Workhop on Power and ndutrial Electronic (NORPE 008).

More information

Emittance limitations due to collective effects for the TOTEM beams

Emittance limitations due to collective effects for the TOTEM beams LHC Project ote 45 June 0, 004 Elia.Metral@cern.ch Andre.Verdier@cern.ch Emittance limitation due to collective effect for the TOTEM beam E. Métral and A. Verdier, AB-ABP, CER Keyword: TOTEM, collective

More information

At the end of this lesson, the students should be able to understand:

At the end of this lesson, the students should be able to understand: Intructional Objective: At the end of thi leon, the tudent hould be able to undertand: Baic failure mechanim of riveted joint. Concept of deign of a riveted joint. 1. Strength of riveted joint: Strength

More information

Analysis and Design of a Third Order Phase-Lock Loop

Analysis and Design of a Third Order Phase-Lock Loop Analyi Deign of a Third Order Phae-Lock Loop DANIEL Y. ABRAMOVITCH Ford Aeropace Corporation 3939 Fabian Way, MS: X- Palo Alto, CA 94303 Abtract Typical implementation of a phae-lock loop (PLL) are econd

More information

SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuits II. Solutions to Assignment 3 February 2005.

SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuits II. Solutions to Assignment 3 February 2005. SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuit II Solution to Aignment 3 February 2005. Initial Condition Source 0 V battery witch flip at t 0 find i 3 (t) Component value:

More information

Physics 2. Angular Momentum. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Physics 2. Angular Momentum. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB Phyic Angular Momentum For Campu earning Angular Momentum Thi i the rotational equivalent of linear momentum. t quantifie the momentum of a rotating object, or ytem of object. To get the angular momentum,

More information

USING NONLINEAR CONTROL ALGORITHMS TO IMPROVE THE QUALITY OF SHAKING TABLE TESTS

USING NONLINEAR CONTROL ALGORITHMS TO IMPROVE THE QUALITY OF SHAKING TABLE TESTS October 12-17, 28, Beijing, China USING NONLINEAR CONTR ALGORITHMS TO IMPROVE THE QUALITY OF SHAKING TABLE TESTS T.Y. Yang 1 and A. Schellenberg 2 1 Pot Doctoral Scholar, Dept. of Civil and Env. Eng.,

More information

Evolutionary Algorithms Based Fixed Order Robust Controller Design and Robustness Performance Analysis

Evolutionary Algorithms Based Fixed Order Robust Controller Design and Robustness Performance Analysis Proceeding of 01 4th International Conference on Machine Learning and Computing IPCSIT vol. 5 (01) (01) IACSIT Pre, Singapore Evolutionary Algorithm Baed Fixed Order Robut Controller Deign and Robutne

More information

Determination of the local contrast of interference fringe patterns using continuous wavelet transform

Determination of the local contrast of interference fringe patterns using continuous wavelet transform Determination of the local contrat of interference fringe pattern uing continuou wavelet tranform Jong Kwang Hyok, Kim Chol Su Intitute of Optic, Department of Phyic, Kim Il Sung Univerity, Pyongyang,

More information

Alternate Dispersion Measures in Replicated Factorial Experiments

Alternate Dispersion Measures in Replicated Factorial Experiments Alternate Diperion Meaure in Replicated Factorial Experiment Neal A. Mackertich The Raytheon Company, Sudbury MA 02421 Jame C. Benneyan Northeatern Univerity, Boton MA 02115 Peter D. Krau The Raytheon

More information

Clustering Methods without Given Number of Clusters

Clustering Methods without Given Number of Clusters Clutering Method without Given Number of Cluter Peng Xu, Fei Liu Introduction A we now, mean method i a very effective algorithm of clutering. It mot powerful feature i the calability and implicity. However,

More information

Annex-A: RTTOV9 Cloud validation

Annex-A: RTTOV9 Cloud validation RTTOV-91 Science and Validation Plan Annex-A: RTTOV9 Cloud validation Author O Embury C J Merchant The Univerity of Edinburgh Intitute for Atmo. & Environ. Science Crew Building King Building Edinburgh

More information

Modelling a 2.5 MW direct driven wind turbine with permanent magnet generator

Modelling a 2.5 MW direct driven wind turbine with permanent magnet generator Modelling a.5 MW direct driven wind turbine with permanent magnet generator Report for the final examination project of the Nordic PhD coure on Wind Power by Thoma P. Fugleth Department of Electrical Power

More information

Linear Motion, Speed & Velocity

Linear Motion, Speed & Velocity Add Important Linear Motion, Speed & Velocity Page: 136 Linear Motion, Speed & Velocity NGSS Standard: N/A MA Curriculum Framework (006): 1.1, 1. AP Phyic 1 Learning Objective: 3.A.1.1, 3.A.1.3 Knowledge/Undertanding

More information

Physics 6A. Angular Momentum. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

Physics 6A. Angular Momentum. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB Phyic 6A Angular Momentum For Campu earning Angular Momentum Thi i the rotational equivalent of linear momentum. t quantifie the momentum of a rotating object, or ytem of object. f we imply tranlate the

More information

Bogoliubov Transformation in Classical Mechanics

Bogoliubov Transformation in Classical Mechanics Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How

More information

Stability regions in controller parameter space of DC motor speed control system with communication delays

Stability regions in controller parameter space of DC motor speed control system with communication delays Stability region in controller parameter pace of DC motor peed control ytem with communication delay Şahin Sönmez, Saffet Ayaun Department of Electrical and Electronic Engineering, Nigde Univerity, 5124,

More information

μ + = σ = D 4 σ = D 3 σ = σ = All units in parts (a) and (b) are in V. (1) x chart: Center = μ = 0.75 UCL =

μ + = σ = D 4 σ = D 3 σ = σ = All units in parts (a) and (b) are in V. (1) x chart: Center = μ = 0.75 UCL = Our online Tutor are available 4*7 to provide Help with Proce control ytem Homework/Aignment or a long term Graduate/Undergraduate Proce control ytem Project. Our Tutor being experienced and proficient

More information

Representation of a Group of Three-phase Induction Motors Using Per Unit Aggregation Model A.Kunakorn and T.Banyatnopparat

Representation of a Group of Three-phase Induction Motors Using Per Unit Aggregation Model A.Kunakorn and T.Banyatnopparat epreentation of a Group of Three-phae Induction Motor Uing Per Unit Aggregation Model A.Kunakorn and T.Banyatnopparat Abtract--Thi paper preent a per unit gregation model for repreenting a group of three-phae

More information

Root Locus Contents. Root locus, sketching algorithm. Root locus, examples. Root locus, proofs. Root locus, control examples

Root Locus Contents. Root locus, sketching algorithm. Root locus, examples. Root locus, proofs. Root locus, control examples Root Locu Content Root locu, ketching algorithm Root locu, example Root locu, proof Root locu, control example Root locu, influence of zero and pole Root locu, lead lag controller deign 9 Spring ME45 -

More information

Assessment of Performance for Single Loop Control Systems

Assessment of Performance for Single Loop Control Systems Aement of Performance for Single Loop Control Sytem Hiao-Ping Huang and Jyh-Cheng Jeng Department of Chemical Engineering National Taiwan Univerity Taipei 1617, Taiwan Abtract Aement of performance in

More information

Introduction to Laplace Transform Techniques in Circuit Analysis

Introduction to Laplace Transform Techniques in Circuit Analysis Unit 6 Introduction to Laplace Tranform Technique in Circuit Analyi In thi unit we conider the application of Laplace Tranform to circuit analyi. A relevant dicuion of the one-ided Laplace tranform i found

More information

Suggested Answers To Exercises. estimates variability in a sampling distribution of random means. About 68% of means fall

Suggested Answers To Exercises. estimates variability in a sampling distribution of random means. About 68% of means fall Beyond Significance Teting ( nd Edition), Rex B. Kline Suggeted Anwer To Exercie Chapter. The tatitic meaure variability among core at the cae level. In a normal ditribution, about 68% of the core fall

More information

LTV System Modelling

LTV System Modelling Helinki Univerit of Technolog S-72.333 Potgraduate Coure in Radiocommunication Fall 2000 LTV Stem Modelling Heikki Lorentz Sonera Entrum O heikki.lorentz@onera.fi Januar 23 rd 200 Content. Introduction

More information

LOW ORDER MIMO CONTROLLER DESIGN FOR AN ENGINE DISTURBANCE REJECTION PROBLEM. P.Dickinson, A.T.Shenton

LOW ORDER MIMO CONTROLLER DESIGN FOR AN ENGINE DISTURBANCE REJECTION PROBLEM. P.Dickinson, A.T.Shenton LOW ORDER MIMO CONTROLLER DESIGN FOR AN ENGINE DISTURBANCE REJECTION PROBLEM P.Dickinon, A.T.Shenton Department of Engineering, The Univerity of Liverpool, Liverpool L69 3GH, UK Abtract: Thi paper compare

More information

THE PARAMETERIZATION OF ALL TWO-DEGREES-OF-FREEDOM SEMISTRONGLY STABILIZING CONTROLLERS. Tatsuya Hoshikawa, Kou Yamada and Yuko Tatsumi

THE PARAMETERIZATION OF ALL TWO-DEGREES-OF-FREEDOM SEMISTRONGLY STABILIZING CONTROLLERS. Tatsuya Hoshikawa, Kou Yamada and Yuko Tatsumi International Journal of Innovative Computing, Information Control ICIC International c 206 ISSN 349-498 Volume 2, Number 2, April 206 pp. 357 370 THE PARAMETERIZATION OF ALL TWO-DEGREES-OF-FREEDOM SEMISTRONGLY

More information

Control Systems Engineering ( Chapter 7. Steady-State Errors ) Prof. Kwang-Chun Ho Tel: Fax:

Control Systems Engineering ( Chapter 7. Steady-State Errors ) Prof. Kwang-Chun Ho Tel: Fax: Control Sytem Engineering ( Chapter 7. Steady-State Error Prof. Kwang-Chun Ho kwangho@hanung.ac.kr Tel: 0-760-453 Fax:0-760-4435 Introduction In thi leon, you will learn the following : How to find the

More information

SERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lectures 41-48)

SERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lectures 41-48) Chapter 5 SERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lecture 41-48) 5.1 Introduction Power ytem hould enure good quality of electric power upply, which mean voltage and current waveform hould

More information

Trajectory Planning and Feedforward Design for High Performance Motion Systems

Trajectory Planning and Feedforward Design for High Performance Motion Systems Trajectory Planning and Feedforward Deign for High Performance Motion Sytem Paul Lambrecht, Matthij Boerlage, Maarten Steinbuch Faculty of Mechanical Engineering, Control Sytem Technology Group Eindhoven

More information

Preemptive scheduling on a small number of hierarchical machines

Preemptive scheduling on a small number of hierarchical machines Available online at www.ciencedirect.com Information and Computation 06 (008) 60 619 www.elevier.com/locate/ic Preemptive cheduling on a mall number of hierarchical machine György Dóa a, Leah Eptein b,

More information

Supplementary Figures

Supplementary Figures Supplementary Figure Supplementary Figure S1: Extraction of the SOF. The tandard deviation of meaured V xy at aturated tate (between 2.4 ka/m and 12 ka/m), V 2 d Vxy( H, j, hm ) Vxy( H, j, hm ) 2. The

More information

Longitudinal automatic control system for a light weight aircraft

Longitudinal automatic control system for a light weight aircraft Longitudinal automatic control ytem for a light eight aircraft Critian VIDAN*,, Silviu Ionut BADEA *Correponding author Military echnical Academy, Faculty of Mechatronic and Integrated Armament Sytem,

More information

Nonlinear Single-Particle Dynamics in High Energy Accelerators

Nonlinear Single-Particle Dynamics in High Energy Accelerators Nonlinear Single-Particle Dynamic in High Energy Accelerator Part 6: Canonical Perturbation Theory Nonlinear Single-Particle Dynamic in High Energy Accelerator Thi coure conit of eight lecture: 1. Introduction

More information

CHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS

CHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS Copyright 22 IFAC 5th Triennial World Congre, Barcelona, Spain CHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS Tritan Pérez Graham C. Goodwin Maria M. Serón Department of Electrical

More information

Adaptive Control of Level in Water Tank: Simulation Study

Adaptive Control of Level in Water Tank: Simulation Study Adaptive Control of Level in Water Tank: Simulation Study Jiri Vojteek and Petr Dotal Abtract An adaptive control i a popular, o called modern control method which could be ued for variou type of ytem

More information

Z a>2 s 1n = X L - m. X L = m + Z a>2 s 1n X L = The decision rule for this one-tail test is

Z a>2 s 1n = X L - m. X L = m + Z a>2 s 1n X L = The decision rule for this one-tail test is M09_BERE8380_12_OM_C09.QD 2/21/11 3:44 PM Page 1 9.6 The Power of a Tet 9.6 The Power of a Tet 1 Section 9.1 defined Type I and Type II error and their aociated rik. Recall that a repreent the probability

More information

A NEW LOAD MODEL OF THE PEDESTRIANS LATERAL ACTION

A NEW LOAD MODEL OF THE PEDESTRIANS LATERAL ACTION A NEW LOAD MODEL OF THE PEDESTRIANS LATERAL ACTION Fiammetta VENUTI PhD Politecnico di Torino Torino, IT Luca Bruno Aociate Profeor Politecnico di Torino Torino, IT Summary Thi paper propoe a new load

More information

UNIT 15 RELIABILITY EVALUATION OF k-out-of-n AND STANDBY SYSTEMS

UNIT 15 RELIABILITY EVALUATION OF k-out-of-n AND STANDBY SYSTEMS UNIT 1 RELIABILITY EVALUATION OF k-out-of-n AND STANDBY SYSTEMS Structure 1.1 Introduction Objective 1.2 Redundancy 1.3 Reliability of k-out-of-n Sytem 1.4 Reliability of Standby Sytem 1. Summary 1.6 Solution/Anwer

More information

HSC PHYSICS ONLINE KINEMATICS EXPERIMENT

HSC PHYSICS ONLINE KINEMATICS EXPERIMENT HSC PHYSICS ONLINE KINEMATICS EXPERIMENT RECTILINEAR MOTION WITH UNIFORM ACCELERATION Ball rolling down a ramp Aim To perform an experiment and do a detailed analyi of the numerical reult for the rectilinear

More information

MODELLING OF FRICTIONAL SOIL DAMPING IN FINITE ELEMENT ANALYSIS

MODELLING OF FRICTIONAL SOIL DAMPING IN FINITE ELEMENT ANALYSIS MODELLING OF FRICTIONAL SOIL DAMPING IN FINITE ELEMENT ANALYSIS S. VAN BAARS Department of Science, Technology and Communication, Univerity of Luxembourg, Luxembourg ABSTRACT: In oil dynamic, the oil i

More information

V = 4 3 πr3. d dt V = d ( 4 dv dt. = 4 3 π d dt r3 dv π 3r2 dv. dt = 4πr 2 dr

V = 4 3 πr3. d dt V = d ( 4 dv dt. = 4 3 π d dt r3 dv π 3r2 dv. dt = 4πr 2 dr 0.1 Related Rate In many phyical ituation we have a relationhip between multiple quantitie, and we know the rate at which one of the quantitie i changing. Oftentime we can ue thi relationhip a a convenient

More information

Digital Control System

Digital Control System Digital Control Sytem Summary # he -tranform play an important role in digital control and dicrete ignal proceing. he -tranform i defined a F () f(k) k () A. Example Conider the following equence: f(k)

More information

Stochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions

Stochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions Stochatic Optimization with Inequality Contraint Uing Simultaneou Perturbation and Penalty Function I-Jeng Wang* and Jame C. Spall** The John Hopkin Univerity Applied Phyic Laboratory 11100 John Hopkin

More information

Journal of Advanced Mechanical Design, Systems, and Manufacturing

Journal of Advanced Mechanical Design, Systems, and Manufacturing Vol., No. 4, 8 Mathematical Model of Linear Motor Stage with Non-Linear Friction Characteritic * Satohi KANEKO **, Ryuta SATO *** and Maaomi TSUTSUMI **** **Pot graduate, Tokyo Univerity of Agriculture

More information

Hylleraas wavefunction for He. dv 2. ,! r 2. )dv 1. in the trial function. A simple trial function that does include r 12. is ) f (r 12.

Hylleraas wavefunction for He. dv 2. ,! r 2. )dv 1. in the trial function. A simple trial function that does include r 12. is ) f (r 12. Hylleraa wavefunction for He The reaon why the Hartree method cannot reproduce the exact olution i due to the inability of the Hartree wave-function to account for electron correlation. We know that the

More information

Performance Improvement of Direct Torque Controlled Interior Permanent Magnet Synchronous Motor Drive by Considering Magnetic Saturation

Performance Improvement of Direct Torque Controlled Interior Permanent Magnet Synchronous Motor Drive by Considering Magnetic Saturation Performance Improvement of Direct Torque Controlled Interior Permanent Magnet Synchronou Motor Drive by Conidering Magnetic Saturation Behrooz Majidi * Jafar Milimonfared * Kaveh Malekian * *Amirkabir

More information

Molecular Dynamics Simulations of Nonequilibrium Effects Associated with Thermally Activated Exothermic Reactions

Molecular Dynamics Simulations of Nonequilibrium Effects Associated with Thermally Activated Exothermic Reactions Original Paper orma, 5, 9 7, Molecular Dynamic Simulation of Nonequilibrium Effect ociated with Thermally ctivated Exothermic Reaction Jerzy GORECKI and Joanna Natalia GORECK Intitute of Phyical Chemitry,

More information