OPTIMAL CONTROL OF PARALLEL HYBRID ELECTRIC VEHICLES BASED ON THEORY OF SWITCHED SYSTEM

Size: px
Start display at page:

Download "OPTIMAL CONTROL OF PARALLEL HYBRID ELECTRIC VEHICLES BASED ON THEORY OF SWITCHED SYSTEM"

Transcription

1 274 Asan Journal of Control, Vol. 8, No. 3, pp , September 26 -Bref Paper- OPTIMAL CONTROL OF PARALLEL HYBRID ELECTRIC VEHICLES BASED ON THEORY OF SWITCHED SYSTEM Lcun Fang and Shyn Qn ABSTRACT The focus of ths paper s the control strategy used to control general parallel hybrd electrc vehcles (HEV). The torque splt control problem of HEV s formulated as the optmal control of a swtched system. A model-based strategy for fuel-optmal control s presented. The optmal control problem of such a swtched system s formulated as a two-stage optmzaton problem. Dynamc programmng s utlzed to determne the optmal control acton that mnmzes the cost functon. Smulated results ndcate that ths method s effectve. KeyWords: Hybrd electrc vehcles (HEV), optmal control, dynamc programmng, swtched system. I. INTRODUCTION The most challengng goals facng the automotve ndustry are ncreasng fuel economy and reducng emssons. There are a seres of constrants that are mposed not only by our socety but also by cooperatve agreements and legslatve efforts [1]. HEVs are begnnng to demonstrate ther capablty to satsfy these requrements, and they wll become a vable alternatve to conventonal vehcles n the future [2,3]. The hybrd powertran s an ntegrated system that may consst of the followng components: an nternal combuston engne (IC engne), a battery pack, and an electrc machne (EM) whch can be utlzed as a tracton motor or generator. In such a system, each sub-system s also a complex system whch has ts own functonalty and desred performance. Moreover, all of the sub-systems need to be coordnated n an optmal manner to acheve the desred objectves. Consequently, an ntegrated vehclelevel s requred to accomplsh such tasks [4]. Recent reports on HEVS have placed great emphass on the control strategy, and several approaches to t have been Manuscrpt receved July 18, 25; revsed January 25, 26, accepted May 15, 26. The authors are wth the School of Automaton Scence and Electrcal Engneerng, Behang Unversty, Bejng 183, Chna (e-mal: lcun_fang@buaa.edu.cn, qsy@buaa.edu.cn). proposed. The power management strateges can be roughly classfed nto three types: statc optmzaton methods [5,6], global optmzaton methods based on dynamc programmng [7,8], and rule-based control approaches [9,1]. HEV systems smultaneously exhbt several knds of dynamc behavor, such as contnuous-tme dynamcs, dscrete-tme dynamcs, jump phenomena, swtchng and the lke. Therefore, an HEV s a typcal hybrd dynamcal system, and many challenges exst n the desgn of a vehcle system for an HEV [4]. The focus of ths paper s the control strategy used to control general parallel HEVs. The remander of ths paper s organzed as follows. Secton 2 wll descrbe the herarchcal archtecture of the HEV system. The dfferent powertran components models used for smulaton and control are presented n secton 3. The optmal control problem s solved n secton 4. In secton 5, a case study nvestgated through smulaton s presented, and the results obtaned n optmzng the control of fuel consumpton based on the optmal strategy proposed n ths paper are gven. Conclusons are drawn n secton 6. II. CONTROL SYSTEM ARCHITECTURE An HEV s a typcal hybrd dynamcal system [4]. Most hybrd dynamcal systems have been desgned usng the concept of centralzed and dstrbuted control. A

2 L. Fang and S. Qn: Optmal Control of Parallel Hybrd Electrc Vehcles Based on Theory of Swtched Systems 275 two-level herarchcal control archtecture [11] s used n the synthess of the control law for the hybrd powertran shown n Fg. 1. The supervsory powertran (SPC) represents a hgh-level vehcle control system that can coordnate the overall powertran to satsfy certan performance ndces. Based on the drver s demand and operatng parameter feedback sgnals of the vehcle and ts components, the SPC must determne the desred propellng torque dstrbuton between the IC engne and EM, or the desred brake torque dstrbuton between the EM and frcton brake system, and the gear poston; therefore, the man functon of the SPC s to determne the desred output of varous sub-systems accordng to drver s commands, system parameter feedback sgnals and a pre-programmed control strategy. Dfferent control strateges wll result n dfferent fuel economy and emsson characterstcs. The focus of ths paper s the control strategy used to control general parallel HEVs. In order to narrow our research scope, only a double shaft parallel powertran s consdered. In ths confguraton, the engne power and electrc motor power can provde torque to the wheels separately through the transmsson. It should be ponted out that the strateges for dealng wth ths stuaton can also be appled to other powertran arrangements by makng approprate amendments. Fgure 2 shows a schematc dagram of the double shaft parallel HEV system consdered n ths paper. CAN bus Engne Drver Supervsory powertran control Motor/Battery Trans. Brake Engne Clutch Motor Trans. Vehcle The followng components are utlzed n the confguraton: spark gnton drect njecton engne: 5kW; permanent magnet motor: 1kW; NMh battery: 2kW, 2kWh; 5speed automated manual transmsson. III. SYSTEM MODEL A qualty dynamc model s paramount for the development of an effectve control strategy, especally when dealng wth a complex coordnated control problem such as a parallel HEV. An HEV powertran dynamc model descrbng key system components can be a sgnfcant ad to the development of the control law and provdes nsght nto the domnant dynamc effects n the control law synthess process. To ths end, the components of a parallel HEV are modeled. However, developng a vehcle s the key objectve of ths paper, so smplfed models are presented. Dynamc models of varous prmary HEV components, such as the ICE, motor, battery, transmsson and vehcle, are descrbed below. 3.1 Engne To avod the dependence on the avalablty of specfc effcency maps for the IC engne and the electrc motor whle optmzng and formulatng the HEV optmal control strategy, a unversal representaton, the Wllans lnes model of an IC engne and motor, s utlzed [12]. A bref summary of the parametrzaton s gven n the followng paragraphs. The scalablty of an IC engne s acheved through the use of two engne szng parameters, the dsplacement V d and the stroke S. Here, some concepts and correspondng varables, such as the mean pston speed c m, the avalable mean effectve pressure p ma and the mean effectve pressure p me, are ntroduced as follows: Sωe c m =, (1) π Battery Fg. 1. Herarchcal control archtecture of an HEV. P ma 4πQLHV m = V ω d e f, (2) Wheel Wheel Mechancal couplng Reducton gears Gear box Motor Clutch Fg. 2. Block dagram of the parallel HEV. Battery ICE P me 4πT = e, (3) V d where ω e s the ICE angular speed, T e s the ICE effectve torque, m f s the fuel mass flow rate, Q LHV s the fuel low heatng value, and e s the nner ablty to transform chemcal to mechancal energy. The engne s effcency can be approxmated by an affne relaton between the torque and fuel mass avalable n one cycle:

