Dynamic interactions of an integrated vehicle–electromagnetic energy harvester–tire system subject to uneven road excitations

Acta Mechanica Sinica, Mar 2017

An investigation is undertaken of an integrated mechanical-electromagnetic coupling system consisting of a rigid vehicle with heave, roll, and pitch motions, four electromagnetic energy harvesters and four tires subject to uneven road excitations in order to improve the passengers’ riding comfort and harvest the lost engine energy due to uneven roads. Following the derived mathematical formulations and the proposed solution approaches, the numerical simulations of this interaction system subject to a continuous sinusoidal road excitation and a single ramp impact are completed. The simulation results are presented as the dynamic response curves in the forms of the frequency spectrum and the time history, which reveals the complex interaction characteristics of the system for vibration reductions and energy harvesting performance. It has addressed the coupling effects on the dynamic characteristics of the integrated system caused by: (1) the natural modes and frequencies of the vehicle; (2) the vehicle rolling and pitching motions; (3) different road excitations on four wheels; (4) the time delay of a road ramp to impact both the front and rear wheels, etc., which cannot be tackled by an often used quarter vehicle model. The guidelines for engineering applications are given. The developed coupling model and the revealed concept provide a means with analysis idea to investigate the details of four energy harvester motions for electromagnetic suspension designs in order to replace the current passive vehicle isolators and to harvest the lost engine energy. Potential further research directions are suggested for readers to consider in the future.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:


Dynamic interactions of an integrated vehicle–electromagnetic energy harvester–tire system subject to uneven road excitations

