轎車膜片彈簧離合器系統(tǒng)設(shè)計(jì)及仿真分析,轎車膜片彈簧離合器系統(tǒng)設(shè)計(jì)及仿真分析,轎車,膜片,彈簧,離合器,系統(tǒng),設(shè)計(jì),仿真,分析 附錄B Energies 2015, 8, 5598-5612; doi:10.3390/en8065598 OPEN ACCESS energies ISSN 1996-1073 www.mdpi.com/journal/energies Article Analysis and Design of a Permanent Magnet Bi-Stable Electro-Magnetic Clutch Unit for In-Wheel Electric Vehicle Drives Wanli Cai *, Chenglin Gu and Xiaodong Hu State Key Laboratory of Advanced Electromagnetic Engineering and Technology, School of Electrical & Electronic Engineering, Huazhong University of Science and Technology, Wuhan 430074, China; E-Mails: clgu@mail.hust.edu.cn (C.G.); xdhu@hust.edu.cn (X.H.) * Author to whom correspondence should be addressed; E-Mail: wlcai@hust.edu.cn; Tel.: +86-134-3718-8225. Academic Editor: Paul Stewart Received: 9 March 2015 / Accepted: 1 June 2015 / Published: 11 June 2015 Abstract: Clutches have been used in internal combustion vehicles and concentrated electric vehicles (EVs) to smoothen impulsion while starting and shifting. This paper proposes a permanent magnet bi-stable electromagnetic clutch unit (PMBECU) which is specially introduced into in-wheel EVs to make the rigid connection between hub and wheel more flexible. Firstly, the operation principle of the PMBECU is illustrated. Then, the basic magnetic circuit model is presented and analyzed, followed by optimal design of the main structural parameters by investigating the PM leakage flux coefficient. Further, according to the basic electromagnetic characteristics of the PMBECU, the current pulse supply is put forward, and the minimum pulse width which enables the operation of the PMBECU and its dynamic characteristics are analyzed by an improved finite element method. Finally, a prototype machine is manufactured and tested to validate all the analysis results. Keywords: clutch unit; dynamics analysis; electromagnetic design; finite element method (FEM); permanent magnet; bi-stable operation Energies 2015, 8 5599 1. Introduction Electric vehicles (EV) have been intensively investigated recently as potential solutions for the growing problems of the energy crisis and environmental pollution [1–4], focusing on the drive form, electric motor, controller, battery, energy system, drive comfort, etc. Compared with centralized drive, the in-wheel EV drive is considered the more competent drive form for EVs in the near future [5–7], because of its merits of direct drive (no-gearbox), more flexible control strategy (torque at each wheel is independently controlled), high mechanical integrity (greatly different from conventional gasoline cars). However, the rigid connection between hub and motor, inevitably introduces mechanical shocks and electromagnetic impulsion during sudden start and stop processes, which can potentially harm the motor and controller and reduce drive comfort [8–10]. Referring to traditional gasoline cars, this electromechanical impulsion in in-wheel EV drives can be ameliorated by introducing a clutch between the hub and motor to make the rigid connection more flexible [11]. The simulation and experimental results of a conventional clutch between motor and load presented in [12,13] show that the starting current and jerk in clutch coupling starts under different idle speeds can be reduced to less than 1/2 compared to direct starting, and the impulsive back electromotive force to the controller can be eliminated by detached braking (the motor stops naturally after being disconnected from the braking load). Besides, in hybrid EVs, the conventional clutch has been used to cut off the engine or electrical machine while idling to avoid spin losses and extend the life cycle of the machine [14]. Moreover, in in-wheel driven EVs, clutches have been used to detach the motor from hub to reduce losses while coasting [15]. However, the conventional mechanical clutch system [16,17] is not suitable for the limited space available in a hub and suffers from a need for regular maintenance which makes it unsuitable for in-wheel EV drives. In addition, electromagnetic clutches [18,19], which can be easily manipulated by current control, are energy-consuming and also suffer from the problem of accommodating their shape in the hub. In other clutches [20] one encounters one or all of the aforementioned problems, thus are also not suitable options. This paper proposes a permanent magnet bi-stable electromagnetic clutch unit (PMBECU), which is controlled by current and held by the PM in a steady state, and thus is energy-saving, and it also has a flat structure that makes its placement in a limited space viable. The clutch system is realized by assembling several PMBECUs around motors, combined with friction or jaw pairs. As key parts of the clutch system, this paper focuses on the electromagnetic design and analysis of the PMBECU. The design and analysis of linear electromagnetic devices, such as electromagnetic valves [21], electric tools [22], oscillators [23,24], and switch gears [25–27], are mainly carried out by the finite element method (FEM). Likewise, aiming to satisfy the need to accommodate the clutch in the limited space available in the hub, the optimal design of the main structure parameters of the PMBECU are carried out by FEM which focuses especially on investigating the leakage flux coefficient of the PM. Moreover, in order to realize simple and reliable control of the operation of the PMBECU, the dynamic characteristics of the PMBECU are calculated by improved FEM, which shows that the low power capacitor pulse supply is very suitable. The influence of the temperature on the dynamic performance is also analyzed. The analysis method and results are finally validated by measurements taken on a prototype machine. Energies 2015, 8 5600 2. Operation Principle The assembly of the PMBECUs to realize the flexible connection between hub and motor is shown in Figure 1a, and the structure of the PMBECU, in which two PMs with opposite polarities are mounted on each side of a rigid E-type ferromagnetic base, is shown in Figure 1b. The ferromagnetic mover is placed in two low-frictional slideways which are non-magnetic. Two coils are connected n series and wound around each slideway. Tyre PMBECU Mover Coil PM Wheel motor PM Slideway Base Slideway (a) (b) Figure 1. (a) Flexible connection of hub and motor; (b) Structure of the PMBECU. The 2-dimensional (2D) analysis model of the PMBECU with its main structure parameters labeled is shown in Figure 2, where the right direction is prescribed as positive for force and movement variables. FEM boundary hm wm fz mg fmx hp fmy x(d1) d2=lt-d1 o Figure 2. 2D analysis model of the PMBECU. The flux line distribution of the PMBECU without current injected into the coils is shown in Figure 3a. Apparently, the mover is held by the left PM in a steady state without energy consumption. When current with a suitable orientation (i.e., the current direction shown in Figure 2) and value accesses the coils, the mover is polarized, and the corresponding flux lines distribution is shown in Figure 3b. The mover will soon be propelled from the left steady state to the right by the resultant electromagnetic force. Meanwhile, the current is switched off automatically by the position sensor, and the mover is held by the right PM, again without any energy consumption, thus it is bi-stable. Energies 2015, 8 5601 (a) (b) Figure 3. Magnetic flux lines distribution. (a) Steady state; (b) Action. It is evident that the PMBECU has a flat structure thus is suitable for placement in a limited space, and the switchover between engagement and disengagement is electrically-controlled thus it can be conveniently manipulated, and only an instant current is required for switchover, but most time it is in a steady state which is held by a PM and thus is energy-saving. 3. Electromagnetic Design 3.1. Magnetic Circuit Model According to the magnetic flux line distribution shown in Figure 3a, assuming the ferromagnetic material has infinite permeability and neglecting the contact air gaps, the magnetic circuit relations of the PMBECU under open circuit of the coils conditions can be expressed by a simplified magnetic network as shown in Figure 4. L d1 L d2 Fm1 F d1 F d2Fm2 Figure 4. Simplified magnetic network. The magnetic network comprises two independent branches, where Φδj are the magneto-motive force furnished to the air gap by the PM, magnetic flux pass through the pole face of the mover at each side, Equations (1)–(3), respectively:  Fmj (j = 1, 2), Λδj, and air gap permeance, and which are calculated by Fmj = H c hm δ j ( kσj hm μr + δ j ) (1) L δj = μ0 Sm δ j (2) Energies 2015, 8 5602 F δj = B r S m ? δ j ? (3) ? kσj + μr ÷ è hm ? where Br, Hc, and μr are remanence, coercivity, and relative permeability of the PM, hm and Sm are the thickness and pole face area of the PM, δj is the air gap length as labeled in Figure 2, μ0 is the permeability of air, and kσj is the leakage flux coefficient which is defined as: kσj = F mj F δj (4) where Φmj is main magnetic flux through bottom face of PM. The Maxwell stress tensors are given by the following equation [27]: t n = (B n2 - Bs2 ) (2μ0 ) (5) t s = B n Bs μ0 where Bn, Bs are the outer normal and tangential components of the flux density on the mover, respectively. Out of an infinite permeable surface, the flux density only has a normal component. Hence, combined with Equation (3), the holding force (horizontal) at steady state (δ1 = 0, δ2 = lt, lt is the travel length of the mover) can be approximately calculated by: ? ? f H = B δ21 S m - B δ22 S m = B r2 Sm ? 1 - 1 ÷ ? ÷ (6) 2 μ 0 2 μ 0 2μ0 ( lt hm ) 2 è ? ? kσ2 + μr ÷ ? ÷ With forces normalized to fb = 0.5Br2Sm/μ0 (the same hereafter), the holding force is: f H = 1- 1 (kσ2 + μr lt hm ) 2 (7) Apparently, the holding force of the PMBECU is determined by kσ2 (leakage coefficient at δ = lt), the ratio of travel length to thickness of the PM lt/hm, and the PM characteristics. Moreover, the leakage flux coefficient kσ2 is a function of the structure parameters, and can be calculated by Equation (4) after the magnetic flux derived from FEM analysis. By increasing the current from 0, the electromagnetic force experienced by the mover can be obtained, and then the ideal threshold current iT which critically enables the action of the mover can be obtained by FEM as well, corresponding to the horizontal electromagnetic force fmx = 0. In this paper, current is all normalized to ib = Hchm/N, where N is the number of turns for one coil. 3.2. Main Structure Parameters Design The PMBECU works at steady state most of the time, which is reliably maintained by the holding force, thus the holding force is the most significant index. According to Equation (7), the leakage flux coefficient kσ2 at the detached side, which is a function of the structure parameters, has a great influence on the holding force. Moreover, the leakage flux coefficient determines the reasonable usage of the PM. Thus, the main structure parameters (as labelled in Figure 2), i.e., the width of the PM wm, the height from the PM to the base hp, and the travel length of mover lt, are optimized by studying kσ2, Energies 2015, 8 5603 combined with accounting for the holding force and threshold current, where, other size ratios (proportioned to hm) remain unchanged while one varies the parameters within the ranges hp/hm = 1.2, wm/hm = 2.5, lt/hm = 2. The variation of kσ2 versus different structure parameters is shown in Figure 5. Figure 5a shows that the leakage flux coefficient increases quite slowly when hp is 1.5 times bigger than hm, hence hp would better be within 1–1.5 times of hm, which also indicates the PMBECU is capable of a flat structure. Likewise, wm would better be around 2.5 times of hm as seen in Figure 5b. Figure 5c shows the leakage flux coefficient kσ2 increases almost linearly with lt, which shows no clear inflection point. But from Figure 5d, the holding force increases very slowly when lt is 2 times larger than hm, meanwhile, the threshold current keeps increasing, which makes the action of the mover harder. Hence, lt within 1.5–2 times the thickness of the PM is more sensible. Leakage flux coefficient kσ2 Leakage flux coefficient kσ2  5 4 3 2 1 0 0 1 2 3 Height of PM to base/thickness of PM hp/hm (a) 7 6 5 4 3 2 1 0 0 1 2 3 4 Travel length/thickness of PM lt/hm (c)  10 σ2 8 k coefficient 6 flux 4 Leakage 2 0 0 1 2 3 4 5 6 7 Width of PM/thickness of PM wm/hm (b) value) 1 1 value) 0.8 0.8 (P.U. (P.U. 0.6 fH T 0.6i H force f 0.4 0.4 iT current Holding 0.2 0.2 0 0 Threshold 1 2 3 0 4 Travel length / thickness of PM lt/hm (d) Figure 5. Opimization. (a) Height from PM to base; (b) Width of PM; (c) Travel length; (d) Travel length. 4. Dynamics Analysis 4.1. Electromagnetic Characteristics Based on the aforementioned analyses, a PMBECU prototype designed with the main parameters listed in Table 1 is shown in Figure 6. Assuming the mover is fixed at different positions, changing the Energies 2015, 8 5604 current (constant DC wave) in the coil, computing the magnetic field by FEM and the forces experienced by mover by (5), then the electromagnetic forces on the mover versus current i and displacement x are obtained as shown in Figure 7. Table 1. Leading design parameters. Parameter Value Parameter Value Thickness of PM hm 2.5 mm Width of base wb 80 mm Width of PM wm 6 mm Remanence of PM Br 0.4 T Length of PM lm 20 mm Coercivity of PM Hc 318 kA/m Height of PM to base hp 3 mm Turns of coil N 60 Travel length lt 4.8 mm Mass of mover m 56 g Figure 6. Prototype. From Figure 7a, for open circuit conditions, the i = 0 horizontal force curve indicates that the PMBECU has two steady states held by the magnetic force from the PM, and an unstable equilibrium point (the half travel length location). When the mover exceeds this unstable point, the mover can be drawn to the other steady state automatically even if the current is switched off. Since the current increases to the ideal threshold current (enabling the action of the mover) i = 0.49, the mover starts moving. The maximum current in the coil is limited by the inflection point of the demagnetizing curve of the PM (critical point of irreversible demagnetization), which is i = 0.77 in this prototype. In fact, because of the asymmetric structure in the vertical direction, the mover experiences a downward vertical electromagnetic force (as shown in Figure 7b) which introduces frictional resistance. Hence, accounting for friction, and other errors (material, model, measuring, etc.), the real threshold current iT is bigger than the calculated value, which is iT = 0.52 for the prototype. Moreover, to guarantee the performance of the PM, the maximum current should be limited to iM = 0.7. Energies 2015, 8 5605 Horizontal force fmx (P.U. value)  3 2 1 0 -1 0  i = 0.7 i = 0.49 i = 0 1.2 2.4 3.6 4.8 Displacement x (mm) (a)  Vertical force fmy (P.U. value)  2 1.5 1 0.5 0 0  i = 0.7 i = 0.49 i = 0 1.2 2.4 3.6 4.8 Displacement x (mm) (b) Figure 7. Electromagnetic forces. (a) Horizontal; (b) Vertical. When the current is larger than the threshold current, the resultant positive horizontal force starts to drive the mover, and the force is a monotonously increasing function of the displacement. After moving through the middle point of the PMBECU, the mover can reach another steady state with the current switched off (i.e., the pulse current sustains only half the travel length width). What’s more, considering the inertial motion and variation of the kinetic friction coefficient, the pulse width of the current can be even smaller. Thus, a dynamics analysis of the PMBECU should be carried out. 4.2. Dynamics Equations and Analysis Method Because of the motion symmetry of the PMBECU, only the movement of the mover from left to right is investigated. Supposing the static friction coefficient is equal to the kinetic friction coefficient, then the magnetic-kinematic coupled mathematic equations which determines the dynamics characteristics are described as: f mx - f z = mdv dt (8) v = dx dt (9) f mx = q ( x, i ), f my = p ( x, i) (10) f z = μs ( f my + mg) (11) where fmx and fmy are the horizontal and vertical electromagnetic forces on the mover, fz is the resisting force, v is the velocity of the mover, μs is the static friction coefficient which is 0.065 in this prototype (measured), and g is the acceleration constant of gravity. The dynamics analysis of the PMBECU is to illustrate the coupling of the magnetic field and the movement. To cope with the varying friction resistance conditions of the PMBECU, and give consideration to the convenience of analysis of varied structure sizes, an improved FEM is proposed. As shown in Figure 8a, two lt length rectangular areas (namely, the material variation area) in proximity to the PMs are established and uniformly meshed into n steps of quadrilateral shape, i.e., the step length is x = lt/n. The initial permeability of the left part and the right are set as iron (μFe) and air Energies 2015, 8 5606 (μ0) respectively. As shown in Figure 8b, if the permeability of the first x meshes in the left material variation area is changed into μ0 and the first x meshes at the right into μFe, a x displacement of the mover is equivalently realized. Thus, a onetime mesh can cover the travel length displacement of the mover [23]. Y Z X (a) n·Dx μFe μ0 μ0 μFe μFe μ0 (b) Figure 8. Onetime mesh technique. (a) Mesh; (b) Principle. Further, by setting displacement as a known quality but time as an unknown variable, and calculating the time, velocity, and current before each time of material variation, the whole PMBECU movement process (i.e., the dynamics characteristics of the PMBECU) can be solved by using only a onetime mesh. This improved FEM analysis flow chart is shown in Figure 9, where both the current change and resistance variation can be taken into account, w