3 276 Asan Journal of Control, Vol. 8, No. 3, September 26 m f Te = e QLHV T ω e loss, (4) where T loss represents the nner losses. The resultng affne functon shows the domnant effect observed n all ICEs,.e., the dramatc decrease n effcency as the torque output decreases. Based on Eqs. (1) ~ (3) and some necessary substtuton n Eq. (4), a dmensonless defnton of engne effcency η can be obtaned: p me η =, (5) pma where Pme = epma Pml, (6) and p ml denotes the nner loss (see [12] for more detals). 3.2 Motor In ths secton, the same concepts that were ntroduced above for an IC engne are developed for an EM. Let p em be the avalable electrcal power. Then, for each tme nstant, the followng equaton defnes the effcency of the EM: ω T = η() p = η() V, (7) em em, e em where ω em and T em, e are the EM angular speed and effectve torque, respectvely; V and are the voltage and current, respectvely. Agan, as a frst approach, the motor s effcency can be approxmated by assumng that an affne relatonshp between the effectve torque and nput energy exsts: T = e V T em, e em em,loss ωem = e T T, (8) em em, a em,loss where e em, T em, a and T em, loss are the nner ablty to transform electrcal energy nto mechancal energy, based on the avalable mean effectve torque of the EM and nner loss (see [12] for more detals). 3.3 Battery The practcal mplementaton of ths control strategy s dependent on the avalablty of an on-lne SOC (State Of Charge) estmator. The battery performance, such as the voltage, current and effcency, from a purely electrc vewpont, s the outcome of thermally dependent electrochemcal processes that are qute complcated. A model based upon open crcut voltage s used n ths paper. Further detals about ths method can be found n [13]. 3.4 Transmsson Two man relatons descrbe the parallel double-shaft arrangement of an HEV. The frst one s the speed relaton: ω e( k) ω em( k) ω w( k) = =. (9) Rk ( ) ρ The second one s the torque relaton: Tw( k) = R( k) Te( k) η gb+ ρ Tm( k)η red, (1) where ω w s the wheel angular speed, R(k) s the gear rato, ρ s the reducton gear rato between the EM and the wheels, T w (k) s the torque of the wheels, and T e (k) and T m (k) are the output torque from the engne and the motor, respectvely. η gb denotes the effcency of the gear box, and η red s that of the reducton gears. 3.5 Vehcle The modelng of the vehcle dynamcs s derved from the basc equaton for sold-body motons (Newton s Second Law): 1 2 F = MgCr + Dar CD Af V 2 dv + M + Mgsn(θ), (11) dt where F s the force requred at the vehcle wheels to reach a certan acceleraton at speed V ; Mg, C r, D ar, C D, A f, M, and θ are the mass of the vehcle, the local acceleraton of gravty, the coeffcent of the rollng resstance between the tres and the road surface, the ar densty, the coeffcent of the aerodynamc drag for the vehcle n the travelng drecton, the frontal area of the vehcle, the nertal mass of the vehcle, ncludng the rotatonal nerta contrbuton, and the gradent of the road, respectvely. The vehcle acceleraton a and speed V are gven by 1 2 η gb RkT ( ) / r MgCr Dar CD Af V mgsn(θ) a = 2, M (12) V = adt, (13) where r s the radus of the wheels. IV. OPTIMAL CONTROL 4.1 Problem Formulaton In the operaton process of a powertran, the transmsson determnes the gear shftng sequence and selects the approprate transmsson gear based on the transmsson output speed, the poston of the acceleraton pedal and brakng pedal, the current gear poston and the clutch state. The powertran system can be formulated as a swtched system. As s well-known, a dscrete-tme swtched system can be descrbed as a set of pece-wse

