Next Article in Journal / Special Issue
Optimal Passive Dynamics for Physical Interaction: Catching a Mass
Previous Article in Journal / Special Issue
Optically Driven Mobile Integrated Micro-Tools for a Lab-on-a-Chip
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis and Modeling of Linear-Switched Reluctance for Medical Application

1
University of Toulouse, INPT, UPS; ENSEEIHT, 2 rue Charles Camichel, BP 7122, F-31071 Toulouse cedex 7, France
2
National Center for Scientific Research (CNRS); LAPLACE; F-31071 Toulouse, France
3
Novatem Simplified joint-stock company (SAS); ENSEEIHT, 27 rue d’Aubuisson F-31071 Toulouse cedex 7, France
4
Cardiology Institute, Pitié-Salpêtrière Hospital, 47-83 Bd de l'Hôpital, 75651 Paris Cedex 13, France
*
Author to whom correspondence should be addressed.
Actuators 2013, 2(2), 27-44; https://doi.org/10.3390/act2020027
Submission received: 18 February 2013 / Revised: 25 March 2013 / Accepted: 7 April 2013 / Published: 22 April 2013
(This article belongs to the Special Issue Human Centered Actuators)

Abstract

:
This paper focuses on the analysis, the modeling and the control of a linear-switched reluctance motor. The application under consideration is medical, and the actuator is to be used as a left ventricular assist device. The actuator has a cylindrical or tubular shape, with a mechanical unidirectional valve placed inside the mover, which provides a pulsatile flow of blood. The analytical expression of the effort based on the linear behavior of the actuator is given. The identification of the characteristics of the prototype actuator and the principle of position control is performed. A modeling of the actuator is proposed, taking into account the variation of inductance with respect to the position. The closed-loop position control of the actuator is performed by simulation. A controller with integral action and anticipatory action is implemented in order to compensate the effects of disturbing efforts and tracking deviations. Moreover, a magic switch is performed in the controller to avoid overshoots. The results show that the closed-loop response of the actuator is satisfactory.

1. Introduction

Heart failure is one of the most common diseases in developed countries. In most of the cases, a left ventricular assist device (LVAD) is used to treat patients. The LVAD assists the failing left ventricle by providing an additional flow of blood in the body [1].
At the beginning of cardiovascular surgery, mechanical circulatory support (MCS), which is a pneumatic pump, mimics the function of the heart by providing a pulsatile flow of blood. The major limitations are the sheer size of the pump, limited mobility, because of a large drive console, and the need for constant anticoagulation solutions [1]. The latest generation of MCS consists of miniature rotary pumps providing a continuous flow of blood. One can cite the Jarvik 2000, the MicroMed DeBakey and the Heart Mate II, which are axial electromagnetic pumps, including a motor and a turbine [2]. They increase the autonomy of the patient and improve overall life comfort. The success of these devices is limited by thromboembolic events and pump thrombosis. Moreover, the implications of continuous flow conditions over long periods remain unclear to the medical and scientific community [3].
To prevent thrombosis, we propose to study a pulsatile flow pump for in vivo tests on pigs. The physiological needs are a pressure of 120 mmHg, a flow of 3 liters/min and a heart rate of 120 bpm (beats per minute). The actuator chosen is a linear-switched reluctance motor (LSRM). The stator and the mover have a cylindrical or tubular shape, with a mechanical unidirectional valve placed inside the mover. This pump, called Pulsamag, incorporates the valve and pump functions.
The aims of the paper are to describe the characteristics and the modeling of this LSRM structure in order to control it and then to show the applicability of the motor for LVAD. The organization of the paper is as follows: the topology of the proposed LSRM and its basic operational principles are described in Section 2; Section 3 gives the dimensional characteristics of the motor; Section 4 presents the modeling of the motor; and Section 5 describes the conclusions and perspectives of the work.

2. Proposed Concept of the LVAD

2.1. Presentation of the Actuator

