Next Article in Journal
Kinematics Analysis and Optimization of a 3-DOF Planar Tensegrity Manipulator under Workspace Constraint
Next Article in Special Issue
Human Gait Data Augmentation and Trajectory Prediction for Lower-Limb Rehabilitation Robot Control Using GANs and Attention Mechanism
Previous Article in Journal
Design, Development and Control of a Therapeutic Robot Incorporating Aquatic Therapy for Ankle Rehabilitation
Previous Article in Special Issue
Design and Analysis of a Lower Limb Rehabilitation Training Component for Bedridden Stroke Patients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrodynamic Bearing Structural Design of Blood Pump Based on Axial Passive Suspension Stability Analysis of Magnetic–Hydrodynamic Hybrid Suspension System

Key Laboratory of Advanced Manufacturing and Intelligent Technology, Harbin University of Science and Technology, Harbin 150080, China
*
Author to whom correspondence should be addressed.
Machines 2021, 9(11), 255; https://doi.org/10.3390/machines9110255
Submission received: 27 September 2021 / Revised: 22 October 2021 / Accepted: 25 October 2021 / Published: 27 October 2021
(This article belongs to the Special Issue Design and Control of Advanced Mechatronics Systems)

Abstract

:
Rotor suspension stability is one of the important performance indexes of a blood pump and the basis of determining whether the blood pump can be used in a clinic. Compared with the traditional magnetic suspension system, a single-winding, bearingless motor has the advantages of a compact structure, simple control system and low power consumption. In this pursuit, the present study aimed to envisage and design the magnetic suspension system coupled with a single-winding bearingless motor and permanent magnet bearings, establish the theoretical models of axial force and electromagnetic torque, and calculate the stiffness of the magnetic suspension system at the equilibrium point. Addressing the problem of the negative axial stiffness of the magnetic suspension system being negative, which leads to the instability of the suspension rotor, the hydrodynamic bearing structure was proposed and designed, and the critical stiffness to realize the stable suspension of the rotor was obtained based on the stability criterion of the rotor dynamics model. The optimal structural parameters of the hydrodynamic bearing are selected by integrating various factors based on the solution of the Reynolds equation and a stiffness analysis. Furthermore, the vibration experiment results proved that the blood pump rotor exhibited a good suspension stability, and the maximum offset under the impact external fluid was no more than 2 μm.

1. Introduction

Heart failure is a common cardiovascular disease and implanting a blood pump is one of the effective ways to treat this kind of disease [1,2]. The earlier developed blood pumps exhibited a pulsatile flow blood pump that simulated cardiac contractions and relaxations, but faced problems such as a large volume, high noise, high failure rate, inability of implantation and a risk of bleeding and thrombosis, which seriously affected the clinical treatment effect [3]. Later, researchers developed the continuous-flow rotary blood pump, which greatly reduced its volume and weight and could be implanted into the abdominal cavity of the patients with heart failure. The typical products, such as the HeartMate II blood pump, are certified by the FDA and widely used in the clinical setting [4], but the mechanical bearing of the blood pump wears off and heats up after long-term use, making it easier to cause hemolysis, thrombosis, and other complications [5]. The rotor of the third-generation blood pump adopts a suspension support system that cancels the mechanical bearing structure, eliminates the blood damage caused by mechanical bearing, and has good blood compatibility [6,7]. As one of the remarkable characteristics of the third-generation blood pump, the suspension system of the rotor directly determines the operation stability and blood compatibility of the blood pump in practical applications.
The suspension system of the blood pump includes the magnetic suspension, hydraulic suspension and magnetic hydrodynamic hybrid suspension. For the magnetic–hydrodynamic hybrid suspension blood pump, the combination of the magnetic bearing and hydraulic bearing is generally used to regulate the movement of the rotor, and the rotating motion is driven by a separate motor. Therefore, in order to realize the active control of the suspension and rotation, more windings, inverters and sensor components are required [8]. For the implantable blood pump, the increased volume of the blood pump increases to a certain extent, the weight and power consumption can adversely affect the implantation of the blood pump and the quality of life of the patients with heart failure. Therefore, reducing the number of power devices is one of the effective ways to reduce power consumption and volume.
Compared with the traditional drive mode of combining magnetic bearings and motor, the axial flux bearingless motor can produce a suspension force and rotating torque at the same time with only one three-phase winding, which greatly simplifies the structure and the control system. In addition, the axial flux motor has the characteristics of a small volume, light weight and high efficiency, when compared to the with radial flux motor. Thus, it has greater advantages in the miniaturization and reduced weight of the blood pump. For example, Toru [9] designed a suspension system with a double offset hybrid magnetic bearing. The upward magnetic pull generated by the magnetic bearing balanced the downward bias magnetic pull of the bearingless motor, and the impeller could be suspended in the axial direction center position of the blood pump through active regulation. The axial and tilt motion of the impeller were actively regulated by the electromagnetic bearing, while the radial motion was constrained by the passive stability. Nevertheless, the magnetic suspension system of the blood pump has problems, such as a large volume and complex structure. If the permanent magnet bearings are used to replace the electromagnetic bearing structure, and the only driving component is the bearingless motor, the number of driving components and the blood pump volume can be greatly reduced. In addition, the permanent magnets of the bearingless motor can adopt the magnetization mode of the Halbach array, which can improve the sinusoidality of the air gap flux and enhance the magnetic density. This kind of motor with a special magnetization direction helps to improve the motor efficiency and reduce the torque ripple on the basis of increasing the average torque [10,11,12]. However, the negative stiffness of the magnetic suspension system makes the impeller unable to suspend stably. The repulsive hydrodynamic bearing can eliminate the negative stiffness of the bearingless motor without adding hardware, ensure the rotor suspension stability, and helps to reduce the vibration of the rotor via damping.
Hydraulic bearings can be categorized into hydrostatic bearings and hydrodynamic bearings. The hydraulic suspension structure design and suspension stability analysis has been extensive studied [13,14]. For example, Wang et al. [15] designed a hydraulic bearing for the radial suspension support of an axial flow blood pump; established a theoretical model for the hydraulic bearing force; studied the influence law of the rotating speed and radial eccentricity on the bearing force through numerical simulation, and measured the magnetic–hydrodynamic, double-levitated impeller stability under different external circulation experiments via the Hall sensor and laser vibration meter. Similarly, Chung et al. [16] studied the rotor dynamic characteristics of the hydrodynamic bearing of the rotary blood pump at different flow rates and impeller speeds; obtained the stiffness coefficient and damping coefficient of the system by applying an external interference, and analyzed the influence of the law of speed and flowrate on the stiffness coefficient and damping coefficient. In order to improve the suspension stability of the blood pump and the hemolysis performance, Kosaka et al. optimized the number of spiral grooves, groove depth and bearing clearance of the hydrodynamic bearing structure and obtained the best parameter combination [17,18]. Kataoka et al. studied the rotor radial movement of the blood pump coupled with a hydrodynamic bearing and permanent magnetic bearing [19]. From the experimental results, it was found that the radial direction motion of the rotor had an unstable circular mode. The occurrence of the unstable rotation mode could be avoided by increasing the speed and adopting a symmetrical structure.
In conclusion, the hybrid suspension system of the axial flux, single-winding, bearingless motor, permanent magnet bearing, and hydrodynamic bearing solves the problems of a complex structure, heavy weight and high energy consumption. The structural design of the hydrodynamic bearing is closely related to the stiffness characteristics of the magnetic suspension system composed of a bearingless motor and permanent magnet bearing. However, there are few studies on the structural design and stability analysis of the hydrodynamic bearing based on the stiffness characteristics of a bearingless motor. Therefore, the present study was envisaged to decipher the axial force and stiffness characteristics of the suspension system by using a combination of theoretical calculations and simulation solutions. This study also aimed to design the hydrodynamic bearing structure according to the relationship between the stiffness and structural parameters and verify the suspension stability of the blood pump through experimental validation. Our study can lay the foundation for the active regulation of the rotor pose of the blood pump and the optimization of blood compatibility.