4 L. Fang and S. Qn: Optmal Control of Parallel Hybrd Electrc Vehcles Based on Theory of Swtched Systems 277 dfference equatons [14]: x( k + 1) = f( x( k), u( k), m( k)) f1 ( x( k), u( k)), f m( k) = 1 = (14) fi ( xk ( ), uk ( )), f mk ( ) = I, where x(k) s a state vector of the system, u(k) s a vector of the contnuous-tme control varables, and m(k) {1, 2, I} s a fnte nteger set n whch every element corresponds to a knd of operatng mode wth event attrbutes. An HEV can be consdered as a specal case of ths model. For a vehcle control system, the contnuous-tme state vector s (V, SOC); the control varables u(k) are T e (k) and T m (k); m(k) s the gear poston, whch s a dscrete-event state varable. If the state varables and control varables are nserted nto Eq. (14), a set of pece-wse dfference equatons for an HEV can be obtaned. It should be ponted out that the motor torque becomes a dependent varable nstead of a control varable due to the equalty constrant on the drvelne torque: Twh _req = Te ( k) + Tm ( k), (15) where T wh_req s the requred torque of the vehcle. For a vehcle control system, once a drvng cycle s gven, the wheel speed profle w wh_req s known, and the torque T wh_req requred to follow the speed profle can be determned n each tme step by nversely solvng the dynamc equatons (11) ~ (13). The practcal wheel speed w wh can be detected by a speed sensor. Problem 1. The power control of an HEV can be formulated as an optmal control problem, n whch the optmzaton goal s to fnd the control nput T e (k) and swtchng sequence of gear shftng g(k) needed to mnmze a cost functon n a gven drvng cycle: N 1 J = ψ( x( N)) + L( x( k), Te ( k), m( k)) k = N 1 + M k( x, x ' ), (16) k = where N s the duraton of the drvng cycle, ψ (x (N)) s the cost functon assocated wth the error n the termnal SOC, L(x (k), T e (k), m (k)) s the nstantaneous fuel consumpton rate, and x, x represent dfferent gear postons. M k (x, x ) s the cost functon of mode swtchng between x and x ; n order to smplfy the problem, the cost functon of mode swtchng s assumed to be a constant set and to take values from ths set. In order to ensure drvablty of the vehcle and safe/smooth operaton of each subsystem, t s necessary to satsfy the followng constrants: T ( w ( k)) T T ( w ( k)), (17) e_mn e e e_max e where T e_mn (w e (k)) and T e_max (w e (k)) are the mnmal and maxmal output torque of the ICE at speed w e (k); w e_mn w e ( k) w e_max, (18) where w e_mn and w e_max are the mnmal and maxmal speeds of the ICE; SOC SOC( k) SOC, (19) mn where SOC mn and SOC max are the lower lmt and upper lmt of the SOC of the battery; m_max max Tm_mn( wm( k), SOC( k)) Tm( k) T ( w ( k),soc( k)), m (2) where T m_mn (w m (k), SOC(k)) and T m_max (w m (k), SOC(k)) are the mnmal and maxmal output torque of the EM at speed w m (k) for the battery SOC(k). 4.2 Soluton For the optmal control problem of a swtched system, the followng condtons must be satsfed [15]. Condton 1. An optmal soluton (σ *, u * ) exsts, where σ * s the optmal swtchng sequence and u * s the optmal nput control. Condton 2. For any gven swtchng sequence σ, there exsts a correspondng u * = u * (σ) such that J σ (u) = J(σ, u) s mnmzed. If the above condtons are satsfed, then, the followng equaton holds: mn J ( σ, u) = mn mn J ( σ, u). (21) σ, u U σ u U Therefore, the optmal control problem of a swtched system can be formulated as a two-stage optmzaton problem. Stage 1. (a) Fx the total number of swtches to be K, fx the order of the actve subsystem, let the mnmum value of J wth respect to u be a functon of the K swtchng nstants,.e., J 1 = J 1 (t 1, t 2,, t k ) for K, and then fnd J 1. (b) Mnmze J 1 wth respect to t 1, t 2,, t k. Stage 2. (a) Vary the order of the actve subsystem to fnd an optmal soluton under K swtches. (b) Vary the number of swtches K to fnd an optmal soluton. Dynamc programmng s a powerful tool for solvng general dynamc optmzaton problems. The man advantage s that t can easly deal wth the constrants and nonlnearty of such problems whle obtanng a globally optmal soluton. The dynamc programmng technque s based on the Bellman s Prncple of Optmalty, whch

5 278 Asan Journal of Control, Vol. 8, No. 3, September 26 states that the optmal polcy can be obtaned f we frst solve a one stage problem nvolvng only the last stage and then gradually extend t to the frst stage. As a matter of fact, dynamc programmng has been used to solve the optmal control of a swtched system [16]. By denotng Σ [k, N] as a subset of Σ [k, N] whch contans all σ Σ [k, N], startng wth subsystem at nstant k, one can compute the cost functons V (x(k), k), I by solvng Eq. (22), where, ' represent a subsystem and V (x (k), k) s the mnmum value of J f the system starts at tme k wth state x(k) and subsystem : V ( x( k), k) = mn mn J(σ, u) ' j {, } σ u U [ kn, ] ψ(( N)) f k = N mn{ Lxk ( ( ), uk ( )) = uk ( ) (22) j + mn V ( x( k + 1), k + 1)} f k < N. Consequently, the optmal soluton (σ, u) of a dscretetme swtched system can be constructed by solvng V (x(k), k) backwards and fnally fndng (x(), ). In a double-shaft parallel hybrd confguraton, condtons 1 and 2 are satsfed, and the gear rato of the transmsson only affects the torque contrbuton of the IC engne; consequently, we can dvde the development of the overall optmal control strategy nto two stages: (1) mplement a gear shftng strategy, whch selects the swtchng sequence of gears n a dscrete set Σ to optmze the operaton of the engne; (2) mplement a power splt strategy, whch defnes the best power splt between the ICE and EM. Wth the above two-stage optmzaton process, the optmal control polcy (T e (k), g(k)) can be constructed. 4.3 Problem Smplfcaton The bg challenge n utlzng the above algorthm to fnd the optmal feedback control polcy s overcomng the dffculty posed by the huge computaton burden. It s, thus, necessary to adopt some approaches to accelerate the computaton speed. In ths paper, we smplfy the above two steps as descrbed below. (1) Shftng Sequence of Gear Poston Determnng the gear shft strategy s crucal to mprovng the fuel economy of hybrd electrc vehcles [17]. The gear shftng sequence of an automatc transmsson can be modeled as a dscrete-tme dynamc system: M, mk ( ) + gk ( ) M mk ( + 1) = 1, mk ( ) + gk ( ) 1 mk ( ) + gk ( ), otherwse, (23) where M s the number of gear postons, and g(k) s constraned to take on values from among 1,, and 1, representng downshft, hold and upshft, respectvely. From the speed profle of a gven drvng cycle, the optmal gear poston can be constructed through statc optmzaton [9], the total swtchng numbers of gear poston can be fxed and the shftng sequence can be determned; consequently, the problem can be transformed nto an optmal control problem for a swtched system wth a fxed swtch number and fxed swtchng sequence. (2) The Optmal Torque Input The state and control varables are frst quantzed nto a fnte grd. For a gven drvng cycle, the vehcle model can be replaced by a fnte set of operatng ponts parameterzed by the requred torque and speed. Pre-computed look-up tables can be constructed for recordng the next states and cost functon as a functon of the quantzed states, control nputs and operatng ponts. Once these tables are bult, they can be used to solve the problem n a very effcent manner. The dynamc programmng procedure produces an * optmal, tme-varyng, state-feedback control polcy ( Te ( k ), g * ()) k for a drvng cycle. Ths control polcy can then be used as a state feedback. V. CASE STUDY To verfy the control polcy determned by the dynamc programmng algorthm, a case was studed and wll be dscussed n ths secton. The UDDS (Urban Dynamometer Drvng Schedule) drvng cycle s a standard test drvng cycle utlzed by the Envronmental Protecton Agency to certfy the fuel economy and emsson performance of vehcles whch are drven n urban areas. The optmal control polcy was found through DP (Dynamc Programmng) n a UDDS drvng cycle. In order to evaluate the fnal fuel economy, the optmal control law was appled to the full-order parallel HEV model. The termnal SOC constrant was selected as.65, and the ntal SOC n the smulaton was chosen as.65 for the purpose of calculatng the fuel economy. The dynamc trajectores of the vehcle under the optmal control polcy for the UDDS drvng cycle are shown n Fg. 3 and Fg. 4. Speed(km/h) Dfference of Tme(s) Fg. 3. Dfference between the requred and acheved speeds.

6 L. Fang and S. Qn: Optmal Control of Parallel Hybrd Electrc Vehcles Based on Theory of Swtched Systems 279 Speed(km/h) SOC Overall rato Fuel(L) B C 5.6 D.4.2 Tme(s) Fg. 4. Dynamc trajectores of the actual speed profle, battery SOC, fuel consumpton and overall rato. The largest dfference between the desred vehcle speed and the actual speed s wthn.4km/h, whch means that the control strategy can trace the objectve speed correctly. The SOC trajectory starts at.65 and ends at around.65 wth a small quantzaton error. The fuel consumpton of the control polcy s.452l, and the average fuel consumpton s 3.8 L/1Km. The results were compared wth the default control strategy n ADVISOR 22, that s, the electrcal peakng energy management strategy. The fuel consumpton of the control polcy s found to be.4695l, and the average fuel consumpton s 3.92 L/1Km, so the optmal control law s effectve n mprovng fuel economy. VI. CONCLUSIONS A The optmal torque splt problem of a parallel HEV has been formulated as the optmal control problem of a swtched system. A control strategy for general parallel HEVs has been developed n ths paper. The proposed optmal torque splt strategy for parallel HEVs, desgned wth the ad of dynamc programmng, can mprove the overall system effcency. A two-stage dynamc programmng algorthm has been used to solve the mnmum fuel optmal control problem for a parallel hybrd electrc vehcle. In order to reduce the computaton complexty, several approaches have been utlzed to accelerate the computng speed. A dynamc optmal soluton to the energy management problem over a drvng cycle has been obtaned. The smulaton results ndcate that mprovement n fuel economy s acheved by optmzng the gear shftng polcy and torque splttng strategy. In general, torque splt algorthms developed through DP are more accurate, so the approach proposed n ths paper can be utlzed to assess other control strateges. REFERENCES 1. Maclean, H.L. and L.B. Lave, Evaluaton Automoble Fuel/Propulson System Technologes, Progr. Energ. Combust. Sc., Vol. 29, pp (23). 2. Btsche, O. and G. Gutmann, System for Hybrd Cars, J. Power Sources, Vol. 127, pp (24). 3. Chau, K.T. and Y.S. Wang, Overvew of Power Management n Hybrd Electrc Vehcles, Energy Convers. Manage., Vol. 43, pp (22). 4. Phllps, A., M. Jankovc, and K. Baley, Vehcle System Controller Desgn for a Hybrd Electrc, Proc. IEEE Int. Conf. Contr. Appl., Alaska, pp (2). 5. Paganell, G., A. Brahma, G. Rzzon, and Y. Guezennec, Control Development for a Hybrd-Electrc Sport-Utlty Vehcles: Strategy, Implementaton and Feld Test Results, Proc. Amer. Contr. Conf., Vol. 6, pp (21). 6. Scarretta, A., M. Back, and L. Guzzella, Optmal Control of Parallel Hybrd Electrc Vehcles, IEEE Trans. Contr. Syst. Technol., Vol. 12, No. 3, pp (24). 7. Paganell, G., G. Ercole, A. Brahma, Y. Guezenne, and G. Rzzon, General Supervsory Control Polcy for the Energy Optmzaton of Charge-Sustanng Hybrd Electrc Vehcles, JSAE Rev., Vol. 22, No. 4, pp (21). 8. Delprat, S., J. Lauber, T.M. Guerra, and J. Rmaux, Control of a Parallel Hybrd Powertran: Optmal Control, IEEE Trans. Veh. Technol., Vol.??, No.??, pp (24). 9. Schouten, N., M. Salman, and N. Kher, Energy Management Strateges for Parallel Hybrd Vehcles Usng Fuzzy Logc, Contr. Eng. Pract., Vol.??, No.??, pp (23). 1. Pusca, R., Y. At-Amrat, A. Berthon, and J.M. Kauffmann, Fuzzy-Logc-Based Control Appled to a Hybrd Electrc Vehcle wth Four Separate Wheel Drves, IEE Proc. Contr. Theory Appl., Vol. 151, Issue 1, pp (24). 11. Ln, C.C., H. Peng, and J.W. Grzzle, Control System Development for an Advanced-Techmology Medum-Duty Hybrd Electrc Truck, SAE Paper, No (23). 12. Rzzon, G., L. Guzzella, and B.M. Baumann, Unfed Modelng of Hybrd Electrc Vehcle Drvetrans,

7 28 Asan Journal of Control, Vol. 8, No. 3, September 26 IEEE/ASME Trans. Mechatron, Vol. 4, No. 3, pp (1999). 13. Gregory, L.P., Hgh-Performance Battery-Pack Power Estmaton Usng a Dynamc Cell Model, IEEE Trans. Veh. Technol., Vol. 5, No. 5, pp (24). 14. Akar, M., U. Mtra, and U. Ozguner, Optmal Control for a Class of Swtchng Systems, Proc. Amer. Contr. Conf., Vol. 4, pp (2). 15. Xu, Xupng, and P.J. Antsakls, A Dynamc Programmng Approach for Optmal Control of Swtched Systems, Proc. 39th IEEE Conf. Decs. Contr., Vol. 2, pp (2). 16. Xu, Xupng and P.J. Antsakls, Optmal Control of Swtched Systems: New Results and Open Problems, Proc. Amer. Contr. Conf., Vol. 4, pp (2). 17. Lee, H.D., S.K. Sul, H.S. Cho, and J.M. Lee, Advanced Gear-Shftng and Clutchng Strategy for a Parallel-Hybrd Vehcle, IEEE Ind. Appl. Mag., Vol. 6, No. 6, pp (2).

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl RECURSIVE SPLINE INTERPOLATION METHOD FOR REAL TIME ENGINE CONTROL APPLICATIONS A. Stotsky Volvo Car Corporaton Engne Desgn and Development Dept. 97542, HA1N, SE- 405 31 Gothenburg Sweden. Emal: astotsky@volvocars.com

More information

Vehicle Propulsion Systems Lecture 2. Energy System Overview. W2M Energy Paths. Evaluating the integral. Mechanical Energy Demand of a Cycle

Vehicle Propulsion Systems Lecture 2. Energy System Overview. W2M Energy Paths. Evaluating the integral. Mechanical Energy Demand of a Cycle Vehcle Propulson Systems Lecture 2 Fuel Consumpton Estmaton & ICE Powertrans Lars Erksson Professor Vehcular Systems Lnköpng Unversty March 21, 2017 2 / 51 3 / 51 W2M Energy Paths Energy System Overvew

More information

Clock-Gating and Its Application to Low Power Design of Sequential Circuits

Clock-Gating and Its Application to Low Power Design of Sequential Circuits Clock-Gatng and Its Applcaton to Low Power Desgn of Sequental Crcuts ng WU Department of Electrcal Engneerng-Systems, Unversty of Southern Calforna Los Angeles, CA 989, USA, Phone: (23)74-448 Massoud PEDRAM

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

EEL 6266 Power System Operation and Control. Chapter 3 Economic Dispatch Using Dynamic Programming

EEL 6266 Power System Operation and Control. Chapter 3 Economic Dispatch Using Dynamic Programming EEL 6266 Power System Operaton and Control Chapter 3 Economc Dspatch Usng Dynamc Programmng Pecewse Lnear Cost Functons Common practce many utltes prefer to represent ther generator cost functons as sngle-

More information

Chapter - 2. Distribution System Power Flow Analysis

Chapter - 2. Distribution System Power Flow Analysis Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

Research Report. Eiichi Ono, Yoshikazu Hattori, Yuji Muragishi. Abstract. Special Issue Estimation and Control of Vehicle Dynamics for Active Safety

Research Report. Eiichi Ono, Yoshikazu Hattori, Yuji Muragishi. Abstract. Special Issue Estimation and Control of Vehicle Dynamics for Active Safety Specal Issue Estmaton and Control of Vehcle Dynamcs for Actve Safety 7 Research Report Estmaton of Tre Frcton Crcle and Vehcle Dynamcs Integrated Control for Four-wheel Dstrbuted Steerng and Four-wheel

More information

Lecture Notes on Linear Regression

Lecture Notes on Linear Regression Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume

More information

Note 10. Modeling and Simulation of Dynamic Systems

Note 10. Modeling and Simulation of Dynamic Systems Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada

More information

EN40: Dynamics and Vibrations. Homework 4: Work, Energy and Linear Momentum Due Friday March 1 st

EN40: Dynamics and Vibrations. Homework 4: Work, Energy and Linear Momentum Due Friday March 1 st EN40: Dynamcs and bratons Homework 4: Work, Energy and Lnear Momentum Due Frday March 1 st School of Engneerng Brown Unversty 1. The fgure (from ths publcaton) shows the energy per unt area requred to

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI 2017 2nd Internatonal Conference on Electrcal and Electroncs: echnques and Applcatons (EEA 2017) ISBN: 978-1-60595-416-5 Study on Actve Mcro-vbraton Isolaton System wth Lnear Motor Actuator Gong-yu PAN,

More information

Modeling of Dynamic Systems

Modeling of Dynamic Systems Modelng of Dynamc Systems Ref: Control System Engneerng Norman Nse : Chapters & 3 Chapter objectves : Revew the Laplace transform Learn how to fnd a mathematcal model, called a transfer functon Learn how

More information

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin Proceedngs of the 007 Wnter Smulaton Conference S G Henderson, B Bller, M-H Hseh, J Shortle, J D Tew, and R R Barton, eds LOW BIAS INTEGRATED PATH ESTIMATORS James M Calvn Department of Computer Scence

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. EE 539 Homeworks Sprng 08 Updated: Tuesday, Aprl 7, 08 DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. For full credt, show all work. Some problems requre hand calculatons.

More information

Army Ants Tunneling for Classical Simulations

Army Ants Tunneling for Classical Simulations Electronc Supplementary Materal (ESI) for Chemcal Scence. Ths journal s The Royal Socety of Chemstry 2014 electronc supplementary nformaton (ESI) for Chemcal Scence Army Ants Tunnelng for Classcal Smulatons

More information

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan Wnter 2008 CS567 Stochastc Lnear/Integer Programmng Guest Lecturer: Xu, Huan Class 2: More Modelng Examples 1 Capacty Expanson Capacty expanson models optmal choces of the tmng and levels of nvestments

More information

Electrical double layer: revisit based on boundary conditions

Electrical double layer: revisit based on boundary conditions Electrcal double layer: revst based on boundary condtons Jong U. Km Department of Electrcal and Computer Engneerng, Texas A&M Unversty College Staton, TX 77843-318, USA Abstract The electrcal double layer

More information

Chapter 6 Electrical Systems and Electromechanical Systems

Chapter 6 Electrical Systems and Electromechanical Systems ME 43 Systems Dynamcs & Control Chapter 6: Electrcal Systems and Electromechancal Systems Chapter 6 Electrcal Systems and Electromechancal Systems 6. INTODUCTION A. Bazoune The majorty of engneerng systems

More information

COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN

COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN Int. J. Chem. Sc.: (4), 04, 645654 ISSN 097768X www.sadgurupublcatons.com COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN R. GOVINDARASU a, R. PARTHIBAN a and P. K. BHABA b* a Department

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

STUDY OF A THREE-AXIS PIEZORESISTIVE ACCELEROMETER WITH UNIFORM AXIAL SENSITIVITIES

STUDY OF A THREE-AXIS PIEZORESISTIVE ACCELEROMETER WITH UNIFORM AXIAL SENSITIVITIES STUDY OF A THREE-AXIS PIEZORESISTIVE ACCELEROMETER WITH UNIFORM AXIAL SENSITIVITIES Abdelkader Benchou, PhD Canddate Nasreddne Benmoussa, PhD Kherreddne Ghaffour, PhD Unversty of Tlemcen/Unt of Materals

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

On the Multicriteria Integer Network Flow Problem

On the Multicriteria Integer Network Flow Problem BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 5, No 2 Sofa 2005 On the Multcrtera Integer Network Flow Problem Vassl Vasslev, Marana Nkolova, Maryana Vassleva Insttute of

More information

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS

More information

Spin-rotation coupling of the angularly accelerated rigid body

Spin-rotation coupling of the angularly accelerated rigid body Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s

More information

Department of Electrical & Electronic Engineeing Imperial College London. E4.20 Digital IC Design. Median Filter Project Specification

Department of Electrical & Electronic Engineeing Imperial College London. E4.20 Digital IC Design. Median Filter Project Specification Desgn Project Specfcaton Medan Flter Department of Electrcal & Electronc Engneeng Imperal College London E4.20 Dgtal IC Desgn Medan Flter Project Specfcaton A medan flter s used to remove nose from a sampled

More information

Simultaneous Optimization of Berth Allocation, Quay Crane Assignment and Quay Crane Scheduling Problems in Container Terminals

Simultaneous Optimization of Berth Allocation, Quay Crane Assignment and Quay Crane Scheduling Problems in Container Terminals Smultaneous Optmzaton of Berth Allocaton, Quay Crane Assgnment and Quay Crane Schedulng Problems n Contaner Termnals Necat Aras, Yavuz Türkoğulları, Z. Caner Taşkın, Kuban Altınel Abstract In ths work,

More information

The Minimum Universal Cost Flow in an Infeasible Flow Network

The Minimum Universal Cost Flow in an Infeasible Flow Network Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran

More information

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed (2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected

More information

An Interactive Optimisation Tool for Allocation Problems

An Interactive Optimisation Tool for Allocation Problems An Interactve Optmsaton ool for Allocaton Problems Fredr Bonäs, Joam Westerlund and apo Westerlund Process Desgn Laboratory, Faculty of echnology, Åbo Aadem Unversty, uru 20500, Fnland hs paper presents

More information

( ) = ( ) + ( 0) ) ( )

( ) = ( ) + ( 0) ) ( ) EETOMAGNETI OMPATIBIITY HANDBOOK 1 hapter 9: Transent Behavor n the Tme Doman 9.1 Desgn a crcut usng reasonable values for the components that s capable of provdng a tme delay of 100 ms to a dgtal sgnal.