Acta Mech. Sin. ( Dynamic interactions of an integrated vehicle-electromagnetic energy harvester-tire system subject to uneven road excitations Jing Tang Xing 0 1 Zhe Sun 0 1 Sulian Zhou 0 1 Mingyi Tan 0 1 0 College of Shipbuilding Engineering, Harbin Engineering University , Harbin 150001 , China 1 Fluid Structure Interactions Group, FEE, University of Southampton , Southampton SO17 1BJ , UK An investigation is undertaken of an integrated mechanical-electromagnetic coupling system consisting of a rigid vehicle with heave, roll, and pitch motions, four electromagnetic energy harvesters and four tires subject to uneven road excitations in order to improve the passengers' riding comfort and harvest the lost engine energy due to uneven roads. Following the derived mathematical formulations and the proposed solution approaches, the numerical simulations of this interaction system subject to a continuous sinusoidal road excitation and a single ramp impact are completed. The simulation results are presented as the dynamic response curves in the forms of the frequency spectrum and the time history, which reveals the complex interaction characteristics of the system for vibration reductions and energy harvesting performance. It has addressed the coupling effects on the dynamic characteristics of the integrated system caused by: (1) the natural modes and frequencies of the vehicle; (2) the vehicle rolling and pitching motions; (3) different road excitations on four wheels; (4) the time delay of a road ramp to impact both the front and rear wheels, etc., which cannot be tackled by an often used quarter vehicle model. The guidelines for engineering applications are given. The developed coupling model and the revealed concept provide a means with analysis idea to investigate the details of four energy harvester motions for electromagnetic suspension designs in order to replace the current passive vehicle isolators and to harvest the lost engine energy. Potential further research directions are suggested for readers to consider in the future. Vibration-energy-harvesters; Electromagnetic suspensions; Mechanical electromagnetic interactions; Vehicle dynamics; Vibration isolations 1 Introduction Environmental concerns and fuel price increases have accelerated the global trend toward the following two important research directions. One is to seek new green energy resources, such as wind/wave energy harvested by different energy converters [1–16], where the fundamental principle includes two main points: ( 1 ) to design an energy harvesting device excited by winds or waves to undergo as large mechanical motions as possible; ( 2 ) to convert the mechanical energy into a type of storable energy, such as electric or chemical ones, etc. To obtain the large motions of an energy harvest device, the mechanisms of resonances in linear dynamic systems and the flutter of aerofoils, which are harmful vibrations to be avoided in traditional dynamic designs, are in a reverse case to be used to obtain the expended large motions of the energy harvesting device [12–16]. The review papers [4,16] provide more details on various wave/wind energy harvesters as well as many available references. The interested readers may refer to them for more information on this topic. Another research direction is to develop low-emission and fuel-efficient designs to reduce energy consumption. In this direction, the majority of current publications aim to harvest the energy of vibrations in various transportation modes: airplanes, trains, and vehicles, etc. These vibrations are caused by different operation environments, for example, the vehicle vibrations excited by uneven roads, airplane cabin noises caused by engines. The energy of this type of vibrations is the output of burning fuels in their engines, but normally is dissipated by the damping of the transport and the environment, which is an energy loss. As indicated in Ref. [17], during the everyday usage of an automobile, only 10%–16% of the fuel energy is used to drive the vehicle to overcome the resistance from the road friction and the air drag, so that about 85% of the fuel energy is dissipated by the shock absorbers of the vehicle suspension system. Therefore, there have been many research papers theoretically, numerically and experimentally dealing with different designs to harvest this type of vibration energy. In vehicle designs, applications of passive viscous shock absorbers in parallel with suspension springs aim to reduce vibrations and to achieve better riding/road handling qualities. Active suspensions have been explored [18], but their operation requires a significant amount of energy, which limits their wide implementations. To overcome this, initial theoretical researches on the feasibility to harvest the vibration energy from vehicle suspensions were reported [19,20], which examined the possibility of using a permanent magnetic motor as mechanical dampers for vehicles to dissipate their vibration energy. Following this initial idea, many publications on this topic have been reported. References such as [21–36] presented some concepts, investigations, numerical simulations, various types of absorber designs, experimental tests, piezoelectric transducers [28], as well as the involved nonlinear behaviour [31]. It has been proposed to adopt regenerative suspensions in harvesting the vehicle lost engine energy while reducing its vibrations for better ride performance. The more details on those researches one can be referred to the comprehensive review papers [29,34] as well as a series of research reports [32–36]. The company, Perpetuum [37], specialises in vibrational energy harvesting devices and has produced a range of energy harvesting nodes on rail applications and several thousand units are currently mounted on UK and European trains. In the available publications to analyse the performance of energy harvesting absorbers, the fundamental concept model is a quarter vehicle one, which has omitted the interactions between vehicle vibration modes, four absorbers as well as different uneven road excitations to four wheels. It has been noticed that a full vehicle schematic model with seven degrees of freedom (DOF) is proposed in Ref. [30], but not for energy harvesting investigations, so that there is no interaction results on vehicle vibration reduction and energy harvesting performance given from this model. This paper develops an integrated mechanical–electrical interaction system consisting of a rigid vehicle with three main important degrees of freedom to describe its heave, roll and pitch motions, four electro-magnetic regenerative suspensions, and four tires, each of them modelled by its stiffness and damper, subjected to uneven road excitations to study its complex coupling dynamic mechanisms. This model is able to reveal the effects of uneven road, the damping of the energy harvester, vehicle natural frequencies/modes, and driving speeds, etc. on the performance of the energy harvester to collect the vehicle vibration energy and at the same time to improve its ride quality. 2 Mathematical model As shown in Fig. 1a, we consider the vehicle as a rigid frame, with the sprung mass M¯ located at its mass center O at which a vehicle body coordinate system O − X Y Z with three unit coordinate vectors iX , iY , and iZ are fixed. This frame is supported by the four suspension units, of which, as shown by Fig. 1b, a representative unit I (I = 1, 2, 3, 4) at point (Xˆ I , YˆI , Zˆ I ) consists of a suspension spring of stiffness K I , an electromagnetic energy converter of damping CI , wheelmagnetic body mass m I , and wheel stiffness kI and damping cI . The vehicle moves on an uneven road surface denoted by a height function h(xI , yI ) measured in a fixed absolute coordinate system o−x yz, of which the origin o is located at point O and its three axes are respectively parallel the three axes of the body coordinate system O − X Y Z when the vehicle is in a stable state moving along an even road in the o − x direction in a constant velocity V . Uneven road excitations to the four wheels will cause the vehicle vibrations, assuming only in the vertical direction, described by a vertical translation Z of the mass center O, a roll rotation angle θx about axis o − x and a pitch rotation angle θy about axis o − y, but neglecting the motions in the horizontal translations and yaw rotations about axis o − z. The roll and pitch inertias of the vehicle about the body system are denoted by Jx and Jy , respectively. The motion z I of mass m I is also considered in the vertical direction with any horizontal motions neglected. The steady position of the system in the even road is chosen as our reference positon to measure the small vibrations of the vehicle, based on which we can derive the dynamic equations of the system as follows. 2.1 Geometrical parameters and dynamical variables We adopt the following vectors and matrices of the geometrical parameters and dynamical variables to derive the governing equations describing the interaction dynamics of the integrated system. Rotation vector of vehicle body θ = θX iX + θY iY , θ = θX θY 2 b M X O Z Y 3 d 4 b ( Xˆ I ,YˆI ) K I 2 CI ˆ RI kI mI FI K I 2 Z I zI cI hI ( x, y) x o a z 1 y a , , , Position vector of suspension unit I Rˆ I = Xˆ I iX + YˆI iY , RˆI = Xˆ I YˆI T . Displacement Z I and velocity Z˙ I of suspension unit I The vehicle body undergoes a small rigid motion consisting of the vertical translation Z and the small rotation denoted by Eq. ( 1 ), so that the displacement Z I of suspension unit I is calculated as θX θY θ˙X θ˙Y θ¨X θ¨Y Z I = Z + {θ × Rˆ I }Z = Z + θX YˆI − θY Xˆ I = Z + YˆI −Xˆ I Zs = Is Z + XsTθ , Zs = X S = Suspension force FI of unit I FI = FI i Z , FI = K I (Z I − z I ) + CI (Z˙ I − z˙ I ), 2.2 Dynamical equations of vehicle body Using Newton’s second law, we obtain the following dynamic equations governing the motions of the system. Vertical translation of vehicle body M¯ Z¨ + FI = M¯ Z¨ +IsTKs (Zs −z)+IsTCs (Z˙ s −z˙) = 0. ( 6 ) T . ( 5 ) Rotation of vehicle body Jθ¨ + (Rˆ I ×FI ) = Jθ¨ +Xs Fs = 0, J = Vertical motion of wheel I mz¨ − Fs + ks (z − h) + cs (z˙ − h˙ ) = 0, m = diag(m I ), ks = diag(kI ), cs = diag(cI ), h = h1 h2 h3 h4 MZ¨ + CZ˙ + KZ = F, ⎡ IsTCs Is IsTCs XsT −IsTCs ⎤ C = ⎣ Xs Cs Is Xs Cs XsT −Xs Cs ⎦ , −Cs Is −Cs XsT Cs + cs K = ⎣ ⎡ IsTKs Is IsTKs XsT −IsTKs ⎤ Xs Ks Is Xs Ks XsT −Xs Ks ⎦ , −Ks Is −Ks XsT Ks + ks M = diag(M¯ , J, ⎡ 0 ⎤ ⎡ 0 ⎤ F = ⎣ 0 ⎦ + ⎣ 0 ⎦ . ks h cs h˙ 2.3 Energy flow equation ⎡ Z ⎤ m), Z = ⎣ θ ⎦ , z Following the energy flow theory [38] and pre-multiplying Eq. ( 9 ) by the velocity vector Z˙ T, we can derive the energy flow equation of the system in the form T˙ + D˙ + Π˙ = P, 1 T T = 2 Z˙ MZ˙ , D = P = Z˙ TF = z˙T(ks h + cs h˙) = PI = z˙ I (kI h I + cI h˙ I ). 0 t Z˙ TCZ˙ dt, 1 ZTKZ, Π = 2 4 Here, T , D, and Π represent the kinetic, dissipated/collected, and potential energy of the system, respectively, while P, the summation of each input power PI from wheel I , gives the total power input from the road excitation. Based on the energy conservation law, this input power, originally produced by the vehicle engine, vanishes if the vehicle is driving along an even road with h I = 0. Against the uneven road, the vehicle engine will output an extra power, which is a nonuseful power to cause vibrations and an uncomfortable riding environment for passengers. In current vehicle designs, a passive suspension system is used to damp these vibrations in order to achieve riding comfort, which physically dissipates this output power by the suspension dampers, so that this power is totally wasted. In this research, we use the four electro-magnetic energy harvesting converters of damping CI to replace the passive suspension dampers in order to collect this part of engine energy and to reduce the vibrations caused by uneven road excitations at the same time. For the characteristic time period Tˆ of the road function, the time averaged energy flow equation is given by T˙ + D˙ + Π˙ = P , ( 11 ) where, for an arbitrary time-variable A(t ), such as A(t ) = P(t ), its time averaged value is defined as 1 A(t ) = Tˆ 0 T ˆ Adt . 2.4 Energy harvesting converter and its collected power As shown in Fig. 1b, the energy harvesting converter I consists of a magnetic body of intensity BI and an electric coil, of effective conduct length lI , electrical inductance L I , internal resistance rI , and capacitance Cˆ I , moving across the electromagnetic field, so that an induced voltage eI at the two ends of the power collection device of resistance Rˆ I is produced. Therefore, an electric charge Qˆ I with the corresponding current IˆI is generated in the moving coil, from which an electromagnetic force FˆI between the moving coil and the magnetic body is further produced to resist their relative motion. Based on the Laplace electromagnetic theorem and electrical equilibrium principle [12,15,16,39,40], the following equations are established. eI = BI lI (Z˙ I − z˙ I ), FˆI = IˆI BI lI . ( 13 ) The electric equilibrium equation of the circuit of total resistance RI = Rˆ I + rI is given by eI = L I Q¨ˆ I + RI Qˆ˙ + Qˆ I /Cˆ I = L I Iˆ˙I + RI IˆI + Qˆ I /Cˆ I , IˆI = Qˆ˙ I . ( 14 ) From Eqs. ( 13 ) and ( 14 ), we obtain the following relationship between the current and the mechanical motions, L I I¨ˆI + RI Iˆ˙I + IˆI /Cˆ I = BI lI (Z¨ I − z¨ I ). Equation ( 15 ) is coupled with Eq. ( 9 ), so that the integrated system is an electric-mechanical interaction system. Considering practical cases for which the capacitance Cˆ I of the moving coil is normally quite large, and the inductance L I is small, especially for the low frequency road excitations, we can reduce Eq. ( 15 ) to RI IˆI = BI lI (Z˙ I − z˙ I ), from which, using Eq. ( 13 ), we have the following electromagnetic force ( 12 ) ( 15 ) ( 16 ) ( 9 ) ( 10 ) FˆI = (BI lI )2 Rˆ I + rI (Z˜˙ I − z˜˙ I ) = CI (Z˜˙ I − z˜˙ I ), CI = (BI lI )2 Rˆ I + rI ( 17 ) As used in Eq. ( 9 ), CI is called the damping coefficient of the energy converter, which is determined by the parameters of the electromagnetic suspension unit. The collected power of the energy converter I is calculated by These vectors give the vibration levels of displacement, velocity, and acceleration at the vehicle mass center and the four points on the vehicle body at which the four energy converters are connected. PˆI = IˆI2 Rˆ I = CI Rˆ I (Z˙ I − z˙ I )2 RI = , Pˆ = PˆI , 4 of which the percentages of energy collection with the corresponding time averaged values are respectively given by λI = PˆI / PI , λ¯I = PˆI / PI , λ = Pˆ / P, λ¯ = Pˆ / P . 2.5 Road uneven functions The uneven road is described by the following function with unit amplitude h(x , y), |h(x , y)| ≤ 1, and this function may be considered as a non-dimensional shape function, of which the positive value implies a convex road while the negative value denotes a concave road. Since the system is considered as linear system, if the amplitude of the uneven road is not a unit, such as H , the dynamic response and the power can be obtained by multiplying them by H and H 2, respectively. Considering our assumption of the vehicle moving along o − x direction in the constant speed V and the distance between the front and rear wheels, we can represent the road functions for each wheel as follows h I = h I (V t, yI ), front: I = 1, 2, h I (V t − a − b, yI ), rear: I = 3, 4, in which the time when the front wheels reach the uneven road surface, is taken as the initial time t = 0 in our analysis. 2.6 Vibration isolation factor vectors To study the performance of the proposed energy harvesting device for reducing the vehicle vibrations due to uneven road excitations, we define the following vibration isolation factors. Based on Eq. ( 9 ), we consider the motions of the vehicle body caused by a unit maximum road height, so that we define the amplifying factor vectors αd = Z ZsT T , αv = T T Z˙ Z˙ s , αa = T T Z¨ Z¨ s . ( 18 ) ( 19 ) ( 20 ) ( 21 ) where q, a time function, is called as the generalized coordinate vector of the mode space. Substituting Eq. ( 29 ) into Eq. ( 9 ), and then pre-multiplying the resultant equation by 3 Solution approach 3.1 Natural vibration To use the mode summation method to solve Eq. ( 9 ), we need to solve its corresponding natural vibrations governed by Eq. ( 9 ) with no damping (C = 0) and no external force (F = 0), i.e., MZ¨ + KZ = 0. Z = ϕejΩt , We consider a solution of this equation in the form which, when substituted into Eq. ( 23 ), gives (−Ω2M + K)ϕ = 0. The condition for a non-zero solution of Eq. ( 25 ) requires its determinant of the coefficient matrix to vanish, which gives the characteristic equation (−Ω2M + K) = 0. The solutions of Eq. ( 26 ) are the seven natural frequencies Ωi (i = 1, 2, . . . , 7), of the system, which is represented by a diagonal matrix Λ = diag(Ωi ). Substituting the natural frequency Ωi into Eq. ( 25 ), we obtain the corresponding natural mode vector ϕi satisfying the orthogonal relationships ΦTMΦ = I, ΦTKΦ = Λ2, Φ = ϕ1 ϕ2 · · · ϕ7 . 3.2 Mode summation Using the mode summation method, we represent the vehicle displacement vector Z = Φq, ( 23 ) ( 24 ) ( 25 ) ( 26 ) ( 27 ) ( 28 ) ( 29 ) ΦT as well as using the orthogonal condition Eq. ( 28 ), we obtain the mode equation Iq¨ + dq˙ + Λ2q = f , d = ΦTCΦ, f = ΦTF. ( 30 ) Here, f is called as the generalized force vector and d is the damping matrix in the mode coordinate space. This damping matrix is not a diagonal matrix in general cases. To judge the efficiency of energy collection and vibration isolation provided by the energy converters, we assume that the uneven road provides a harmonic excitation of unit amplitude and frequency ω˜ , represented by a complex force f = Fejω˜ t˜, of which its imaginary part denotes the real force, i.e., f = Im{f } = Im{Fejω˜ t˜}. In this representation form, the dynamic response of the system can also be represented in a complex form q = Im{q}, q = Qejω˜ t˜, which transforms Eq. ( 30 ) into a corresponding complex form (Λ2 − ω˜ 2I + jω˜ d)Q = F. Pre-multiplying this equation by Λ−2, we obtain a nondimensional equation (I − η2 + 2jηζ )Q = F, F = Λ−2F, η = diag(ω/Ωi ), ζ = Λ−1d/2. Here, η is a non-dimensional diagonal matrix of frequency ratio and ζ is a non-dimensional mode damping matrix. From this equation, we obtain the complex amplitude of the generalized coordinate q in the form Q = (I − η2 + 2jηζ )−1F, ( 31 ) ( 32 ) ( 33 ) ( 34 ) ( 36 ) h I = sin(β V t ) = sin(V t ), front wheels: I = 1, 2, sin[β(V t − a − b, )] = sin(V t − a − b, ), rear wheels: I = 3, 4, and their corresponding imaginary parts are the physical variables. As indicated in the book [38], when calculating powers by Eqs. ( 10 ) and ( 18 ), you have to use the real physical variables to get their correct results. 4 Case study 4.1 Vehicle parameters We use the vehicle parameters given in Ref. [41]: sprung mass M¯ = 1440 kg with its roll inertia JX = 900 kg · m2 and pitch inertia JY = 2000 kg · m2; distance of mass center from front axle a = 1.016 m and from rear axle b = 1.524 m; front/rear track width d = 1.5 m; front suspension stiffness K I = 35, 000 N · m−1 with energy harvester damping coefficient CI = 2500 N · s · m−1 (I = 1, 2); rear suspension stiffness K I = 30, 000 N · m−1 with damping coefficient CI = 2000 N · s · m−1 (I = 3, 4); front/rear un-sprung mass m I = 80 kg with tire stiffness kI = 200,000 N · m−1 and damping cI = 10 N · s · m−1 (I = 1, 2, 3, 4). These parameters and their geometrical layout shown in Fig. 1, when substituted into Eq. ( 3 ), give the related prescribed matrices in Eq. ( 9 ), such as X S = −d/2 −a d/2 −a −d/2 b d/2 b We use the maximum height of the uneven road to measure the length, so that the road function will provide a non-dimensional unit height road excitation. You can assume any types of road functions, continues or discrete, to use the developed model for the corresponding simulation. Here, for an illustration purpose, we choose the following two road functions. Sinusoidal road function ( 37 ) ( 38 ) ( 39 ) which, when substituted into Eq. ( 29 ), gives the complex/real physical displacement vector h(x , y) = sin β x , Z = ΦQejω˜ t˜, Z˜ = Im{Z}. It should be mentioned herein, we use the complex representations to solve Eq. ( 9 ), of which the solutions are complex ( 35 ) A sinusoidal road function with unit amplitude is defined as where β is a non-dimensional parameter involving the spatial period of the road function. The big value of parameter β corresponds a short spatial period of road function. For simplicity without losing the generality, we choose β = 1 in the following derivations, so that Eq. ( 21 ) becomes of which the time derivative Therefore, from Eq. ( 9 ) we find that the vehicle moving speed V involves a road excitation frequency, i.e., ω = β V = V . In this equation, the vehicle speed is measured by the road amplitude, so that its physical dimension, m·s−1·m−1 = s−1, is just the dimension of angular frequency ω.The corresponding road vector of this road function is denoted by a complex vector as follows h = Im{h}, h = 1 1 e−j(a+b) e−j(a+b) T ejV t . (41b) (41a) A sine convex or concave road function To consider a type of ramp on roads, we assume a sine convex or concave road function defined in the range 0 ≤ x ≤ x0 in the form h˙ I = h = h˙ = h¨ = h I = ± sin x , 0, ±V cos x , 0, ∓V 2 sin x , 0, 0 ≤ x < x0, x0 ≤ x , 0 ≤ x < x0, x0 ≤ x , 0 ≤ x < x0, x0 ≤ x . β V cos(V t ) = V cos(V t ), β V cos[β(V t − a − b, )] = V cos(V t − a − b, ), front wheels: I = 1, 2, back wheels: I = 3, 4. ( 40 ) shown in Fig. 2. Table 1 gives the detailed information of frequencies and mode shape characteristics of these seven natural modes. For each mode shape in Fig. 2, the left side is the front of vehicle, the dashed line denotes the stationary position of the vehicle while the real line represents the mode shape displacement, i.e., shift from the dashed position. The first three modes are observed with very small motions of the energy harvesters, while the last four modes are with the large energy harvester motions. The mode 1 of frequency 1.3274 Hz shows a rolling mode of the vehicle while the modes 2 and 3 represent two pitch modes of the vehicle in which mode 2 of frequency 1.3424 Hz and mode 3 of frequency 1.5682 Hz relate to the large motion amplitudes at the vehicle rear and front sides, respectively. Mode 4 of frequency 8.5404 Hz and mode 5 of frequency 8.556 Hz show respectively that the harvester masses 3 and 4, at the vehicle rear side, undergo the rolling and heave motions, while mode 6 of 8.6366 Hz and mode 7 of 8.6450 Hz behave the opposite motion shapes of the harvester masses 1 and 2, connected to the vehicle front side, with the rolling and heave motions. In the energy harvesting suspension unit designs for vehicles to isolate the vibrations caused by road excitations as well as to collect this vibration energy, the following general rules involving natural modes may be followed. ( 1 ) The energy collected from the energy harvester is proportional to the square of the relative velocity between the vehicle and harvester mass motions, as given by Eq. ( 18 ), so that the modes with large relative velocities will provide the large contribution for energy collection and vibration reduction. ( 2 ) In principle, considering the symmetry of the system, for symmetrical road excitations, the dynamic responses of the anti-symmetrical modes vanish, therefore, they will not provide any contributions for energy collection and vibration isolation caused by the symmetrical road excitations. Similarly, for anti-symmetrical road excitations, symmetrical modes have also no contributions for harvesting energy with vibration isolation. As shown in Fig. 2 and Table 1, the modes 2, 3, 5, 7 provide symmetric motions about the central line of the vehicle, while the modes 1, 4, 6 are antisymmetrical. Therefore, for a single ramp in Eq. (42a), the road excitations applied to the four wheels take the form for the front wheels sin(V t ), 0, 0 ≤ t < x0/ V , x0/ V ≤ t, I = 1, 2, (42b) and for the rear wheels ⎧ 0, 0 ≤ t < (a + b)/ V , h I = ⎨ sin(V t ), (a + b)/ V ≤ t < x0/ V , ⎩ 0, x0/ V ≤ t. 4.3 Simulation results 4.3.1 Natural frequencies and modes (42a) I = 3, 4, (42c) The natural vibrations of the investigated system are governed by Eqs. ( 22 )–( 28 ), from which, when solving Eq. ( 25 ) for eigenvalues/eigenvectors, the following seven natural frequencies with the corresponding modes are obtained as -2 -1 f5: 8.5560 Hz 0 -1 f2: 1.3414 Hz 4.3.2 Energy collection and vibration isolation Based on the mode summation approach and choosing different road functions with the unit amplitude as given by Eq. ( 41 ) or ( 42 ), we solve the dynamic equations of the system to investigate the performance of energy collections and vibration isolations of the proposed electromagnetic energy harvester for cases 1, 2, 3, 4, respectively with the energy harvester damping coefficients 25%, 50%, 75%, and 100% CI . The simulation results are given in the frequency response spectrums or the time histories as follows. Symmetrical continuous sinusoidal road excitation For a continuous sinusoidal road excitation defined by Eq. ( 41 ), where the vehicle speed links the angular frequency of the road excitation. After solving the dynamic equations subject to this road excitation, we obtain the following frequency response spectrums of displacements, velocities, accelerations at chosen points of the system as well as the road input energy and the harvested energy for the four cases with different energy harvester damping coefficients. In these figures, the horizontal axes denote the road excitation frequency ω = V , i.e., the vehicle speed, while the vertical axes represent the related dynamic response variables excited by a unit amplitude (m) road excitation. These curves provide the base to examine the performance of the energy harvester suspension system in reducing vibrations and harvesting energy. Figure 3 shows the response curves of vertical displacement Z , velocity, and acceleration at the mass center and the roll/pitch angles θx and θy with the corresponding angular velocities and angular accelerations of the vehicle, which describes the integrated motion of rigid vehicle. The energy Mode number Vehicle (V) /E-converter (E) motion V. Roll V. Rear pitch V. Front pitch E. Rear roll E. Rear pitch E. Front roll E. Front pitch collection of each harvester involves the relative motion between the vehicle and the harvester body, so that as an example, Fig. 4 presents the dynamic response spectrums of the displacement z1, velocity z˙1, and acceleration z¨1 at the mass m1 of energy converter 1, the corresponding ones Z1, Z˙ 1, Z¨ 1 at vehicle point 1, where converter 1 fixed, as well as their differences z1 − Z1, z˙1 − Z˙ 1, z¨1 − Z¨ 1. Figure 6 gives the input power by the road excitation, the collection power by the energy harvester as well as the energy collection efficiency, in which for the lower frequency range the curves have been locally amplified to get its more details. Observing these three Figs. 3–5, we can clearly see that the investigated symmetrical road excitation excites only symmetrical motions consisting of four symmetrical modes shown in Table 1 with the curve peaks located around 8−9 rad · s−1 and 53–54 rad · s−1. Especially, as shown in Fig. 4, the dynamic response of anti-symmetric rolling motion θx of the vehicle about x -axis is extremely small in the calculation error range. With increasing the energy harvester damping values, the vibration peaks of vehicle motions are much reduced, so that the power input due to vibrations and the harvesting energy are also reduced. Since in the simulations, we assume that the tire damping is quite small and neglect other types of physical damping, therefore, the nearly full vibration energy is collected. For practical cases, due to existence of various physical damping, a part of vibration energy will be dissipated, so that the energy collection efficiency is less than 1. Non-symmetrical continuous sinusoidal road excitation on the right side of vehicle As discussed above, the symmetrical road excitation cannot excite the responses of the anti-symmetrical modes of the vehicle. Here, we consider a non-symmetrical continuous sinusoidal road excitation defined also by Eq. ( 41 ), but only acting on the right side wheels of the vehicle, i.e. the left side wheels of vehicle running along an even road with no road excitations. Figures 6–9 give the dynamic response spectrums in the similar structures of Figs. 3–5, respectively. Since the road excitation is non-symmetrical, so that all modes of the system contribute to the dynamic responses. As shown in Fig. 6 for this excitation, the roll rotation motion θx , θ˙x , θ¨x , of the vehicle does not vanish. Furthermore, a comparison of the curves in Figs. 7 and 8 shows that the dynamic responses of the variables at positions 1 and 2 are obviously different for this case of the non-symmetrical road excitation applied only on the right side of vehicle. The response values on the left side position 2 shown in Fig. 8 are much smaller than the ones on the right side position 1 as given in Fig. 7. For this non-symmetrical continuous sinusoidal road excitation, we consider the road excitation of amplitude 0.2 m and use the standard Newmark time integration method [42] with time step t = 0.0005 s in the simulations to reveal the characteristics of time histories of the integrated system, which are shown in Figs. 10 and 11. Figure 10 shows the time histories of the displacements of vehicle points and suspension masses at positions 1 and 2 as well as the ones of the displacement at vehicle mass center and its roll/pitch rotation angles, while Fig. 11 gives the time histories of input and collection powers. It is observed that: ( 1 ) The displacement at excitation side 1 is larger than the non-excited side 2. ( 2 ) At the initial road impact instance, although the vehicle body suffers big motions but its order is much lower than the suspension mass motions, which absorbs the impact energy and reduces the road impact transmitting to the vehicle body. ( 3 ) Since the system is linear, the stable response shows a sinusoidal form as a function of the frequency (vehicle speed). ( 4 ) The input and output collection powers are positive, of which the values at vehicle speed 50 are maximum. This is because the speed (ω = 50) is near to the resonance peak shown in Fig. 7, so that there is large relative motion z1 − Z1 from which the large road excitation energy is collected. A sine convex ramp Here, we consider the investigated vehicle, moving in different speeds V = 30, 50, 70 m · s−1, is impacted by a half sine convex ramp (amplitude 0.2 m, x0 = π in Eqs. (42b, 42c)) only at positions 1 and 3 on its right-side wheels in succession. Based on the standard Newmark approach [42], Fig. 12 gives the time histories of the displacements of four vehicle points (Z1 − Z4) and the 4 suspension masses points (z1 − z4) a b 15 while Fig. 13 shows the time histories of the displacement Z at mass center and roll θx /pitch θy rotation angles of the vehicle. Figures 14 and 15 provide the instant time histories of road input powers and output collection powers at four suspension positions and the total ones of the integrated vehicle, respectively. These time histories show that: C C C C a a b b c c C C C C a b 0.6 0.4 m0.2 0 -0.2 -0.40 0 0.5 2 2.5 0.5 2 Fig. 9 (Color online) Frequency spectrums of powers. a Road excitation input. b Power collection. c Energy collection efficiency = collection power/input power Potential further research directions The following directions may be considered for future research. ( 1 ) This paper has conducted research only on a fundamental concept but has not presented an idea on how practically to collect the harvested energy. The best possible way is to design an electrical device to transfer the collected energy into the battery of vehicle. However, since the electrical device operation needs power, a proper evaluation of final economic efficiency is required. Another option is to apply the collected energy to some vehicle facilities, such as window deforester, etc. To realize any 0.5 0.5 0.1 0.05 -0.1 0 a 15 1 Time (s) 1 Time (s) 1 Time (s) 2 1 Time (s) 1 Time (s) 1 Time (s) Fig. 12 (Color online) Time histories. a Displacements (Z1 – Z4) of four vehicle points. b Displacements (z1 − z4) at four suspension masses points possible practical designs, analysis and tests are necessary. ( 2 ) The other three motions of rigid vehicle, sway in y direction, surge in x -direction and yaw about z-axis may be considered for a full rigid vehicle model. Especially, the sway and yaw motions are important for vehicle control performance. ( 3 ) The further detailed research may consider the vehicle as an elastic body with more practical structure, for which finite element method may be used, but the solution approaches discussed in this paper are still valid. The main difference of this case is the order of the final dynamic matrix equation is higher than 7. ( 4 ) The harvester damping depends on the parameters of electromagnetic body as well as its moving electric coil, therefore practical designs must be analyzed and tested to obtain the required damping value of energy harvester. 5 Conclusions This paper has developed an integrated numerical model consisting of a rigid vehicle, four electromagnetic energy V = 30 m s-1 V = 50 m s-1 V = 70 m s-1 2 2 Z 4 0.1 -0.01 -0.02 0.4 Time (s) 0.6 0 0.2 a b 15 10 1 - s 5 m N 0 -5 7 6 5 1 0 -1 s4 m 3 N 2 2 2 Z 1 Time (s) 1 Time (s) No. 1 No. 2 No. 3 harvesters, and four vehicle tires to investigate and reveal the characteristics of the system in harvesting energy while reducing vibrations caused by uneven road excitations. The mathematical formulations with their solution approaches to the mechanical–electromagnetic interaction system subject to various uneven road excitations are presented. The theoretical analysis and numerical simulations demonstrate that this type of energy harvesting provides a damping mechanism in the integrated coupling system. The numerical examples for the continuous sine wave road excitation and the single ramp impact are simulated, from which the following guidelines are obtained: ( 1 ) the current passive vehicle suspension units can be replaced by the investigated electromagnetic energy harvesters in reducing vibrations 2 2 2 0.5 1 Time (s) 0.5 1 Time (s) V=30mm/ss-1 V=50m/ss-1 VV=70mm/ss-1 3 3 3 No. 4 V = 30mm/s-1 V = 50mm/s-1 V = 70mm/s-1 1 Time (s) 1 Time (s) 2 2 b 10 -1 s m5 N 0 0 0.5 0.5 caused by uneven roads and at same time harvesting the lost engine energy; ( 2 ) the model with seven DOF reveals the more complex interaction characteristics of the system than a quarter vehicle model, so that it can address the efficiency of vibration isolation/energy harvesting affected by the natural modes/frequencies of the vehicle, the possible different road excitations on four wheels, vehicle speeds and road impact time delays for the front and rear wheels, etc. This model provides a means for vehicle electromagnetic suspension unit designs in engineering. Potential further research directions are suggested for readers to consider in future. Acknowledgements The authors acknowledge CSC for supporting S. Zhou to visit University of Southampton for one year to engage in this research, as well as Harbin Engineering University for supporting J. T. Xing to visit Harbin Engineering University (Grant HEUCF160104). 1. Thorpe , T.W.: A brief review of wave energy. A report produced for The UK Department of Trade and Industry . ETSU R-120, AEA Technology Plc. ( 1999 ) 2. Office of Naval Research: Environmental assessment for proposed wave energy technology project in Kaneohe Bay, Report by Office of Naval Research and US Department of the Navy, Hawaii . ( 2003 ) 3. Bedard , R. , Hagerman , G. , Previsic , M. , et al: Offshore Wave Power Feasibility Demonstration Project: Final Summary ReportProject Definition Study , EPRI Global WP 009 - US Rev 1 , The Electric Power Research Institute ( 2005 ) 4. Paasch , R. , Ruehl , K. , Howland , J. , et al.: Wave energy: a pacific perspective . Philos. Trans. R. Soc. A 370 , 481 - 501 ( 2012 ). doi: 10 . 1098/rsta. 2011 .0225 5. U.S. Department of the Interior: Technology white paper on wave energy potential on the U.S. Outer Continental Shelf . http:// ocsenergy.anl.gov ( 2006 ). (Accessed 3 Oct 2016 ) 6. Washio , Y. : Wave energy research and development at JAMSTEC, offshore floating wave energy device, Mighty Whale . http://www. takesteps.org/empower/exhibition ( 2006 ). (Accessed 3 Oct 2016 ) 7. Ingram , D. : Edinburgh makes waves in power generation research . http://www.ed.ac.uk/files/atoms/files/engineering-wavemakers. pdf ( 2016 ). (Accessed 3 Oct 2016 ) 8. Pelamis Wave Power: World's first commercial wave energy project , Agucadoura, Portugal. http://www.power-technology. com/projects/pelamis/ ( 2016 ). (Accessed 3 Oct 2016 ) 9. Wave Dragon : Wave Dragon has started the development of a 1.5 MW North Sea Demonstrator , March 2011 . http://www. wavedragon.net ( 2016 ). (Accessed 3 Oct 2016 ) 10. Wave Plane Productions: Wave Plane-the only device in the world . http://www.waveplane.com ( 2016 ). (Accessed 3 Oct 2016 ) 11. Falnes , J.: Ocean Waves and Oscillating Systems , Linear Interactions Including Wave-Energy Extraction . Cambridge University Press, London ( 2002 ) 12. Xing , J.T. , Xiong , Y.P. , Tan , M. , et al.: A numerical investigation of a wave energy harness device-water interaction system subject to the wave maker excitation in a towing tank . In: Proceedings of the ASME 28th International Conference on Ocean, Offshore and Arctic Engineering , Honolulu, Hawaii, 1 - 10 ( 2009 ) 13. Xu , X. , Pavlovskaia , E. , Wiercigroch , M. , et al.: Dynamic interactions between parametric pendulum and electro-dynamical shaker . Z. Angew. Math. Mech . 87 , 172 - 186 ( 2007 ) 14. Yang , J. , Xiong , Y.P. , Xing , J.T. : Investigations on a nonlinear energy harvesting system consisting of a flapping foil and an electro-magnetic generator using power flow analysis . In: Proceedings of the ASME 2011 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference (IDETC/CIE 2011 ), Washington, DC, 1 - 8 ( 2011 ) 15. Xing , J.T. , Xiong , Y.P. , Wiercigroch , M. , et al.: A mathematical modelling for an integrated electric converter-nonlinear oscillatorwater interaction system to harness wave energy . In: Proceedings of 7th European Nonlinear Dynamics Conference (ENOC 2011 ), Rome, Italy, 1 - 6 ( 2011 ) 16. Xing , J.T.: Developments of numerical methods for linear and nonlinear fluid-solid interaction dynamics with applications . Adv. Mech . 46 , 201602 ( 2016 ). doi: 10 .6052/ 1000 -0992-15-038 17. Zuo , L. , Scully , B. , Shestani , J. , et al.: Design and characterization of an electromagnetic energy harvester for vehicle suspensions . Smart Mater. Struct . 19 , 045003 ( 2010 ). doi: 10 .1088/ 0964 -1726/ 19/4/045003 18. Sharp , R. , Crolla , D. : Road vehicle suspension system design-a review . Vehicle Syst. Dyn . 16 , 167 - 192 ( 1987 ) 19. Karnopp , D. : Permanent magnet linear motors used as variable mechanical dampers for vehicle suspensions . Veh. Syst. Dyn . 18 , 187 - 200 ( 1989 ) 20. Karnopp , D. : Power requirement for vehicle suspension systems . Veh. Syst. Dyn . 21 , 65 - 71 ( 1992 ) 21. Segel , L. , Lu , X.P. : Vehicular resistance to motion as influenced by road roughness and highway alignment . Aust. Road Res . 12 , 211 - 222 ( 1982 ) 22. Hsu , P. : Power recovery property of electrical active suspension systems . In: Proceedings of the 31st Intersociety Energy Conversion Engineering Conference (IECEC 96) , Washington, DC, 1899 - 1904 ( 1996 ) 23. Goldner , R. , Zerigian , P. , Hull , J.: A preliminary study of energy recovery in vehicles by using regenerative magnetic shock absorbers . SAE Technical Paper , 2001 - 01 - 2071 ( 2001 ). doi: 10 . 4271/2001-01-2071 24. Abouelnour , A. , Hammad , N.: Electric utilization of vehicle damper dissipated energy . In: Al-Azhar Engineering Seventh International Conference (AEIC), Cairo Egypt, April 7 - 10 ( 2003 ) 25. Zheng , X. , Yu , F. : Study on the potential benefits of an energyregenerative active suspension for vehicles . SAE Technical Paper , 2005 - 01 - 3564 ( 2005 ). doi: 10 .4271/2005-01-3564 26. Kawamoto , Y. , Suda , Y. , Inoue , H. , et al.: Modelling of electromagnetic damper for automobile suspension . J. Syst. Des. Dyn . 1 , 524 - 535 ( 2007 ) 27. Zhang , Y. , Huang , K. , Yu , F. , et al.: Experimental verification of energy-regenerative feasibility for an automotive electrical suspension system . In: Proceedings of the IEEE International Conference on Vehicular Electronics and Safety , Beijing, December 13- 15 ( 2007 ) 28. Priya , S.: Advances in energy harvesting using low profile piezoelectric transducers . J. Electroceram . 19 , 165 - 182 ( 2007 ) 29. Zhang , Y. , Yu , F. , Huang , K. : A state of art review on regenerative vehicle active suspension . In: Proceedings of the 3rd International Conference on Mechanical Engineering and Mechanics (ICMEM) , Beijing, October 21- 23 ( 2009 ) 30. Montazeri-Gh , M. , Soleymani , M. : Investigation of the energy regeneration of active suspension system in hybrid electric vehicles . IEEE Trans. Ind. Electron . 57 , 918 - 925 ( 2010 ) 31. Ando , B. , Baglio , S. , Trigona , C. , et al.: Nonlinear mechanism in MEMS devices for energy harvesting applications . J. Micromech. Microeng . 20 , 125020 ( 2010 ). doi: 10 .1088/ 0960 -1317/20/12/ 125020 32. Li , Z. , Brindak , Z. , Zuo , L. : Modeling of an electromagnetic vibration energy harvester with motion magnification . In: ASME 2011 International Mechanical Engineering Congress and Exposition , Denver, Colorado, 285 - 293 ( 2011 ) 33. Li , Z. , Zuo , L. , Kuang , J. , et al.: Energy-harvesting shock absorber with a mechanical motion rectifier . Smart Mater. Struct . 22 , 025008 ( 2013 ) 34. Zuo , L. , Tang , X. : Large-scale vibration energy harvesting . J. Intell. Mater. Syst. Struct . 24 , 1405 - 1430 ( 2013 ) 35. Li , Z. , Zuo , L. , Luhrs , G. , et al.: Electromagnetic energy-harvesting shock absorbers: design, modeling, and road tests . IEEE Trans. Vehicul. Technol . 62 , 1065 ( 2013 ) 36. Zuo , L. , Zhang , P.: Energy harvesting, ride comfort and road handling of regenerative vehicle suspensions . J. Vib. Acoust . 133 , 1 - 8 ( 2013 ) 37. Perpetuum : PMG17 vibration energy harvester-IEEE Xplore Digital Library . doi: 10 .1109/ MCS . 2007 .911267 38. Xing , J.T.: Energy Flow Theory of Nonlinear Dynamical Systems with Applications . Springer, Heidelberg ( 2015 ) 39. Kittel , C. : Berkeley Physics Course. McGraw-Hill , New York ( 1967 ) 40. Xing , J.T. , Xiong , Y.P. , Price , W.G.: A generalized mathematical model and analysis for integrated multi-channel vibration structure-control interaction systems . J. Sound Vib . 322 , 584 - 616 ( 2009 ) 41. Shim , T. , Ghike , C. : Understanding the limitations of different vehicle models for roll dynamics studies . Veh. Syst. Dyn . 45 , 191 - 216 ( 2007 ) 42. Bathe , K.J.: Finite Element Procedures. Prentice Hall , New Jersey ( 1996 )

This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1007%2Fs10409-017-0647-x.pdf

Jing Tang Xing, Zhe Sun, Sulian Zhou, Mingyi Tan. Dynamic interactions of an integrated vehicle–electromagnetic energy harvester–tire system subject to uneven road excitations, Acta Mechanica Sinica, 2017, 440-456, DOI: 10.1007/s10409-017-0647-x