Abstract. Keywords. 1. Introduction. Eric S. Sanches, José A. Santisteban

Size: px
Start display at page:

Download "Abstract. Keywords. 1. Introduction. Eric S. Sanches, José A. Santisteban"

Transcription

1 Engineering, 2014, 6, Published Online October 2014 in SciRes. Comprtive Study of Conventionl, Fuzzy Logic nd Neurl PID Speed Controllers with Torque Ripple Minimiztion for n Axil Mgnetic Flux Switched Reluctnce Motor Eric S. Snches, José A. Sntistebn Electricl Engineering Deprtment, Mechnicl Engineering Postgrdute Progrm (PGMEC), Fluminense Federl University, Niterói, Brzil Emil: ericsnches@vm.uff.br, jsntistebn@vm.uff.br Received 31 July 2014; revised 25 August 2014; ccepted 8 September 2014 Copyright 2014 by uthors nd Scientific Reserch Publishing Inc. This work is licensed under the Cretive Commons Attribution Interntionl License (CC BY). Abstrct Three speed controllers for n xil mgnetic flux switched reluctnce motor with only one sttor, re described nd experimentlly tested. As it is known, when current pulses re imposed in their windings, high ripple torque is obtined. In order to reduce this ripple, control strtegy with modified current shpes is proposed. A workbench consisting of mchine prototype nd the control system bsed on microcontroller ws built. These controllers were: conventionl PID, fuzzy logic PID nd neurl PID type. From experimentl results, the effective reduction of the torque ripple ws confirmed nd the performnce of the controllers ws compred. Keywords Axil Flux SRM, PID Speed Controller, Fuzzy Logic PID Speed Controller, Neurl PID Speed Controller, Torque Ripple Minimiztion, Current Shpe Control Strtegy 1. Introduction The rdil mgnetic flux version of switched reluctnce motors (SRM) hs been widely used in mny vrible speed industril pplictions nd some dvntges hve been reported: high torque output, wide rnge of operting speed, geometricl simplicity, relibility nd robustness [1] [2]. In the cse of n xil flux SRM (AFSRM) the ir gp mgnetic flux is prllel to the rottion xis. As the AFSRM cn hve smller xis length thn the rdil mgnetic flux motor, it is believed tht it is good solution for use in pplictions where size is importnt, How to cite this pper: Snches, E.S. nd Sntistebn, J.A. (2014) Comprtive Study of Conventionl, Fuzzy Logic nd Neurl PID Speed Controllers with Torque Ripple Minimiztion for n Axil Mgnetic Flux Switched Reluctnce Motor. Engineering, 6,

2 E. S. Snches, J. A. Sntistebn s in electric crs, for instnce. The switched reluctnce motor (SRM) nlyses re complex becuse of their doubly slient pole structure nd nonliner mgnetic chrcteristics. The developed torque is nonliner function of the currents pplied to the sttor windings nd their inductnces, which depend on the rotor position. Nevertheless, with n pproprite control system, minimum ripple torque cn be obtined. From literture, it cn be observed tht mny works re relted to rdil flux SRM; however few works re relted to xil flux SRM [3]. Some of them refer to xil flux SRM with two sttors [4]-[6], but in this pper the motor hs only one sttor. With respect to the control strtegies, severl strtegies re reported [7]-[9]. For ordinry rdil flux SRM, the mutul inductnces of the sttor windings re considered smll [10] nd most reserchers do not tke them into ccount. For AFSRM this is not true becuse the mgnetic flux pths re quite different to those found in rdil flux SRM. In this work, the AFSRM self nd mutul inductnces were estimted bsed on the three-dimensionl finite element method nd then they were used in the motor electromechnicl model. When voltges or current pulses re imposed in the sttor windings, the level of torque ripple is high, if compred with other kind of motors, which is the primry disdvntge of SRM s it contributes to the vibrtion nd coustic noise. For this reson, the imposition of different current shpes to reduce the torque ripple ppers s n interesting solution. In this work, the procedure to obtin lterntives shpes is described nd the effectiveness of this strtegy is experimentlly tested in n AFSRM prototype. Even more, in order to evlute the performnce of different kinds of speed controllers, three PID types were implemented: conventionl, fuzzy logic nd neurl bsed pproch. Although in different levels the experimentl results confirm the reduction of torque ripple, there re lso differences with respect to the execution time nd the speed response. 2. The Axil Flux SRM Prototype The prototype hs six poles in the sttor, corresponding to three phses (, b nd c ), nd four poles in the rotor. It is 3-phse 6/4 poles AFSRM, shown in Figure 1. The chrcteristics of the 3-phse 6/4 AFSRM prototype re shown in Tble 1. Tble 1. The chrcteristics of the 3-phse 6/4 xil flux SRM. Prmeter Outer dimeter of rotor nd sttor Inner dimeter of rotor nd sttor poles Shft dimeter Air gp Sttor pole width Rotor pole width Sttor nd rotor poles rc Poles rdil length Sttor yoke thickness Vlue 126 mm 63 mm 40 mm 1.9 mm 34 mm 26 mm mm 5 mm Sttor pole re (xil cross-section) 1039 mm 2 Rotor yoke thickness 17 mm Number of turns per sttor pole 175 Turn wire Coil resistnce Sttor nd rotor cores mteril Motor shft mteril 24 AWG 2.3 Ω Steel SAE-1020 Stinless steel 656

3 E. S. Snches, J. A. Sntistebn Figure 1. The AFSRM prototype. Ech sttor pole hs one coil, so there re six coils nmed 1, 2, b1, b2, c1 nd c2. The sttor nd rotor cores were solid becuse the objective of this work ws the torque ripple minimiztion nd not the loss reduction. 3. Torque Ripple Minimiztion Strtegy Considering tht the motor core is operting in the liner region, the net electromgnetic torque T e for the three phse AFSRM prototype is found s: ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) L11 θ L12 θ L13 θ I θ 1 Te( θ, I, Ib, Ic) = I( θ) Ib( θ) Ic( θ) L21 θ L22 θ L23 θ Ib θ 2 θ L31 θ L32 θ L33 θ Ic θ L = + L + M + M + L (2) L = + M M M + M (3) 12 b 1 1 b21 b 1 2 b22 L = M + M + M M (4) 13 c 1 1 c21 c 1 2 c22 L = + M M M + M (5) 21 b 1 1 2b1 b 1 2 2b2 L = + L + M + M + L (6) 22 bb 1 1 b b2b1 b bb 1 2 b b2b2 b L = + M M M + M (7) 23 cb 1 1 b c2b1 b cb 1 2 b c2b2 b L = M + M + M M (8) 31 c 1 1 2c1 c 1 2 2c2 L = + M M M + M (9) 32 bc 1 1 b b2c1 b bc 1 2 b b2c2 b L = + L + M + M + L (10) 33 cc 1 1 c c2c1 c cc 1 2 c c2c2 c where: L nd M stnd for the self nd mutul inductnces. The rotor position ngle θ, for phse, is considered s 45 when the rotor pole ws completely overlpping with the sttor pole, s depicted in Figure 2, nd this ngle is mesured clockwise. The θ b nd θ c ngles correspond to rotor ngulr positions for phses b nd c, respectively. The reltionships between these ngles re: for θ 30 θ = θ + 30 b (1) (11) θ = θ + 60 c (12) 657

