THREE AXIS CON TROL OF THE HUB BLE SPACE TELE SCOPE USING TWO RE AC TION WHEELS AND MAG NETIC TORQUER BARS FOR SCI ENCE OB SER VA TIONS

Size: px
Start display at page:

Download "THREE AXIS CON TROL OF THE HUB BLE SPACE TELE SCOPE USING TWO RE AC TION WHEELS AND MAG NETIC TORQUER BARS FOR SCI ENCE OB SER VA TIONS"

Transcription

1 AAS HREE AXIS CON ROL OF HE HUB BLE SPACE ELE SCOPE USING WO RE AC ION WHEELS AND MAG NEIC ORQUER BARS FOR SCI ENCE OB SER A IONS Sn Hr-Diaz, * John Wirzbrger and Dan Smith he Hb ble Space ele scope (HS) is re nowned for its s perb point ing ac cracy of less than 10 milli-arcseconds ab so lte point ing er ror. o ac com plish this, the HS re lies on its com ple ment of for re ac tion wheel as sem blies (RWAs) for at ti tde con trol and for mag netic torqer bars (MBs) for momen tm man age ment. As with most sat el lites with re ac tion wheel con trol, the forth RWA pro ides for falt tol er ance to main tain three-axis point ing ca pabil ity shold a fail re oc cr and a wheel is lost from op er a tions. If an ad ditional fail re is en con tered, the abil ity to main tain three-axis point ing is jeopar dized. In or der to pre pare for this po ten tial sit a tion, HS Pointing Con trol Sb sys tem (PCS) eam de el oped a wo Re ac tion Wheel Sci ence (RS) control mode. his mode ti lizes two RWAs and for mag netic torqer bars to achiee three-axis sta bi li za tion and point ing ac c racy nec es sary for a con tin ed sci ence ob ser ing pro gram. his pa per pres ents the de sign of the RS mode and op er a tional con sid er ations nec es sary to pro tect the space craft while al lowing for a sb stan tial sci ence pro gram. INRODUCION he Hbble Space elescope (HS) was lanched in April 1990 and since that time has been renowned for its sperb pointing accracy of less than 10 milli-arcseconds absolte pointing error. o accomplish this, the HS relies on its complement of for reaction wheel assemblies (RWAs) for attitde control and for magnetic torqer bars (MBs) for momentm management. As with most satellites with reaction wheel control, the forth RWA proides for falt tolerance to maintain three-axis pointing capability shold a failre occr and a wheel is lost from operations. If an additional failre is encontered, the ability to maintain three-axis pointing is jeopardized. Oer the corse of the 18 year operational life of HS, there hae been two RWA anomalies on-orbit. De to the niqe ability for HS to spport sericing by the Space Shttle, these two RWAs were replaced in Febrary 1997 and March With the decommissioning of the Shttle planned prior to the HS end of life, the ability to serice the telescope will be greatly diminished. In order to prepare for this potential sitation, the HS Pointing Control Sbsystem (PCS) eam deeloped a wo Reaction Wheel Science (RS) control mode. his mode tilizes two * Emer gent Space ech nol ogies, Inc., Greenbelt, Mary land sn.hr-diaz@emergentspace.com. Honeywell echnology Soltions, Inc., Colmbia, Maryland jwirzbrger@hst.nasa.go. Lockheed Mar tin Mis sion Ser ices, Greenbelt, Mary land dsmith@hst.nasa.go. 335

2 RWAs and for magnetic torqer bars to achiee three-axis stabilization and pointing accracy necessary for a contined science obsering program. With jst two operational reaction wheels, there exists one axis where no wheel control is aailable. Magnetic torqer bars can proide control abot the wheel-less axis, bt they mst also contine to dmp momentm from the wheels to preent wheel speed satration. Becase of the redced wheel momentm management capability of the bars, the wheel speeds can exceed those from normal operations and mst be tolerated by the ehicle and safing systems. Since the MBs can only generate a torqe in the plane perpendiclar to the Earth s magnetic field, no control torqe can be generated abot the wheel-less axis if the magnetic field aligns with the wheel-less axis. his means that certain maneers or attitdes may not be feasible and can lead to an ncontrolled ehicle state. herefore, carefl planning of science attitdes and maneers is necessary, and significant modification of the crrent planning and schedling tools is foreseen. Redced-wheel pointing capability had been considered in other missions as well, notably the Far Ultraiolet Spectroscopic Explorer (FUSE) 1,2,3. Unlike FUSE, HS can experience significant aerodynamic torqe which can be on the same order of magnitde as the graity gradient torqe depending on atmospheric density conditions. HS s pointing reqirement is also more stringent which places stricter constraints for the redced-wheel operations. his paper presents the design of the RS mode and the operation considerations necessary to protect the spacecraft while allowing for a sbstantial science program. his mode can operate with any pair of reaction wheels and assmes for magnetic torqer bars and at least three gyros. he paper begins with a brief backgrond followed by the algorithms and some simlation reslts. hen key isses in planning and schedling are discssed. BACKGROUND he location and direction of the for RWAs on HS along with the definition of the ehicle frame are shown in Figre 1. he boresight of the telescope is along the +1 axis. Figre 1 Reaction wheel assembly configration here are six possible pairs of wheels that may be in RS operation. he control athority, hence the performance, depends on the particlar pair of reaction wheels sed as well as the MB capability at arios leels of aerodynamic torqe. he maximm torqe from a single wheel is 0.82 N-m for wheel speeds below 3200 RPM. As the wheel speed increases aboe this ale, the torqe otpt begins to diminish. here is also a software limit on the wheel speed for safety. 336

3 he torqe otpt from the MBs depends on their orientation relatie to the Earth s magnetic field strength and direction. he maximm dipole moment from a single bar is 3500 amp-meter 2. he maximm torqe that can be prodced by all for bars is 0.34 N-m. Figre 2 shows the distribtion of maximm graity gradient and aerodynamic torqes in the body frame as well as the maximm MB torqe and acceleration capability. he wheel-less axis directions of the six possible pairs of operating wheels are also shown Graity Gradient Aerodynamic (=9.6e-12 kg/m 3 ) otal External Bar Acceleration Capability Bar orqe Capability Figre 2 External torqe and MB capability distribtion he MB acceleration capability distribtion shows that wheel pairs (1,2) and (3,4) hae the worst control athority abot the wheel-less axis as well as the worst external torqes. Performance for these two-wheel configrations is expected to be worse than the other for twowheel configrations. For adeqate control in all three axes, there needs to be enogh control athority abot the wheel-less axis by the MBs. his means that certain maneers or attitdes may not be feasible and can lead to an ncontrolled ehicle state. herefore, carefl planning of science attitdes and maneers is necessary. In addition to the maximm torqe realizable by the MBs, their dynamics are also important. Figre 3 shows the time response from a PSPICE model of PSPICE ME Response amp-m 2 to 1800 amp-m 2 (-450 ma to 450 ma) 600 Crrent (milliamps) deg C 24 deg C 71 deg C ime (seconds) Figre 3 Magnetic orqer Electronics ime Response 337