For the application considered, the cylindrical LSRM is a good candidate. Its structure is simple and robust, and the mover part has no windings, as in [4,5]. It consists of a hollow, mobile ferromagnetic cylinder with transversal slots and an outer cylindrical stator containing a number of identical phase sets. Each set consists of a magnetic core, formed by two discs and a ring and an enclosed solenoid winding. In the same manner, a six phase cylindrical LSRM has been developed for circulatory heart assistance [6].
The actuator studied is called Pulsamag and is presented in [7]; it has a mechanical unidirectional valve fixed inside the cylindrical mover. This aortic valve, made by St. Jude Medical, has a diameter of 25 mm and allows blood circulation. As shown in Figure 1, the cylindrical LSRM studied is composed of a stator with four phases and eight coils (two coils per phase) and a tubular mover. Each set of the stator is separated from its neighbor by a non-magnetic ring. Phase number 1 is composed of coils, numbers 1 and 5. The other coils and phases are numbered in the same way.
The coils are sequentially excited in order to create a magnetic field moving from one end to the other. Then the cylindrical magnetic mover performs an oscillatory movement. The closed aortic valve pushes the fluid in the flow direction, and when the cylinder comes back to its initial position, the valve is opened, so as not to interact with the fluid.
Figure 1. (a) Basic structure of Pulsamag;(b) Design of Pulsamag.
Figure 1. (a) Basic structure of Pulsamag;(b) Design of Pulsamag.
Actuators 02 00027 g001
In this low speed application, because of low eddy current losses, the magnetic core does not have to be made of laminated steel. Moreover, due to the cylindrical shape, the normal forces are neutralized. The other advantages of this kind of actuator are its reliability, its simplicity of implementation, its ability to generate a linear motion directly (without mechanical processing), etc. Its disadvantages include low efficiency, due to a loss of joules that increases proportionally to the desired mechanical effort generated, and the difficulty of guiding the mover part to ensure a constant air gap.

2.2. Kinetic Characteristics of the Motor

According to [7], the stroke of the St. Jude Medical valve and the mover is:
Actuators 02 00027 i001
where L is the stroke of the mover in mm, Q is the desired flow in L/min, S is the section of the valve in cm² and f is the desired frequency in Hz.
The thrust of the mover is:
Actuators 02 00027 i002
where F is the thrust in N and P is the pressure in kPa.
The physiological needs for a ventricular assist device are a flow of 3 L/min and a pressure of 120 mmHg (16kPa) under a frequency of 2 Hz (120 bpm). Therefore, with a 25 mm diameter valve, the stroke is L ≅ 50mm and the thrust is F ≅ 8 N.
According to Equation (1), the control of the position of the mover at a given frequency will fix the blood flow. By setting a value of supply voltage, a current is injected in each phase. Thus, the driving force required is produced and, in accordance with Equation (2), blood pressure is imposed.

3. Dimensional Characteristics of the LVAD

3.1. Prototype Dimensions and Characteristics