More information

Uncertainty in measurements of power and energy on power networks

Uncertainty in measurements of power and energy on power networks Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:

More information

Design and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm

Design and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm Desgn and Optmzaton of Fuzzy Controller for Inverse Pendulum System Usng Genetc Algorthm H. Mehraban A. Ashoor Unversty of Tehran Unversty of Tehran h.mehraban@ece.ut.ac.r a.ashoor@ece.ut.ac.r Abstract:

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

829. An adaptive method for inertia force identification in cantilever under moving mass

829. An adaptive method for inertia force identification in cantilever under moving mass 89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,

More information

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS) Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998

More information

coordinates. Then, the position vectors are described by

coordinates. Then, the position vectors are described by Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009 College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:

More information

PERFORMANCE OF HEAVY-DUTY PLANETARY GEARS

PERFORMANCE OF HEAVY-DUTY PLANETARY GEARS THE INTERNATIONAL CONFERENCE OF THE CARPATHIAN EURO-REGION SPECIALISTS IN INDUSTRIAL SYSTEMS 6 th edton PERFORMANCE OF HEAVY-DUTY PLANETARY GEARS Attla Csobán, Mhály Kozma 1, 1 Professor PhD., Eng. Budapest

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Annexes. EC.1. Cycle-base move illustration. EC.2. Problem Instances