4 the magnetic torqer electronics at three different temperatres. he general response is that of a second-order system. o characterize the torqe capability abot the wheel-less axis, the following parameter is defined: AWX w x B ( eh / M M ) w x GG AERO ( I) IaMR (1) w ( h where a w B x MR eh / RWA h RWA x eh / RWA I ehicle inertia RWA magnetic field ector in ehicle frame maneer acceleration command nit ector in the wheel - less axis in ehicle frame eh / M transformation matrix from MB frame to ehicle frame M bar dipole moments in MB frame GG graity gradient torqe AERO aerodynamic torqe ehicle anglar rate in ehicle frame transformation matrix from wheel to ehicle frame nominal wheel momentm in wheel frame A positie AWX assres controllability abot the wheel-less axis. For planning and schedling prposes, some of the ariables are assmed nominal ales. For maximm torqe from the bars abot the wheel-less axis, the bar moments in the aboe eqation are set to sign( w B (2) M max x ) where [ ] X represents the cross-prodct matrix of its argment. For maximm wheel gyroscopic torqe abot the wheel-less axis for a gien maneer axis, wheel momentm is set to h h sign w (3) RWA RWA _ nom X M ) ( x X eh / RWA Dring non-maneer periods, the acceleration command is zero and the two gyroscopic terms containing ehicle anglar elocity are negligible. Dring maneers, the acceleration command and the ehicle anglar rate are the expected ales compted by the maneer command generator. he expected maneer attitde profile is also sed in the comptation of the magnetic field in the ehicle frame as well as the graity gradient torqe and the aerodynamic torqe. ALGORIHMS One of the main considerations in the design of the RS mode for HS was to minimize changes to the existing flight software (FSW). herefore, the existing proportional-integralderiatie control architectre for the normal 3- and 4-wheel modes was left intact as mch as possible. A major deiation is the comptation of the feedback control in a wheel-orthogonalframe (WOF) to maximize controllability in the wheel plane. Feedback control is compted in 338

5 WOF with gains in the wheel-plane generally higher than the gains in the wheel-less axis since the wheels hae higher actator capability than the bars. Figre 4 illstrates WOF defined by a pair of any two wheels A and B. Wheel Plane Normal (Wheel-less axis) ŵ X ŵ N rb Wheel B w Wheel A A r A Figre 4 Wheel-Orthogonal-Frame (WOF) he transformation from the wheel orthogonal frame to the ehicle frame is gien by / WOF w A w N w X r A r B r A r B r A r A r A r B r A r B where the parameters with the ^ symbol indicate nit ectors in the ehicle frame and the other ectors are defined in Figre 4. RS Control Wheel Plane Becase the MBs are needed for both ehicle control and wheel momentm management, the primary problem to sole is how to distribte these two potentially competing tasks. hree possible approaches were explored and are described below. Method 1. he first method takes the feedback control that was transformed to the ehicle frame and sms it with the feedforward control torqe compensation comprised of momentm management torqe and other known torqes. his total torqe is sed for compting the wheel command as well as for compting the deficit torqe, which is the reqired control torqe not achieable by the wheels. he total magnetic dipole command for the torqer bars is compted from a combination of the deficit torqe and momentm management torqe. he dipole moments of the torqer bars are compted from a modified cross-prodct law described in a later section. o inclde in the feedforward control torqe, the momentm management torqe is processed throgh the cross-prodct law, and the actal momentm management torqe that can be prodced from the bars is compted as follows: B MM MM (5) 2 B MM B where MM is the dipole moments of the bars also gien in the ehicle frame. he negatie of this momentm management torqe is fed forward along with other known feedforward torqes so that the total desired control torqe is gien by (4) _ ACUAL MM (6) 339

6 C FB FB FF MM _ ACUAL FF _ OHERS (7) he control torqe that is realizable by the wheels is gien by / lim / (8) C eh RWA LIMRW RWA eh where eh / RWA is transformation from the wheel frame to the ehicle frame, and RWA/ eh is the transformation from the ehicle frame to the wheel frame. hese transformations are shown below: where C eh / RWA ra rb (9) 1 RWA / eh eh / RWA eh / RWA eh / RWA r A r B (10) nit ector nit ector of wheel A in ehicle frame of wheel B in ehicle frame Note that the total control torqe transformed to the wheel frame is limited by LIMRW before transforming back to the ehicle frame. he control torqe that is not realizable by the wheels is the deficit torqe gien by DEFICI eh / WOF lim WOF / eh C C (11) RSLMM which is to be compensated by the bars as mch as possible. he limit RSLMM, applied in WOF, controls how mch deficit torqe is compensated. he total torqe to be prodced by the bars is then a combination of the deficit control torqe and the momentm management torqe. he corresponding bar moment command is then gien by: lim M MAX lim M / eh M / eh MAX t DEFICI B B MM DEFICI 2 where M / eh is the (4x3) transformation from the ehicle frame to the MB frame, Max is the maximm moment limit, and is a scaling factor applied to achiee the fll desired control torqe (within the bar limit) and is gien by DEFICI MM (12) 2 2 B DEFICI (13) 2 B he scaling of the bar torqe is illstrated in Figre 5. With scaling, the projection of the magnetic torqe along the deficit torqe direction prodces the desired magnitde proided that 340

7 the bars are not satrated. If the magnetic field ector and the deficit torqe are aligned, then is set to 1. DEFICI orqe from nominal cross-prodct law orqe from cross-prodct law with scaling Figre 5 Scaling of the bar torqe command for control Any extraneos torqe in the wheel plane prodced by the bars is compensated by the wheels. his torqe is the desired magnetic torqe mins the actal: RWA _ M _ COMP DES ( eh / M M ) B (14) where the desired magnetic torqe is gien by (15) DES DEFICI MM _ ACUAL he total command to the wheels is pdated with the compensation torqe: (16) C C RWA _ M _ COMP he total torqe to be prodced by the wheels is gien by the pdated command lim (17) RWA RWA _ Max RWA / eh Method 2. In the second method, the desired torqe is decomposed into a non-orthogonal UN frame defined by the nit ector n along the intersection of the plane of magnetic control athority and the plane of the wheel control athority which is gien by w X B n, (18) w B X the nit ector in the plane of wheel control athority gien by w X n, (19) w n X and the nit ector in the plane of magnetic control athority gien by B n. (20) B n Figre 6 illstrates the UN coordinate frame. C 341

8 342 Figre 6 UN Coordinate Frame he control torqe that is originally compted in the wheel-orthogonal-frame is conerted to the ehicle frame and then conerted to this non-orthogonal UN frame as follows: n n n 2 2 (21) he commands to the wheels and the bars are then allocated as follows. n n eh RWA RWA Max RWA lim / _ (22) MM eh M M B B MAX 2 / lim (23) he control torqe component abot the n direction (intersection of the wheel plane and the magnetic control athority plane) is allocated flly to the wheels since the wheels hae a higher torqe capability. here is a singlarity when the magnetic field ector is aligned with the wheel-less axis. In this case, torqe abot the wheel-less axis is zero and we only command the control torqe in the wheel-plane. Method 3. he third method inoles a least-sqares soltion to optimally combine the wheel torqe and the bar torqe to prodce the desired control torqe. Since momentm management torqe mst be proided by the bars in order to redce the system momentm, it is not combined with the control torqe when sing this concrrent method.

9 he system of eqations to sole for is: C B eh / RWA B X eh / M RWA M _ C (24) where is the nknown control ariable. he aboe eqation represents an nder-determined system since there are 6 nknowns and 3 eqations, so there is no niqe soltion. One possible soltion is the minimm norm soltion which can be soled by forming an objectie fnction agmented with the aboe eqality constraint sing a Lagrange mltiplier: min J W B (25) where W is a sqare matrix of weights on the nknowns. We take the partial of J with respect to and, set each partial to zero, and sole for. After some matrix algebra, we get C 1 BW B C 1 1 W B (26) which is essentially a psedo-inerse of the eqality constraint eqation. he wheel torqe command for control and the bar moment command for control are then gien by RWA lim (1: 2) RWA _ Max (27) (4 : 6) M _ C he bar command is limited after the momentm management component is added: lim (28) M MAX t M _ C M / eh he matrix BW 1 B is not inertible when the magnetic field is aligned with the wheel-less axis since the rank of the matrix wold then be only two. In this case, we can only command the torqe in the wheel plane. he psedo-inerse soltion becomes no longer optimal when the control ariables begin to satrate, and control limits wold hae to be incorporated. For example, the problem can be posed as a qadratic programming problem with linear and ineqality constraints and soled sing established nonlinear programming methods. Comparison of Methods. he three methods of torqe allocation presented aboe were simlated in Matlab. In a case where there is no actator satration, all methods performed identically. In a case where there is a period of little or no torqe athority abot the wheel-less axis, methods 2 and 3 performed worse than method 1. he torqe athority for the second case and the attitde errors for all three methods are shown in Figre 7(a) and (b), respectiely. MM 343

10 0.2 orqe Athority Abot Wheel-Less Axis (Wx) 4.5 x 104 re Error Magnitde Nm re Error, arc-seconds Methods 2 and 3 Method t (sec) (a) t (sec) (b) Figre 7 (a) AWX, (b) Attitde Error he lighter cre in Figre 7(b) corresponds to methods 2 and 3 whose attitde control error is abot 10 times worse than that of method 1. he reason is becase the wheels satrate in methods 2 and 3 and ths lose controllability. Figre 8(a) shows the wheel torqe history for Method 1, and Figre 8(b) shows the wheel torqe history for Methods 2 and Wheel orqe RWA 1 RWA 2 RWA 3 RWA Wheel orqe RWA 1 RWA 2 RWA 3 RWA 4 N-m 0 N-m t (sec) t (sec) (a) (b) Figre 8 Wheel orqe for (a) Method 1 (b) Methods 2 and 3 For Method 2 or 3 to be iable, actator limits will hae to be addressed. his will be the sbject of a ftre paper. Momentm Management For minimal FSW change, the nominal momentm management (MM) proportionalderiatie control law is maintained except that the MM command is limited to be only along the direction that wold not corrpt control. If the control torqe can be flly achieed, i.e., there is sfficient torqe athority abot the wheel-less axis, then sch restriction is not necessary. Howeer, to aoid redcing this torqe athority with nnecessary corrption from the momentm management torqe, only the n -component of the nominal momentm management torqe,, is applied: MM 344

11 MM n n (29) When the wheel-less axis is aligned with the magnetic field ector within a tolerance, the intersection line is not defined. In this case, most of the realizable momentm management torqe is already in the wheel plane and only the component of momentm management torqe that is in the wheel plane is commanded. Maneer Planning For maneer planning in RS, the same 3- and 4-wheel command generator algorithm is sed which is based on an eigen-axis maneer from an initial attitde to a final attitde. he command generator determines the ehicle rate profile for the control system to follow that satisfies the maximm rate, maximm acceleration, and maximm jerk parameters. For 3- and 4- wheel operations, these parameters are fixed for all maneers. For 2-wheel operation, on the other hand, these parameters hae to be specified per maneer becase the amont of control athority abot the wheel-less axis is dependent on the wheel configration, the ehicle attitde and orbit, atmospheric density, as well as the magnetic field strength and alignment. he comptation of these parameters, which will be performed on the grond by planning and schedling, is discssed in the next section. OPERAIONAL CONSIDERAIONS With jst two reaction wheels, there exists one axis where no wheel control is aailable. Magnetic torqe bars can proide control abot the wheel-less axis, bt they mst also contine to dmp momentm from the wheels. Becase of the redced wheel momentm management capability of the bars, the wheel speeds can get high. o allow for operation in this mode, the wheel speed safing limit is increased bt to a leel that wold still allow enogh margin for safe mode. For arios operational modes of RS, the minimm torqe athority has to be specified AWX AWX Min (30) Becase of the dependency on the aerodynamic drag torqe, the minimm torqe athority aries with the expected leel of the atmospheric density. Attitde Hold and Science For attitde hold and science modes, AWX mst be sfficient to preent nallowable attitde errors and maintain RWA momentm adeqately to allow for ehicle maneers and preent safemodes. It is obios that AWX mst be greater than 0 N-m at all times dring science interals. Once AWX becomes negatie, the external torqes become larger than the control torqes and the ehicle wold incr attitde errors and loss of lock on the target. At 0 N-m, there is no additional torqe aailable for the feedback control loop as all the control torqe aailable is dedicated to controlling expected external torqes. Some attitde error is acceptable dring attitde hold portions, as long as the schedling system protects for the largest error possible. Howeer, once AWX becomes negatie, momentm management of the RWAs is no longer garanteed and wheel speeds increase beyond acceptable limits. wo major performance criteria are attitde errors and wheel speeds. Figre 9 shows the maximm attitde errors relatie to minimm AWX from a Monte Carlo simlation of 1000 random initial conditions for wheel configration (1,2) with each case simlated for 3.5 orbits and MM 345

12 with an aero density of 3.5e-12 kg/m 3. Figre 10 shows the wheel speeds relatie to the minimm AWX for all 1000 cases. (a) (b) Figre 9 Maximm attitde error (a) all 1000 cases shown, (b) Close-p 346

13 Figre 10 Maximm Wheel Speeds Based on the simlations, setting AWX Min to N-m wold maintain attitde errors below the science acqisition search radii and maintain RWA speeds below the proposed safing speed limit. Once HS transitions to RS science mode, the absolte pointing error is maintained below 20 milli-arcseconds for the worst case wheel geometry with the inclsion of attitde information from the Fine Gidance Sensors. his increase in pointing error oer the nominal science pointing error with three RWAs of less than 10 milli-arcseconds is primarily de to the lagged response of the MB as compared to the RWAs. Figre 11 shows the probability of haing a certain torqe margin for the entire dration of a 3.5 orbit science obseration at a specific time for wheel configration (1,2) and (1,3) for aero density ale of 3.5e-12 kg/m 3. As expected, the efficiency is lower for the (1,2) configration becase the torqe athority abot the wheel-less axis is less. Positie torqe margin occrs approximately 40 percent of the time for wheel configration (1,2) and 60 percent of the time for wheel configration (1,3). When performing the calclation of AWX oer attitde hold or science interals, the selection of atmospheric density is critical. he density chosen mst be the maximm density expected to be seen on-orbit dring the interal or higher, not a mean ale. his ensres torqe athority and bilds in a slight pad for times when the actal density is lower than the schedling density. If a mean density ale is sed, AWX Min mst be reealated. 347

14 Wheel Configration (1,3) Wheel Configration (1,2) Figre 11 Science Planning Efficiency One of HS s defining performance characteristics is its jitter leel below 7 milli-arcseconds allowing for clear scientific images. Since HS is relatiely insensitie to jitter abot the boresight axis, this metric is defined as the Root-Smmed Sqared of the 60-second standard deiation of the pointing error in 2 and 3. With the increase in actator rise time incrred by sing the MB as opposed to a RWA, the fine control of HS is degraded. he amont of degradation of jitter is dependent on the remaining two RWAs. he RWA 1-2 and the RWA 3-4 configrations place the wheel-less axis in the 2/3 plane. his reslts in the maximm jitter transferred to the boresight as the total affect of the increased jitter de to the MBs is captred in the jitter calclation. Anticipated jitter leels for the worst-case wheel configrations are on the order of 12 milli-arcseconds, while jitter for the other configrations meets the 7 milli-arcseconds goal. he impact on the science mission de to the increase jitter wold be minimal shold HS fall into a worst-case wheel configration. Similarly, the ability for HS to track a moing target, sch as a planet or comet, cold be affected. his moing target track makes se of feedforward acceleration to keep the target in the field of iew. hrogh carefl gain selection, no significant degradation was obsered in moing target tracking. hrogh carefl design of the attitde hold and science modes, there is minimal performance degradation in RS mode. Howeer, a transition to RS is not withot cost. he biggest impact to these modes is the ability to schedle attitdes that hae sfficient torqe margin in the wheelless axis to control attitde errors and RWA speeds. Maneers For a gien maneer, the command generator parameters max, maximm elocity; Amax, maximm acceleration; and start, maneer start time, that satisfy the minimm AWX 348

15 reqirement will hae to be determined. A combination of these command generator parameters that optimizes an objectie fnction, e.g., minimm end-of-maneer time, can be fond: min start slew (31), max,amax start sbject to AWX AWX Min (32) he parameter start is to allow delays in the maneer for a possible improement in the magnetic field and orbit, and slew is the maneer dration compted in the command generator. he lower and the pper limits of the sole-for command generator parameters in the search space will hae to be specified for the arios 2-wheel configrations. In addition, the minimm torqe-athority abot the wheel-less axis AWX Min has to be specified. Becase aerodynamic torqe can be a significant torqe on the HS, all of these parameters will hae to be specified for different leels of atmospheric density. A method for soling the optimization problem is a global parametric search based on discrete ales of the sole-for parameters where each combination of the sole-for parameters is simlated oer the maneer dration. Cases that iolate the AWX constraint are rejected, and the combination that minimizes the objectie fnction is selected. Recall that gyroscopic torqe in the comptation of AWX in Eqn. (1) is negligible except dring maneers. Since the momentm of the wheels at the start of a maneer cannot be predicted, AWX is compted sing a nominal wheel momentm ale as shown in Eqn. (3). Figre 12 shows plots of AWX from a sample simlation. he top plot shows AWX withot the wheel gyroscopic term sing nominal ales as well as the actal ales. he bottom plot shows AWX with the gyroscopic term sing a nominal momentm in each wheel of 250 N-m-s as well as the actal wheel momentm. Figre 12 shows the attitde error corresponding to this case. 0.2 Nominal and Actal abot Wx Nm 0 M -abs(gg+aero-wcxiwc) M -abs(gg+aero-wactxiwact) t (sec) Nominal and Actal AWX abot Wx 0.2 Nm 0 M - wcxhnom -abs(gg+aero-wcxiwc) M - wactxhact -abs(gg+aero-wactxiwact) t (sec) Figre 12 Nominal AWX sing h RWA_nom =250 N-m-s and the actal AWX for a sample case 349

16 6 otal attitde error 5 4 deg ime, seconds Figre 13 Attitde error corresponding to the case shown in Figre 12 he gyroscopic torqe from the npredictable momentm of the wheels can hae a significant effect on the control athority abot the wheel-less axis. For this particlar case shown, the actal wheel momentm reaches a maximm ale of abot 350 N-m-s per wheel. o assre that cases sch as this are not planned, the optimization problem gien by Eqns. (31) and (32) is soled sing a larger ale of the nominal wheel momentm or a larger ale of AWX Min. It is possible that a feasible soltion does not exist for a gien desired maneer. If no sitable maneer exists, the science timeline wold hae to be modified accordingly. CONCLUSION A preliminary design of the two-reaction wheel science (RS) mode for the Hbble Space elescope was presented. hree methods of combining wheel torqe and magnetic bar torqe for RS control were deried and compared. he first method was chosen for the preliminary RS design becase it has the least impact to the existing flight software code. Withot accommodation of actator limits, the latter two methods perform worse than the first when there is a period of negatie torqe athority abot the wheel-less axis. Consideration of actator limits in the latter two methods is the sbject of a ftre paper. Operational considerations in terms of the torqe-athority abot the wheel-less axis were also presented. Simlation reslts of both attitde hold and maneers show that acceptable performance can be had when there is sfficient torqe athority. ACKNOWLEDGMEN he athors wold like to thank Landis Markley and Peiman Maghami of NASA Goddard Space Flight Center for their peer reiew and sggestions of the RS design. REFERENCES 1. Roberts, Bryce A., et al, hree-axis Attitde Control with wo Reaction Wheels and Magnetic orqer Bars, AIAA Gidance, Naigation, and Control Conference and Exhibit, Paper 5245, Proidence, Rhode Island, Agst Krk, Jeffrey W., et al, FUSE In-Orbit Attitde Control with wo Reaction Wheels and No Gyroscopes, SPIE, ol. 4854, Paper 72, Sahnow, Daid J., Operations with the new FUSE obseratory: three-axis control with one reaction wheel, SPIE, ol. 6266, Paper 2,

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2 MATH 307 Vectors in Rn Dr. Neal, WKU Matrices of dimension 1 n can be thoght of as coordinates, or ectors, in n- dimensional space R n. We can perform special calclations on these ectors. In particlar,

More information

The Real Stabilizability Radius of the Multi-Link Inverted Pendulum

The Real Stabilizability Radius of the Multi-Link Inverted Pendulum Proceedings of the 26 American Control Conference Minneapolis, Minnesota, USA, Jne 14-16, 26 WeC123 The Real Stabilizability Radis of the Mlti-Link Inerted Pendlm Simon Lam and Edward J Daison Abstract

More information

We automate the bivariate change-of-variables technique for bivariate continuous random variables with

We automate the bivariate change-of-variables technique for bivariate continuous random variables with INFORMS Jornal on Compting Vol. 4, No., Winter 0, pp. 9 ISSN 09-9856 (print) ISSN 56-558 (online) http://dx.doi.org/0.87/ijoc.046 0 INFORMS Atomating Biariate Transformations Jeff X. Yang, John H. Drew,

More information

Adaptive Fault-tolerant Control with Control Allocation for Flight Systems with Severe Actuator Failures and Input Saturation

Adaptive Fault-tolerant Control with Control Allocation for Flight Systems with Severe Actuator Failures and Input Saturation 213 American Control Conference (ACC) Washington, DC, USA, Jne 17-19, 213 Adaptive Falt-tolerant Control with Control Allocation for Flight Systems with Severe Actator Failres and Inpt Satration Man Wang,

More information

u P(t) = P(x,y) r v t=0 4/4/2006 Motion ( F.Robilliard) 1

u P(t) = P(x,y) r v t=0 4/4/2006 Motion ( F.Robilliard) 1 y g j P(t) P(,y) r t0 i 4/4/006 Motion ( F.Robilliard) 1 Motion: We stdy in detail three cases of motion: 1. Motion in one dimension with constant acceleration niform linear motion.. Motion in two dimensions

More information

Technical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty

Technical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty Technical Note EN-FY160 Revision November 30, 016 ODiSI-B Sensor Strain Gage Factor Uncertainty Abstract Lna has pdated or strain sensor calibration tool to spport NIST-traceable measrements, to compte

More information

1 The space of linear transformations from R n to R m :

1 The space of linear transformations from R n to R m : Math 540 Spring 20 Notes #4 Higher deriaties, Taylor s theorem The space of linear transformations from R n to R m We hae discssed linear transformations mapping R n to R m We can add sch linear transformations

More information

LPV control of an active vibration isolation system

LPV control of an active vibration isolation system 29 American Control Conference Hyatt Regency Rierfront, St. Lois, MO, USA Jne 1-12, 29 ThC15.1 LPV control of an actie ibration isolation system W.H.T.M. Aangenent, C.H.A. Criens, M.J.G. an de Molengraft,

More information

Nonparametric Identification and Robust H Controller Synthesis for a Rotational/Translational Actuator

Nonparametric Identification and Robust H Controller Synthesis for a Rotational/Translational Actuator Proceedings of the 6 IEEE International Conference on Control Applications Mnich, Germany, October 4-6, 6 WeB16 Nonparametric Identification and Robst H Controller Synthesis for a Rotational/Translational

More information

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University 9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)

