Quaternion-Based Tracking Control law Design for Tracking Mode

Size: px
Start display at page:

Download "Quaternion-Based Tracking Control law Design for Tracking Mode"

Transcription

1 Quatenion-ased aking Contol law Design fo aking Mode A. M. Elbeltagy Aiaft Mehanis Dep. MC, Caio, Egypt; Y. Z. Elhalwagy Guidane and Navigation Dep. MC, Caio, Egypt; A. M. ayoumy Aiaft Mehanis Dep. MC, Caio, Egypt; +767 A. M. Youssef Aiaft Eleti system Dep. MC, Caio, Egypt; SSC6-V-8 ASRAC Vaious ontol design tehniues ae model dependent. hey typially euie knowledge of the inetia mati. hee ae majo hallenges fo eah poposed ontolle to ope with spaeaft mission objetive in tems of pointing and jitte euiements. hese hallenges inlude sensitivity to noise effets ando modeling eos, while othes ae sensitive to etenal toue distubanes, suh as toues indued by sola adiation pessue. Robust ontolles have been developed to mitigate these sensitivities. n this pape, a obust nonlinea taking ontol algoithm intodued peviously in the open liteatue is modified and toleated to be utilized with ehange momentum atuatos, e.g. eation wheels. he ontol law is using the ommanded attitude ate, ommanded attitude aeleation, attitude eo uatenion and gyosopi tems. aking eo dynamis euivalent to ellite losed-loop time-vaying nonlinea dynami system is used altenatively to onfim that a globally stable taking ontolle always eists. he poposed ontolle is applied to meet euiements of a taking omple mode. Geneation of the needed taget attitude and attitude ate ae deived in details. he motion and kinematis of the gound taget elative to the in-obit ellite is analyzed and desibed in obitefeened oodinates. he ellite dynamis ae deived fom fist piniples and efomulated also in obit efeened oodinates.. A taking sheme fo the pointing ais along the body z-ais of ellite is highlighted. Consideing attitude and obit ontol system (AOCS with ideal attitude and obit detemination sensos with symmeti ellite inetia, the validity of poposed ontolle and taget data geneato is demonstated unde MALASMULNK envionment. MALA optimization tool is used fo optimal gains seletion. Robustness of the globally stable modified ontol law to spaeaft inetia mati unetainty is also disussed. Simulation esults show that the poposed ontol law an be used suessfully onboad fo fast taking and is obust enough to keep the pointing auay within aeptable limits with onsideable inetia unetainty. Elbeltagy 3 th Annual AAAUSU Confeene on Small Satellites

2 NRODUCON Fo seveal spae missions, Eath-pointing ellites ae euied to point a payload suh as a high-gain antenna, amea, and telesope to tak a fied taget on the Eath fo a etain peiod of time in ode to povide impoved long-peiod up-down-link ommuniation, low-distotion imaging and auate obsevation. hese missions euie gound taget taking ontol. -3 A linea uatenion feedbak egulato with open loop deoupling ontol toue fo gyosopi foes to ensue inetially efeened eigenais otations is poposed in Ref. 3. A taking ontolle whih deals with muh moe ompliated ontol task onsides the taget desied ate duing ealizing the desied taget attitude. A taking ontolle is ealized though a feedbak ontol sheme ontaining an inne veloity loop that taks the desied ate ommand and oute attitude loop that taks the desied attitude ommand to ahieve taking via an instantaneous eigenais otation. Goeee and Shuke developed a taking ontolle with a feedbak elasti tem fo taking attitude, a feedbak visous tem fo taking a desied angula ate, and a feed-fowad model-based ompenion tem. he aim of thei algoithm is to enable the telesope mounted along the z-ais of ellite body to tak a gound station while passing ove and to eeive the lase ommuniation signals fom the gound station. A new non-linea taking ontol algoithm based on an attitude eo uatenion uses the ommanded attitude ate without tansfomation into the body fame is also intodued. 4 he diet use of the ommanded attitude ate simplifies the alulation of its deivative, whih is used in the ontol law. his algoithm is poposed to be used with etenal atuatos fo geneal taking task. Petubation dynamis with seven state vaiables is used to analyze the stability of the system and the taking pefomane. Simulation esults show that the spaeaft an tak the ommanded attitude and ate uikly fo a non-zeo aeleation ate ommand. his poposed taking ontolle was the motivation to popose a modified non-linea ontolle fo fied gound taget taking puposes. Most of the ontol algoithms intodued peviously in the open liteatue ae about pointing auay and ignoe the time elapsed to eah these stit euiements. Atually thee is an ineasing demand fom ustomes fo a fast gound-taget taking even on the epense of pointing auay itself. When the taget loations and ellite gound station ae both within the ellite s footpint, it may be muh moe impotant fo the gound station to stee the main image boe sight towads the immediate euied taget aeas befoe aptuing the imaging data. Quik ellite esponse allows fo the militay intelligene data gatheing and downloading in the same ommuniation session. n this pape, a non-linea taking ontol algoithm intodued in Ref. 4 is modified and toleated to be utilized with ehange momentum atuatos, e.g. eation wheels, to tak any desied taget on the Eath fo a nominally Eath- pointing ellite as fast as possible. he poposed taking ontolle uses attitude eo epesented in uatenion, ate eo without the need to be epesented in ellite body oodinate fame, taget aeleation, gyosopi tems. he pape is oganized as following: Setion define the oodinates fo attitude tansfomation though the pape. he ellite non-linea ombined model euipped with thee-ais eation wheels (RW is pesented in setion 3. he RW non-linea attitude taking ontolle fo fied gound taget is pesented in setion 4 followed by stability analysis using taking eo dynamis in setion 5. he kinemati of fied gound taget taking inluding the taget desied ate and aeleation is totally oveed in details in setion 6. A design eample with its simulation esults using MALA Simulink softwae is pesented in setion 7 and finally a onlusion is given in setion in setion 8. DEFNON OF COORDNAE FRAMES Seveal oodinate fames ae used oodinates fo attitude tansfomation though the pape as follows Obit Refeened Coodinate fame ( he z-ais points towads the ente of the Eath stating fom the ellite in obit. he y-ais points to the negative obit nomal. he -ais is hosen to fom a ight-handed othogonal efeene fame. heefoe, fo a iula obit, the -ais will be along the veloity veto of the ellite. Spaeaft ody Fied Coodinate fame (SC he SC oodinate fame oiginates fom the ente of mass of ellite. he -ais points to the ente of the haness side of the ellite. he z-ais is othogonal to the -ais. he y-ais is hosen to fom a ight-handed othogonal efeene system. When the ellite is pefet nadi pointing without any otation aound the z-ais, the SC and oodinate fames ae assumed to be aligned. SC oodinates oiginate fom the ente of mass of ellite. he -ais points to the ente of the haness side of the ellite. he z-ais is othogonal to the -ais. he y-ais is hosen to fom a ight-handed othogonal efeene system. When the ellite is pefet nadi pointing without any otation aound the z-ais, the SC and oodinates ae assumed to be aligned. Elbeltagy 3 th Annual AAAUSU Confeene on Small Satellites

3 Eath Fied Coodinate fame (EFC he EFC fame is Eath-enteed. he -ais will point towads the pime meidian passing though the ente of the Royal Geenwih Obsevatoy in London. he z ais points to the noth elestial pole. he y-ais is hosen to fom a ight-handed othogonal efeene fame. Eath Centeed netial Coodinates (EC he EC fame is Eath-enteed. he -ais points towads the mean euino. he z-ais points to the elestial pole. he y-ais is hosen to fom a ight handed othogonal efeene fame. SPACECRAF COMNED NONLNEAR MODEL Assuming a igid ellite euipped with thee ais eation wheels (RW as intenal toue atuatos and by ealling the Eule s moment euation, then the dynami model of an Eath-pointing ellite is given by 5 (. J ω ( = N et ω Jω + h h ( v v ωo = [ ] (8 ( 4 ( ω = J ( N Ω ( ω ( Jω + h h et Whee J is the spaeaft inetia mati, ω = ( ω, ω, ω3 is the body angula veloity veto with espet to the inetial efeene fame, N et is the etenal distubane toue, h = ( h, h, h3 is the RW angula momentum veto, h is RW applied toue and Ω ( ω is a skew-symmeti mati defined by ω ω Ω ( ω ω ω 3 = 3 ω ω (3 Futhe assuming the ellite is 3-ais stabilized, then the absolute angula veloity in the inetial spae esolved in the fame is given by 5 ω = ω + A ω (4 SC O Whee ω O = obital angula ate veto of the ellite motion in fame. Substitute E. (4 into E. (, then the dynami euation of the Eath-pointing ellite beomes ω = J ( N + N Ω ( ω ( Jω + h h ωo et With SC SC N ωo = JA ωo JA ωo SC SC = JΩ( ω A ω JA ω Whee ω and O O O Whee N eff = effetive RW toue input. he nonlinea attitude kinematis euations of motion of an Eath-pointing ellite an be epesented by using vaious attitude paametes. Repesentation though uatenion paamete has the popety of (. nonsingulaity and it is fee fom the tigonometi omponent. heefoe, this epesentation is widely used to study the attitude behaviou of spaeaft. he kinematis of the ellite model is the pat whih epesses the elation between the attitude and angula veloities of the body and an be desibed by 6 (5 (6 ω O ae obit angula ate and obit angula aeleation epesented in fame and deived as follows v ωo = [ ] (7 An euivalent but moe useful fom of E. (5 fo ontol puposes is given by ω = J ( N ωo + N et Ω ( ω ( Jω + N eff h = Ω( ω ( h N eff = Ω ( ω + 4ω 4 ( = ω (9 ( Elbeltagy 3 3 th Annual AAAUSU Confeene on Small Satellites

4 he ellite attitude w..t fame as a efeene fame is detemined by its uatenion. oth and 4 of the uatenion ae defined as ϕ e sin( ϕ = esin( = 3 ϕ e 3 sin( ( ϕ 4 = os ( ( = 4 (3 Whee e, e and e 3 ae the omponents of the otation ais unit veto along the efeene fame; ϕ is the otation angle. he ombined dynami, E. (5 and kinemati, E. ( give the geneal nonlinea model fo the spaeaft angula motion with ten state vaiables d J ( N ωo + N et Ω ( ω ( Jω + N eff ω Ω ( ω + 4ω = 4 ( ω h ( ω ( h N Ω eff (4 RW NONLNEAR AUDE RACKNG CONROLLER ased on the non-linea taking ontol law intodued in Ref. 4 whih use attitude eo uatenion and desied ate without tansfomation in the body fame, a modified poposed RW nonlinea taking ontolle fo gound taget taking puposes is witten as follows h = N ω + N Ω ( ω ( J ω + h o et Ω ( ω J ω +Ω( ω J ω Ω( Jω ( ω ω J ω + D ( ω ω + K δ (5 Whee δ = attitude eo uatenion veto of δ attitude eo uatenion, whih is defined as follows δ δ (. δ 4 ( δ (7 Whee = the invese of ; (. = uatenion multipliation. Knowing thatδ is still a uatenion whih has a physial meaning as it epesents the attitude diffeene between and. n suh ase effetive RW toue input is N eff = N ωo N et +Ω( ω J ω + J ω D ( ω ω K δ +Ω( ω J ω +Ω( Jω ( ω ω Ω( ω J ω Substitute E. (5 into E. (5 gives ω ω = J ( D ( ω ω + K δ Ω( ω J ω +Ω( Jω ( ω ω (8 +Ω( ω J ω E. (9 and E. ( fom the losed-loop time-vaying nonlinea dynami system of the spaeaft attitude with the ontol law (5. SALY ANALYSS (9 Closed-loop solution t is lea that the following is a solution of euation (9 and euation ( O ω = ω = ( Elbeltagy 4 3 th Annual AAAUSU Confeene on Small Satellites

5 Solution ( means theoetially that it is possible fo the spaeaft to follow the ommanded attitude and attitude ate without estition although in eality thee ae always etain limits in using diffeent atuatos. he tems N ω, Ω ( ω ( J ω + h and o J ω ae essential to be inluded in the poposed RW ontolle dediated fo taking task of nominal eath pointing ellite. nluding these tems ensue that the solution ( always eist and hene an be applied in diffeent maneuve onditions (e.g. est to est maneuve, taking ommanded attitude with onstant ate. he stability analysis an daw a lea onlusion onening the ability to tak the ommanded attitude when solution ( eists. he ommanded ate and its deivative ω and ω espetively will estimated late fo fied gound taget with aid of GPS data. heefoe, the nonlinea ontolle is appliable fo spaeaft attitude ontol and ommand taking. aking eo dynamis Geneally asymptoti stability analysis is euied to ensue that the ellite an follow the ommanded attitude with ommanded attitude ate fom any ondition of both ω and. he stability analysis is stated by studying the loal asymptoti stability, at whih the initial values of both ω and ae nea the solution (. One this loal stability is onfimed, a global asymptoti stability, at whih any values of both ω and ae assumed, is futhe seahed. One of the objetives of this pape is to find the onstant gain maties D and K having a globally asymptoti stable solution ( and guaanteeing the needed solution onvege pefomane. hese gains ae used as initial design values to the Matlab optimization tools. aking eo may be used as paamete to be minimized in limited time fame.he taking eo dynamis intodued in Ref. 4 ae employed fo the indiet stability analysis altenatively to the nononvenient diet analysis of solution (. he taking eo dynamis defined in the following ae always valid whateve the diffeenes between the ommanded ates and the atual ates ae lage o small. We define the following state vaiable to epesent the taking eo dynamis as follows Whee ( ω 3 ω = ω ω = ω = ω 3 = = = = ( ( = ( ( 3 3 (3 46= With (4 Substituting (-4 into E. (9 and E. (, getting simila fomula dedued in Ref. 4 fo the taking eo dynamis with mino hanges as 46 3 = f3 = D3 K[ ] 7 46 = f 46 = Q Ω ( ω 46 + ω 7 X ( = = 3 + ω f (5 Elbeltagy 5 3 th Annual AAAUSU Confeene on Small Satellites

6 Whee D = J D (6 K = J K (7 [ ] Q 3 3 X 3 3 (8 (9 (3 he taking eo dynamis E. (5 is euivalent to E. (9 and E. (. E. (5 has an euilibium point =, whih oesponds to solution (. heefoe, instead of studying the deviation of ω and fom ω and fo E. (9 and E. (, we may simply study the stability of the taking eo E. (5 with espet to the euilibium point =. Global stability analysis A Lyapunov funtion based global stability analysis of the time-vaying non-linea taking eo dynamis E. (5 is employed. he states vaiable in this eo dynamis epesents the deviation fom obit efeened ellite body ate, uatenion and desied taget ate and uatenion espetively. Although the states vaiable of the nonlinea dynami system E.(9 and E. ( with the poposed ontol law (5 is diffeent fom those used in Ref. 5, similaity in the stutue of taking eo dynamis is enough to use theoem onditions n Ref. 4 to guaantee stability. t is poven that when K is symmeti and positive definite and that K D is positive definite, the euilibium point = of taking eo dynamis is globally stable. he andidate Lyapunov funtion V, whih is independent of time and is adially unbound, is defined by V( Whee P P = P ( = (3 P is a symmeti and positive definite mati. 4 4is an identity mati of ode 4. he total time deivative of V along the tajetoies of the taking eo dynamis is poved to be negative semidefinite when P is seleted P = K (33 hen (34 V( = PD 3 3 heefoe the euilibium point = of E. (5 is globally stable. Futhemoe using theoem 8.4 of Ref. appoahes zeo as t. 7, we onlude that 3 Assuming D and K ae diagonal maties epessed as D = diag ( d, d, d 3 K = diag ( k, k, k 3 (35 he dynami behaviou of the nonlinea system E.(5 an be ontolled aound the euilibium point in the dietion of eigenvalues with nonzeo eal pats of the Jaobian mati. 4 he haateistis euation of the Jaobian mati fo the speial ase when ω = an be epanded as 4 k k λ( λ + d λ + ( λ + dλ + k 3 ( λ + d 3λ + = (36 Elbeltagy 6 3 th Annual AAAUSU Confeene on Small Satellites

7 Epessing E. (36 by the desied damping atio ξ and natual feuenyω, then the gains maties ae n elated dietly to the desied damping atio ξ and natual feuenyω as 4 n d= ξω n; k= ωn d= ξω n; k = ωn d3= ξω 3 n3; k3= ωn3 (37 Sine atuato onstaints ae onsideed in this pape, these gain maties figued using E. (37 an be used as initial values guide. A MatlabSimulink Optimization tool is used to adjust the dynami behaviou of the nonlinea system E. (5 aound the euilibium point. aking eo is the seleted paamete to be optimized within shot time. KNEMAC OF FXED GROUND ARGE RACKNG aget aking Mode will most likely be ativated fom Nominal Mode as soon as a pedefined gound taget loation (e.g. Gound station is in ange. Knowledge of the ellite position is theefoe euied by the on-boad ompute. he method that will be used to alulate the desied attitude and angula ates fo taking a taget was pesented by Chen et al. in. 8 he ellite will need to be able to detemine whethe the taget is in ange (i.e. in the ellite s field of view of the Eath. he distane D ma fom the ellite to the futhest point on Eath that the ellite an "see" an be alulated as D = ( R + h R (38 ma E E Whee R E = adius of the eath and altitude. is ellite An illustation of the above-mentioned geomety an be seen in Fig.. Compaing the uent distane to the taget with maimum taget taking distane allow the on-boad ontol algoithms to know if the taget is within the visibility zone o not and hene taking task an be stated. Figue : maimum taget taking distane Sine the nominal mode of eath obiting ellite usually maintains the ellite body fame aligned with obital ( fame, it is pefeable to desibe the kinematis of the gound taget with espet to the inobit ellite in fame. he veto fom the ellite to the gound taget epesented in will speify the taking dietion of the seleted pointing ais whih may be the positive mounting ais of a amea o antenna. Fistly the veto fom the ente of the Eath to the gound taget with espet to EFC fame fom the given geoenti latitude (assuming a spheial Eath and longitude will be deived as EFC = R E ( λ os( ϕ ( λ os( ϕ os sin sin( ϕ (39 Using the tansfomation mati fom EFC fame to EC fame,, the loation of the taget with espet to the EC fame is given by EC EC EFC AEFC = (4 Assuming that the Eath has a onstant angula ate aound its otation ais duing the taking maneuve, then the attitude tansfomation mati fom the EFC fame to EC fame an be epesented by os( ωet + α sin( ωet + α EC AEFC = sin( ωet α os( ωet α + + (4 Whee t is time sala and is the initial phase between the -aes of both EFC fame and EC fame. EC and (the ellite s position veto in EC fame, as obtained by onboad GPS eeive o Elbeltagy 7 3 th Annual AAAUSU Confeene on Small Satellites

8 podued by any obit popagato fo simulation puposes an then be used to alulate S, the veto fom the ellite to the taget epesented in obit fame, thus EC S = AEC ( (4 A EC As, is the tansfomation mati fom EC fame to fame and defined by A EC Whee wˆ [ uvw ˆ ˆ ˆ ] = (43 = (44 v ˆ = v v (45 uˆ = vˆ wˆ (46 f the magnitude of S is lage than ma, then the taget is not in ange and the ellite will ontinue opeating in nominal mode. Howeve, if < D, taget taking mode will be S ma ativated. atually the dietion of S is of use to the ellite s ADCS, hene the unit veto u S must be alulated by nomalised, thus u S S S S = (47 aget taking task euies that the mounting ais of the ommissioned payload being ontolled to point towads the dietion of u S. n most ases, the ommissioned payload, suh as a amea, u SC om antenna o telesope is usually mounted along the z- ais of the ellite body. his onfiguation theoetially allows the ADCS designes to ahieve the task with minimum esoues (i.e. task an be ahieved by only using the - and y- eation wheels. Nevetheless this onfiguation pemits speifially ontolling the angle aound the mounting ais. D SC Moeove, defining A _ d to epesent the desied attitude fom fame to SC fame, then A u u SC SC _ d S om = (48 Let the --3 seuene of otations is used fo the SC desiption of the desied attitude mati A _ d fo taking. Substituting suh attitude mati 9 in tems of φθ, and ψ espetively epesent the oll, pith and yaw angles in E. (48 leads to us = [sin θ osθsin φ osθosφ] (49 ased upon E. 49 the oll and pith angles φ and fo taking an be easily solved. Assuming a onstant yaw ψ is euied duing taking, the desied attitude SC mati A _ d an be omputed. heefoe, the efeened uatenion ommand an be deived. Aodingly, δ an be solved using E. (7. he ellite s angula ate fo taget taking will not be zeo, sine the taget is fied with espet to the Eath and will thus be otating elative to the obiting ellite. he ommanded angula ate ω an be alulated as epesented in 8 ω = u u (5 S S Due to the assumption of a onstant yaw angle duing taking, the desied angula ate along the z-ais will be zeo. ased upon E. (5, ω an be edued to 8 S _ y S _ us _ z us _ z ω = u u [ ] (5 he ommanded angula ate ω an be deived by taking the time deivative of E. (5 ω = u u + u u S S S S = us u S (5 he u S tem an be deived by fist taking the time deivative of E. (47 yields Elbeltagy 8 3 th Annual AAAUSU Confeene on Small Satellites

9 u = [ u ( u ] (53 S 3 S S S S u So S an be deived by taking the time deivative of E. 53 yields d EC EC S = ( A EC ( AEFC d EC EC + ( A EC ( A EFC v (58 u d( S S S = S S S d u ( u S S S ( S (54 As d ( A EC EC EC ( A EFC = A ( A + A ( A EC EC EC EC EC EFC EC EFC A EC EFC EC (59 ut geneally fo any veto Α d ( Α = ( Α Α Α (55 Using E. (55 and substitute the oesponding elations in E. (54 and make some manipulation yields u S S = ( ( 3 S S S S + ( ( u S 3 S S S S ( u ( u ( u S S S S S + u ( u + u S S S S S S (56 Futhemoe fom E. (4, S an be detemined as = A ( A + A EC EC ( A v EC EC S EC EFC EC EFC (57 d ( A EC EC EC ( A EFC v = EC EC EC EC A EC ( AEFC v + AEC ( A EFC (6 a Whee A EC [ ˆ ˆ ˆ = uvw ] (3 ˆ w = ( ˆˆ 3 ww v (4 v ˆ = (5 uˆ = vˆ [( ˆˆ 3 ww v ] (6 sin( ωet + α os( ωet + α EC A EFC = ω E os( ωet α -sin( ωet α + + (7 [ ˆ ˆ ˆ A = uvw ] (8 EC Diffeentiating E and using E. 55 to get the olumn vetos of A EC as follows aking the deivative of E. (57 get Elbeltagy 9 3 th Annual AAAUSU Confeene on Small Satellites

10 ˆ a w = + ( v v ˆ ˆ ˆ ww wv ˆ + w a + v + ( ( ˆ ˆ v w v w (67 v fom a GPS eeive. Matlab Simulink Optimization tools is used to seah almost best gains to allow the taking poess meets the euied apid pointing to the taget with a little elaed auay. his yields the following gains v ˆ = (68 ˆ u = ˆ ( ˆ 3 v v v + v a + ( ˆ ˆˆ 3 v v ww v ˆ ( vˆ ( ww ˆ v + wˆˆ w a (69 Simulation is made fo one of the ovepass flight. Fig. 4-5 shows that the pointing eo of the poposed taking ontolle of E. (5 sueeded to maneuve ellite though a lage angle to a pedetemined attitude with the euied elaed auay (.8 within only (4 se whih is well below the euiements. Fig. 6-7 shows the uatenion eo and ate eo simulation esults. t is seen that both obit efeened ellite attitude and desied attitude ae almost oinide within only (3 se. Whee a is ellite aeleation veto as 5 µ asa t = = (7 3 So finally using E to get u S and substitute in E. 5 to obtain ω A MatlabSimulink softwae ode has been built to simulate the ellite fied gound taget taking task. he oveall system simulato inludes many subsystems o modules. DESGN EXAMPLE n this setion, a typial design eample has been used in simulation tests to veify the pefomane of the modified poposed taking ontolle pesented above. An imaginay ellite in a low-eath-eenti obit is used as an eample duing these simulations. n ode to investigate the attitude hange of the ellite, we define the oll, pith and yaw angles espetively to epesent the otations of the ellite body, y and z aes with espet to the oodinates. he simulation paametes ae given in able. Duing simulation, we assume pefet attitude knowledge and. We also assume pefet measuements of the vetos and Reuied aking auay aget position Unit veto of ommissioned ais in SC fame Moment of able : Simulation paametes (.3 within ( min of being in ellite FOV EFC [ ] = u = SC om [ ] inetia Reation wheel Obit paametes nitial phase α = kgm Ma toue = (. Nm, Ma momentum = (4 Nms, nitial Momentum = [ ] Nms peigee altitude = (65 km ; inlination (64.5 ; eentiity (.3, agument of peigee ( ; ight asension of asending node ( ; initial mean anomaly ( Elbeltagy 3 th Annual AAAUSU Confeene on Small Satellites

11 he obit efeened angula veloity ω of the ellite body does tak the desied angula ate ommand ω with aeptable limits duing the taking peiod at whih the taget is within the FOV of the ellite. he ativities of the thee-ais eation wheels epesented in wheel ealized toue and wheel momenta ae shown in Figs Duing the taking peiod, the toues of thee ais eation wheels is fa fom uation. he eation wheels will apply maimum toue (. Nm. to povide a bang-bang ontol manoeuve to meet the euied taget taking within speified limited time with elaed pointing auay. he wheel momentums ae kept well away below the pemissible limits (4 Nms. Although the gavity gadient distubane is onsideed in simulation, this safe patten an guaantee allowable gadually momentum build up aused by othe seula etenal distubane toues suh as aeodynami and sola pessue foes whih ae not onsideed in simulation in this pape. Figue 4: Quatenion Eo Figue 5: Rate Diffeene Eo Figue : Gound aget aking Eo Figue 6: Wheel oue Figue 3: Gound aget aking Eo (zoom in Elbeltagy 3 th Annual AAAUSU Confeene on Small Satellites

12 Figue 7: Wheel Momenta Due to the feedbak natue of poposed ontolle, this ontolle is obust against model unetainties. Simulations wee also done to investigate the obust behaviou of the taking ontolle against a ±% eo in the moment of inetia tenso of the ellite. hese simulations show in Fig. that the taking eo an still be maintained within euied limits duing the fast gound taget taking. Figue 9: Gound aget aking Eo Figue 8: Gound aget aking Eo with up to ±% inetia unetainty (Zoom in he gain seletion of the poposed RW non-linea ontolle an also be seleted to guaantee fo muh bette taking auay. Matlab Simulink Optimization tools is used again to seah fo almost best gains to allow the taking poess meets bette taking auay. his yields the following gain Figue : Gound aget aking Eo (Zoom in he simulation esults Figs - shows that the poposed RW non-linea ontolle is not limited to fast taking with elaed pointing auay. he pointing auay an eah (. and futhe enhanement an be aied out. he enhanement in taking auay will be on the epense of the time elapsed to tak the fied gound taget and atuato effot. CONCLUSONS n this pape a non-linea taking ontol algoithm intodued peviously in the open liteatue is modified and toleated to be utilized with ehange momentum atuatos, e.g. eation wheels, fo speifially fied gound taget taking task. he poposed ontolle aes about fast gound-taget taking. Quik ellite esponse allows fo the militay intelligene data gatheing and downloading in the same ommuniation session. he poposed RW non-linea taking ontolle is using the ommanded attitude ate, ommanded attitude aeleation, attitude eo uatenion and gyosopi tems. A systemati method fo detemine a omplete kinematis of fied gound taget elative to the obit fame has been pesented. he simulation esults show that RW non-linea Elbeltagy 3 th Annual AAAUSU Confeene on Small Satellites

13 ontolle is apable to ahieve a fast gound taget taking maneuve one the taget is being in the ellite FOV. he eation wheels will apply maimum toue to povide bang-bang ontol maneuves to meet the euied taget taking within speified limited time with elaed pointing auay. he wheel momentum is kept well away below the pemissible limits. he simulation esults also show well and onsideable obustness behavio of the poposed ontolle.he enhanement in taking auay is possible but it will be on the epense of the time elapsed to tak the fied gound taget and atuato effot. Refeenes. Weiss, H., "Quatenion-ased RateAttitude aking System with Appliation to Gimbal Attitude," Jounal of Guidane, Contol and Dynamis, Vol. 6, No.4, July-August 993, pp Goeee,., and Shuke,., "Geometi Attitude Contol of a Small Satellite fo Gound aking Maneuves," Poeedings of the AAAUSU Confeenes on Small Satellites, Utah State Univesity, Sept Wie,., and aba, P., "Quatenion Feedbak fo Spaeaft Lage Angle Maneuves," Jounal of Guidane, Contol, and Dynamis, Vol. 8, No. 3, May-June 985, pp Z. Zhou and R. COLGREN, "NONLNEAR AUDE CONROL FOR LARGE AND FAS MANEUVERS," AAA Guidane, Navigation, and Contol Confeene and Ehibit, 5-8 August 5 5. MARCEL J. SD, " Spaeaft Dynamis and Contol". 6. F. Landis Makley, Jhon L. Cassidis, "Fundamentals of Spaeaft Attitude Detemination and Contol", 4s 7. Hassan Khalil, "Nonlinea Systems," (3ed, Xiaojiang Chen, Willem H. Steyn, Yoshi Hashida"Gound-aget aking Contol Of Eath-Pointing Satellites," AAA Guidane, Navigation, and Contol Confeene and Ehibit,. 9. A. M. Elbeltagy, A. M. ayoumy Ali, A. M. Youssef, Y. Z. Elhalwagy, "Modeling And Simulation of Spaeaft Pointing Modes Using Quatenion-ased Nonlinea Contol Laws," AAASieh Confeene and Ehibit, 5. Elbeltagy 3 3 th Annual AAAUSU Confeene on Small Satellites

From E.G. Haug Escape Velocity To the Golden Ratio at the Black Hole. Branko Zivlak, Novi Sad, May 2018

From E.G. Haug Escape Velocity To the Golden Ratio at the Black Hole. Branko Zivlak, Novi Sad, May 2018 Fom E.G. Haug Esape eloity To the Golden Ratio at the Blak Hole Banko Zivlak, bzivlak@gmail.om Novi Sad, May 018 Abstat Esape veloity fom the E.G. Haug has been heked. It is ompaed with obital veloity

More information

Discrete-Time Immersion and Invariance Adaptive Control of a Slider-crank Mechanism

Discrete-Time Immersion and Invariance Adaptive Control of a Slider-crank Mechanism Pepints of the 9th Wold Congess he Intenational Fedeation of Automati Contol Cape own South Afia August -9 0 Disete-ime Immesion and Invaiane Adaptive Contol of a Slide-an Mehanism Yapa Yalçın Dept of

More information

(conservation of momentum)

(conservation of momentum) Dynamis of Binay Collisions Assumptions fo elasti ollisions: a) Eletially neutal moleules fo whih the foe between moleules depends only on the distane between thei entes. b) No intehange between tanslational

More information

Answers to Coursebook questions Chapter 2.11

Answers to Coursebook questions Chapter 2.11 Answes to Couseook questions Chapte 11 1 he net foe on the satellite is F = G Mm and this plays the ole of the entipetal foe on the satellite, ie mv mv Equating the two gives π Fo iula motion we have that

More information

Extra Examples for Chapter 1

Extra Examples for Chapter 1 Exta Examples fo Chapte 1 Example 1: Conenti ylinde visomete is a devie used to measue the visosity of liquids. A liquid of unknown visosity is filling the small gap between two onenti ylindes, one is

More information

OBSTACLE DETECTION USING RING BEAM SYSTEM

OBSTACLE DETECTION USING RING BEAM SYSTEM OBSTACLE DETECTION USING RING BEAM SYSTEM M. Hiaki, K. Takamasu and S. Ozono Depatment of Peision Engineeing, The Univesity of Tokyo 7-3-1 Hongo, Bunkyo-ku, Tokyo, Japan Abstat: In this pape, we popose

More information

Vision Sensor. Vision. (Phase 1) pre-shaping. Actuator. Tactile Sensor. Vision. (Phase 2) shaping. Actuator. Tactile Sensor.

Vision Sensor. Vision. (Phase 1) pre-shaping. Actuator. Tactile Sensor. Vision. (Phase 2) shaping. Actuator. Tactile Sensor. Optimal Gasping using Visual and Tatile Feedbak Akio NAMIKI Masatoshi ISHIKAWA Depatment of Mathematial Engineeing and Infomation Physis Univesity of Tokyo Tokyo 3, Japan namik@k.t.u-tokyo.a.jp Abstat

More information

Rigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018

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

Time Dilation in Gravity Wells

Time Dilation in Gravity Wells Time Dilation in Gavity Wells By Rihad R. Shiffman Digital Gaphis Asso. 038 Dunkik Ave. L.A., Ca. 9005 s@isi.edu This doument disusses the geneal elativisti effet of time dilation aused by a spheially

More information

1 Fundamental Solutions to the Wave Equation

1 Fundamental Solutions to the Wave Equation 1 Fundamental Solutions to the Wave Equation Physial insight in the sound geneation mehanism an be gained by onsideing simple analytial solutions to the wave equation One example is to onside aousti adiation

More information

DARK MATTER AND THE DYNAMICS OF GALAXIES: A NEWTONIAN APPROACH 1. INTRODUCTION

DARK MATTER AND THE DYNAMICS OF GALAXIES: A NEWTONIAN APPROACH 1. INTRODUCTION DARK MATTER AND THE DYNAMICS OF GALAXIES: A NEWTONIAN APPROACH Mugu B. RĂUŢ Coesponding autho: Mugu RĂUŢ, E-mail: m_b_aut@yahoo.om Abstat In this pape I popose a oetion to the well-known Newtonian gavitational

More information

8.022 (E&M) Lecture 13. What we learned about magnetism so far

8.022 (E&M) Lecture 13. What we learned about magnetism so far 8.0 (E&M) Letue 13 Topis: B s ole in Mawell s equations Veto potential Biot-Savat law and its appliations What we leaned about magnetism so fa Magneti Field B Epeiments: uents in s geneate foes on hages

More information

1 Fundamental Solutions to the Wave Equation

1 Fundamental Solutions to the Wave Equation 1 Fundamental Solutions to the Wave Equation Physial insight in the sound geneation mehanism an be gained by onsideing simple analytial solutions to the wave equation. One example is to onside aousti adiation

More information

Numerical Modeling in Biomedical Systems

Numerical Modeling in Biomedical Systems Numeial Modeling in Biomedial Systems BME 15:35 Letue 7 9/6/17 Nonlinea Systems Dunn Chapte 5 Nonlinea equations Root finding Baketing methods Open methods Gaphial Bisetion False Position Newton s method

More information

Red Shift and Blue Shift: A realistic approach

Red Shift and Blue Shift: A realistic approach Red Shift and Blue Shift: A ealisti appoah Benhad Rothenstein Politehnia Uniesity of Timisoaa, Physis Dept., Timisoaa, Romania E-mail: benhad_othenstein@yahoo.om Coina Nafonita Politehnia Uniesity of Timisoaa,

More information

PHYS 110B - HW #7 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased

PHYS 110B - HW #7 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased PHYS B - HW #7 Fall 5, Solutions by David Pae Equations efeened as Eq. # ae fom Giffiths Poblem statements ae paaphased [.] Poblem.4 fom Giffiths Show that Eq..4, V, t an be witten as Eq..44, V, t q t

More information

Photographing a time interval

Photographing a time interval Potogaping a time inteval Benad Rotenstein and Ioan Damian Politennia Univesity of imisoaa Depatment of Pysis imisoaa Romania benad_otenstein@yaoo.om ijdamian@yaoo.om Abstat A metod of measuing time intevals

More information

Physics 218, Spring March 2004

Physics 218, Spring March 2004 Today in Physis 8: eleti dipole adiation II The fa field Veto potential fo an osillating eleti dipole Radiated fields and intensity fo an osillating eleti dipole Total satteing oss setion of a dieleti

More information

In electrostatics, the electric field E and its sources (charges) are related by Gauss s law: Surface

In electrostatics, the electric field E and its sources (charges) are related by Gauss s law: Surface Ampee s law n eletostatis, the eleti field E and its soues (hages) ae elated by Gauss s law: EdA i 4πQenl Sufae Why useful? When symmety applies, E an be easily omputed Similaly, in magnetism the magneti

More information

AVS fiziks. Institute for NET/JRF, GATE, IIT-JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES

AVS fiziks. Institute for NET/JRF, GATE, IIT-JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES ELECTROMAGNETIC THEORY SOLUTIONS GATE- Q. An insulating sphee of adius a aies a hage density a os ; a. The leading ode tem fo the eleti field at a distane d, fa away fom the hage distibution, is popotional

More information

Special Relativity in Acoustic and Electromagnetic Waves Without Phase Invariance and Lorentz Transformations 1. Introduction n k.

Special Relativity in Acoustic and Electromagnetic Waves Without Phase Invariance and Lorentz Transformations 1. Introduction n k. Speial Relativit in Aousti and Eletomagneti Waves Without Phase Invaiane and Loentz Tansfomations Benhad Rothenstein bothenstein@gmail.om Abstat. Tansfomation equations fo the phsial quantities intodued

More information

Mass- and light-horizons, black holes' radii, the Schwartzschild metric and the Kerr metric

Mass- and light-horizons, black holes' radii, the Schwartzschild metric and the Kerr metric 006-010 Thiey De Mees Mass- and light-hoizons, blak holes' adii, the Shwatzshild meti and the Ke meti mpoved alulus. (using gavitomagnetism) T. De Mees - thieydm@pandoa.be Abstat Blak holes geneally ae

More information

Revised Newtonian Formula of Gravity and Equation of Cosmology in Flat Space-Time Transformed from Schwarzschild Solution

Revised Newtonian Formula of Gravity and Equation of Cosmology in Flat Space-Time Transformed from Schwarzschild Solution Intenational Jounal of Astonomy and Astophysis,,, 6-8 http://dx.doi.og/.46/ijaa.. Published Online Mah (http://www.sip.og/jounal/ijaa) evised Newtonian Fomula of Gavity and Equation of Cosmology in Flat

More information

Chapter 4. Sampling of Continuous-Time Signals

Chapter 4. Sampling of Continuous-Time Signals Chapte 4 Sampling of Continuous-Time Signals 1 Intodution Disete-time signals most ommonly ou as epesentations of sampled ontinuous-time signals. Unde easonable onstaints, a ontinuous-time signal an be

More information

Relativity for Global Navigation Satellite Systems

Relativity for Global Navigation Satellite Systems Relativity fo Global Navigation Satellite Systems Notes by Anna Heffenan based on the Living eviews atile, Relativity in the Global Positioning Systems, Neil Ashby, Living Rev. Relativity 6, (003),1 whih

More information

Khmelnik S.I. Mathematical Model of Dust Whirl

Khmelnik S.I. Mathematical Model of Dust Whirl Khmelnik S.I. Mathematial Model of Dust Whil Abstat The question of the soue of enegy in a dust whil is onsideed. Atmosphei onditions annot be the sole soue of enegy, as suh dust whils exist on Mas, whee

More information

Experiment 1 Electric field and electric potential

Experiment 1 Electric field and electric potential Expeiment 1 Eleti field and eleti potential Pupose Map eleti equipotential lines and eleti field lines fo two-dimensional hage onfiguations. Equipment Thee sheets of ondutive papes with ondutive-ink eletodes,

More information

TORSIONAL VIBRATIONS IN THE SAW UNIT OF A KIND OF CIRCULAR SAW. NUMERICAL INVESTIGATIONS OF THE NATURAL FREQUENCIES AND MODE SHAPES

TORSIONAL VIBRATIONS IN THE SAW UNIT OF A KIND OF CIRCULAR SAW. NUMERICAL INVESTIGATIONS OF THE NATURAL FREQUENCIES AND MODE SHAPES TRIESKOVÉ A BEZTRIESKOVÉ OBRÁBANIE DREVA, 8(): 7 78, Zvolen, Tehniká univezita vo Zvolene, ISBN 978-8-8-85- 7 TORSIONAL VIBRATIONS IN THE SAW UNIT OF A KIND OF CIRCULAR SAW NUMERICAL INVESTIGATIONS OF

More information

The Kerr-metric, mass- and light-horizons, and black holes' radii.

The Kerr-metric, mass- and light-horizons, and black holes' radii. 006 Thiey De Mees The Ke-meti, mass- and light-hoizons, and blak holes' adii. (using the Analogue Maxwell theoy) T. De Mees - thieydm @ pandoa.be Abstat Blak holes an geneally be defined as stella objets

More information

dp p v= = ON SHOCK WAVES AT LARGE DISTANCES FROM THE PLACE OF THEIR ORIGIN By Lev D. Landau J. Phys. U.S.S.R. 9, 496 (1945).

dp p v= = ON SHOCK WAVES AT LARGE DISTANCES FROM THE PLACE OF THEIR ORIGIN By Lev D. Landau J. Phys. U.S.S.R. 9, 496 (1945). ON SHOCK WAVES AT LARGE DISTANCES FROM THE PLACE OF THEIR ORIGIN By Lev D. Landau J. Phys. U.S.S.R. 9, 496 (1945). It is shown that at lage distanes fom the body, moving with a. veloity exeeding that of

More information

SAMPLE LABORATORY SESSION FOR JAVA MODULE B. Calculations for Sample Cross-Section 2

SAMPLE LABORATORY SESSION FOR JAVA MODULE B. Calculations for Sample Cross-Section 2 SAMPLE LABORATORY SESSION FOR JAVA MODULE B Calulations fo Sample Coss-Setion. Use Input. Setion Popeties The popeties of Sample Coss-Setion ae shown in Figue and ae summaized below. Figue : Popeties of

More information

FORMATION FLYING DYNAMICS OF MICRO- SATELLITES NEAR EQUATORIAL LOW ORBITS UNDER THE INFLUENCE OF J2 AND J3

FORMATION FLYING DYNAMICS OF MICRO- SATELLITES NEAR EQUATORIAL LOW ORBITS UNDER THE INFLUENCE OF J2 AND J3 IAA-AAS-DyCoSS-14-08-09 FORMATION FLYING DYNAMICS OF MICRO- SATELLITES NEAR EQUATORIAL LOW ORBITS UNDER THE INFLUENCE OF J AND J Haijono Djojodihadjo 1 The pesent wok assess the effet of the Eath s oblateness

More information

Study of the Endface Friction of the Revolving Vane Mechanism

Study of the Endface Friction of the Revolving Vane Mechanism Pudue Univesity Pudue e-pubs Intenational Compesso Engineeing Confeene Shool of Mehanial Engineeing 010 Study of the Endfae Fition of the Revolving Vane Mehanism Alison Subiantoo Shool of Mehanial and

More information

Non-Ideal Gas Behavior P.V.T Relationships for Liquid and Solid:

Non-Ideal Gas Behavior P.V.T Relationships for Liquid and Solid: hemodynamis Non-Ideal Gas Behavio.. Relationships fo Liquid and Solid: An equation of state may be solved fo any one of the thee quantities, o as a funtion of the othe two. If is onsideed a funtion of

More information

Suppose you have a bank account that earns interest at rate r, and you have made an initial deposit of X 0

Suppose you have a bank account that earns interest at rate r, and you have made an initial deposit of X 0 IOECONOMIC MODEL OF A FISHERY (ontinued) Dynami Maximum Eonomi Yield In ou deivation of maximum eonomi yield (MEY) we examined a system at equilibium and ou analysis made no distintion between pofits in

More information

Recitation PHYS 131. must be one-half of T 2

Recitation PHYS 131. must be one-half of T 2 Reitation PHYS 131 Ch. 5: FOC 1, 3, 7, 10, 15. Pobles 4, 17, 3, 5, 36, 47 & 59. Ch 5: FOC Questions 1, 3, 7, 10 & 15. 1. () The eloity of a has a onstant agnitude (speed) and dietion. Sine its eloity is

More information

On the indirect e ect in the Stokes±Helmert method of geoid determination

On the indirect e ect in the Stokes±Helmert method of geoid determination Jounal of Geodesy (1999) 7: 87±9 On the indiet e et in the Stokes±Helmet method of geoid detemination L. E. SjoÈ beg, H. Nahavandhi oyal Institute of Tehnology, Depatment of Geodesy and Photogammety, S-100

More information

Generalized Vapor Pressure Prediction Consistent with Cubic Equations of State

Generalized Vapor Pressure Prediction Consistent with Cubic Equations of State Genealized Vapo Pessue Pedition Consistent with Cubi Equations of State Laua L. Petasky and Mihael J. Misovih, Hope College, Holland, MI Intodution Equations of state may be used to alulate pue omponent

More information

Eddy Currents and Magnetic Calibrations in LDX using a Copper Plasma. D.P. Boyle, PPPL M.E. Mauel, D.T. Garnier, Columbia J.

Eddy Currents and Magnetic Calibrations in LDX using a Copper Plasma. D.P. Boyle, PPPL M.E. Mauel, D.T. Garnier, Columbia J. Eddy Cuents and Magneti Calibations in LDX using a Coppe Plasma D.P. Boyle PPPL M.E. Mauel D.T. Ganie Columbia J. Kesne MIT PSFC Coppe Plasma Oveview LDX Magnetis Goals Calibate magneti diagnostis positions

More information

On the Optimal Multirate Control of Networked Control Systems

On the Optimal Multirate Control of Networked Control Systems WSEAS RANSACIONS on SYSES On the Optimal ultiate Contol of Netwoked Contol Systems ZHANG XIANG, XIAO JIAN Key Laboatoy of agneti Suspension ehnology and aglev Vehile inisty of Eduation Shool of Eletial

More information

Physics 181. Assignment 4

Physics 181. Assignment 4 Physics 181 Assignment 4 Solutions 1. A sphee has within it a gavitational field given by g = g, whee g is constant and is the position vecto of the field point elative to the cente of the sphee. This

More information

7.2. Coulomb s Law. The Electric Force

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

Modeling of Vector Controlled Induction Motor in Different Reference Frames

Modeling of Vector Controlled Induction Motor in Different Reference Frames ISSN (Online) : 39-8753 ISSN (Pint) : 347-67 Intenational Jounal of Innovative Reseah in Siene, Engineeing and Tehnology Volume 3, Speial Issue 3, Mah 4 4 Intenational Confeene on Innovations in Engineeing

More information

IMPLEMENTATION OF MUR S ABSORBING BOUNDARIES WITH PERIODIC STRUCTURES TO SPEED UP THE DESIGN PROCESS USING FINITE-DIFFERENCE TIME-DOMAIN METHOD

IMPLEMENTATION OF MUR S ABSORBING BOUNDARIES WITH PERIODIC STRUCTURES TO SPEED UP THE DESIGN PROCESS USING FINITE-DIFFERENCE TIME-DOMAIN METHOD Pogess In Eletomagnetis Reseah, PIER 58, 101 114, 006 IMPLEMENTATION OF MUR S ABSORBING BOUNDARIES WITH PERIODIC STRUCTURES TO SPEED UP THE DESIGN PROCESS USING FINITE-DIFFERENCE TIME-DOMAIN METHOD G.

More information

E(r,t) = e 3. r 3. (b) Show that the transverse current, J t,is 3n(n e 3 ) e 3

E(r,t) = e 3. r 3. (b) Show that the transverse current, J t,is 3n(n e 3 ) e 3 Polem Set 3 (Jakson 6.20).. An example of the pesevation of ausality and finite speed of popagation in spite of the use of the Coulomg gauge is affoded y a unit stength dipole soue that is flashed on and

More information

Vibrational Modes and Instabilities of a Dust Particle Pair in a Complex Plasma

Vibrational Modes and Instabilities of a Dust Particle Pair in a Complex Plasma > REPLACE THIS LINE WITH YOUR PAPER IDENTIFICATION NUMBER (DOUBLE-CLICK HERE TO EDIT) < Vibational Modes and Instabilities of a Dust Patile Pai in a Comple Plasma K. Qiao, L. S. Matthews, and T. W. Hyde,

More information

Electric Anisotropy, Magnetic Anisotropy, Uniaxial and Biaxial Materials, Bianisotropic Media (Definitions)

Electric Anisotropy, Magnetic Anisotropy, Uniaxial and Biaxial Materials, Bianisotropic Media (Definitions) leti nisotop agneti nisotop Uniaial and iaial ateials ianisotopi edia efinitions medium is alled eletiall anisotopi if tenso Note that and ae no longe paallel medium is magnetiall anisotopi if tenso Note

More information

Gain-Scheduled Controller Design: An Analytic Framework Directly Incorporating Non-Equilibrium Plant Dynamics

Gain-Scheduled Controller Design: An Analytic Framework Directly Incorporating Non-Equilibrium Plant Dynamics Gain-Sheduled Contolle Design: An Analyti Fameok Dietly Inopoating Non-Equilibium Plant Dynamis D.J.Leith W.E.Leithead Abstat Depatment of Eletoni & Eletial Engineeing, Univesity of Stathlyde, GLASGOW

More information

Mass Transfer (Stoffaustausch)

Mass Transfer (Stoffaustausch) Mass Tansfe (Stoffaustaush) Examination 3. August 3 Name: Legi-N.: Edition Diffusion by E. L. Cussle: none nd 3 d Test Duation: minutes The following mateials ae not pemitted at you table and have to be

More information

e sin cos i sin sin j cos k [2 POINTS] (c) Hence, determine expressions for sin sin i sin cos j sin e

e sin cos i sin sin j cos k [2 POINTS] (c) Hence, determine expressions for sin sin i sin cos j sin e EN: Continuum Mehanis Homewok : Kinematis Due : noon Fiday Febuay 4th Shool of Engineeing Bown Univesity. To analyze the defomation of a onial membane, it is poposed to use a two-dimensional onial-pola

More information

PROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.

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

7.2.1 Basic relations for Torsion of Circular Members

7.2.1 Basic relations for Torsion of Circular Members Section 7. 7. osion In this section, the geomety to be consideed is that of a long slende cicula ba and the load is one which twists the ba. Such poblems ae impotant in the analysis of twisting components,

More information

arxiv: v4 [physics.class-ph] 14 Jul 2018

arxiv: v4 [physics.class-ph] 14 Jul 2018 Noname manusipt No. will be inseted by the edito Long-Range Longitudinal Eleti Wave in Vauum Radiated by Eleti Dipole: Pat I Altay Zhakatayev, Leila Tlebaldiyeva axiv:7.v4 [physis.lass-ph] 4 Jul 8 Reeived:

More information

3.1 Random variables

3.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 information

matschek (ccm2548) Ch17-h3 chiu (57890) 1

matschek (ccm2548) Ch17-h3 chiu (57890) 1 matshek m2548) Ch17-h3 hiu 5789) 1 This pint-out should have 16 questions. Multiple-hoie questions may ontinue on the next olumn o page find all hoies efoe answeing. 1 1. points A student said, The eleti

More information

Anisotropic 2-D Wavelet Packets and Rectangular Tiling: Theory and Algorithms

Anisotropic 2-D Wavelet Packets and Rectangular Tiling: Theory and Algorithms Anisotopi -D Wavelet Pakets and Retangula Tiling: Theoy and Algoithms Dan Xu and Minh N. Do Depatment of Eletial and Compute Engineeing and Bekman Institute Univesity of Illinois at Ubana-Champaign Email:

More information

How to Obtain Desirable Transfer Functions in MIMO Systems Under Internal Stability Using Open and Closed Loop Control

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

not to be republished NCERT ELECTROMAGNETIC WAVES Chapter Eight MCQ I

not to be republished NCERT ELECTROMAGNETIC WAVES Chapter Eight MCQ I Chapte Eight ELECTROMAGNETIC WAVES MCQ I 8 One equies ev of enegy to dissoiate a abon monoxide moleule into abon and oxygen atoms The minimum fequeny of the appopiate eletomagneti adiation to ahieve the

More information

Stress, Cauchy s equation and the Navier-Stokes equations

Stress, Cauchy s equation and the Navier-Stokes equations Chapte 3 Stess, Cauchy s equation and the Navie-Stokes equations 3. The concept of taction/stess Conside the volume of fluid shown in the left half of Fig. 3.. The volume of fluid is subjected to distibuted

More information

Design of Brushless DC motor Drive System for Electric Vehicle Applications Yueying ZHU1,2, a, Xu CAO1,b,Shihai CUI1,2

Design of Brushless DC motor Drive System for Electric Vehicle Applications Yueying ZHU1,2, a, Xu CAO1,b,Shihai CUI1,2 nd Intenational Confeene on Advanes in Mehanial Engineeing and Industial Infomatis (AMEII 016) Design of Bushless DC moto Dive System fo Eleti Vehile Appliations Yueying ZHU1,, a, Xu CAO1,b,Shihai CUI1,

More information

Macroelement Modelling of Laterally Loaded Piles and Pile-groups

Macroelement Modelling of Laterally Loaded Piles and Pile-groups 1 st Intenational Confeene on Natual Hazads & Infastutue 8-30 June, 016, Chania, Geee Maoelement Modelling of Lateally Loaded Piles and Pile-goups Nikos Geolymos 1 National Tehnial Univesity of Athens

More information

A NEW VARIABLE STIFFNESS SPRING USING A PRESTRESSED MECHANISM

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 information

APPENDIX D COMPRESSIBILITY FACTOR EQUATIONS D.1 THE REDLICH KWONG EQUATION

APPENDIX D COMPRESSIBILITY FACTOR EQUATIONS D.1 THE REDLICH KWONG EQUATION AENDIX D COMRESSIBILIY FACOR EQUAIONS D.1 HE REDLICH KWONG EQUAION he Redlih-Kwong equation is atually an equation of state. It was fomulated by Otto Redlih and Joseph N. S. Kwong in 1949 [Chemial Review

More information

THEORETICAL AND EXPERIMENTAL STUDY ON DROPWISE CONDENSATION IN PLATE HEAT EXCHANGERS

THEORETICAL AND EXPERIMENTAL STUDY ON DROPWISE CONDENSATION IN PLATE HEAT EXCHANGERS Abstat THEORETICAL AND EXPERIMENTAL STUDY ON DROPWISE CONDENSATION IN PLATE HEAT EXCHANGERS V. Bendt, S. Zunft and H. Mülle-Steinhagen Geman Aeospae Cente (DLR), Stuttgat, Gemany This pape desibes the

More information

Quaternion-Based Tracking Control Law Design For Tracking Mode

Quaternion-Based Tracking Control Law Design For Tracking Mode A. M. Elbeltagy Egyptian Armed forces Conference on small satellites. 2016 Logan, Utah, USA Paper objectives Introduction Presentation Agenda Spacecraft combined nonlinear model Proposed RW nonlinear attitude

More information

A Theory of the Podkletnov Effect based on General Relativity: Anti-Gravity Force due to the Perturbed Non-Holonomic Background of Space

A Theory of the Podkletnov Effect based on General Relativity: Anti-Gravity Force due to the Perturbed Non-Holonomic Background of Space July, 007 PROGRESS IN PHYSICS Volume 3 SPECIAL REPORT A Theoy of the Podkletnov Effet based on Geneal Relativity: Anti-Gavity Foe due to the Petubed Non-Holonomi Bakgound of Spae Dmiti Rabounski and Laissa

More information

Newton s Laws, Kepler s Laws, and Planetary Orbits

Newton s Laws, Kepler s Laws, and Planetary Orbits Newton s Laws, Keple s Laws, and Planetay Obits PROBLEM SET 4 DUE TUESDAY AT START OF LECTURE 28 Septembe 2017 ASTRONOMY 111 FALL 2017 1 Newton s & Keple s laws and planetay obits Unifom cicula motion

More information

SIMPLE MODELING AND CONTROL OF PLASMA CURRENT PROFILE. François Delebecque Jean-Pierre Quadrat Sylvain Brémond Emmanuel Witrant Jean-François Artaud

SIMPLE MODELING AND CONTROL OF PLASMA CURRENT PROFILE. François Delebecque Jean-Pierre Quadrat Sylvain Brémond Emmanuel Witrant Jean-François Artaud SIMPLE MOELING AN CONTROL OF PLASMA CURRENT PROFILE Fançois elebeue Jean-Piee Quadat Sylvain Bémond Emmanuel Witant Jean-Fançois Ataud INRIA-Rouenout, BP105 78150 Le Chesnay Cedex (Fane) CEA-Cadaahe, 13108

More information

Lab #9: The Kinematics & Dynamics of. Circular Motion & Rotational Motion

Lab #9: The Kinematics & Dynamics of. Circular Motion & Rotational Motion Reading Assignment: Lab #9: The Kinematics & Dynamics of Cicula Motion & Rotational Motion Chapte 6 Section 4 Chapte 11 Section 1 though Section 5 Intoduction: When discussing motion, it is impotant to

More information

Chapter 13 Gravitation

Chapter 13 Gravitation Chapte 13 Gavitation In this chapte we will exploe the following topics: -Newton s law of gavitation, which descibes the attactive foce between two point masses and its application to extended objects

More information

Physics 235 Chapter 5. Chapter 5 Gravitation

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

On the integration of the equations of hydrodynamics

On the integration of the equations of hydrodynamics Uebe die Integation de hydodynamischen Gleichungen J f eine u angew Math 56 (859) -0 On the integation of the equations of hydodynamics (By A Clebsch at Calsuhe) Tanslated by D H Delphenich In a pevious

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Depatment Physics 8.033 Decembe 5, 003 Poblem Set 10 Solutions Poblem 1 M s y x test paticle The figue above depicts the geomety of the poblem. The position

More information

Graphs of Sine and Cosine Functions

Graphs of Sine and Cosine Functions Gaphs of Sine and Cosine Functions In pevious sections, we defined the tigonometic o cicula functions in tems of the movement of a point aound the cicumfeence of a unit cicle, o the angle fomed by the

More information

2. Equation of generalized Dynamics. Let rectangular right hand coordinate triple is fixed in three-dimensional Euclidian space.

2. Equation of generalized Dynamics. Let rectangular right hand coordinate triple is fixed in three-dimensional Euclidian space. Genealized Dynamis about Foes Ating on Chage Moving in Capaito and Solenoid. J.G. Klyushin, Ph. D. Aademy of Civil Aviation, hai of applied mathematis; e-mail: klyushin@shaping.og; mail: Intenational Club

More information

SKP-2 ALGORITHM: ON FORMING PART AND MACHINE CLUSTERS SEPARATELY

SKP-2 ALGORITHM: ON FORMING PART AND MACHINE CLUSTERS SEPARATELY Poeedings of the 1998 Paifi Confeene on Manufatuing, August 18-20, 1998, Bisbane, Queensland, Austalia SKP-2 ALGORITHM: ON FORMING PART AND MACHINE CLUSTERS SEPARATELY Susanto,S., Kennedy,R.D. and Pie,

More information

EEO 401 Digital Signal Processing Prof. Mark Fowler

EEO 401 Digital Signal Processing Prof. Mark Fowler EEO 41 Digital Signal Pocessing Pof. Mak Fowle Note Set #31 Linea Phase FIR Design Optimum Equiipple (Paks-McClellan) Reading: Sect. 1.2.4 1.2.6 of Poakis & Manolakis 1/2 Motivation The window method and

More information

4) Magnetic confinement of plasma

4) Magnetic confinement of plasma 4) Magneti onfineent of plasa Due to the shielding in the plasa, thee is alost no ontol with eleti fields. A ontol is possible with agneti fields, as patiles ae bound to the field lines. This is alled

More information

DonnishJournals

DonnishJournals DonnishJounals 041-1189 Donnish Jounal of Educational Reseach and Reviews. Vol 1(1) pp. 01-017 Novembe, 014. http:///dje Copyight 014 Donnish Jounals Oiginal Reseach Pape Vecto Analysis Using MAXIMA Savaş

More information

Analysis of Spacecraft Attitude Control

Analysis of Spacecraft Attitude Control Lehigh Univesity Lehigh Peseve Theses and Dissetations 2016 Analysis of Spacecaft Attitude Contol Shiqiao Zhu Lehigh Univesity Follow this and additional woks at: http://peseve.lehigh.edu/etd Pat of the

More information

Classical Approach to the Theory of Elementary Particles

Classical Approach to the Theory of Elementary Particles Classial Appoah to the Theoy of Elementay Patiles By Yui N. Keilman Abstat: Pesented hee is an attempt to modify /extend lassial eletodynamis (CED) in ode to enable the lassial appoah (the appoah based

More information

ASTR415: Problem Set #6

ASTR415: 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 information

Solving Problems of Advance of Mercury s Perihelion and Deflection of. Photon Around the Sun with New Newton s Formula of Gravity

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

AST 121S: The origin and evolution of the Universe. Introduction to Mathematical Handout 1

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

Circular Motion Problem Solving

Circular Motion Problem Solving iula Motion Poblem Soling Aeleation o a hange in eloity i aued by a net foe: Newton nd Law An objet aeleate when eithe the magnitude o the dietion of the eloity hange We aw in the lat unit that an objet

More information

Kinematic Motion Control of Wheeled Mobile Robots Considering Curvature Constraints

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

Light Time Delay and Apparent Position

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

COMPARING MORE THAN TWO POPULATION MEANS: AN ANALYSIS OF VARIANCE

COMPARING MORE THAN TWO POPULATION MEANS: AN ANALYSIS OF VARIANCE COMPARING MORE THAN TWO POPULATION MEANS: AN ANALYSIS OF VARIANCE To see how the piniple behind the analysis of vaiane method woks, let us onside the following simple expeiment. The means ( 1 and ) of

More information

Interaction of Feedforward and Feedback Streams in Visual Cortex in a Firing-Rate Model of Columnar Computations. ( r)

Interaction of Feedforward and Feedback Streams in Visual Cortex in a Firing-Rate Model of Columnar Computations. ( r) Supplementay mateial fo Inteaction of Feedfowad and Feedback Steams in Visual Cotex in a Fiing-Rate Model of Columna Computations Tobias Bosch and Heiko Neumann Institute fo Neual Infomation Pocessing

More information

INTEGRATION OF THE SELF-ORGANIZING MAP AND NEURAL GAS WITH MULTIDIMENSIONAL SCALING

INTEGRATION OF THE SELF-ORGANIZING MAP AND NEURAL GAS WITH MULTIDIMENSIONAL SCALING ISSN 39 4X INFORMATION TECHNOLOGY AND CONTROL, 0, Vol.40, No. INTEGRATION OF THE SELF-ORGANIZING MAP AND NEURAL GAS WITH MULTIDIMENSIONAL SCALING Olga Kuasova, Alma Molytė Vilnius Univesity, Institute

More information

Quantum theory of angular momentum

Quantum theory of angular momentum Quantum theoy of angula momentum Igo Mazets igo.mazets+e141@tuwien.ac.at (Atominstitut TU Wien, Stadionallee 2, 1020 Wien Time: Fiday, 13:00 14:30 Place: Feihaus, Sem.R. DA gün 06B (exception date 18 Nov.:

More information

ME 3600 Control Systems Frequency Domain Analysis

ME 3600 Control Systems Frequency Domain Analysis ME 3600 Contol Systems Fequency Domain Analysis The fequency esponse of a system is defined as the steady-state esponse of the system to a sinusoidal (hamonic) input. Fo linea systems, the esulting steady-state

More information

Physics 2B Chapter 22 Notes - Magnetic Field Spring 2018

Physics 2B Chapter 22 Notes - Magnetic Field Spring 2018 Physics B Chapte Notes - Magnetic Field Sping 018 Magnetic Field fom a Long Staight Cuent-Caying Wie In Chapte 11 we looked at Isaac Newton s Law of Gavitation, which established that a gavitational field

More information

Dissolution of Solid Particles in Liquids: A Shrinking Core Model

Dissolution of Solid Particles in Liquids: A Shrinking Core Model Wold Aademy of Siene, Engineeing and Tehnology 5 9 Dissolution of Solid Patiles in Liquids: A Shining oe Model Wei-Lun Hsu, Mon-Jyh Lin, and Jyh-Ping Hsu Astat The dissolution of spheial patiles in liquids

More information

STABILITY AND PARAMETER SENSITIVITY ANALYSES OF AN INDUCTION MOTOR

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

APPLICATION OF MAC IN THE FREQUENCY DOMAIN

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

Between any two masses, there exists a mutual attractive force.

Between any two masses, there exists a mutual attractive force. YEAR 12 PHYSICS: GRAVITATION PAST EXAM QUESTIONS Name: QUESTION 1 (1995 EXAM) (a) State Newton s Univesal Law of Gavitation in wods Between any two masses, thee exists a mutual attactive foce. This foce

More information

Planck Quantization of Newton and Einstein Gravitation

Planck Quantization of Newton and Einstein Gravitation Plank Quantization of Newton and Einstein Gavitation Espen Gaade Haug Nowegian Univesity of Life Sienes Mah 0, 06 Abstat In this pape we ewite the gavitational onstant based on its elationship with the

More information

An analytic calculation method on air gap flux in permanent magnet. brushless DC motor with ironless rotor

An analytic calculation method on air gap flux in permanent magnet. brushless DC motor with ironless rotor Intenational Confeene on Enegy and Envionmental Potetion ICEEP 6 An analyti alulation method on ai gap flux in pemanent magnet bushless DC moto with ionless oto Xinghua Wang,Yaolong Sheng andshugang Zhao,,

More information

Tidal forces. m r. m 1 m 2. x r 2. r 1

Tidal forces. m r. m 1 m 2. x r 2. r 1 Tidal foces Befoe we look at fee waves on the eath, let s fist exaine one class of otion that is diectly foced: astonoic tides. Hee we will biefly conside soe of the tidal geneating foces fo -body systes.

More information