Annexes. EC.1. Cycle-base move illustration. EC.2. Problem Instances ec Annexes Ths Annex frst llustrates a cycle-based move n the dynamc-block generaton tabu search. It then dsplays the characterstcs of the nstance sets, followed by detaled results of the parametercalbraton

More information

SWITCHING PROCESS IN LIMITED SLIP DIFFERENTIAL

SWITCHING PROCESS IN LIMITED SLIP DIFFERENTIAL 6th Internatonal DAAAM Baltc Conference INDUSTRIAL ENGINEERING 24-26 Aprl 2008, Tallnn, Estona SWITCHING PROCESS IN LIMITED SLIP DIFFERENTIAL Resev, J.; Roosmölder, L.; Stas, M. Abstract: The energy flow

More information

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc

More information

= z 20 z n. (k 20) + 4 z k = 4

= z 20 z n. (k 20) + 4 z k = 4 Problem Set #7 solutons 7.2.. (a Fnd the coeffcent of z k n (z + z 5 + z 6 + z 7 + 5, k 20. We use the known seres expanson ( n+l ( z l l z n below: (z + z 5 + z 6 + z 7 + 5 (z 5 ( + z + z 2 + z + 5 5

More information

Constitutive Modelling of Superplastic AA-5083

Constitutive Modelling of Superplastic AA-5083 TECHNISCHE MECHANIK, 3, -5, (01, 1-6 submtted: September 19, 011 Consttutve Modellng of Superplastc AA-5083 G. Gulano In ths study a fast procedure for determnng the constants of superplastc 5083 Al alloy

More information

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k) ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of