4 E. S. Snches, J. A. Sntistebn STATOR 90 ROTOR STATOR STATOR ROTOR SHAFT ROTOR STATOR ROTOR STATOR STATOR for 30 < θ 60 for 60 < θ 90 Figure 2. Position of rotor pole 45 reltive to sttor pole of the AFSRM prototype. θ = θ + 30 b θ = θ 30 c θ = θ 60 b (13) (14) (15) θ = θ 30 c The reltionship mong the three phses self nd mutul inductnces re: L 22 (θ ) = L 11 (θ ); L b1b1 (θ b ) = L 11 (θ b ); L b2b2 (θ b ) = L b1b1 (θ b ); L c1c1 (θ c ) = L 11 (θ c ); L c2c2 (θ c ) = L c1c1 (θ c ); M 22 (θ ) = M 12 (θ ); M 2b1 (θ ) = M 1b2 (θ ); M 2b2 (θ ) = M 1b1 (θ ); M 2c1 (θ ) = M 1c2 (θ ); M 2c2 (θ ) = M 1c1 (θ ), M b11 (θ ) = M 1b1 (θ ); M b12 (θ ) = M 2b1 (θ ); M b1b2 (θ b ) = M 12 (θ b ); M b1c1 (θ b ) = M 1b1 (θ b ); M b1c2 (θ b ) = M 1b2 (θ b ); M b21 (θ ) = M 1b2 (θ ); M b22 (θ ) = M 2b2 (θ ); M b2b1 (θ b ) = M b1b2 (θ b ); M b2c1 (θ b ) = M b1c2 (θ b ); M b2c2 (θ b ) = M b1c1 (θ b ); M c11 (θ ) = M 1c1 (θ ); M c12 (θ ) = M 2c1 (θ ); M b2b1 (θ b ) = M b1b2 (θ b ); M c1b1 (θ b ) = M b1c1 (θ b ); M c1b2 (θ b ) = M b2c1 (θ b ); M c1c2 (θ c ) = M 12 (θ c ); M c21 (θ ) = M 1c2 (θ ); M c22 (θ ) = M 2c2 (θ ); M c2b1 (θ b ) = M b1c2 (θ b ); M c2b2 (θ b ) = M b2c2 (θ b ); nd M c2c1 (θ c ) = M c1c2 (θ c ). The self nd mutul inductnces were estimted bsed on simultion model using three dimensionl finite element method of the ANSYS Multiphysics softwre [11] with mgneto sttic nlysis. The 96 type solid elements were used in the simultion softwre. In Figure 3, for eleven rotor positions, the self nd mutul inductnces estimted, for phse, re shown. The self-inductnce is represented by L 11 nd the mutul inductnces re nmed M 12, M 1b1, M 1b2, M 1c1 nd M 1c2, respectively. As shown in Figure 3, the self-inductnce hs similr profile to the rdil flux SRM cse, which mens tht the minimum inductnce corresponds to the unligned position (no overlp between rotor nd sttor poles) while the mximum inductnce corresponds to the full lignment position between rotor nd sttor poles. The inductnces profiles repet every 90. From Figure 3, it cn be noted tht the mutul inductnces between the coil 1 nd the other five coils re not negligible. For instnce, the mximum reltionship between them (M 1b1 /L 11 ) is round 0.25 ner 20. In this wy, in order to provide constnt electromgnetic torque for ll rotor positions, pproprited reference currents should be designed. From previous tests in open loop [12], the best current feed strtegy consists of the following sequence of energizing the windings: (16) 658

5 E. S. Snches, J. A. Sntistebn Figure 3. AFSRM self nd mutul inductnces. phses nd b from rotor position 0 to 15 ; only phse from rotor position 15 to 30 ; phses nd c from rotor position 30 to 45 ; only phse c from rotor position 45 to 60 ; phses b nd c from rotor position 60 to 75 ; only phse b from rotor position 75 to 90. This cycle repets ech 90. On the other hnd, to define reference torque without hrmonics, it ws considered s tht generted when current of 3A is imposed in phse, while the rotor position ws 30. In this sitution net electromgnetic torque of Nm ws clculted. Next, for every one degree, the current vlues of the phse were clculted s: from rotor position 15 to 30 I ( θ ) = L + θ ( θ ) L ( θ ) from rotor position 0 to 15 function tht strts in zero nd it then ssumes the vlue clculted with (17) in the rotor position of 15, ws used ( θ ) I θ ( θ 4) (17) = e (18) from rotor position 30 to 45 the I current vlues re clculted from qudrtic Eqution (19) considering tht I c ssumes the vlues of I currents obtined from (18). For exmple, the I c vlue for 35 is equl to the I vlue for 5 ( ) ( ) ( θ ) L ( θ ) ( ) ( ) L θ L θ 2 L θ 2 L θ ( ) ( ) ( ) c 11 c 12 + I θ + + Ic θ I θ θ θ θ θ Lcc 11 cc Ic ( θ) = 0 θ θ The reference currents of phses c nd b re shifted by 30 nd 60 reltive to phse current. All these reference currents re shown in Figure Controllers Description In order to evlute the performnce of different kinds of speed controllers, three PID types were implemented: conventionl, fuzzy logic nd neurl bsed pproch. (19) 659

6 E. S. Snches, J. A. Sntistebn In Figure 5, the conventionl PID structure is shown. As noted, the current references re obtined multiplying the controller output by tble with three outputs, which depend on the ngulr position of the rotor. In our experiments, two tbles were evluted. In the first one, the outputs hve pulsed shpes. The finl effect is torque with high ripple components. In the second cse, the outputs hve modified shpes ccording to Figure 4. As the experimentl results show, the reduction of torque ripple is importnt. In this cse, the PID gins were djusted by simultions using MtLb Simulink softwre, resulting in gins proportionl, integrl nd derivtive, of 16, 3 nd 1, respectively. The controller output ws limited to the intervl of [0, 1]. In Figure 6, the fuzzy or neurl control structure is shown. As noted, the outputs of the controller re directly the current references for the AFSRM. Figure 4. Reference currents wveforms. +. θ REF -. θ PID CONTROLLER CURRENT WAVEFORM TABLE Reference Currents POWER CONVERTER CIRCUIT Current Sensor ω d/dt θ POSITION SENSOR Figure 5. Block digrm of the PID speed control system. Figure 6. Block digrm of the fuzzy or neurl PID speed control system. 660