More information

Optimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance

Optimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance Optimal Control of a Heterogeneos Two Server System with Consideration for Power and Performance by Jiazheng Li A thesis presented to the University of Waterloo in flfilment of the thesis reqirement for

More information

Relativity II. The laws of physics are identical in all inertial frames of reference. equivalently

Relativity II. The laws of physics are identical in all inertial frames of reference. equivalently Relatiity II I. Henri Poincare's Relatiity Principle In the late 1800's, Henri Poincare proposed that the principle of Galilean relatiity be expanded to inclde all physical phenomena and not jst mechanics.

More information

EE2 Mathematics : Functions of Multiple Variables

EE2 Mathematics : Functions of Multiple Variables EE2 Mathematics : Fnctions of Mltiple Variables http://www2.imperial.ac.k/ nsjones These notes are not identical word-for-word with m lectres which will be gien on the blackboard. Some of these notes ma

More information

Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS

Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS 3. System Modeling Mathematical Modeling In designing control systems we mst be able to model engineered system dynamics. The model of a dynamic system

More information

UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL

UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL 8th International DAAAM Baltic Conference "INDUSTRIAL ENGINEERING - 19-1 April 01, Tallinn, Estonia UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL Põdra, P. & Laaneots, R. Abstract: Strength analysis is a