More information

Problem Set 9 Solutions

Problem Set 9 Solutions Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

Predictive Energy Management of a Power-Split Hybrid Electric Vehicle

Predictive Energy Management of a Power-Split Hybrid Electric Vehicle 29 Amercan Control Conference Hyatt Regency Rverfront, St. Lous, MO, USA June 1-12, 29 FrA1.6 Predctve Energy Management of a Power-Splt Hybrd Electrc Vehcle H. Al Borhan Ardalan Vahd Anthony M. Phllps

More information

System in Weibull Distribution

System in Weibull Distribution Internatonal Matheatcal Foru 4 9 no. 9 94-95 Relablty Equvalence Factors of a Seres-Parallel Syste n Webull Dstrbuton M. A. El-Dacese Matheatcs Departent Faculty of Scence Tanta Unversty Tanta Egypt eldacese@yahoo.co

More information

4DVAR, according to the name, is a four-dimensional variational method.

4DVAR, according to the name, is a four-dimensional variational method. 4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The

More information

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve

More information

Curve Fitting with the Least Square Method

Curve Fitting with the Least Square Method WIKI Document Number 5 Interpolaton wth Least Squares Curve Fttng wth the Least Square Method Mattheu Bultelle Department of Bo-Engneerng Imperal College, London Context We wsh to model the postve feedback