7 E. S. Snches, J. A. Sntistebn In the cse of the fuzzy controller, the rules were bsed on the observtion of the simultion results of the PID controller tht use the modified current shpes. The inputs of the fuzzy logic controller were speed error, derivtive of speed error, integrl of speed error nd ngulr rotor position. The output vribles re the three phse s reference currents. In the ctul implementtion, this strtegy ws trnslted to tble. The fuzzy logic chrcteristics used in MtLb Simulink simultions re: minimum And Method, mximum Or Method, minimum Impliction, mximum Aggregtion nd centroid Defuzzifiction. The input vrible speed error ( E ) hs three linguistic vlues: negtive, zero nd positive. The corresponding membership functions re: trpezoidl type [ ] for negtive ( N ); tringulr type [ 1 0 1] zero ( Z ); trpezoidl type [ ] for positive ( P ). The input vrible derivtive of speed error ( CE ) hs three linguistic vlues: negtive, zero nd positive. The corresponding membership functions re: tringulr type [ ] for negtive ( N ); tringulr type [ ] for zero ( Z ); trpezoidl type [ ] for positive ( P ). The input vrible integrl of speed error ( IE ) hs three linguistic vlues: negtive, zero nd positive. The corresponding membership functions re: tringulr type [ ] for negtive ( N ); for zero ( Z ); tringulr type [ ] tringulr type [ ] for positive ( P ). The input vrible ngulr rotor position ( Pos ) hs ten linguistic vlues: P1, P2, P3, P4, P5, P6, P7, P8, P9 nd P10, whose corresponding membership functions re shown in Figure 7. The output vribles re the three phse s reference currents. They hve four linguistic vlues: zero ( Z ), low ( L ), medium ( M ) nd high ( H ). The corresponding membership functions for phse current re shown in Figure 8. Tke into ccount the number of membership functions, there re 270 ( ) fuzzy rules. Its structure is similr to: if (E is P ) nd (CE is P ) nd (IE is N ) nd (Pos is P5 ) then (I is M ) (I b is Z ) (I c is H ). Figure 7. Membership functions of input vribles of the fuzzy speed controller. 661

8 E. S. Snches, J. A. Sntistebn Figure 8. Membership functions of output vribles of the fuzzy speed controller. For the neurl PID controller design, the usul procedure of off-line trining ws dopted. This ws bsed in the dt, with smpling frequency of 1 khz, from the simultion results of the PID controller tht use the modified current shpes. The simultion with MtLb Simulink used 12,000 input dt vlues (speed error, derivtive of speed error, integrl of speed error nd ngulr rotor position) distributed s 6000 for trining, 3000 for test nd 3000 for vlidtion. On the other hnd, 12,000 output dt vlues (currents in three phses) were used, obtined with the simultion of the PID conventionl controller. For neurl network trining the MtLb softwre ws used ( nntool ). The neurl network properties were: feed-forwrd bck propgtion network type, network trining function tht updtes weight nd bis vlues ccording to Levenberg-Mrqurdt optimiztion (TRAINLM), grdient descent with momentum weight nd bis dption lerning function (LEARNGDM), men squred error performnce function (MSE), two lyers: lyer 1 using 36 neurons with hyperbolic tngent sigmoid trnsfer function (TANSIG) nd lyer 2 using 3 neurons with liner trnsfer function (PURELIN). Figure 9 shows the neurl rchitecture nd Figure 10 shows the performnce of the neurl network trining used in the simultion. The trining generted weight nd bis vlues of the two lyers tht were used in the neurl PID controller to clculte the references currents of the simultion. In the ctul implementtion, correspondent tble ws prepred. 5. Controllers Implementtion The hrdwre to control the AFSRM consists of controller circuit, three power converter circuits, one for ech phse, position sensor circuit, three Hll Effect current mesurement circuits, lso one for ech phse, torque meter connected between the AFSRM shft nd the mechnicl lod shft. The shft of DC mchine ws used s n inertil lod. These components re rrnged s shown in the block digrm of Figure 11. The experimentl hrdwre is shown in Figure 12, including the AFSRM. In this photo the prts of the hrdwre used in the experiments cn be identified: the controller circuit (1), the power converter circuits (2), the position sensor circuit (3), Hll Effect current mesurement circuits (4), torque sensor nd disply (5), the AFSRM (6) nd the DC mchine (7). The position sensor consists of n infrred opticl circuit nd n luminum dish with 180 holes, fixed to the rotor structure. The signl from the sensor is trnsformed into rectngulr pulse strem with mplitude of 5 V, 662

9 E. S. Snches, J. A. Sntistebn Figure 9. Neurl rchitecture generted by nntool of the MtLb softwre. Figure 10. Neurl network trining of the MtLb softwre. MEASURED CURRENTS REFERENCE SPEED CONTROLLER REFERENCE CURRENTS POWER CONVERTER AFSRM CONTROLLER REAL SPEED Lod DC MACHINE POSITION SENSOR Figure 11. AFSRM drive block digrm. which is sent to the controller circuit. The ctul phse currents re mesured using commercil Hll Effect circuits. Their gin nd zero djustments were done in order to obtin 1 V for 1 A of phse current. The power circuits consist of three symmetric bridge converters nd uxiliry circuits. In Figure 13, one of the power circuits is depicted. As illustrted, the ctul mesured current is compred with the reference current coming from the controller circuit. If the rel current is greter thn the reference, the power circuit switches off the Mosfets of this phse in order to decrese the current in the coils. The controller circuit consists of one microcontroller PIC 18F4680, progrmmed using the C lnguge, tht receives the rectngulr pulses from the position sensor nd clcultes the reference currents of the three phses 663

10 E. S. Snches, J. A. Sntistebn Figure 12. Experimentl hrdwre nd AFSRM prototype. Power Converter Circuit Reference Current Clock +v Controller Circuit Q Isoltor Mesurement Current Circuit Comprtor _ Q Isoltor SRM Phse Coil Actul Current Figure 13. Power converter block digrm for one phse. in ccordnce with the controller type used, sending them to the power converter. Due to the 180 holes in the luminum dish, the controller imposed the reference currents every Experimentl Results In order to evlute the torque ripple minimiztion strtegy, four experiments were performed: current pulses were imposed nd conventionl PID speed controller ws used; proposed current wveforms, shown in Figure 4, were imposed nd conventionl PID speed controller ws used; fuzzy logic PID speed controller ws used; neurl PID speed controller ws used Imposing Current Pulses In this cse the current pulses were pplied in the following sequence: phse from rotor position 0 to 30 ; phse c from rotor position 30 to 60 ; phse b from rotor position 60 to 90. The Figure 14 shows the reference current imposed by the controller nd the ctul current in the coils of one phse. The Figure 15 shows the speed response nd the Figure 16 shows the torque response, when current pulses re imposed nd conventionl PID speed controller is used. The reference speed is 350 rpm. In this cse, the torque ripple in stedy stte ws round 17.1% of the verge torque. 664

