Next Article in Journal
SRV Method: Lubricating Oil Screening Test for FZG
Previous Article in Journal
Review on Torque Distribution Scheme of Four-Wheel In-Wheel Motor Electric Vehicle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Response and Synchronizing Characteristic for the Dual-Motor Driving System in Non-Inertial System

College of Mechanical Engineering, Chongqing University of Technology, Chongqing 400054, China
*
Author to whom correspondence should be addressed.
Machines 2022, 10(8), 620; https://doi.org/10.3390/machines10080620
Submission received: 24 June 2022 / Revised: 21 July 2022 / Accepted: 26 July 2022 / Published: 28 July 2022
(This article belongs to the Section Electrical Machines and Drives)

Abstract

:
As one of the typical propulsion systems of the new-energy vehicles (NEVs), the dual-motor driving system (DMDS) is installed and fixed on the bodywork and moves together with the bodywork in space, that is, the DMDS works in non-inertial systems. The previous researches on the DMDS are based on the assumption that the bodywork is stationary. In fact, since the DMDS moves with the bodywork, besides its own excitation forces, it is inevitably affected by the additional inertial terms (AITs) caused by the change in the operating state of the NEVs. In order to investigate the dynamic response and synchronizing characteristic of the DMDS under different motion forms of NEVs, the dynamic model of the DMDS in a non-inertial system is established, considering the permanent magnet synchronous motors (PMSMs), time-varying meshing stiffness and transmission error of gears. Subsequently, the expressions of the AITs are deduced under different non-inertial conditions. The translational motion and circular motion of the vehicle are selected to analyze the dynamic response and synchronizing characteristics of the DMDS in a non-inertial system. The results show that the acceleration has a significant influence on the displacement response of gears, torque and speed of PMSMs, but the torque synchronization and speed synchronization between PMSMs are minimally influenced by both the acceleration and the AITs. Meanwhile, the AITs that affect the displacement response are analyzed and quantified.

1. Introduction

It is generally believed that vehicles that use fossil fuels are one of the top sources of air pollution, and the depletion of oil is another serious problem that people have to face. As a consequence, new-energy vehicles (NEVs) are increasingly needed and have been listed as a development priority for all countries. Normally, the NEVs are driven by electric drive gear systems, such as the combination of a single motor or multiple motors with a gear transmission system. As one of the typical propulsion systems of NEVs, the dual-motor driving system (DMDS) has become the focus of more and more people, and has been reported in many literatures in the past. For example, Yang et al. designed a high power dual motor drive system to resolve the limitation of permanent magnet synchronous motor (PMSMs) application (such as temperature characteristics of permanent magnet materials and flux-weakening control problems) and utilize its good performances [1] and Tang provided a method and apparatus for optimizing the torque applied by each motor of a dual motor drive system of an all-electric vehicle [2]. In the traditional research on the dynamic response or synchronization property of DMDS, the DMDS is considered to be stationary relative to the earth reference system. In fact, the DMDS is installed and fixed on the bodywork, and moves together with the bodywork in space, that is, the DMDS operates in a non-inertial system. Besides its own excitation forces, the DMDS is inevitably affected by additional inertial forces and gyroscopic moments, which affects the dynamic response and synchronization property of DMDS. Therefore, the effects of non-inertia should be considered when the NEVs perform a spatial motion.
In the past few decades, many scholars across the world have carried out meaningful research on the combination of a single motor or multiple motors with a gear transmission system attached to a moving body, including modeling, dynamic response and synchronization property, etc. In modeling and dynamic property terms, Inalpolat et al. established a transverse-torsional dynamic mode of a spur gear reducer, including periodically time-varying meshing stiffness and nonlinearities caused by tooth separations in resonance regions, and the influence of gear tooth indexing errors on the dynamic response was investigated [3]. Giagopulos et al. proposed a dynamic model of a gear-pair system considering the gear backlash and bearing stiffness non-linearities, and then the effect of non-linearities in the identification and fault detection of gear-pair system were investigated [4]. Setiawan et al. developed a mathematical model for multispeed transmissions in electric vehicles, and the transient dynamic response of the transmission was studied. The established model was subsequently validated by comparing its results with experiment results [5]. Sun et al. developed a coupling dynamic model of a multi-motor gear system by employing the proposed hierarchical modeling method and obtained the influence rules of system parameters on the natural frequency. Meanwhile, dynamic behaviors under different load conditions were also obtained [6]. Liu et al. used the proposed hybrid user-defined element method to establish the dynamic model of a gear transmission system in an electric vehicle, and investigated the effect of meshing impact force and meshing stiffness on the dynamic behavior of the system [7]. Tang et al. established a novel torsional vibration dynamic model of a hybrid power train, and the effect of a dual-mass flywheel on the torsional vibration characteristics of a power-split hybrid power train was studied [8]. Zhang et al. established a multi-degree-of freedom (M-D-O-F) dynamic model of a two-stage gear set that included the time-variant stiffness and gear backlash feature, and then the nonlinear response of the two-stage gear set under the combined periodic and stochastic excitations was investigated [9]. Yan et al. presented a new design by integrating a planetary gear train (PGT) within a DC motor to form a compact structure assembly, and then investigated the effects of the gear tooth on the dynamic responses of the DC motor associated with the kinematics of the PGT [10]. Mashimo et al. proposed a micro geared ultrasonic motor and its dynamic model was established to analyze the dynamic characteristics [11]. Kim et al. proposed a dual actuator unit (DAU) composed of two actuators and a planetary gear train; the mathematic model of the DAU was developed to investigate the performance of position tracking, stiffness variation, etc. [12]. In addition, Lee et al. also proposed a similar actuator system using a planetary gear and two motors; subsequently, a dynamic model of the actuator system was established to analyze the output velocity and torque, and obtain frequency response characteristics [13]. In addition, regarding synchronization characteristics, Wei et al. established a coupling dynamic model of a multi-source driving transmission system considering the flexible shafts, support bearing and meshing gear pair; the speed and torque synchronization properties under an impact load at different load change rates were studied [14]. Zhao et al. investigated the real-time speed synchronization of multiple induction motors during load variation and speed acceleration [15]. Xiong et al. addressed the system dynamic response synchronization and real-time problems by employing a dynamic simulation test system for electric vehicles with a dual drive circuit [16]. Shi et al. investigated the coupling characteristics and synchronization of multi-motors by using the presented synchronization control algorithm for a multi-motor speed synchronous driving system. The algorithm was employed to control the motion of the electronic circuit axis and virtual transmission shaft. The results showed that it improved the running response speed and synchronization accuracy between the motors [17]. Yang et al. established a translation-torsion coupling dynamic model of a multi-motor torque coupling system; the influence of system speed synchronization performance and torque synchronization performance on the load sharing behavior of the system were discussed [18].
At present, the related research on the combination of a single motor or multiple motors with a gear transmission system attached to a moving body mainly focuses on the dynamic response and synchronizing characteristics under inertial conditions, while related research under non-inertial conditions has rarely been reported. In fact, the gear transmission system is installed and fixed on the moving body and moves with the moving body. Obviously, under non-inertial conditions, the gear transmission system is not only affected by its own excitation forces, but is also affected by the additional inertial forces generated by the change in the operating state of the moving body. Therefore, it is necessary and meaningful to investigate the dynamic response and synchronizing characteristics of the combination of a single motor or multiple motors with a gear transmission system attached to a moving body. This paper will take the DMDS installed on the NEVs as an example to carry out the related research, and the organization of the paper is as follows. In Section 2, the dynamic model of DMDS in a non-inertial system is introduced. Then, the response and synchronizing characteristic of the DMDS under the translational motion and circular motion, two typical forms of NEV motion, are discussed in Section 3. Finally, some conclusions are drawn in Section 4.