More information

Sources of Non Stationarity in the Semivariogram

Sources of Non Stationarity in the Semivariogram Sorces of Non Stationarity in the Semivariogram Migel A. Cba and Oy Leangthong Traditional ncertainty characterization techniqes sch as Simple Kriging or Seqential Gassian Simlation rely on stationary

More information

Lecture 3. (2) Last time: 3D space. The dot product. Dan Nichols January 30, 2018

Lecture 3. (2) Last time: 3D space. The dot product. Dan Nichols January 30, 2018 Lectre 3 The dot prodct Dan Nichols nichols@math.mass.ed MATH 33, Spring 018 Uniersity of Massachsetts Janary 30, 018 () Last time: 3D space Right-hand rle, the three coordinate planes 3D coordinate system:

More information

Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers

Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers D.R. Espinoza-Trejo and D.U. Campos-Delgado Facltad de Ingeniería, CIEP, UASLP, espinoza trejo dr@aslp.mx Facltad de Ciencias,

More information

Reduction of over-determined systems of differential equations

Reduction of over-determined systems of differential equations Redction of oer-determined systems of differential eqations Maim Zaytse 1) 1, ) and Vyachesla Akkerman 1) Nclear Safety Institte, Rssian Academy of Sciences, Moscow, 115191 Rssia ) Department of Mechanical