11 E. S. Snches, J. A. Sntistebn Figure 14. Reference pulsed current nd rel current. Figure 15. Speed response for current pulses nd PID conventionl speed controller. Figure 16. Torque response for current pulses nd conventionl PID speed controller Imposing Proposed Current Wveforms The current wveform proposed s reference nd the ctul current in phse re shown in Figure 17. Note tht the ctul nd reference currents re prcticlly overlpping. The Figure 18 shows the speed response nd Figure 19 shows the torque response for current wveforms proposed nd conventionl PID speed controller. In this cse, the torque ripple in stedy stte ws round 2.1% of the verge torque. 665

12 E. S. Snches, J. A. Sntistebn Figure 17. Proposed current reference nd ctul current. Figure 18. Speed response for current wveforms proposed nd conventionl PID speed controller. Figure 19. Torque response for current wveforms proposed nd conventionl PID speed controller Using Fuzzy Logic PID Speed Controller The ctul currents in two phses of the AFSRM re shown in Figure 20. These were very close to their references. Figure 21 shows the speed response nd Figure 22 shows the torque response for the fuzzy logic PID controller. In this cse, the torque ripple in stedy stte ws round 6.5% of the verge torque. 666

13 E. S. Snches, J. A. Sntistebn Figure 20. Actul currents of two phses using the fuzzy logic PID controller. Figure 21. Speed response for the fuzzy logic PID controller. Figure 22. Torque response for the fuzzy logic PID controller Using Neurl PID Speed Controller The ctul currents in two phses of the AFSRM re shown in Figure 23. The Figure 24 shows the speed response nd the Figure 25 shows the torque response for neurl controller. In this cse, the torque ripple in stedy stte ws round 2.0% of the verge torque. 7. Conclusions In this work, prticulr reference wveforms for the currents of n ASFRM in closed loop speed control 667

14 E. S. Snches, J. A. Sntistebn Figure 23. Actul currents of two phses using the neurl PID controller. Figure 24. Speed response for neurl PID controller. Figure 25. Torque response for neurl PID controller. system hve been proposed nd experimentlly tested. Compring the performnce of conventionl PID speed controller using current pulses with conventionl PID speed controller, using the current wveforms proposed, it is confirmed tht the level of torque ripple is reduced. All the controllers were djusted through simultions, using motor model tht is composed by self nd mutul inductnces obtined through 3D FEM simultions. Compring the performnce of conventionl PID speed controller using the proposed current wveforms with fuzzy logic PID speed controller, it ws noted tht the torque ripple for the first one is lower, which cn be explined becuse the currents obtined with fuzzy logic controller re not exctly similr to those tht would 668

15 E. S. Snches, J. A. Sntistebn produce constnt torque. Nevertheless, in respect of the speed response, the fuzzy logic PID controller reches the stedy stte in less time thn the conventionl PID controller. Even more, with respect to the implementtion, the execution time of the fuzzy logic PID is lower thn tht used with the conventionl PID controller but similr to tht spent by the neurl PID controller. Finlly, with respect to the neurl PID controller, it notes similr torque ripple to the conventionl PID controller but lower thn tht obtined using the fuzzy logic pproch, which is expected s this ws obtined from trining bsed on the conventionl PID controller results. However, this hs cost, which is the time spent in off-line trining. References [1] Krishnn, R. (2001) Switched Reluctnce Motor Drives: Modeling, Simultion, Anlysis, Design nd Applictions. CRC Press, Boc Rton. [2] Miller, T.J.E. (1993) Switched Reluctnce Motors nd Their Control. Oxford Science, Oxford. [3] Vijykumr, K., Krthikeyn, R., Prmsivm, S., Arumugm, R. nd Srinivs, K.N. (2008) Switched Reluctnce Motor Modeling, Design, Simultion, nd Anlysis: A Comprehensive Review. IEEE Trnsctions on Mgnetics, 44, [4] Hung, S., Luo, J., Leonrdi, F. nd Lipo, T.A. (2000) A Comprison of Power Density for Axil Flux Mchines Bsed on Generl Purpose Sizing Equtions. IEEE Trnsctions on Energy Conversion, 14, [5] Mdhvn, R. nd Fernndes, B.G. (2013) Axil Flux Segmented SRM with Higher Number of Rotor Segments for Electric Vehicles. IEEE Trnsctions on Energy Conversion, 28, [6] Mdhvn, R. nd Fernndes, B.G. (2012) Comprtive Anlysis of Axil Flux SRM Topologies for Electric Vehicle Appliction. IEEE Interntionl Conference on Power Electronics, Drives nd Energy Systems, Bengluru, December 2012, 1-6. [7] Asghr, K. (2013) Anlysis of Switched Reluctnce Motor Drives for Reduced Torque Ripple Using FPGA Bsed Simultion Technique. Americn Journl of Informtion Sciences, 6, [8] Wdnerkr, V.S., Bhskr, M.M., Ds, T.R. nd RjKumr, A.D. (2010) A New Fuzzy Logic Bsed Modeling nd Simultion of Switched Reluctnce Motor. Journl of Electricl Engineering & Technology, 5, [9] Rj, E.F.I. nd Kmrj, V. (2013) Neurl Network Bsed Control for Switched Reluctnce Motor Drive. IEEE Interntionl Conference on Emerging Trends in Computing, Communiction nd Nnotechnology, [10] Pul, P.P., Silv, W.M., Crdoso, J.R. nd Nbet, S.I. (2003) Assessment of the Influences of the Mutul Inductnces on Switched Reluctnce Mchines Performnce. Interntionl Electric Mchine nd Drives Conference, 3, [11] ANSYS Inc. (2005) Low-Frequency Electromgnetics Anlysis Guide: Relese ANSYS Inc., Cnonsburg. [12] Sss, F., Sntistebn, J.A. nd Snches, E. (2009) Design nd Implementtion of Digitl Control System for n Axil Flux Switched Reluctnce Motor. Brzilin Power Electronics Conference, Bonito, 27 September-1 October 2009,

16

ELE B7 Power Systems Engineering. Power System Components Modeling

ELE 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 information

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

Fig. 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 information

FEM ANALYSIS OF ROGOWSKI COILS COUPLED WITH BAR CONDUCTORS

FEM 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 information

DESIGN OF FERRITE MAGNET CFSM TO IMPROVE SPEED RESPONSE OF EPS SYSTEM

DESIGN 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 information

New Expansion and Infinite Series

New 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 information

#6A&B Magnetic Field Mapping

#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 information

Lec 3: Power System Components

Lec 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 information

ADVANCEMENT 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 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 information