From the actuator presented in Figure 1, it can be noted that each phase set forms an independent magnetic circuit thanks to the non-magnetic ring. The phase set or basic pattern of the actuator is presented in Figure 2(a). It is formed of a stator ferromagnetic part with one slot (with a coil) and two teeth and a ferromagnetic mover part with one slot (without a coil) and two teeth. When the coil is excited, the stator teeth and the mover teeth tend to be aligned in order to maximize the magnetic flux.
The radial dimensions (y-axis) are: stator yoke thickness, esy, stator tooth length, hst, mover tooth length, hmt, mover yoke thickness, emy, external radius, Rext, air gap radius, Rg, and valve radius, Rv. The axial dimensions (x-axis) are: tooth width, a, and slot width, b.
Figure 2. (a) Basic pattern;(b) Picture of Pulsamag prototype.
Figure 2. (a) Basic pattern;(b) Picture of Pulsamag prototype.
Actuators 02 00027 g002
The sizing and the two-dimensional finite element analysis of the actuator are presented in [7,8]. In order to ensure smooth functioning, we chose the tooth pitch of the stator, λs and the tooth pitch of the mover, λm, to be equal, the tooth width, a, and the slot width, b, to be equal, too:
Actuators 02 00027 i003
We can notice in Figure 1(a) that to separate two basic patterns, a non-magnetic ring has to be inserted. Its thickness value, c, depends on the number of phases and on the tooth pitch values. For the Pulsamag prototype, we chose to make the thickness of the non-magnetic ring, c, half that of the tooth width, a.
The stroke of the actuator depends on: the tooth pitch length, λ, the non-magnetic ring thickness, c, the phase number, m, and the coils per phase number, n. To increase the stroke for a given elementary pattern dimensions, we can increase m and/or n. By supplying several coils in series (n ≠ 1), both the stroke and the thrust increase.
The total length of the stator, LsT, is given by:
Actuators 02 00027 i004
The total length of the mover, LmT, depends on the desired stroke, L, and is given by:
Actuators 02 00027 i005
By assuming a linear model of the basic pattern, the analytical expression of the traction or propulsion force can be obtained using the principle of conversion of energy [9,10] and by neglecting the magnetic field in the iron relative to that in the air gap [5,8]. According to [8], the expression of the total expected thrust FxT is given by:
Actuators 02 00027 i006
where Ni is the magnetomotive force (mmf) and g the length of the air gap.
It is to be noted that the lateral force will be high if the air gap length, g, is low, and the air gap radius, Rg, is high. Moreover, increasing the mmf, Ni leads to an increase of the thrust, but also a proportional increase in copper loss (square increase).
The actuator presented in Figure 1 has m = 4 phases and n = 2 coils per phase. The first phase is composed of coils 1 and 5 in series, and the other phases are formed and numbered in the same way. Supplying the phases in order 1-2-3-4 by a four-phase current system leads to a leftward shift, whereas the reverse sequence. 4-3-2-1. leads to a rightward shift. The distance covered by the mover when switching from one coil to the next is the mechanical step:
Actuators 02 00027 i007
The prototype dimensions and a picture of the LSRM are given, respectively, in Table 1 and in Figure 2(b).
Table 1. Prototype dimensions.
Table 1. Prototype dimensions.
NameAbbreviationValue (mm)
Valve radiusRv12.5
Mover yoke thicknessemy2
Mover tooth lengthhmt1.5
Air gap lengthg0.2
Air gap radiusRg16.1
Stator tooth lengthhst9
Stator yoke thicknessesy1.5
External radiusRext26.7
Tooth widtha2.9
Slot widthb2.9
Non-magnetic ring thicknessc1.45
Total length of the statorLsT80
Total length of the moverLmT140
For the application considered, the desired stroke and the thrust are, respectively, 50 mm and 8 N. According to the above dimensions, the stroke is available.
The winding characteristics of the prototype are given in Table 2.
Table 2. Winding characteristics.
Table 2. Winding characteristics.
NameAbbreviationValue
Number of phasesm4
Number of coils per phasen2
Number of turns per slotN155
Turn diameterФturn0.335 mm
Slot areaAslot = bhst26.1 mm²
Coil areaAcoil = NπФturn2/413.66 mm²
Slot fill factorkfill = Acoil/Aslot52.30%
The expression of the mmf, Ni, in terms of the previous areas expressed in mm² and the current density, J, expressed in A/mm², is:
Actuators 02 00027 i008
Therefore, for a current supply of i = 1 A corresponding to a mmf of Ni = 155 At and a current density in the coil, J = 11.3 A/mm², the analytical thrust using Equation (6) is equal to 7.6 N. This value is close to that desired for the application. However, it should be remembered that the analytical expression assumes a linear operation, an independence of each magnetic circuit associated with each coil (due to the non-magnetic ring) and a null magnetic field in the iron relative to that in the air gap.
By measurements, it was determined that the prototype had the parameters and characteristics given in Table 3.
Table 3. Prototype characteristics.
Table 3. Prototype characteristics.
CharacteristicValue
Mass of the moverm = 270.8 g
Static dry friction force fs = 1.75 N
Winding resistance (average value)R = 8.5 Ω
Unaligned Position inductance (average value)Lu = 34.1 mH
Aligned Position inductance (average value)La = 44.6 mH
Winding inductance (average value)L = 39.4 mH
Rate of change of inductance (average value)∆L/∆x = 10/2.9