More information

Numerical Model for Studying Cloud Formation Processes in the Tropics

Numerical Model for Studying Cloud Formation Processes in the Tropics Astralian Jornal of Basic and Applied Sciences, 5(2): 189-193, 211 ISSN 1991-8178 Nmerical Model for Stdying Clod Formation Processes in the Tropics Chantawan Noisri, Dsadee Skawat Department of Mathematics

More information

Lesson 81: The Cross Product of Vectors

Lesson 81: The Cross Product of Vectors Lesson 8: The Cross Prodct of Vectors IBHL - SANTOWSKI In this lesson yo will learn how to find the cross prodct of two ectors how to find an orthogonal ector to a plane defined by two ectors how to find

More information

Chapter 6 Momentum Transfer in an External Laminar Boundary Layer

Chapter 6 Momentum Transfer in an External Laminar Boundary Layer 6. Similarit Soltions Chapter 6 Momentm Transfer in an Eternal Laminar Bondar Laer Consider a laminar incompressible bondar laer with constant properties. Assme the flow is stead and two-dimensional aligned

More information

MATH2715: Statistical Methods

MATH2715: Statistical Methods MATH275: Statistical Methods Exercises III (based on lectres 5-6, work week 4, hand in lectre Mon 23 Oct) ALL qestions cont towards the continos assessment for this modle. Q. If X has a niform distribtion

More information

The Cross Product of Two Vectors in Space DEFINITION. Cross Product. u * v = s ƒ u ƒƒv ƒ sin ud n

The Cross Product of Two Vectors in Space DEFINITION. Cross Product. u * v = s ƒ u ƒƒv ƒ sin ud n 12.4 The Cross Prodct 873 12.4 The Cross Prodct In stdying lines in the plane, when we needed to describe how a line was tilting, we sed the notions of slope and angle of inclination. In space, we want

More information

Low-emittance tuning of storage rings using normal mode beam position monitor calibration

Low-emittance tuning of storage rings using normal mode beam position monitor calibration PHYSIAL REVIEW SPEIAL TOPIS - AELERATORS AND BEAMS 4, 784 () Low-emittance tning of storage rings sing normal mode beam position monitor calibration A. Wolski* Uniersity of Lierpool, Lierpool, United Kingdom

More information

MAT389 Fall 2016, Problem Set 6

MAT389 Fall 2016, Problem Set 6 MAT389 Fall 016, Problem Set 6 Trigonometric and hperbolic fnctions 6.1 Show that e iz = cos z + i sin z for eer comple nmber z. Hint: start from the right-hand side and work or wa towards the left-hand

More information

Online Stochastic Matching: New Algorithms and Bounds

Online Stochastic Matching: New Algorithms and Bounds Online Stochastic Matching: New Algorithms and Bonds Brian Brbach, Karthik A. Sankararaman, Araind Sriniasan, and Pan X Department of Compter Science, Uniersity of Maryland, College Park, MD 20742, USA

More information

6.4 VECTORS AND DOT PRODUCTS

6.4 VECTORS AND DOT PRODUCTS 458 Chapter 6 Additional Topics in Trigonometry 6.4 VECTORS AND DOT PRODUCTS What yo shold learn ind the dot prodct of two ectors and se the properties of the dot prodct. ind the angle between two ectors

More information

A New Approach to Direct Sequential Simulation that Accounts for the Proportional Effect: Direct Lognormal Simulation

A New Approach to Direct Sequential Simulation that Accounts for the Proportional Effect: Direct Lognormal Simulation A ew Approach to Direct eqential imlation that Acconts for the Proportional ffect: Direct ognormal imlation John Manchk, Oy eangthong and Clayton Detsch Department of Civil & nvironmental ngineering University

More information

PROBABILISTIC APPROACHES TO STABILITY AND DEFORMATION PROBLEMS IN BRACED EXCAVATION

PROBABILISTIC APPROACHES TO STABILITY AND DEFORMATION PROBLEMS IN BRACED EXCAVATION Clemson Uniersity TigerPrints All Dissertations Dissertations 12-2011 PROBABILISTIC APPROACHES TO STABILITY AND DEFORMATION PROBLEMS IN BRACED EXCAVATION Zhe Lo Clemson Uniersity, jerry8256@gmail.com Follow

More information

The Brauer Manin obstruction

The Brauer Manin obstruction The Braer Manin obstrction Martin Bright 17 April 2008 1 Definitions Let X be a smooth, geometrically irredcible ariety oer a field k. Recall that the defining property of an Azmaya algebra A is that,

More information

Lecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2

Lecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 BIJU PATNAIK UNIVERSITY OF TECHNOLOGY, ODISHA Lectre Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 Prepared by, Dr. Sbhend Kmar Rath, BPUT, Odisha. Tring Machine- Miscellany UNIT 2 TURING MACHINE

More information

Discontinuous Fluctuation Distribution for Time-Dependent Problems

Discontinuous Fluctuation Distribution for Time-Dependent Problems Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation

More information

1. Solve Problem 1.3-3(c) 2. Solve Problem 2.2-2(b)

1. Solve Problem 1.3-3(c) 2. Solve Problem 2.2-2(b) . Sole Problem.-(c). Sole Problem.-(b). A two dimensional trss shown in the figre is made of alminm with Yong s modls E = 8 GPa and failre stress Y = 5 MPa. Determine the minimm cross-sectional area of

More information

Modelling, Simulation and Control of Quadruple Tank Process

Modelling, Simulation and Control of Quadruple Tank Process Modelling, Simlation and Control of Qadrple Tan Process Seran Özan, Tolgay Kara and Mehmet rıcı,, Electrical and electronics Engineering Department, Gaziantep Uniersity, Gaziantep, Trey bstract Simple

More information

Safe Manual Control of the Furuta Pendulum

Safe Manual Control of the Furuta Pendulum Safe Manal Control of the Frta Pendlm Johan Åkesson, Karl Johan Åström Department of Atomatic Control, Lnd Institte of Technology (LTH) Box 8, Lnd, Sweden PSfrag {jakesson,kja}@control.lth.se replacements

More information

Direct linearization method for nonlinear PDE s and the related kernel RBFs

Direct linearization method for nonlinear PDE s and the related kernel RBFs Direct linearization method for nonlinear PDE s and the related kernel BFs W. Chen Department of Informatics, Uniersity of Oslo, P.O.Box 1080, Blindern, 0316 Oslo, Norway Email: wenc@ifi.io.no Abstract