2. Magnetic Suspension System Modeling

2.1. Structure and Working Principle

The structure composition of the magnetic–hydrodynamic hybrid suspension centrifugal blood pump designed in this study is shown in Figure 1a. The blood pump shell was composed of the inlet, upper casing, lower casing, and motor casing, as shown in Figure 1b. The blood pump impeller structure was in the form of a closed straight blade cover plate. The upper end face was embedded with a concentric, movable, permanent magnetic ring, which was matched with the concentric, fixed magnetic ring of the upper casing to form a permanent magnetic bearing. The permanent magnets with different polarities were embedded in the lower end face of the impeller to cooperate with the motor stator, in order to realize the rotating motion of the impeller. Therefore, the impeller can also be regarded as the motor rotor. The impeller adopted the support scheme of the combination of passive and active suspension, and the rotating motion of the impeller was consistent with the driving mode of the ordinary, permanent, magnet-synchronous motor. The stator coil and iron core of the motor were placed in the motor casing to form a single-winding bearingless motor with the permanent magnets on the impeller. The hydrodynamic bearing was located on the upper and lower pump casing and exerted an axial fluid force on the impeller under the fluid action. Therefore, the impeller can be suspended in the blood pump cavity under the coupling action of the permanent magnet bearing, the motor electromagnetic force and the hydrodynamic bearing force.
As shown in Figure 2a, the structural design scheme of a single-winding bearingless motor consisted of a rotor backplate, magnetic pole, coils, iron core and stator yoke. The motor of the blood pump was an axial flux, bearingless, permanent, magnet-synchronous motor, which realized the rotation of the impeller in the form of single-winding structure, and cooperated with the permanent magnet bearing to ensure the stable suspension of the impeller in the pump casing. The structure composition and size of the bearingless motor stator are shown in Figure 2b. The Halbach array structure was adopted for a magnetic pole arrangement to improve the blood pump efficiency and reduce motor power consumption. Its structure is shown in Figure 2c.

2.2. Modeling of the Bearingless Motor

