Research on Ride Comfort Model of Wheel Motor Driving Vehicle Based on Matlab / Simulink

Simulink is a visual simulation tool in MATLAB; through Simulink software, to establish a model can reduce the amount of programming workload, and improve the efficiency of the establishment of automotive models. The ride comfort of the vehicle is a measure of the most basic indicators of a car performance. By establishing a ride comfort model in Matlab/Simulink, the wheel motor electric vehicle mainly affects the smoothness of the car mainly in the following aspects: pavement, tire, suspension, motor and so on. Through the establishment of the above model, we can effectively study the wheel motor drive electric vehicle ride comfort research.


Introduction
With the development of society and scientific progress, people's requirements on the ride comfort (that is, the car ride comfort) are getting higher and higher; the car's ride comfort is mainly to keep the car in the process of vibration and shock environment.The impact on the comfort of the occupants is within a certain limit [1].Due to the limited number of experimental equipment, a lot of experiments need to simulate the simulation to establish a virtual simulation environment, as well as virtual conditions.However, the most important point of simulation is the establishment of simulation model.The closer the simulation model is to the real conditions, the more conducive the experiment will be.In 2009, Zhang Lu of Wuhan University of Technology put forward different me-

The Establishment of Random Pavement
Now that the road surface roughness of the time domain model methods are harmonic superposition method, white noise method and other methods [6] [7] These methods have made some assumptions on the road, assuming that the road contour in line with the Gaussian distribution, due to white The noise method is realized by the differential equation, and the calculation is relatively small, using MATLAB/Simunlink to facilitate modeling.As shown in Figure 1.In order to reflect the authenticity of the road, a white noise model with cutoff frequency is adopted.The model can be used to reflect the actual characteristics of the road spectrum closer to the horizontal line at low frequencies.
The following formula can be used to analyze the pavement in time domain for different pavement grades and vehicle speeds.
This formula, ( ) x t is a random displacement of the front wheel, unit m; oo n is the cut-off space frequency,

Reference Coordinate System Selection
Before the simulation model is built, the reference coordinate system should be developed so that the movement state of the vehicle can be better described.
Commonly used reference coordinate system mainly has tire coordinate system, vehicle coordinate system and earth coordinate system.The vehicle coordinate system is a moving coordinate system, used to describe the movement of the car state, which is consolidated in the movement of the car [9].
Figure 2 shows the consolidation of the car in the oyez right-angled coordinate system is the vehicle coordinate system.When the vehicle is in a stationary state on the horizontal road, the x-axis is parallel to the ground, and the y-axis points to the left side of the driver.The z-axis passes through the center of the center, and the origin of the coordinate system.Can make it coincide with the center of mass [10].
Figure 3 for the tire coordinate system, perpendicular to the wheel rotation axis of the tire in the sub-plane known as the wheel plane, Z axis perpendicular to the ground plane, the provisions of the above point is positive.Y axis in the ground plane, the provisions of the direction of the direction of the wheel when the direction of the left is positive, the wheel plane and the ground plane of the intersection of the X axis, the provisions of forward is [11].The origin of the coordinate system is the intersection of the intersection of the wheel plane and the ground plane and the projection line of the wheel rotation axis on the ground plane.
The geodetic coordinate system is used to describe the absolute motion of the vehicle.Its origin is fixed in a certain coordinate point on the ground.The positive direction of the Z axis is the direction perpendicular to the ground, the positive direction of the X axis is the direction of the initial state when the vehicle is running.The positive direction of the Y axis is the left side of the vehicle's initial state direction.

Body Dynamics Model
According to the vehicle in Figure 4 in the XY plane within the force analysis and coordinate system requirements, in the horizontal plane, the car's motion equation is: ( ) ( ) sin ( ) ( ) ( ) ( ) Let This formula, i is the wheel number ( 1, 2,3, 4 i = ); s M is the vehicle's sprung mass; i M is the total mass of the car; Ti δ is the wheel steering angle, rad; f F is the tire rolling resistance; a ρ is the air density, kg•m 3 ; xx I , yy I , zz I Respectively, for the vehicle around the x, y, z axis of inertia; g is the gra- vitational acceleration; ui w is the vertical speed of the first wheel, 1, 2,3, 4 i = ; i B is the vehicle suspension related parameters; g a is the front projection area of the car; g a is the vertical velocity of the body; g a is the slope of the road; yi F is the tire lateral force, kg•m 2 ; xxs I , zzs I , yys I are the moment of inertia of the sprung mass around their respective axes of rotation, kg•m 2 ; For the corresponding one of the wheel suspension geometry parameters; si K The action of the first wheel on the suspension; si K is the stiffness of each suspension; si C is the damping of each suspension; D C is the air resistance coefficient; θ is the pitch angle; ϕ is the roll angle; i d Distance from the center of the sus- pension to the center of mass of the sprung mass; ui z The center of mass of the non-sprung mass; v is the lateral velocity of the body; s z is the centroid height of the sprung mass; p is the roll angle; q is the pitch angular velocity; r is the yaw rate; u is the longitudinal velocity of the body; xi F is the tire longitudinal force [12].
According to the equation to build out the body model shown in Figure 5.

