Lazhar Roubache 1, *, Kamel Boughrara 1,Frédéric Dubas 2,andRachidIbtiouen 1
|
|
- Marjory Knight
- 5 years ago
- Views:
Transcription
1 Progress In Electromgnetics Reserch B, Vol. 77, 85 11, 17 Semi-Anlyticl Modeling of Spoke-Type Permnent-Mgnet Mchines Considering the Iron Core Reltive Permebility: Subdomin Technique nd Tylor Polynomil Lzhr Roubche 1, *, Kmel Boughrr 1,Frédéric Dubs,ndRchidIbtiouen 1 Abstrct This rticle presents novel contribution to the improvement of the nlytic modeling of electricl mchines using two-dimensionl -D) subdomin technique with Tylor polynomil. To vlidte this novel method, the semi-nlyticl model hs been implemented for spoke-type permnentmgnet PM) mchines STPMM). Mgnetosttic Mxwell s equtions re solved in polr coordintes nd in ll prts of the mchine. The globl solution is obtined by using the trditionl boundry conditions BCs), in ddition to new rdil BCs e.g., between the PMs nd the rotor teeth) which re trduced into system of liner equtions ccording to Tylor series expnsion. The mgnetic field clcultions re performed for two different vlues of iron core reltive permebility viz., 1 nd 1,) nd compred to finite-element method FEM) predictions. The results show tht very good greement is obtined. 1. INTRODUCTION The full clcultion of mgnetic field in electricl mchines is the first step for their design nd optimiztion. The methods of mgnetic field prediction cn be clssified in vrious ctegories [1]: Lehmnn s grphicl, numericl equivlent circuit, Schwrz-Christoffel mpping, Mxwell-Fourier Some comprehensive reviews of mgnetic field prediction in electricl mchine cn be found in [1 8], nd their references. Numericl methods i.e., the finite-element, finite-difference, or boundry-element nlysis) [9 11], which cn be clssified s the very ccurte method compred to the rel results with lrge flexibility to vrious geometries, include nonliner nd nonhomogeneous mterils. The most ccurte models re three-dimensionl 3-D) numericl methods. However, these pproches re very time-consuming nd not very suitble for optimiztion nlysis. Nevertheless, in [1, 13], it is possible to optimize electromgnetic systems from numericl methods. Nowdys, in order to reduce the computtion time, hybrid numericl methods cn be developed [14]. The ctul design works re minly bsed on semi-)nlyticl models i.e., equivlent circuit, SC mpping nd Mxwell-Fourier methods) [1]. These pproches hve been proposed under some geometricl nd physicl ssumptions. Mxwell-Fourier methods re one of the most resent semi-nlytic pproches with good ccurte results in D or 3-D electromgnetic performnce clcultion. These models re bsed on the forml resolution of Mxwell s equtions pplied in ech subdomin of the electromgnetic devices. These subdomins cn be divided into two types of subdomins, viz., periodic e.g., ir-gp) nd non-periodic e.g., slots, teeth, tooth-tips... ). In the second type, the generl solution is derived by pplying homogeneous Neumnn BCs. These conditions re obtined from the pproximtion tht the iron prts re considered to be infinitely permeble. In the literture, we find the pplictions of these methods for the nlysis of Received 1 My 17, Accepted 8 July 17, Scheduled 7 July 17 * Corresponding uthor: Lzhr Roubche roubche.lzhr@gmil.com). 1 Ecole Ntionle Polytechnique LRE-ENP), Algiers, 1, Av. Psteur, El Hrrch, BP 18, 16, Algeri. Déprtement ENERGIE, FEMTO-ST, CNRS, Univ. Bourgogne Frnche-Comté, Belfort, Frnce.
2 86 Roubche et l. severl types of electricl mchines, such s [, 7, 8, 15 19] for PM synchronous mchines, [ ] for solid or cge rotor induction motors, nd [3 5] for reluctnce mchines. It is interesting to note tht n overview on the existing semi-)nlyticl models in Mxwell-Fourier methods with the locl/globl sturtion hs been relized in [1], whose some detils nd the dis)dvntges of these techniques cn be found. In [6,3 5], the uthors pproximte the djcent regions e.g., rotor slots/teeth)in the hrmonic modeling technique s one homogeneous region with reltive permebility developed s Fourier series expnsion. In [1], Dubs nd Boughrr developed the first model introducing the iron prts in the mgnetic field clcultion using subdomin technique, where the uthors solved prtil differentil equtions EDPs) of mgnetic potentil vector in Crtesin coordintes in which the subdomins connection is performed directly in both directions i.e., x-ndy-edges). However, to pply this contribution in the modultion of electricl mchines, it should trnsform it to polr coordintes. Thus, the work in this pper tkes prt in the development nd improvement of the subdomin technique on the scientific topic. The im of this pper is to propose new contribution improving -D subdomin technique with Tylor polynomil by focusing on the considertion of iron. The proposed technique involves solution in polr coordintes of Lplce s nd Poisson s equtions in the sttor yoke, sttor slots nd teeth, ir-gp, rotor teeth, PMs, nd rotor nonmgnetic regions. For the rotor i.e., PMs/teeth) nd sttor i.e., slots/teeth), we tke the generl solution considering the nonhomogeneous Neumnn BCs. In ddition to trditionl interfce conditions between two djcent regions represented with rdius vlue nd intervl of ngles e.g., between the ir-gp nd sttor slots/teeth), we dd the considertion of the interfce conditions between two djcent regions represented with ngle vlue nd intervl of rdius e.g., between the PMs nd the rotor teeth). The first type of interfce conditions permits to get system of liner equtions ccording to Fourier series expnsion, nd the second type permits to obtin system of equtions ccording to Tylor series expnsion. All results from the developed semi-nlyticl model re then compred to those found by FEM [6].. PROBLEM DEFINITIONS..1. Mgnetic Field Solution Figure 1 represents the topology of frctionl-slot STPM mchines with the definition of region where Region I represents the ir-gp, Region II the nonmgnetic mteril under rotor mgnetic i.e., under PMs nd rotor teeth), Region III the sttor yoke, Region IV the buried PMs, Region V the rotor Figure 1. Exmple of frctionl-slot STPM mchines i.e., 4-slots/-poles) with the definition of regions.
3 Progress In Electromgnetics Reserch B, Vol. 77, teeth, Region VI the sttor slots, nd Region VII the sttor teeth. The -D semi-nlyticl model is formulted in mgnetic vector potentil A z nd polr coordintes r, θ) with the following ssumptions: the end-effects re neglected i.e., tht the mgnetic vribles re independent of z); the sttor tooth-tips nd rotor bridge re not considered. However, they cn be introduced esily; the sttor slots/teeth, rotor teeth nd PMs hve rdil sides; the current density hs only one component long the z-xis; the electricl conductivities of mterils re ssumed to be null i.e., the eddy-currents induced in the copper/iron/pms re neglected); the PMs demgnetiztion curve is ssumed to be liner; the direction of PMs mgnetiztion is supposed purely tngentil, i.e., M = {,M θ, } In this rticle, we tke into ccount the iron core reltive permebility in sttor i.e., yoke nd teeth) nd rotor i.e., teeth) subdomin regions. The generl EDP issued from mgnetosttic Mxwell s equtions in continuous nd isotopic region cn be expressed, in term of mgnetic potentil vector A z,by ΔA z r, θ) =, in Region I, II, III, V, nd VII 1) ΔA z r, θ) = μ M, in Region IV ) ΔA z r, θ) = μ J z, in Region VI 3) where M is the PMs mgnetiztion, J z the current density in the sttor slots, nd μ the vcuum permebility. The field vectors B = {B r ; B θ ;} nd H = {H r ; H θ ;} in the different regions re coupled by: B = μ H, in Region I, II, nd VI 4) B = μ μ rm H + μ M, in Region IV 5) B = μ μ rch, in Region III, V, nd VII 6) where μ rm is the reltive recoil permebility of PMs, nd μ rc is the reltive permebility of the iron core. Using B = A,ther- ndθ-components of mgnetic flux density re deduced from A z by B r = 1 r A θ nd B θ = A r.1. Generl Solution of Lplce s Eqution with Nonhomogeneous Neumnn BCs In generl cse of electricl mchine nlysis using subdomin technique, the iron prts re considered to hve infinite reltive permebility which leds to homogeneous Neumnn BCs in slots/pms/teeth subdomins), nd Lplce s, Poisson s or Helmholtz s equtions should be solved only in slots/pms/airgp regions [8, ]. In order to dd rotor/sttor teeth in the mgnetic field prediction in electricl mchines, it is necessry to consider the generl solution of Mxwell s equtions with nonhomogeneous Neumnn BCs. In slots/pms/teeth subdomins, we hve to solve the Lplce s eqution in polr coordintes A z r, θ) r + 1 A z r, θ) + 1 A z r, θ) r r r θ = 8) The electricl mchines hve cylindricl form, nd the following periodicity condition between nd π should be dded to solution A z r, θ) θ= = A z r, θ) θ=π 9) Considering the periodicity condition given in Eq. 9), the generl solution of Eq. 8) cn be written s A z r, θ)=d1 +D lnr)+ D1 k r k +D k r k) sink θ θ 1 ))+ D3 k r k +D4 k r k) cosk θ θ 1 )) 1) 7)
4 88 Roubche et l. A Z D ) r θ = Δ A Z = A Z D ) 1 r θ = θ 1 θ R1 R Figure. Slots/PMs/teeth subdomins with nonhomogeneous Neumnn BCs. Figure 3. Description: Air-gp i.e., Region I), nonmgnetic mteril under rotor mgnetic i.e., Region II) nd sttor yoke Region III). where k is positive integer, nd D1 D4 k re the integrtion constnts. As shown in Figure, the slots/pms/teeth subdomins re delimited by r [R 1 ; R ] θ [θ 1 ; θ ] nd re ssocited with the following BCs: Ar, θ) D 1 r) = θ = k D1 k r k + D k r k) 11) θ=θ1 D r) = Ar, θ) θ = θ=θ k D1 k r k + D k r k) cos k)+ k D3 k r k +D4 k r k) sink) 11b) In the cse where D 1 r) = D r) = i.e., homogeneous Neumnn BCs), the derivtion of Eqs. 11) nd 11b) gives k D1 k r k + D k r k) = D1 k = D k = 1) k D3 k r k + D4 k r k) sin k)= k = mπ 1b) where m is positive integer, nd the solution of Eq. 1) cn be simplified s ) A z r, θ) =D1 + D lnr)+ D3 m r mπ + D4m r mπ mπ ) cos θ θ 1) 13) m=1 As explicted in [1], D1 D4 m cn be determined by Fourier series expnsions ssocited with BCs for the rdii R 1 nd R. In the generl cse where the vrition of A z versus θ for θ = θ 1 nd θ = θ is not necessrily null, it should consider the generl solution of Eq. 1), which cn be written with nother form s ) A z = D1 + D lnr)+ B1 m r mπ + B m r mπ mπ ) cos θ θ 1) + m=1 D1 k r k + D k r k) sin kθ θ 1 )) + k mπ D3 k r k + D4 k r k) cos kθ θ 1 )) 14)
5 Progress In Electromgnetics Reserch B, Vol. 77, where B1 m & B m re new integrtion constnts which respectively represent D3 k & D4 k for k = mπ/. It is interesting to note tht the dded Fourier constnts D1 k D4 k cn be determined by Tylor series expnsions ssocited with BCs for θ 1 nd θ, nd some explntions of the obtined equtions re explined in Section.6. The solution of Poisson s equtions is derived by dding the corresponding prticulr solution... Solution of Lplce s Eqution in Regions I, II nd III 1) Air-gp Region I): In the ir-gp, which is n nnulr domin t r [R ; R 3 ] see Figure 3), the solution of Eq. 1) is defined by A zi r, θ) =A1 + A lnr)+ A1n r n + A n r n) sinnθ)+ A3n r n + A4 n r n) cosnθ) 15) n=1 where n is positive integer, nd A1 A4 n re the integrtion constnts in Region I. ) Nonmgnetic Mteril under Rotor Mgnetic Region II): The generl solution in this subdomin hs the sme form s the solution in Region I. Nevertheless, in dding the condition of the finite vlue of A z for r =, the solution cn be written s A zii r, θ) = A5 n r n sinnθ)+ A6 n r n cosnθ) 16) n=1 where A5 n & A6 n re the integrtion constnts in Region II. 3) Sttor Yoke Region III): In the sttor yoke, which is n nnulr domin t r [R 4 ; R ext ]see Figure 3), the generl solution hs the sme form s the solution in Region I. Nevertheless, in dding the Dirichlet s BC for r = R ext, the solution cn be written s ) n r ) n ) r A ziii r, θ) = A7 n sinnθ)+ R n=1 ext R ext n=1 where A7 A8 n re the integrtion constnts in Region III. n=1 n=1 A8 n r R ext ) n ) ) r n cosnθ) 17) R ext.3. Solution of Poisson s Eqution with Nonhomogeneous Neumnn in Regions IV nd VI 1) Buried PMs Region IV): In ech buried PM subdomin j) of Region IV see Figure 1), we hve to solve Eq. ), i.e., Poisson s eqution. Becuse the direction of PMs Mgnetiztion is purely tngentil, Eq. ) cn be reduced to A z r, θ) r + 1 r A z r, θ) r + 1 r A z r, θ) θ = μ M θ r where M θ = M j = 1) j B rm /μ with j vrying from 1 to p poles, nd B rm is the remnent flux density of PMs. As shown in Figure 4), the jth PM is ssocited with nonhomogeneous Neumnn BCs, nd using the method explicted in Section., the solution of Eq. 18) cn be written s A ziv j r, θ) = B1 j + B j lnr) μ M j r k mπ D1 jk r k + D jk r k) sin m=1 k B1 jm r mπ θ α j + 18) ) + B jm r mπ mπ cos θ α j + )) )) D3 jk r k + D4 jk r k) cos k θ α j + )) 19)
6 9 Roubche et l. ) b) Figure 4. Nonhomogeneous Neumnn BCs for ) jth buried PMs nd b) ith sttor slot. where m nd k re the positive integers; B1 j B jm nd D1 jk D4 jk re the integrtion constnts; α j is the ngulr position of the jth PMs; is the PM-opening. ) Sttor Slots Region VI):In ech slot subdomin i) of Region VI, we hve to solve Eq. ), i.e., Poisson s eqution, A z r, θ) r + 1 A z r, θ) + 1 A z r, θ) r r r θ = μ J i ) where J zi is the current density in the ith sttor slot with i vrying from 1 to Q s in which Q s represents the number of sttor slots. Tking into ccount the BCs shown in Figure 4b), the solution of Eq. ) cn be written s A zv Ii r, θ) =C1 i + C i lnr) 1 4 μ J i r + + E1 ik r k + E ik r k) sin k θ γ i + c m=1 )) + C1 im r mπ c k mπ c ) + C im r mπ mπ c cos θ γ i + c )) c E3 ik r k + E4 ik r k) cos k θ γ i + c )) 1) where m nd k re the positive integers; C1 i C im nd E1 ik E4 ik re the integrtion constnts; γ i is the ngulr position of the ith sttor slot; c is the sttor slot-opening Solution of Lplce s Eqution with Nonhomogeneous Neumnn BCs in Regions V nd VII 1) Rotor Teeth Region V): The jth rotor tooth is subdomin region with nonhomogeneous Neumnn BCs s shown in Figure 5). As explicted in Section., the solution of Eq. 1), i.e., Lplce s eqution, cn be written s ) A zv j r, θ) = B3 j + B4 j lnr)+ B3 jl r lπ b + B4 jl r lπ lπ b cos θ β j + b )) b l=1 + D5 jk r k + D6 jk r k) sin k θ β j + b )) + k lπ b D7 jk r k + D8 jk r k) cos k θ β j + b )) )
7 Progress In Electromgnetics Reserch B, Vol. 77, ) b) Figure 5. Nonhomogeneous Neumnn BCs for ) jth rotor nd b) ith sttor teeth. where l nd k re the positive integers; B3 j B4 jl nd D5 jk D8 jk re the integrtion constnts; β j is the ngulr position of the j th rotor teeth; b is the rotor tooth-opening. ) Sttor Teeth Region VII): The ith sttor tooth is subdomin region with nonhomogeneous Neumnn BCs shown in Figure 5b). As explined in Section., the solution of Eq. 1), i.e., Lplce s eqution, cn be written s A zv IIi r, θ) = C3 i + C4 i lnr)+ + + l=1 C3 il r lπ d + C4 il r lπ d ) lπ cos d E5 ik r k + E6 ik r k) sin k θ δ i + d )) E7 ik r k + E8 ik r k) cos k θ δ i + d )) k lπ d θ δ i + d )) where l nd k re the positive integers; C3 i C4 il nd E5 ik E8 ik re the integrtion constnts; δ i is the ngulr position of the ith sttor teeth; d is the sttor tooth-opening..3.. Interfces Conditions between Regions To determine the integrtion constnts in Eqs. 15) 17), 19), nd 1) 3), the BCs t the interfce between the vrious regions should be introduced. The interfces conditions in this model cn be divided into two types. One is over ngle intervl for given rdius vlue {R 1,R,R 3,R 4 }, nd the other is over rdius intervl for given ngle {α j ± /, β j ± b/, γ i ± c/, δ i ± d/}. As the obtined equtions from ech type of these interfces conditions hve the sme form, it is sufficient to show one exmple for ech type. The interfce conditions re: between Region IV, V nd I t r = R : 3) A zi R,θ)=A ziv j R,θ) for α j θ α j + A zi R,θ)=A zv j R,θ) for β j b θ β j + b H θi R,θ)=H θiv j R,θ) for α j θ α j + 4) 5) 6)
8 9 Roubche et l. H θi R,θ)=H θv j R,θ) for β j b θ β j + b 7) between Region IV nd V t α j + / =β j b/ ndα j+1 / =β j + b/ forr [R 1 ; R ]: A ziv j r, α j + ) = A zv j r, β j b ) 8) H rivj r, α j + ) = H rvj r, β j b ) 9) A ziv j+1) r, α j+1 ) = A zv j r, β j + b ) 3) H rivj+1) r, α j+1 ) = H rvj r, β j + b ) 31) The interfce conditions in Eqs. 4) 7) concern regions with different subdomin frequencies which need Fourier series expnsions to stisfy equlities of mgnetic potentil vector nd tngentil mgnetic field. According to Fourier series expnsion, the interfce condition in Eq. 4) gives = B1 j + B j lnr ) μ M j R k mπ D3 jk R k + D4 jkr k B1 jm R mπ +B jm R mπ + + k mπ α j + α j D3 jk R k + D4 jk R k ) α j+ α j D1 jk R k + D jkr k ) α j+ α j cos k θ α j + )) dθ = 1 D1 jk R k +D jk R k ) α j+ α j ) α j+ α j sin k θ α j + )) dθ α j + α j sin k θ α j + )) cos A I R,θ)dθ 3) mπ cos k θ α j + )) mπ cos θ α j + )) dθ θ α j + )) dθ mπ A I R,θ)cos θ α j + )) dθ 33) The interfce condition in Eq. 5) gives B3 j + B4 j lnr ) b k mπ D7 jk R k + D8 jkr k B3 jm R mπ +B4 jm R mπ + b D5 jk R k + D6 jkr k ) β j+ b β j b ) β j+ b β j b cos k θ β j + b )) dθ = 1 b D5 jk R k +D6 jk R k ) β j+ b β j + b sin k θ β j + b )) dθ β j + b β j b A I R,θ)dθ 34) sin k θ β j + b )) mπ cos θ β j + b )) dθ
9 Progress In Electromgnetics Reserch B, Vol. 77, = b + b k mπ b β j + b β j + b D7 jk R k + D8 jkr k ) β j+ b β j + b cos k θ β j + b )) mπ cos θ β j + b )) dθ mπ A I R,θ)cos θ β j + b )) dθ 35) The interfces conditions in Eqs. 6) nd 7) give A R = 1 π n A1n R n 1 A n R n 1 ) 1 = μ π Q r α j + j=1 α j Q r α j + j=1 α j Q r n A3n R n 1 A4 n R n 1 ) 1 = μ π α j + j=1 α j H θiv j R,θ)dθ + 1 π Q r H θiv j R,θ)sinnθ)dθ+ 1 π β j + b j=1 β j b H θiv j R,θ)cosnθ)dθ+ 1 π Q r j=1 β j b Q r H θv j R,θ)dθ 36) β j + b β j + b j=1 β j b H θv j R,θ)sinnθ)dθ 37) H θv j R,θ)cosnθ)dθ 38) The interfce conditions in Eqs. 8) 31) concern two polynomils equtions with different degree, which need Tylor series expnsions round R t such s R t ]R 1 ; R ] to stisfy equlities of mgnetic potentil vector nd rdil mgnetic field. The Tylor polynomil degree is chosen to stisfy the lck of equtions obtined by Fourier series expnsion of previous interfce conditions compred to the number of unknowns in system, so for k vrying from 1 to K r it is enough to tke Tylor polynomil of degree K r 1. Therefore, the interfce condition in Eq. 8) gives A IV j R t,α j + ) = A Vj R t,β j b k A IV j R r k t,α j + ) = k A Vj r k R t,β j b ) From the interfce condition in Eq. 9), we get Hr IV j R t,α j + ) = Hr Vj R t,β j b ) k Hr IV j R r k t,α j + ) = k Hr Vj r k R t,β j b ) From the interfce condition in Eq. 3), we get A IV j+1 R t,α j+1 ) k A IV j+1 R r k t,α j + ) From the interfce condition in Eq. 31), we get Hr IV j+1 R t,α j+1 ) k Hr IV j+1 R r k t,α j + ) ) = A Vj R t,β j + b ) = k A Vj r k R t,β j b ) = Hr Vj R t,β j + b ) = k Hr Vj r k R t,β j b ) 39) 4) 41) 4) 43) 44) 45) 46)
10 94 Roubche et l. There re no simple expressions to define the system obtined by Tylor series expnsions 39) 46). However, the coefficients of the unknowns in this system cn be defined recursively. For exmple, the development of interfce conditions Eq. 39)) nd Eq. 4)) permits to define new coefficients C B1k for the vector potentil of region IV s k = k =1 k = B1 j B1 jmr mπ t cosmπ)+...= A zv j R t,β j b ) Bj mπ B1 jm R mπ t cosmπ)+...= AzV j R t,β j b ) R t R t r Bj mπ mπ B1 1 Rt jm R mπ t cosmπ)+...= A zv j R t R t r R t,β j b ) C B1 = R mπ t cosmπ) C B11 = mπ C B1 R t C B1 = The coefficients C B1 k cn be defined by the recursive formul s: { C B1 = R mπ t cosmπ) C B1 k+1 = mπ +1 k 48) R t C B1 k The other coefficients re deduced with the sme resoning. It is interesting to note tht the interfce conditions for the rdii {R 1,R 3,R 4 } nd for the ngels {γ i ± c/, δ i ± d/} re developed with the sme previous method 3. RESULTS AND VALIDATION The developed semi-nlyticl method for frctionl-slot STPM mchines tking into ccount the iron core reltive permebility is used to determine electromgnetic performnces viz., the mgnetic flux density, the mgnetic flux linkge, the bck EMF, the electromgnetic/cogging/ripple torques, the nonintrinsic unblnced mgnetic forces,... ) whose vrious formuls hve been clrified in [19]. The min dimensions nd prmeters of the frctionl-slot STPM mchines e.g., 4-slots/-poles) with the buried PMs nd single lyer winding re given in Tble 1. Then, semi-nlytic results re verified by -D FEM [6]. For the finite-elements simultion, we hve used 99,199 nodes nd 194,796 elements. Tble 1. Prmeters of the studied mchine. mπ 1 R t Symbol Prmeters Vlue Units) B rm Remnence flux density of PMs 1. T) μ rm Reltive permebility of PMs 1 N c Number of conductors per sttor slot 18 I m Pek phse current 18 A) Q s Number of sttor slots 4 c Sttor slot-opening 7.5 deg) PM-opening 3.7 deg) p Number of pole pirs 11 R ext Rdius of the externl sttor surfce 11 mm) R 4 Outer rdius of sttor slot 1 mm) R 3 Rdius of the sttor outer surfce 79 mm) R Rdius of the rotor inner surfce t the PM surfce 78 mm) R 1 Rdius of the rotor inner surfce t the PM bottom 57 mm) g Air-gp length 1 mm) L u Axil length of the mchine 63 mm) Ω Mechnicl pulse of synchronism 1,1 rpm) C B11 47)
11 Progress In Electromgnetics Reserch B, Vol. 77, FEM Anlytic μ rc = 1. FEM Anlytic μ rc = 1 B r T) B θ T) Angle mech. degrees) 1 ) Angle mech. degrees)..5.1 B r T) B θ T) FEM Anlytic μrc = 1 -. FEM Anlytic μrc = Angle mech. degrees) b) Angle mech. degrees).4 1 FEM Anlytic μrc = 1. FEM Anlytic μrc = 1 B r T) B θ T) Angle mech. degrees) c) Angle mech. degrees) Figure 6. Rdil nd tngentil components of the mgnetic flux density in the middle of the ir-gp i.e., Region I): ) PMs lone, b) sttor current lone nd c) on lod condition.
12 96 Roubche et l. The wveforms of rdil nd tngentil components of the mgnetic flux density in the vrious regions re computed with finite number of hrmonic terms, viz., N =, M = L =ndk =. The nlytic clcultion of mgnetic flux distribution in ll regions is done with two different vlues of the iron core reltive permebility viz., 1 nd 1,). In ech cse, the reltive permebility is supposed the sme for ll iron prts i.e., sttor yoke, sttor nd rotor teeth). However, it should be noted tht it is possible to tke specil permebility vlue for ech region. In Figure 6, comprison between the numericl results nd semi-nlyticl predictions is shown in term of the mgnetic flux density in the middle of the ir-gp i.e., Region I) for three operting conditions : i) PMs lone, ii) sttor currents lone, nd iii) on lod condition. One cn see tht the diminution of the iron core reltive permebility is ccompnied by decrese in mgnetic flux mplitude. It cn be seen tht very good greement for the rdil nd tngentil components of the mgnetic flux density is obtined. The proposed model is chrcterized by the bility to know the mgnetic flux density in ll mchine regions, such s sttor/rotor yoke, sttor slots, sttor/rotor teeth nd PMs. Figure 7) shows the rdil nd tngentil components of the mgnetic flux density in the middle of the PMs/rotor teeth under pole-pitch for on lod condition. Very good greement is obtined for the tngentil component. For B r T) Mgnet Anlytic Rotor tooth Femm μ rc = Angle mech. degrees).3 B θ T) Mgnet Anlytic Rotor tooth Femm μ rc = Angle mech. degrees) ).3 Anlytic Femm μ rc = 1 Anlytic Femm μ rc = 1. Sttor tooth B r T).1 Sttor Slot B θ T) Sttor Slot Sttor tooth Angle mech. degrees) Angle mech. degrees) b) Figure 7. Rdil nd tngentil components of the mgnetic flux density under tooth-pitch for on lod condition in the middle of ) PMs/rotor tooth nd b) sttor slot/tooth.
13 Progress In Electromgnetics Reserch B, Vol. 77, Shft Mgnet Sttor tooth Yoke Shft Mgnet Sttor tooth Yoke B r T) Anlytic Femm μ rc = Rdius mm) Shft Rotor tooth Sttor slot Yoke B θ T) ) Anlytic Femm μ rc = Rdius mm).6.4. Shft Rotor tooth Sttor slot Yoke B r T). B θ T) Anlytic Femm μ rc = Rdius mm) b) Anlytic Femm μ rc = Rdius mm) Figure 8. Rdil nd tngentil components of the mgnetic flux density for on lod condition t: ) θ = α 1 nd b) θ = γ 3. the rdil component, the results re good with greeble error, the sme order of ccurcy is seen for sttor slot/teeth under tooth-pitch shown in Figure 7b). Another presenttion of obtining results is shown in Figure 8, where the rdil nd tngentil components of the mgnetic flux density for on lod condition re plotted in function of the rdius t given ngle. For the semi-nlytic results, the mgnetic flux density is clculted for ech intervl of rdius by the corresponding solution. The results re obtined for θ = α 1 nd θ = γ 3, nd very good greement is obtined. The sttic torque versus rotor position is presented in Figure 9). The sttor currents re supposed constnt, nd the rotor position is updted to hve nominl rottion speed. As obtined, the mechnicl ngle hs Θ rs =π/p fundmentl period, nd the mximum torque is obtined for Θ rs = π/p. The mximum torque decreses for smll vlue of reltive permebility. It cn be seen tht very good greement is obtined for the two proposed vlues of reltive permebility. Figure 9b) shows the electromgnetic torque wveform versus electricl ngle. At ech instnt, the sttor currents re updted to hve sinusoidl current wveform. The initil rotor position is tken for giving mximum torque vlue, i.e., Θ rs = π/p. The induced bck EMF for no-lod condition is shown in Figure 1. The simultion is done for
14 98 Roubche et l. Torque Nm) Anlytic Femm μ rc = Angle mech. degrees) ) Torque Nm) Anlytic Femm μ rc = Electricl ngle degrees) b) Figure 9. ) Sttic torque nd b) electromgnetic torque t Θ rs = π/p for on lod condition with I =18A). 3 Bck EMF V) Anlytic Femm μ rc = Electricl ngle degrees) Figure 1. Induced bck EMF for no-lod condition with I =A). two different vlues of iron core reltive permebility, nd the obtined results confirm the ccurcy of the proposed semi-nlyticl model. Figure 11) shows the evolution of simultion time versus the number of hrmonics. To reduce the computtion time of the semi-nlyticl model, the number of sptil hrmonics N, M, L nd K cn be reduced. However, this leds to reduction in ccurcy. It is seen tht for giving vlue of sptil hrmonics, the RMS error will not be ffected very much by incresing the number of sptil hrmonics. This vlue cn be chosen s n optiml vlue. Figure 11b) shows the evolution of RMS error for the rdil component of the mgnetic flux density in the middle of the ir-gp versus hrmonic order of Regions I, II, nd III i.e., N) for severl vlues of hrmonic order of the other regions. The RMS error is clculted by error = Npc m=1 B Numeric Ir,m / ) B Anlytic Ir,m N pc 49) It cn be seen from Figure 11b) tht the error strts to converge round N = 18, M = L =6nd K =1.
15 Progress In Electromgnetics Reserch B, Vol. 77, Simultion time s) M=6 L=6 K=3 M=6 L=6 K=1 M=4 L=4 K=3 RMS error T) M=6 L=6 K=3 M=6 L=6 K=1 M=4 L=4 K= Number of hrmonics ) Number of hrmonics b) Figure 11. Influence of the number of hrmonics ssocited with Regions I, II nd III i.e., N) onthe rdil component of the mgnetic flux density in the middle of the ir-gp versus numeric simultions for severl vlues of the number of hrmonics ssocited with the other regions: ) the clcultion time nd b) the RMS error. 4. CONCLUSION In this pper, we hve proposed n improved semi-nlyticl model bsed on subdomin technique with Tylor polynomil for prediction open-circuit, rmture rection nd on lod mgnetic field distribution in STPM mchines. The proposed model tkes into ccount the reltive permebility in ll mchine regions including the iron prts. One of the lrgest dvntges of the discussed pproch is its bility to ccount the finite permebility in rotor/sttor teeth nd yoke, with the possibility to get the solution of mgnetic flux in these regions. Currently, only liner PM mteril proprieties re considered. However, nonliner PM mterils could be ccounted by mens of n itertive lgorithm, nd lso the proposed model offers the possibility to tke especil finite permebility for ech PM or rotor/sttor teeth, sme for rotor nd sttor yoke. Semi-nlytic results re in excellent greement with the ones obtined by -DFEM. From these results, we show the ccurcy of the semi-nlyticl model for severl vlues of iron core reltive permebility. APPENDIX A. NOMENCLATURE H rx,h θx Rdil nd tngentil components of the mgnetic flux intensity in x domin where x cn be I, II, III, IV, V, VI, or VII. A zx Mgnetic potentil vector in x domin. α j, Position nd opening width of j th PMs. β j,b Position nd opening width of j th rotor tooth. γ i,c Position nd opening width of i th sttor slot. δ i,d Position nd opening width of i th sttor tooth. n Hrmonic order for ir-gp, nonmgnetic mteril, nd sttor yoke. k, m, l Hrmonic order for slots/pms/teeth. J z Current density. μ Vcuum permebility. μ rc Reltive permebility of iron core. Reltive permebility of PMs. μ rm
16 1 Roubche et l. REFERENCES 1. Dubs, F. nd K. Boughrr, New scientific contribution on the -D subdomin technique in crtesin coordintes: Tking into ccount of iron prts, Mth. Comput. Appl., Vol., No. 1, 17, 17, DOI: 1.339/mc117.. Dubs, F. nd C. Espnet, Anlyticl solution of the mgnetic field in permnent-mgnet motors tking into ccount slotting effect: no-lod vector potentil nd flux density clcultion, IEEE Trns. on Mgn., Vol. 45, No. 5, 97 19, Devillers, E., J. Le Besneris, T. Lubin, M. Hecquet, nd J. P. Lecointe, A review of subdomin modeling techniques in electricl mchines: Performnces nd pplictions, Proc. ICEM, Lusnne, Switzerlnd, Sep. 4 7, Tiegn, H., Y. Amr, nd G. Brkt, Overview of nlyticl models of permnent mgnet electricl mchines for nlysis nd design purposes, Mthemtics nd Computer in Simultion, Vol. 9, , Curti, M., J. J. H. Pulides, nd E. A. Lomonov, An overview of nlyticl methods for mgnetic field computtion, Proc. EVER, Grimldi Forum, Monco, Mr. 31 Apr., Sprngers, R. L. J., J. J. H. Pulides, B. L. J. Gysen, nd E. A. Lomonov, Mgnetic sturtion in semi-nlyticl hrmonic modeling for electric mchine nlysis, IEEE Trns. on Mgn., Vol. 5, No., Art. ID 8141, Pfister, P.-D., X. Yin, nd Y. Fng, Slotted permnent-mgnet mchines: Generl nlyticl model of mgnetic fields, torque, eddy currents, nd permnent-mgnet power losses including the Diffusion effect, IEEE Trns. on Mgn., Vol. 5, No. 5, Art. ID 81313, Dubs, F. nd A. Rhideh, Two-dimensionl nlyticl permnent-mgnet eddy-current loss clcultions in slotless PMSM equipped with surfce-inset mgnets, IEEE Trns. on Mgn., Vol. 5, No. 3, Art. ID 633, Yilmz, M. nd P. T. Krein, Cpbilities of finite element nlysis nd mgnetic equivlent circuits for electricl mchine nlysis nd design, Proc. PESC, Rhodes, Greece, Jun , Sulimn, E. B., F. Khn, nd T. Kosk, Field-excited flux switching motor design, optimiztion nd nlysis for future hybrid electric vehicle using finite element nlysis, Progress In Electromgnetics Reserch B, Vol. 71, , Konwr, R. S., K. Klit, A. Bnerjee, nd W. K. S. Khoo, Electromgnetic nlysis of bridge configured winding cge induction mchine using finite element method, Progress In Electromgnetics Reserch B, Vol. 48, , Schutte, J. nd J. M. Struss, Optimiztion of trnsverse flux liner PM genertor using 3D finite element nlysis, Proc. ICEM, Rome, Itly, Sep. 6 8, Espnet, C., C. Kieffer, A. Mir, S. Giurge, nd F. Gustin, Optiml design of specil permnent mgnet synchronous mchine for mgnetocloric refrigertion, Proc. ECCE, Denver, CO, USA, Sep , Benlmine, R., F. Dubs, S.-A. Rndi, D. Lhotellier, nd C. Espnet, 3-D numericl hybrid method for PM eddy-current losses clcultion: Appliction to xil-flux PMSMs, IEEE Trns. on Mgn., Vol. 51, No. 7, Art. ID 81611, Rhideh, A., H. Moyed-Jhromi, M. Mrdneh, F. Dubs, nd T. Korkinitis, Anlyticl clcultions of electromgnetic quntities for slotted Brushless mchines with surfce-inset mgnets, Progress In Electromgnetics Reserch B, Vol. 7, 49 65, Teymoori, S., A. Rhideh, H. Moyed-Jhromi, nd M. Mrdneh, -D nlyticl mgnetic field prediction for consequent-pole permnent mgnet synchronous mchines, IEEE Trns. on Mgn., Vol. 5, No. 6, Art. ID 8114, Boughrr, K., T. Lubin, R. Ibtiouen, nd N. Benlll, Anlyticl clcultion of prllel double excittion nd spoke-type permnent-mgnet motors; simplified versus exct model, Progress In Electromgnetics Reserch B, Vol. 47, , Lubin, T., S. Mezni, nd A. Rezzoug, Two-dimensionl nlyticl clcultion of mgnetic field nd electromgnetic torque for surfce-inset permnent-mgnet motors, IEEE Trns. on Mgn.,
17 Progress In Electromgnetics Reserch B, Vol. 77, Vol. 48, No. 6, 8 91, Boughrr, K., R. Ibtiouen, nd F. Dubs, Anlyticl prediction of electromgnetic performnces nd unblnced mgnetic forces in frctionl-slot spoke-type permnent-mgnet mchines, Proc. ICEM, Lusnne, Switzerlnd, Sep. 4 7, 16.. Boughrr, K., N. Tkorbet, R. Ibtiouen, O. Touhmi, nd F. Dubs, Anlyticl nlysis of cge rotor induction motors in helthy, defective, nd broken brs conditions, IEEE Trns. on Mgn., Vol. 51, No., Art. ID 8317, Roubche, L., K. Boughrr, nd R. Ibtiouen, Anlyticl electromgnetic nlysis of multi-phses cge rotor induction motors in helthy, broken brs nd open phses conditions, Progress In Electromgnetics Reserch B, Vol. 7, , 16.. Boughrr, K., F. Dubs, nd R. Ibtiouen, -D nlyticl prediction of eddy currents, circuit model prmeters, nd stedy-stte performnces in solid rotor induction motors, IEEE Trns. on Mgn., Vol. 5, No. 1, Art. ID 7814, Sprngers,R.L.J.,J.J.H.Pulides,B.L.J.Gysen,J.Wrm,ndE.A.Lomonov, Semi nlyticl frmework for synchronous reluctnce motor nlysis including finite soft-mgnetic mteril Permebility, IEEE Trns. on Mgn., Vol. 51, No. 11, Art. ID 81154, Djelloul, K. Z., K. Boughrr, R. Ibtiouen, nd F. Dubs, Nonliner nlyticl clcultion of mgnetic field nd torque of switched reluctnce mchines, Proc. CISTEM, Mrrkech-Benguérir, Mroc, Oct. 6 8, Djelloul, K. Z., K. Boughrr, F. Dubs, nd R. Ibtiouen, Nonliner nlyticl prediction of mgnetic field nd electromgnetic performnces in switched reluctnce mchines, IEEE Trns. on Mgn., 17, DOI: 1.11/TMAG Meeker, D. C., Finite Element Method Mgnetics ver. 4., [Online], Avilble: Apr. 1, 9.
ANALYTICAL CALCULATION OF PARALLEL DOU- BLE EXCITATION AND SPOKE-TYPE PERMANENT- MAGNET MOTORS; SIMPLIFIED VERSUS EXACT MODEL
Progress In Electromgnetics Reserch B, Vol. 47, 145 178, 13 ANALYTICAL CALCULATION OF PARALLEL DOU- BLE EXCITATION AND SPOKE-TYPE PERMANENT- MAGNET MOTORS; SIMPLIFIED VERSUS EXACT MODEL Kmel Boughrr 1,
More informationELE B7 Power Systems Engineering. Power System Components Modeling
Power Systems Engineering Power System Components Modeling Section III : Trnsformer Model Power Trnsformers- CONSTRUCTION Primry windings, connected to the lternting voltge source; Secondry windings, connected
More informationHigh speed machines using advanced magnetic materials analyzed by appropriate finite element models
High speed mchines using dvnced mgnetic mterils nlyzed y pproprite finite element models G. D. KALOKIRIS, A. G. KLADAS Lortory of Electricl Mchines nd Power Electronics, Electric Power Division Fculty
More informationA New Method for Estimating Permanent Magnet Synchronous Machine Parameters
J. Bsic. Appl. Sci. Res., (9)9145-9151, 01 01, TextRod Publiction ISSN 090-4304 Journl o Bsic nd Applied Scientiic Reserch www.textrod.com A New Method or Estimting Permnent Mgnet Synchronous Mchine Prmeters
More informationIRON LOSSES EVALUATION IN HIGH SPEED DRIVES. G. D. Kalokiris 1, P. S. Georgilakis 2 and A. G. Kladas 1. (
IRON LOSSES EVALUATION IN HIGH SPEED DRIVES G. D. Klokiris 1, P. S. Georgilkis 2 nd A. G. Klds 1 1 Ntionl Technicl University of Athens, GR-1578, Greece (e-mil: kldsel@centrl.ntu.gr) 2 Technicl University
More informationFEM ANALYSIS OF ROGOWSKI COILS COUPLED WITH BAR CONDUCTORS
XIX IMEKO orld Congress Fundmentl nd Applied Metrology September 6 11, 2009, Lisbon, Portugl FEM ANALYSIS OF ROGOSKI COILS COUPLED ITH BAR CONDUCTORS Mirko Mrrcci, Bernrdo Tellini, Crmine Zppcost University
More informationStudy on the Calculation of Magnetic Force Based on the Equivalent Magnetic Charge Method
Avilble online t www.sciencedirect.com Physics Procedi 4 () 9 97 Interntionl Conference on Applied Physics nd Industril Engineering Study on the Clcultion of Mgnetic Force Bsed on the Equivlent Mgnetic
More informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
More informationSolutions of Klein - Gordan equations, using Finite Fourier Sine Transform
IOSR Journl of Mthemtics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 13, Issue 6 Ver. IV (Nov. - Dec. 2017), PP 19-24 www.iosrjournls.org Solutions of Klein - Gordn equtions, using Finite Fourier
More informationA ROTATING DISC IN CONSTANT PURE SHEAR BY S. KUMAR AND C. V. JOGA RAO
A ROTATING DISC IN CONSTANT PURE SHEAR BY S. KUMAR AND C. V. JOGA RAO (Deprtment of Aeronuticl Engineering, Indin Institute of Science, Bnglore-3) Received April 25, 1954 SUMMARY The disc of constnt pure
More informationThis final is a three hour open book, open notes exam. Do all four problems.
Physics 55 Fll 27 Finl Exm Solutions This finl is three hour open book, open notes exm. Do ll four problems. [25 pts] 1. A point electric dipole with dipole moment p is locted in vcuum pointing wy from
More information#6A&B Magnetic Field Mapping
#6A& Mgnetic Field Mpping Gol y performing this lb experiment, you will: 1. use mgnetic field mesurement technique bsed on Frdy s Lw (see the previous experiment),. study the mgnetic fields generted by
More informationDESIGN OF FERRITE MAGNET CFSM TO IMPROVE SPEED RESPONSE OF EPS SYSTEM
Interntionl Journl of Automotive Technology, Vol. 18, No. 3, pp. 499 504 (2017) DOI 10.1007/s12239 017 0049 x Copyright 2017 KSAE/ 096 13 pissn 1229 9138/ eissn 1976 3832 DESIGN OF FERRITE MAGNET CFSM
More informationPartial Differential Equations
Prtil Differentil Equtions Notes by Robert Piché, Tmpere University of Technology reen s Functions. reen s Function for One-Dimensionl Eqution The reen s function provides complete solution to boundry
More informationIntroduction to Finite Element Method
Introduction to Finite Element Method Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pn.pl/ tzielins/ Tble of Contents 1 Introduction 1 1.1 Motivtion nd generl concepts.............
More informationMath 100 Review Sheet
Mth 100 Review Sheet Joseph H. Silvermn December 2010 This outline of Mth 100 is summry of the mteril covered in the course. It is designed to be study id, but it is only n outline nd should be used s
More informationModification Adomian Decomposition Method for solving Seventh OrderIntegro-Differential Equations
IOSR Journl of Mthemtics (IOSR-JM) e-issn: 2278-5728, p-issn: 239-765X. Volume, Issue 5 Ver. V (Sep-Oct. 24), PP 72-77 www.iosrjournls.org Modifiction Adomin Decomposition Method for solving Seventh OrderIntegro-Differentil
More informationWELCOME TO THE LECTURE
WELCOME TO THE LECTURE ON DC MOTOR Force on conductor If conductor is plced in mgnetic field nd current is llowed to flow through the conductor, the conductor will experience mechnicl force. N S Electric
More informationPlates on elastic foundation
Pltes on elstic foundtion Circulr elstic plte, xil-symmetric lod, Winkler soil (fter Timoshenko & Woinowsky-Krieger (1959) - Chpter 8) Prepred by Enzo Mrtinelli Drft version ( April 016) Introduction Winkler
More informationLec 3: Power System Components
Lec 3: Power System Components Dr. Mlbik Bsu 8/0/2009 Lesson pln 3 nd L.O. Sequence nlysis exmple ( detil fult nlysis next sem) Trnsformer model recp, tp chnge nd phse chnge, 3-phse Modeling of Synchronous
More informationA Modified ADM for Solving Systems of Linear Fredholm Integral Equations of the Second Kind
Applied Mthemticl Sciences, Vol. 6, 2012, no. 26, 1267-1273 A Modified ADM for Solving Systems of Liner Fredholm Integrl Equtions of the Second Kind A. R. Vhidi nd T. Dmercheli Deprtment of Mthemtics,
More informationPhysics 116C Solution of inhomogeneous ordinary differential equations using Green s functions
Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner
More informationMAC-solutions of the nonexistent solutions of mathematical physics
Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences - Finite Elements - Finite Volumes - Boundry Elements MAC-solutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE
More informationNUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.
NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More informationKINEMATICS OF RIGID BODIES
KINEMTICS OF RIGID ODIES Introduction In rigid body kinemtics, e use the reltionships governing the displcement, velocity nd ccelertion, but must lso ccount for the rottionl motion of the body. Description
More informationADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS
ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS F. Tkeo 1 nd M. Sk 1 Hchinohe Ntionl College of Technology, Hchinohe, Jpn; Tohoku University, Sendi, Jpn Abstrct:
More informationEnergy creation in a moving solenoid? Abstract
Energy cretion in moving solenoid? Nelson R. F. Brg nd Rnieri V. Nery Instituto de Físic, Universidde Federl do Rio de Jneiro, Cix Postl 68528, RJ 21941-972 Brzil Abstrct The electromgnetic energy U em
More informationMultiple Integrals. Review of Single Integrals. Planar Area. Volume of Solid of Revolution
Multiple Integrls eview of Single Integrls eding Trim 7.1 eview Appliction of Integrls: Are 7. eview Appliction of Integrls: Volumes 7.3 eview Appliction of Integrls: Lengths of Curves Assignment web pge
More informationMath& 152 Section Integration by Parts
Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible
More informationPressure Wave Analysis of a Cylindrical Drum
Pressure Wve Anlysis of Cylindricl Drum Chris Clrk, Brin Anderson, Brin Thoms, nd Josh Symonds Deprtment of Mthemtics The University of Rochester, Rochester, NY 4627 (Dted: December, 24 In this pper, hypotheticl
More informationTravelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing
Applied Mthemtics E-Notes 8(8) - c IN 67-5 Avilble free t mirror sites of http://www.mth.nthu.edu.tw/ men/ Trvelling Profile olutions For Nonliner Degenerte Prbolic Eqution And Contour Enhncement In Imge
More informationFrobenius numbers of generalized Fibonacci semigroups
Frobenius numbers of generlized Fiboncci semigroups Gretchen L. Mtthews 1 Deprtment of Mthemticl Sciences, Clemson University, Clemson, SC 29634-0975, USA gmtthe@clemson.edu Received:, Accepted:, Published:
More informationIn the diagram below, the rotation continues until N-S alignment, resulting in lock-up that is, if nothing is done to prevent it.
25-1 DC motors DC motors re importnt in mny pplictions. In portble pplictions using bttery power, DC motors re nturl choice. DC mchines re lso used in pplictions where high strting torque nd ccurte speed
More informationFig. 1. Open-Loop and Closed-Loop Systems with Plant Variations
ME 3600 Control ystems Chrcteristics of Open-Loop nd Closed-Loop ystems Importnt Control ystem Chrcteristics o ensitivity of system response to prmetric vritions cn be reduced o rnsient nd stedy-stte responses
More informationThe Velocity Factor of an Insulated Two-Wire Transmission Line
The Velocity Fctor of n Insulted Two-Wire Trnsmission Line Problem Kirk T. McDonld Joseph Henry Lbortories, Princeton University, Princeton, NJ 08544 Mrch 7, 008 Estimte the velocity fctor F = v/c nd the
More informationPhysics 3323, Fall 2016 Problem Set 7 due Oct 14, 2016
Physics 333, Fll 16 Problem Set 7 due Oct 14, 16 Reding: Griffiths 4.1 through 4.4.1 1. Electric dipole An electric dipole with p = p ẑ is locted t the origin nd is sitting in n otherwise uniform electric
More informationChapter 6 Electrostatic Boundary Value Problems. Dr. Talal Skaik
Chpter 6 Electrosttic Boundry lue Problems Dr. Tll Skik 1 1 Introduction In previous chpters, E ws determined by coulombs lw or Guss lw when chrge distribution is known, or potentil is known throughout
More informationarxiv: v1 [math.gm] 30 Dec 2015
A NEW FRACTIONAL DERIVATIVE WITHOUT SINGULAR KERNEL: APPLICATION TO THE MODELLING OF THE STEADY HEAT FLOW rxiv:161.1623v1 [mth.gm] 3 Dec 215 by Xio-Jun YANG, H. M. SRIVASTAVA b,c, J. A. Tenreiro MACHADO
More informationSummary: Method of Separation of Variables
Physics 246 Electricity nd Mgnetism I, Fll 26, Lecture 22 1 Summry: Method of Seprtion of Vribles 1. Seprtion of Vribles in Crtesin Coordintes 2. Fourier Series Suggested Reding: Griffiths: Chpter 3, Section
More informationMultiple Integrals. Review of Single Integrals. Planar Area. Volume of Solid of Revolution
Multiple Integrls eview of Single Integrls eding Trim 7.1 eview Appliction of Integrls: Are 7. eview Appliction of Integrls: olumes 7.3 eview Appliction of Integrls: Lengths of Curves Assignment web pge
More informationAn approximation to the arithmetic-geometric mean. G.J.O. Jameson, Math. Gazette 98 (2014), 85 95
An pproximtion to the rithmetic-geometric men G.J.O. Jmeson, Mth. Gzette 98 (4), 85 95 Given positive numbers > b, consider the itertion given by =, b = b nd n+ = ( n + b n ), b n+ = ( n b n ) /. At ech
More informationTheoretische Physik 2: Elektrodynamik (Prof. A.-S. Smith) Home assignment 4
WiSe 1 8.1.1 Prof. Dr. A.-S. Smith Dipl.-Phys. Ellen Fischermeier Dipl.-Phys. Mtthis Sb m Lehrstuhl für Theoretische Physik I Deprtment für Physik Friedrich-Alexnder-Universität Erlngen-Nürnberg Theoretische
More informationPhysics 202, Lecture 14
Physics 202, Lecture 14 Tody s Topics Sources of the Mgnetic Field (Ch. 28) Biot-Svrt Lw Ampere s Lw Mgnetism in Mtter Mxwell s Equtions Homework #7: due Tues 3/11 t 11 PM (4th problem optionl) Mgnetic
More information3.4 Numerical integration
3.4. Numericl integrtion 63 3.4 Numericl integrtion In mny economic pplictions it is necessry to compute the definite integrl of relvlued function f with respect to "weight" function w over n intervl [,
More informationNumerical Analysis: Trapezoidal and Simpson s Rule
nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =
More informationTHE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.
THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem
More informationPOLYPHASE CIRCUITS. Introduction:
POLYPHASE CIRCUITS Introduction: Three-phse systems re commonly used in genertion, trnsmission nd distribution of electric power. Power in three-phse system is constnt rther thn pulsting nd three-phse
More informationJackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jckson 2.26 Homework Problem Solution Dr. Christopher S. Bird University of Msschusetts Lowell PROBLEM: The two-dimensionl region, ρ, φ β, is bounded by conducting surfces t φ =, ρ =, nd φ = β held t zero
More informationFredholm Integral Equations of the First Kind Solved by Using the Homotopy Perturbation Method
Int. Journl of Mth. Anlysis, Vol. 5, 211, no. 19, 935-94 Fredholm Integrl Equtions of the First Kind Solved by Using the Homotopy Perturbtion Method Seyyed Mhmood Mirzei Deprtment of Mthemtics, Fculty
More informationINTRODUCTION. The three general approaches to the solution of kinetics problems are:
INTRODUCTION According to Newton s lw, prticle will ccelerte when it is subjected to unblnced forces. Kinetics is the study of the reltions between unblnced forces nd the resulting chnges in motion. The
More informationAP Calculus Multiple Choice: BC Edition Solutions
AP Clculus Multiple Choice: BC Edition Solutions J. Slon Mrch 8, 04 ) 0 dx ( x) is A) B) C) D) E) Divergent This function inside the integrl hs verticl symptotes t x =, nd the integrl bounds contin this
More informationPhys 6321 Final Exam - Solutions May 3, 2013
Phys 6321 Finl Exm - Solutions My 3, 2013 You my NOT use ny book or notes other thn tht supplied with this test. You will hve 3 hours to finish. DO YOUR OWN WORK. Express your nswers clerly nd concisely
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More informationA REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007
A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus
More informationOn the Decomposition Method for System of Linear Fredholm Integral Equations of the Second Kind
Applied Mthemticl Sciences, Vol. 2, 28, no. 2, 57-62 On the Decomposition Method for System of Liner Fredholm Integrl Equtions of the Second Kind A. R. Vhidi 1 nd M. Mokhtri Deprtment of Mthemtics, Shhr-e-Rey
More informationA Bernstein polynomial approach for solution of nonlinear integral equations
Avilble online t wwwisr-publictionscom/jns J Nonliner Sci Appl, 10 (2017), 4638 4647 Reserch Article Journl Homepge: wwwtjnscom - wwwisr-publictionscom/jns A Bernstein polynomil pproch for solution of
More informationApplication of Exp-Function Method to. a Huxley Equation with Variable Coefficient *
Interntionl Mthemticl Forum, 4, 9, no., 7-3 Appliction of Exp-Function Method to Huxley Eqution with Vrible Coefficient * Li Yo, Lin Wng nd Xin-Wei Zhou. Deprtment of Mthemtics, Kunming College Kunming,Yunnn,
More informationNumerical Integration
Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the
More informationReview of Gaussian Quadrature method
Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge
More informationProblem Set 3 Solutions
Msschusetts Institute of Technology Deprtment of Physics Physics 8.07 Fll 2005 Problem Set 3 Solutions Problem 1: Cylindricl Cpcitor Griffiths Problems 2.39: Let the totl chrge per unit length on the inner
More informationAPPLICATIONS OF THE DEFINITE INTEGRAL
APPLICATIONS OF THE DEFINITE INTEGRAL. Volume: Slicing, disks nd wshers.. Volumes by Slicing. Suppose solid object hs boundries extending from x =, to x = b, nd tht its cross-section in plne pssing through
More informationChapter 4 Contravariance, Covariance, and Spacetime Diagrams
Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz
More informationCzechoslovak Mathematical Journal, 55 (130) (2005), , Abbotsford. 1. Introduction
Czechoslovk Mthemticl Journl, 55 (130) (2005), 933 940 ESTIMATES OF THE REMAINDER IN TAYLOR S THEOREM USING THE HENSTOCK-KURZWEIL INTEGRAL, Abbotsford (Received Jnury 22, 2003) Abstrct. When rel-vlued
More informationQuantum Physics II (8.05) Fall 2013 Assignment 2
Quntum Physics II (8.05) Fll 2013 Assignment 2 Msschusetts Institute of Technology Physics Deprtment Due Fridy September 20, 2013 September 13, 2013 3:00 pm Suggested Reding Continued from lst week: 1.
More informationStudies on Nuclear Fuel Rod Thermal Performance
Avilble online t www.sciencedirect.com Energy Procedi 1 (1) 1 17 Studies on Nucler Fuel od herml Performnce Eskndri, M.1; Bvndi, A ; Mihndoost, A3* 1 Deprtment of Physics, Islmic Azd University, Shirz
More informationAdomian Decomposition Method with Green s. Function for Solving Twelfth-Order Boundary. Value Problems
Applied Mthemticl Sciences, Vol. 9, 25, no. 8, 353-368 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/.2988/ms.25.486 Adomin Decomposition Method with Green s Function for Solving Twelfth-Order Boundry
More informationTHE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM
ROMAI J., v.9, no.2(2013), 173 179 THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM Alicj Piseck-Belkhyt, Ann Korczk Institute of Computtionl Mechnics nd Engineering,
More informationResearch Article Numerical Treatment of Singularly Perturbed Two-Point Boundary Value Problems by Using Differential Transformation Method
Discrete Dynmics in Nture nd Society Volume 202, Article ID 57943, 0 pges doi:0.55/202/57943 Reserch Article Numericl Tretment of Singulrly Perturbed Two-Point Boundry Vlue Problems by Using Differentil
More informationEuler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )
Euler, Iochimescu nd the trpezium rule G.J.O. Jmeson (Mth. Gzette 96 (0), 36 4) The following results were estblished in recent Gzette rticle [, Theorems, 3, 4]. Given > 0 nd 0 < s
More informationJournal of Inequalities in Pure and Applied Mathematics
Journl of Inequlities in Pure nd Applied Mthemtics GENERALIZATIONS OF THE TRAPEZOID INEQUALITIES BASED ON A NEW MEAN VALUE THEOREM FOR THE REMAINDER IN TAYLOR S FORMULA volume 7, issue 3, rticle 90, 006.
More informationKINEMATICS OF RIGID BODIES
KINEMTICS OF RIGI OIES Introduction In rigid body kinemtics, e use the reltionships governing the displcement, velocity nd ccelertion, but must lso ccount for the rottionl motion of the body. escription
More informationTHREE DIMENSIONAL MAGNETOSTATIC ANALYSIS OF A 3-PHASE, 12/8 POLE 250 W SWITCHED RELUCTANCE MOTOR USED IN A WASHING MACHINE DRIVE
THREE DIMENSIONAL MAGNETOSTATIC ANALYSIS OF A 3-PHASE, /8 POLE 50 W SWITCHED RELUCTANCE MOTOR USED IN A WASHING MACHINE DRIVE Ahmet FENERCİOĞLU Ünl KURT Gziosmnpş University, Turhl Higher Voctionl School,
More informationAPPROXIMATE LIMIT CYCLES FOR THE RAYLEIGH MODEL
ROMAI J, 4, 228, 73 8 APPROXIMATE LIMIT CYCLES FOR THE RAYLEIGH MODEL Adelin Georgescu, Petre Băzăvn, Mihel Sterpu Acdemy of Romnin Scientists, Buchrest Deprtment of Mthemtics nd Computer Science, University
More information(See Notes on Spontaneous Emission)
ECE 240 for Cvity from ECE 240 (See Notes on ) Quntum Rdition in ECE 240 Lsers - Fll 2017 Lecture 11 1 Free Spce ECE 240 for Cvity from Quntum Rdition in The electromgnetic mode density in free spce is
More informationCIT and QIT with a Periodic Impulsional Potential
Adv. Studies Theor. Phys., Vol. 6, 2012, no. 6, 285-291 nd with Periodic Impulsionl Potentil S. Seddighi Chhrborj srkhosh@mth.upm.edu.my Islmic Azd University, Science nd Reserch Brnch Deprtment of Mthemtics,
More informationAn improvement to the homotopy perturbation method for solving integro-differential equations
Avilble online t http://ijimsrbiucir Int J Industril Mthemtics (ISSN 28-5621) Vol 4, No 4, Yer 212 Article ID IJIM-241, 12 pges Reserch Article An improvement to the homotopy perturbtion method for solving
More informationOn the Uncertainty of Sensors Based on Magnetic Effects. E. Hristoforou, E. Kayafas, A. Ktena, DM Kepaptsoglou
On the Uncertinty of Sensors Bsed on Mgnetic Effects E. ristoforou, E. Kyfs, A. Kten, DM Kepptsoglou Ntionl Technicl University of Athens, Zogrfou Cmpus, Athens 1578, Greece Tel: +3177178, Fx: +3177119,
More informationElectromagnetism Answers to Problem Set 10 Spring 2006
Electromgnetism 76 Answers to Problem Set 1 Spring 6 1. Jckson Prob. 5.15: Shielded Bifilr Circuit: Two wires crrying oppositely directed currents re surrounded by cylindricl shell of inner rdius, outer
More informationProblems for HW X. C. Gwinn. November 30, 2009
Problems for HW X C. Gwinn November 30, 2009 These problems will not be grded. 1 HWX Problem 1 Suppose thn n object is composed of liner dielectric mteril, with constnt reltive permittivity ɛ r. The object
More information1 Bending of a beam with a rectangular section
1 Bending of bem with rectngulr section x3 Episseur b M x 2 x x 1 2h M Figure 1 : Geometry of the bem nd pplied lod The bem in figure 1 hs rectngur section (thickness 2h, width b. The pplied lod is pure
More informationExact solutions for nonlinear partial fractional differential equations
Chin. Phys. B Vol., No. (0) 004 Exct solutions for nonliner prtil frctionl differentil equtions Khled A. epreel )b) nd Sleh Omrn b)c) ) Mthemtics Deprtment, Fculty of Science, Zgzig University, Egypt b)
More informationAn inverse steady state thermal stresses in a thin clamped circular plate with internal heat generation
Americn Journl of Engineering Reserch (AJER) e-issn : 2320-0847 p-issn : 2320-0936 Volume-02, Issue-10, pp-276-281 www.jer.org Reserch Pper Open Access An inverse stedy stte therml stresses in thin clmped
More informationODE: Existence and Uniqueness of a Solution
Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =
More informationConsequently, the temperature must be the same at each point in the cross section at x. Let:
HW 2 Comments: L1-3. Derive the het eqution for n inhomogeneous rod where the therml coefficients used in the derivtion of the het eqution for homogeneous rod now become functions of position x in the
More informationTangent Line and Tangent Plane Approximations of Definite Integral
Rose-Hulmn Undergrdute Mthemtics Journl Volume 16 Issue 2 Article 8 Tngent Line nd Tngent Plne Approximtions of Definite Integrl Meghn Peer Sginw Vlley Stte University Follow this nd dditionl works t:
More informationIntegral points on the rational curve
Integrl points on the rtionl curve y x bx c x ;, b, c integers. Konstntine Zeltor Mthemtics University of Wisconsin - Mrinette 750 W. Byshore Street Mrinette, WI 5443-453 Also: Konstntine Zeltor P.O. Box
More informationArithmetic Mean Derivative Based Midpoint Rule
Applied Mthemticl Sciences, Vol. 1, 018, no. 13, 65-633 HIKARI Ltd www.m-hikri.com https://doi.org/10.1988/ms.018.858 Arithmetic Men Derivtive Bsed Midpoint Rule Rike Mrjulis 1, M. Imrn, Symsudhuh Numericl
More informationOn the Linear Stability of Compound Capillary Jets
ILASS Americs, th Annul Conference on Liquid Atomiztion nd Spry Systems, Chicgo, IL, My 7 On the Liner Stbility of Compound Cpillry Jets Mksud (Mx) Ismilov, Stephen D Heister School of Aeronutics nd Astronutics,
More information63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1
3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =
More informationMath 8 Winter 2015 Applications of Integration
Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl
More informationSynchronous Motor Design using Particle Swarm Optimization Technique
Proceedings o the 4 th Interntionl Middle Est Power Systems Conerence (MEPCON ), Ciro University, Egypt, December 9-,, Pper ID 87. Synchronous Motor Design using Prticle Swrm Optimiztion Technique R. A.
More informationCandidates must show on each answer book the type of calculator used.
UNIVERSITY OF EAST ANGLIA School of Mthemtics My/June UG Exmintion 2007 2008 ELECTRICITY AND MAGNETISM Time llowed: 3 hours Attempt FIVE questions. Cndidtes must show on ech nswer book the type of clcultor
More informationDisclaimer: This Final Exam Study Guide is meant to help you start studying. It is not necessarily a complete list of everything you need to know.
Disclimer: This is ment to help you strt studying. It is not necessrily complete list of everything you need to know. The MTH 33 finl exm minly consists of stndrd response questions where students must
More informationMATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1
MATH34032: Green s Functions, Integrl Equtions nd the Clculus of Vritions 1 Section 1 Function spces nd opertors Here we gives some brief detils nd definitions, prticulrly relting to opertors. For further
More informationKirchhoff and Mindlin Plates
Kirchhoff nd Mindlin Pltes A plte significntly longer in two directions compred with the third, nd it crries lod perpendiculr to tht plne. The theory for pltes cn be regrded s n extension of bem theory,
More informationBernoulli Numbers Jeff Morton
Bernoulli Numbers Jeff Morton. We re interested in the opertor e t k d k t k, which is to sy k tk. Applying this to some function f E to get e t f d k k tk d k f f + d k k tk dk f, we note tht since f
More information- 5 - TEST 2. This test is on the final sections of this session's syllabus and. should be attempted by all students.
- 5 - TEST 2 This test is on the finl sections of this session's syllbus nd should be ttempted by ll students. Anything written here will not be mrked. - 6 - QUESTION 1 [Mrks 22] A thin non-conducting
More information5.7 Improper Integrals
458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the
More informationExam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1
Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution
More information