More information

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion Assessment of Ste Amplfcaton Effect from Input Energy Spectra of Strong Ground Moton M.S. Gong & L.L Xe Key Laboratory of Earthquake Engneerng and Engneerng Vbraton,Insttute of Engneerng Mechancs, CEA,

More information

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations Quantum Physcs 量 理 Robert Esberg Second edton CH 09 Multelectron atoms ground states and x-ray exctatons 9-01 By gong through the procedure ndcated n the text, develop the tme-ndependent Schroednger equaton

More information

Over-Temperature protection for IGBT modules

Over-Temperature protection for IGBT modules Over-Temperature protecton for IGBT modules Ke Wang 1, Yongjun Lao 2, Gaosheng Song 1, Xanku Ma 1 1 Mtsubsh Electrc & Electroncs (Shangha) Co., Ltd., Chna Room2202, Tower 3, Kerry Plaza, No.1-1 Zhongxns

More information

Design and Analysis of Landing Gear Mechanic Structure for the Mine Rescue Carrier Robot

Design and Analysis of Landing Gear Mechanic Structure for the Mine Rescue Carrier Robot Sensors & Transducers 214 by IFSA Publshng, S. L. http://www.sensorsportal.com Desgn and Analyss of Landng Gear Mechanc Structure for the Mne Rescue Carrer Robot We Juan, Wu Ja-Long X an Unversty of Scence

More information

Physics 181. Particle Systems

Physics 181. Particle Systems Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system

More information

VQ widely used in coding speech, image, and video

VQ widely used in coding speech, image, and video at Scalar quantzers are specal cases of vector quantzers (VQ): they are constraned to look at one sample at a tme (memoryless) VQ does not have such constrant better RD perfomance expected Source codng

More information

Integrals and Invariants of Euler-Lagrange Equations

Integrals and Invariants of Euler-Lagrange Equations Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,

More information

NON LINEAR ANALYSIS OF STRUCTURES ACCORDING TO NEW EUROPEAN DESIGN CODE

NON LINEAR ANALYSIS OF STRUCTURES ACCORDING TO NEW EUROPEAN DESIGN CODE October 1-17, 008, Bejng, Chna NON LINEAR ANALYSIS OF SRUCURES ACCORDING O NEW EUROPEAN DESIGN CODE D. Mestrovc 1, D. Czmar and M. Pende 3 1 Professor, Dept. of Structural Engneerng, Faculty of Cvl Engneerng,

More information

10.34 Numerical Methods Applied to Chemical Engineering Fall Homework #3: Systems of Nonlinear Equations and Optimization

10.34 Numerical Methods Applied to Chemical Engineering Fall Homework #3: Systems of Nonlinear Equations and Optimization 10.34 Numercal Methods Appled to Chemcal Engneerng Fall 2015 Homework #3: Systems of Nonlnear Equatons and Optmzaton Problem 1 (30 ponts). A (homogeneous) azeotrope s a composton of a multcomponent mxture

More information

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.

More information

Supplementary Notes for Chapter 9 Mixture Thermodynamics

Supplementary Notes for Chapter 9 Mixture Thermodynamics Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects

More information

Comparative Analysis of SPSO and PSO to Optimal Power Flow Solutions

Comparative Analysis of SPSO and PSO to Optimal Power Flow Solutions Internatonal Journal for Research n Appled Scence & Engneerng Technology (IJRASET) Volume 6 Issue I, January 018- Avalable at www.jraset.com Comparatve Analyss of SPSO and PSO to Optmal Power Flow Solutons

More information

Some modelling aspects for the Matlab implementation of MMA

Some modelling aspects for the Matlab implementation of MMA Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton

More information

Week 9 Chapter 10 Section 1-5

Week 9 Chapter 10 Section 1-5 Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,

More information

Section 8.3 Polar Form of Complex Numbers

Section 8.3 Polar Form of Complex Numbers 80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

Pivot-Wheel Drive Crab with a Twist! Clem McKown Team November-2009 (eq 1 edited 29-March-2010)

Pivot-Wheel Drive Crab with a Twist! Clem McKown Team November-2009 (eq 1 edited 29-March-2010) Pvot-Wheel Drve Crab wth a Twst! Clem McKown Team 1640 13-November-2009 (eq 1 edted 29-March-2010) 4-Wheel Independent Pvot-Wheel Drve descrbes a 4wd drve-tran n whch each of the (4) wheels are ndependently

More information

Lecture 14: Forces and Stresses

Lecture 14: Forces and Stresses The Nuts and Bolts of Frst-Prncples Smulaton Lecture 14: Forces and Stresses Durham, 6th-13th December 2001 CASTEP Developers Group wth support from the ESF ψ k Network Overvew of Lecture Why bother? Theoretcal

More information

An Admission Control Algorithm in Cloud Computing Systems

An Admission Control Algorithm in Cloud Computing Systems An Admsson Control Algorthm n Cloud Computng Systems Authors: Frank Yeong-Sung Ln Department of Informaton Management Natonal Tawan Unversty Tape, Tawan, R.O.C. ysln@m.ntu.edu.tw Yngje Lan Management Scence

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS

A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS HCMC Unversty of Pedagogy Thong Nguyen Huu et al. A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS Thong Nguyen Huu and Hao Tran Van Department of mathematcs-nformaton,

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

ELE B7 Power Systems Engineering. Power Flow- Introduction

ELE B7 Power Systems Engineering. Power Flow- Introduction ELE B7 Power Systems Engneerng Power Flow- Introducton Introducton to Load Flow Analyss The power flow s the backbone of the power system operaton, analyss and desgn. It s necessary for plannng, operaton,

More information

A Robust Method for Calculating the Correlation Coefficient

A Robust Method for Calculating the Correlation Coefficient A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal

More information

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,

More information

Physics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.

Physics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding. Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays

More information

En Route Traffic Optimization to Reduce Environmental Impact

En Route Traffic Optimization to Reduce Environmental Impact En Route Traffc Optmzaton to Reduce Envronmental Impact John-Paul Clarke Assocate Professor of Aerospace Engneerng Drector of the Ar Transportaton Laboratory Georga Insttute of Technology Outlne 1. Introducton

More information

AERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY

AERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY The Bot-Savart Law The velocty nduced by the sngular vortex lne wth the crculaton can be determned by means of the Bot- Savart formula

More information

Adaptive sliding mode reliable excitation control design for power systems

Adaptive sliding mode reliable excitation control design for power systems Acta Technca 6, No. 3B/17, 593 6 c 17 Insttute of Thermomechancs CAS, v.v.. Adaptve sldng mode relable exctaton control desgn for power systems Xuetng Lu 1, 3, Yanchao Yan Abstract. In ths paper, the problem

More information

Dynamic Programming. Preview. Dynamic Programming. Dynamic Programming. Dynamic Programming (Example: Fibonacci Sequence)

Dynamic Programming. Preview. Dynamic Programming. Dynamic Programming. Dynamic Programming (Example: Fibonacci Sequence) /24/27 Prevew Fbonacc Sequence Longest Common Subsequence Dynamc programmng s a method for solvng complex problems by breakng them down nto smpler sub-problems. It s applcable to problems exhbtng the propertes

More information

FUZZY GOAL PROGRAMMING VS ORDINARY FUZZY PROGRAMMING APPROACH FOR MULTI OBJECTIVE PROGRAMMING PROBLEM

FUZZY GOAL PROGRAMMING VS ORDINARY FUZZY PROGRAMMING APPROACH FOR MULTI OBJECTIVE PROGRAMMING PROBLEM Internatonal Conference on Ceramcs, Bkaner, Inda Internatonal Journal of Modern Physcs: Conference Seres Vol. 22 (2013) 757 761 World Scentfc Publshng Company DOI: 10.1142/S2010194513010982 FUZZY GOAL

More information

The Finite Element Method

The Finite Element Method The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of

More information

6.3.7 Example with Runga Kutta 4 th order method

6.3.7 Example with Runga Kutta 4 th order method 6.3.7 Example wth Runga Kutta 4 th order method Agan, as an example, 3 machne, 9 bus system shown n Fg. 6.4 s agan consdered. Intally, the dampng of the generators are neglected (.e. d = 0 for = 1, 2,

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

More information

Scroll Generation with Inductorless Chua s Circuit and Wien Bridge Oscillator

Scroll Generation with Inductorless Chua s Circuit and Wien Bridge Oscillator Latest Trends on Crcuts, Systems and Sgnals Scroll Generaton wth Inductorless Chua s Crcut and Wen Brdge Oscllator Watcharn Jantanate, Peter A. Chayasena, and Sarawut Sutorn * Abstract An nductorless Chua

More information

A Simple Strategy for Optimal Operation of Heat Exchanger Networks

A Simple Strategy for Optimal Operation of Heat Exchanger Networks A Smple Strategy for Optmal Operaton of Heat Exchanger Networks V. LERSBAMRUNGSUK 1, S. SKOGESTAD 2 AND T. SRINOPHAKUN 1 1 Department of Chemcal Engneerng, Kasetsart Unversty, Bangkok, Thaland 2 Department

More information

A SEPARABLE APPROXIMATION DYNAMIC PROGRAMMING ALGORITHM FOR ECONOMIC DISPATCH WITH TRANSMISSION LOSSES. Pierre HANSEN, Nenad MLADENOVI]

A SEPARABLE APPROXIMATION DYNAMIC PROGRAMMING ALGORITHM FOR ECONOMIC DISPATCH WITH TRANSMISSION LOSSES. Pierre HANSEN, Nenad MLADENOVI] Yugoslav Journal of Operatons Research (00) umber 57-66 A SEPARABLE APPROXIMATIO DYAMIC PROGRAMMIG ALGORITHM FOR ECOOMIC DISPATCH WITH TRASMISSIO LOSSES Perre HASE enad MLADEOVI] GERAD and Ecole des Hautes

More information