Study on the Calculation of Magnetic Force Based on the Equivalent Magnetic Charge Method

Study 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 information

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

NUMERICAL 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 information

THREE 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, 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 information

MAC-solutions of the nonexistent solutions of mathematical physics

MAC-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 information

A New Method for Estimating Permanent Magnet Synchronous Machine Parameters

A 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 information

CBE 291b - Computation And Optimization For Engineers

CBE 291b - Computation And Optimization For Engineers The University of Western Ontrio Fculty of Engineering Science Deprtment of Chemicl nd Biochemicl Engineering CBE 9b - Computtion And Optimiztion For Engineers Mtlb Project Introduction Prof. A. Jutn Jn

More information

Measuring Electron Work Function in Metal

Measuring Electron Work Function in Metal n experiment of the Electron topic Mesuring Electron Work Function in Metl Instructor: 梁生 Office: 7-318 Emil: shling@bjtu.edu.cn Purposes 1. To understnd the concept of electron work function in metl nd

More information

WELCOME TO THE LECTURE

WELCOME 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 information

Modelling of the near infra-red radiation pulse propagation in biological tissues for medical imaging application

Modelling of the near infra-red radiation pulse propagation in biological tissues for medical imaging application JOURNAL OF INTENSE PULSED LASERS AND APPLICATIONS IN ADVANCED PHYSICS Vol. 3, No. 4, p. 4-45 Modelling of the ner infr-red rdition pulse propgtion in biologicl tissues for medicl imging ppliction A. SAOULI

More information

Estimation of Binomial Distribution in the Light of Future Data

Estimation of Binomial Distribution in the Light of Future Data British Journl of Mthemtics & Computer Science 102: 1-7, 2015, Article no.bjmcs.19191 ISSN: 2231-0851 SCIENCEDOMAIN interntionl www.sciencedomin.org Estimtion of Binomil Distribution in the Light of Future

More information

Driving Cycle Construction of City Road for Hybrid Bus Based on Markov Process Deng Pan1, a, Fengchun Sun1,b*, Hongwen He1, c, Jiankun Peng1, d

Driving Cycle Construction of City Road for Hybrid Bus Based on Markov Process Deng Pan1, a, Fengchun Sun1,b*, Hongwen He1, c, Jiankun Peng1, d Interntionl Industril Informtics nd Computer Engineering Conference (IIICEC 15) Driving Cycle Construction of City Rod for Hybrid Bus Bsed on Mrkov Process Deng Pn1,, Fengchun Sun1,b*, Hongwen He1, c,

More information

Modeling and Simulation of Voltage Source Inverter with Voltage Drop and Its Application for Direct Torque Control of Induction Motors

Modeling and Simulation of Voltage Source Inverter with Voltage Drop and Its Application for Direct Torque Control of Induction Motors Modeling nd Simultion of Voltge Source Inverter with Voltge Drop nd Its Appliction for Direct Torque Control of Induction Motors H. L. Bui,*, Shoudo Hung, D. C. Phm College of Electricl nd Informtion Engineering,

More information

Week 10: Line Integrals

Week 10: Line Integrals Week 10: Line Integrls Introduction In this finl week we return to prmetrised curves nd consider integrtion long such curves. We lredy sw this in Week 2 when we integrted long curve to find its length.

More information

Studies on Nuclear Fuel Rod Thermal Performance

Studies 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 information

Motion of Electrons in Electric and Magnetic Fields & Measurement of the Charge to Mass Ratio of Electrons

Motion of Electrons in Electric and Magnetic Fields & Measurement of the Charge to Mass Ratio of Electrons n eperiment of the Electron topic Motion of Electrons in Electric nd Mgnetic Fields & Mesurement of the Chrge to Mss Rtio of Electrons Instructor: 梁生 Office: 7-318 Emil: shling@bjtu.edu.cn Purposes 1.

More information

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Physics 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 information

4.4 Areas, Integrals and Antiderivatives

4.4 Areas, Integrals and Antiderivatives . res, integrls nd ntiderivtives 333. Ares, Integrls nd Antiderivtives This section explores properties of functions defined s res nd exmines some connections mong res, integrls nd ntiderivtives. In order

More information

In the diagram below, the rotation continues until N-S alignment, resulting in lock-up that is, if nothing is done to prevent it.

In 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 information

Review of Calculus, cont d

Review 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 information

I1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3

I1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3 2 The Prllel Circuit Electric Circuits: Figure 2- elow show ttery nd multiple resistors rrnged in prllel. Ech resistor receives portion of the current from the ttery sed on its resistnce. The split is

More information

LINEAR ALGEBRA APPLIED

LINEAR ALGEBRA APPLIED 5.5 Applictions of Inner Product Spces 5.5 Applictions of Inner Product Spces 7 Find the cross product of two vectors in R. Find the liner or qudrtic lest squres pproimtion of function. Find the nth-order

More information

On the Uncertainty of Sensors Based on Magnetic Effects. E. Hristoforou, E. Kayafas, A. Ktena, DM Kepaptsoglou

On 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 information

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response

More information

Problem Set 3 Solutions

Problem Set 3 Solutions Chemistry 36 Dr Jen M Stndrd Problem Set 3 Solutions 1 Verify for the prticle in one-dimensionl box by explicit integrtion tht the wvefunction ψ ( x) π x is normlized To verify tht ψ ( x) is normlized,

More information

APPLICATIONS OF THE DEFINITE INTEGRAL

APPLICATIONS 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 information

1 Bending of a beam with a rectangular section

1 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 information

POLYPHASE CIRCUITS. Introduction:

POLYPHASE 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 information

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite

More information

High speed machines using advanced magnetic materials analyzed by appropriate finite element models

High 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 information

13.4 Work done by Constant Forces

13.4 Work done by Constant Forces 13.4 Work done by Constnt Forces We will begin our discussion of the concept of work by nlyzing the motion of n object in one dimension cted on by constnt forces. Let s consider the following exmple: push

More information

Math 8 Winter 2015 Applications of Integration

Math 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 information

We partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.

We partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b. Mth 255 - Vector lculus II Notes 4.2 Pth nd Line Integrls We begin with discussion of pth integrls (the book clls them sclr line integrls). We will do this for function of two vribles, but these ides cn

More information

PART 1 MULTIPLE CHOICE Circle the appropriate response to each of the questions below. Each question has a value of 1 point.

PART 1 MULTIPLE CHOICE Circle the appropriate response to each of the questions below. Each question has a value of 1 point. PART MULTIPLE CHOICE Circle the pproprite response to ech of the questions below. Ech question hs vlue of point.. If in sequence the second level difference is constnt, thn the sequence is:. rithmetic

More information

Summary of equations chapters 7. To make current flow you have to push on the charges. For most materials:

Summary of equations chapters 7. To make current flow you have to push on the charges. For most materials: Summry of equtions chpters 7. To mke current flow you hve to push on the chrges. For most mterils: J E E [] The resistivity is prmeter tht vries more thn 4 orders of mgnitude between silver (.6E-8 Ohm.m)

More information

Cross-section section of DC motor. How does a DC Motor work? 2 Commutator Bars N X. DC Motors 26.1

Cross-section section of DC motor. How does a DC Motor work? 2 Commutator Bars N X. DC Motors 26.1 DC Motors 26.1 How does DC Motor work? Crosssection section of DC motor Mgnetic field vector, B oft Iron Core (otor) Wire length vector, dl Force vector, df Current, i Permnent Mgnet (ttor) Crosssection

More information

Acceptance Sampling by Attributes

Acceptance Sampling by Attributes Introduction Acceptnce Smpling by Attributes Acceptnce smpling is concerned with inspection nd decision mking regrding products. Three spects of smpling re importnt: o Involves rndom smpling of n entire

More information

Math& 152 Section Integration by Parts

Math& 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 information

DISTRIBUTION OF SUB AND SUPER HARMONIC SOLUTION OF MATHIEU EQUATION WITHIN STABLE ZONES

DISTRIBUTION OF SUB AND SUPER HARMONIC SOLUTION OF MATHIEU EQUATION WITHIN STABLE ZONES Fifth ASME Interntionl Conference on Multibody Systems, Nonliner Dynmics nd Control Symposium on Dynmics nd Control of Time-Vrying nd Time-Dely Systems nd Structures September 2-2, 05, Long Bech, Cliforni,

More information

Undergraduate Research

Undergraduate Research Undergrdute Reserch A Trigonometric Simpson s Rule By Ctherine Cusimno Kirby nd Sony Stnley Biogrphicl Sketch Ctherine Cusimno Kirby is the dughter of Donn nd Sm Cusimno. Originlly from Vestvi Hills, Albm,

More information

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn

More information

Partial Derivatives. Limits. For a single variable function f (x), the limit lim

Partial Derivatives. Limits. For a single variable function f (x), the limit lim Limits Prtil Derivtives For single vrible function f (x), the limit lim x f (x) exists only if the right-hnd side limit equls to the left-hnd side limit, i.e., lim f (x) = lim f (x). x x + For two vribles

More information

Fourier Series and Their Applications

Fourier Series and Their Applications Fourier Series nd Their Applictions Rui iu My, 006 Abstrct Fourier series re of gret importnce in both theoreticl nd pplied mthemtics. For orthonorml fmilies of complex vlued functions {φ n}, Fourier Series

More information

Maejo International Journal of Science and Technology

Maejo International Journal of Science and Technology Mejo nt. J. Sci. Technol. 4, 8(), 9- Full Pper Mejo nterntionl Journl of Science nd Technology SSN 95-787 vilble online t www.mijst.mju.c.th Protective coordintion of min nd bckup overcurrent relys with

More information

UNIT 1 FUNCTIONS AND THEIR INVERSES Lesson 1.4: Logarithmic Functions as Inverses Instruction

UNIT 1 FUNCTIONS AND THEIR INVERSES Lesson 1.4: Logarithmic Functions as Inverses Instruction Lesson : Logrithmic Functions s Inverses Prerequisite Skills This lesson requires the use of the following skills: determining the dependent nd independent vribles in n exponentil function bsed on dt from

More information

1 Part II: Numerical Integration

1 Part II: Numerical Integration Mth 4 Lb 1 Prt II: Numericl Integrtion This section includes severl techniques for getting pproimte numericl vlues for definite integrls without using ntiderivtives. Mthemticll, ect nswers re preferble

More information

Chaos in drive systems

Chaos in drive systems Applied nd Computtionl Mechnics 1 (2007) 121-126 Chos in drive systems Ctird Krtochvíl, Mrtin Houfek, Josef Koláčný b, Romn Kříž b, Lubomír Houfek,*, Jiří Krejs Institute of Thermomechnics, brnch Brno,

More information

7.2 The Definite Integral

7.2 The Definite Integral 7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where

More information

IRON 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. 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 information

Thermal Performance of Electrocaloric Refrigeration using Thermal Switches of Fluid Motion and Changing Contact Conductance

Thermal Performance of Electrocaloric Refrigeration using Thermal Switches of Fluid Motion and Changing Contact Conductance Americn Journl of Physics nd Applictions 2016; 4(5): 134-139 http://www.sciencepublishinggroup.com/j/jp doi:.11648/j.jp.20160405.12 ISSN: 2330-4286 (Print); ISSN: 2330-4308 (Online) Therml Performnce of

More information

Multi-objective optimization of dielectric layer photonic crystal filter

Multi-objective optimization of dielectric layer photonic crystal filter Optic Applict, Vol. XLVII, No. 1, 017 DOI: 10.577/o170103 Multi-objective optimiztion of dielectric lyer photonic crystl filter HONGWEI YANG *, CUIYING HUANG, SHANSHAN MENG College of Applied Sciences,

More information

Practical exercise 7. Surge tank

Practical exercise 7. Surge tank Prcticl exercise 7. Surge tnk Introduction Surge tnk is used t hydro power plnts for reduction of wter hmmer tht occurs t closing of the turbine inlet vlve. In this exercise we will mesure mss osciltions

More information

u( t) + K 2 ( ) = 1 t > 0 Analyzing Damped Oscillations Problem (Meador, example 2-18, pp 44-48): Determine the equation of the following graph.

u( t) + K 2 ( ) = 1 t > 0 Analyzing Damped Oscillations Problem (Meador, example 2-18, pp 44-48): Determine the equation of the following graph. nlyzing Dmped Oscilltions Prolem (Medor, exmple 2-18, pp 44-48): Determine the eqution of the following grph. The eqution is ssumed to e of the following form f ( t) = K 1 u( t) + K 2 e!"t sin (#t + $

More information

The Trapezoidal Rule

The Trapezoidal Rule _.qd // : PM Pge 9 SECTION. Numericl Integrtion 9 f Section. The re of the region cn e pproimted using four trpezoids. Figure. = f( ) f( ) n The re of the first trpezoid is f f n. Figure. = Numericl Integrtion

More information

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

Exam 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

MATH 144: Business Calculus Final Review

MATH 144: Business Calculus Final Review MATH 144: Business Clculus Finl Review 1 Skills 1. Clculte severl limits. 2. Find verticl nd horizontl symptotes for given rtionl function. 3. Clculte derivtive by definition. 4. Clculte severl derivtives

More information

Songklanakarin Journal of Science and Technology SJST R1 Thongchan. A Modified Hyperbolic Secant Distribution

Songklanakarin Journal of Science and Technology SJST R1 Thongchan. A Modified Hyperbolic Secant Distribution A Modified Hyperbolic Secnt Distribution Journl: Songklnkrin Journl of Science nd Technology Mnuscript ID SJST-0-0.R Mnuscript Type: Originl Article Dte Submitted by the Author: 0-Mr-0 Complete List of

More information

Application of Exact Discretization for Logistic Differential Equations to the Design of a Discrete-Time State-Observer

Application of Exact Discretization for Logistic Differential Equations to the Design of a Discrete-Time State-Observer 5 Proceedings of the Interntionl Conference on Informtion nd Automtion, December 58, 5, Colombo, Sri Ln. Appliction of Exct Discretiztion for Logistic Differentil Equtions to the Design of Discrete-ime

More information

Predict Global Earth Temperature using Linier Regression

Predict Global Earth Temperature using Linier Regression Predict Globl Erth Temperture using Linier Regression Edwin Swndi Sijbt (23516012) Progrm Studi Mgister Informtik Sekolh Teknik Elektro dn Informtik ITB Jl. Gnesh 10 Bndung 40132, Indonesi 23516012@std.stei.itb.c.id

More information

The momentum of a body of constant mass m moving with velocity u is, by definition, equal to the product of mass and velocity, that is

The momentum of a body of constant mass m moving with velocity u is, by definition, equal to the product of mass and velocity, that is Newtons Lws 1 Newton s Lws There re three lws which ber Newton s nme nd they re the fundmentls lws upon which the study of dynmics is bsed. The lws re set of sttements tht we believe to be true in most

More information

Simulated Performance of Packed Bed Solar Energy Storage System having Storage Material Elements of Large Size - Part I

Simulated Performance of Packed Bed Solar Energy Storage System having Storage Material Elements of Large Size - Part I The Open Fuels & Energy Science Journl, 2008, 1, 91-96 91 Open Access Simulted Performnce of Pcked Bed Solr Energy Storge System hving Storge Mteril Elements of Lrge Size - Prt I Rnjit Singh *,1, R.P.

More information

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl

More information

Section 4.8. D v(t j 1 ) t. (4.8.1) j=1

Section 4.8. D v(t j 1 ) t. (4.8.1) j=1 Difference Equtions to Differentil Equtions Section.8 Distnce, Position, nd the Length of Curves Although we motivted the definition of the definite integrl with the notion of re, there re mny pplictions

More information

DIRECT CURRENT CIRCUITS

DIRECT CURRENT CIRCUITS DRECT CURRENT CUTS ELECTRC POWER Consider the circuit shown in the Figure where bttery is connected to resistor R. A positive chrge dq will gin potentil energy s it moves from point to point b through

More information

Realistic Method for Solving Fully Intuitionistic Fuzzy. Transportation Problems

Realistic Method for Solving Fully Intuitionistic Fuzzy. Transportation Problems Applied Mthemticl Sciences, Vol 8, 201, no 11, 6-69 HKAR Ltd, wwwm-hikricom http://dxdoiorg/10988/ms20176 Relistic Method for Solving Fully ntuitionistic Fuzzy Trnsporttion Problems P Pndin Deprtment of

More information

13: Diffusion in 2 Energy Groups

13: Diffusion in 2 Energy Groups 3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups

More information

Credibility Hypothesis Testing of Fuzzy Triangular Distributions

Credibility Hypothesis Testing of Fuzzy Triangular Distributions 666663 Journl of Uncertin Systems Vol.9, No., pp.6-74, 5 Online t: www.jus.org.uk Credibility Hypothesis Testing of Fuzzy Tringulr Distributions S. Smpth, B. Rmy Received April 3; Revised 4 April 4 Abstrct

More information

MIXED MODELS (Sections ) I) In the unrestricted model, interactions are treated as in the random effects model:

MIXED MODELS (Sections ) I) In the unrestricted model, interactions are treated as in the random effects model: 1 2 MIXED MODELS (Sections 17.7 17.8) Exmple: Suppose tht in the fiber breking strength exmple, the four mchines used were the only ones of interest, but the interest ws over wide rnge of opertors, nd

More information

x = b a N. (13-1) The set of points used to subdivide the range [a, b] (see Fig. 13.1) is

x = b a N. (13-1) The set of points used to subdivide the range [a, b] (see Fig. 13.1) is Jnury 28, 2002 13. The Integrl The concept of integrtion, nd the motivtion for developing this concept, were described in the previous chpter. Now we must define the integrl, crefully nd completely. According

More information

CONTRIBUTIONS TO THE STUDY OF TUBULAR

CONTRIBUTIONS TO THE STUDY OF TUBULAR TECHNICAL UNIVERSITY OF CLUJ-NAPOCA IOANA-CORNELIA VESE CONTRIBUTIONS TO THE STUDY OF TUBULAR LINEAR ELECTRIC ACTUATORS PH.D. THESIS Ph.D. Supervisor: Prof.Dr.ing. Mirce M. RADULESCU ABSTRACT OCTOBER 2010

More information

SOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY 5 NOVEMBER 2014

SOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY 5 NOVEMBER 2014 SOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY 5 NOVEMBER 014 Mrk Scheme: Ech prt of Question 1 is worth four mrks which re wrded solely for the correct nswer.

More information

potentials A z, F z TE z Modes We use the e j z z =0 we can simply say that the x dependence of E y (1)

potentials A z, F z TE z Modes We use the e j z z =0 we can simply say that the x dependence of E y (1) 3e. Introduction Lecture 3e Rectngulr wveguide So fr in rectngulr coordintes we hve delt with plne wves propgting in simple nd inhomogeneous medi. The power density of plne wve extends over ll spce. Therefore

More information

SOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD. Mehmet Pakdemirli and Gözde Sarı

SOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD. Mehmet Pakdemirli and Gözde Sarı Mthemticl nd Computtionl Applictions, Vol., No., pp. 37-5, 5 http://dx.doi.org/.99/mc-5- SOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD Mehmet Pkdemirli nd Gözde

More information

The Velocity Factor of an Insulated Two-Wire Transmission Line

The 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 information

Section 5.1 #7, 10, 16, 21, 25; Section 5.2 #8, 9, 15, 20, 27, 30; Section 5.3 #4, 6, 9, 13, 16, 28, 31; Section 5.4 #7, 18, 21, 23, 25, 29, 40

Section 5.1 #7, 10, 16, 21, 25; Section 5.2 #8, 9, 15, 20, 27, 30; Section 5.3 #4, 6, 9, 13, 16, 28, 31; Section 5.4 #7, 18, 21, 23, 25, 29, 40 Mth B Prof. Audrey Terrs HW # Solutions by Alex Eustis Due Tuesdy, Oct. 9 Section 5. #7,, 6,, 5; Section 5. #8, 9, 5,, 7, 3; Section 5.3 #4, 6, 9, 3, 6, 8, 3; Section 5.4 #7, 8,, 3, 5, 9, 4 5..7 Since

More information

Designing Information Devices and Systems I Anant Sahai, Ali Niknejad. This homework is due October 19, 2015, at Noon.

Designing Information Devices and Systems I Anant Sahai, Ali Niknejad. This homework is due October 19, 2015, at Noon. EECS 16A Designing Informtion Devices nd Systems I Fll 2015 Annt Shi, Ali Niknejd Homework 7 This homework is due Octoer 19, 2015, t Noon. 1. Circuits with cpcitors nd resistors () Find the voltges cross

More information

Linear and Non-linear Feedback Control Strategies for a 4D Hyperchaotic System

Linear and Non-linear Feedback Control Strategies for a 4D Hyperchaotic System Pure nd Applied Mthemtics Journl 017; 6(1): 5-13 http://www.sciencepublishinggroup.com/j/pmj doi: 10.11648/j.pmj.0170601.1 ISSN: 36-9790 (Print); ISSN: 36-981 (Online) Liner nd Non-liner Feedbck Control

More information

Calculus - Activity 1 Rate of change of a function at a point.

Calculus - Activity 1 Rate of change of a function at a point. Nme: Clss: p 77 Mths Helper Plus Resource Set. Copright 00 Bruce A. Vughn, Techers Choice Softwre Clculus - Activit Rte of chnge of function t point. ) Strt Mths Helper Plus, then lod the file: Clculus

More information

LABYRINTH SEALS DYNAMIC COEFFICIENTS AND CRITICAL SPEED RESEARCH BASED ON CHILDS MODEL

LABYRINTH SEALS DYNAMIC COEFFICIENTS AND CRITICAL SPEED RESEARCH BASED ON CHILDS MODEL Abstrct LABYINTH SALS DYNAMIC COFFICINTS AND CITICAL SPD SACH BASD ON CHILDS MODL M Wensheng, Hung Hi, Feng Guoqun AVIC Shenyng Aeroengine eserch Institute, Shenyng, Chin Keywords:lbyrinth sels, dynmics

More information

Freely propagating jet

Freely propagating jet Freely propgting jet Introduction Gseous rectnts re frequently introduced into combustion chmbers s jets. Chemicl, therml nd flow processes tht re tking plce in the jets re so complex tht nlyticl description

More information

Chapter 5 : Continuous Random Variables

Chapter 5 : Continuous Random Variables STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 216 Néhémy Lim Chpter 5 : Continuous Rndom Vribles Nottions. N {, 1, 2,...}, set of nturl numbers (i.e. ll nonnegtive integers); N {1, 2,...}, set of ll

More information

Purpose of the experiment

Purpose of the experiment Newton s Lws II PES 6 Advnced Physics Lb I Purpose of the experiment Exmine two cses using Newton s Lws. Sttic ( = 0) Dynmic ( 0) fyi fyi Did you know tht the longest recorded flight of chicken is thirteen

More information

Supplementary Information for Directional Reflective Surface Formed via Gradient- Impeding Acoustic Meta-surfaces

Supplementary Information for Directional Reflective Surface Formed via Gradient- Impeding Acoustic Meta-surfaces Supplementry Informtion for Directionl Reflective Surfce Formed vi Grdient- Impeding Acoustic Met-surfces Kyungjun Song 1*, Jedo Kim 2, Hur Shin 1, Jun-Hyuk Kwk 1, Seong-Hyun Lee 3,Tesung Kim 4 1 Deprtment

More information

1 Online Learning and Regret Minimization

1 Online Learning and Regret Minimization 2.997 Decision-Mking in Lrge-Scle Systems My 10 MIT, Spring 2004 Hndout #29 Lecture Note 24 1 Online Lerning nd Regret Minimiztion In this lecture, we consider the problem of sequentil decision mking in

More information

expression simply by forming an OR of the ANDs of all input variables for which the output is

expression simply by forming an OR of the ANDs of all input variables for which the output is 2.4 Logic Minimiztion nd Krnugh Mps As we found ove, given truth tle, it is lwys possile to write down correct logic expression simply y forming n OR of the ANDs of ll input vriles for which the output

More information

Entropy ISSN

Entropy ISSN Entropy 006, 8[], 50-6 50 Entropy ISSN 099-4300 www.mdpi.org/entropy/ ENTROPY GENERATION IN PRESSURE GRADIENT ASSISTED COUETTE FLOW WITH DIFFERENT THERMAL BOUNDARY CONDITIONS Abdul Aziz Deprtment of Mechnicl

More information

Department of Mechanical Engineering MECE 551 Final examination Winter 2008 April 16, 9:00 11:30. Question Value Mark

Department of Mechanical Engineering MECE 551 Final examination Winter 2008 April 16, 9:00 11:30. Question Value Mark Deprtment of Mechnicl Engineering MECE 55 Finl exmintion Winter 8 April 6, 9: :3 Notes: You my hve your text book nd one pge formul sheet Electronic devices re not llowed except n pproved clcultor NAME:

More information

Math 131. Numerical Integration Larson Section 4.6

Math 131. Numerical Integration Larson Section 4.6 Mth. Numericl Integrtion Lrson Section. This section looks t couple of methods for pproimting definite integrls numericlly. The gol is to get good pproimtion of the definite integrl in problems where n

More information

1 Which of the following summarises the change in wave characteristics on going from infra-red to ultraviolet in the electromagnetic spectrum?

1 Which of the following summarises the change in wave characteristics on going from infra-red to ultraviolet in the electromagnetic spectrum? Which of the following summrises the chnge in wve chrcteristics on going from infr-red to ultrviolet in the electromgnetic spectrum? frequency speed (in vcuum) decreses decreses decreses remins constnt

More information

1. a) Describe the principle characteristics and uses of the following types of PV cell: Single Crystal Silicon Poly Crystal Silicon

1. a) Describe the principle characteristics and uses of the following types of PV cell: Single Crystal Silicon Poly Crystal Silicon 2001 1. ) Describe the principle chrcteristics nd uses of the following types of PV cell: Single Crystl Silicon Poly Crystl Silicon Amorphous Silicon CIS/CIGS Gllium Arsenide Multijunction (12 mrks) b)

More information

Industrial Electrical Engineering and Automation

Industrial Electrical Engineering and Automation CODEN:LUTEDX/(TEIE-719)/1-7/(7) Industril Electricl Engineering nd Automtion Estimtion of the Zero Sequence oltge on the D- side of Dy Trnsformer y Using One oltge Trnsformer on the D-side Frncesco Sull

More information

4.6 Numerical Integration

4.6 Numerical Integration .6 Numericl Integrtion 5.6 Numericl Integrtion Approimte definite integrl using the Trpezoidl Rule. Approimte definite integrl using Simpson s Rule. Anlze the pproimte errors in the Trpezoidl Rule nd Simpson

More information

INTRODUCTION. The three general approaches to the solution of kinetics problems are:

INTRODUCTION. 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 information