3.2. Input-Output Characteristics

The output power, Pout, is given as the product of the kinetic thrust, FxT, and average speed, Actuators 02 00027 i024, as:
Actuators 02 00027 i009
By supplying one phase after the other, the average speed can be expressed according to the desired frequency, f (i.e., desired heart rate) and the desired stroke, L, as:
Actuators 02 00027 i010
In this low speed application (the frequency is low, typically under 3 Hz), the iron losses can be neglected compared to the copper losses. Therefore, the efficiency, η, is given by:
Actuators 02 00027 i011
where PJ are the Joule losses.
The Joule losses for one energized phase (constituted of n coils) are:
Actuators 02 00027 i012
where R is the winding resistance, which can be expressed according to the electrical resistivity of copper, ρ, the winding average length of turn, lt, and its cross section, At:
Actuators 02 00027 i013
where:
Actuators 02 00027 i014
and:
Actuators 02 00027 i015
Finally, the Joule losses can be expressed by:
Actuators 02 00027 i016
One dimensioning factor for switched reluctance motors is the maximization of the mean torque per copper loss. In [5], this optimization gives the following result (by neglecting the dimension, g/2, to the Rg dimension in the previous expression):
Actuators 02 00027 i017
If we neglect the stator yoke thickness, esy, compared to the external radius, Rext, the optimum ratio is 0.4 between the air gap radius, Rg, and the external radius, Rext. For our prototype, the external radius, Rext, according to Equation (17), will be larger than 31 mm, considering the radius valve value, Rv, which is equal to 12.5 mm. In order to limit the sheer size of the actuator (external radius inferior to 30 mm) and the mass of the mover, we have chosen the dimensions given in Table 1.
By supplying one phase after the other, we can calculate the output power, the input power and the efficiency at different frequencies (cf.Figure 3). As mentioned previously, because the frequency is low, iron losses are negligible. The Joule losses are considered constant and equal for each phase supplied. We can notice that the efficiency of the LSRM is very low, as in [6]. In this medical application, efficiency is not the main design criterion. The choice of this architecture is guided by the application.
Figure 3. Input power, output power and efficiency versus frequency.
Figure 3. Input power, output power and efficiency versus frequency.
Actuators 02 00027 g003

4. Control Principle

4.1. Open Loop Control of the Motor

4.1.1. Presentation

The motor is controlled using a dSPACE system with a DS1005 PWM card and a DS2202 acquisition card. A Mitutoyo LGF-550L position sensor with the EH10P signal conditioner is used to obtain the real position of the mover. The power supply of the prototype has been developed by the Novatem SAS Company and consists of four independent H-bridges (one H-bridge per phase). Figure 4 shows the test bench of the Pulsamag prototype.
In this pump application, the functioning of the motor requires a back and forth motion. The motor is to be controlled in position with a sinusoidal variation of the position:
Actuators 02 00027 i018
where Actuators 02 00027 i025 is equal to half the desired stroke and f is the operating frequency equal to the desired cardiac frequency (cf. Equation (1)).
Figure 4. Test bench of Pulsamag prototype.
Figure 4. Test bench of Pulsamag prototype.
Actuators 02 00027 g004
The control of the motor has been developed using Matlab Simulink software in discrete mode with a sampling period: Te = 100 µs. The real position of the mover, called xm, measured by the sensor position, will be compared to the reference sinusoidal signal, x, of Equation (18). Figure 5 shows the open loop control of the motor. Concerning the use of Matlab Simulink and dSPACE tools for the control of SRM, we can cite the works [11,12] for variable speed pumping applications.
Figure 5. Open loop control of the motor. PWM, pulse-width modulation.
Figure 5. Open loop control of the motor. PWM, pulse-width modulation.
Actuators 02 00027 g005
The elements that are inside the dotted box in Figure 5 are the material elements of the control. In the “PWM & power supply” block, the DS1005 pulse-width modulation (PWM) card generates the PWM signals needed for the Novatem SAS power supply. This provides voltages, U1 to U4, to the motor. In the “Sensor position” block, the DS2202 card is used to retrieve the real position of the mover, xm, measured by the Mitutoyo sensor position (analog to digital conversion). The other blocks of Figure 5 (presented in next section) are the control blocks, which are located in the PC (in the form of Matlab Simulink blocks).