Tire Model
In this paper, the "magic formula" tire model to carry out dynamic simulation analysis.Because the tire model can accurately predict the tire performance, it is widely used, it is a combination of triangular function to fit the tire test data, obtained a set of the same form and can simultaneously express longitudinal force, lateral force and back Positive torque of the tire model [13], in the form of the following: where, y is the positive moment; Lateral force or lateral force; C represents ( ) G x is the weight function of the slip and side-side com- pound conditions; B represents the stiffness factor; D represents the peak coefficient; x can be expressed as the tire side angle or the longitudinal slip rate, respectively; y S represents the vertical offset; E represents the curvature factor; h S represents the horizontal offset.The magic formula is illustrated by Figure 6, and the curve shown in the figure may be a positive moment, a lateral force or a longitudinal force curve, where, when determined, the product of the coefficient, which is the slope of the curve, that is, the shape of the curve, which determines the shape of the curve.The value can be determined by the peak and steady-state values [14].The coefficient is used to control the curvature at the peak of the curve,: the amount

Motor Model
In this paper, the four-wheel hub motor is a permanent magnet brushless DC motor.This kind of motor has the advantages of long service life, strong overload capacity, high power density, strong reliability and large starting torque compared with other kinds of motors, with a wide range of speed and good mechanical properties of the advantages [15].
The maximum braking torque of the hub motor has a great influence on the regenerative braking distribution strategy of the wheel motor-driven electric vehicle, which in turn affects the effect of braking energy recovery [16].So the choice of wheel motor to meet the requirements of the electric car, the following is the selected motor, the motor parameters shown in Table 1.
Figure 7 for the Matlab/Simulink platform to build the motor model.
Four wheel hub motor drive electric vehicle motor speed and speed between the following relationship:  0.377 where, r is the wheel rolling radius; n is the motor speed, (r/min); v is the driving speed, km/h.The relationship between the wheel brake force and the wheel motor braking torque is: where, m T is the wheel motor braking torque; w r is the wheel rolling radius; m F is the hub motor braking force.Since the hub motor has its own peak torque limit, the required braking force for the braking force distribution by fuzzy identification cannot be dispatched to the hub motor unconditionally [17].In order to extend the service life of the hub motor, the generator torque should be two the smaller the braking torque of the brake.The peak torque of the hub motor at any speed can be measured by the test of the motor peak torque data set [14]. Figure 8 is the motor torque calculation of the module diagram

Suspension Model
The suspension of the car suspension system is non-linear and is multi-degree free.As shown in Figure 9.The degree of freedom will be more difficult to establish the model, so this article on the car suspension system two degree of freedom 1/4 model for simplified analysis [18] [19] [20].To get a simplified model, the following assumptions were made: 1) body and frame are rigid, and in the  course of the study to ignore the body suspension of the elasticity and damping; 2) the car on the road evenly and always with its contact, no impact occurred; 3) the stiffness of the suspension and the stiffness of the tire are linear functions of displacement, and ignore the tire damping; 4) The input function for road displacement can be assumed to act at the center of the tire and ground contact point.
Figure 10 2 m is the quality of the body, B is a damper, 1 m is the frame, 1 k represents the elastic element in the tire system.
By the Newton's second law we can see that the kinematic equation of the suspension system is ( ) ( ) ( ) By introducing some variables where, 1 Build the above constructed model as shown in Figure 11.

Conclusion
With MATLAB/Simulink software, you can easily and easily build a variety of automotive dynamics models, reduce the time in the programming, and quickly
the speed of the ve- hicle, unit m/s; o n is the reference space frequency, the reference space frequency under the road power density, which is the road roughness coefficient, unit m 3 ; ( ) w t is the Gaussian white noise with mean zero[8].

Figure 4 .
Figure 4.The force analysis of vehicle in the XY plane.

Figure 6 .Figure 7 .
Figure 6.Illustrates the "magic formula" in the parameters of the tire characteristic curve.

Figure 8 .
Figure 8. Module diagram of motor torque calculation.

Figure 10 .
Figure 10.Simplify the system model.