Minimizing torque ripple in a brushless DC motor with fuzzy logic: applied to control rod driving mechanism of MNSR

Minimizing torque ripple in a brushless DC motor with fuzzy logic: applied to control rod driving mechanism of MNSR 1 corresp first-author Corresponding author, mohammad.divandari@gmail.com Corresponding author, mohammad.divandari@gmail.com mohammad.divandari@gmail.com 1 1 1 2 Department of Electrical Engineering, Aliabad Katoul branch, Islamic Azad University, Aliabad Katoul 93451-49417, Iran Department of Electrical Engineering, Aliabad Katoul branch, Islamic Azad University, Aliabad Katoul 93451-49417, Iran Department of Mechanical Engineering, Babol Noshirvani University of Technology, Babol 71167-47148, Iran Department of Mechanical Engineering, Babol Noshirvani University of Technology, Babol 71167-47148, Iran A permanent magnet BLDC (brushless direct current) motor is used to move the control rod of a miniature neutron source reactor (MNSR). The BLDC motor drive is modeled using MATLAB/SIMULINK. Two main parts of the modeling are the inverter switching and the current control. Current control with chopping used to minimize the torque ripple of the MNSR control rod drive. Fuzzy logic current control together with soft chopping control shows the best response of all the three strategies. The prototype drive mechanism has an ATmega32 controller and power MOSFET switches. The simulation results are compared with experimental drive mechanism. Torque ripple BLDC motor Fuzzy logic Control rod MNSR I. INTRODUCTION 1 The permanent magnet brushless direct current (BLDC) motor has been widely used in industrial applications of variable speed drives. The BLDC motor produces torque from the interaction of magnetic flux of the rotor magnets and the current carrying stator conductors of the stator core assembly. Its benefits include high power density, high efficiency, low maintenance, high reliability, simplicity, lower price, and fast dynamic response. Stator current commutation torque ripples, occurring due to the loss of ideal phase current commutation and the methods to minimize ripples, have been reported in literatures. Ashabani et al . reduced torque ripple of BLDC motor by control of the DC link voltage during the commutation time [ 1 1 ]. The intelligence method was the basis of their research, where the magnitude of voltage and commutation time were estimated by a neural network and optimized with particle swarm optimization (PSO) algorithm. In 2011 and 2014, the brushless linear DC motor was proposed to move the nuclear reactor control rod-assemblies as primary and auxiliary driving motor, respectively [ 2 2 , 3 3 ]. In this study, the BLDC motor was employed to move the only control rod of a miniature neutron source reactor (MNSR). Torque ripple of the BLDC motor minimizes with fuzzy logic to satisfy the requirements of fast response and precision movement of the control rod. II. REACTOR DESCRIPTION 1 The MNSR developed by the Chinese Institute of Atomic Energy, is a compact research reactor based on the Canadian SLOWPOKE reactor design. It is of a low power tank-in-pool type and uses highly enriched uranium as fuel, light water as moderator and coolant, and metallic beryllium plates as reflector [ 3 3 ]. The core consists of the fuel cage, one central cadmium absorber with stainless steel cladding control rod, fuel pins (fuel rods), dummy rods, and tie rods, with the generated heat removed through natural circulation [ 4 4 ]. A. Reactor description 2 The control rod drive mechanism (CRDM) is driven by a servo motor, type SDE-45 through a mechanical gear system. The servo motor ensures position control of the control rod, and hence stability of the neutron-flux. The CRDM assembly is mounted on top of the reactor vessel, about 0.6 m above the volume of reactor pool water. The pool top is covered with a Perspex material to protect the pool water from environmental contamination and to reduce evaporation. Although most of the materials in the control rod drive mechanism assembly are made of stainless steel, the servo motor contains corrodible materials [ 5 5 , 6 6 ]. In 2014, an electromagnetic levitation system together with synchronous motor was used to navigate the control rod. The control system was programmed in MATLAB through the open-loop system, closed-loop with state feedback, and closed-loop with state feedback integral tracking. The results from the prototype system showed that the proposed method was useful in such a sensitive drive mechanism [ 4 4 ]. III. MODELING OF BLDC 1 The BLDC motor back electromotive force (back-EMF) induced in stator windings has a trapezoidal shape. In this study, due to high efficiency, cost effective control, optimum number of power electronic devices (number of MOSFETs), and low torque ripple, the three phases of BLDC motor have been used. Fig. 1 Fig. 1 shows the cross section of the three-phase BLDC motor construction. The BLDC drive modeling is implemented based on the following steps: • The phase currents and position are measured; • Differential equations of the BLDC consist of speed, position, and the current of phases a and b; • Trapezoidal back-EMF of BLDC drive simulated using a look-up table; • By means of torque ripple minimization an inverter with soft chopping has been used; • The current control unit provides a DC voltage level of the inverter. Figure 2 Figure 2 shows the schematic of the BLDC drive modeling. A. BLDC differential equations 2 Fundamental of differential equations in BLDC motor can be written as (refer to Fig. 2 Fig. 2 for the variants) [ 7 7 ] <math display="block" id="M1"><mrow><msub><mi>v</mi><mrow><mtext>ab</mtext></mrow></msub><mo>=</mo><mi>R</mi><msub><mi>i</mi><mrow><mtext>ab</mtext></mrow></msub><mo>+</mo><mi>R</mi><mtext>d</mtext><msub><mi>i</mi><mrow><mtext>ab</mtext></mrow></msub><mo>/</mo><mtext>d</mtext><mi>t</mi><mo>+</mo><msub><mi>E</mi><mrow><mtext>ab</mtext></mrow></msub><mn>,</mn></mrow></math> <math display="block" id="M2"><mrow><msub><mi>v</mi><mrow><mtext>bc</mtext></mrow></msub><mo>=</mo><mi>R</mi><msub><mi>i</mi><mrow><mtext>bc</mtext></mrow></msub><mo>+</mo><mi>R</mi><mtext>d</mtext><msub><mi>i</mi><mrow><mtext>bc</mtext></mrow></msub><mo>/</mo><mtext>d</mtext><mi>t</mi><mo>+</mo><msub><mi>E</mi><mrow><mtext>bc</mtext></mrow></msub><mn>,</mn></mrow></math> <math display="block" id="M3"><mrow><msub><mi>v</mi><mrow><mtext>ca</mtext></mrow></msub><mo>=</mo><mi>R</mi><msub><mi>i</mi><mrow><mtext>ca</mtext></mrow></msub><mo>+</mo><mi>R</mi><mtext>d</mtext><msub><mi>i</mi><mrow><mtext>ca</mtext></mrow></msub><mo>/</mo><mtext>d</mtext><mi>t</mi><mo>+</mo><msub><mi>E</mi><mrow><mtext>ca</mtext></mrow></msub><mn>,</mn></mrow></math> <math display="block" id="M4"><mrow><msub><mi>i</mi><mrow><mtext>ab</mtext></mrow></msub><mo>=</mo><msub><mi>i</mi><mtext>a</mtext></msub><mo>−</mo><msub><mi>i</mi><mtext>b</mtext></msub><mn>,</mn><mtext> </mtext><msub><mi>E</mi><mrow><mtext>ab</mtext></mrow></msub><mo>=</mo><msub><mi>e</mi><mtext>a</mtext></msub><mo>−</mo><msub><mi>e</mi><mtext>b</mtext></msub><mn>,</mn></mrow></math> <math display="block" id="M5"><mrow><msub><mi>i</mi><mrow><mtext>bc</mtext></mrow></msub><mo>=</mo><msub><mi>i</mi><mtext>b</mtext></msub><mo>−</mo><msub><mi>i</mi><mtext>c</mtext></msub><mn>,</mn><mtext> </mtext><msub><mi>E</mi><mrow><mtext>bc</mtext></mrow></msub><mo>=</mo><msub><mi>e</mi><mtext>b</mtext></msub><mo>−</mo><msub><mi>e</mi><mtext>c</mtext></msub><mn>,</mn></mrow></math> <math display="block" id="M6"><mrow><msub><mi>i</mi><mrow><mtext>ca</mtext></mrow></msub><mo>=</mo><msub><mi>i</mi><mtext>c</mtext></msub><mo>−</mo><msub><mi>i</mi><mtext>a</mtext></msub><mn>,</mn><mtext> </mtext><msub><mi>E</mi><mrow><mtext>ca</mtext></mrow></msub><mo>=</mo><msub><mi>e</mi><mtext>c</mtext></msub><mo>−</mo><msub><mi>e</mi><mtext>a</mtext></msub><mn>,</mn></mrow></math> <math display="block" id="M7"><mrow><msub><mi>i</mi><mtext>a</mtext></msub><mo>+</mo><msub><mi>i</mi><mtext>b</mtext></msub><mo>+</mo><msub><mi>i</mi><mtext>c</mtext></msub><mo>=</mo><mn>0,</mn></mrow></math> where v ab , v bc and v ca are line-to-line voltage in the three phases a, b and c; i ab , i bc and i ca are line-to-line current in the three phases; i a , i b and i c are current of the three phases; E ab , E bc and E ca are line-to-line back-EMF in the three phases; e a , e b and e c are voltage of the three phases; R is the phase resistance ( Ω ). The back-EMF for three-phases are [ 6 6 ] <math display="block" id="M8"><mrow><msub><mi>e</mi><mtext>a</mtext></msub><mo>=</mo><msub><mi>k</mi><mrow><mi>e</mi><mi>ω</mi></mrow></msub><mi>K</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><mn>,</mn></mrow></math> <math display="block" id="M9"><mrow><msub><mi>e</mi><mtext>b</mtext></msub><mo>=</mo><msub><mi>k</mi><mrow><mi>e</mi><mi>ω</mi></mrow></msub><mi>K</mi><mo stretchy="false">(</mo><mi>θ</mi><mo>−</mo><mn>120</mn><mo stretchy="false">)</mo><mn>,</mn></mrow></math> <math display="block" id="M10"><mrow><msub><mi>e</mi><mtext>c</mtext></msub><mo>=</mo><msub><mi>k</mi><mrow><mi>e</mi><mi>ω</mi></mrow></msub><mi>K</mi><mo stretchy="false">(</mo><mi>θ</mi><mo>−</mo><mn>240</mn><mo stretchy="false">)</mo><mn>,</mn></mrow></math> <math display="block" id="M11"><mrow><mi>K</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><mo>=</mo><mrow><mo>{</mo><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mn>1</mn></mtd><mtd columnalign="left"><mrow><mtext> </mtext><mn>0</mn><mo>≤</mo><mi>θ</mi><mo><</mo><mn>120</mn></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mn>1</mn><mo>−</mo><mn>6</mn><mo stretchy="false">(</mo><mi>θ</mi><mo>−</mo><mn>120</mn><mo stretchy="false">)</mo><mo>/</mo><mi>π</mi></mrow></mtd><mtd columnalign="left"><mrow><mtext> </mtext><mn>120</mn><mo>≤</mo><mi>θ</mi><mo><</mo><mn>180</mn></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mo>−</mo><mn>1</mn></mrow></mtd><mtd columnalign="left"><mrow><mtext> </mtext><mn>180</mn><mo>≤</mo><mi>θ</mi><mo><</mo><mn>300</mn></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mn>1</mn><mo>−</mo><mn>6</mn><mo stretchy="false">(</mo><mi>θ</mi><mo>−</mo><mn>300</mn><mo stretchy="false">)</mo><mo>/</mo><mi>π</mi></mrow></mtd><mtd columnalign="left"><mrow><mtext> </mtext><mn>300</mn><mo>≤</mo><mi>θ</mi><mo><</mo><mn>360</mn></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow/></mtd></mtr></mtable><mn>,</mn></mrow></mrow></mrow></math> <math display="block" id="M12"><mrow><msub><mi>k</mi><mrow><mi>e</mi><mi>ω</mi></mrow></msub><mo>=</mo><msub><mi>k</mi><mi>e</mi></msub><msub><mi>ω</mi><mtext>m</mtext></msub><mo>/</mo><mn>2,</mn></mrow></math> where, k e is back-EMF constant (V s/rad), ω is angular velocity (rad/s), θ is the angular displacement (radian) and K ( θ ) is trapezoidal waveform of the back-EMF. The electro-mechanical torque can be expressed as [ 7 7 ] <math display="block" id="M13"><mrow><msub><mi>T</mi><mtext>e</mtext></msub><mo>−</mo><msub><mi>T</mi><mtext>L</mtext></msub><mo>=</mo><msub><mi>k</mi><mtext>f</mtext></msub><msub><mi>ω</mi><mtext>m</mtext></msub><mo>+</mo><mi>J</mi><mtext>d</mtext><msub><mi>ω</mi><mtext>m</mtext></msub><mo>/</mo><mtext>d</mtext><mi>t</mi><mn>,</mn></mrow></math> where, T e is electromagnetic torque (N m), T L is load torque (N m), k f is motor friction constant (N s/rad) and J is rotor inertia (kg m 2 ). The BLDC torque is a function of the current and the back-EMF <math display="block" id="M14"><mrow><msub><mi>T</mi><mtext>e</mtext></msub><mo>=</mo><msub><mi>k</mi><mtext>t</mtext></msub><mo stretchy="false">[</mo><mi>K</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><msub><mi>i</mi><mtext>a</mtext></msub><mo>+</mo><mi>K</mi><mo stretchy="false">(</mo><mi>θ</mi><mo>−</mo><mn>120</mn><mo stretchy="false">)</mo><msub><mi>i</mi><mtext>b</mtext></msub><mo>+</mo><mi>K</mi><mo stretchy="false">(</mo><mi>θ</mi><mo>−</mo><mn>240</mn><mo stretchy="false">)</mo><msub><mi>i</mi><mtext>c</mtext></msub><mo stretchy="false">]</mo><mo>/</mo><mn>2,</mn></mrow></math> where k t is torque constant (N m/A). From Eqs. (2) and (7), one has <math display="block" id="M15"><mrow><msub><mi>v</mi><mrow><mtext>bc</mtext></mrow></msub><mo>=</mo><mi>R</mi><mo stretchy="false">(</mo><msub><mi>i</mi><mtext>a</mtext></msub><mo>+</mo><mn>2</mn><msub><mi>i</mi><mtext>b</mtext></msub><mo stretchy="false">)</mo><mo>+</mo><mi>L</mi><mtext>d</mtext><mo stretchy="false">(</mo><msub><mi>i</mi><mtext>a</mtext></msub><mo>+</mo><mn>2</mn><msub><mi>i</mi><mtext>b</mtext></msub><mo stretchy="false">)</mo><mtext>d</mtext><mi>t</mi><mo>+</mo><msub><mi>E</mi><mrow><mtext>bc</mtext></mrow></msub><mn>,</mn></mrow></math> where L is phase inductance (H). According to Eqs. (1), (2) and (15), the BLDC motor differential equations for three-phase can be reduced to two-phase as below <math display="block" id="M16"><mrow><mtext>d</mtext><msub><mi>i</mi><mtext>a</mtext></msub><mo>/</mo><mtext>d</mtext><mi>t</mi><mo>=</mo><mo>−</mo><mi>R</mi><msub><mi>i</mi><mtext>a</mtext></msub><mo>/</mo><mi>L</mi><mo>+</mo><mn>2</mn><mo stretchy="false">(</mo><msub><mi>v</mi><mrow><mtext>ab</mtext></mrow></msub><mo>−</mo><msub><mi>E</mi><mrow><mtext>ab</mtext></mrow></msub><mo stretchy="false">)</mo><mo>/</mo><mo stretchy="false">(</mo><mn>3</mn><mi>L</mi><mo stretchy="false">)</mo><mo>+</mo><mo stretchy="false">(</mo><msub><mi>v</mi><mrow><mtext>bc</mtext></mrow></msub><mo>−</mo><msub><mi>E</mi><mrow><mtext>bc</mtext></mrow></msub><mo stretchy="false">)</mo><mo>/</mo><mo stretchy="false">(</mo><mn>3</mn><mi>L</mi><mo stretchy="false">)</mo><mn>,</mn></mrow></math> <math display="block" id="M17"><mrow><mtext>d</mtext><msub><mi>i</mi><mtext>b</mtext></msub><mo>/</mo><mtext>d</mtext><mi>t</mi><mo>=</mo><mo>−</mo><mi>R</mi><msub><mi>i</mi><mtext>b</mtext></msub><mo>/</mo><mi>L</mi><mo>−</mo><mo stretchy="false">(</mo><msub><mi>v</mi><mrow><mtext>ab</mtext></mrow></msub><mo>−</mo><msub><mi>E</mi><mrow><mtext>ab</mtext></mrow></msub><mo stretchy="false">)</mo><mo>/</mo><mo stretchy="false">(</mo><mn>3</mn><mi>L</mi><mo stretchy="false">)</mo><mo>+</mo><mo stretchy="false">(</mo><msub><mi>v</mi><mrow><mtext>bc</mtext></mrow></msub><mo>−</mo><msub><mi>E</mi><mrow><mtext>bc</mtext></mrow></msub><mo stretchy="false">)</mo><mo>/</mo><mo stretchy="false">(</mo><mn>3</mn><mi>L</mi><mo stretchy="false">)</mo><mn>,</mn></mrow></math> Equation (13) can be simplified to <math display="block" id="M18"><mrow><mtext>d</mtext><msub><mi>ω</mi><mtext>m</mtext></msub><mo>/</mo><mtext>d</mtext><mi>t</mi><mo>=</mo><mo stretchy="false">(</mo><msub><mi>T</mi><mtext>e</mtext></msub><mo>−</mo><msub><mi>T</mi><mtext>L</mtext></msub><mo>−</mo><msub><mi>k</mi><mtext>f</mtext></msub><msub><mi>ω</mi><mtext>m</mtext></msub><mo stretchy="false">)</mo><mo>/</mo><mi>J</mi><mn>,</mn></mrow></math> where <math display="block" id="M19"><mrow><mtext>d</mtext><msub><mi>θ</mi><mtext>m</mtext></msub><mo>/</mo><mtext>d</mtext><mi>t</mi><mo>=</mo><msub><mi>ω</mi><mtext>m</mtext></msub><mn>.</mn></mrow></math> B. Inverter of BLDC drive 2 The BLDC motor drive was actuated by a six-switch inverter as shown in Fig. 3 Fig. 3 . The switches are metal oxide semi-conductor field effect transistor (MOSFET) [ 7 7 ]. In an inverter, the current in active phase was turned on and off by the electronic devices, usually at a high frequency. It is called chopping. Chopping can be done by hard chopping and soft chopping. In hard chopping, the upper and lower switches are driven by the same signal; while in soft chopping, only the upper switch is driven by the chopping signal and the lower switch is left on during the whole interval. The advantage of soft chopping is that it creates less switching losses and current ripple than the hard chopping, hence the minimized torque ripple. Proposed soft chopping circuit is shown in Fig. 4 Fig. 4 . The inverter voltages in soft chopping mode applied to six-switch inverter are given in Table 1 Table 1 . 1 1 Rotor position 1 1 Phase current (A) 1 1 v ab - e ab 1 1 v bc - e bc 1 1 0°–60° 1 1 i a ≠0, i c ≠0 1 1 - e a + e b 1 1 - e b + e c 1 1 1 1 i a ≠0, i c =0 1 1 - e a + e b 1 1 0.5( e a - e b ) 1 1 1 1 i a =0, i c ≠0 1 1 0.5( e b - e c ) 1 1 - e b + e c 1 1 1 1 i a =0, i c =0 1 1 0 1 1 0 1 1 60°–120° 1 1 i a ≠0, i b ≠0 1 1 -V s - e a + e b 1 1 V s - e b + e c 1 1 1 1 i a ≠0, i b =0 1 1 0.5(- e a + e c ) 1 1 0.5(- e a + e c ) 1 1 1 1 i a =0, i b ≠0 1 1 0.5(- V s + e b - e c ) 1 1 V s - e b + e c 1 1 1 1 i a =0, i b =0 1 1 0 1 1 0 1 1 120°–180° 1 1 i b ≠0, i a ≠0 1 1 - e a + e b 1 1 - e b + e c 1 1 1 1 i b ≠0, i a =0 1 1 0.5( e b - e c 1 1 - e b + e c 1 1 1 1 i b =0, i a ≠0 1 1 0.5(- e a + e c 1 1 0.5(- e a + e c 1 1 1 1 i b =0, i a =0 1 1 0 1 1 0 1 1 180°–240° 1 1 i b ≠0, i c ≠0 1 1 - e a + e b 1 1 -V s - e b + e c 1 1 1 1 i b ≠0, i c =0 1 1 - e a + e b 1 1 0.5( e a - e b 1 1 1 1 i b =0, i c ≠0 1 1 0.5(- V s - e a + e c ) 1 1 0.5(- V s - e a + e c ) 1 1 1 1 i b =0, i c =0 1 1 0 1 1 0 1 1 240°–300° 1 1 i c ≠0, i b ≠0 1 1 - e a + e b 1 1 - e b + e c 1 1 1 1 i c ≠0, i b =0 1 1 0.5(- e a + e c ) 1 1 0.5 (- e a + e c ) 1 1 1 1 i c =0, i b ≠0 1 1 - e a + e b 1 1 0.5( e a - e b 1 1 1 1 i c =0, i b =0 1 1 0 1 1 0 1 1 300°–360° 1 1 i c ≠0, i a ≠0 1 1 V s - e a + e b 1 1 - e b + e c 1 1 1 1 i c ≠0, i a =0 1 1 0.5( e b - e c 1 1 - e b + e c 1 1 1 1 i c =0, i a ≠0 1 1 V s - e a + e b 1 1 0.5(- V s + e a - e b ) 1 1 1 1 i c =0, i a =0 1 1 0 1 1 0 C. Current control of BLDC drive 2 There are two main sources of torque ripple generated by the BLDC drives: chopping in inverter and current control. By means of torque ripple minimization, a new soft chopping control has been presented in this study, and the current control in the BLDC drive via hysteresis current control (HCC) and fuzzy logic current control (FLCC) was investigated. 1. Hysteresis current control 3 Reviewing researches in the past decade about torque ripple reduction BLDC drive, most authors worked on the torque ripple reduction by HCC. Current control of hysteresis band is a simple and usual current closed-loop control. In HCC, the value of controlled variable is forced to stay within limits (hysteresis band) around a reference value. HCC technique generates switching frequencies when a narrow hysteresis is used, and the wider the hysteresis band is, the larger the ripple is produced. Fig. 5 Fig. 5 is the block diagram of HCC. The uncertain switching frequency makes filtering of electromagnetic noise difficult. Moreover, current chattering produces noise in electronic devices and generates high mechanical stress. 2. Fuzzy logic current control 3 The torque ripple of BLDC drive could be minimized by using fuzzy logic current control. The block diagram of FLCC for minimization of torque ripple of BLDC drive is shown in Fig. 6 Fig. 6 . Current error is the differential between current reference and root mean square (RMS) of three phase current. In FLCC method, fuzzy rules act on current errors (five member function for -0.5–0.5 A) and current error derivative (seven member function for -0.06–0.06 A/ms) as inputs, and also act on inverter voltage level (five members function for 0–1 V) as output. Fig. 7 Fig. 7 shows the fuzzy logic member functions for proposed FLCC. Table 2 Table 2 shows the fuzzy logic rules used to reduce the torque ripple. 1 2 Current Error 7 1 Current error derivative 1 1 NB 1 1 NM 1 1 NS 1 1 Z 1 1 PS 1 1 PM 1 1 PB 1 1 NM 1 1 NS 1 1 NM 1 1 NM 1 1 NM 1 1 NM 1 1 NM 1 1 NM 1 1 NS 1 1 NS 1 1 NS 1 1 Z 1 1 Z 1 1 NS 1 1 NS 1 1 NM 1 1 Z 1 1 NM 1 1 NS 1 1 NS 1 1 Z 1 1 Z 1 1 NS 1 1 NS 1 1 PS 1 1 PM 1 1 PM 1 1 PS 1 1 PS 1 1 Z 1 1 NS 1 1 NS 1 1 PM 1 1 PM 1 1 PM 1 1 PM 1 1 PM 1 1 PM 1 1 PM 1 1 PS IV. SIMULATION RESULTS AND DISCUSSION 1 The proposed BLDC drive which applied to the CRDM of MNSR has been simulated by MATLAB/ SIMULINK. The parameters of the model are listed in table 3 table 3 . 1 1 Parameters 1 1 Value 1 1 Number of poles 1 1 2 1 1 Nominal Power (W) 1 1 1.2 1 1 Nominal voltage (V) 1 1 6 1 1 No load speed (RPM) 1 1 8000 1 1 Terminal resistance phase to phase (Ω) 1 1 12.50 1 1 Terminal inductance phase to phase (H) 1 1 0.091×10 -3 1 1 No load current (A) 1 1 0.5 1 1 Back-EMF constant (V s/rad) 1 1 1.05×10 -3 1 1 Torque constant (N m A) 1 1 1.05×10 -3 1 1 Rotor inertia (kg m 2 ) 1 1 0.005×10 -7 1 1 Friction constant (N s/rad) 1 1 1.38×10 -8 1 1 MNSR control rod mass (g) 1 1 50 1 1 Control rod length (mm) 1 1 230 1 1 Control rod speed limit (cm/s) 1 1 1.06 1 1 Gearbox ratio 1 1 1/22 1 1 Pulley diameter (cm) 1 1 10 1 1 MOSFET model of power switches 1 1 IRF840 1 1 Micro-controller 1 1 ATmega32 The BLDC drive has been operated by three cases: a) HCC with hard chopping, b) HCC with soft chopping and c) FLCC with soft chopping. Fig. 8 Fig. 8 shows the fast dynamic response of the BLDC motor where the speed of the motor was regulated at 360 RPM. The simulation results of the BLDC drive operation are shown Fig. 9 Fig. 9 . From Fig. 9 Fig. 9 , the total torque of the BLDC motor is affected by current control. The hard chopping control has a good response to HCC rather than soft chopping control with HCC. However, the best performance can be obtained using soft chopping with FLCC. This torque ripple reduction technique can be used in such a precision movement of the MNSR control rod drive. A prototype of the optimized drive mechanism with BLDC motor is shown in Fig. 10 Fig. 10 . This setup includes the BLDC motor, drive, gearbox, and the moving rod. The drive consists of the three parts; power switches (6 MOSFETs), micro-controller, and interface circuits. V. CONCLUSION 1 The BLDC motor drive mechanism has been investigated to be used in a research reactor CRDM. The MATLAB/SIMULINK was used to simulate the three understudy cases: HCC with hard chopping, HCC with soft chopping, and FLCC with soft chopping. The results showed the combination of FLCC method with soft chopping as an improved technique for reduction of torque ripple has a good performance. This minimized torque ripple of the BLDC motor drive will directly effect on the movement of the MNSR control rod. In nuclear industrial application, frequent control rod movements for load-following operation induce xenon-oscillation. Therefore, such an FLCC controller together with the soft chopping that can subdue this phenomenon effectively is needed. References M Ashabani , Kaviani A Kashefi , J Milimonfared , et al . Minimization of commutation torque ripple in brushless DC motors with optimized input voltage control . International Symposium on Power Electronics, Electrical Drives, Automation and Motion, SPEEDAM , Italy , 2008 . Minimization of commutation torque ripple in brushless DC motors with optimized input voltage control K Paul . US patent 8811562 , Control rod drive mechanism for nuclear reactor , 2011 . Control rod drive mechanism for nuclear reactor P K David , L T Stafford , C C Brian . US patent 88634513 , System for exchanging a component of a nuclear reactor , 2012 . System for exchanging a component of a nuclear reactor M Divandari , F Shahsavari , M Rezaeimoghadam , et al . A novel control-rod drive mechanism via electromagnetic levitation in MNSR . Nukleonika , 2014 , 59 : 73 - 79 . DOI: 10.2478/nuka-2014-0008 http://doi.org/10.2478/nuka-2014-0008 A novel control-rod drive mechanism via electromagnetic levitation in MNSR I Yusuf , A A Mati . Effects of extended shutdown on the control rod drive mechanism of Nigeria research reactor-1(NIRR-1) . African conference on research reactor utilization and safety , Nigeria , 2010 . Effects of extended shutdown on the control rod drive mechanism of Nigeria research reactor-1(NIRR-1) International atomic energy agency technology and use of low power research reactors. IAEA-TECDOC-384 , 1986 . International atomic energy agency technology and use of low power research reactors. IAEA-TECDOC-384 M Divandari . et al . A novel fuzzy dynamic observer for high speed BLDC motor . IEEE XXXIII International Scientific Conference Electronics and Nanotechnology (ELNANO) , Ukraine , 2013 . A novel fuzzy dynamic observer for high speed BLDC motor 10.13538/j.1001-8042/nst.26.010601 Supported by Research Contract of the Islamic Azad University’s Aliabad Katoul branch 2014-08-11 2014-11-01 2015-02-20 2015-02-20 2021-04-11 2015 Nuclear Science and Techniques 1 26