The magnetic poles of the motor rotor were fixed on the rotor backplate. Therefore, the motor was a non-salient pole structure, and the magnetic poles were magnetized according to the Halbach array arrangement. The specific composition of the magnetic suspension system and permanent magnet bearing are shown in Figure 3a,b.
The two concentric, permanent magnet rings embedded in the upper-end face of the impeller and the two concentric, permanent magnet rings embedded in the upper casing cooperated with each other, and the permanent magnet bearing generated a magnetic pull force in the z axial direction. Therefore, the component of the magnetic force generated by the motor stator and the permanent magnet rings in the z axial direction jointly determined the axial position of the impeller. The x, y axial translation and ϕx, ϕy rotation of the impeller were passively constrained by the magnetic coupling between the permanent magnet bearing and the bearingless motor. Therefore, the translation and rotation of the impeller in x and y axis direction were limited. Let the upper and lower air gaps of the rotor be represented as g1 and g0, respectively. When the rotor is translated z in the positive axial direction due to disturbance, the upper and lower air gaps becomes g1z and g0 + z. Therefore, in order to maintain the stable suspension and rotation of the impeller at the balance position, it is crucial to regulate the motor stator current. The specifications of the bearingless motor and permanent magnet bearing are shown in Table 1.
The bearingless motor stator adopted the structure of the three-phase concentrated windings. The three-phase windings, A, B and C, were arranged in a space at an interval of 120°, and the number of turns was Ns. The rotor and stator structure of the bearingless motor, as shown in Figure 4a, is used to analyze the rotation motion of the motor, expand the spatial structure of the motor along the circumferential direction and establish the coordinate system, as shown in Figure 4b. The direct axis and the quadrature axis are synchronous, dynamic coordinate systems fixed on the motor rotor, which are used to describe the rotating motion of the motor. The magnetic field direction of the motor is the d axis and the axis perpendicular to the magnetic field direction is the q axis. Therefore, the spatial phase difference between the d axis and the q axis is π/2. Permanent magnets with different magnetic pole directions are pasted on the lower surface of the rotor of the bearingless motor, and the number p of rotor pole pairs was four. The three-phase coils of the stator were fed with a sinusoidal current to drive the motor. The magnetic flux distribution ψ formed by the magnetic field interaction of the permanent magnet and the motor energizing coil is shown by the black dotted arrow in Figure 4b. In order to study the rotor dynamic characteristics of blood pump, the magnetic field of the bearingless motor was theoretically analyzed using the inductance matrix method. The air gap magnetic energy was calculated by the magnetic flux distribution and inductance matrix, and the torque and axial electromagnetic force were obtained by deriving the rotation direction θ and z-axis direction. In order to calculate the inductance and deduce the electromagnetic force and torque model through the effective permeability, the magnetic field generated by the permanent magnet on the rotor was equivalent to a direct current excitation winding with the turns Nf and current if.
The electromagnetic force and the torque of the bearingless motor were theoretically modeled, and the assumptions were as follows:
(a)
The spatial distribution of the magnetomotive force was sinusoidal.
(b)
The stator surface was smooth and the slot harmonics could be ignored.
(c)
Magnetic saturation was ignored.
(d)
The permeability of the stator core was infinite, and the permeability of the permanent magnet and air was considered equal to that of vacuum.
(e)
The iron loss was ignored.
For the non-salient rotor, the self-inductance of the stator coil was independent of the rotor position angle θ. Therefore, the self-inductance of the stator coil was constant, and the d axis inductance and q axis inductance were equal. The self-inductances of the three-phase stator coils are given by Equation (1):
L a a = L b b = L c c   = N s ψ i a + L a l   = L a a 0   + L a l  
Therefore, the self-inductance of the stator was composed of Laa0 and Lal, where Laa0 is the self-inductance component generated in the coil by the fundamental flux in the air gap; and Lal is the additional component generated for the leakage flux of armature winding, which is also called the leakage inductance Lal. The self-inductance of a winding is generally a function of the air gap g [20]. Therefore, the self-inductance of stator winding can be approximately expressed by Equation (2) with air gap g:
L a a 0   = L a a 0 g  
where L′aa0 is the stator magnetizing inductance multiplied by the air gap length, and the unit is Hm. L′aa0 and Lal are constants, which are obtained by the motor parameter identification and the phase inductance measurement. The mutual-inductance between the two windings is the product of the self-inductance and cosα, whose phase angle difference is the electrical angle α. Since the angle between the different phases of armature winding was 120°, the mutual inductances between the different phases of armature winding are equal, which is expressed by Equation (3):
L a b = L b a   = L a c = 1 2 L a a 0   = 1 2 L a a 0 g  
Similar to the self-inductance of stator, the self-inductance Lf of the rotor excitation winding is expressed in Equation (4):
L f f   = L f   = L f f 0   + L f l = L f g + L f l  
where Lff0 is the inductance generated by the fundamental component of magnetic flux in the air gap. The leakage inductance Lfl of the rotor excitation winding is used to characterize the influence of permanent magnet leakage flux. The mutual inductance between the stator and rotor is related to the rotor position and changes periodically with θ, so the mutual inductance between the excitation winding f and phase winding a is given by Eqaution (5):
L a f = L f a   = L m c o s θ = L m g c o s θ  
where L′m is the constant part of the mutual inductance of the motor stator and rotor. It was observed that the relationship between the self-inductance and the mutual inductance of the stator and rotor that the relationship   L m = L a a 0 L f was satisfied. When the windings of the motor stator were three-phase symmetrical, ia + ib + ic = 0, and thus ib + ic = −ia. On substituting the self-inductance of the above phase A and the mutual inductance of the excitation winding f, the phases B and C into the flux linkage expression were obtained as shown in Equation (6):
λ a = ( 3 2 L a a 0 + L a l ) i a + L a f i f  
If the synchronous inductance Ls = 3/2Laa0 + Lal is defined, the corresponding flux linkage equation of phase A is:
λ a = L s i a + L a f i f  
Therefore, the phase voltage of the motor is composed of the voltage on the phase resistance and the back electromotive force, and its expression is:
{ u d = R s i d + L s d i d d t ω L s i q u q = R s i q + L s d i q d t + ω L s i d + ω λ m  
The flux linkage was obtained by the product of the inductance matrix of the motor and the current, as shown in Equation (9):
λ = ( λ f λ a λ b λ c ) T = [ L f f L a f L b f L c f L a f L a a L a b L c a L b f L a b L b b L b c L c f L c a L b c L c c ] [ i f i a i b i c ]
where L is the inductance matrix of the motor; the conversion matrix between the 3-phase current and d axis and q axis currents is given by Equation (10):
C = 2 3 [ 3 2 0 0 0 c o s θ s i n θ 0 c o s ( θ 2 π 3 ) s i n ( θ 2 π 3 ) 0 c o s ( θ + 2 π 3 ) s i n ( θ + 2 π 3 ) ]
Then, the motor flux linkage equation was converted into:
[ λ f λ d λ q ] = [ L f 3 2 L m 0 L m L s 0 0 0 L s ] [ i f i d i q ]  
where, λm = Lmif is the flux linkage generated by the magnetic field of the rotor permanent magnet in the stator winding. The magnetic energy W stored in the coil was calculated by Equation (12):
W = 1 2 ( i f i d i q ) ( λ f λ d λ q ) T
If the motor rotor deviated slightly from the balance position g0 in the axial direction, the magnetic energy changed accordingly. Therefore, the axial magnetic force and torque were calculated by the partial derivative of the magnetic energy in the axial displacement g and rotation direction θ [21]. The axial electromagnetic force between the motor stator and rotor is:
F m 0 = W g = 1 2 g 2 [ L f i f 2 + 5 2 L m i d i f + 3 2 L a a 0 ( i d 2 + i q 2 ) ]  
Additionally, the torque Ts was expressed as:
T s = 3 2 p ( λ d i q λ q i d ) = 3 p L m i f i q 2 g  

2.3. Axial Force Modeling of the Permanent Magnet Bearing

The magnetic field of permanent magnetic bearing was solved by the static magnetic field. The permanent magnetic force was obtained by integrating the Maxwell stress tensor on the surface of a permanent magnetic bearing [22]. The calculation formula used was:
f p = n · T = 1 2 n ( H · B ) + ( n · H ) B T  
where n is the outward normal vector of the outer surface of the permanent magnet bearing and T is the Maxwell stress tensor on the outer surface of the permanent magnet bearing. The axial component fz of the permanent magnetic force fp is the axial force of the permanent magnetic bearing on the rotor, and the direction is vertical upward along the z axial direction. Its derivative relative to the axial displacement is the magnetic stiffness [23], which was defined as:
k z 1 = d f z d z = f ( z k ) f ( z k + 1 ) z k + 1 z k  
The three-dimensional model of the permanent magnet bearing assembly was established, and the schematic diagram is shown in Figure 5. The upper concentric magnetic rings were arranged on the blood pump casing, which were static magnetic rings. The lower concentric magnetic rings were arranged on the upper-end face of the impeller and were movable magnetic rings.
The material of the permanent magnetic bearing was the permanent, magnetic, rare earth material NdFeB35, which adopted an axial magnetization. The specific parameters of the permanent magnetic bearing are shown in Table 2.
The magnetic field of the permanent magnet bearing was solved by the static field module in the COMSOL software, and the magnetic force was obtained by integrating the Maxwell stress surface. When simulating the static magnetic field of the permanent magnet bearing, the axial air gap scanning range was set to 0.1~3.1 mm. In the permanent magnet bearing, the movable magnetic ring was offset under the action of the external load and interference. Figure 6a shows the magnetic density and magnetic induction line distribution when the air gap g1 was 0.4 mm and 1.72 mm (Δz = 1.32 mm). The air gap was inversely proportional to the magnetic density and the strength of the magnetic field. The variation law of the axial restoring force Fm1 is shown in Figure 6b. There was a strong nonlinearity between the axial restoring force Fm1 and the air gap g1. The smaller the air gap between the movable magnetic ring and the static magnetic ring, the greater the axial magnetic force. The equilibrium position of the rotor was g1 = 1.6 mm and the movable range was 1.5–1.7 mm. Within this range, the axial restoring force Fm1 and axial offset were approximately linear. Therefore, the axial force was linearized near the equilibrium position to obtain the linearization model of the permanent magnet bearing.
The axial force at the different axial displacements was fitted, and the mathematical model between the axial force Fm1 and the axial displacement g1 was established, as shown in Equation (17). The root mean square error RMSE of the fitting results was 0.042 and the correlation coefficient R was 0.999:
F m 1 = 2.1 g 1 3 46.88 g 1 2 + 135.32 g 1 1.5 127.65 g 1 + 37.72  

2.4. Model Linearization of Magnetic Suspension System

Without considering the effects of fluid force, gravity and external disturbance, the total axial force on the rotor is given by Equation (18):
F z = F m 0   + F m 1  
When the motor rotor is in a balanced and stable position, the air gap between the stator and the rotor surface is g0. Assuming that the rotor has an axial offset z and the rotor deviates from the blood pump inlet along the axial direction, the air gap becomes g0 z. The axial force and electromagnetic torque applied to the rotor are linearized at the equilibrium point, and Equations (13) and (14) are expanded by the Lagrange linearization method, then the respective linearization models of axial force and torque are:
F z   = 1 2 g 0 2 [ L f i f 2 + 5 2 L m i d i f + 3 2 L a a 0 ( i d 2 + i q 2 ) 9.92 g 0 2 ] z 2 g 0 3 [ L f i f 2 + 5 2 L m i d i f + 3 2 L a a 0 ( i d 2 + i q 2 ) 9.6 g 0 3 ]
T l = 3 p L m i f i q 2 g 0 ( 1 z g 0 )
The air gap g0 of the motor was 1.8 mm. Due to structural constraints, the maximum displacement on one side of the rotor was 0.1 mm and the axial displacement z was relatively small during the actual axial position control of the rotor. Therefore, z/g0 was approximately 0. In addition, kt = 3pL′mif/2g0 instead of the constant term in Equation (20), Equation (20) is simplified as:
T l = k t i q  
According to Equation (19), the axial force Fz was coupled by the direct axis current id and the quadrature axis current iq. If the coupling term is ignored, the axial force is only related to id and can be controlled linearly. According to Equation (21), the electromagnetic torque is only related to the quadrature axis current iq. After linearization, the axial force shown in Equation (19) can be expressed as in the form of z and id as:
F z = K i i d + K z z + f  
where, K i = 5 L m i f / 4 g 0 2 is the force gain; K z = [ L f i f 2 + 2.5 L m i f i d + 1.5 L a a 0 ( i d 2 + i q 2 ) 9.6 g 0 3 ] / 2 g 0 3 is the axial displacement stiffness; other terms of the model can be regarded as disturbance terms, that is, f = [ L f i f 2 + 1.5 L a a 0 ( i d 2 + i q 2 )   9.92 g 0 2 ] / 2 g 0 2 . When the impeller was passively suspended in the equilibrium position, f is 0 N. Without considering gravity and fluid force, according to Newton’s second law of motion, the axial motion of rotor can be expressed as:
m z ¨ K z z = K i i d + f  
According to Equation (23), it can be seen that the rotor of the bearingless motor exhibited a negative stiffness in the axial direction, which is unstable for the rotor suspension state and seriously affects the stable operation of the blood pump.

2.5. Parameter Identification of the Magnetic Suspension System

In order to establish an accurate motor suspension force and torque model, the motor parameters were obtained by parameter identification, and the parameters in the suspension force and torque model were calculated. Figure 7 shows the bearingless motor parameter identification platform, which adjusted the air gap between the rotor and stator by changing the position of the sliding table. The sliding table displacement was measured by the spectral confocal displacement sensor: the sensor probe was fixed on the fine-tuning platform; the probe was fixed through the fine-tuning platform and adjusted to the measurement range. Subsequently, the change in the air gap was calculated by measuring the sliding distance of the sliding table in real time. The displacement sensor mode was IFC2461(Micro-Epsilon Co. Ltd., Ortenburg, Germany), the measurement frequency was 25 kHz, the absolute error was less than 0.05% and the resolution was 0.004%, which helped to realize the micron precision measurement of displacement. The motor, displacement sensor and motor controller were powered by a DC-regulated power supply. The power supply model was ITECH6322 and the voltage accuracy was 1 mV. The voltage of the motor and displacement sensor was 24 V and the power supply voltage of motor controller was 5 V.
The position of the sliding table was adjusted, so that the air gap between the motor stator and rotor was 1.8 mm, and the rotor position was fixed. The power supply was turned on and the voltage was adjusted. The motor parameters were identified through the motor controller and multiple experiments were conducted. The motor parameter identification results are shown in Table 3.
Through the motor parameters, the axial force and torque of the motor was calculated, and the model parameters are L′aa0 = 8 × 10−7 Hm, the leakage inductance Lal = 5.3 × 10−5 H; L′mif = 4.88 × 10−6 Ahm and the torque parameter kt = 1.63 × 10−2 N·m/A.

2.6. Verification of the Magnetic Suspension System Theoretical Model

The transient electromagnetic field of the bearingless motor was solved using the ANSYS Electronics software. The rotor backplate was made of titanium alloy. The permanent magnet material was NdFeB35, the residual magnetic flux density Br was set to 1.3 T, and the different magnetization directions were set according to the magnetic pole orientation. The material of the stator winding was copper and the number of turns was 115. The material of the iron core and yoke was silicon steel. The number of turns and sinusoidal current excitation were set for the coil. The AC current excitation amplitude im was 0.4 A and the spatial phase difference was 120 degrees. The electromagnetic field solution time was set to 10 ms, and the simulation solution time step was the time corresponding to each 1° rotation of the motor.
The axial force received by the rotor of the bearingless motor during rotation is shown in Figure 8a. The axial force obtained by the FEM simulation fluctuated all the time during the rotation process, and the amplitude was essentially within 0.5 N. The axial force was negative, indicating that the motor rotor was always subject to the suction of the stator. The average axial force was calculated by FEM to be −5.2 N. Assuming that the rotor had no axial offset, the direct axis current id was 0 A and the quadrature axis current iq was 0.4 A; the axial force calculated by the theoretical model was −4.63 N. The two results were close, and the error of the theoretical model was 11%. The variation process of the electromagnetic torque during the motor rotation is shown in Figure 8b. There is an obvious pulsation phenomenon of the rotor torque during rotation. The potential causes of torque ripple were the interactions of the stator and rotor magnetic field, current commutation, the cogging effect, and the armature reaction. The average value of the electromagnetic torque in the FEM simulation results was 6.6 N·mm. According to the theoretical Equation (21), the theoretical torque at 3000 rpm was 6.5 N·mm, and the error was 1.5%. The two values were very close, which verified the accuracy of the motor torque theoretical model.
When the motor rotor was offset under the action of external disturbance, the changes in the permanent magnetic force on the rotor are shown in Figure 9a. Because the motor structure and the permanent magnet bearings were spatially symmetrical, the permanent magnetic force in the x and y direction was always 0 N, when the impeller was offset, and the axial force Fz changed with the offset. When the impeller was offset to the inlet of the blood pump, the axial force on the impeller was upward, that is, the axial force direction was consistent with the offset direction, and its amplitude increased with an increase in the offset. Therefore, the impeller was unstable in the axial direction. The resultant axial force at the balance point was close to 0 N.
Figure 9b shows the theoretical value and the FEM simulation results of the axial force of the rotor within the axial movable range. It can be seen from Figure 9b that the theoretical value of axial force was very close to the FEM simulation results, which verifies the accuracy of the theoretical model. Thus, the stiffness calculated by the theoretical model can be used as a reference index for the structural design of the hydrodynamic bearing. The position where the resultant force was 0 N was located below the theoretical equilibrium point. Furthermore, the linear relationship between the axial force and the axial displacement at the equilibrium point also verified the rationality of the linearization of the axial force theoretical model. According to the definition of stiffness, the axial stiffness was negative, which indicated that the impeller was unstable under the action of external disturbance.

2.7. Hydrodynamic Force and Stiffness Model of Hydrodynamic Bearing

The hydrodynamic bearing was located on the inner cavity surface of the upper and lower casings of the blood pump. The fluid force generated by the dynamic pressure effect acted on the rotor to realize the axial stable suspension. Due to the secondary flow field at the upper and lower clearance of the blood pump, the blood flow field was in a spiral form. Although the herringbone groove hydrodynamic bearing had a better stability, the groove type was inconsistent with the blood flow form. When the blood flowed through the hydrodynamic bearing, flow separation and turbulence occurred, which could easily damage the blood. Therefore, in order to obtain a better blood compatibility, the hydrodynamic bearing was designed in the form of a spiral groove. The specific structure is shown in Figure 10a. The groove type of the hydrodynamic bearing is spiral. According to the rotation direction of the impeller and the flow direction of blood, the hydrodynamic bearing on the lower pump casing was designed as the “Pump-in” type, and the hydrodynamic bearing on the upper casing was designed as the “Pump-out” type. The structural parameters of the two hydrodynamic bearings were consistent. The inlet of the hydrodynamic bearing was located outside the spiral groove with radius rs. When the impeller rotated in the direction shown in Figure 10a, the high-pressure fluid emitted by the impeller flowed into the hydrodynamic bearing from the inlet of the spiral groove. As the groove width of the spiral groove decreased in the flow direction, the wedge-shaped spiral groove formed a dynamic pressure effect, resulting in the bearing capacity to realize the suspension of the impeller. Figure 10b shows the cross-sectional view of the hydrodynamic bearing. The vertical distance between the lower surface of the impeller and the pump casing is the liquid film thickness h1 of the hydrodynamic bearing, that is, the suspension clearance of the impeller. The groove depth of the spiral groove of the hydrodynamic bearing is h2, and the sum of h1 and h2 is the distance h between the bottom surface of the impeller and the groove depth. The liquid film thickness h1 did not only affect the suspension performance of the hydrodynamic bearing, but also affected the hemolysis performance. Therefore, the liquid film thickness was designed after comprehensively considering the performance of the hemolysis and suspension. Here, this factor is not considered temporarily in the structural design of the hydrodynamic bearing, and only the influence of the groove depth h2 on the suspension performance was analyzed. In addition to the groove depth, the groove number N, the groove width ratio δ = Bg/B and spiral angle β were the key parameters affecting the hydrodynamic bearing performance.
The spiral groove of the hydrodynamic bearing was in the form of a logarithmic spiral, and its structural shape was closely related to the spiral angle β, which is expressed as:
r = r g e θ t a n β  
where, rg is the base circle radius of spiral line, θ is the spiral line control angle, and β is the spiral angle. If the pressure change in the clearance direction was ignored, the liquid film pressure distribution of the hydrodynamic bearing could be solved based on the Reynolds equation in lubrication theory. The steady-state Reynolds equation is shown in Equation (25):
t · ( ρ h 1 v a v e ) = 0  
where, ρ is the fluid density, h1 is the clearance, h1 = hc + hb, hc = −uc·nref, hb = hb1 + ub·nrefub·thb1, uc is the impeller displacement, nref is the reference normal direction, hb1 is the initial clearance, vave is the average speed, and the calculation formula is:
v a v e = 1 2 ( v c , t + v b , t ) h 1 12 μ t p f  
where μ is blood viscosity; vc,t and vb,t are the velocity vectors of the impeller surface and hydrodynamic bearing wall at the clearance, respectively. The specific calculation method is shown in Equations (27) and (28):
v c , t = ( I n r e f n r e f ) { u ˙ c . + Ω × ( X X c ) }
v b , t = ( 1 n r e f n r e f ) u ˙ b
The absolute pressure pA at any point on the bearing surface was calculated from the sum of the reference pressure pref and the relative pressure pf, which is pA = pref + pf. The pressure distribution on the whole hydrodynamic bearing surface was obtained through an iterative calculation. The pressure of each node on the hydrodynamic bearing surface was double integrated along the radial direction and circumferential direction to obtain the bearing capacity Ff of the liquid film:
F f = r i r o p A r d r d θ
where, ri and ro are the inner radius and outer radius of the hydrodynamic bearing, respectively. When the axial displacement of the impeller occurred, the clearance between the impeller and the hydrodynamic bearing changed correspondingly, as did the bearing capacity imposed by the hydrodynamic bearing on the impeller. The change rate is the stiffness of the hydrodynamic bearing, which represents the repulsion of the hydrodynamic suspension bearing to the disturbance, which is defined as:
K f = d F f d Δ h 1  

2.8. Axial Suspension Stability Analysis of Magnetic–Hydrodynamic Hybrid Suspension System

In order to eliminate the negative stiffness of the magnetic suspension system and realize the stable suspension in the axial direction, the hydrodynamic bearing is required to provide stable support for the rotor, eliminate the negative stiffness and ensure that the rotor has a certain passive suspension stability under the action of external disturbance. The stiffness and damping characteristics of the hydrodynamic bearing can effectively limit the axial movement and amplitude of the rotor vibration. Therefore, based on the dynamic model of the magnetic suspension system, it was also necessary to include the model parameters related to the action of the hydrodynamic bearing and the flow field, in order to establish the dynamic model of the rotor under the action of the magnetic field, flow field, gravity and external disturbance. The hydrodynamic performance was optimized through passive suspension stability analysis, in order to improve the suspension stability of the blood pump and reduce the power consumption of the blood pump motor. The axial dynamic model of the rotor is shown in Equation (31):
m z ¨ + C f z   ˙ + ( K f K z ) z = 0
where m is the rotor mass, Cf is the damping of the hydrodynamic bearing and Kf is the stiffness of the hydrodynamic bearing. The stiffness Kf and damping Cf were similar to the active regulation of impeller spatial position through the PD controller, and their values determined the stability of the suspension system. According to Equation (31), the damping coefficient Cf > 0; therefore, the prerequisite for the system stability was KfKz > 0; the axial stiffness amplitude of the hydrodynamic bearing was greater than that of the magnetic suspension system, and the greater the difference between the two stiffnesses, the more stable the system. For the passive suspension stability, the motor current id and iq were both 0 A. Therefore, according to the calculation formula of kz, the axial stiffness kz of magnetic suspension system was −2547 N/m. Since the hydrodynamic bearing is symmetrically arranged on the inner cavity surface of the upper and lower pump casings, when analyzing the suspension stability, the suspension stability was guaranteed only if the stiffness of the one-sided hydrodynamic bearing was greater than −1274 N/m.

3. Stiffness Analysis and Structural Design of the Hydrodynamic Bearing

The parameters of the hydrodynamic bearing mainly included the groove number N, spiral angle β, groove width ratio δ and groove depth h2. The effects of different structural parameters on the stiffness and hydrodynamic force were studied to obtain the best combination of the structural parameters. The inlet and outlet pressures of the hydrodynamic bearing were 82 mmHg and 48 mmHg, respectively, and the impeller speed was 2600 rpm. The blood viscosity was 3.5 mPa·s and the density was 1055 kg·m3. Figure 11 shows the solution results of the pressure distribution on the surface of the hydrodynamic bearing. The pressure distribution was uniform and the blood pressure in the spiral groove was relatively high. When the blood from the impeller flowed through the spiral groove of the hydrodynamic bearing, its pressure decreased continuously, and the bearing capacity was formed through the dynamic pressure effect in the wedge gap to realize the suspension of the impeller.
The radial pressure distribution of the hydrodynamic bearing is shown in Figure 12a. When the blood flowed from the pad edge of the hydrodynamic bearing into the spiral groove, the blood pressure increased continuously, the blood maintained a large pressure in the spiral groove, and the dynamic pressure effect was strong. When the blood flowed out of the spiral groove to the edge of the pad, the blood pressure decreased continuously. It can be seen from Figure 12b that when the groove width ratio δ was 0.5, the blood pressure first increased and later decreased along the rotation direction in the circumferential direction. However, the blood pressure at the spiral groove was higher than that in the pad area, indicating that there was an obvious dynamic pressure effect of the spiral groove, which helped to provide a large bearing capacity for the impeller. In addition, according to the pressure distribution diagram of the hydrodynamic bearing, the pressure clearly changed at the junction of the spiral groove and the pad.
It can be seen from Figure 13a that the greater the number of grooves, the greater the dynamic pressure effect on the spiral groove. Accordingly, the hydrodynamic force of the hydrodynamic bearing was higher, and its range was from 7.25–7.36 N. The stiffness Kf also increased with the increase in the number of grooves, and the stiffness range was from 2300–3900 N/m, greater than the critical stiffness value of −1274 N/m. Therefore, in order to achieve a stable suspension, a higher number of grooves in the hydrodynamic bearing offered a better suspension performance. The influence of the groove depth on the stiffness and hydrodynamic force is shown in Figure 13b. The hydrodynamic force increased with an increase in the groove depth h2. The stiffness first increased and then decreased with an increase in the groove depth h2, and the stiffness range was 2450–4000 N/m. It can be seen from the stiffness curve that, in order to obtain a high axial stiffness Kf, the range of the groove depth h2 was around 0.06 mm. If it was considered to obtain high stiffness Kf, while ensuring high axial hydrodynamic force, the groove depth was set to about 0.1 mm.
The influence of the groove width ratio δ on the stiffness and hydrodynamic force is shown in Figure 14a. When the groove width ratio is 0 and 1, the hydrodynamic bearing becomes a full pad and groove, and the dynamic pressure effect disappears. Therefore, a groove width ratio in the range from 0–1 will maximize the hydrodynamic force. The selection range of groove width ratio was from 0.3–0.9. With the increase in the groove width ratio, the hydrodynamic force initially increased and later a decrease was observed. The groove width ratio with a peak hydrodynamic force was near 0.9, but the stiffness decreased with an increase in the groove width ratio, and its value decreased from 3500 N/m to 2060 N/m. Figure 14b shows the variation law of the stiffness and hydrodynamic force with spiral angle. With a gradual increase in the spiral angle, the hydrodynamic force initially increased and then gradually decreased. The stiffness first increased followed by a decrease accompanied by an increase in the spiral angle, and its range decreased from 6000 N/m to 5600 N/m. A high hydrodynamic force and axial stiffness was achieved when the spiral angle β was 65°.
From the above analysis, it can be seen that under different structural parameters, the hydrodynamic force was essentially between 7.17 and 7.43 N with little changes, while the stiffness range corresponding to different structural parameters of hydrodynamic bearing changed significantly. According to the variation law of each structural parameter pair, in order to obtain a relatively high axis stiffness, the number of grooves and the spiral angle should be as large as possible, and the groove depth and groove width ratio should be as small as possible. In order to ensure the stable suspension of the blood pump impeller under an external disturbance, the axial stiffness should be much higher than the axial negative stiffness of the bearingless motor. According to the hydraulic characteristics and the flow field distribution of the blood pump, based on the above structural parameters and stiffness analysis of the hydrodynamic bearing and considering the constraints of the structural dimensions of the blood pump, the structural parameters of the hydrodynamic bearing were designed as follows: the number of grooves was 20, the spiral angle was 45°, the groove width ratio was 0.5 and the groove depth was 0.1 mm. The titanium alloy casing machined by CNC (Computerized Numerical Control) with the hydrodynamic bearing is shown in Figure 15. After calculation, the axial stiffness at 1800 rpm was found to be 3819 N/m, which was much greater than the critical negative stiffness −1274 N/m of the magnetic suspension system. Therefore, the blood pump could also ensure the ability of stable suspension at high speeds.

4. Experimental Validation of the Axial Passive Suspension Stability of the Blood Pump

In order to simulate the interaction between the blood pump and the cardiovascular system and study the suspension stability of the blood pump during the auxiliary blood pumping, an experimental platform, used to imitate the circulatory loop, was built, as shown in Figure 16. The system consisted of a left atrial simulator, left ventricular simulator, aortic simulator, venous circuit pipeline, mitral valve, aortic valve and damping valve. The measurement module included a pressure sensor (PX409-015CGUSBH, ±1 bar, 0.06% accuracy, Omega Engineering Inc., Norwalk, CA, USA) and a displacement sensor IFC2461. The control module was composed of the left ventricular simulator controller and the blood pump controller. The measured pressure signal was transmitted to the computer through the communication interface. By controlling the motion of the left ventricular simulator, the contraction and relaxation of the ventricle were simulated to produce a pulsatile physiological blood flow. The blood pump controller controlled the rotation of the impeller according to the set target speed to realize the physiological blood perfusion. The mock circulatory loop could simulate the response of the blood pump to cardiac unloading and cardiovascular system at a constant speed and was used to study the influence of the time-varying disturbance generated by the cardiovascular system on the suspension state of the blood pump impeller.
The axial relative displacement of the impeller was measured by the spectral confocal displacement sensor, and the signal with high-frequency noise was filtered. The impeller position was set to 0 mm, when the blood pump was not started, and the measured value of the sensor truly reflected the change in the axial suspension position of the impeller. It can be seen from Figure 17a that the impeller of the blood pump could realize the axial suspension at a low speed of 1800 rpm; the higher the speed, the greater the axial displacement.
Figure 17b shows the id and iq current waveforms corresponding to different rotating speeds of the blood pump. Since the study focused on the passive suspension stability of the blood pump, the id was always maintained at 0 A through a closed-loop control, and the iq current increased with an increase in the rotating speed. When no active control was applied in the axial direction, the higher the speed, the more clear the dynamic pressure effect, and the greater the hydrodynamic force of the bearing. In addition, with the increase in the impeller speed, the flowrate of the blood pump and the impact force of the fluid on the impeller also increased. The above reasons made the axial displacement of the impeller larger, that is, the suspension clearance became larger. It was observed that the different speed changes exhibited little effect on the axial displacement, and the change range was 0.2 μm.
Figure 18 shows the hydraulic efficiency curve of the blood pump at different speeds and flowrates. It can be seen from the hydraulic efficiency curve that the speed and flowrate had a great impact on the hydraulic efficiency of the blood pump. The smaller the flow, the higher the speed, and the greater the hydraulic efficiency of the blood pump. When the blood pump impeller rotates at high speed, the dynamic pressure effect is more clear, the suspension stability is better and the hydraulic efficiency is higher.
In the experimental validation of the mock circulatory loop works, the influence of the blood pumped by the failed heart on the suspension rotor of the blood pump was simulated and the output blood pressure waveform is shown in Figure 19a. The peak blood pressure is 70 mmHg, which is lower than the normal human blood pressure. At this time, the impeller speed was set at different values to study the change law of the axial suspension position when the impeller speed changed under specific disturbance. It can be seen from Figure 19b that the impeller has achieved a stable suspension under the action of external disturbance. Under the same external disturbance, the axial displacement of the impeller changes greatly compared to that without the external disturbance, indicating that the external fluid impact can lead to the displacement of the impeller suspension position, and the impeller could move axially with a certain frequency with the action of the external disturbance. In addition, the change of impeller speed could also affect the impeller suspension position. Thus, it was found that an increase in the speed led to a higher suspension clearance and a better suspension stability.
Figure 20a shows the disturbance waveform applied to the blood pump by the extracorporeal mock circulatory loop. The pulsation frequency of Case 1 to Case 4 continues to decrease, but the pulsation amplitude of the blood pressure shows a persistent increase. When the blood pump impeller was run at a constant speed of 2600 rpm, the extracorporeal mock circulatory loop gave an output for the four blood pressure waveforms shown above to simulate the suspension stability of the blood pump under the action of different frequencies and the blood pressure waveforms pumped by the failed heart. Under the external fluid disturbance, the blood pump impeller floated up and down in the axial direction, and its frequency was consistent with the external disturbance frequency. The greater the external fluid pressure, the greater the inlet pressure of the hydrodynamic bearing, and the more obvious the dynamic pressure effect. It was found that a larger the hydrodynamic force increased the suspension clearance of the blood pump. The results are shown in Figure 20b. The results show that the fluid pressure significantly affected the axial suspension position of the impeller. Compared to the axial clearance, the axial displacement is small. Before the blood pump starts, the magnetic suspension system is unstable, and the impeller is generally adsorbed on the inlet side casing. When the rotating speed of the blood pump increases, the dynamic pressure effect increases, the hydrodynamic force counteracts the attraction of the permanent magnet bearing and is in the force balance under the action of the inlet fluid, so that the impeller is separated from the contact with the pump casing and is in a suspended state. Additionally, the hydrodynamic force is slightly greater than that of the permanent magnet bearing. Therefore, the axial displacement is small relative to the gap.

5. Conclusions

The blood pump with a suspension system eliminates the problems of wear, heat, and hemolysis caused by the mechanical bearings and greatly improves the blood compatibility. However, the realization of a stable suspension of the rotor plays a vital role in improving the durability and blood compatibility of the blood pump. In the present study, the magnetic suspension system composed of a blood pump with single-winding bearingless motor and permanent magnet bearing was designed, the axial dynamic model of bearingless motor rotor was established, and the accuracy of the electromagnetic force and the torque theoretical model was verified using a simulation. The radial suspension stability of the rotor was analyzed, the axial dynamics model of the rotor with a magnetic–hydrodynamic double suspension system was established, the suspension stability was analyzed from the perspective of stiffness, and a hydrodynamic bearing was designed to eliminate the influence of axial negative stiffness on the suspension stability. By solving the Reynolds equation, the influence law of the spiral groove structural parameters on stiffness and hydrodynamic force was studied. When the number of grooves was 20, the spiral angle was 45°, the groove width ratio was 0.5, and the groove depth was 0.1 mm, the system exhibited an optimal axial support performance. Finally, the suspension stability under the fluid impact was studied by building an extracorporeal mock circulatory loop. The experimental results showed that the bearingless rotor could float stably and ensured the suspension stability under the impact of external fluid, which laid a foundation for the follow-up active suspension control and practical clinical application. In addition, the structural parameters and axial clearance control of the hydrodynamic bearing affected the blood compatibility of blood pump. Therefore, the structure of magnetic suspension system and hydrodynamic bearing needs to be optimized by integrating the suspension performance and blood compatibility.

Author Contributions

Y.W. and P.S. conceptualized the study; P.F. and L.Z. performed the formal analysis; P.S. worked on the methodology and software-related evaluations; L.L. and Y.C. wrote, reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the National Natural Science Foundation of China (grant number 51875143).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Han, J.; Trumble, D.R. Cardiac Assist Devices: Early Concepts, Current Technologies, and Future Innovations. Bioengineering 2019, 6, 18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Simon, M.A.; Bachman, T.N.; Watson, J.; Baldwin, J.T.; Wagner, W.R.; Borovetz, H.S. Current and Future Considerations in the Use of Mechanical Circulatory Support Devices: An Update, 2008–2018. In Annual Review of Biomedical Engineering; Yamush, M.L., Ed.; Annual Reviews: San Mateo, CA, USA, 2019; Volume 21, pp. 33–60. [Google Scholar]
  3. Ward, R.A.; Wellhausen, S.R.; Dobbins, J.J.; Johnson, G.S.; Devries, W.C. Thromboembolic and Infectious Complications of Total Artificial Heart Implantation. Ann. N. Y. Acad. 2010, 516, 638–650. [Google Scholar] [CrossRef] [PubMed]
  4. Sheikh, F.H.; Russell, S.D. HeartMate II continuous-flow left ventricular assist system. Expert Rev. Med. Devices 2011, 8, 11–21. [Google Scholar] [CrossRef] [PubMed]
  5. Kounis, N.G.; Soufras, G.D.; Davlouros, P.; Tsigkas, G.; Hahalis, G. Thrombus Formation Patterns in HeartMate II Continuous-Flow Left Ventricular Assist Devices: A Multifactorial Phenomenon Involving Kounis Syndrome? Asaio J. 2014, 60, 369–371. [Google Scholar] [CrossRef] [PubMed]
  6. Bourque, K.; Gernes, D.B.; Loree, H.M.; Richardson, J.S.; Poirier, V.L.; Barletta, N.; Fleischli, A.; Foiera, G.; Gempp, T.M.; Schoeb, R. HeartMate III: Pump design for a centrifugal LVAD with a magnetically levitated rotor. Am. Soc. Artif. Intern. Organs J. 2001, 47, 401–405. [Google Scholar] [CrossRef] [PubMed]
  7. Hoshi, H.; Katakoa, K.; Ohuchi, K.; Asama, J.; Shinshi, T.; Shimokohbe, A.; Takatani, S. Magnetically suspended centrifugal blood pump with a radial magnetic driver. ASAIO J. 2005, 51, 60–64. [Google Scholar] [CrossRef] [PubMed]
  8. Hoshi, H.; Shinshi, T.; Takatani, S. Third-generation blood pumps with mechanical noncontact magnetic bearings. Artif. Organs 2010, 30, 324–338. [Google Scholar] [CrossRef] [PubMed]
  9. Masuzawa, T. Magnetically suspended motor system applied to artificial hearts and blood pumps. Proc. Inst. Mech. Eng. Part I J. Syst. Control. Eng. 2017, 231, 330–338. [Google Scholar] [CrossRef]
  10. Chen, Z.F.; Li, Z.X.; Ma, H.Z. Design and performance characterization of a novel salient-pole surface-mounted permanent magnet synchronous machine. Int. J. Appl. Electromagn. Mech. 2018, 56, 347–355. [Google Scholar] [CrossRef]
  11. Chou, P.; Yang, S.; Jhong, C.; Huang, J.; Chen, J. High-Precision Permaent Magent Motor Design for Satellite Attitude Control with High Torque Density and Low Torque Ripple. In Proceedings of the IECON 2019—45th Annual Conference of the IEEE Industrial Electronics Society, Lisbon, Portugal, 14–17 October 2019; pp. 802–805. [Google Scholar]
  12. Nasiri-Zarandi, R.; Ghaheri, A.; Abbaszadeh, K. Cogging Torque Reduction in U-Core TFPM Generator Using Different Halbach-Array Structures. In Proceedings of the International Symposium on Power Electronics, Electrical Drives, Automation and Motion (SPEEDAM), Amalfi, Italy, 20–22 June 2018; pp. 1153–1158. [Google Scholar]
  13. Fu, Y.; Hu, Y.; Huang, F.; Zhou, M. The Impact of Pulsatile Flow on Suspension Force for Hydrodynamically Levitated Blood Pump. J. Healthc. Eng. 2019, 2019, 8065920. [Google Scholar] [CrossRef] [PubMed]
  14. Kosaka, R.; Yoshida, F.; Nishida, M.; Maruyama, O.; Yamane, T. Bearing gap adjustment for improvement of levitation performance in a hydrodynamically levitated centrifugal blood pump. In Proceedings of the Engineering in Medicine & Biology Society, Milan, Italy, 25–29 August 2015. [Google Scholar]
  15. Wang, L.; Tang, X.; Yun, Z.; Xiang, C.; Lu, Z. Analysis of a Novel Magnetic-Hydrodynamic Double Levitated Motor for an Implantable Axial Flow Blood Pump. Shock. Vib. 2020, 2020, 8833994. [Google Scholar] [CrossRef]
  16. Chung, M.; Nong, Z.; Tansley, G.D.; Yi, Q. Experimental Determination of Dynamic Characteristics of the VentrAssist Implantable Rotary Blood Pump. Artif. Organs 2015, 28, 1089–1094. [Google Scholar] [CrossRef] [PubMed]
  17. Kosaka, R.; Yasui, K.; Nishida, M.; Kawaguchi, Y.; Maruyama, O.; Yamane, T. Optimal bearing gap of a multiarc radial bearing in a hydrodynamically levitated centrifugal blood pump for the reduction of hemolysis. Artif. Organs 2014, 38, 818–822. [Google Scholar] [CrossRef] [PubMed]
  18. Yasui, K.; Kosaka, R.; Nishida, M.; Maruyama, O.; Kawaguchi, Y.; Yamane, T. Optimal design of the hydrodynamic multi-arc bearing in a centrifugal blood pump for the improvement of bearing stiffness and hemolysis level. Artif. Organs 2013, 37, 768–777. [Google Scholar] [CrossRef] [PubMed]
  19. Kataoka, H.; Kimura, Y.; Fujita, H.; Takatani, S. Measurement of the Rotor Motion and Corresponding Hemolysis of a Centrifugal Blood Pump With a Magnetic and Hydrodynamic Hybrid Bearing. Artif. Organs 2015, 29, 547–556. [Google Scholar] [CrossRef] [PubMed]
  20. Nguyen, Q.D.; Ueno, S. Analysis and Control of Nonsalient Permanent Magnet Axial Gap Self-Bearing Motor. IEEE Trans. Ind. Electron. 2011, 58, 2644–2652. [Google Scholar] [CrossRef]
  21. Asama, J.; Hamasaki, Y.; Oiwa, T.; Chiba, A. Proposal and Analysis of a Novel Single-Drive Bearingless Motor. IEEE Trans. Ind. Electron. 2013, 60, 129–138. [Google Scholar] [CrossRef]
  22. Cheng, S.; Olles, M.W.; Burger, A.F.; Day, S.W. Optimization of a Hybrid Magnetic Bearing for a Magnetically Levitated Blood Pump via 3-D FEA. Mechatronics 2011, 21, 1163–1169. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Pohlmann, A.; Hameyer, K. A study on permanent magnet topologies for hybrid bearings for medical drives applied in Ventricular Assist Devices. Arch. Electr. Eng. 2011, 60, 371–380. [Google Scholar] [CrossRef]
Figure 1. Magnetic–hydrodynamic hybrid suspension blood pump: (a) Structural composition of the blood pump; (b) Three-dimensional casing model.
Figure 1. Magnetic–hydrodynamic hybrid suspension blood pump: (a) Structural composition of the blood pump; (b) Three-dimensional casing model.
Machines 09 00255 g001
Figure 2. Bearingless motor of the blood pump: (a) Single-winding bearingless motor; (b) Motor stator; (c) Rotor with Halbach array.
Figure 2. Bearingless motor of the blood pump: (a) Single-winding bearingless motor; (b) Motor stator; (c) Rotor with Halbach array.
Machines 09 00255 g002
Figure 3. Hybrid magnetic suspension system with a single-winding bearingless motor and a permanent magnet bearing: (a) Hybrid magnetic suspension system; (b) permanent magnet bearing.
Figure 3. Hybrid magnetic suspension system with a single-winding bearingless motor and a permanent magnet bearing: (a) Hybrid magnetic suspension system; (b) permanent magnet bearing.
Machines 09 00255 g003
Figure 4. Bearingless motor structure and working principle diagram of the blood pump: (a) Rotor and stator of the bearingless motor; (b) Schematic representation of the working principle diagram.
Figure 4. Bearingless motor structure and working principle diagram of the blood pump: (a) Rotor and stator of the bearingless motor; (b) Schematic representation of the working principle diagram.
Machines 09 00255 g004
Figure 5. Structure composition of the permanent magnet bearing.
Figure 5. Structure composition of the permanent magnet bearing.
Machines 09 00255 g005
Figure 6. Variation law of the magnetic field and axial force during the axial offset: (a) Magnetic flux density and magnetic induction line distribution of the permanent magnet bearing; (b) Relationship curve between axial force and axial displacement.
Figure 6. Variation law of the magnetic field and axial force during the axial offset: (a) Magnetic flux density and magnetic induction line distribution of the permanent magnet bearing; (b) Relationship curve between axial force and axial displacement.
Machines 09 00255 g006
Figure 7. Parameter identification platform of the bearingless motor.
Figure 7. Parameter identification platform of the bearingless motor.
Machines 09 00255 g007
Figure 8. Variation curve of the electromagnetic torque and axial force during motor rotation: (a) Axial force of the bearingless motor; (b) Electromagnetic torque of the bearingless motor.
Figure 8. Variation curve of the electromagnetic torque and axial force during motor rotation: (a) Axial force of the bearingless motor; (b) Electromagnetic torque of the bearingless motor.
Machines 09 00255 g008
Figure 9. Relationship curve between the permanent magnetic force and the impeller displacement: (a) Variation law of the permanent magnetic force under axial displacement; (b) Variation curve of axial stiffness with axial displacement.
Figure 9. Relationship curve between the permanent magnetic force and the impeller displacement: (a) Variation law of the permanent magnetic force under axial displacement; (b) Variation curve of axial stiffness with axial displacement.
Machines 09 00255 g009
Figure 10. Structure composition of the hydrodynamic bearing: (a) Structural parameters of the hydrodynamic bearing; (b) Cross-sectional view of the hydrodynamic bearing.
Figure 10. Structure composition of the hydrodynamic bearing: (a) Structural parameters of the hydrodynamic bearing; (b) Cross-sectional view of the hydrodynamic bearing.
Machines 09 00255 g010
Figure 11. Pressure distribution of the hydrodynamic bearing surface.
Figure 11. Pressure distribution of the hydrodynamic bearing surface.
Machines 09 00255 g011
Figure 12. Pressure distribution of the hydrodynamic bearing: (a) Radial pressure distribution; (b) Circumferential pressure distribution.
Figure 12. Pressure distribution of the hydrodynamic bearing: (a) Radial pressure distribution; (b) Circumferential pressure distribution.
Machines 09 00255 g012
Figure 13. Influence of the groove number and groove depth on stiffness and hydrodynamic force: (a) Number of the spiral grooves; (b) Spiral groove depth.
Figure 13. Influence of the groove number and groove depth on stiffness and hydrodynamic force: (a) Number of the spiral grooves; (b) Spiral groove depth.
Machines 09 00255 g013
Figure 14. Influence of the groove width ratio and spiral angle on the stiffness and hydrodynamic force: (a) Groove width ratio; (b) Spiral angle.
Figure 14. Influence of the groove width ratio and spiral angle on the stiffness and hydrodynamic force: (a) Groove width ratio; (b) Spiral angle.
Machines 09 00255 g014
Figure 15. Processed blood pump casing with a hydrodynamic bearing.
Figure 15. Processed blood pump casing with a hydrodynamic bearing.
Machines 09 00255 g015
Figure 16. Extracorporeal mock circulatory loop.
Figure 16. Extracorporeal mock circulatory loop.
Machines 09 00255 g016
Figure 17. Axial displacement and current waveform corresponding to different rotating speeds of the blood pump impeller: (a) Axial displacement; (b) Current waveform.
Figure 17. Axial displacement and current waveform corresponding to different rotating speeds of the blood pump impeller: (a) Axial displacement; (b) Current waveform.
Machines 09 00255 g017
Figure 18. Hydraulic efficiency of blood pump corresponding to different speeds and flowrates.
Figure 18. Hydraulic efficiency of blood pump corresponding to different speeds and flowrates.
Machines 09 00255 g018
Figure 19. Pressure corresponding to external fluid disturbance at constant speeds: (a) Pressure waveform; (b) Axial displacement.
Figure 19. Pressure corresponding to external fluid disturbance at constant speeds: (a) Pressure waveform; (b) Axial displacement.
Machines 09 00255 g019
Figure 20. Effects of different external disturbances on suspension performance of the rotor at constant speeds: (a) Pressure waveform; (b) Axial displacement.
Figure 20. Effects of different external disturbances on suspension performance of the rotor at constant speeds: (a) Pressure waveform; (b) Axial displacement.
Machines 09 00255 g020
Table 1. Specifications of the magnetic suspension system.
Table 1. Specifications of the magnetic suspension system.
SymbolExplanationsSize
RsoOuter radius of stator20 (mm)
RpoOuter radius of pole shoe15 (mm)
RpiInner radius of pole shoe8 (mm)
RcCoil radius6 (mm)
RiIron core radius2 (mm)
γAngle of pole shoe15 (°)
RriInner radius of rotor7 (mm)
RroOuter radius of rotor16 (mm)
αrAngle of magnet pole22.5 (°)
hpsThickness of static magnetic ring1.3 (mm)
hpmThickness of movable magnetic ring0.8 (mm)
hrRotor backplate thickness1(mm)
hchCore height7 (mm)
hmMagnetic pole thickness1 (mm)
hsThickness of pole shoe0.5 (mm)
hyYoke thickness2 (mm)
R1iInner radius of magnetic bearing8 (mm)
R1oOuter radius of magnetic bearing10 (mm)
R2iInner radius of magnetic bearing11.5 (mm)
R2oOuter radius of magnetic bearing13.5 (mm)
Table 2. Material parameters of the permanent magnet bearing.
Table 2. Material parameters of the permanent magnet bearing.
ParameterSymbolValueUnit
Residual flux densityBr1.3T
CoercivityHc960kA/m
Magnetic energyBHmax40MGOe
Table 3. Motor parameter identification results.
Table 3. Motor parameter identification results.
ParameterSymbolValueUnit
Rated voltageVbus24V
Rated currentI1.8A
Rated speedω4500rpm
Armature windingRs1.49Ω
Winding inductanceLs0.72mH
Polar pairsp4
Flux linkage of permanent magnetλm2.81 × 10−3Wb
Moment of inertiaJ4.5 × 10−6kg·m2
Damping coefficientB2.05 × 10−5N·s/m
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shen, P.; Wang, Y.; Chen, Y.; Fu, P.; Zhou, L.; Liu, L. Hydrodynamic Bearing Structural Design of Blood Pump Based on Axial Passive Suspension Stability Analysis of Magnetic–Hydrodynamic Hybrid Suspension System. Machines 2021, 9, 255. https://doi.org/10.3390/machines9110255

AMA Style

Shen P, Wang Y, Chen Y, Fu P, Zhou L, Liu L. Hydrodynamic Bearing Structural Design of Blood Pump Based on Axial Passive Suspension Stability Analysis of Magnetic–Hydrodynamic Hybrid Suspension System. Machines. 2021; 9(11):255. https://doi.org/10.3390/machines9110255

Chicago/Turabian Style

Shen, Peng, Yiwen Wang, Yun Chen, Pengqiang Fu, Lijie Zhou, and Lijia Liu. 2021. "Hydrodynamic Bearing Structural Design of Blood Pump Based on Axial Passive Suspension Stability Analysis of Magnetic–Hydrodynamic Hybrid Suspension System" Machines 9, no. 11: 255. https://doi.org/10.3390/machines9110255

APA Style

Shen, P., Wang, Y., Chen, Y., Fu, P., Zhou, L., & Liu, L. (2021). Hydrodynamic Bearing Structural Design of Blood Pump Based on Axial Passive Suspension Stability Analysis of Magnetic–Hydrodynamic Hybrid Suspension System. Machines, 9(11), 255. https://doi.org/10.3390/machines9110255

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop