Motion Sensorless Control of BLDC PM Motor with Offline FEM Info Assisted State Observer

Size: px
Start display at page:

Download "Motion Sensorless Control of BLDC PM Motor with Offline FEM Info Assisted State Observer"

Transcription

1 2, 2th Internatonal Conference on Optmzaton of Electrcal and Electronc Equpment, OPTIM 2 Moton Sensorless Control of BLDC PM Motor wth Offlne FEM Info Asssted State Observer Aln Ştrban, Ion Boldea, Fellow, IEEE, Gheorghe-Danel Andreescu, Senor Member, IEEE, Dorn Iles, Member, IEEE, and Frede Blaabjerg, Fellow, IEEE Dept. of Electrcal Machnes and Drves, Dept. of Automaton and Appled Informatcs, Unversty Poltehnca of Tmsoara, Tmsoara, Romana, e-mal: strban.aln@yahoo.com, boldea@lselnux.upt.ro, danel.andreescu@aut.upt.ro ebm-papst St. Georgen GmbH&Co. KG., 782 St. Georgen, Germany, e-mal: les@eee.org Insttute of Energy Technology, Aalborg Unversty, Aalborg, Denmark, e-mal: fbl@et.aau.dk Abstract Ths paper descrbes a new offlne FEM asssted poston and speed observer, for brushless dc (BLDC) PM motor drve sensorless control, based on the lne-to-lne PM flux lnkage estmaton. The zero-crossng of the lne-to-lne PM flux lnkage occurs rght n the mddle of two commutaton ponts, whch s used as a bass for the poston and speed observer. For performance applcatons, the poston between commutaton ponts s obtaned by comparng the estmated lne-to-lne PM flux lnkage wth the FEM calculated lne-to-lne PM flux lnkage. Even f the proposed observer reles on the fundamental model of the machne, a safe startng strategy under heavy load torque, called I-f control, s used, wth seamless transton to the proposed sensorless control. The I-f startng method allows lowspeed sensorless control, wthout knowng the ntal poston, and wthout machne parameters dentfcaton. Dgtal smulatons and expermental results are shown, demonstratng the relablty of the FEM asssted poston and speed observer for BLDC PM motor sensorless control operaton. Index Terms Brushless dc (BLDC) motor, sensorless control, startng method, varable speed range. F I. INTRODUCTION OR brushless dc permanent magnet (BLDCPM) motors, many sensorless control strateges for obtanng the rotor poston and speed have been proposed n the lterature [- 6]. One category, the back-emf sensng technque () [3-8], s a scheme estmatng the rotor poston ndrectly, by usng the zero-crossng pont detecton from the termnal voltage of the unenergzed phase wndng. It s the most popular sensorless control method, and has been mproved for wde range speed control [6]. Other sensorless methods are: back- EMF ntegraton technques (2) [9]; flux lnkage-based technque (3) [9]; and freewheelng dode conducton (4) []. All mentoned methods have advantages (easy mplementaton and low computatonal burden) as well as dsadvantages (need of an external hardware crcutry [, 2] and [4], error accumulaton problem at low speeds [2] and [3], poston error n transent state ([] and [4]). References [2] and [3] ntroduce a sensorless poston detecton technque based on a speed-ndependent poston functon for BLDC motors, whch proved good relablty n a wde speed range. Reference [4] proposes the calculaton of the commutaton nstants based on the slope varatons of the common dc current I MAX, whch s obtaned by takng the absolute values of two of the three phase currents. Ths last technque ntroduces unavodable delay n the swtchng status at hgher frequences. Ths paper descrbes a new offlne FEM asssted poston and speed observer for BLDC PM motor drve sensorless control, based on the lne-to-lne PM flux lnkage. Usng measured phase currents and lne-to-lne voltages that can be measured or calculated (by multplcaton of the dc bus voltage by swtchng status, whch are known n the controller), the lne-to-lne PM flux lnkage can be estmated. The zero-crossng of the lne-to-lne PM flux occurs rght n the mddle of two commutatons ponts. At constant or slowly varyng speed, the tme perod from one commutaton pont to zero-crossng, and the tme perod from zero-crossng to the next commutaton pont are equal to each other. Ths s used as a bass for the poston and speed observer. Snce the shape of lne-to-lne PM flux lnkage s dentcal at all speeds, t provdes rather precse commutaton ponts. For BLDC PM motors where the objectve s to acheve quas-square current waveforms, only the poston of commutaton ponts s necessary. But, for specal purpose control strateges (e.g. advanced angle control for the feld weakenng operaton), the poston between commutaton ponts s obtaned by comparng the estmated lne-to-lne PM flux lnkage wth the FEM calculated one (or the one measured under no-load operaton). As the proposed observer reles on the fundamental model of the machne, a safe startng strategy under heavy load torque s stll needed. For the startng method from standstll to a certan low speed, varous methods were proposed: nductance varaton method [] or the algn and go method wth the peak currents lmted durng the startng perod [6]. In ths paper, a safe, robust sensorless start-up method, called I-f control [7], wthout ntal poston detecton, s used. Seamless transton from I-f control to the proposed sensorless control, based on the lne-to-lne PM flux lnkage s acheved. Dgtal smulatons and expermental results demonstrate the relablty of the offlne FEM asssted poston and speed observer for BLDC PM motor sensorless control operaton. The paper s organzed as follows: Secton II: BLDC PM Machne model, Secton III: I-F control startng method, Secton IV: FEM asssted poston and speed observer, Secton V: Transton strateges, Secton VI: Smulaton results, Secton VII: Expermental results, Secton VIII: Concluson //$26. '2 IEEE 32 Authorzed lcensed use lmted to: Gheorghe Andreescu. Downloaded on August 3,2 at 2:6:7 UTC from IEEE Xplore. Restrctons apply.

2 II. BLDC PM MACHINE MODEL Fg. a shows the cross-secton of the nteror permanent magnet brushless dc motor wth 8 poles and (6+6) nonunform slots wth concentrated wndngs, chosen as a case study [8]. The BLDC PM specfcaton data and machne parameters are summarzed n the Appendx. (It has been shown [9] that, wth segmental skewng the full torque pulsatons wth trapezodal shape currents s less than 3%, despte of the IPM character of the rotor). B C B C A A A A Fg.. BLDC PM motor cross-secton. Generally, BLDC PM motors are drven by a three-phase nverter wth what s called sx-step commutaton. Each phase s conductng 2 electrcal degrees. Therefore, for ths knd of machne, only two of the three phases are conductng at any tme, leavng the thrd (floatng) phase, open. Fg. 2 shows an equvalent crcut of the three-phase BLDC PM motor drve. /2V DC /2V DC T T 3 D D 3 V a a V b V c T 4 D 4 b T 6 D 6 T D c T 2 D 2 a V ab b V bc V ca c C B C B R s V bn L bb Fg. 2. Schematc nverter system dagram. R s L E a aaeb n Each phase n an ac motor can be descrbed by a frst order dfferental equaton. The general equatons of the phase voltages can be represented n the matrx form, as: ( θ,,, ) ( θ,,, ) ( θ,,, ) Va a Ψ A r a b c Vn d Vb = Rs b + ΨB r a b c Vn dt + Vc c ΨC r a b c V n where V x s the phase voltage, R s s the phase resstance, x s the phase current, θ r s the rotor poston, Ψ x s the total flux lnkage of the wndng, and V n s the voltage between motor and nverter neutral ponts. The flux lnkage n a phase ncludes both self and mutual flux lnkages, and may be defned n the followng form: Ψ A Laa( θr, Lab( θr, Lac( θr, a λpma( θr) Ψ B Lba( θr, Lbb( θr, Lbc( θr, = b + λpmb( θr) Ψ L ( θ, ) L ( θ, ) L ( θ, ) λ ( θ ) C ca r a cb r b cc r c c PMc r V an Ec V cn L cc R s () (2) where, the term, L xx (θ r, x ) x, represents the self flux lnkage of phase x; the term, L xy (θ r, y ) y stands for the mutual flux lnkage between phases x and y; and the term, λ PMx (θ r ) s the phase flux lnkage n phase x, due to the permanent magnets. The dependence of machne nductances on currents has been thoroughly treated n [9]; here, however, for control desgn, approxmated constant values are consdered (see Appendx). The phase flux lnkages λ PMx (θ r ), are modeled through Fourer seres, based on FEM-analyss of the BLDC IPM motor: λ θ = A cos k p θ ( ) ( ) PMa r k r k =,3, 2 PMb ( r ) = Ak cos k p r 3 k =,3, 2 PMc ( r ) = Ak cos k p r + 3 k =,3, λ θ θ λ θ θ π π where p s the number of pole pars, and the terms A k, are the frst three coeffcents of Fourer seres (n Appendx, Table II). The back-emf equaton s gven by dλpmx dλpmx dθ e r x = = = kex ωm (4) dt dθr dt where k ex stands for back-emf constants, ω m s the mechancal rotor speed. Substtutng (2) and (4) nto (), the phase voltage equatons become: Va a Laa Lab Lac a ea Vn d Vb = Rs b + Lba Lbb Lbc b eb Vn dt + + () Vc b Lca Lcb Lcc b e c V n Snce the neutral pont of the motor s not avalable, the lne-to-lne voltage equatons are used as follows: Vab a Laa Lab Lac a eab d Vbc = Rs K b + K Lba Lbb Lbc b + ebc dt Vca c Lca Lcb Lcc c e ca where K s the matrx transformaton from phase quanttes to lne-to-lne quanttes, and has the followng form, K = From (6) and knowng the fact that the lne-to-lne back EMF s the lne-to-lne PM flux lnkage dervatve, the lne-tolne PM flux lnkages can be evaluated n the followng form: λ t PMab Vab a λpmbc = Vbc Rs K b dτ λ PMca Vca c Laa Lab Lac a + K Lba Lbb Lbc b Lca Lcb L cc c Drect measurement of the phase currents and lne-to-lne voltages, can allow the estmaton of the lne-to-lne PM flux (3) (6) (7) (8) 322 Authorzed lcensed use lmted to: Gheorghe Andreescu. Downloaded on August 3,2 at 2:6:7 UTC from IEEE Xplore. Restrctons apply.

3 lnkages. The lne-to-lne voltages can also be calculated, by multplcaton of the dc bus voltage and swtchng status. From these lne-to-lne PM flux lnkages, the poston of the rotor can be estmated as explaned n the next sectons. III. I-F CONTROL STARTING METHOD I-f control method conssts n rampng the stator current frequency whle the current reference s kept constant by the current controller [7]. In ths paper, durng the experments, n order to have a fnte acceleraton durng start-up, the frequency tme-varaton form was chosen as: fmax ( cos ( ω f t) ), ωf t π f = 2 (9) fmax, ωf > t π ωf tf = π 2 π ff tf = π ff = () 2 t f where t f s the tme perod of the reference frequency ramp (t f =.2 [sec], n our case). The frequency s used to calculate a current reference angle θ er, whch s the feedforward poston-angle correspondng to the mposed synchronous reference frame: IV. θ () er = 2 π f dt FEM ASSISTED POSITION AND SPEED OBSERVER A. Commutaton Strategy and Poston Observer Fg. 3 shows the back-emfs, phase currents, lne-to-lne PM flux lnkages of the BLDC PM motor n [8], and the detecton of the zero crossngs of the PM lne-to-lne flux lnkages. Back-EMFs [V] Phase Currents [A] Lne-to-lne PM flux lnkages [Wb] Zero-crossng detecton of the lne PM flux lnakges a b c x -3 Mode Mode 2 Mode 3 Mode 4 Mode Mode 6 λ PMbc λ PMca e a e a e b e b e c e c tme [s] Fg. 3. Bass of the sensorless operaton of the BLDC PM: back-emfs; phase currents; lne-to-lne PM flux lnkages; zero-crossng detecton of the lne-to-lne PM flux lnkages. As can be seen n Fg. 3, the zero-crossngs of the lne-tolne PM flux lnkages occur rght n the mddle of two commutatons ponts. At constant or slowly varyng speed, the tme perod from one commutaton pont to zero-crossng, and the tme perod from zero-crossng to the next commutaton pont are equal to each other. Ths s used as a bass for the mplementaton of the poston and speed sensorless observer. In each mode (sector), only two of the three phases are conductng at any tme, leavng the thrd floatng phase open. The phase voltage n the floatng phase s equal wth the Back-EMF. Only one lne-to-lne PM flux lnkage (from Eq. (8)) s used n each mode: ths s the PM flux lnkage between the two actve phases. Based on ths dea, n each mode, the proper estmated lne-to-lne PM flux lnkage can be used to detect the zero-crossng, and then a delay of 3 s appled to swtch-on the current n the correspondng phases (next commutaton pont). After commutaton s produced, a new couple of phases wll be conductng, and the correspondng lne-to-lne PM flux lnkage wll gve the nformaton for the next commutaton. Table I shows the proper lne-to-lne PM flux lnkage at each mode. TABLE I. LINE-TO-LINE PM FLUX LINKAGES AT EACH MODE Mode and 4 Mode 2 and Mode 3 and 6 Postve speed λ PMca λ PMbc The lne-to-lne PM flux lnkage estmator (Fg. 4) based on the voltage model n stator reference frame, employs an equvalent ntegrator n close-loop, wth a speed-adaptve proportonal controller to compensate the dc offset and phasedelay. The transfer functon of ths equvalent ntegrator s: T H( s) = ; T = (2) + st ω where ω s the corner frequency. For BLDC PM motor general drves, where the objectve s just to acheve quas-square current waveforms, knowng the poston of commutaton ponts suffces. But, for specal purpose drves wth advance angle control for the feld weakenng operaton, the poston between commutaton ponts can be obtaned by comparng the estmated lne-tolne PM flux lnkage wth a look-up table, whch contans the poston versus characterzed lne-to-lne PM flux lnkage. Ths PM flux lnkage can be calculated by usng FEMmethod [8], f the motor was desgned usng ths method, or measured by runnng the machne as a generator wth open phases at constant speed. The lne voltage sgnals are recorded and by ntegratng these voltages, the lne-to-lne PM flux lnkages can be calculated. The poston between commutaton ponts can be estmated because the correspondng lne-to-lne PM flux lnkage at each mode s ether ncreasng or decreasng. 323 Authorzed lcensed use lmted to: Gheorghe Andreescu. Downloaded on August 3,2 at 2:6:7 UTC from IEEE Xplore. Restrctons apply.

4 I abc R s V ab, V bc,v ca Κ I abc [L] X Κ ω ω r λ PM lneto-lne Fg. 4. Lne-to-lne PM flux estmator (K s from Eq. (7)). λpm(θer) θer Lookup Table B. Speed Estmaton Wth no sensors avalable, the speed must be calculated (for speed feedback) usng the nformaton from the commutaton controller. The mechancal model of the drve s approxmated as: θer θer d ω r p J ω r p J Te dt = + TL TL θˆer (3) where θ er s the rotor poston, ω r s the rotor electrcal speed and the load torque T L was assumed constant. The electromagnetc torque T e s regarded as known nput and s calculated by the followng equaton: ea a + eb b + ec c Te = = kea a + keb b + kec c (4) ωm From the mechancal model an observer can be constructed: ˆ θ ˆ er θ er k d ˆ ω ˆ ˆ r p J ω r p J = Te k2 ε dt + + () Tˆ ˆ k L T L 3 The PLL observer error s: ( ˆ er ˆ r ) ( ˆ er ) ( ˆ r ) ( ˆ er ) ( ˆ r ) ε = sn θ θ = sn θ cos θ cos θ sn θ (6) The PLL speed observer block dagram s shown n Fg.. Ths observer estmates the rotor speed and the load torque. However, the load torque estmaton s poor. TˆL Tˆe p/j ˆr ω θˆr J k = S, k2 = SP2, k3 = P (7) p where: S = p + p 2 + p 3, SP 2 = p p 2 + p 2 p 3 + p 3 p, P = p p 2 p 3. Gan values are lmted only by nose consderaton. If the nose level s small, t s possble to select all three poles real, negatve. However, for better dynamcs, two of them can be selected complex conjugates wth not too large magnary parts. V. TRANSITION STRATEGIES A. Transton from I-f Control to FEM Asssted Poston and Speed Observer I-f control method s used only for start-up and low speeds. When the reference frequency exceeds a certan level (f mn ), the system automatcally ntates the transton to the sensorless control operaton of the BLDC, based on FEM asssted poston and speed observer (Fg. 6. In order to acheve a smooth transton, the angle used for control has to er go smoothly from θ to ˆer θ : where: ( k) k ˆ θ = θ + θ (8) < ; for f f mn f fmn k = ; for fmn f f fmax fmn ; for f > fmax er er max (9) When the reference frequency reaches f max, the control swtches to FEM asssted poston and speed observer moton sensorless control (k=). Because the I-f feedforward poston-angle and the estmated poston are kept between and 2 π, based on the fact that Δθ sn ( Δθ) for small Δθ, (3) s mplemented n the followng form: θer ; for f < f mn ˆ = er + k sn ( er er ); for fmn f fmax ˆ θer ; for f > f max θ θ θ θ (2) k 3 k 2 k ε sn Δθ θˆer Fg.. Equvalent PLL state-observer structure The observer gans are determned by pole placement method. Let p, p 2, and p 3 be the observer poles, allocated such that the observer s fast and stable. In ths case, the observer gans are: Fg. 6. Control strategy: ) I-f control; 2) FEM asssted poston and speed observer moton sensorless control; 3) I-f control. 324 Authorzed lcensed use lmted to: Gheorghe Andreescu. Downloaded on August 3,2 at 2:6:7 UTC from IEEE Xplore. Restrctons apply.

5 B. Transton from FEM Asssted Poston and Speed Observer to I-f Control The transton from FEM asssted poston and speed observer moton sensorless control, to I-f control, s done automatcally by the system when the reference speed s below a certan level (Fg. 6. When the above condton s met, the control system generates an angle ( θ er - feedforward poston angle): f abs (2) er= dt= 2 f dt θ ω π 2 π f mn back θˆer u x θ er Fg. 7. Transton strategy from FEM asssted poston and speed observer moton sensorless control to I-f method. The transton to I-f control s performed by settng the ntal condton of the feedforward poston angle θ er, equal to the last value of the estmated rotor poston θ er (Fg. 7). VI. SIMULATION RESULTS The FEM asssted poston and speed observer for a brushless dc PM motor drve sensorless control, have been frst verfed through dgtal smulaton. A BLDC PM wth the characterstcs summarzed n the Appendx, was mplemented n Matlab/Smulnk, usng the general equatons of a brushless dc PM motor. In the mplementaton, the phase commutaton phenomenon was consdered [2]. In order to obtan accurate and effcent smulaton [2], the back-emf and the nductances (IPM rotor) were modeled through Fourer seres, and through look-up tables, based on FEM-method calculaton of the BLDC IPM motor, respectvely. Fg. 8 llustrates the performance of the proposed FEM asssted poston and speed observer moton sensorless control at RPM. The phase currents are controlled by a PI controller. The commutaton ponts from Fg. 8, represent n fact, the zero-crossng detecton of the lne-to-lne PM flux lnkages, delayed wth 3. Phase A Commutaton ponts Fg. 8. (contnue. Lne-to-lne PM flux lnkages [Wb] Poston [deg.] Poston error [deg.] Speed [rpm] x -3 estmated from FEM estmated encoder estmated speed actual speed f) Fg. 8. Smulaton results of moton-sensorless control of BLDC PM motor wth the proposed poston and speed FEM asssted observer at RPM, at. Nm load: phase A current; commutaton ponts; lne-to-lne PM flux lnkage; electrcal poston; e) poston error; f) speed. The numercal smulatons demonstrate the relablty of ths observer. Next, expermental results are presented, n order to valdate the FEM asssted poston and speed observer moton sensorless control n real-tme applcatons. VII. EXPERIMENTAL RESULTS The expermental setup conssts of BLDC IPM motor, wth the parameters as n smulatons (Appendx), loaded wth an dentcal nteror permanent magnet synchronous generator wth constant load. The dspace CLP3 platform s used to control the entre drve system (Fg. 9). Encoder BLDC PM motor 3-phase nverter & drver crcut IPMSG Current measurement devce Fg. 9. Expermental set-up. e) 32 Authorzed lcensed use lmted to: Gheorghe Andreescu. Downloaded on August 3,2 at 2:6:7 UTC from IEEE Xplore. Restrctons apply.

6 V dc ω r ˆr ω k pw, k w s s k p, k V dc Duty cycle drecton COMMUTATION STRATEGY ˆer θ CP S a, S b, S c : swtchng status PWM INVERTER a b Vab Vbc BLDC Motor s drecton - z Duty cycle MAX c ˆr ω Speed Estmaton (Fg. ) CP θˆer Lne-to-lne PM flux estmator (Fg. 4) ω r V ca Fg.. Block dagram of the proposed sensorless drve (CP - commutaton ponts). Fg. shows an overall system block dagram of the proposed sensorless drve. In ths experment, the lne-to-lne voltages are measured, but they also can be calculated by multplcaton of the dc bus voltage by swtchng status, whch are known n the controller. The use of the measured lne-to-lne voltages nstead of the calculated ones was done n order to acheve better stablty and dynamc responses durng the sensorless operaton. The performance of the proposed moton sensorless control system s nvestgated n four cases: constant speed, speed varatons, start-up wth I-f method and transton to the FEM asssted poston and speed observer, and back when the reference speed s below a certan level. A. Constant Speed Operaton The sensorless operaton of the BLDC motor wth the proposed FEM asssted poston and speed observer at constant speed (Fg. rpm, loaded at. Nm), shows a rather satsfactory agreement between encoder poston and estmated poston, certfyng good estmaton results. In both cases there are some dfferences between the estmated and measured speeds. Phase A Lne-to-lne PM flux lnkages [Wb] - - x -3 estmated - from FEM Fg.. (contnue. Poston [deg.] Poston error [deg.] Speed [rpm] estmated encoder estmated speed measured speed 9 e) Fg.. Expermental results of moton-sensorless control of BLDC PM motor wth the proposed poston and speed FEM asssted observer at RPM, at. Nm load: phase A current; lne-to-lne PM flux lnkage; electrcal poston; poston error; e) spee. B. Dynamcs Fg. 2 shows f the FEM asssted poston and speed observer s relable under acceleraton and deceleraton, and that the load s proportonal to speed. Durng these transent tests, the agreement between encoder poston and estmated poston s rather satsfactory. The error between these two postons may be mnmzed by a more carefully tunng of the poston and speed observer. From Fg. 2 t s obvous that the speed estmated wth the PLL observer (Fg. ) s n agreement wth the measured speed, n both steady state and transents. 326 Authorzed lcensed use lmted to: Gheorghe Andreescu. Downloaded on August 3,2 at 2:6:7 UTC from IEEE Xplore. Restrctons apply.

7 Phase A Speed [rpm] estmated speed measured speed Lne-to-lne PM flux lnkages [Wb] Poston error [deg.] Speed [rpm] x -3 estmated from FEM estmated speed measured speed reference speed Fg. 2. Expermental results of moton-sensorless control of BLDC PM motor wth the proposed poston and speed FEM asssted observer for acceleraton from to RPM and deceleraton from to RPM: phase A current; lne-to-lne PM flux lnkage; poston error; speed. C. Transton from I-F Control to FEM Asssted Poston and Speed Observer BLDC PM Moton-Sensorless Control In ths paper, for a smooth start-up, the frequency tmevaraton form was chosen as n (9), up to Hz ( rpm). The control strategy automatcally swtches from I-f control method to the FEM asssted poston and speed observer moton sensorless control, when the reference frequency s between 7 and Hz (Fg. 3). Phase A Poston [deg.] estmated encoder Fg. 3. (contnue. Speed error [rpm] Poston [deg.] Fg. 3. Transton from I-f control method to FEM asssted poston and speed observer: phase A current; poston; speed; and speed error I-f estmated w Start transton End transton Fg. 4. I-f poston, estmated rotor poston, and the angle used for control [deg.], n transton to FEM asssted poston and speed observer. It s obvous from Fg. 4 that when the system automatcally swtches between the two control methods, the angle used for control swtches smoothly (Eq. (2)) from I-f feedforward poston to estmated poston D. Transton from FEM Asssted Poston and Speed Observer BLDC PM Moton-Sensorless Control to I-F Control When the reference frequency s below a certan level (2 Hz), the control strategy automatcally swtches from FEM asssted poston and speed observer moton sensorless control, to I-f control (Fg. ) wth currents frequency equal wth the reference frequency. A seamless transton from the estmated poston to I-f feedforward poston s vsble n Fg. 6 by usng the soluton from Fg. 7. Phase A Poston [deg.] estmated encoder Fg.. (contnue. 327 Authorzed lcensed use lmted to: Gheorghe Andreescu. Downloaded on August 3,2 at 2:6:7 UTC from IEEE Xplore. Restrctons apply.

8 Speed [rpm] Poston [deg.] estmated speed measured speed reference speed Fg.. Transton FEM asssted poston and speed observer to I-f control method: phase A current; poston; and speed To I-f I-f estmated Fg. 6. I-f poston and estmated rotor poston [deg.] n transton to I-f control. VIII. CONCLUSION The paper has ntroduced a FEM asssted poston and speed observer for brushless dc PM motor drve sensorless control, based on the lne-to-lne PM flux lnkage. The lne-to-lne PM flux lnkage can be estmated from measured phase currents and measured (or calculate lneto-lne voltages. The zero-crossng of the lne-to-lne PM flux occurs rght n the mddle of two commutatons ponts. Ths s used as a bass for the poston and speed observer. Ths sensorless control method s used together wth the I- f sensorless control, for start-up and for low speed control, wth seamless transtons between them. Dgtal smulatons and expermental results demonstrate that the FEM asssted poston and speed observer for BLDC PM motor sensorless control operaton, provdes rather precse commutaton ponts, even durng speed or load transent states. APPENDIX BLDC PM parameters: 2 p = 8 rotor poles, (6+6) slots, rated torque T emax =.4 [Nm], rated speed n n = [rpm], dc bus voltage V DC = 2 [V], very low coggng torque T cogg <% T emax, phase resstance R s = [mω], nerta: J= [kg m 2 ], constant values for nductances: L s = L aa = L bb = L cc = [uh], M = L ab = L ac = L ba = L bc = L ca = L cb = 7. [uh] [8]. TABLE II. FOURIER SERIES COEFFICIENTS A A 3 A [Wb] [Wb] [Wb] ACKNOWLEDGEMENT Ths work was partally supported by EU-FP7 EE-VERT Project (Grant agreement no. SCP7-GA EE- VERT). REFERENCES [] J.P. Johnson, M. Ehsan, and Y. Güzelgünler, Revew of Sensorless Methods for Brushless DC, IAS Annual Meetng 999, vol., pp.43-, October 999. [2] I. Boldea and S.A. Nasar, Electrc Drves, 2 nd ed., CRC Press, Taylor & Francs, New York, 2. [3] AVR444: Sensorless Control of 3-Phase Brushless DC Motors, 2. [4] J. Shao, D. Nalon, and T. Hopkns, A Novel Drect Back EMF Detecton for Sensorless Brushless DC (BLDC) Motor Drves, Proc. of IEEE APEC, pp.33-37, 22. [] Y.-S. La and Y.-K. Ln, A Unfed Approach to Back-EMF Detecton for Brushless DC Motor Drves wthout Current and Hall Sensors, IECON 26, pp , November 26. [6] Y.-S. La and Y.-K. Ln, Back-EMF Detecton Technque of Brushless DC Motor Drve for Wde Range Control, IECON 26, pp.6-, November 26. [7] L. Mngyao, Z. Zhyao, and L. Keman, A Novel and Easy-Realzng Intal Rotor Poston Detecton Method and Speedup Algorthm for Sensorless BLDC Motor Drves, ICEMS 28, pp , October 28. [8] C.-T. Ln, C.-W. Hung, and C.-W. Lu, Sensorless Control for Four- Swtch Three-Phase Brushless DC Motor Drve, Conference Record of the IEEE on Industry Applcatons, vol. 4, pp , October 26. [9] J. X. Shen and S. Iwasak, Sensorless Control of Ultrahgh-Speed PM Brushless Motor Usng PLL and Thrd Harmonc Back EMF, IEEE Trans. on Industral Electroncs, vol. 3, no. 2, pp , Aprl 26. [] N. Ertugrul and P. Acarnely, A New Algorthm for Sensorless Operaton of Permanent Magnets Motors, IEEE Trans. on Industry Applcatons, vol. 3, no., pp.26-33, January/February 994. [] S. Ogasawara and H. Akag, An Approach to Poston Sensorless Drve for Brushless DC Motors, IEEE Trans. Ind. Applcatons, vol. 27, pp , September/October 99. [2] T.-H. Km and M. Ehsan, Sensorless Control of the BLDC Motors from Near-Zero to Hgh Speeds, IEEE Trans. on Power Electroncs, vol. 9, no. 6, pp.63-64, November 24. [3] T.-S. Km, B.-G. Park, D.-M. Lee, J.-S. Ryu, and D.-S. Hyun, A New Approach to Sensorless Control Method for Brushless DC Motors, Internatonal Journal of Control, Automaton, and Systems, vol. 6, no. 4, pp , August 28. [4] J. Dxon, M. Rodrguez, and R. Huerta, Poston Estmator and Smplfed Current Control Strategy for Brushless-DC Motors, Usng DSP Technology, IECON 2, vol., pp.9-96, November 22. [] K.-W. Lee, D.-K. Km, B.-T. Km, and B.-I. Kwon, A Novel Startng Method of the Surface Permanent-Magnet BLDC Motors Wthout Poston Sensor for Recprocatng Compressor, IEEE Trans. on Industry Applcatons, vol. 44, no., pp.8-92, January/February 28. [6] W.-J. Lee and S.-K. Sul, A New Startng Method of BLDC Motors Wthout Poston Sensor, IEEE Trans. on Industry Applcatons, vol. 42, no. 6, pp.32-38, November/December 26. [7] M. Fătu, R. Teodorescu, I. Boldea, G.-D. Andreescu, and F. Blaabjerg, I-F Startng Method wth Smooth Transton to EMF Based Moton-Sensorless Vector Control of PM Synchronous Motor/Generator, PESC 28, pp , June 28. [8] A. Ştrban, D. Iles-Klumpner, M. Rstcevc, and I. Boldea, (6+6) Slot/ 8 Pole 3 Phase IPM Brushless D.C. Automotve Actuator: Torque, EMF and Inductances Characterzaton by FEM vs. Experments, IEEE- IEMDC 29, May 29. [9] A. Ştrban, L. Tutelea, D. Iles-Klumpner, and I. Boldea, FEM analyss of concentrated cols nonunform slot (6+6/8) IPMSM fed wth trapezodal current, OPTIM 28, pp. 4-2, May 28. [2] W. Hong, W. Lee, and B.-K. Lee, Dynamc Smulaton of Brushless DC Motor Drves Consderng Phase Commutaton for Automotve Applcatons, IEMDC 7, vol. 2, pp , May 27. [2] O.A. Mohammed, S. Lu, and Z. Lu, A Phase Varable Model of Brushless DC Motors Based on Fnte Element Analyss and Its Couplng wth External Crcuts, IEEE Trans. on Magnetcs, vol. 4, no., pp.76-79, May Authorzed lcensed use lmted to: Gheorghe Andreescu. Downloaded on August 3,2 at 2:6:7 UTC from IEEE Xplore. Restrctons apply.

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

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

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

Analysis of the Magnetomotive Force of a Three-Phase Winding with Concentrated Coils and Different Symmetry Features

Analysis of the Magnetomotive Force of a Three-Phase Winding with Concentrated Coils and Different Symmetry Features Analyss of the Magnetomotve Force of a Three-Phase Wndng wth Concentrated Cols and Dfferent Symmetry Features Deter Gerlng Unversty of Federal Defense Munch, Neubberg, 85579, Germany Emal: Deter.Gerlng@unbw.de

More information

On-line Parameter Estimation of Interior Permanent Magnet Synchronous Motor Using an Extended Kalman Filter

On-line Parameter Estimation of Interior Permanent Magnet Synchronous Motor Using an Extended Kalman Filter J Electr Eng echnol Vol. 9, No.?: 742-?, 2014 http://dx.do.org/10.5370/jee.2014.9.?.742 ISSN(Prnt) 1975-0102 ISSN(Onlne) 2093-7423 On-lne Parameter Estmaton of Interor Permanent Magnet Synchronous Motor

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

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

Modeling and Control of SRM in Electromechanic Brake System

Modeling and Control of SRM in Electromechanic Brake System Modelng and Control of SRM n Electromechanc Brake System Wenzhe u, Ph.D. student Al Keyhan, Advsor Mechatroncs aboratory Department of Electrcal and Computer Engneerng The Oho State Unversty Agenda I.

More information

Equivalent Circuit Analysis of Interior Permanent Magnet Synchronous Motor Considering Magnetic saturation

Equivalent Circuit Analysis of Interior Permanent Magnet Synchronous Motor Considering Magnetic saturation Page 0114 World Electrc Vehcle Journal Vol. 3 - ISSN 2032-6653 - 2009 AVERE EVS24 Stavanger, Norway, May 13-16, 2009 Euvalent Crcut Analyss of Interor Permanent Magnet Synchronous Motor Consderng Magnetc

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

( ) = ( ) + ( 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

(b) i(t) for t 0. (c) υ 1 (t) and υ 2 (t) for t 0. Solution: υ 2 (0 ) = I 0 R 1 = = 10 V. υ 1 (0 ) = 0. (Given).

(b) i(t) for t 0. (c) υ 1 (t) and υ 2 (t) for t 0. Solution: υ 2 (0 ) = I 0 R 1 = = 10 V. υ 1 (0 ) = 0. (Given). Problem 5.37 Pror to t =, capactor C 1 n the crcut of Fg. P5.37 was uncharged. For I = 5 ma, R 1 = 2 kω, = 5 kω, C 1 = 3 µf, and C 2 = 6 µf, determne: (a) The equvalent crcut nvolvng the capactors for

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

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

6.3.4 Modified Euler s method of integration

6.3.4 Modified Euler s method of integration 6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from

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

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

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

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient 58:080 Expermental Engneerng 1 OBJECTIVE Lab 2e Thermal System Response and Effectve Heat Transfer Coeffcent Warnng: though the experment has educatonal objectves (to learn about bolng heat transfer, etc.),

More information

Appendix B. The Finite Difference Scheme

Appendix B. The Finite Difference Scheme 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton

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

PHY2049 Exam 2 solutions Fall 2016 Solution:

PHY2049 Exam 2 solutions Fall 2016 Solution: PHY2049 Exam 2 solutons Fall 2016 General strategy: Fnd two resstors, one par at a tme, that are connected ether n SERIES or n PARALLEL; replace these two resstors wth one of an equvalent resstance. Now

More information

CONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING INTRODUCTION

CONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING INTRODUCTION CONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING N. Phanthuna 1,2, F. Cheevasuvt 2 and S. Chtwong 2 1 Department of Electrcal Engneerng, Faculty of Engneerng Rajamangala

More information

PHYSICS - CLUTCH CH 28: INDUCTION AND INDUCTANCE.

PHYSICS - CLUTCH CH 28: INDUCTION AND INDUCTANCE. !! www.clutchprep.com CONCEPT: ELECTROMAGNETIC INDUCTION A col of wre wth a VOLTAGE across each end wll have a current n t - Wre doesn t HAVE to have voltage source, voltage can be INDUCED V Common ways

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

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

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

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

Pop-Click Noise Detection Using Inter-Frame Correlation for Improved Portable Auditory Sensing

Pop-Click Noise Detection Using Inter-Frame Correlation for Improved Portable Auditory Sensing Advanced Scence and Technology Letters, pp.164-168 http://dx.do.org/10.14257/astl.2013 Pop-Clc Nose Detecton Usng Inter-Frame Correlaton for Improved Portable Audtory Sensng Dong Yun Lee, Kwang Myung Jeon,

More information

Digital Signal Processing

Digital Signal Processing Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over

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

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

INRUSH CURRENT. By Mario Célio Contin. R&D Department WEG Equipamentos Elétricos S.A. - Motores.

INRUSH CURRENT. By Mario Célio Contin. R&D Department WEG Equipamentos Elétricos S.A. - Motores. 1 INUSH UENT y Maro élo ontn &D Department WEG Equpamentos Elétrcos S.A. - Motores. 1. INTODUTION Any nductve system, whch can be represented by a seres-connected col wth a resstor (g.1), when energzed,

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

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

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

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

APPLICATION OF A SLIDING MODE OBSERVER FOR SENSORLESS OPERATION OF SWITCHED RELUCTANCE MOTORS. Stefan Brock

APPLICATION OF A SLIDING MODE OBSERVER FOR SENSORLESS OPERATION OF SWITCHED RELUCTANCE MOTORS. Stefan Brock Preprnt of the paper publshed n: Archves of Electrcal Engneerng Vol. 56, no 2 (2007) pp: 163--172 http://www.aee.put.poznan.pl/ DOI: http://dx.do.org/10.6084/m9.fgshare.729075 Abstract APPLICATION OF A

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

Boise State University Department of Electrical and Computer Engineering ECE 212L Circuit Analysis and Design Lab

Boise State University Department of Electrical and Computer Engineering ECE 212L Circuit Analysis and Design Lab Bose State Unersty Department of Electrcal and omputer Engneerng EE 1L rcut Analyss and Desgn Lab Experment #8: The Integratng and Dfferentatng Op-Amp rcuts 1 Objectes The objectes of ths laboratory experment

More information

Adiabatic Sorption of Ammonia-Water System and Depicting in p-t-x Diagram

Adiabatic Sorption of Ammonia-Water System and Depicting in p-t-x Diagram Adabatc Sorpton of Ammona-Water System and Depctng n p-t-x Dagram J. POSPISIL, Z. SKALA Faculty of Mechancal Engneerng Brno Unversty of Technology Techncka 2, Brno 61669 CZECH REPUBLIC Abstract: - Absorpton

More information

Numerical Transient Heat Conduction Experiment

Numerical Transient Heat Conduction Experiment Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use

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

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

Electrical Circuits 2.1 INTRODUCTION CHAPTER

Electrical Circuits 2.1 INTRODUCTION CHAPTER CHAPTE Electrcal Crcuts. INTODUCTION In ths chapter, we brefly revew the three types of basc passve electrcal elements: resstor, nductor and capactor. esstance Elements: Ohm s Law: The voltage drop across

More information

CS 468 Lecture 16: Isometry Invariance and Spectral Techniques

CS 468 Lecture 16: Isometry Invariance and Spectral Techniques CS 468 Lecture 16: Isometry Invarance and Spectral Technques Justn Solomon Scrbe: Evan Gawlk Introducton. In geometry processng, t s often desrable to characterze the shape of an object n a manner that

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

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

Electric Machinery and Apparatus 1 AE1B14SP1. Miroslav Chomát room B3-248

Electric Machinery and Apparatus 1 AE1B14SP1. Miroslav Chomát room B3-248 Electrc Machnery and Apparatus 1 AE1B14SP1 Mroslav Chomát chomat@fel.cvut.cz room B3-48 Inducton Machne (IM) Applcatons Constructon Prncple Equatons Equvalent crcut Torque-speed characterstc Crcle dagram

More information

Outline. Communication. Bellman Ford Algorithm. Bellman Ford Example. Bellman Ford Shortest Path [1]

Outline. Communication. Bellman Ford Algorithm. Bellman Ford Example. Bellman Ford Shortest Path [1] DYNAMIC SHORTEST PATH SEARCH AND SYNCHRONIZED TASK SWITCHING Jay Wagenpfel, Adran Trachte 2 Outlne Shortest Communcaton Path Searchng Bellmann Ford algorthm Algorthm for dynamc case Modfcatons to our algorthm

More information

Gravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11)

Gravitational Acceleration: A case of constant acceleration (approx. 2 hr.) (6/7/11) Gravtatonal Acceleraton: A case of constant acceleraton (approx. hr.) (6/7/11) Introducton The gravtatonal force s one of the fundamental forces of nature. Under the nfluence of ths force all objects havng

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016 ME140 - Lnear rcuts - Wnter 16 Fnal, March 16, 2016 Instructons () The exam s open book. You may use your class notes and textbook. You may use a hand calculator wth no communcaton capabltes. () You have

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

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 6. Operational Amplifier. inputs can be defined as the average of the sum of the two signals.

Chapter 6. Operational Amplifier.  inputs can be defined as the average of the sum of the two signals. 6 Operatonal mpler Chapter 6 Operatonal mpler CC Symbol: nput nput Output EE () Non-nvertng termnal, () nvertng termnal nput mpedance : Few mega (ery hgh), Output mpedance : Less than (ery low) Derental

More information

Chapter 4 The Wave Equation

Chapter 4 The Wave Equation Chapter 4 The Wave Equaton Another classcal example of a hyperbolc PDE s a wave equaton. The wave equaton s a second-order lnear hyperbolc PDE that descrbes the propagaton of a varety of waves, such as

More information

PHYSICS - CLUTCH 1E CH 28: INDUCTION AND INDUCTANCE.

PHYSICS - CLUTCH 1E CH 28: INDUCTION AND INDUCTANCE. !! www.clutchprep.com CONCEPT: ELECTROMAGNETIC INDUCTION A col of wre wth a VOLTAGE across each end wll have a current n t - Wre doesn t HAVE to have voltage source, voltage can be INDUCED V Common ways

More information

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification E395 - Pattern Recognton Solutons to Introducton to Pattern Recognton, Chapter : Bayesan pattern classfcaton Preface Ths document s a soluton manual for selected exercses from Introducton to Pattern Recognton

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

THE CURRENT BALANCE Physics 258/259

THE CURRENT BALANCE Physics 258/259 DSH 1988, 005 THE CURRENT BALANCE Physcs 58/59 The tme average force between two parallel conductors carryng an alternatng current s measured by balancng ths force aganst the gravtatonal force on a set

More information

Analytical Gradient Evaluation of Cost Functions in. General Field Solvers: A Novel Approach for. Optimization of Microwave Structures

Analytical Gradient Evaluation of Cost Functions in. General Field Solvers: A Novel Approach for. Optimization of Microwave Structures IMS 2 Workshop Analytcal Gradent Evaluaton of Cost Functons n General Feld Solvers: A Novel Approach for Optmzaton of Mcrowave Structures P. Harscher, S. Amar* and R. Vahldeck and J. Bornemann* Swss Federal

More information

Temperature. Chapter Heat Engine

Temperature. Chapter Heat Engine Chapter 3 Temperature In prevous chapters of these notes we ntroduced the Prncple of Maxmum ntropy as a technque for estmatng probablty dstrbutons consstent wth constrants. In Chapter 9 we dscussed the

More information

Dynamics of a Superconducting Qubit Coupled to an LC Resonator

Dynamics of a Superconducting Qubit Coupled to an LC Resonator Dynamcs of a Superconductng Qubt Coupled to an LC Resonator Y Yang Abstract: We nvestgate the dynamcs of a current-based Josephson juncton quantum bt or qubt coupled to an LC resonator. The Hamltonan of

More information

2010 Black Engineering Building, Department of Mechanical Engineering. Iowa State University, Ames, IA, 50011

2010 Black Engineering Building, Department of Mechanical Engineering. Iowa State University, Ames, IA, 50011 Interface Energy Couplng between -tungsten Nanoflm and Few-layered Graphene Meng Han a, Pengyu Yuan a, Jng Lu a, Shuyao S b, Xaolong Zhao b, Yanan Yue c, Xnwe Wang a,*, Xangheng Xao b,* a 2010 Black Engneerng

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

More information

Odd/Even Scroll Generation with Inductorless Chua s and Wien Bridge Oscillator Circuits

Odd/Even Scroll Generation with Inductorless Chua s and Wien Bridge Oscillator Circuits Watcharn Jantanate, Peter A. Chayasena, Sarawut Sutorn Odd/Even Scroll Generaton wth Inductorless Chua s and Wen Brdge Oscllator Crcuts Watcharn Jantanate, Peter A. Chayasena, and Sarawut Sutorn * School

More information

Boise State University Department of Electrical and Computer Engineering ECE 212L Circuit Analysis and Design Lab

Boise State University Department of Electrical and Computer Engineering ECE 212L Circuit Analysis and Design Lab Bose State Unersty Department of Electrcal and omputer Engneerng EE 1L rcut Analyss and Desgn Lab Experment #8: The Integratng and Dfferentatng Op-Amp rcuts 1 Objectes The objectes of ths laboratory experment

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

Computational Modelling of the Unbalanced Magnetic Pull by Finite Element Method

Computational Modelling of the Unbalanced Magnetic Pull by Finite Element Method Avalable onlne at www.scencedrect.com Proceda Engneerng 48 (2012 ) 83 89 MMaMS 2012 Computatonal Modellng of the Unbalanced Magnetc Pull by Fnte Element Method Martn Donát a * a Brno Unversty of Technology,

More information

Resource Allocation with a Budget Constraint for Computing Independent Tasks in the Cloud

Resource Allocation with a Budget Constraint for Computing Independent Tasks in the Cloud Resource Allocaton wth a Budget Constrant for Computng Independent Tasks n the Cloud Wemng Sh and Bo Hong School of Electrcal and Computer Engneerng Georga Insttute of Technology, USA 2nd IEEE Internatonal

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

Canonical transformations

Canonical transformations Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,

More information

Chapter 9: Statistical Inference and the Relationship between Two Variables

Chapter 9: Statistical Inference and the Relationship between Two Variables Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,

More information

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors An Algorthm to Solve the Inverse Knematcs Problem of a Robotc Manpulator Based on Rotaton Vectors Mohamad Z. Al-az*, Mazn Z. Othman**, and Baker B. Al-Bahr* *AL-Nahran Unversty, Computer Eng. Dep., Baghdad,

More information

Estimating the Fundamental Matrix by Transforming Image Points in Projective Space 1

Estimating the Fundamental Matrix by Transforming Image Points in Projective Space 1 Estmatng the Fundamental Matrx by Transformng Image Ponts n Projectve Space 1 Zhengyou Zhang and Charles Loop Mcrosoft Research, One Mcrosoft Way, Redmond, WA 98052, USA E-mal: fzhang,cloopg@mcrosoft.com

More information

Formulas for the Determinant

Formulas for the Determinant page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use

More information

Field computation with finite element method applied for diagnosis eccentricity fault in induction machine

Field computation with finite element method applied for diagnosis eccentricity fault in induction machine Proceedngs of the Internatonal Conference on Recent Advances n Electrcal Systems, Tunsa, 216 Feld computaton wth fnte element method appled for dagnoss eccentrcty fault n nducton machne Moufd Mohammed,

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

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

Parameter Estimation for Dynamic System using Unscented Kalman filter

Parameter Estimation for Dynamic System using Unscented Kalman filter Parameter Estmaton for Dynamc System usng Unscented Kalman flter Jhoon Seung 1,a, Amr Atya F. 2,b, Alexander G.Parlos 3,c, and Klto Chong 1,4,d* 1 Dvson of Electroncs Engneerng, Chonbuk Natonal Unversty,

More information

Harmonic Detection Algorithm based on DQ Axis with Fourier Analysis for Hybrid Power Filters

Harmonic Detection Algorithm based on DQ Axis with Fourier Analysis for Hybrid Power Filters Harmonc Detecton Algorthm based on DQ Axs wth Fourer Analyss for Hybrd Power Flters K-L. AREERAK Power Qualty Research Unt, School of Electrcal Engneerng Insttute of Engneerng, Suranaree Unversty of Technology

More information

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

A Fast Computer Aided Design Method for Filters

A Fast Computer Aided Design Method for Filters 2017 Asa-Pacfc Engneerng and Technology Conference (APETC 2017) ISBN: 978-1-60595-443-1 A Fast Computer Aded Desgn Method for Flters Gang L ABSTRACT *Ths paper presents a fast computer aded desgn method

More information

JAB Chain. Long-tail claims development. ASTIN - September 2005 B.Verdier A. Klinger

JAB Chain. Long-tail claims development. ASTIN - September 2005 B.Verdier A. Klinger JAB Chan Long-tal clams development ASTIN - September 2005 B.Verder A. Klnger Outlne Chan Ladder : comments A frst soluton: Munch Chan Ladder JAB Chan Chan Ladder: Comments Black lne: average pad to ncurred

More information

Coupling Element and Coupled circuits. Coupled inductor Ideal transformer Controlled sources

Coupling Element and Coupled circuits. Coupled inductor Ideal transformer Controlled sources Couplng Element and Coupled crcuts Coupled nductor Ideal transformer Controlled sources Couplng Element and Coupled crcuts Coupled elements hae more that one branch and branch oltages or branch currents

More information

The optimal delay of the second test is therefore approximately 210 hours earlier than =2.

The optimal delay of the second test is therefore approximately 210 hours earlier than =2. THE IEC 61508 FORMULAS 223 The optmal delay of the second test s therefore approxmately 210 hours earler than =2. 8.4 The IEC 61508 Formulas IEC 61508-6 provdes approxmaton formulas for the PF for smple

More information

SRM control algorithm for electric vehicles drives applications

SRM control algorithm for electric vehicles drives applications Krzysztof TOMCZUK, Marcn PARCHOMIUK Electrotechncal Insttute, Power Converters Department SRM control algorthm for electrc vehcles drves applcatons Abstract. In the artcle author descrbes a measurements

More information

FTCS Solution to the Heat Equation

FTCS Solution to the Heat Equation FTCS Soluton to the Heat Equaton ME 448/548 Notes Gerald Recktenwald Portland State Unversty Department of Mechancal Engneerng gerry@pdx.edu ME 448/548: FTCS Soluton to the Heat Equaton Overvew 1. Use

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

Polynomial Regression Models

Polynomial Regression Models LINEAR REGRESSION ANALYSIS MODULE XII Lecture - 6 Polynomal Regresson Models Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Test of sgnfcance To test the sgnfcance

More information

FE REVIEW OPERATIONAL AMPLIFIERS (OP-AMPS)( ) 8/25/2010

FE REVIEW OPERATIONAL AMPLIFIERS (OP-AMPS)( ) 8/25/2010 FE REVEW OPERATONAL AMPLFERS (OP-AMPS)( ) 1 The Op-amp 2 An op-amp has two nputs and one output. Note the op-amp below. The termnal labeled l wth the (-) sgn s the nvertng nput and the nput labeled wth

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

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

Consistency & Convergence

Consistency & Convergence /9/007 CHE 374 Computatonal Methods n Engneerng Ordnary Dfferental Equatons Consstency, Convergence, Stablty, Stffness and Adaptve and Implct Methods ODE s n MATLAB, etc Consstency & Convergence Consstency

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

R. W. Erickson. Department of Electrical, Computer, and Energy Engineering University of Colorado, Boulder

R. W. Erickson. Department of Electrical, Computer, and Energy Engineering University of Colorado, Boulder R. W. Erckson Department of Electrcal, Computer, and Energy Engneerng Unersty of Colorado, Boulder 3.5. Example: ncluson of semconductor conducton losses n the boost conerter model Boost conerter example

More information

Grover s Algorithm + Quantum Zeno Effect + Vaidman

Grover s Algorithm + Quantum Zeno Effect + Vaidman Grover s Algorthm + Quantum Zeno Effect + Vadman CS 294-2 Bomb 10/12/04 Fall 2004 Lecture 11 Grover s algorthm Recall that Grover s algorthm for searchng over a space of sze wors as follows: consder the

More information

Lecture 5 Decoding Binary BCH Codes

Lecture 5 Decoding Binary BCH Codes Lecture 5 Decodng Bnary BCH Codes In ths class, we wll ntroduce dfferent methods for decodng BCH codes 51 Decodng the [15, 7, 5] 2 -BCH Code Consder the [15, 7, 5] 2 -code C we ntroduced n the last lecture

More information

FFT Based Spectrum Analysis of Three Phase Signals in Park (d-q) Plane

FFT Based Spectrum Analysis of Three Phase Signals in Park (d-q) Plane Proceedngs of the 00 Internatonal Conference on Industral Engneerng and Operatons Management Dhaka, Bangladesh, January 9 0, 00 FFT Based Spectrum Analyss of Three Phase Sgnals n Park (d-q) Plane Anuradha

More information