More information

The Linear Quadratic Regulator

The Linear Quadratic Regulator 10 The Linear Qadratic Reglator 10.1 Problem formlation This chapter concerns optimal control of dynamical systems. Most of this development concerns linear models with a particlarly simple notion of optimality.

More information

SECTION 6.7. The Dot Product. Preview Exercises. 754 Chapter 6 Additional Topics in Trigonometry. 7 w u 7 2 =?. 7 v 77w7

SECTION 6.7. The Dot Product. Preview Exercises. 754 Chapter 6 Additional Topics in Trigonometry. 7 w u 7 2 =?. 7 v 77w7 754 Chapter 6 Additional Topics in Trigonometry 115. Yo ant to fly yor small plane de north, bt there is a 75-kilometer ind bloing from est to east. a. Find the direction angle for here yo shold head the

More information

Reducing Conservatism in Flutterometer Predictions Using Volterra Modeling with Modal Parameter Estimation

Reducing Conservatism in Flutterometer Predictions Using Volterra Modeling with Modal Parameter Estimation JOURNAL OF AIRCRAFT Vol. 42, No. 4, Jly Agst 2005 Redcing Conservatism in Fltterometer Predictions Using Volterra Modeling with Modal Parameter Estimation Rick Lind and Joao Pedro Mortaga University of

More information

FOUNTAIN codes [3], [4] provide an efficient solution