4.1.2. Determination of the Control Blocks

4.1.2.1. “Position to Speed” Block

By setting a reference signal for the position with a sinusoidal variation, the desired speed has a sinusoidal variation. It is obtained simply by using a derivation block.

4.1.2.2. “Motor Simplified Electric Motor” Block

To control the motor, we must provide duty cycles for the four H-bridges of the inverter to the “PWM & power supply” block. To do this, the “Motor simplified electric model” block receives the currents, i1 to i4, as input and generates the duty cycles, k1 to k4. The content of this block is based on Equation (19) for phase number k:
Actuators 02 00027 i019
We assume that the magnetic circuit of the motor has a linear behavior without saturation and that the phases are independent (no mutual inductance). This equation is commonly used in the works dealing with modeling and control of SRM [11,13].
The values of R, L and dL/dx are given in Section 3.1. The DC bus voltage, Ubus, corresponding to the supply voltage of the H-bridges is fixed at 30 V.

4.1.2.3. “Current Signal Generator” Block

The generation of currents, i1 to i4, is performed by the “Current signal generator” block. This block is an “S-function” of Matlab Simulink. The inputs are the desired total force, F, and the reference signal position, x. Figure 6 illustrates the current supply sequence required for a reference signal position.
Figure 6. Generation of the phase currents according to the reference position.
Figure 6. Generation of the phase currents according to the reference position.
Actuators 02 00027 g006
We can see that when the reference position increases (positive speed), the current supply sequence generated is in the order i1, i2, i3, i4 allowing movement of the mover in the direction of increasing x. On the contrary, when the reference position decreases (negative speed), the current supply sequence generated is in the order i4, i3, i2, i1, allowing movement of the mover in the direction of decreasing x.
The current amplitude is function of the desired force and is determined by:
Actuators 02 00027 i020
It is assumed here that only one phase is energized at a time, and only the phase supplied contributes to generating the desired force. Previous assumptions on the linearity, the non-saturation of the magnetic material and the independence of the phases are maintained. This equation is also used in work [14,15], which are related to the modeling of SRM.

4.1.2.4. “Speed to Force” Block

The desired total force, F, is generated from the reference speed, v, by the “Speed to force” block, which is chosen as a simple gain, KvF, in order to simplify the control. Thus, the speed, v, and the desired total force, F, have the same evolution, which explains why the current waveform shown in Figure 6 is sinusoidal. A preliminary test was performed in order to determine the gain, KvF. During this experiment, one phase of the motor was supplied by a square wave signal of 1 A, and the position of the mover was recorded (cf. Figure 7).
Figure 7. Preliminary test in order to determine KvF.
Figure 7. Preliminary test in order to determine KvF.
Actuators 02 00027 g007
The determination of the slope of the curve, x(t) (position vs. time curve zoomed in Figure 7), and the use of Equation (20) allowed us to determine the gain, KvF = 0.05, where v is expressed in mm/s and F in N.

4.1.4. Results of the Open Loop Control

The motor being modeled, Figure 13 illustrates the open loop response of the motor, xm*, to the desired position, x (10 mm amplitude and 2 Hz frequency).
Figure 13. Open loop configuration: input (position x) and output (position xm*) waveforms.
Figure 13. Open loop configuration: input (position x) and output (position xm*) waveforms.
Actuators 02 00027 g013
The steady state error requires setting up a closed loop control with an integral and derivative terms controller.

4.2. Closed Loop Control of the Motor

From Figure 8, the block diagram of the closed loop control of the motor can be established in Figure 14.
Figure 14. Closed loop control of the motor.
Figure 14. Closed loop control of the motor.
Actuators 02 00027 g014
Thus, a block called “Force Controller” is implemented. From the reference position, x, and the estimated measured position, xm*, this block will generate the appropriate control effort. This corrector includes an integral action and an anticipatory action, called feed-forward, which allows for compensation of the tracking deviations during the speed ramps, and the integral action allows for compensation of the effects of disturbing efforts. To implement all these elements, we used the reference speed, v, and the estimated measured velocity, vm*.
The output of the “Force Controller” block is called the control force, Fc, which is provided to the “Current signal generator” block. We recall that this block must generate the current supply sequence, which allows the mover to achieve the desired position. It has the estimated measured position, xm*, as the second input. Figure 15 shows the implemented “Force Controller” block.
Figure 15. “Force Controller” block.
Figure 15. “Force Controller” block.
Actuators 02 00027 g015
K1 defines the proportional gain, which acts on the position error. The output of this gain is called the control speed, Vc. The integrator gain is chosen equal to ¼ of K1 to obtain a proper integral action of the controller. Then, the control speed is corrected by an anticipatory (feed-forward) action. The gain, K2, defines the proportional gain, which acts on the speed error. The output of the corrector is the control force, Fc. Upstream of the integrator, a “magic-switch” has been inserted. When the switch is open (out = 0), the integral action is held. This block is an “S-function” of Matlab Simulink and has x, xm* and the control speed, Vc, for inputs.
Figure 16. (a) Flow chart of the “magic-switch”;(b) Closed loop configuration: input (position x) and output (position xm*) waveforms.
Figure 16. (a) Flow chart of the “magic-switch”;(b) Closed loop configuration: input (position x) and output (position xm*) waveforms.
Actuators 02 00027 g016
The flow chart of the S-function is presented in Figure 16(a). In principle, the speed error should not continue to be integrated to inflate the command when the dynamics of this error are such that it tends to disappear naturally, otherwise the extra control contributes to overshooting. To do this, signal Vc was compared to the estimated speed, Vest. When the estimated speed, Vest, is too rapid compared to Vc, then the switch is open (out = 0); otherwise, it will be closed (out = Vc).
Figure 16(b) shows the results of the closed loop control for the same input waveform, x, as in Section 4.1.4. We can see that the corrector cancels the steady state error. We can see also that the phase shift is very small between the two waveforms. The corrector thus ensures its function, and the measured output is very close to the input reference.

5. Conclusions

In this paper, we have presented a tubular linear-switched reluctance motor for a medical application. This actuator must provide a pulsatile pump operation for a ventricular assistance device. After describing the operating principle of this actuator, we presented the characteristics of the prototype actuator and the principle of position control. This control is based on a dSPACE control system associated with Matlab Simulink. A dedicated power supply based on a H-bridge inverter for the four phases of the machine has been developed.
A simplified model of the actuator was proposed, which can be generalized to a linear-switched reluctance motor for all applications. This modeling takes into account the variation of inductance with respect to the position and assumes an independence of the phases (and a linear behavior of the magnetic materials).
After this modeling, we present the open- and closed-loop position control of the actuator in Matlab Simulink. For the closed-loop control, a controller with an integral action and anticipatory action, called feed-forward, is implemented. We provide the implementation of a magic switch that prevents the integral action to inflate the command and that avoids overshoots of the output. The results obtained show that the use of a simplified model is satisfactory.
In the perspective of this work, a real control of the prototype under dSPACE based on the Simulink models already developed must be carried out.

