THE steady-state operation of rotating electrical machines
|
|
- Samantha King
- 5 years ago
- Views:
Transcription
1 IEEE TRANSACTIONS ON MAGNETICS, DATE 1 Generalized Slip Transformaions and Air Gap Harmonics in Field Models of Elecrical Machines Ville Räisänen, Saku Suuriniemi, Lauri Keunen Absrac In numerical field analysis of seady-saes of elecrical machines, frequency-domain mehods are ofen much faser han ime-domain mehods. Frequency-domain echniques ofen involve simplificaions ha impac heir accuracy on harmonic effecs due o sloing, sauraion and ime harmonics in winding currens. The naure of hese inaccuracies and heir relaionship o air gap field harmonics is no ha well covered in he lieraure. Correc predicion of air gap field harmonics can be based on careful use of specral Dirichle-o-Neumann (D-o-N) mappings and generalized slip ransformaions (GST). We show how nonzero harmonics in he air gap fields can be reliably prediced and explain he naure of inaccuracies in cerain common seady-sae soluion echniques. Index Terms Boundary value problems, Fourier series, harmonics, roaing machines. I. INTRODUCTION THE seady-sae operaion of roaing elecrical machines can be prediced in a number of ways. Approximae soluion by Finie Elemen Mehods (FEM) in ime domain is a sandard approach, bu full soluion of hese compuaionally inensive problems is no always feasible. Then frequencydomain or saic mehods are ofen used. Many moor analyses depend on a mehod ha accouns for he relaive movemen beween he roor and he saor. In he case of inducion machines, ofen only a single imefrequency is used in he roor equaions. This roor frequency is he single saor frequency muliplied by he roor slip he simples example of a slip ransformaion. The res of problem is solved as usual a he synchronous frequency 1]. While he single slip ransform approach may provide a useful firs orque esimae, i is inadequae for mos oher asks: Inverers higher harmonics and he roaion of a sloed roor induce addiional ime-frequencies ino he saor fields. When he roor slip is no an ineger, hese ime-frequencies are no always ineger-muliples of he synchronous frequency 5]. The fundamenal slip ransformaion is incorrec for mos oher space and ime harmonics 2], 3]: boh he ime and space harmonic indices influence he appropriae slip ransformaion, as found in previous work 4] on seady-sae problems. The generalized slip ransformaions (GST) arise from his need. GSTs describe how he ime-frequencies of he air gap field harmonics ransform beween he roor and saor coordinae sysems in relaive moion. They are he key o accuraely accoun for he air gap field harmonics in wo-dimensional V. Räisänen was wih he Deparmen of Elecrical Engineering, Tampere Universiy of Technology, Tampere 3311, Finland. S. Suuriniemi and L. Keunen are wih he Deparmen of Elecrical Engineering, Tampere Universiy of Technology, Tampere 3311, Finland. Corresponding auhor: Saku Suuriniemi ( saku.suuriniemi@u.fi) frequency-domain BVPs for elecrical machines. The applicaion of he ransformaion needs aenion o deail, no well covered in he lieraure on frequency-domain soluion mehods. The roor and he saor can be seen o respond o he air gap field harmonics. Cerain air gap sream funcion harmonic excies a unique combinaion of harmonics in he normal derivaive of he sream funcion. This is reducible o air gap harmonic coefficiens, and compleely described in erms of specral Dirichle-o-Neumann (D-o-N) maps 7], 8], 9]. The marix represenaions of linear specral D-o-N maps exhibi sparsiy paerns ha sem from he roaional symmery of sloed geomeries. The sparsiy paerns, ogeher wih he GSTs, explain how roor and saor subproblems creae harmonics ino he air gap fields 8]. The non-zero air gap field harmonics can be hen accuraely prediced. The specral D- o-n mappings and GSTs also provide a solid foundaion for accurae and efficien frequency-domain mehods. Prior frequency-domain echniques ha he auhors know of suffer from over-simplified slip ransformaions. The permeance wave heory, a popular air gap harmonics analysis mehod, describes he effec of slos by a permeance relaion beween he magneomoive force over he air gap and he air gap flux densiy 1]. While some permeance relaions can be derived direcly from simple boundary value problems, hey rely on very resricive assumpions abou slo shapes, curren densiies, and maerial relaions 6], 1], 11]. In conras, he sparsiy paerns for specral D-o-N mappings follow formally from roaional symmery of maerial parameers in sloed geomeries and apply in any such BVP, regardless of he shape or maerial parameers of he slos. The Muli-Layer Transfer Marix (MLTM) mehods build on analyical soluions of roor subproblems in erms of individual space harmonics 12]. However, since he higher space harmonics in MLTM are due o he space harmonics in he curren shee expansion of he saor winding currens, he saor models in MLTM fail o accoun for he influence of he saor sloing 4], 13]. The roor and saor subproblems can be coupled also in FE models by a space harmonic analyical model of he air gap 14]. Such an air gap model can be inerpreed as a specral D-o-N mapping from he viewpoin of he roor and he saor subproblems 15]. In his paper, he specral D-o-N mappings of he roor and saor subproblems are expressed in co-moving coordinae sysems, and hey do no couple differen ime-frequencies. Consequenly, he effec of moion is isolaed o he slip ransformaions.
2 IEEE TRANSACTIONS ON MAGNETICS, DATE 2 Ω s Ω r Figure 1. Half of he problem geomery. The air gap widh is exaggeraed. Λ s Λ s Λ s θ Γ dr Λ r Γ sr Λ r Γ ds Λ r θ Figure 2. An example of he ransformaion of harmonics beween he saor and he roor coordinae sysems. II. AIR GAP FIELD HARMONICS Le us denoe by Ω he cross secion of a machine. This is furher divided ino he saor and he roor regions Ω s and Ω r, whose common inerface Γ sr = Ω s Ω r is a he cener of he air gap (see Fig. 1). The ouer boundary of he saor is denoed by Γ ds and he inner boundary of he roor region by Γ dr. We assume ha he magneic flux densiy has no componen in he axial direcion and ha he curren densiy is purely axial. An A z -formulaion wih a sream funcion A corresponding o he magneic vecor poenial is used 9]. As in publicaion 9], magneic sream funcions A s and A r are used in he regions Ω s and Ω r in heir respecive comoving coordinae sysems. The change of coordinaes from he saor o he roor sysem is u : R 2 R 2 : rs θ s ] rr θ r ] := r s θ s (1 s) ωs p where s is he roor slip, p is he number of pole pairs and ω s is he synchronous frequency. If A is he magneic sream funcion in he whole problem domain expressed in he saor coordinae sysem, he sream funcion expressed in he roor coordinae sysem is A r := A Ωr u. A. Dirichle and Neumann daa The magneic sream funcion and is radial derivaive on Γ sr are expanded in boh coordinae sysems as a Fourier series expansion of wo variables: ime and respecive mechanical angle θ k. Le us denoe he se of non-zero air gap harmonics in he subproblem k {s,r} by Λ k and index he harmonics in Λ k wih naural numbers. In order o faciliae he use of marix algebra and numerical compuaions, he series expansions for boh subproblems k mus be runcaed ] (1) o respecive combinaions over finie ses of harmonics. Tha is, we wrie Λ k := { ( k n, k n) : n = 1,..., Λ k } Z 2, (2) where n k and n k are he ime and space harmonic indices of he n:h non-zero harmonic of he subproblem k. The selecion of he ses Λ k will be discussed in he Subsecion IIIc. The se of space and ime harmonics in he se Λ k are denoed by Λ k θ and Λk, respecively. The se of space harmonics Λ k θ is independen of he coordinae sysem whereas he relaion beween Λ s and Λ r is dependen on he roor slip (see Fig. 2). We assume ha he magneic sream funcions in boh coordinae sysems are ime-periodic a some fundamenal imefrequency ω f. An appropriae ω f can be always deermined, and is choice is discussed in secion IIID. The ime-harmonic index of he synchronous frequency is denoed by h = ω s /ω f and he harmonics ±(h, p) are called he working harmonics. Thus, we expand A k (θ k,) = D k e j( ω f θ k ) (3) (,) Λ k A k r (θ k,) = Nη k ηe j(η ω f ηθ k ). (4) (η,η) Λ k The expansions (3) (4) are called he Dirichle and Neumann daa in coordinae sysem k, respecively. The coefficiens D k and Nk η η are he Dirichle and Neumann coefficiens, respecively. Since A k and is radial derivaive mus be real, conjugae symmery D( k, ) = Dk (,) mus hold for all (,) Λ k, where denoes he complex conjugae. Le us define inner producs, θ and, on θ,2π] and,2π/ω f ], respecively, so ha f,g θ := 1 2π f,g := ω f 2π 2π 2π/ωf f(θ)g(θ) dθ (5) f()g() d (6) hold. Thereafer, he Dirichle and Neumann coefficiens can be obained from he Dirichle and Neumann daa (3) (4) Ak D k = k,e jθ ω f (7) θ,ej Nη k A k η = r,e jηθ k,e jη ω f. (8) For convenience, he harmonic indices for he harmonics in he Dirichle and Neumann daa are denoed wih differen symbols so ha ( k n, k n) = (η k n, η k n) (9) hold for all n = 1,2,..., Λ k. Then, he coefficiens are assembled ino he vecors ] D k = D k, D k... D k T 1 kk 1 2 kk 2 k (1) Λ k Λ ] N k = N k, N k... N k T η 1 kηk 1 η 2 kηk 2 η k Λ ηk Λ. (11) θ
3 IEEE TRANSACTIONS ON MAGNETICS, DATE 3 Roor Slip s (, 37) (, 35) (, 5) (, ) (1,1) (1,5) (, 7) (1,7) Opposie (1, 35) (1, 37) Locked Synchronous Generalized Slip s 1 Figure 3. Slip ransformaions for he ime-frequencies in he Dirichle daa for some harmonics (,) when ω f = ω s and p = 1. Change in he ime harmonic number corresponds o a horizonal shif in he generalized slip s p. B. Slip Transformaions and Inerface Condiions When he roor is no locked, he ime-frequencies of he space harmonics in he Dirichle and Neumann daa are no equal in differen coordinae sysems. The relaionships beween he ime-frequencies of he fields of he roor and he saor subproblems are known as slip ransformaions, which we shall derive here. Equivalen ransformaions have been widely discussed in he lieraure 2], 4], 5]. Slip ransformaions arise from he inerface condiions on Γ sr for he Dirichle and Neumann daa: ( A s θ s + 1 s ) p ω s, = A r (θ s,) (12) A s r ( θ s + 1 s ) p ω s, = A r r (θ s,). (13) Dirichle daa in he roor coordinae sysem can be obained by he subsiuion of (3) ino (12) ( A r = A s θ r + 1 s ) p ω s (14) = D s expj(s p ω f θ r ) (,) Λ s where we have defined he generalized slip s p := ω s ω f 1 s p. (15) In Fig. 3, he ransformaion (15) is demonsraed for a se of harmonics. From he inerface condiions (12) (13), i follows ha he number of non-zero air gap harmonics is independen of he coordinae sysem. The harmonics of differen subproblems are relaed by D s = D r s p (16) N s = N r s p. (17) The ses Λ r and Λ s used for he expansions (3) (4) should be seleced so ha(,) Λ s holds if and only if(s p,) Λr holds. The Dirichle and Neumann coefficiens for boh subproblems k are assembled ino vecors D k,n k C Λk so ha roor and saor coefficiens relaed by (16) (17) share common rows. Thereafer, (16) (17) are equivalen o C. Winding Currens D s = D r (18) N s = N r. (19) The curren densiy in each subproblem can consis of winding currens and eddy-currens in massive conducors. Tha is, he curren densiy is expanded as J k = Q k q=1 Λ k γ= Q k /2 I k qj q k ej ω f σ A k (2) where J q k is he source curren densiy corresponding o uni curren in he winding of he slo q in he subproblem k. For simpliciy, we assume ha slos conain only a single layer. Wihou loss of generaliy, we assume ha he slos are ordered in ani-clockwise direcion according o mechanical angle. Thereafer, we expand he slo winding currens as a discree Fourier series Q k /2 ( I k q = O k γ exp j 2πγq ) (21) Q k where O γ is called he harmonic (,γ) in he slo winding currens of subproblem k. The equaion (21) can wrien in he marix form I k = ζ k] O k. (22) D. Specral Dirichle-o-Neumann Mappings From he viewpoin of he inerface Γ sr a he cener of he air gap, he elecromagneics of each subproblem is compleely characerized by a mapping from he Dirichle daa o he Neumann daa. Such mappings are known as Dirichle-o- Neumann (D-o-N) mappings. We call represenaions of D-o- N mappings in erms of Dirichle and Neumann coefficiens specral D-o-N mappings 8], 9]. They compacly express e.g. differen roor varians in a design siuaion. If he fundamenal frequency ω f and he ses Λ r,λ s have been chosen appropriaely, he specral D-o-N mapping for he subproblem k in seady sae can be expressed as a funcion N k = S(D k,i k ). (23) If sauraion can be negleced, here exis marices B k] C Λk Λ k and F k] C Λk Q k such ha N k = B k] D k + F k] I k (24) holds (see Appendix A). The marix B k characerizes he behavior of he subproblem in he absence of winding currens. The marix F k characerizes he behavior of he subproblem when no flux crosses Γ sr. Combinaion of (24) wih he inerface condiions (18) (19), yields B s B r] D s = F r ζ r] O r F s ζ s] O s. (25)
4 IEEE TRANSACTIONS ON MAGNETICS, DATE 4 III. PREDICTION OF AIR GAP HARMONICS This secion explains how he sparsiy paerns for he marices B k] and F k ζ k] ogeher wih he GSTs (15) can applied o (25) in order o predic he air gap field harmonics, which can be non-zero. A. Sloed Geomeries and Sparsiy Paerns When sauraion can be negleced, he maerial parameers of sloed geomeries of he roor and he saor subproblems are roaionally symmeric wih respec o roaions by a single slo. Le e k be he sandard basis vecor corresponding o he harmonic(,) in he subproblemk. The roaional symmery leads o sparsiy paern in he marix B k] such ha ] B k e k = ] B k (η,η)(,) ek η η (26) (η,η) Z 2 = n Z B k ] (,+nq k )(,) ek (+nq k ) (,) (s p,) Dirichle (s p,) (,) Dirichle D-o-N D-o-N ( +n r, +nq r ) η (s p, +nq r) Neumann (s p, +nq s) ( n s, +nq s ) Neumann Figure 4. The influence of roor moion o he specral Dirichle-o-Neumann map for he roor (above) and he saor (below) subproblems from he viewpoin of he oher subproblem. η η η holds for every (,) 8]. The erms in he series (26) ouside Λ k are runcaed. Tha is, he marix B k] maps harmonic (,) in he Dirichle daa o he harmonics {(η,η) = (, +nq k ) : n Z} in he Neumann daa. In a similar way, each harmonic (,γ) in he slo winding currens is mapped by he marix F k ζ k] o he harmonics {(,γ +nq k ) : n Z} in he Neumann daa. Tha is, F k ζ k] e k γ = F k ζ k] (,γ+nq k )(,γ) ek (γ+nq k ) (27) holds. n Z B. Coordinae Sysems and Movemen The sparsiy paerns (26) (27) imply ha specral D- o-n maps consruced for he roor and saor subproblems in heir own coordinae sysems do no couple differen ime-frequencies in he Dirichle and Neumann daa. Insead, differen ime-frequencies become coupled due o he relaive movemen of he roor wih respec o he saor. In order o undersand his coupling, le us invesigae he combined effec of he slip ransformaions and he specral D- o-n mapping for he roor from he viewpoin of he saor. Consider a single harmonic (,) in he saor Dirichle daa. In he roor coordinae sysem, he corresponding harmonic is (s p,), where he change in he ime-frequency depends on he space harmonic index according o (15). In Fig. 4, his change in ime-frequency is depiced as a verical shif from he saor line o he roor line. From (26), i follows ha he specral D-o-N mapping for he roor maps he harmonic (s p,) in he roor Dirichle daa o he harmonics {(s p, +n rq r ) : n r Z} in he roor Neumann daa. The corresponding harmonics in he saor sysem are obained wih (15) s p (s p )( (+n rq r)) = s p +( +n rq r ) ω s 1 s ω f p = +n r r (28) r s Q r +( r,q r Q s) +( r,q r) +( r,q r +Q s) +( r,q r +2Q s) (,Q s) (,γ) +(,Q s) +(,2Q s) ( r,q r) ( r,q r Q s) ( r,q r Q s) Q s +( s, Q r Q s) +( s, Q s) +( s,q r Q s) +( s,2q r Q s) (,Q r) (,γ) +(,Q r) +(,2Q r) +( s,q s Q r) Q s +( s,q s) Q r +( s,q r +Q s) +( s,2q r +Q s) Figure 5. Equivalence classes,γ] s and,γ] r of he harmonic (,γ) in he slo winding currens in he saor (above) and roor (below) coordinae sysems, respecively. where we have defined he laice cell heigh k := Q k ω s ω f 1 s p. (29) Thus, he combined effec of he specral D-o-N map for he roor and he relaive roaion can couple differen imefrequencies in he saor subproblem. In a symmeric way, from he viewpoin of he roor, differen ime-frequencies can be coupled due o he combined effec of relaive roaion and he specral D-o-N mapping for he saor (see Fig. 4). C. Dirichle and Neumann Laices Expansion of he row (η,η) of (25) in erms of he saor Dirichle daa wih(28) and applicaion of he sparsiy paerns
5 IEEE TRANSACTIONS ON MAGNETICS, DATE 5 (26) (27), yields ] B s (η,η)(η,η+nq Ds s) (η,η+nq s) (3) n Z Λ k n Z B r ] (s p η η,η)(sp η η,η+nqr)ds (η +n r,η+nq r) = n Z F s ζ s] (η,η)(η,η+nq s) Os (η,η+nq s) n Z F r ζ r] (s p η η,η)(sp η η,η+nqr)or (s p η η,η+nqr). The sysem of equaions (3) is linear wih respec o he source curren vecors O s and O r. Thus, he Dirichle coefficien vecor D s solved from (3) wih O s and O r is he sum of he Dirichle coefficien vecors solved from (3) wih individual rows of he source curren vecors. The se of non-zero Dirichle coefficiens is herefore a subse of he union of non-zero Dirichle coefficiens obained from (3) wih individual winding curren harmonics. If O s = e s γ and Or = hold, i follows from (27) ha on he righ-hand side of (3), only he rows (η,η) = (,γ +nq s ) : n s Z can be non-zero. Each row (η,η) on he lef-hand side of (3) is dependen on he saor Dirichle coefficiens (η,η + n s Q s ) : n s Z and (η + n r r,η + n r Q r ) : n r Z. Thus, he harmonic (,γ) in he saor slo winding currens and he harmonic (s p γ,γ) in he roor winding currens is coupled only o he harmonics,γ] s := {( +n r r, γ +n s Q s +n r Q r ) : n s,n r Z} (31) in he saor Dirichle daa, where he laice cell heigh r = ω Q s 1 s r ω f p (see Fig. 5). In a similar way, harmonic (,γ) in he roor slo winding currens and he harmonic (s p γ,γ) in he saor winding currens are coupled o he harmonics,γ] r := {( n s s, γ +n r Q r +n s Q s ) : n s,n r Z} (32) in he roor Dirichle daa, where he laice cell heigh s = ω Q s 1 s s ω f p (see Fig. 5). The same harmonics can be non-zero in he Neumann daa. Harmonics in he Dirichle and Neumann daa can be pariioned ino equivalence classes (31) (32). Le us denoe he se of non-zero winding curren harmonics in he subproblem k wih Λ k I. Harmonics in each equivalence class,γ] k form a laice (see Fig. 5) and he ses of non-zero harmonics Λ s and Λ r saisfy Λ s,γ] s s p ( γ),γ] s (33) Λ r respecively. (,γ) Λ s I (,γ) Λ r I,γ] r (,γ) Λ r I (,γ) Λ s I s p γ,γ] r (34) Figure 6. Pars of equivalence classes for wo coil curren ime harmonics (black and whie nodes). Time harmonics in orque arise for nodes wih equal space harmonic indices (a few illusraed by dashed lines). Heigh of he line corresponds o he ime harmonic index γ of he orque (righ). ω f should be chosen o make boh h = ω s /ω f and h(1 s)/p inegers in(15) heoreically possible for raional values of slip only. However, any slip can be approximaed o arbirary precision by a raional number. When he roor slip is an ineger and he number of slos in boh subproblems is divisible by he number of pole pairs, k Z holds for boh subproblems k and he synchronous frequency is suiable as he fundamenal frequency. E. Immediae Applicaions The air gap harmonics uniquely correspond o roor and saor fields. Therefore, loss esimaes can be compued from he Dirichle coefficiens and he soluions of he field problems ha are needed o consruc he D-o-N map 4]. The laice (33) (34) of non-zero air gap harmonics relaes simply o he orque ime harmonics. Le r δ be he radius of he circle Γ sr and L he lengh of he acive area of he roor. Subsiuion of (3) (4) o he orque formula obained wih he Maxwell Sress Tensor in 9] is T = Lr 2π δ D k µ N k η η e j( +η )ω f (+η)θ k dθ k (, )(µ,µ ) (35) The inegral is nonzero only for η =, and applicaion of he complex conjugae propery N η k η = Nk η η yields T = 2π Lr δ µ ( ) η Λ k D k N k η e j( +η )ω f. (36) The ime harmonic coefficien of index γ in he orque is T,e jγ ω f = 2π Lr δ D k µ N ( k γ ). (37) (,) Λ k Thus, orque harmonics arise from Dirichle and Neumann coefficiens wih equal space harmonic indices. The orque harmonic index γ = η is he difference of he ime harmonic indices of he coefficiens (see Fig. 6). If, for example, a orque frequency componen causes an objecionable noise, one can idenify he air gap harmonics ha cause i. D. Selecion of he Fundamenal Frequency In order o Fourier series expansions (3) (4) of he Dirichle and Neumann daa o exis, he fundamenal frequency ω f has o be seleced so ha he ime-frequencies of all non-zero harmonics are ineger-muliples of ω f. IV. EXAMPLES This secion applies he framework discussed in secions II- III o wo realisic examples. Fourier series expansion (21) of he winding currens is used o predic a se of laices. The unions (33) (34) of hese laices are compared o a Fourier
6 IEEE TRANSACTIONS ON MAGNETICS, DATE 6 log 1 D s, s= Hz 42 Hz Hz 335 Hz Hz 25 Hz Hz 165 Hz Hz 8 Hz 1 5 Hz 5 Hz log 1 D s, s= Hz Hz Hz Hz Hz Hz Hz Hz Hz Hz 1 5 Hz 5 Hz log 1 D s, s = 1 log 1 D r, s = 9 45 Hz Hz Hz Hz 18 9 Hz Hz log 1 D r, s= Hz 445 Hz Hz 3235 Hz Hz 2425 Hz Hz 1615 Hz Hz Hz 1 5 Hz 5 Hz log 1 D r, s = Hz 5 Hz Hz 5 Hz Figure 7. The Dirichle coefficiens prediced and compued from a ime-domain FE soluion. Solid lines are drawn for he he equivalence classes ±h,p] s and ±s p hp,p]r of he working harmonic ±(h,p) in he saor and roor coordinae sysems, respecively. The remaining equivalence classes in Λs I and Λr I are drawn wih dashed lines. Circles are drawn for he Dirichle coefficiens compued from he FE soluion. Circles darkness correspond o ampliudes of he harmonics wihin he displayed range. series coefficiens exraced from ime-domain FE soluions wih (7) (8). In all es cases, 4-pole inegral-slo saor windings are driven wih balanced 3-phase currens. Thus, he slo winding currens conain harmonics 8] Λ s I = {±(h,2),±(h,14),±(h,)}. (38) By subsiuion of (38) ino (33) (34), we obain he predicions Λ s ±h,2] s ±h,14] s ±h,] s (39) Λ r ±s p h2,2] r ±s p h14,14] r ±s p h,] r. (4) A. Squirrel-Cage Inducion Moor Consider a 4-pole squirrel-cage inducion moor wih Q s = 36 saor slos and Q r = 34 roor slos depiced in Fig. 1. The es case is a simplified version of he moor presened in 17]. A no-load, he slip s = and we selec ω f := ω s. Thereafer, h = 1, r = 17, s = 18 and h,] s = {(h+17n r, +34n r +36n s ) : n r,n s Z} (41) h,] r = {(h8n s, +34n r +36n s ) : n r,n s Z} (42) hold. A s =.1, we selec ω f := ω s /1. Thereafer, h = 1, r = 153, s = 162 and h,] s = {(h+153n r, +34n r +36n s ) : n r,n s Z} (43) h,] r = {(h62n s, +34n r +36n s ) : n r,n s Z} (44) hold. When he roor is locked, he slip s = 1 and we selec ω f := ω s. Thereafer, h = 1, r = s = and h,] r = h,] s = {(h, +2n) : n Z} (45) hold. In Fig. 7, he predicions (38) (45) are compared o resuls obained wih (7) (8) from a ime-domain FE simulaion using he code discussed in 9] wih he assumpion of linear maerials. B. Synchronous Relucance Moor Consider a 4-pole synchronous relucance moor wih Q s = 36 saor slos and Q r = 4 roor poles depiced in Fig. 8. Since he roor roaes a he synchronous speed, ω f = ω s, r = Q r /p = 2 and s = Q s /p = 18 hold. Thus,,] s = {( +2n r, +4n r +36n s : n r,n s Z)} (46),] r = {( 8n s, +4n r +36n s : n r,n s Z)} (47) hold. In Fig. 9, he predicions (46) (47) are compared o resuls obained wih (7) (8) from a ime-domain FE simulaion using he code discussed in 9]. V. INACCURACIES IN PRIOR METHODS This secion uses he approach of Secions II-III o gain insigh ino inaccuracies in wo formulaions used wih prior
7 IEEE TRANSACTIONS ON MAGNETICS, DATE 7 d-axis Figure 8. Geomery of he synchronous relucance machine wih he roor and he saor drawn separaely. log 1 D s Hz Hz Hz Hz Hz Hz Hz Hz Hz Hz 9 45 Hz 7 35 Hz 5 25 Hz 3 15 Hz 1 5 Hz 1 2 log 1 D r Hz 3 15 Hz Hz 18 9 Hz 12 6 Hz 6 3 Hz Hz 2 4 Figure 9. Dirichle coefficiens prediced wih (46) (47) and compued from a ime-domain FE soluion for he synchronous relucance machine wih (7) (8). Lines and circles are used as in Fig. 7. seady-sae compuaional mehods. Boh cases follow he approach used in Secion IIIB. The formulaion for he roor subproblem is inerpreed as a specral D-o-N map and he coupling echnique used o couple he roor and saor subproblems is inerpreed as wo disinc slip ransformaion rules for he Dirichle and Neumann daa, respecively. By comparison o he correc consrucion of he specral D-o-N maps and slip ransformaions developed in Secions II-III, we can undersand error sources in he discussed mehods. A. Single-Slip Approach The mos widely applied approach for inducion machines is presened in 1]. The problem is formulaed as a imeharmonic eddy-curren problem a he synchronous frequency. To accoun for he relaive movemen of he roor wih respec o he saor, he ime-frequency used for he roor equaions is muliplied wih he slip. Muliplicaion of he ime-frequency used in he roor subproblem wih he slip corresponds o applicaion of he same slip ransformaion s p hp = s o all space harmonics ±(h,) in he Dirichle daa of he saor subproblem. In Fig. 1, his is depiced as a verical shif from he saor line o he roor line. Thus, if he slip is small, slo harmonics wih large space harmonic indices are given he usually smaller ime-frequency s p hp ω f = sω s of he working space harmonic = p. This can lead o severe inaccuracies in eddy-curren loss esimaes 4]. The roor subproblem can be inerpreed as a specral D-o- N map, which maps harmonics (sh,) in he roor Dirichle daa o harmonics (sh, + nq r ) in he roor Neumann daa (see Fig. 1). Coninuiy of he Neumann daa on Γ sr implies ha hese harmonics are ransformed o harmonics (h, + nq r ) in he saor Neumann daa wih he slip ransformaion appropriae only for (sh, p). (h,) (s p hp,) D-o-N Dirichle η (h, +nq r) (s p hp, +nqr) η Neumann Figure 1. Inerpreaion of he common approach of using a single slip ransformaion. B. One-Way Generalized Slip In 2] 4], separae ime-harmonic BVPs are formulaed for differen space harmonics in he roor subproblem wih he appropriae slip ransformaion of ime-frequency used in each roor BVP. Only he synchronous frequency is considered in he saor subproblem. The saor Dirichle daa is expanded as a Fourier series wih respec o he mechanical angle and differen space harmonics are coupled o appropriae roor models. Each roor BVP can be inerpreed as par of a specral D-o-N map for he roor, which maps harmonic(s p h,) in he roor Dirichle daa o he harmonics {(η,η) = (s p h, +n rq r ) : n r Z} in he roor Neumann daa (see Fig. 11). In order o saisfy he coninuiy of he Neumann daa on Γ sr, Neumann daa for he saor subproblem is assumed o be equal o he sum of he Neumann daa from individual roor models. However, his corresponds o he assumpion ha he slip ransformaion for n r = is appropriae for all n r Z (compare Fig. 11 and Fig. 4). Dirichle (h,) (s p h,) D-o-N Neumann η (h, +nq r ) Figure 11. Inerpreaion of he approach discussed in 2] 4]. (s p h, +nq r) η Thus, all space harmonics are coupled back o he saor model a he synchronous ime-frequency. This is incorrec for he new space harmonics wihn r, creaed o he Neumann daa due o sloing of he roor. In ime-domain simulaions of inducion machines, effecs of hese ime harmonics creaed by roor sloing are ofen observed as addiional ime harmonics in he saor winding currens and as oscillaions in he orque. VI. CONCLUSION Analysis of air gap field harmonics in seady-sae analysis of wo-dimensional elecromagneic BVPs for elecrical machines was presened. From he viewpoin of he air gap, he
8 IEEE TRANSACTIONS ON MAGNETICS, DATE 8 roor and he saor responses can be compleely represened by specral Dirichle-o-Neumann maps in co-moving coordinae frames. The effec of moion is hen isolaed o he GSTs ha describe how ime-frequencies ransform beween he roor and saor subproblems. The main limiaion of he presened framework is he assumpion of linear maerial relaions. Nonnegligible sauraion in a subproblem requires a nonlinear specral D-o-N mapping, and he convenien lineariy of equaions (24) (25) is los. Addiional sauraion harmonics hen also appear. The roaional symmery of maerial parameers in sloed geomeries leads o saic sparsiy paers in he marix represenaions of he specral D-o-N maps. The analysis exhausively explains which air gap and orque ripple harmonics arise and enables efficien and generally applicable frequencydomain compuaion mehods. APPENDIX A CONSTRUCTION OF SPECTRAL D-TO-N MAPPINGS Consider he boundary value problem: Find A k defined in Ω k ha saisfies µ A k = σ A k Q k I k qj q ω f k ej Λ k q=1 A k = D k e j( ω f θ k ) (,) Λ k (48) (49) A k Γdk = (5) wih σ he elecric conduciviy and µ he magneic permeabiliy. The soluion o (48) (5) is expanded as he sum A k (r,) = + Λ k q=1 Q k I k qa q ω f ik (r)ej (,) Λ k D k à dk (r)e j ω f where A q ik for each q = 1,2,..,Q k is he soluion o (51) µ A q ik = J k q (52) A q ik Γ sr = (53) A q ik Γ dk = (54) and à dk for each (,) Λ k is he soluion o µ à dk = j ω f σã dk (55) à dk = e jθ k (56) =. (57) Γdk à dk The radial derivaive of A k on Γ sr can be expanded as A k = Q k I k A q ik r q r e j ω f + Λ k q=1 (,) Λ k D k à dk r e j ω f. (58) Le δ denoe he Kronecker dela. Applicaion of he inner producs(5) (6) o(58) and definiion of he marices B k] C Λk Λ k and F k] C Λk Q k such ha (24) follows: B k ] (η,η)(,) F k ] (η,η)(,q) à dk,e jηθ k r A q ik := δ η r,e jηθ k := δ η REFERENCES θ θ (59) (6) 1] E. Vassen, G. Meunier, A. Foggia, Simulaion of Inducion Machines Using Complex Magneodynamic Finie Elemen Mehod Coupled Wih he Circui Equaions, IEEE Trans. Magn., vol. 27, no. 5, pp , ] H. De Gersem, K. Hameyer, Air-Gap Flux Spliing for he Time- Harmonic Finie-Elemen Simulaion of Single-Phase Inducion Machines, IEEE Trans. Magn., vol. 38, no. 2, pp , 22. 3] Y. Ouazir, N. Takorabe, R. Ibiouen, O. Touhami, S. Mezani, Consideraion of Space Harmonics in Complex Finie Elemen Analysis of Inducion Moors Wih an Air-Gap Inerface Coupling, IEEE Trans. Magn., vol. 42, no. 4, pp , 26. 4] V. Räisänen, S. Suuriniemi, S. Kurz, L. Keunen, Rapid Compuaion of Harmonic Eddy-Curren Losses in High-Speed Solid-Roor Inducion Machines, IEEE Trans. Energy Convers., vol. 28, no. 3, pp , Sep ] K. Oberrel, Die Oberfeldheorie des Käfigmoors uner Berücksichigung der durch die Ankerrückwirkung verursachen Saorobersröme und der parallelen Wicklungszweige, Archiv für Elekroechnik, vol. 49, no. 6, pp , ] K. Oberrel, Losses, orques and magneic noise in inducion moors wih saic converer supply, aking muliple armaure reacion and slo openings ino accoun IET Elecric Power Applicaions, vol. 1, issue 4, pp , 27. 7] D. Givoli, Recen advances in he DN FE Mehod, Arch. Compu. Meh. Eng. vol. 6, no. 2, pp , ] V. Räisänen, Air Gap Fields in Elecrical Machines: Harmonics and Modeling of Movemen, PhD hesis, Tampere Universiy of Technology, ] V. Räisänen, S. Suuriniemi, S. Kurz, L. Keunen, Subdomain Reducion by Dirichle-o-Neumann Mappings in Time-Domain Elecrical Machine Modeling, IEEE Trans. Energy Convers., vol. 3, no. 1, pp , ] B. Heller, V. Hamaa, Harmonic Field Effecs in Inducion Machines, 1977: Elsevier. 11] W. Gibbs, Conformal Transformaions in Elecrical Engineering, Chapman & Hall, London, ] E. Freeman, Travelling waves in inducion machines: inpu impedance and equivalen circuis, Proceedings of he IEEE, vol. 115, no.12, ] J. Gieras, Analysis of mulilayer roor inducion moor wih higher space harmonics aken ino accoun, IEE Proceedings B: Elecric Power Applicaions, vol. 138, no. 2, pp , ] A. Razek, J.L. Coulomb, M. Feliachi, J. Sabonnadiere, Concepion of an air-gap elemen for he dynamic analysis of he elecromagneic field in elecric machines, IEEE Trans. Magn., vol. 18, no. 2, pp , ] H. De Gersem, T. Weiland, A compuaionally efficien air-gap elemen for 2D FE machine models, IEEE Trans. Magn., vol. 41, no. 5, pp , ] R. Merens, U. Pahner, K. Hameyer, R. Belmans, Force Calculaion Based on a Local Soluion of Laplace s Equaion. 17] A. Arkkio, Analysis of inducion moors based on he numerical soluion of he magneic field and circui equaions,, Ph.D. disseraion, Helsinki Univ. Technol., Helsinki, Finland, 1987.
Computation of the Effect of Space Harmonics on Starting Process of Induction Motors Using TSFEM
Journal of elecrical sysems Special Issue N 01 : November 2009 pp: 48-52 Compuaion of he Effec of Space Harmonics on Saring Process of Inducion Moors Using TSFEM Youcef Ouazir USTHB Laboraoire des sysèmes
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More information23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes
Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals
More informationd 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3
and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or
More informationTime Domain Transfer Function of the Induction Motor
Sudies in Engineering and Technology Vol., No. ; Augus 0 ISSN 008 EISSN 006 Published by Redfame Publishing URL: hp://se.redfame.com Time Domain Transfer Funcion of he Inducion Moor N N arsoum Correspondence:
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More information4. Electric field lines with respect to equipotential surfaces are
Pre-es Quasi-saic elecromagneism. The field produced by primary charge Q and by an uncharged conducing plane disanced from Q by disance d is equal o he field produced wihou conducing plane by wo following
More informationEXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE
Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,
More information14 Autoregressive Moving Average Models
14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class
More informationGuest Lectures for Dr. MacFarlane s EE3350 Part Deux
Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A
More informationSome Basic Information about M-S-D Systems
Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,
More information( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is
UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires
More informationSecond quantization and gauge invariance.
1 Second quanizaion and gauge invariance. Dan Solomon Rauland-Borg Corporaion Moun Prospec, IL Email: dsolom@uic.edu June, 1. Absrac. I is well known ha he single paricle Dirac equaion is gauge invarian.
More informationSTATE-SPACE MODELLING. A mass balance across the tank gives:
B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing
More informationLinear Response Theory: The connection between QFT and experiments
Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and
More informationLAPLACE TRANSFORM AND TRANSFER FUNCTION
CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions
More informationProblem Set #1. i z. the complex propagation constant. For the characteristic impedance:
Problem Se # Problem : a) Using phasor noaion, calculae he volage and curren waves on a ransmission line by solving he wave equaion Assume ha R, L,, G are all non-zero and independen of frequency From
More informationNavneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi
Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec
More informationLet us start with a two dimensional case. We consider a vector ( x,
Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our
More informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationA DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS
A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS Xinping Guan ;1 Fenglei Li Cailian Chen Insiue of Elecrical Engineering, Yanshan Universiy, Qinhuangdao, 066004, China. Deparmen
More informationAnalyze patterns and relationships. 3. Generate two numerical patterns using AC
envision ah 2.0 5h Grade ah Curriculum Quarer 1 Quarer 2 Quarer 3 Quarer 4 andards: =ajor =upporing =Addiional Firs 30 Day 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 andards: Operaions and Algebraic Thinking
More informationMagnetic Field Synthesis in an Open-Boundary Region Exploiting Thévenin Equivalent Conditions
Magneic Field Synhesis in an Open-Boundary Region Exploiing Thévenin Equivalen Condiions Paolo DI BARBA, PhD Dep of Indusrial and Informaion Engineering Universiy of Pavia, Ialy paolo.dibarba@unipv.i Ouline
More informationt is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...
Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger
More information15. Vector Valued Functions
1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,
More informationEchocardiography Project and Finite Fourier Series
Echocardiography Projec and Finie Fourier Series 1 U M An echocardiagram is a plo of how a porion of he hear moves as he funcion of ime over he one or more hearbea cycles If he hearbea repeas iself every
More informationThe motions of the celt on a horizontal plane with viscous friction
The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More information23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes
Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion
More informationFinal Spring 2007
.615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o
More informationLecture Outline. Introduction Transmission Line Equations Transmission Line Wave Equations 8/10/2018. EE 4347 Applied Electromagnetics.
8/10/018 Course Insrucor Dr. Raymond C. Rumpf Office: A 337 Phone: (915) 747 6958 E Mail: rcrumpf@uep.edu EE 4347 Applied Elecromagneics Topic 4a Transmission Line Equaions Transmission These Line noes
More informationClass Meeting # 10: Introduction to the Wave Equation
MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion
More informationGENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT
Inerna J Mah & Mah Sci Vol 4, No 7 000) 48 49 S0670000970 Hindawi Publishing Corp GENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT RUMEN L MISHKOV Received
More informationTwo Coupled Oscillators / Normal Modes
Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own
More informationOrdinary Differential Equations
Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a
More informationLecture Notes 2. The Hilbert Space Approach to Time Series
Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship
More informationR.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder#
.#W.#Erickson# Deparmen#of#Elecrical,#Compuer,#and#Energy#Engineering# Universiy#of#Colorado,#Boulder# Chaper 2 Principles of Seady-Sae Converer Analysis 2.1. Inroducion 2.2. Inducor vol-second balance,
More informationThen. 1 The eigenvalues of A are inside R = n i=1 R i. 2 Union of any k circles not intersecting the other (n k)
Ger sgorin Circle Chaper 9 Approimaing Eigenvalues Per-Olof Persson persson@berkeley.edu Deparmen of Mahemaics Universiy of California, Berkeley Mah 128B Numerical Analysis (Ger sgorin Circle) Le A be
More informationRC, RL and RLC circuits
Name Dae Time o Complee h m Parner Course/ Secion / Grade RC, RL and RLC circuis Inroducion In his experimen we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors.
More informationChapter 2. First Order Scalar Equations
Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.
More informationCHAPTER 2 Signals And Spectra
CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par
More informationSolutions from Chapter 9.1 and 9.2
Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is
More informationBasic Circuit Elements Professor J R Lucas November 2001
Basic Circui Elemens - J ucas An elecrical circui is an inerconnecion of circui elemens. These circui elemens can be caegorised ino wo ypes, namely acive and passive elemens. Some Definiions/explanaions
More informationTurbulent Flows. Computational Modelling of Turbulent Flows. Overview. Turbulent Eddies and Scales
School of Mechanical Aerospace and Civil Engineering Turbulen Flows As noed above, using he mehods described in earlier lecures, he Navier-Sokes equaions can be discreized and solved numerically on complex
More information11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu
ON EQUATIONS WITH SETS AS UNKNOWNS BY PAUL ERDŐS AND S. ULAM DEPARTMENT OF MATHEMATICS, UNIVERSITY OF COLORADO, BOULDER Communicaed May 27, 1968 We shall presen here a number of resuls in se heory concerning
More informationFinish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!
MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his
More informationMath 315: Linear Algebra Solutions to Assignment 6
Mah 35: Linear Algebra s o Assignmen 6 # Which of he following ses of vecors are bases for R 2? {2,, 3, }, {4,, 7, 8}, {,,, 3}, {3, 9, 4, 2}. Explain your answer. To generae he whole R 2, wo linearly independen
More informationChapter 3 Boundary Value Problem
Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le
More informationINDEX. Transient analysis 1 Initial Conditions 1
INDEX Secion Page Transien analysis 1 Iniial Condiions 1 Please inform me of your opinion of he relaive emphasis of he review maerial by simply making commens on his page and sending i o me a: Frank Mera
More informationEECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits
EEE25 ircui Analysis I Se 4: apaciors, Inducors, and Firs-Order inear ircuis Shahriar Mirabbasi Deparmen of Elecrical and ompuer Engineering Universiy of Briish olumbia shahriar@ece.ubc.ca Overview Passive
More informationThe Maxwell Equations, the Lorentz Field and the Electromagnetic Nanofield with Regard to the Question of Relativity
The Maxwell Equaions, he Lorenz Field and he Elecromagneic Nanofield wih Regard o he Quesion of Relaiviy Daniele Sasso * Absrac We discuss he Elecromagneic Theory in some main respecs and specifically
More information- If one knows that a magnetic field has a symmetry, one may calculate the magnitude of B by use of Ampere s law: The integral of scalar product
11.1 APPCATON OF AMPEE S AW N SYMMETC MAGNETC FEDS - f one knows ha a magneic field has a symmery, one may calculae he magniude of by use of Ampere s law: The inegral of scalar produc Closed _ pah * d
More informationThe Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales
Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions
More informationAnalysis of Microstrip Coupling Gap to Estimate Polymer Permittivity
Analysis of Microsrip Couplin Gap o Esimae Polymer Permiiviy Chanchal Yadav Deparmen of Physics & Elecronics Rajdhani Collee, Universiy of Delhi Delhi, India Absrac A ap in he microsrip line can be modeled
More informationSub Module 2.6. Measurement of transient temperature
Mechanical Measuremens Prof. S.P.Venkaeshan Sub Module 2.6 Measuremen of ransien emperaure Many processes of engineering relevance involve variaions wih respec o ime. The sysem properies like emperaure,
More informationHW6: MRI Imaging Pulse Sequences (7 Problems for 100 pts)
HW6: MRI Imaging Pulse Sequences (7 Problems for 100 ps) GOAL The overall goal of HW6 is o beer undersand pulse sequences for MRI image reconsrucion. OBJECTIVES 1) Design a spin echo pulse sequence o image
More informationSymmetry and Numerical Solutions for Systems of Non-linear Reaction Diffusion Equations
Symmery and Numerical Soluions for Sysems of Non-linear Reacion Diffusion Equaions Sanjeev Kumar* and Ravendra Singh Deparmen of Mahemaics, (Dr. B. R. Ambedkar niversiy, Agra), I. B. S. Khandari, Agra-8
More informationWaves are naturally found in plasmas and have to be dealt with. This includes instabilities, fluctuations, waveinduced
Lecure 1 Inroducion Why is i imporan o sudy waves in plasma? Waves are naurally found in plasmas and have o be deal wih. This includes insabiliies, flucuaions, waveinduced ranspor... Waves can deliver
More informationRobert Kollmann. 6 September 2017
Appendix: Supplemenary maerial for Tracable Likelihood-Based Esimaion of Non- Linear DSGE Models Economics Leers (available online 6 Sepember 207) hp://dx.doi.org/0.06/j.econle.207.08.027 Rober Kollmann
More informationSection 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients
Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous
More informationPENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD
PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.
More informationVoltage/current relationship Stored Energy. RL / RC circuits Steady State / Transient response Natural / Step response
Review Capaciors/Inducors Volage/curren relaionship Sored Energy s Order Circuis RL / RC circuis Seady Sae / Transien response Naural / Sep response EE4 Summer 5: Lecure 5 Insrucor: Ocavian Florescu Lecure
More information6.003 Homework #9 Solutions
6.00 Homework #9 Soluions Problems. Fourier varieies a. Deermine he Fourier series coefficiens of he following signal, which is periodic in 0. x () 0 0 a 0 5 a k sin πk 5 sin πk 5 πk for k 0 a k 0 πk j
More informationOn Multicomponent System Reliability with Microshocks - Microdamages Type of Components Interaction
On Mulicomponen Sysem Reliabiliy wih Microshocks - Microdamages Type of Componens Ineracion Jerzy K. Filus, and Lidia Z. Filus Absrac Consider a wo componen parallel sysem. The defined new sochasic dependences
More informationTHE MATRIX-TREE THEOREM
THE MATRIX-TREE THEOREM 1 The Marix-Tree Theorem. The Marix-Tree Theorem is a formula for he number of spanning rees of a graph in erms of he deerminan of a cerain marix. We begin wih he necessary graph-heoreical
More informationModal identification of structures from roving input data by means of maximum likelihood estimation of the state space model
Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix
More informationElectrical and current self-induction
Elecrical and curren self-inducion F. F. Mende hp://fmnauka.narod.ru/works.hml mende_fedor@mail.ru Absrac The aricle considers he self-inducance of reacive elemens. Elecrical self-inducion To he laws of
More information6.2 Transforms of Derivatives and Integrals.
SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.
More information2. Nonlinear Conservation Law Equations
. Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear
More informationElectrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit
V() R L C 513 Elecrical Circuis Tools Used in Lab 13 Series Circuis Damped Vibraions: Energy Van der Pol Circui A series circui wih an inducor, resisor, and capacior can be represened by Lq + Rq + 1, a
More informationProgram: RFEM 5, RSTAB 8, RF-DYNAM Pro, DYNAM Pro. Category: Isotropic Linear Elasticity, Dynamics, Member
Verificaion Example Program: RFEM 5, RSTAB 8, RF-DYNAM Pro, DYNAM Pro Caegory: Isoropic Linear Elasiciy, Dynamics, Member Verificaion Example: 0104 Canilever Beam wih Periodic Exciaion 0104 Canilever Beam
More information4.1 Other Interpretations of Ridge Regression
CHAPTER 4 FURTHER RIDGE THEORY 4. Oher Inerpreaions of Ridge Regression In his secion we will presen hree inerpreaions for he use of ridge regression. The firs one is analogous o Hoerl and Kennard reasoning
More informationSystem of Linear Differential Equations
Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y
More informationApplication of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing
Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology
More informationDYNAMIC ECONOMETRIC MODELS vol NICHOLAS COPERNICUS UNIVERSITY - TORUŃ Józef Stawicki and Joanna Górka Nicholas Copernicus University
DYNAMIC ECONOMETRIC MODELS vol.. - NICHOLAS COPERNICUS UNIVERSITY - TORUŃ 996 Józef Sawicki and Joanna Górka Nicholas Copernicus Universiy ARMA represenaion for a sum of auoregressive processes In he ime
More informationEE650R: Reliability Physics of Nanoelectronic Devices Lecture 9:
EE65R: Reliabiliy Physics of anoelecronic Devices Lecure 9: Feaures of Time-Dependen BTI Degradaion Dae: Sep. 9, 6 Classnoe Lufe Siddique Review Animesh Daa 9. Background/Review: BTI is observed when he
More informationMatrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality
Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]
More informationNew effective moduli of isotropic viscoelastic composites. Part I. Theoretical justification
IOP Conference Series: Maerials Science and Engineering PAPE OPEN ACCESS New effecive moduli of isoropic viscoelasic composies. Par I. Theoreical jusificaion To cie his aricle: A A Sveashkov and A A akurov
More informationNon-uniform circular motion *
OpenSax-CNX module: m14020 1 Non-uniform circular moion * Sunil Kumar Singh This work is produced by OpenSax-CNX and licensed under he Creaive Commons Aribuion License 2.0 Wha do we mean by non-uniform
More informationA New Perturbative Approach in Nonlinear Singularity Analysis
Journal of Mahemaics and Saisics 7 (: 49-54, ISSN 549-644 Science Publicaions A New Perurbaive Approach in Nonlinear Singulariy Analysis Ta-Leung Yee Deparmen of Mahemaics and Informaion Technology, The
More informationChapter 2. Models, Censoring, and Likelihood for Failure-Time Data
Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based
More informationLinear Control System EE 711. Design. Lecture 8 Dr. Mostafa Abdel-geliel
Linear Conrol Sysem EE 7 MIMO Sae Space Analysis and Design Lecure 8 Dr. Mosafa Abdel-geliel Course Conens Review Sae Space SS modeling and analysis Sae feed back design Oupu feedback design Observer design
More informationComparative study between two models of a linear oscillating tubular motor
IOSR Journal of Elecrical and Elecronics Engineering (IOSR-JEEE) e-issn: 78-676,p-ISSN: 3-333, Volume 9, Issue Ver. IV (Feb. 4), PP 77-83 Comparaive sudy beween wo models of a linear oscillaing ubular
More information!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)
"#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5
More informationPrinciple and Analysis of a Novel Linear Synchronous Motor with Half-Wave Rectified Self Excitation
Principle and Analysis of a Novel Linear Synchronous Moor wih Half-Wave Recified Self Exciaion Jun OYAMA, Tsuyoshi HIGUCHI, Takashi ABE, Shoaro KUBOTA and Tadashi HIRAYAMA Deparmen of Elecrical and Elecronic
More information2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS
Andrei Tokmakoff, MIT Deparmen of Chemisry, 2/22/2007 2-17 2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS The mahemaical formulaion of he dynamics of a quanum sysem is no unique. So far we have described
More informationComputer-Aided Analysis of Electronic Circuits Course Notes 3
Gheorghe Asachi Technical Universiy of Iasi Faculy of Elecronics, Telecommunicaions and Informaion Technologies Compuer-Aided Analysis of Elecronic Circuis Course Noes 3 Bachelor: Telecommunicaion Technologies
More information0.1 MAXIMUM LIKELIHOOD ESTIMATION EXPLAINED
0.1 MAXIMUM LIKELIHOOD ESTIMATIO EXPLAIED Maximum likelihood esimaion is a bes-fi saisical mehod for he esimaion of he values of he parameers of a sysem, based on a se of observaions of a random variable
More information6.003 Homework #9 Solutions
6.003 Homework #9 Soluions Problems. Fourier varieies a. Deermine he Fourier series coefficiens of he following signal, which is periodic in 0. x () 0 3 0 a 0 5 a k a k 0 πk j3 e 0 e j πk 0 jπk πk e 0
More informationDetermination of the Sampling Period Required for a Fast Dynamic Response of DC-Motors
Deerminaion of he Sampling Period Required for a Fas Dynamic Response of DC-Moors J. A. GA'EB, Deparmen of Elecrical and Compuer Eng, The Hashemie Universiy, P.O.Box 15459, Posal code 13115, Zerka, JORDAN
More informationDifferential Equations
Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding
More informationRobust estimation based on the first- and third-moment restrictions of the power transformation model
h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,
More informationSOLUTIONS TO ECE 3084
SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no
More informationON THE BEAT PHENOMENON IN COUPLED SYSTEMS
8 h ASCE Specialy Conference on Probabilisic Mechanics and Srucural Reliabiliy PMC-38 ON THE BEAT PHENOMENON IN COUPLED SYSTEMS S. K. Yalla, Suden Member ASCE and A. Kareem, M. ASCE NaHaz Modeling Laboraory,
More informationChapter 8 The Complete Response of RL and RC Circuits
Chaper 8 The Complee Response of RL and RC Circuis Seoul Naional Universiy Deparmen of Elecrical and Compuer Engineering Wha is Firs Order Circuis? Circuis ha conain only one inducor or only one capacior
More informationChapter 2: Principles of steady-state converter analysis
Chaper 2 Principles of Seady-Sae Converer Analysis 2.1. Inroducion 2.2. Inducor vol-second balance, capacior charge balance, and he small ripple approximaion 2.3. Boos converer example 2.4. Cuk converer
More informationFrequency-Domain Steady-State Analysis of Circuits with Mem-Elements
ELECO 211 7h Inernaional Conference on Elecrical and Elecronics Engineering, 1-4 December, Bursa, TURKEY Frequency-Domain Seady-Sae Analysis of Circuis wih em-elemens Zdene Kola 1, Dalibor Biole 2, Viera
More informationWavelet Methods for Time Series Analysis. What is a Wavelet? Part I: Introduction to Wavelets and Wavelet Transforms. sines & cosines are big waves
Wavele Mehods for Time Series Analysis Par I: Inroducion o Waveles and Wavele Transforms waveles are analysis ools for ime series and images as a subjec, waveles are relaively new (983 o presen) a synhesis
More informationTraveling Waves. Chapter Introduction
Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from
More informationChapter 4 AC Network Analysis
haper 4 A Nework Analysis Jaesung Jang apaciance Inducance and Inducion Time-Varying Signals Sinusoidal Signals Reference: David K. heng, Field and Wave Elecromagneics. Energy Sorage ircui Elemens Energy
More information