FOUNTAIN codes [3], [4] provide an efficient solution Inactivation Decoding of LT and Raptor Codes: Analysis and Code Design Francisco Lázaro, Stdent Member, IEEE, Gianligi Liva, Senior Member, IEEE, Gerhard Bach, Fellow, IEEE arxiv:176.5814v1 [cs.it 19 Jn

More information

BLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students

BLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students BLOOM S TAXONOMY Topic Following Bloom s Taonomy to Assess Stdents Smmary A handot for stdents to eplain Bloom s taonomy that is sed for item writing and test constrction to test stdents to see if they

More information

Model (In-)Validation from a H and µ perspective

Model (In-)Validation from a H and µ perspective Model (In-)Validation from a H and µ perspectie Wolfgang Reinelt Department of Electrical Engineering Linköping Uniersity, S-581 83 Linköping, Seden WWW: http://.control.isy.li.se/~olle/ Email: olle@isy.li.se

More information

Prandl established a universal velocity profile for flow parallel to the bed given by

Prandl established a universal velocity profile for flow parallel to the bed given by EM 0--00 (Part VI) (g) The nderlayers shold be at least three thicknesses of the W 50 stone, bt never less than 0.3 m (Ahrens 98b). The thickness can be calclated sing Eqation VI-5-9 with a coefficient

More information

ANOVA INTERPRETING. It might be tempting to just look at the data and wing it

ANOVA INTERPRETING. It might be tempting to just look at the data and wing it Introdction to Statistics in Psychology PSY 2 Professor Greg Francis Lectre 33 ANalysis Of VAriance Something erss which thing? ANOVA Test statistic: F = MS B MS W Estimated ariability from noise and mean

More information

Thermal balance of a wall with PCM-enhanced thermal insulation

Thermal balance of a wall with PCM-enhanced thermal insulation Thermal balance of a wall with PCM-enhanced thermal inslation E. Kossecka Institte of Fndamental Technological esearch of the Polish Academy of Sciences, Warsaw, Poland J. Kośny Oak idge National aboratory;

More information

Visualisations of Gussian and Mean Curvatures by Using Mathematica and webmathematica

Visualisations of Gussian and Mean Curvatures by Using Mathematica and webmathematica Visalisations of Gssian and Mean Cratres by Using Mathematica and webmathematica Vladimir Benić, B. sc., (benic@grad.hr), Sonja Gorjanc, Ph. D., (sgorjanc@grad.hr) Faclty of Ciil Engineering, Kačićea 6,

More information

University of Bristol

University of Bristol Uniersity of Bristol AEROGUST Workshop 27 th - 28 th April 2017, Uniersity of Lierpool Presented by Robbie Cook and Chris Wales Oeriew Theory Nonlinear strctral soler copled with nsteady aerodynamics Gst

More information

ON TRANSIENT DYNAMICS, OFF-EQUILIBRIUM BEHAVIOUR AND IDENTIFICATION IN BLENDED MULTIPLE MODEL STRUCTURES

ON TRANSIENT DYNAMICS, OFF-EQUILIBRIUM BEHAVIOUR AND IDENTIFICATION IN BLENDED MULTIPLE MODEL STRUCTURES ON TRANSIENT DYNAMICS, OFF-EQUILIBRIUM BEHAVIOUR AND IDENTIFICATION IN BLENDED MULTIPLE MODEL STRUCTURES Roderick Mrray-Smith Dept. of Compting Science, Glasgow Uniersity, Glasgow, Scotland. rod@dcs.gla.ac.k

More information

Math 144 Activity #10 Applications of Vectors

Math 144 Activity #10 Applications of Vectors 144 p 1 Math 144 Actiity #10 Applications of Vectors In the last actiity, yo were introdced to ectors. In this actiity yo will look at some of the applications of ectors. Let the position ector = a, b

More information

Control Performance Monitoring of State-Dependent Nonlinear Processes

Control Performance Monitoring of State-Dependent Nonlinear Processes Control Performance Monitoring of State-Dependent Nonlinear Processes Lis F. Recalde*, Hong Ye Wind Energy and Control Centre, Department of Electronic and Electrical Engineering, University of Strathclyde,

More information

The Open Civil Engineering Journal

The Open Civil Engineering Journal Send Orders for Reprints to reprints@benthamscience.ae 564 The Open Ciil Engineering Jornal, 16, 1, 564-57 The Open Ciil Engineering Jornal Content list aailable at: www.benthamopen.com/tociej/ DOI: 1.174/187414951611564

More information

Creating a Sliding Mode in a Motion Control System by Adopting a Dynamic Defuzzification Strategy in an Adaptive Neuro Fuzzy Inference System

Creating a Sliding Mode in a Motion Control System by Adopting a Dynamic Defuzzification Strategy in an Adaptive Neuro Fuzzy Inference System Creating a Sliding Mode in a Motion Control System by Adopting a Dynamic Defzzification Strategy in an Adaptive Nero Fzzy Inference System M. Onder Efe Bogazici University, Electrical and Electronic Engineering

More information

Simulation Based Analysis of Two Different Control Strategies for PMSM

Simulation Based Analysis of Two Different Control Strategies for PMSM International Jornal of Engineering Trends and Technology (IJETT) - Volme4Isse4- April 23 Simlation Based Analysis of Two Different Control Strategies for PMSM Lindita Dhamo #, Aida Spahi #2 # Department

More information

Essentials of optimal control theory in ECON 4140

Essentials of optimal control theory in ECON 4140 Essentials of optimal control theory in ECON 4140 Things yo need to know (and a detail yo need not care abot). A few words abot dynamic optimization in general. Dynamic optimization can be thoght of as

More information

Minimizing Intra-Edge Crossings in Wiring Diagrams and Public Transportation Maps

Minimizing Intra-Edge Crossings in Wiring Diagrams and Public Transportation Maps Minimizing Intra-Edge Crossings in Wiring Diagrams and Pblic Transportation Maps Marc Benkert 1, Martin Nöllenbrg 1, Takeaki Uno 2, and Alexander Wolff 1 1 Department of Compter Science, Karlsrhe Uniersity,

More information

Study on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom

Study on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom EPJ Web of Conferences 80, 0034 (08) EFM 07 Stdy on the implsive pressre of tank oscillating by force towards mltiple degrees of freedom Shigeyki Hibi,* The ational Defense Academy, Department of Mechanical

More information

Cosmic Microwave Background Radiation. Carl W. Akerlof April 7, 2013

Cosmic Microwave Background Radiation. Carl W. Akerlof April 7, 2013 Cosmic Microwave Backgrond Radiation Carl W. Akerlof April 7, 013 Notes: Dry ice sblimation temperatre: Isopropyl alcohol freezing point: LNA operating voltage: 194.65 K 184.65 K 18.0 v he terrestrial

More information

Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introduction The transmission line equations are given by

Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introduction The transmission line equations are given by Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introdction The transmission line eqations are given by, I z, t V z t l z t I z, t V z, t c z t (1) (2) Where, c is the per-nit-length

More information

An Investigation into Estimating Type B Degrees of Freedom

An Investigation into Estimating Type B Degrees of Freedom An Investigation into Estimating Type B Degrees of H. Castrp President, Integrated Sciences Grop Jne, 00 Backgrond The degrees of freedom associated with an ncertainty estimate qantifies the amont of information

More information

Complex Tire-Ground Interaction Simulation: Recent Developments Of An Advanced Shell Theory Based Tire Model

Complex Tire-Ground Interaction Simulation: Recent Developments Of An Advanced Shell Theory Based Tire Model . ozdog and W. W. Olson Complex Tire-Grond Interaction Simlation: ecent eelopments Of n danced Shell Theory ased Tire odel EFEECE: ozdog. and Olson W. W. Complex Tire-Grond Interaction Simlation: ecent

More information

System identification of buildings equipped with closed-loop control devices

System identification of buildings equipped with closed-loop control devices System identification of bildings eqipped with closed-loop control devices Akira Mita a, Masako Kamibayashi b a Keio University, 3-14-1 Hiyoshi, Kohok-k, Yokohama 223-8522, Japan b East Japan Railway Company

More information

Change of Variables. (f T) JT. f = U

Change of Variables. (f T) JT. f = U Change of Variables 4-5-8 The change of ariables formla for mltiple integrals is like -sbstittion for single-ariable integrals. I ll gie the general change of ariables formla first, and consider specific

More information

Pulses on a Struck String

Pulses on a Struck String 8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a

More information

Graphs and Networks Lecture 5. PageRank. Lecturer: Daniel A. Spielman September 20, 2007

Graphs and Networks Lecture 5. PageRank. Lecturer: Daniel A. Spielman September 20, 2007 Graphs and Networks Lectre 5 PageRank Lectrer: Daniel A. Spielman September 20, 2007 5.1 Intro to PageRank PageRank, the algorithm reportedly sed by Google, assigns a nmerical rank to eery web page. More

More information

Primary dependent variable is fluid velocity vector V = V ( r ); where r is the position vector

Primary dependent variable is fluid velocity vector V = V ( r ); where r is the position vector Chapter 4: Flids Kinematics 4. Velocit and Description Methods Primar dependent ariable is flid elocit ector V V ( r ); where r is the position ector If V is known then pressre and forces can be determined

More information

2 Faculty of Mechanics and Mathematics, Moscow State University.

2 Faculty of Mechanics and Mathematics, Moscow State University. th World IMACS / MODSIM Congress, Cairns, Astralia 3-7 Jl 9 http://mssanz.org.a/modsim9 Nmerical eamination of competitie and predator behaior for the Lotka-Volterra eqations with diffsion based on the

More information

ENGINEERING COUNCIL DYNAMICS OF MECHANICAL SYSTEMS D225 TUTORIAL 2 LINEAR IMPULSE AND MOMENTUM

ENGINEERING COUNCIL DYNAMICS OF MECHANICAL SYSTEMS D225 TUTORIAL 2 LINEAR IMPULSE AND MOMENTUM ENGINEERING COUNCIL DYNAMICS OF MECHANICAL SYSTEMS D5 TUTORIAL LINEAR IMPULSE AND MOMENTUM On copletion of this ttorial yo shold be able to do the following. State Newton s laws of otion. Define linear

More information

Simplified Identification Scheme for Structures on a Flexible Base

Simplified Identification Scheme for Structures on a Flexible Base Simplified Identification Scheme for Strctres on a Flexible Base L.M. Star California State University, Long Beach G. Mylonais University of Patras, Greece J.P. Stewart University of California, Los Angeles

More information

III. Demonstration of a seismometer response with amplitude and phase responses at:

III. Demonstration of a seismometer response with amplitude and phase responses at: GG5330, Spring semester 006 Assignment #1, Seismometry and Grond Motions De 30 Janary 006. 1. Calibration Of A Seismometer Using Java: A really nifty se of Java is now available for demonstrating the seismic

More information

NUCLEATION AND SPINODAL DECOMPOSITION IN TERNARY-COMPONENT ALLOYS

NUCLEATION AND SPINODAL DECOMPOSITION IN TERNARY-COMPONENT ALLOYS NUCLEATION AND SPINODAL DECOMPOSITION IN TERNARY-COMPONENT ALLOYS COLLEEN ACKERMANN AND WILL HARDESTY Abstract. The Cahn-Morral System has often been sed to model the dynamics of phase separation in mlti-component

More information

A Fully-Neoclassical Finite-Orbit-Width Version. of the CQL3D Fokker-Planck code

A Fully-Neoclassical Finite-Orbit-Width Version. of the CQL3D Fokker-Planck code A Flly-Neoclassical Finite-Orbit-Width Version of the CQL3 Fokker-Planck code CompX eport: CompX-6- Jly, 6 Y. V. Petrov and. W. Harvey CompX, el Mar, CA 94, USA A Flly-Neoclassical Finite-Orbit-Width Version

More information

Simulation investigation of the Z-source NPC inverter

Simulation investigation of the Z-source NPC inverter octoral school of energy- and geo-technology Janary 5 20, 2007. Kressaare, Estonia Simlation investigation of the Z-sorce NPC inverter Ryszard Strzelecki, Natalia Strzelecka Gdynia Maritime University,

More information

Model Explaining the Gravitational Force

Model Explaining the Gravitational Force Model Explaining the Gravitational Force Clade Mercier eng., Agst nd, 015 Rev. October 17 th, 015 clade.mercier@gctda.com The niversal gravitation eqation of Newton is widely sed in the calclations in

More information

Discussion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli

Discussion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli 1 Introdction Discssion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli Søren Johansen Department of Economics, University of Copenhagen and CREATES,

More information

Analysis of Enthalpy Approximation for Compressed Liquid Water

Analysis of Enthalpy Approximation for Compressed Liquid Water Analysis of Entalpy Approximation for Compressed Liqid Water Milioje M. Kostic e-mail: kostic@ni.ed Nortern Illinois Uniersity, DeKalb, IL 60115-2854 It is cstom to approximate solid and liqid termodynamic

More information

Concept of Stress at a Point

Concept of Stress at a Point Washkeic College of Engineering Section : STRONG FORMULATION Concept of Stress at a Point Consider a point ithin an arbitraril loaded deformable bod Define Normal Stress Shear Stress lim A Fn A lim A FS

More information

Flood flow at the confluence of compound river channels

Flood flow at the confluence of compound river channels Rier Basin Management VIII 37 Flood flow at the conflence of compond rier channels T. Ishikawa 1, R. Akoh 1 & N. Arai 2 1 Department of Enironmental Science and Technology, Tokyo Institte of Technology,

More information

3.3 Operations With Vectors, Linear Combinations

3.3 Operations With Vectors, Linear Combinations Operations With Vectors, Linear Combinations Performance Criteria: (d) Mltiply ectors by scalars and add ectors, algebraically Find linear combinations of ectors algebraically (e) Illstrate the parallelogram

More information

Sareban: Evaluation of Three Common Algorithms for Structure Active Control

Sareban: Evaluation of Three Common Algorithms for Structure Active Control Engineering, Technology & Applied Science Research Vol. 7, No. 3, 2017, 1638-1646 1638 Evalation of Three Common Algorithms for Strctre Active Control Mohammad Sareban Department of Civil Engineering Shahrood

More information

ON THE PERFORMANCE OF LOW

ON THE PERFORMANCE OF LOW Monografías Matemáticas García de Galdeano, 77 86 (6) ON THE PERFORMANCE OF LOW STORAGE ADDITIVE RUNGE-KUTTA METHODS Inmaclada Higeras and Teo Roldán Abstract. Gien a differential system that inoles terms

More information

Chapter 5 Computational Methods for Mixed Models

Chapter 5 Computational Methods for Mixed Models Chapter 5 Comptational Methods for Mixed Models In this chapter we describe some of the details of the comptational methods for fitting linear mixed models, as implemented in the lme4 package, and the

More information

Linear Strain Triangle and other types of 2D elements. By S. Ziaei Rad

Linear Strain Triangle and other types of 2D elements. By S. Ziaei Rad Linear Strain Triangle and other tpes o D elements B S. Ziaei Rad Linear Strain Triangle (LST or T6 This element is also called qadratic trianglar element. Qadratic Trianglar Element Linear Strain Triangle

More information

Elements of Coordinate System Transformations

Elements of Coordinate System Transformations B Elements of Coordinate System Transformations Coordinate system transformation is a powerfl tool for solving many geometrical and kinematic problems that pertain to the design of gear ctting tools and

More information

A Model-Free Adaptive Control of Pulsed GTAW

A Model-Free Adaptive Control of Pulsed GTAW A Model-Free Adaptive Control of Plsed GTAW F.L. Lv 1, S.B. Chen 1, and S.W. Dai 1 Institte of Welding Technology, Shanghai Jiao Tong University, Shanghai 00030, P.R. China Department of Atomatic Control,

More information

arxiv: v1 [cs.sy] 22 Nov 2018

arxiv: v1 [cs.sy] 22 Nov 2018 AAS 18-253 ROBUST MYOPIC CONTROL FOR SYSTEMS WITH IMPERFECT OBSERVATIONS Dantong Ge, Melkior Ornik, and Ufk Topc arxiv:1811.09000v1 [cs.sy] 22 Nov 2018 INTRODUCTION Control of systems operating in nexplored

More information

THE DIFFERENTIAL GEOMETRY OF REGULAR CURVES ON A REGULAR TIME-LIKE SURFACE

THE DIFFERENTIAL GEOMETRY OF REGULAR CURVES ON A REGULAR TIME-LIKE SURFACE Dynamic Systems and Applications 24 2015 349-360 TH DIFFRNTIAL OMTRY OF RULAR CURVS ON A RULAR TIM-LIK SURFAC MIN OZYILMAZ AND YUSUF YAYLI Department of Mathematics ge Uniersity Bornoa Izmir 35100 Trkey

More information

Department of Industrial Engineering Statistical Quality Control presented by Dr. Eng. Abed Schokry

Department of Industrial Engineering Statistical Quality Control presented by Dr. Eng. Abed Schokry Department of Indstrial Engineering Statistical Qality Control presented by Dr. Eng. Abed Schokry Department of Indstrial Engineering Statistical Qality Control C and U Chart presented by Dr. Eng. Abed

More information

Effects of MHD Laminar Flow Between a Fixed Impermeable Disk and a Porous Rotating Disk

Effects of MHD Laminar Flow Between a Fixed Impermeable Disk and a Porous Rotating Disk Effects of MHD Laminar Flow Between a Fixed Impermeable Disk and a Poros otating Disk Hemant Poonia * * Asstt. Prof., Deptt. of Math, Stat & Physics, CCSHAU, Hisar-54.. C. Chadhary etd. Professor, Deptt.

More information

The Faraday Induction Law and Field Transformations in Special Relativity

The Faraday Induction Law and Field Transformations in Special Relativity Apeiron, ol. 10, No., April 003 118 The Farada Indction Law and Field Transformations in Special Relatiit Aleander L. Kholmetskii Department of Phsics, elars State Uniersit, 4, F. Skorina Aene, 0080 Minsk

More information

SUBJECT:ENGINEERING MATHEMATICS-I SUBJECT CODE :SMT1101 UNIT III FUNCTIONS OF SEVERAL VARIABLES. Jacobians

SUBJECT:ENGINEERING MATHEMATICS-I SUBJECT CODE :SMT1101 UNIT III FUNCTIONS OF SEVERAL VARIABLES. Jacobians SUBJECT:ENGINEERING MATHEMATICS-I SUBJECT CODE :SMT0 UNIT III FUNCTIONS OF SEVERAL VARIABLES Jacobians Changing ariable is something e come across er oten in Integration There are man reasons or changing

More information

EVALUATION OF GROUND STRAIN FROM IN SITU DYNAMIC RESPONSE

EVALUATION OF GROUND STRAIN FROM IN SITU DYNAMIC RESPONSE 13 th World Conference on Earthqake Engineering Vancover, B.C., Canada Agst 1-6, 2004 Paper No. 3099 EVALUATION OF GROUND STRAIN FROM IN SITU DYNAMIC RESPONSE Ellen M. RATHJE 1, Wen-Jong CHANG 2, Kenneth

More information

FRTN10 Exercise 12. Synthesis by Convex Optimization

FRTN10 Exercise 12. Synthesis by Convex Optimization FRTN Exercise 2. 2. We want to design a controller C for the stable SISO process P as shown in Figre 2. sing the Yola parametrization and convex optimization. To do this, the control loop mst first be

More information

Discussion Papers Department of Economics University of Copenhagen

Discussion Papers Department of Economics University of Copenhagen Discssion Papers Department of Economics University of Copenhagen No. 10-06 Discssion of The Forward Search: Theory and Data Analysis, by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli Søren Johansen,

More information

New Regularized Algorithms for Transductive Learning

New Regularized Algorithms for Transductive Learning New Reglarized Algorithms for Transdctie Learning Partha Pratim Talkdar and Koby Crammer Compter & Information Science Department Uniersity of Pennsylania Philadelphia, PA 19104 {partha,crammer}@cis.penn.ed

More information

4 Exact laminar boundary layer solutions

4 Exact laminar boundary layer solutions 4 Eact laminar bondary layer soltions 4.1 Bondary layer on a flat plate (Blasis 1908 In Sec. 3, we derived the bondary layer eqations for 2D incompressible flow of constant viscosity past a weakly crved

More information

Model Discrimination of Polynomial Systems via Stochastic Inputs

Model Discrimination of Polynomial Systems via Stochastic Inputs Model Discrimination of Polynomial Systems via Stochastic Inpts D. Georgiev and E. Klavins Abstract Systems biologists are often faced with competing models for a given experimental system. Unfortnately,

More information

DEPTH CONTROL OF THE INFANTE AUV USING GAIN-SCHEDULED REDUCED-ORDER OUTPUT FEEDBACK 1. C. Silvestre A. Pascoal

DEPTH CONTROL OF THE INFANTE AUV USING GAIN-SCHEDULED REDUCED-ORDER OUTPUT FEEDBACK 1. C. Silvestre A. Pascoal DEPTH CONTROL OF THE INFANTE AUV USING GAIN-SCHEDULED REDUCED-ORDER OUTPUT FEEDBACK 1 C. Silvestre A. Pascoal Institto Sperior Técnico Institte for Systems and Robotics Torre Norte - Piso 8, Av. Rovisco

More information

Linear System Theory (Fall 2011): Homework 1. Solutions

Linear System Theory (Fall 2011): Homework 1. Solutions Linear System Theory (Fall 20): Homework Soltions De Sep. 29, 20 Exercise (C.T. Chen: Ex.3-8). Consider a linear system with inpt and otpt y. Three experiments are performed on this system sing the inpts

More information

Frequency Estimation, Multiple Stationary Nonsinusoidal Resonances With Trend 1

Frequency Estimation, Multiple Stationary Nonsinusoidal Resonances With Trend 1 Freqency Estimation, Mltiple Stationary Nonsinsoidal Resonances With Trend 1 G. Larry Bretthorst Department of Chemistry, Washington University, St. Lois MO 6313 Abstract. In this paper, we address the

More information