2. Dynamic Model of Dual-Motor Driving System in Non-Inertial System

The DMDS, as shown in Figure 1, is composed of PMSMs, pinions and a bull wheel, etc. The two pinions are driven by two PMSMs (the rated power of each one is 100 kW) separately and the wheel is driven by the two pinions together; the shaft is connected to the wheel via the spline. Based on the finite element theory, the DMDS can be divided into the following four basic coupling elements: meshing element, shafting element, connecting element and bearing element. The dynamic model and node diagram of DMDS are shown in Figure 2.

2.1. Dynamic Model of Basic Coupling Elements

Meshing elements. The meshing elements can be regarded as a bridge to couple two or more rotating shafts together and transmit motion and force, and can be described by a pair of rigid disks connected by a stiffness-damping spring system on the plane of action, as shown in Figure 3. φ, β and α are the relative position angle, helix angle and pressure angle of gears, respectively. rp and rg are the base radius of gear p and gear g, respectively. km is the meshing stiffness, which can be represented by rectangular ware and regarded as the time-varying function of gear rotational angular displacement, as shown in Figure 4. em and cm are the transmission error and meshing damping, respectively, and they are the time-varying variables. Herein, em can be solved by referring to Ref. [18] and viewed as the time-varying function of gear rotational angular displacement. Therefore, the time-varying variables km, cm and em can be expressed as Equation (2). xi, yi and zi (i = p and g) are the horizontal, vertical and axial direction displacements of gears, respectively. θi (i = xp, yp and xg, yg) is the swing angular displacement of gears around the x-axis or y-axis, θj (j = zp and zg) is the rotational angular displacement of gears around the z-axis. The dynamic model of the meshing elements can be expressed as
M p g X ¨ p g + C p g X ˙ p g + K p g X p g = F p g
where
{ M p g = d i a g [ m p   m p   m p   I x p   I y p   I z p   m g   m g   m g   I x g   I y g   I z g ] ;     X p g = [ x p   y p   z p   θ x p   θ y p   θ z p   x g   y g   z g   θ x g   θ y g   θ z g ] T K p g = k m V p g T V p g ; F p g = ( k m e m + c m e ˙ m ) V p g T ; C p g = c m V p g T V p g V p g = [ cos β sin ( α φ ) ,   cos β cos ( α φ ) ,   sin β ,   r p sin β sin ( α φ ) ,   r p sin β cos ( α φ ) ,   r p cos β ,   cos β sin ( α φ ) ,   cos β cos ( α φ ) ,   sin β ,   r g sin β sin ( α φ ) ,   r g sin β cos ( α φ ) , r g cos β ]
Mpg, Cpg and Kpg are the mass, damping and stiffness matrices of the gear pair, respectively. Xpg and Vpg are the displacement vector and meshing matrix of the gear pair. Fpg is the exciting force vector, which is a non-constant vector and is related to the time-varying variables km, cm and em.
{ k m = k ¯ p g + l = 1 7 k ^ p g l cos ( l 2 π r p θ p p b p + γ i )     c m = 2 ξ k m / ( 1 / m e q , 1 + 1 / m e q , 2 )   e m = e ¯ p g + l = 1 7 e ^ p g l cos ( l 2 π r p θ p p b p + γ i )
where k ^ p g l and k ¯ p g are the lth harmonic term amplitude and mean mesh stiffness, respectively. e ^ p g l and e ¯ p g are the lth harmonic term amplitude and mean transmission error, respectively. γi is the phase angle. θp and pbp denote the rotational angular displacement around the z-axis of gear p, and the base pitch of gear p, respectively. ξ is the damping ratio. meq,1 and meq,2 are the equivalent qualities of the gears.
Shafting elements. The rotating components (such as shafts and gears) of the DMDS can be equivalent to the shafting elements with different cross sections, including hollow circle cross-section, circle cross-section, etc. The force state and deformation state of the shafting elements can be analyzed by employing the Timoshenko beam theory [19,20], as shown in Figure 5. The deformation of the shafting elements can be described by the displacement of the nodes and the dynamic model of the two-node shafting element can be express as
M i , i + 1 X ¨ i , i + 1 + C i , i + 1 X ˙ i , i + 1 + K i , i + 1 X i , i + 1 = 0
Mi,i+1 and Xi,i+1 are the mass matrix and displacement vector, respectively. Ki,i+1 is the stiffness matrix with a 12 × 12 order form and can be solved by referring to Ref. [20]. Ci,i+1 is the damping matrix and can be solved by Rayleigh damping, Ci,i+1 = λ1Mi,i+1 + λ2Ki,i+1, where λ1 and λ2 represent the proportionality coefficients.
Connecting elements. The connecting elements can be employed to describe the coupling relationship between two connecting components, such as couplings and splines. For example, as shown in Figure 6, the coupling relationship between nodes i on the internal spline and node j on the external spline can be described by means of a spring-damper system. Therefore, the dynamic model of the connecting elements can be expressed as
[ M i O O M j ] [ X ¨ i X ¨ j ] + [ C c o C c o C c o C c o ] [ X ˙ i X ˙ j ] + [ K c o K c o K c o K c o ] [ X j X ˙ j ] = 0
Mi and Mj are the mass matrices of nodes i and j, respectively. Xi and Xj represent the displacement vectors of nodes i and j, respectively. Kco and Cco are the coupling stiffness and damping matrices between the internal and external splines, respectively.
Bearing elements. The bearings are used to support the rotating shafts and constrain its horizontal and vertical movement, which can be equivalent to a spring-damper system, as shown in Figure 7. Kbi is the support stiffness matrix of the support node i on the rotating shaft, and the Ref. [21] can be referred to in order to provide the specific form of Kbi. Cbi represents the support damping matrix of the support node i on the rotating shaft, and Cbi has the same structure as Kbi. The dynamic model of the bearing element’s corresponding support node i on the rotating shaft can be expressed as
M b i X ¨ b i + C b i X ˙ b i + K b i X b i = 0
Xbi and Mbi represent the displacement vector and mass matrix of the support node i on the rotating shaft, respectively.

2.2. Dynamic Model of Dual-Motor Driving System in Non-Inertial System

The rotating shafts of the DMDS are simulated by the shafting elements, and they can be assembled together by means of the meshing elements, connecting elements and bearing elements according to the node distribution diagram illustrated in Figure 2. The coupling matrix assemble method described in the Ref. [21] can be used to establish the overall system dynamic model of DMDS, which contains 44 nodes. The overall system dynamic model can be expressed as
M X ¨ ( t ) + C ( t ) X ˙ ( t ) + K ( t ) X ( t ) = F ( t ) + G ( t ) + M ( t ) + T ( t )
M and X(t) are the overall mass matrix and displacement vector, and their dimensions are 264 × 264 and 264 × 1, respectively. C(t) and K(t) represent the overall damping and stiffness matrices, and the dimensions of these matrices are 264 × 264, respectively. F(t), G(t) and M(t) and T(t) are the additional inertial force vector, gravity and gyroscopic moment vector and external load excitation vector, and the dimensions of these vectors are 264 × 1, respectively. T(t) can be expressed as
T ( t ) = [ 0 ,   0 ,   0 ,   0 ,   0 ,   0 ,   ,   0 ,   0 ,   0 ,   0 ,   0 ,   T i n p o w e r   i n p u t   n o d e   6 ,   ,   0 ,   0 ,   0 ,   0 ,   0 ,   0 n o d e   i ,   ,   0 ,   0 ,   0 ,   0 ,   0 ,   T o u t p o w e r   o u p u t   n o d e   28 ,   ,   0 ,   0 ,   0 ,   0 ,   0 ,   T i n p o w e r   i n p u t   n o d e   35 ,   ,   0 ,   0 ,   0 ,   0 ,   0 ,   0 n o d e   N ] T
Tout is the load torque that acts on the power output node. Tin is the output torque of PMSMs that acts on the power input node. Tin can be obtained by establishing the simulation model of PMSMs, and the modeling process is described in detail in Ref. [22], and the key parameters of PMSMs are listed in Table 1.
Additional force and moment. The DMDS is installed and fixed to the bodywork and moves together with the bodywork in space. Therefore, the effects of non-inertia on the DMDS operating in moving bodywork should be considered, and the thermal changes between PMSM and the gear transmission system is ignored. As shown in Figure 8a, a non-inertial system is employed to describe the absolute motion of point Q in space, where OXYZ is a static coordinate system fixed on earth, and oxyz is a moving coordinate system that is relative to the OXYZ and fixed on moving bodywork. Ei and ei (i = 1, 2 and 3) are the base vectors of OXYZ and oxyz, respectively. ro and rQ are the position vectors of origin o and moving point Q, respectively and they are measured in OXYZ. ro/Q is the relative position vector of moving point Q, which is measured in oxyz. Ω is the angular velocity vector of oxyz and measured in OXYZ. The position vector rQ can be expressed as
r Q = r o + r o / Q = X o E 1 + Y o E 2 + Z o E 3 + x Q e 1 + y Q e 2 + z Q e 3
where Xo, Yo, Zo and xQ, yQ, zQ are the variables. Ei and ei (i = 1, 2 and 3) are the constant vector and variable vector, respectively. The absolute velocity vector vQ can be obtained by taking the derivative of the position vector rQ with respect to time, and it can be expressed as
{ v Q = d r Q d t = X ˙ o E 1 + Y ˙ o E 2 + Z ˙ o E 3 + x ˙ Q e 1 + y ˙ Q e 2 + z ˙ Q e 3 + x Q d e 1 d t + y Q d e 2 d t + z Q d e 3 d t      = v o + v o / Q + Ω × r o / Q d e i d t = Ω × e i ,   i = 1 ,   2   a n d   3
Then, the absolute acceleration vector aQ can be obtained by taking the derivative of the absolute velocity vector vQ with respect to time, and it can be expressed as
{ a Q = d v Q d t = a o + d v o / Q d t + d Ω d t × r o / Q + Ω × d r o / Q d t = a o + a o / Q + Ω ˙ × r o / Q + 2 Ω × v o / Q + Ω × ( Ω × r o / Q ) d r o / Q d t = x ˙ Q e 1 + y ˙ Q e 2 + z ˙ Q e 3 + x Q d e 1 d t + y Q d e 2 d t + z Q d e 3 d t
(1) When the NEVs perform a variable speed translational motion, from Equation (10), it can be observed that the absolute acceleration term of the moving point Q can be derived as aQ = ao + ao/Q, where ao is the convected acceleration term and ao/Q is the relative acceleration term. Since there is no relative movement between DMDS and the bodywork, this means that ao/Q = 0. Therefore, when the NEVs perform a variable speed translational motion, the additional inertial force term F(t) can be expressed as
F ( t ) = m i a i / o = m i [ a 1 x / o ,   a 1 y / o ,   a 1 z / o ,   0 ,   0 ,   0 n o d e   1 ,   ,   a i x / o ,   a i y / o ,   a i z / o ,   0 ,   0 ,   0 n o d e   i ,   ,   a N x / o ,   a N y / o ,   a N z / o ,   0 ,   0 ,   0 n o d e   N ] T
where mi is the mass of node i on the rotating shaft, and ai/o (i = 1, 2, …, N) represents the convected acceleration vector of node i on the rotating shaft and is measured in OXYZ. aix/o, aiy/o and aiz/o are the components of ai/o in x, y and z directions, respectively. Besides F(t), the node i on the rotating shaft is also affected by gravity G(t), and it can be expressed as
G ( t ) = m i a i g = m i [ a 1 g x ,   a 1 g y ,   a 1 g z ,   0 ,   0 ,   0 n o d e   1 ,   ,   a i g x ,   a i g y ,   a i g z ,   0 ,   0 ,   0 n o d e   i ,   ,   a N g x ,   a N g y ,   a N g z ,   0 ,   0 ,   0 n o d e   N ] T
where aig (i = 1, 2, …, N) represents the gravitational acceleration vector of node i on the rotating shaft. aigx, aigy and aigz are the components of aig in x, y and z directions, respectively.
(2) When the NEVs performs a variable speed circular motion, from Equation (10), the absolute acceleration term of moving point Q can be derived as follows: aQ = ao/Q + ×ro/Q + Ω × (Ω × ro/Q) + 2Ω × vo/Q, where Ω ˙ × ro/Q and Ω × (Ω × ro/Q) are the first convected acceleration term and second convected acceleration term, respectively, and are denoted as ar1 and ar2, respectively. Herein, ar = ar1 + ar2. 2Ω × vo/Q is the Coriolis acceleration term and is denoted as ac. Similarly, since there is no relative movement between DMDS and the bodywork, the relative acceleration term ao/Q = 0. Thus, when the NEVs performs a variable speed circular motion, the additional inertial force term F(t) can be expressed as
F ( t ) = m i ( a i r + a i c ) = m i [ a 1 r x + a 1 c x ,   a 1 r y + a 1 c y ,   a 1 r z + a 1 c z ,   0 ,   0 ,   0 n o d e   1 ,   ,     a i r x + a i c x ,   a i r y + a i c y ,   a i r z + a i c z ,   0 ,   0 ,   0 n o d e   i ,   ,   a N r x + a N c x ,   a N r y + a N c y ,   a N r z + a N c z ,   0 ,   0 ,   0 n o d e   N ] T
where air and aic (i = 1, 2, …, N) represent the convected acceleration and Coriolis acceleration vectors of node i on the rotating shaft, respectively. airx, airy, airz and aicx, aicy, aicz are the components of air and aic in x, y and z directions, respectively. The gravity term G(t) can be calculated by Equation (12). Moreover, the rotating shaft will be subjected to an additional gyroscopic moment, due to the rotating angular velocity ω when the NEVs perform a circular motion. As shown in Figure 8b, the moment of momentum can be described as H = J·ω and J is the polar moment of inertia. Hence, under the interaction of Ω and ω, the additional gyroscopic moment M(t) can be expressed as
M ( t ) = H × Ω = J ω × Ω = [ 0 ,   0 ,   0 ,   J 1 x ω 1 x Ω 1 x ,   J 1 y ω 1 y Ω 1 y ,   J 1 z ω 1 z Ω 1 z n o d e   1 ,   ,       0 ,   0 ,   0 ,   J i x ω i x Ω i x ,   J i y ω i y Ω i y ,   J i z ω i z Ω i z n o d e   i ,   ,   0 ,   0 ,   0 ,   J N x ω N x Ω N x ,   J N y ω N y Ω N y ,   J N z ω N z Ω N z n o d e   N ] T
Jij and ωij (j = x, y and z) are the moment of inertia and angular velocity of node i around the j-axis, respectively, and Ωij (j = x, y and z) is the angular velocity of node i around the axes other than the j-axis.

3. Dynamic Response and Synchronizing Characteristic Analysis in Non-Inertial System

The software MATLAB/Simulink (license: 40927247) is used to build the simulation model of the DMDS, and the translational motion and circular motion of the NEVs are selected to investigate the response and synchronizing characteristics of the DMDS in a non-inertial system. The key parameters of the DMDS are listed in Table 2.

3.1. Case I: The NEVs Perform a Variable Speed Translational Motion

If one supposes that the NEVs travel in a straight line with variable speed in the x-direction, the force condition is shown as Figure 9, and the basic parameters of NEVs are listed in Table 3. For the gravity of node j, Gj = −mj·ag·(0·e1 + e2 + 0·e3), as shown in Figure 9a. For the additional convected inertial force (ACIF) of node j, Fj = −mj·ajx·(e1 + 0·e2 + 0·e3), as shown in Figure 9b. The total gravity and ACIF of a shaft are equal to the sum of Gj and Fj of all the nodes on the shaft, respectively.
The curves of the translation acceleration ax and translational speed v in the x-direction shown in Figure 10 are employed to investigate the response and synchronizing characteristics of the DMDS under different non-inertial conditions. The initial translational speed of the vehicle is 80 km/h and the rotating speeds of PMSMs are 1900 rpm, and the load imposed on the DMDS is 600 Nm. If one takes as an example the displacement responses of node 14 (gear node) under different non-inertial conditions, as shown in Figure 11, it can be clearly observed that the displacement response curves in the x- and y-directions are significantly associated with the ax-t curve; the difference between both of them is that they exhibit a full reversal. The displacement responses increase with the increase in ax, but the magnitude in the y-direction is markedly greater than that in the x-direction.
According to the theoretical analysis, when the NEVs perform a translational motion, the nodes are not only affected by the torque of PMSMs (Tin), but also by ACIF (F) in the x-direction and gravity (G) in the y-direction. Figure 12a shows the displacement responses with or without F when ax changes from 10 to −10 m·s−2. As a consequence, the offset of the red line is significantly greater than that of blue line. Similarly, due to the effect of gravity G, the offset of the displacement response in the y-direction has increased, but is not obvious, as shown in Figure 12b.
In addition, the relative influence proportion γi is used to analyze the degree of influence of F or G on the displacement response, and it can be calculated by Equation (15). As shown in Figure 13a, the influence of F on the displacement response in the x-direction is obviously greater than that of G on the displacement response in the y-direction. When ax changes from −15 to 15 m·s−2, the relative influence proportion of F has a numerically negative correlation with ax, and then has a numerically positive correlation with ax. When ax = −15 or 15 m·s−2, the relative influence proportion reaches the maximum value, which is about 24%. However, when ax changes from −15 to 15 m·s−2, the relative influence proportion of G has a numerically positive correlation with ax, and then has a numerically negative correlation with ax, and the maximum relative influence proportion only reaches 2.3% when ax = 0 m·s−2, as shown in Figure 13b. Therefore, the gravity G is not the primary influence factor, and its influence can be ignored whether it is in inertial or non-inertial conditions.
γ i = ( δ i δ t ) / ( δ i δ t ) × 100 %
δi is the offset when the factor i is not considered and δt is the offset when all factors are considered.
The vibration of gears can be transmitted to the PMSMs and affect the torque and speed of PMSMs. Considering that the gravity G has very little effect on the displacement responses of gears, we next focus on the effect of ACIF on the torque and speed of PMSMs. Herein, the synchronization error (SE) is applied to measure the synchronization of the output variables (torque and speed) between PMSMs. The greater the SE, the greater the difference in the output variables between PMSMs, which can be expressed as
E λ = ( λ l λ m e a n ) max λ m e a n × 100 % , λ m e a n = 1 N i = 1 N λ l
Eλ is the SE, λ = Tin and v, l = 1, 2, …, N. Tin and v represent the torque and speed of PMSMs, respectively. N is the number of PMSMs and N = 2 in this study.
The torque of PMSMs and its torque synchronization error (TSE) are shown in Figure 14. Except for the fact that there is some fluctuation at 4.0–4.5 s and 5.0–5.5 s, the curve of torque is very close to the axt curve, as shown in Figure 14a. Obviously, the torque remains almost unchanged before and after considering the ACIF. Similarly, the TSE between PMSMs has hardly changed, and their amplitudes fluctuate between −2.1% and 2.3% during acceleration and fluctuate between −14.2% and 18.8% during deceleration. This phenomenon is firstly because the growth of the numerator in Equation (16) is lesser than that of the denominator, and second because the mean value of torque (namely, λmean) during acceleration is significantly greater than that during deceleration. Figure 14c shows the curve of the root mean square (RMS) value of TSE. When ax changes in the range of −15 to 15 m·s−2, with the decrease in the numerical value of ax, the RMS value of TSE increases and then decreases and the maximum value appears when ax = −2.5 m·s−2 and is about 2.8%. The reason for this is that the torque of PMSMs is at its minimum value when ax = −2.5 m·s−2, namely when λmean in Equation (16) is at its minimum value, so the calculated Eλ is at its maximum value. It indicates that the change in translation acceleration has a small effect on the synchronization of torque between PMSMs. Moreover, through the analysis and comparison, it is found that the change in the RMS value of TSE is minor before and after considering the ACIF, and the maximum difference between the two is only 0.26%, indicating that the ACIF also has a small effect on the synchronization of torque between PMSMs.
The speed of PMSM and its speed synchronization error (SSE) are shown in Figure 15. The curve of speed exhibits good correlation with the vt curve, as shown in Figure 15a. There is no obvious change in speed at any of the times, regardless of whether the ACIF is considered or not. This can also be verified by Figure 15b and compared with the case without considering F, the change in SSE is not noticeable, and its amplitude fluctuates between −2.2‰ and 2.3‰. In addition, the curve of the RMS value of SSE under different translation acceleration conditions is plotted, as shown in Figure 15c. When ax changes in the range of −15 to 15 m·s−2, with the decrease in the numerical value of ax, the RMS value of SSE increases and then decreases and the maximum value appears when ax = −0 m·s−2 and is about 5.03‰. This is because the speed of PMSMs is at its minimum value when ax = 0 m·s−2, namely when λmean in Equation (16) is at its minimum value, so the calculated Eλ is at its maximum value. This means that the changes in ax have a small effect on speed synchronization between PMSMs. Meanwhile, it can also be clearly observed that the increase in the RMS value of SSE is not obvious after considering ACIF, and the maximum difference between the two is only 0.27‰, indicating that the ACIF has a very small influence on speed synchronization between PMSMs.

3.2. Case II: The NEVs Perform a Variable Speed Circular Motion

If one supposes that the NEVs perform a variable-speed circular motion with radius RΩ around the Y-axis, θY is the angle displacement and ΩY is the rotating velocity of the vehicle around the Y-axis and they are measured in OXYZ, as shown in Figure 16. The gravity Gj in case II is calculated in the same way as that in case I. For the additional inertial force of node j, Fj = −mj·(ajr +ajc); ajr and ajc are the convected acceleration vector and Coriolis acceleration vector of node j, respectively. For the additional gyroscopic moment of node j, Mj = Jj·ω × ΩY.
As shown in Figure 17a, the convected acceleration vector ajr = Ω ˙ Y × ro/j + ΩY × (ΩY × ro/j) = Ω ˙ Y × (rj/z + rj/x + rj/y) + ΩY × (ΩY × (rj/z + rj/x + rj/y)) = Ω ˙ Y × (rj/z + rj/x) + ΩY × (ΩY × (rj/z + rj/x)). One must suppose that ajr1 = Ω ˙ Y × (rj/z + rj/x) is the additional first convected acceleration vector and ajr2 = ΩY × (ΩY × (rj/z + rj/x)) is the additional second convected acceleration vector. Hence, the complex three-dimension spatial problem can be solved in a two-dimensional plane, i.e., the xoz plane. As shown in Figure 17b, (m, n) and (m’, n’) represent the initial position coordinate and instantaneous position coordinate of node j in the xoz plane and x1o1z1 plane, respectively, and the convected acceleration vector ajr can be derived as
{ a j r 1 = | Ω ˙ Y × r o / j | cos ( γ j ) e 1 | Ω ˙ Y × r o / j | sin ( γ j ) e 3 ; a j r 2 = | Ω Y × ( Ω Y × r o / j ) | sin ( γ j ) e 1 | Ω Y × ( Ω Y × r o / j ) | cos ( γ j ) e 3 | Ω ˙ Y × r o / j | = | Ω ˙ Y | ( m + m ) 2 + ( n + n ) 2 ;   | Ω Y × ( Ω Y × r o / j ) | = | Ω Y | 2 ( m + m ) 2 + ( n + n ) 2 ;   γ j = arctan ( ( n + n ) / ( m + m ) )
In Figure 17b, vj is the velocity vector of node j in the xoz plane, and the Coriolis acceleration vector ajc can be derived as
{ v j = v j x e 1 + v j z e 3 a j c = 2 Ω Y × v j = 2 Ω Y × ( v j x e 1 + v j z e 3 ) = 2 Ω Y v j x e 2 × e 1 + 2 Ω Y v j z e 2 × e 3 = 2 Ω Y v j z e 1 2 Ω Y v j x e 3
where ΩY is an angular velocity vector in the y-axis, ΩY is an angular velocity of node i around the y-axes and it is a number.
The initial positions and instantaneous positions of nodes i and j are shown in Figure 17c and Lij represents the initial axial distance between i and j. ( x i , y i , z i ) and ( x j , y j , z j ) are the instantaneous position coordinates of nodes i and j. According to the Ref. [23], the additional gyroscopic moment M(t) of any node can be derived as Equation (19), and the terms with a value of zero are not given in Equation (19).
M ( t ) = J ω Ω Y [ L i j + ( z j z i ) ] ( x j x i ) 2 + ( y j y i ) 2 + [ L i j + ( z j z i ) ] 2 e 1 + J ω Ω Y ( x j x i ) ( x j x i ) 2 + ( y j y i ) 2 + [ L i j + ( z j z i ) ] 2 e 3
The Ω ˙ Y t, ΩYt and θYt curves shown in Figure 18 are used to study the response and synchronizing characteristics of the DMDS when the NEVs perform a variable speed circular motion. Similarly, the initial speed of the vehicle, rotating speeds of PMSMs and load imposed on the DMDS are also 80 km/h, 1900 rpm and 600 Nm, respectively. The turning radius RΩ is set to 5 m, and the rotating acceleration Ω ˙ Y is set to 5, 10 and 15 m·s−2, respectively. The displacement responses of node 14 (gear node) are shown in Figure 19. The displacement response curves in the y- and z-directions are significantly associated with the Ω ˙ Y t and ΩYt curves, respectively, and the difference is that the directions of the displacement response curves is completely opposite to that of the Ω ˙ Y t and ΩYt curves. However, due to the coupling effect of Ω ˙ Y and ΩY, the displacement response curves in the x-direction are neither completely related to Ω ˙ Y t nor ΩYt. The displacement responses in three directions increase with the increase in Ω ˙ Y , but the increases are different.
Based on the above theoretical analysis, under the given non-inertial condition, except for the torque of PMSM (Tin), the nodes are also affected by the following five additional terms: first convected inertial force Fr1, second convected inertial force Fr2, Coriolis inertial force Fc, gyroscopic moment M and gravity G. The displacement responses of node 14 with or without additional inertial terms (AITs) are shown in Figure 20. In the x-direction, after considering AIFs, the displacement increases gradually during the II and III stages, and then decreases during the IV stage, and eventually reaches stabilization during the V stage. The displacement increases from about 0.62 µm to 1.18 µm, an increase of 90.3%, as shown in Figure 20a. In the y-direction, the displacement shows a small increase during the II-V stages, and increases from about 25 µm to 28 µm with an increase of 12% during the V stage. In the z-direction, the difference by not considering AITs is that the displacement increases gradually during the II-IV stages and then tends to become steady during the V stage. Finally, the displacement increases from approximately 3.6 µm to 7.5 µm, with an increase of 108.3%.
Herein, the Equation (14) is also used to calculate the relative influence proportion (γi) of the five additional terms and to analyze which terms are the major influence factors and which are the weakest ones, as shown in Figure 21. As demonstrated by the results in Figure 21a, during the I stage, the relative influence proportions of Fr1, Fc and M are greater, reaching 2.9–76%, 9.1–48.6% and 10.6–48.6%, respectively. Then, their relative influence proportions decrease gradually during the II-V stages. However, during the II-V stages, the relative influence proportion of Fr2 starts to increase quickly and becomes the main influence factor, and the maximum proportion can reach 93.9%. In Figure 21b, during the I stage, the relative proportions of G, Fr1, Fc and M are greater, and there is a little difference in their average proportions. Conversely, the relative influence proportion of M tends to become steady after the increase during the II-V stages, the maximum proportion reaches 99.3%, and the relative influence proportions of other additional terms are very small and their influences can be ignored. By comparing Figure 21a,b, it is found that the relative influence proportions of the five additional terms have a similar variation rule, but the shapes of the relative influence proportions of Fc and M in the two figures are exactly opposite. Likewise, the proportions of Fr1, Fc and M are still greater during the I stage and then begin to decrease rapidly during the II-V stages. When it comes to stage II, the proportion of Fr2 begins to increase quickly and then becomes the primary influence factor.
The torque and RMS value of TSE of the PMSMs are shown in Figure 22. Obviously, the curves of torque are very close to the Ω ˙ Y t curve, and increase with the increase in rotating acceleration, as shown in Figure 22a. The results in Figure 22b show that the curves of RMS value of TSE are associated with the Ω ˙ Y t curve, and exhibit a full reversal. Unlike torque, the RMS value of TSE decreases with the increase in rotating acceleration. When Ω ˙ Y changes from 5 to 15 m·s−2, the RMS value of TSE is only reduced from about 1.0% to 0.9%, and the reduction is no more than 10%, indicating that the rotating acceleration also has a small influence on the synchronization of torque between PMSMs. Moreover, the RMS values of TSE with or without AITs are shown in Figure 22c and it can be very clearly observed that the red solid line almost coincides with the blue solid line at all times whether the AITs are considered or not, and the maximum difference does not exceed 0.05%, which means that the influence of AITs on the synchronization of torque between the PMSMs can almost be ignored when the NEVs perform a circular motion.
The speed and RMS value of SSE of the PMSMs are shown in Figure 23. The curves of speed and RMS values of SSE are very close to the Ω Y t curve, but the curves of RMS values of SSE and Ω Y t have a diametrically opposite changing tendency. The speed increases as the rotating acceleration increases; conversely, the RMS value of SSE decreases as the rotating acceleration increases, with a maximum reduction of no more than 1.5‰, as shown in Figure 23a,b. Combining the results from Figure 23c, it is found that the influence of both AITs and rotating acceleration on the synchronization of speed between PMSMs can also be ignored when the vehicle performs a circular motion.

4. Conclusions

The dynamic model of the DMDS installed on the NEVs was established considering the effects of the non-inertial system in the form of additional inertial terms. The response and synchronizing characteristics of the DMDS under different non-inertial conditions were investigated and compared, and the following conclusions were obtained:
(1) When the NEVs perform a variable speed translational motion, in addition to its own excitation forces and gravity, the components of the DMDS are subjected to an additional convected inertial force caused by the variable speed translational motion. There is a numerically positive correlation between the displacement response of the gears and the translation acceleration, and the contribution of the additional convected inertial force to displacement offset in the x-direction is greater than that of gravity in the y-direction.
(2) Due to the influence of multiple additional inertial terms’ superimposition, the variable speed circular motion of the NEVs results in greater displacement responses. For the displacement responses in the x and z directions, the additional second relative inertial term (namely, the Fr2 term) is the main influence factor among all the additional inertial terms that change the displacement offset. In the y-direction, the additional gyroscopic moment term (namely, the M term) is the dominant influence factor, and its influence is more prominent under high-speed operating conditions.
(3) Regardless of whether the NEVs perform a variable speed translation motion or variable speed circular motion, the acceleration (translation acceleration or rotating acceleration) has a significant influence on the torques and speeds of PMSMs, and they increase with the increase in the acceleration. However, the torque synchronization and speed synchronization between PMSMs is minimally influenced by both the additional inertial terms and the acceleration.

Author Contributions

Writing—review and editing, Z.X.; Writing—original draft, Investigation, R.S.; Conceptualization, J.H.; Methodology, B.F.; Validation, Z.Z.; Data curation, R.T. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (No. 51905060/51905062/51905064/51805060/51875068); China Postdoctoral Science Foundation (No. 2021M700619); Natural Science Foundation Project of Chongqing Science and Technology Commission (No. cstc2020jcyj-msxmX0346/cstc2021jcyj-msxmX0361); Science and Technology Research Program of Chongqing Municipal Education Commission (No. KJQN201901107); National Key Research and Development Project of China (No. 2019YFB1703700). We also want to thank the support from the Scientific Research Foundation of Chongqing University of Technology and the State Key Laboratory of Mechanical Transmissions of Chongqing University, China.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no conflict of interest to report regarding the present study.

References

  1. Yang, H.; Zhang, J.; Wen, X. High Power Dual Motor Drive System Used in Fuel Cell Vehicles. In Proceedings of the IEEE Vehicle Power and Propulsion Conference, Harbin, China, 3–5 September 2008. [Google Scholar]
  2. Tang, Y. Dual Motor Drive and Control System for an Electric Vehicle. U.S. Patent 8453770B2, 4 June 2013. [Google Scholar]
  3. Inalpolat, M.; Handschuh, M.; Kahraman, A. Influence of indexing errors on dynamic response of spur gear pairs. Mech. Syst. Signal Process. 2015, 60, 391–405. [Google Scholar] [CrossRef]
  4. Giagopulos, D.; Salpistis, C.; Natsiavas, S. Effect of non-linearities in the identification and fault detection of gear-pair systems. Int. J. Non-Linear Mech. 2006, 41, 213–230. [Google Scholar] [CrossRef]
  5. Setiawan, Y.D.; Roozegar, M.; Zou, T.; Angeles, J. A Mathematical Model of Multi-Speed Transmissions in Electric Vehicles in the Presence of Gear-Shifting. IEEE Trans. Veh. Technol. 2018, 67, 397–408. [Google Scholar] [CrossRef]
  6. Sun, W.; Ding, X.; Wei, J.; Wang, X.; Zhang, A. Hierarchical modeling method and dynamic characteristics of cutter head driving system in tunneling boring machine. Tunn. Undergr. Space Technol. 2016, 52, 99–110. [Google Scholar] [CrossRef]
  7. Liu, C.; Fang, Z.D.; Liu, X.; Hu, S.Y. Multibody dynamic analysis of a gear transmission system in electric vehicle using hybrid user-defined elements. Proc. Inst. Mech. Eng. Part K J. Multi-Body Dyn. 2018, 233, 30–42. [Google Scholar] [CrossRef]
  8. Tang, X.; Hu, X.; Yang, W.; Yu, H. Novel Torsional Vibration Modeling and Assessment of a Power-Split Hybrid Electric Vehicle Equipped with a Dual Mass Flywheel. IEEE Trans. Veh. Technol. 2018, 67, 1990–2000. [Google Scholar] [CrossRef]
  9. Zhang, Y.; Spanos, P.D. Efficient response determination of a M-D-O-F gear model subject to combined periodic and stochastic excitations. Int. J. Non-Linear Mech. 2019, 120, 103378. [Google Scholar] [CrossRef]
  10. Yan, H.S.; Wu, Y.C. A Novel Design of a Brushless DC Motor Integrated with an Embedded Planetary Gear Train. IEEE/ASME Trans. Mechatron. 2005, 11, 551–557. [Google Scholar] [CrossRef]
  11. Mashimo, T.; Urakubo, T.; Shimizu, Y. Micro Geared Ultrasonic Motor. IEEE/ASME Trans. Mechatron. 2018, 23, 781–787. [Google Scholar] [CrossRef]
  12. Kim, B.S.; Song, J.B.; Park, J.J. A Serial-Type Dual Actuator Unit With Planetary Gear Train: Basic Design and Applications. IEEE/ASME Trans. Mechatron. 2009, 15, 108–116. [Google Scholar]
  13. Lee, H.; Choi, Y. A New Actuator System Using Dual-Motors and a Planetary Gear. IEEE/ASME Trans. Mechatron. 2012, 17, 192–197. [Google Scholar] [CrossRef]
  14. Wei, J.; Shu, R.; Qin, D.; Lim, T.C.; Zhang, A.; Meng, F. Study of synchronization characteristics of a multi-source driving transmission system under an impact load. Int. J. Precis. Eng. Manuf. 2016, 17, 1157–1174. [Google Scholar] [CrossRef]
  15. Zhao, D.Z.; Li, C.W.; Ren, J. Speed synchronisation of multiple induction motors with adjacent cross-coupling control. Iet Control Theory Appl. 2009, 4, 119–128. [Google Scholar] [CrossRef]
  16. Xiong, H.; Zhang, M.; Zhang, R.; Zhu, X.; Yang, L.; Guo, X.; Cai, B. A new synchronous control method for dual motor electric vehicle based on cognitive-inspired and intelligent interaction. Future Gener. Comput. Syst. 2019, 94, 536–548. [Google Scholar] [CrossRef]
  17. Shi, T.; Liu, H.; Geng, Q.; Xia, C. Improved Relative Coupling Control Structure for Multi-motor Speed Synchronous Driving System. IET Electr. Power Appl. 2016, 10, 451–457. [Google Scholar]
  18. Yang, Y.; Li, M.; Hu, M.; Qin, D. Influence of controllable parameters on load sharing behavior of torque coupling gear set. Mech. Mach. Theory 2018, 121, 286–298. [Google Scholar] [CrossRef]
  19. Kahrobaiyan, M.H.; Asghari, M.; Ahmadian, M.T. A Timoshenko beam element based on the modified couple stress theory. Int. J. Mech. Sci. 2014, 79, 75–83. [Google Scholar] [CrossRef]
  20. Zhang, A.; Wei, J.; Qin, D.; Hou, S.; Lim, T.C. Coupled Dynamic Characteristics of Wind Turbine Gearbox Driven by Ring Gear Considering Gravity. J. Dyn. Syst. Meas. Control 2018, 140, 1–15. [Google Scholar] [CrossRef]
  21. Hu, Z.; Tang, J.; Zhong, J.; Chen, S. Frequency spectrum and vibration analysis of high speed gear-rotor system with tooth root crack considering transmission error excitation. Eng. Fail. Anal. 2016, 60, 405–441. [Google Scholar] [CrossRef]
  22. Lu, Z.; Sheng, H.; Hess, H.L.; Buck, K.M. The modeling and simulation of a permanent magnet synchronous motor with Direct Torque Control based on Matlab/Simulink. In Proceedings of the Electric Machines and Drives, San Antonio, TX, USA, 15 May 2005; pp. 1150–1156. [Google Scholar]
  23. Zhang, A.; Wei, J.; Shi, L.; Qin, D.; Lim, T.C. Modeling and dynamic response of parallel shaft gear transmission in non-inertial system. Nonlinear Dyn. 2019, 98, 997–1017. [Google Scholar] [CrossRef]
Figure 1. Structure of DMDS.
Figure 1. Structure of DMDS.
Machines 10 00620 g001
Figure 2. Dynamic model and node distribution diagram of DMDS.
Figure 2. Dynamic model and node distribution diagram of DMDS.
Machines 10 00620 g002
Figure 3. Equivalent model of meshing element.
Figure 3. Equivalent model of meshing element.
Machines 10 00620 g003
Figure 4. Meshing stiffness.
Figure 4. Meshing stiffness.
Machines 10 00620 g004
Figure 5. Timoshenko beam.
Figure 5. Timoshenko beam.
Machines 10 00620 g005
Figure 6. Equivalent model of connecting elements.
Figure 6. Equivalent model of connecting elements.
Machines 10 00620 g006
Figure 7. Bearing element.
Figure 7. Bearing element.
Machines 10 00620 g007
Figure 8. (a) Kinematic analysis in non-inertial system; (b) analysis of additional gyroscopic moment.
Figure 8. (a) Kinematic analysis in non-inertial system; (b) analysis of additional gyroscopic moment.
Machines 10 00620 g008
Figure 9. Force analysis diagram: (a) force analysis when NEV travels in a straight line; (b) force analysis of node j.
Figure 9. Force analysis diagram: (a) force analysis when NEV travels in a straight line; (b) force analysis of node j.
Machines 10 00620 g009
Figure 10. The curves of translation acceleration ax and translational speed v.
Figure 10. The curves of translation acceleration ax and translational speed v.
Machines 10 00620 g010
Figure 11. Displacement responses of node 14 under different non-inertial conditions: (a) x-direction; (b) y-direction.
Figure 11. Displacement responses of node 14 under different non-inertial conditions: (a) x-direction; (b) y-direction.
Machines 10 00620 g011
Figure 12. Displacement responses with or without additional inertial forces of node 14 when ax changes from 10 to −10 m·s−2: (a) x-direction; (b) y-direction.
Figure 12. Displacement responses with or without additional inertial forces of node 14 when ax changes from 10 to −10 m·s−2: (a) x-direction; (b) y-direction.
Machines 10 00620 g012
Figure 13. Relative proportion of each influence factor: (a) displacement response in x-direction; (b) displacement response in y-direction.
Figure 13. Relative proportion of each influence factor: (a) displacement response in x-direction; (b) displacement response in y-direction.
Machines 10 00620 g013
Figure 14. Torque and TSE when ax changes from 10 to −10 m·s−2: (a) torque of PMSMs; (b) TSE; (c) the changes in RMS of TSE with ax.
Figure 14. Torque and TSE when ax changes from 10 to −10 m·s−2: (a) torque of PMSMs; (b) TSE; (c) the changes in RMS of TSE with ax.
Machines 10 00620 g014
Figure 15. Speed and SSE when ax changes from 10 to −10 m·s−2: (a) speed of PMSMs; (b) SSE; (c) the changes in RMS of SSE with ax.
Figure 15. Speed and SSE when ax changes from 10 to −10 m·s−2: (a) speed of PMSMs; (b) SSE; (c) the changes in RMS of SSE with ax.
Machines 10 00620 g015
Figure 16. Schematic diagram of circular motion.
Figure 16. Schematic diagram of circular motion.
Machines 10 00620 g016
Figure 17. Relative motion: (a) displacement vector of node j; (b) convected acceleration vector and Coriolis acceleration vector; (c) rotational angular velocity.
Figure 17. Relative motion: (a) displacement vector of node j; (b) convected acceleration vector and Coriolis acceleration vector; (c) rotational angular velocity.
Machines 10 00620 g017
Figure 18. The curves of rotating acceleration Ω ˙ Y , rotating speed ΩY and angle displacement θY in variable-speed circular motion.
Figure 18. The curves of rotating acceleration Ω ˙ Y , rotating speed ΩY and angle displacement θY in variable-speed circular motion.
Machines 10 00620 g018
Figure 19. Displacement responses of node 14 under different rotating accelerations: (a) x-direction; (b) y-direction; (c) z-direction.
Figure 19. Displacement responses of node 14 under different rotating accelerations: (a) x-direction; (b) y-direction; (c) z-direction.
Machines 10 00620 g019
Figure 20. Displacement responses with or without AIF of node 14 when R = 5 m and Ω ˙ Y = 20 m·s−2: (a) x-direction; (b) y-direction; (c) z-direction.
Figure 20. Displacement responses with or without AIF of node 14 when R = 5 m and Ω ˙ Y = 20 m·s−2: (a) x-direction; (b) y-direction; (c) z-direction.
Machines 10 00620 g020
Figure 21. Relative influence proportion of each factor: (a) displacement response in x-direction; (b) displacement response in y-direction; (c) displacement response in z-direction.
Figure 21. Relative influence proportion of each factor: (a) displacement response in x-direction; (b) displacement response in y-direction; (c) displacement response in z-direction.
Machines 10 00620 g021
Figure 22. Torque and RMS value of TSE: (a) torque of PMSMs under different rotating accelerations; (b) RMS values of TSE under different rotating accelerations; (c) RMS values of TSE with or without AITs.
Figure 22. Torque and RMS value of TSE: (a) torque of PMSMs under different rotating accelerations; (b) RMS values of TSE under different rotating accelerations; (c) RMS values of TSE with or without AITs.
Machines 10 00620 g022
Figure 23. Speed and RMS value of SSE: (a) speed of PMSM under different rotating accelerations; (b) RMS values of SSE under different rotating accelerations; (c) RMS values of SSE with or without AITs.
Figure 23. Speed and RMS value of SSE: (a) speed of PMSM under different rotating accelerations; (b) RMS values of SSE under different rotating accelerations; (c) RMS values of SSE with or without AITs.
Machines 10 00620 g023
Table 1. Key parameters of PMSMs.
Table 1. Key parameters of PMSMs.
StatorRotorUnit
Resistance8.296 × 10−37.11 × 10−2Ω
Rated speed2400rpm
Rated voltage540VDC
Leakage inductance5.05 × 10−55.05 × 10−5H
Rotational inertia8.86 × 10−1kg·m2
Magnetizing inductance8.49 × 10−58.49 × 10−5H
Dimensionsφ380 × L399mm
Mass130kg
Rated power100kW
Rated torque400Nm
Table 2. Key parameters of DMDS.
Table 2. Key parameters of DMDS.
ParametersValuesParametersValues
Tooth number (pinion/wheel)23/69Length of shaft 1 (or 3) L1 (mm)32
Module (mm)4Diameter of shaft 2 D/d (mm)65/45
Width (mm)35Length of shaft 2 L2 (mm)20
Pressure angle α (°)20Bearing translational stiffness (Nm−1)kx = ky = 1.2 × 108; kz = 9.5 × 107
Helix angle β (°)15Bearing torsional stiffness (Nm−1)kθx = kθy = 1.0 × 107
Diameter of shaft 1 (or 3) D/d (mm)60/45Transmission error amplitude (µm)20
Table 3. Basic parameters of NEVs.
Table 3. Basic parameters of NEVs.
ParametersValuesParametersValuesParametersValues
Length (mm)4675High (mm)1500Front/rear track (mm)1525/1520
Width (mm)1770Wheelbase (mm)2670Weight (kg)1531
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xie, Z.; Shu, R.; Huang, J.; Fu, B.; Zou, Z.; Tan, R. Dynamic Response and Synchronizing Characteristic for the Dual-Motor Driving System in Non-Inertial System. Machines 2022, 10, 620. https://doi.org/10.3390/machines10080620

AMA Style

Xie Z, Shu R, Huang J, Fu B, Zou Z, Tan R. Dynamic Response and Synchronizing Characteristic for the Dual-Motor Driving System in Non-Inertial System. Machines. 2022; 10(8):620. https://doi.org/10.3390/machines10080620

Chicago/Turabian Style

Xie, Zhengqiu, Ruizhi Shu, Jin Huang, Benyuan Fu, Zheng Zou, and Rulong Tan. 2022. "Dynamic Response and Synchronizing Characteristic for the Dual-Motor Driving System in Non-Inertial System" Machines 10, no. 8: 620. https://doi.org/10.3390/machines10080620

APA Style

Xie, Z., Shu, R., Huang, J., Fu, B., Zou, Z., & Tan, R. (2022). Dynamic Response and Synchronizing Characteristic for the Dual-Motor Driving System in Non-Inertial System. Machines, 10(8), 620. https://doi.org/10.3390/machines10080620

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