References

  1. DiGiorgi, P.L.; Naka, Y.; Oz, M.C. Left Ventricular Assist Devices. Contemporary Cardiology: Surgical Management of Congestive Heart Failure; Springer Humana Press: Totawa, NJ, USA, 2005; Chapter 7; pp. 155–189. [Google Scholar]
  2. Blitz, A.; Fang, J.C. Ventricular Assist Devices and Total Artificial Hearts. In Contemporary Cardiology: Device Therapy in Heart Failure; Springer Humana Press: Boston, MA, USA, 2009; Chapter 13; pp. 339–371. [Google Scholar]
  3. Allen, G.S.; Murray, K.D; Olsen, D.B. The importance of pulsatile and non pulsatile flow in the design of blood pumps. Artif. Organs 1997, 21, 922–928. [Google Scholar] [CrossRef]
  4. Corda, J.; Skopljak, E. Linear Switched Reluctance Actuator. In Proceedings of 6th International Conference on Electrical Machines and Drives; Inspec/IEE: London, UK, 1993; pp. 535–539. [Google Scholar]
  5. El Amraoui, L.; Gillon, F.; Vivier, S.; Brochet, P.; Benrejeb, M. Optimal Design Approach for Linear Tubular Machines. In Proceedings of IEEE International Conference on Systems Man and Cybernetics; IEEE: Danvers, MA, USA, 2002; Volume 5, pp. 1–6. [Google Scholar]
  6. Yamada, H.; Hamajima, T.; Xiang, S.; Nishizawa, N. Six-phase linear pulse motor as linear oscillatory actuator. IEEE Trans. Mag. 1987, MAG-23, 2841–2843. [Google Scholar]
  7. Llibre, J.F.; Martinez, N.; Nogarede, B.; Leprince, P. Linear Tubular Switched Reluctance Motor for heart assistance circulatory: Analytical and finite element modeling. In Proceedings of 10th International Workshop on Electronics, Control, Measurement and Signals (ECMS); IEEE: Danvers, MA, USA, 2011; pp. 1–6. [Google Scholar]
  8. Llibre, J.F.; Martinez, N.; Leprince, P.; Nogarède, B. Innovative Linear Pulsatile Pump for heart assistance circulatory. In Proceedings of The Eighth International Symposium on Linear Drives for Industry Applications (LDIA 2011), Eindhoven, The Netherlands, 3–6 July 2011.
  9. Lee, B.-S.; Bae, H.-K.; Vijayraghavan, P.; Krishnan, R. Design of a linear switched reluctance machine. IEEE Trans. Ind. Appl. 2000, 36, 1571–1580. [Google Scholar] [CrossRef]
  10. Viorel, I.-A.; Hameyer, K.; Strete, L. Transverse flux tubular switched reluctance motor. In Proceedings of 11th International Conference on Optimization of Electrical and Electronic Equipment, OPTIM'08; IEEE: Danvers, MA, USA, 2008; pp. 131–136. [Google Scholar]
  11. Vasquez, H.; Parker, J.K. A new simplified mathematical model for a switched reluctance motor in a variable speed pumping application. Mechatronics 2004, 14, 1055–1068. [Google Scholar] [CrossRef]
  12. Vasquez, H.; Parker, J.K.; Haskew, T. Control of a 6/4 switched reluctance motor in a variable speed pumping application. Mechatronics 2005, 15, 1061–1071. [Google Scholar] [CrossRef]
  13. Krishnan, R. Switched Reluctance Motor Drives, Modeling, Simulation, Analysis, Design, and Applications; CRC Press LCC: Boca Raton, FL, USA, 2001; pp. 21–22. [Google Scholar]
  14. Bae, H.-K.; Lee, B.-S.; Vijayraghavan, P.; Krishnan, R. A linear switched reluctance motor: converter and control. IEEE Trans. Ind. Appl. 2000, 36, 1351–1359. [Google Scholar] [CrossRef]
  15. Lim, H.S.; Krishnan, R.; Lobo, N.S. Design and control of a linear propulsion system for an elevator using linear switched reluctance motor drives. IEEE Trans. Ind. Appl. 2008, 55, 534–542. [Google Scholar]

Share and Cite

MDPI and ACS Style

Llibre, J.-F.; Martinez, N.; Leprince, P.; Nogarede, B. Analysis and Modeling of Linear-Switched Reluctance for Medical Application. Actuators 2013, 2, 27-44. https://doi.org/10.3390/act2020027

AMA Style

Llibre J-F, Martinez N, Leprince P, Nogarede B. Analysis and Modeling of Linear-Switched Reluctance for Medical Application. Actuators. 2013; 2(2):27-44. https://doi.org/10.3390/act2020027

Chicago/Turabian Style

Llibre, Jean-Francois, Nicolas Martinez, Pascal Leprince, and Bertrand Nogarede. 2013. "Analysis and Modeling of Linear-Switched Reluctance for Medical Application" Actuators 2, no. 2: 27-44. https://doi.org/10.3390/act2020027

Article Metrics

Back to TopTop