Integrated self-powered active suspension actuator control system and control method thereof

文档序号:1497661 发布日期:2020-02-07 浏览:14次 中文

阅读说明:本技术 一种集成式自供能主动悬架作动器控制系统及其控制方法 (Integrated self-powered active suspension actuator control system and control method thereof ) 是由 梁津 李旭林 王印 马相飞 任全 孟令建 赵鹏飞 王俊娇 王星斌 曹洲 于 2019-10-14 设计创作,主要内容包括:本发明公开了一种集成式自供能主动悬架作动器控制系统,包括主动控制系统和能量回收系统,主动控制系统包括作动器、车速传感器、路面不平度传感器、非簧载质量位移传感器、簧载质量位移传感器、作动器速度传感器、路面不平度探测器,能量回收系统包括驱动模块、整流滤波模块、升压模块、第一MOS开关模块、第二MOS开关模块、能量存储模块、OCV检测模块、控制模块和控制电路供电模块。本发明应用于重型商用车,集成了空气弹簧、螺旋弹簧与直线电机,可有效减小悬架系统刚度,减小悬架系统偏频,并可以进行主动控制,根据行使工况调整悬架系统阻尼,提高平顺性与操纵稳定性。另外,可进行振动能量回收,实现系统自供能,不需要消耗外部能量。(The invention discloses an integrated self-powered active suspension actuator control system which comprises an active control system and an energy recovery system, wherein the active control system comprises an actuator, a vehicle speed sensor, a road surface irregularity sensor, an unsprung mass displacement sensor, a sprung mass displacement sensor, an actuator speed sensor and a road surface irregularity detector, and the energy recovery system comprises a driving module, a rectifying and filtering module, a boosting module, a first MOS (metal oxide semiconductor) switch module, a second MOS switch module, an energy storage module, an OCV (open circuit voltage) detection module, a control module and a control circuit power supply module. The invention is applied to heavy commercial vehicles, integrates the air spring, the spiral spring and the linear motor, can effectively reduce the rigidity of a suspension system, reduce the offset frequency of the suspension system, can carry out active control, adjusts the damping of the suspension system according to the driving working condition, and improves the smoothness and the operation stability. In addition, vibration energy recovery can be carried out, self-energy supply of the system is realized, and external energy does not need to be consumed.)

1. An integrated self-powered active suspension actuator control system, comprising: the pneumatic actuator comprises an actuator (23), the actuator (23) comprises an actuator mounting support (2), a first actuator mounting hole (1) is formed in the actuator mounting support (2), the actuator mounting support (2) is fixed to the top end of an actuator outer shell (10), a top fixed connection limiting block fixing seat (4) in the actuator outer shell (10), a rubber limiting block (3) used for limiting the maximum compression amount of the actuator (23) is installed on the limiting block fixing seat (4), an air spring bag skin (5) installed on the actuator outer shell (10) is arranged on the periphery of the rubber limiting block (3), the lower end of the air spring bag skin (5) is connected with an air spring piston (9) in an interference fit manner, a limiting buffer seat support (8) is installed inside the air spring piston (9), the bottom of the limiting buffer seat support (8) is fixedly connected with the bottom of the air spring piston (9), the top fixed connection limit buffer seat (6) of limit buffer seat support (8), limit buffer seat (6) are located the outside of air spring piston (9), the inside of air spring bag skin (5), the top of limit buffer seat support (8) is equipped with limit buffer seat support inner chamber air vent (7), air spring piston (9) bottom is connected with coil spring (11), coil spring (11) lower part links to each other with actuator shell body (10), first permanent magnet fixing base (12) is installed to air spring piston (9) lower part, the bottom of air spring piston (9) is equipped with air spring additional chamber air vent (13), first permanent magnet fixing base (12) bottom is connected with secondary iron core (14), S utmost point permanent magnet array (15) and N utmost point permanent magnet array (16) are installed to secondary iron core (14) outside, install primary iron core (17) between actuator shell body (10) and S utmost point permanent magnet array (15), N utmost point permanent magnet array (16), the inside winding of primary iron core (17) has winding coil (18), second permanent magnet fixing base (19) are installed to secondary iron core (14) lower part, first permanent magnet fixing base (12) and second permanent magnet fixing base (19) all link to each other with S utmost point permanent magnet array (15), N utmost point permanent magnet array (16), secondary iron core (14) lower part is equipped with second actuator mounting hole (20), actuator shell body (10) top is equipped with gas pressure sensor (21) and blow vent (22).

2. The integrated self-powered active suspension actuator control system of claim 1, wherein: the actuator control system further comprises an active control system and an energy recovery system, wherein the active control system comprises a vehicle speed sensor (29), a road surface irregularity detector (30), an unsprung mass displacement sensor (31), a sprung mass displacement sensor (32), an actuator speed sensor (33), a road surface irregularity displacement sensor (34), a vehicle-mounted computer (35) and a control module (36); the energy recovery system comprises a driving module (24), a first MOS switch module (25), a rectifying and filtering module (26), a boosting module (27), a second MOS switch module (28), a 0CV detection module (37), a control circuit power supply module (38) and an energy storage module (39).

3. The integrated self-powered active suspension actuator control system of claim 1, wherein: the actuator mounting bracket (2) is connected with the whole vehicle frame in a bolt connection mode through the first actuator mounting hole (1); the upper plane of the actuator outer shell (10) is attached to the lower plane of the whole vehicle frame; and the secondary iron core (14) is connected with the whole vehicle axle in a bolt connection mode through a second actuator mounting hole (20).

4. The integrated self-powered active suspension actuator control system of claim 1, wherein: the air spring piston (9), the actuator outer shell (10), the first permanent magnet fixing seat (12), the secondary iron core (14) and the second permanent magnet fixing seat (19) are all made of magnetism isolating materials, and the primary iron core (17) is made of magnetic conducting materials.

5. The integrated self-powered active suspension actuator control system of claim 1, wherein: the air spring bag skin (5) is internally provided with a main air chamber which is communicated with a middle air chamber formed by a limit buffer seat support (8) through an air vent hole (7) in an inner cavity of the limit buffer seat support, and the middle air chamber formed by the limit buffer seat support (8) is communicated with an additional air chamber formed by a secondary iron core (14) through an air spring additional cavity air vent hole (13); a gap exists between the top of the air spring piston (9) and the limiting buffer seat (6), and a main air chamber inside the air spring bag skin (5) is communicated with an inner air chamber inside the air spring piston (9).

6. The integrated self-powered active suspension actuator control system of claim 1, wherein: the winding coil (18) is a double-layer winding, namely two coils are arranged in each tooth slot, the corresponding difference of each phase voltage is 120 degrees, and the distribution sequence of each phase voltage is as follows: c-, A +; a-, B +; b-, C +; c-, A + … … where A, B, C represents three phases and +, -represents current direction.

7. The integrated self-powered active suspension actuator control system of claim 2, wherein:

the control module (36) is responsible for receiving sensor signals and performing circuit control;

the vehicle-mounted computer (35) is responsible for calculating a control algorithm;

the vehicle speed sensor (29) is responsible for periodically collecting vehicle speed data;

the road surface unevenness detector (30) is arranged at the front part of the vehicle and is responsible for detecting the road surface unevenness of a road at a certain distance in front;

the unsprung mass displacement sensor (31) is responsible for periodically collecting unsprung mass displacement;

the sprung mass displacement sensor (32) is responsible for periodically collecting the sprung mass displacement;

the actuator speed sensor (33) is responsible for periodically acquiring the movement speed of the actuator (23);

the road surface irregularity displacement sensor (34) is responsible for periodically collecting road surface irregularity information of the ground where the front axle is located;

the gas pressure sensor (21) is responsible for periodically collecting the gas pressure in the actuator;

the driving module (24) is responsible for controlling and driving the motor according to the control signal;

the first MOS switch module (25) and the second MOS switch module (28) are responsible for switching an energy recovery mode and an active control mode in a hardware circuit, and when the first MOS switch module (25) is connected and disconnected with the second MOS switch module (28), the actuator (23) is in the active control mode; when the first MOS switch module (25) is off and the second MOS switch module (28) is on, the actuator (23) is in an energy recovery mode;

the rectification filter module (26) is responsible for converting the direct-current voltage stored by the energy storage module (39) into the alternating-current voltage required by the motor, or converting the alternating-current voltage generated by the motor into the direct-current voltage which can be stored by the energy storage module (39);

the boosting module (27) is responsible for boosting the electric energy generated by the actuator (23) in the energy recovery mode so as to store the electric energy by the energy storage module (39);

the OCV detection module (37) is responsible for detecting an open circuit voltage of the energy storage module (39) and transmitting an open circuit voltage signal to the control module (36);

the control circuit power supply module (38) is responsible for supplying required electric energy to the control module (36);

the energy storage module (39) is responsible for storing electrical energy.

8. A control method of an integrated self-powered active suspension actuator control system is characterized by comprising the following steps: the method comprises the following steps:

step I, detection data acquisition, which comprises the following specific processes: when the wheels are excited by road unevenness to generate vertical acting force, the lower end of the secondary iron core (14) is connected with the axle through the second actuator mounting hole (20), the acting force is transmitted to the actuator (23), and at the moment, the S pole permanent magnet array (15) and the N pole permanent magnet array (16) which are fixed on the secondary iron core (14) and the air spring piston (9) which is fixedly connected with the secondary iron core (14) vibrate vertically; the actuator mounting support (2) is fixedly connected with an actuator outer shell (10), the actuator mounting support (2) is fixedly connected with a frame through a first actuator mounting hole (1), a primary iron core (17) is fixedly connected with the actuator outer shell (10), a winding coil (18) is wound in the primary iron core (17), an S pole permanent magnet array (15) and an N pole permanent magnet array (16) can generate relative displacement with the winding coil (18), a conductor cuts a magnetic induction line to generate electric energy, and the electric energy is in an energy recovery state; after the winding coil (18) is electrified, a traveling wave magnetic field is generated and moves along the axial direction of the actuator (23), when the running speeds of the secondary iron core (14) and the primary iron core (17) are different, a slip rate is generated, so that electromagnetic force in the linear direction is generated, the electromagnetic force is the active control force of the actuator (23), and the active control force is adjusted in real time according to different running conditions of a vehicle, so that the actuator (23) can be kept in the optimal damping force state at any time; the expression of the slip ratio isWhen the slip ratio is between 0 and 1, the actuator (23) is in a normal working state; when the slip ratio is more than 1, the moving direction of the secondary iron core (14) is opposite to the traveling magnetic field directionWhen the actuator (23) is in a braking state; when the slip ratio is less than 0, the movement speed of the secondary iron core (14) is larger than the speed of the traveling wave magnetic field, and the actuator (23) is in an energy recovery state;

the vehicle speed sensor (29) detects the vehicle speed in real time, a road unevenness detector (30) arranged at the front part of the vehicle detects the road unevenness of a road at a certain distance in front, and the road input displacement x at the front wheel of the vehicle after a certain time delay is obtained0f(ii) a The road surface irregularity displacement sensor (34) detects the road surface irregularity of the front wheels of the vehicle in real time to obtain the road surface input displacement x at the rear wheels of the vehicle after a certain time delay0r(ii) a The unsprung mass displacement sensor (31) detects unsprung mass displacement in real time, the sprung mass displacement sensor (32) detects sprung mass displacement in real time, the actuator speed sensor (33) detects the motion speed of a suspension actuator in real time, and the control module (36) respectively obtains vehicle speed V and unsprung mass displacement z of a front wheel of a vehicle, which are obtained through periodic samplingwfUnsprung mass displacement z of the rear wheel of the vehiclewrVehicle front wheel sprung mass displacement zbfVehicle rear wheel sprung mass displacement zbrAnd the moving speed u of an actuator (23) mounted on the front wheel of the vehicleafAnd the moving speed u of an actuator (23) mounted on the rear wheel of the vehiclearAnd transmitted to the vehicle computer (35);

step II, establishing a system state equation, which comprises the following specific processes:

step 201, the vehicle-mounted computer (35) establishes a half vehicle model motion equation:

Figure FDA0002232833130000051

Figure FDA0002232833130000052

Figure FDA0002232833130000053

Figure FDA0002232833130000054

wherein the content of the first and second substances,

Figure FDA0002232833130000055

when the pitch angle thetabSmaller, there is

Figure FDA0002232833130000059

The equation of motion of the 4-degree-of-freedom semi-vehicle model is expressed in a matrix form:

Figure FDA00022328331300000510

in the formula, the quality matrix [ M ] is:

[M]=diag(mb,Ib,mwf,mwr)

the damping matrix [ C ] is:

the stiffness matrix [ K ] is:

Figure FDA0002232833130000062

tire stiffness matrix [ K ]t]Comprises the following steps:

Figure FDA0002232833130000063

road surface input vector x0Comprises the following steps:

x0=[x0f,x0r]T

the system input vector x is:

x=[z θbzwfzwr]T

step 202, when the value of the suspension mass distribution coefficient epsilon is close to 1, 2-degree-of-freedom models can be respectively established for the front axis and the rear axis, and at the moment, a state variable is taken as follows:

Figure FDA0002232833130000064

wherein z iswIs the unsprung mass displacement of the vehicle, and zwIs taken as zwfOr zwr,zbIs the sprung mass displacement of the vehicle, and zbIs taken as zbfOr zbr(ii) a If the road speed input is white noise, thenThe system state equation at this time is:

Figure FDA0002232833130000071

wherein u is the main power of the actuator, and the value of u is the main power u of the front wheel actuatorfOr rear wheel actuator main power ur(ii) a A is a constant matrix relating to state variables of the vehicle model and

Figure FDA0002232833130000072

Step III, the working mode switching rule of the integrated self-powered active suspension actuator control system is as follows:

step 301, defining an energy recovery mode of the integrated self-powered active suspension actuator control system as MeDefining the active control mode as Mc

Step 302, set moveDeflection threshold fdyAnd duration threshold tyWherein dynamic deflection threshold fdyThe mean square root value of the dynamic deflection of the suspension under the condition of a good road surface and a low vehicle speed is smaller than that of a specific vehicle type; when the actual road surface dynamic deflection fd<fdyWhen the actual road surface dynamic deflection f is larger than the actual road surface dynamic deflection f, the vehicle is shown to run on a better road surfaced>fdyWhen, it indicates that the vehicle is traveling on a poor road surface; duration threshold tyThe influence on the working mode switching frequency of the integrated self-powered active suspension actuator control system is great, if the time threshold t is setyToo small, which may result in frequent switching of the operating mode of the integrated self-powered active suspension actuator control system, if the time threshold t is setyToo large, which may result in slow response of the integrated self-powered active suspension actuator control system to switching modes, time threshold tyThe specific numerical value of the road is determined according to the actual vehicle type and the common road condition;

step 303, the OCV detection module (37) collects the open-circuit voltage of the energy storage module (39), the control module (36) adopts an AH integration method to estimate the SOC, and the rated capacity of the energy storage module (39) is set to be QeReal-time remaining capacity of QcInitial state capacity of Q0η is coulombic efficiency, then at tsReal-time remaining capacity at time QcThe state, using the integral operation, can be expressed as:

Figure FDA0002232833130000081

the remaining energy of the energy storage module (39) is set to 3 levels, namely a state of charge (SOC)1Normal state of charge SOC2State of sufficient charge SOC3

Step 304, the main function of the integrated self-powered active suspension actuator control system is to ensure the smoothness and the handling stability of the vehicleSelf-energizing features of an integrated self-energized active suspension actuator control system are integrated such that when the energy storage module (39) is at SOC2、SOC3In the state or when the running road surface condition is poor and the whole vehicle has large vibration, the vehicle is in McAn active control mode; the energy storage module (39) is in SOC1The state or the running road surface condition is good, and when the whole vehicle vibrates less, the vehicle is in MeAn energy recovery mode; the mode switching rule is therefore:

step IV, performing pre-aiming optimal control on the unevenness of the front axle road surface of the integrated self-powered active suspension actuator control system, wherein the specific process is as follows:

step 401, the vehicle computer (35) defines the vehicle speed v at any moment iiLower, the pre-aiming distance LiWith the speed v of the vehicleiIs the aiming time t of the ith aiming pointiIs formulated as:

wherein, the value of i is a natural number between 1 and n, and n is the total sampling times of the road surface unevenness detector (30);

step 402, the vehicle-mounted computer (35) inputs the pavement input displacement x at the ith pre-aiming pointiInput displacement x from road surface at front wheel of vehicle0fThe Laplace functional relation is expressed as:

Figure FDA0002232833130000093

and e is-tisThe approximate calculation of (a) is expressed as:

then, the second order Pade approximation is taken and is used in advanceWhen the aiming time is not 0, the numerator and denominator of the formula (F14) are divided by the numerator and denominator simultaneously

Figure FDA0002232833130000095

Figure FDA0002232833130000101

wherein, ai0And ai1Are all a coefficient of months

Figure FDA0002232833130000102

The formula (F12) is then rewritten as:

[x0f(s)-zi(s)][ai0+ai1s+s2]=-2ai1sxi(s) (F16)

step 403, the vehicle computer (35) defines x0f(s)-xi(s)=yi(s),b1=-2ai1Performing an inverse Laplace transform on the formula (F16) yields:

redefining the state variables as:

ηi1=yi0xt=yt

Figure FDA0002232833130000104

β0=0

β1=b1-ai1β0=b1

and rewriting the state equation of the ith pre-aiming point as follows:

Figure FDA0002232833130000106

wherein the content of the first and second substances,

Figure FDA0002232833130000107

the state equation of the n preview points is obtained as follows:

Figure FDA0002232833130000108

wherein the content of the first and second substances,

Figure FDA0002232833130000109

Figure FDA0002232833130000111

xp=[x1x2… xn]T

η=[η11η12η21η22… ηn1ηn2]T

step 404, the vehicle-mounted computer (35) defines the road surface input displacement variable at the wheel as:

x0f=Dηη+Eηxp(F20)

wherein D isηAnd EηAre all constant matrices and

Figure FDA0002232833130000113

step 405, substituting the formula (F20) into the formula (F8) to obtain a system state equation:

step 406, the vehicle computer (35) sends q1fDefined as the front axle spring mass displacement weighting coefficient representing ride comfort, q2fDefining a front axle suspension dynamic stroke weighting coefficient which influences the vehicle body posture and is related to the structural design and the arrangement, and q3fThe front axle tire dynamic displacement weighting coefficient representing the steering stability is defined, and the performance index of the LQR controller design in the T moment is defined as:

Figure FDA0002232833130000115

wherein z iswfFor unsprung mass displacement of the front axle, zbfFor the displacement of the front axle spring-loaded mass,

Figure FDA0002232833130000116

step 407, the onboard computer (35) brings the system state equation (F8) into the equation (F22) to obtain:

Figure FDA0002232833130000117

wherein the state variable gain matrix

Figure FDA0002232833130000121

Control variable gain matrix

Figure FDA0002232833130000122

Step 408, the vehicle computer (35) outputs an optimal control feedback gain matrix KfWherein

Figure FDA0002232833130000124

Step 409, the vehicle computer (35) calculates the formula Uf=-KfxfObtaining the optimal control force U of the front axlef

Step 410, the linear motor electromagnetic thrust expression is

Faf=BδfLfIf(F24)

Wherein B isδfIs the average air gap magnetic density, L, at which the moving coil is locatedfIs the total length of a linear motor moving coil winding conductor, IfThe coil current is led into the front axle actuator;

by passing a coil current I to the front-axle actuatorfSo that Faf=UfThe front axle actuator can output the optimal active control force of the front axle;

v, performing pre-aiming optimal optimization control on the rear axle base of the integrated self-powered active suspension actuator control system, wherein the specific process is as follows:

step 501, the road surface irregularity displacement sensor (34) collects the road surface irregularity impact on the front wheels of the vehicle, and the road surface irregularity input displacement x at the front wheels at the moment0fWill be delayed by a time of

Figure FDA0002232833130000127

Figure FDA0002232833130000131

step 502, let q1r、q2r、q3rThe vertical displacement of a vehicle body representing riding comfort, the suspension dynamic stroke which influences the posture of the vehicle body and is related to structural design and arrangement and the tire dynamic displacement representing operation stability are respectively controlled by 3 weighting coefficients, so that the performance index of the LQR controller design in the T moment is as follows:

Figure FDA0002232833130000132

the system state equation (F25) is brought into formula (F26)

Wherein the state variable gain matrix

Figure FDA0002232833130000134

Control variable gain matrixVariable correlation matrix

Figure FDA0002232833130000136

Step 503, using adaptive multi-forgetting factor recursive least square method to Qr、Rr、NrOptimizing, and setting a coefficient optimization matrix theta ═ QrRrNr]TSetting a covariance matrix Pi=106E, wherein E is a third-order identity matrix, and a data matrix phi is set:

Figure FDA0002232833130000137

calculate time T

Figure FDA0002232833130000138

Figure FDA0002232833130000141

Calculating the value of the variable forgetting factor:

Figure FDA0002232833130000142

where l is the step size, calculate the gain K of a single parameteri,k

Figure FDA0002232833130000143

Where i is 1, 2 … n, where n is the number of identification parameters, calculating gain matrix Lk

Figure FDA0002232833130000144

And (3) performing parameter estimation:

Figure FDA0002232833130000145

updating the covariance matrix:

Figure FDA0002232833130000146

step 504, repeating step 503 until the identification is finished, and obtaining the optimal performance index J meeting the expectationyQ of (2)r、Rr、NrFurther, 3 weighting coefficients q are obtained1r、q2r、q3r

505, outputting an optimal control feedback gain matrix K by the vehicle-mounted computer (35)rWherein

Figure FDA0002232833130000147

Step 506, the vehicle-mounted computer (35) according to a formula Ur=-KrxrCalculating the optimal control force U of the rear axler

Step 507, the linear motor electromagnetic thrust expression is

Far=BδrLrIr(F35)

Wherein B isδrIs the average air gap magnetic density, L, at which the moving coil is locatedrIs the total length of a linear motor moving coil winding conductor, IrThe coil current is led into the front axle actuator;

by passing a coil current I to the rear-axle actuatorrSo that Far=UrThe rear axle actuator can output the optimal active control force of the rear axle;

step VI, energy recovery of the integrated self-powered active suspension actuator control system, which comprises the following specific processes: when the integrated self-powered active suspension actuator control system is in MeIn the energy recovery mode, the first MOS switch module (25) is disconnected, the second MOS switch module (28) is connected, the wheel is excited by road surface unevenness, when vertical acting force is generated, the lower end of the secondary iron core (14) is connected with the axle through the second actuator mounting hole (20), the acting force is transmitted to the actuator (23), and the S pole permanent magnet array (15) and the N pole permanent magnet array (16) which are fixed on the secondary iron core (14) and the air spring piston (9) fixedly connected with the secondary iron core (14) vibrate vertically; the actuator mounting supportThe frame (2) is fixedly connected with an actuator outer shell (10), the actuator mounting bracket (2) is fixedly connected with a frame through a first actuator mounting hole (1), the primary iron core (17) is fixedly connected with the actuator outer shell (10), and a winding coil (18) is wound in the primary iron core (17); the S pole permanent magnet array (15) and the N pole permanent magnet array (16) can generate relative displacement with the winding coil (18) to cut the magnetic induction line to generate electric energy, and as the unevenness excitation of the road surface is changed and the movement speed of the actuator (23) is also changed, the speed of the winding coil (18) cutting the magnetic induction line is also changed, so that the energy feedback voltage is also changed, the energy feedback voltage needs to be processed by the rectifier filter module (26), and the energy feedback voltage is boosted to a voltage range which can be received by the energy storage module (39) through the boosting module (27), and the vibration energy recovery is completed.

Technical Field

The invention belongs to the technical field of suspension systems of heavy commercial vehicles, and particularly relates to an integrated self-powered active suspension actuator control system and a control method thereof.

Background

When the vehicle runs, the sprung mass and the unsprung mass of the vehicle can generate relative displacement under the influence of the unevenness of the road surface, so that the vehicle vibrates. The suspension system is an important part playing a role in damping vibration in an automobile chassis, is positioned between a frame and an axle, and has the functions of buffering and attenuating impact and vibration caused by road unevenness and transmitting various forces and moments acting between wheels and the axle to influence the smoothness and the operation stability of the automobile. At present, the power performance and the like of heavy commercial vehicles are better and better, and the market also puts forward higher requirements on the transportation timeliness of the heavy commercial vehicles, so when the heavy commercial vehicles with large full load and mass are excited by road unevenness at a high speed, the vehicles are strongly impacted and generate severe vibration, and the condition is very unfavorable for the cargo integrity, the vehicle smoothness and the operation stability. At present, the conventional passive suspension is generally adopted by the heavy commercial vehicle, and after the parameters such as rigidity, damping and the like of the heavy commercial vehicle are determined during design, the parameters are fixed and unchanged when the vehicle runs on different road conditions, so that the vibration damping effect of the conventional passive suspension is greatly limited. The active suspension can generate controllable active force, so that different optimal active force can be obtained when the vehicle runs on different road surfaces, and the smoothness and the operation stability are improved. However, the energy consumption problem of the active suspension is always one of the main factors restricting the market popularization, and the integrated self-powered active suspension actuator provided by the invention has a self-powered function and does not need to consume external energy. For example, a ball screw type shock absorber related in chinese patent publication No. 201410650452.4, "ball screw type inerter with mechanically variable inerter coefficient", is mainly to alleviate the instantaneous impact on the road surface through a mechanical structure, but once manufactured, various parameters thereof are determined, and the parameters cannot be adjusted according to different road conditions and different vehicle types, so that there is a certain limitation. Chinese patent publication No. CN 102261417B, "additional air chamber of constant internal pressure air spring with continuously changing volume and working method", discloses an additional air chamber of constant internal pressure air spring with continuously changing volume and working method, wherein a partition plate for air spring is divided into 3 air chambers, and a linear motor is used to control the movement of the partition plate, so that the volumes of three chambers of the air chamber cylinder of the air spring are changed, thereby controlling the stiffness of the suspension system, but the method can only control the stiffness of the suspension system, and does not solve the energy consumption problem.

Therefore, at present, an active suspension actuator integrating an air spring, a spiral spring and a linear motor of a heavy commercial vehicle with large full load mass does not exist, the active suspension actuator can be in an active control mode to carry out vibration control or energy recovery mode to carry out vibration energy recovery according to different driving conditions of the heavy commercial vehicle, external energy does not need to be consumed, and the active suspension actuator has a self-powered function.

Disclosure of Invention

The invention aims to solve the technical problem of providing an integrated self-powered active suspension actuator control system and a control method thereof aiming at the defects in the prior art. The suspension system is applied to heavy commercial vehicles, integrates the air spring, the spiral spring and the linear motor, can effectively reduce the rigidity of the suspension system, reduces the offset frequency of the suspension system, can be actively controlled, adjusts the damping of the suspension system according to the driving working condition, and improves the smoothness and the operation stability. In addition, vibration energy recovery can be carried out, self-energy supply of the system is realized, and external energy does not need to be consumed.

In order to solve the technical problems, the invention adopts the technical scheme that:

an integrated self-powered active suspension actuator control system comprises an actuator, wherein the actuator comprises an actuator mounting support, a first actuator mounting hole is formed in the actuator mounting support, the actuator mounting support is fixed at the top end of an actuator outer shell, a top fixed connection limiting block fixing seat in the actuator outer shell is provided with a rubber limiting block for limiting the maximum compression amount of the actuator, an air spring bag skin arranged on the actuator outer shell is arranged on the periphery of the rubber limiting block, the lower end of the air spring bag skin is connected with an air spring piston in an interference fit manner, a limiting buffer seat support is arranged in the air spring piston, the bottom of the limiting buffer seat support is fixedly connected with the bottom of the air spring piston, the top of the limiting buffer seat support is fixedly connected with a limiting buffer seat, the air spring piston comprises an air spring piston, a limiting buffer seat, an air spring bag cover, a limiting buffer seat support, an air spring piston, a coil spring, a primary iron core, a winding coil, a secondary iron core, a first permanent magnet fixing seat, a second permanent magnet fixing seat, a first permanent magnet fixing seat, an S-pole permanent magnet array, an N-pole permanent magnet array, a winding coil, a second permanent magnet fixing seat, an air spring support inner cavity air vent, a spiral spring, a first permanent magnet fixing seat, a second permanent magnet fixing seat, an air spring bag cover and a first permanent magnet fixing seat, wherein the limiting buffer seat is located outside the air spring piston and inside the air spring bag cover, the top of the limiting buffer seat support is provided with the limiting buffer seat support inner cavity air vent, the, The N pole permanent magnet arrays are connected, a second actuator mounting hole is formed in the lower portion of the secondary iron core, and a gas pressure sensor and a vent are arranged at the top of the outer shell of the actuator.

Preferably, the actuator control system further comprises an active control system and an energy recovery system, wherein the active control system comprises a vehicle speed sensor, a road surface unevenness detector, an unsprung mass displacement sensor, a sprung mass displacement sensor, an actuator speed sensor, a road surface unevenness displacement sensor, a vehicle-mounted computer and a control module; the energy recovery system comprises a driving module, a first MOS switch module, a rectifying and filtering module, a boosting module, a second MOS switch module, an OCV detection module, a control circuit power supply module and an energy storage module.

Preferably, the actuator mounting bracket is connected with the whole vehicle frame through a first actuator mounting hole in a bolt connection mode; the upper plane of the actuator shell is attached to the lower plane of the whole vehicle frame; and the secondary iron core is connected with the whole vehicle axle in a bolt connection mode through a second actuator mounting hole.

Preferably, the air spring piston, the actuator outer shell, the first permanent magnet fixing seat, the secondary iron core and the second permanent magnet fixing seat are made of magnetism isolating materials, and the primary iron core is made of magnetic conducting materials.

Preferably, the air spring bag skin is internally provided with a main air chamber which is communicated with a middle air chamber formed by the limit buffer seat bracket through an air vent hole in an inner cavity of the limit buffer seat bracket, and the middle air chamber formed by the limit buffer seat bracket is communicated with an additional air chamber formed by the secondary iron core through an air vent hole in an additional cavity of the air spring; a gap exists between the top of the air spring piston and the limiting buffer seat, and a main air chamber inside the air spring bag skin is communicated with an inner air chamber inside the air spring piston.

Preferably, the winding coil is a double-layer winding, that is, there are two coils in each slot, each phase voltage has a corresponding difference of 120 °, and the distribution sequence of each phase voltage is: c-, A +; a-, B +; b-, C +; c-, A + … … where A, B, C represents three phases and +, -represents current direction.

Furthermore, the control module is responsible for receiving the sensor signal and performing circuit control;

the vehicle-mounted computer is responsible for calculating a control algorithm;

the vehicle speed sensor is responsible for periodically collecting vehicle speed data;

the road surface irregularity detector is arranged at the front part of the vehicle and is responsible for detecting the road surface irregularity of a road at a certain distance ahead;

the unsprung mass displacement sensor is responsible for periodically collecting unsprung mass displacement;

the sprung mass displacement sensor is responsible for periodically collecting sprung mass displacement;

the actuator speed sensor is used for periodically collecting the movement speed of the actuator;

the road surface irregularity displacement sensor is responsible for periodically collecting road surface irregularity information of the ground where the front axle is located;

the gas pressure sensor is used for periodically collecting the gas pressure in the actuator;

the driving module is used for controlling and driving the motor according to the control signal;

the first MOS switch module and the second MOS switch module are responsible for switching an energy recovery mode and an active control mode in a hardware circuit, and when the first MOS switch module is connected and disconnected with the second MOS switch module, the actuator is in the active control mode; when the first MOS switch module is disconnected and the second MOS switch module is connected, the actuator is in an energy recovery mode;

the rectification filter module is used for converting the direct-current voltage stored by the energy storage module into alternating-current voltage required by the motor or converting the alternating-current voltage generated by the motor into direct-current voltage capable of being stored by the energy storage module;

the boosting module is responsible for boosting the electric energy generated by the actuator in the energy recovery mode so as to store the electric energy by the energy storage module;

the OCV detection module is used for detecting the open-circuit voltage of the energy storage module and transmitting an open-circuit voltage signal to the control module;

the control circuit power supply module is responsible for supplying required electric energy to the control module;

the energy storage module is responsible for storing electrical energy.

A method of controlling an integrated self-powered active suspension actuator control system, comprising the steps of:

step I, detection data acquisition, which comprises the following specific processes: when the wheels are excited by road surface unevenness to generate vertical acting force, the lower end of the secondary iron core is connected with the axle through a second actuator mounting hole, the acting force is transmitted to the actuator, and the S pole permanent magnet array, the N pole permanent magnet array and the air spring piston fixedly connected with the secondary iron core which are fixed on the secondary iron core vibrate vertically; the actuator mounting bracket is fixedly connected with the actuator outer shell, the actuator mounting bracket is fixedly connected with the frame through a first actuator mounting hole, the primary iron core is fixedly connected with the actuator outer shell, a winding coil is wound in the primary iron core, the S-pole permanent magnet array and the N-pole permanent magnet array generate relative displacement with the winding coil, a conductor cuts a magnetic induction line to generate electric energy, and the electric energy is in an energy recovery state; after the winding coil is electrified, a traveling wave magnetic field is generated and moves along the axial direction of the actuator, when the running speeds of the secondary iron core and the primary iron core are different, the slip rate is generated, so that the electromagnetic force in the linear direction is generated, the electromagnetic force is the active control force of the actuator, and the active control force is adjusted in real time according to different running working conditions of the vehicle, so that the actuator can be kept in the optimal damping force state at any time; the expression of the slip ratio is

Figure BDA0002232833140000051

When the slip ratio is between 0 and 1, the actuator is in a normal working state; when the slip ratio is larger than 1, the movement direction of the secondary iron core is opposite to the direction of the traveling wave magnetic field, and the actuator is in a braking state at the moment; when the slip ratio is less than 0, the movement speed of the secondary iron core is greater than the speed of the traveling wave magnetic field, and the actuator is in an energy recovery state;

the vehicle speed sensorThe device detects the speed of the vehicle in real time, a road unevenness detector arranged at the front part of the vehicle detects the road unevenness of a road at a certain distance in front, and the road input displacement x at the front wheel of the vehicle after a certain time delay is obtained0f(ii) a The road surface irregularity displacement sensor detects the road surface irregularity of the front wheels of the vehicle in real time to obtain the road surface input displacement x of the rear wheels of the vehicle after a certain time delay0r(ii) a The unsprung mass displacement sensor detects unsprung mass displacement in real time, the sprung mass displacement sensor detects sprung mass displacement in real time, the actuator speed sensor detects the motion speed of the suspension actuator in real time, and the control module respectively obtains vehicle speed V and unsprung mass displacement z of a front wheel of the vehicle, which are obtained by periodic samplingwfUnsprung mass displacement z of the rear wheel of the vehiclewrVehicle front wheel sprung mass displacement zbfVehicle rear wheel sprung mass displacement zbrAnd a moving speed u of an actuator mounted on a front wheel of the vehicleafAnd the speed of movement u of an actuator mounted on the rear wheel of the vehiclearAnd transmitting to the vehicle-mounted computer;

step II, establishing a system state equation, wherein the specific process is as follows:

step 201, the vehicle-mounted computer establishes a half vehicle model equation of motion:

Figure BDA0002232833140000061

Figure BDA0002232833140000063

Figure BDA0002232833140000064

wherein the content of the first and second substances,vertical acceleration of the sprung mass of the front wheel of the vehicle,

Figure BDA0002232833140000066

The vertical acceleration of the sprung mass of the rear wheel of the vehicle, the displacement of the sprung mass, and mbIs a sprung mass, IbIs the moment of inertia, theta, of the vehicle bodybIs the pitch angle, a is the distance from the front axle to the center of mass of the vehicle, b is the distance from the rear axle to the center of mass of the vehicle, ksfFor the stiffness, k, of the front active suspension actuatorsrFor the stiffness, m, of the rear active suspension actuatorwfFor unsprung mass, m, of the front wheels of the vehiclewrIs the unsprung mass of the rear wheel of the vehicle,

Figure BDA0002232833140000067

Is the vertical acceleration of the unsprung mass of the front wheel of the vehicle,

Figure BDA0002232833140000068

For unsprung mass vertical acceleration, k, of the rear wheels of the vehicletfAs the rigidity, k, of the front wheels of the vehicletrAs the stiffness, x, of the rear wheels of the vehicle0fInputting displacement, x, to the road surface of the front wheel0rInputting displacement for the road surface of the rear wheel;

when the pitch angle thetabSmaller, there is

Figure BDA0002232833140000069

The equation of motion of the 4-degree-of-freedom semi-vehicle model is expressed in a matrix form:

Figure BDA00022328331400000610

in the formula, the quality matrix [ M ] is:

[M]=diag(mb,Ib,mwf,mwr)

the damping matrix [ C ] is:

Figure BDA0002232833140000071

the stiffness matrix [ K ] is:

tire stiffness matrix [ K ]t]Comprises the following steps:

road surface input vector x0Comprises the following steps:

x0=[x0f,x0r]T

the system input vector x is:

x=[z θbzwfzwr]T

step 202, when the value of the suspension mass distribution coefficient epsilon is close to 1, 2-degree-of-freedom models can be respectively established for the front axis and the rear axis, and at the moment, a state variable is taken as follows:

Figure BDA0002232833140000074

wherein z iswIs the unsprung mass displacement of the vehicle, and zwIs taken as zwfOr zwr,zbIs the sprung mass displacement of the vehicle, and zbIs taken as zbfOr zbr(ii) a If the road speed input is white noise, thenThe system state equation at this time is:

Figure BDA0002232833140000075

wherein u is the main power of the actuator, and the value of u is the main power u of the front wheel actuatorfOr rear wheel actuator main power ur(ii) a A is a constant matrix relating to state variables of the vehicle model and

Figure BDA0002232833140000081

Figure BDA0002232833140000082

b is a constant matrix relating to the variation of the speed of movement of the actuator (23) and

Figure BDA0002232833140000083

g is a constant matrix related to the input displacement variable of the road surface at the wheel and

Figure BDA0002232833140000084

csis damping of vehicle suspension and csIs taken as the front wheel actuator damping csfOr rear wheel actuator damping csr、mbIs the mass of the vehicle body, ksIs the suspension stiffness and ksIs taken as the front wheel actuator stiffness ksfOr rear wheel actuator stiffness ksr、ktIs tire stiffness and ktIs taken as the front wheel tire stiffness ktfOr rear wheel tire stiffness ktr、mwIs the unsprung mass of the vehicle and mwIs taken as the unsprung mass m of the front wheelwfOr unsprung mass m of the rear wheelwr

Step III, the working mode switching rule of the integrated self-powered active suspension actuator control system is as follows:

step 301, defining an energy recovery mode of the integrated self-powered active suspension actuator control system as MeDefining the active control mode as Mc

Step 302, setting a dynamic deflection threshold fdyAnd duration threshold tyWherein dynamic deflection threshold fdyThe mean square root value of the dynamic deflection of the suspension under the condition of a good road surface and a low vehicle speed is smaller than that of a specific vehicle type; when the actual road surface dynamic deflection fd<fdyWhen the actual road surface dynamic deflection f is larger than the actual road surface dynamic deflection f, the vehicle is shown to run on a better road surfaced>fdyWhen, it indicates that the vehicle is traveling on a poor road surface; duration threshold tyControl system for integrated self-powered active suspension actuatorThe frequency of switching the working modes is greatly influenced if the time threshold tyToo small, which may result in frequent switching of the operating mode of the integrated self-powered active suspension actuator control system, if the time threshold t is setyToo large, which may result in slow response of the integrated self-powered active suspension actuator control system to switching modes, time threshold tyThe specific numerical value of the road is determined according to the actual vehicle type and the common road condition;

303, collecting the open-circuit voltage of the energy storage module by the OCV detection module, estimating the SOC by the control module by adopting an AH integration method, and setting the rated capacity of the energy storage module to be QeReal-time remaining capacity of QcInitial state capacity of Q0η is coulombic efficiency, then at tsReal-time remaining capacity at time QcThe state, using the integral operation, can be expressed as:

Figure BDA0002232833140000091

the remaining capacity of the energy storage module is set to be divided into 3 grades, namely a power-lack state SOC1Normal state of charge SOC2State of sufficient charge SOC3

Step 304, the main function of the integrated self-powered active suspension actuator control system is to ensure the smoothness and the handling stability of the vehicle, and the self-powered characteristic of the integrated self-powered active suspension actuator control system is combined, so that when the energy storage module is in the SOC state2、SOC3In the state or when the running road surface condition is poor and the whole vehicle has large vibration, the vehicle is in McAn active control mode; energy storage module at SOC1The state or the running road surface condition is good, and when the whole vehicle vibrates less, the vehicle is in MeAn energy recovery mode; the mode switching rule is therefore:

Figure BDA0002232833140000093

step IV, performing pre-aiming optimal control on the unevenness of the front axle road surface of the integrated self-powered active suspension actuator control system, wherein the specific process is as follows:

step 401, the vehicle computer defines the vehicle speed v at any moment iiLower, the pre-aiming distance LiWith the speed v of the vehicleiIs the aiming time t of the ith aiming pointiIs formulated as:

Figure BDA0002232833140000101

wherein, the value of i is a natural number from 1 to n, and n is the total sampling times of the road surface unevenness detector;

step 402, the vehicle-mounted computer inputs the road surface displacement x at the ith pre-aiming pointiInput displacement x from road surface at front wheel of vehicle0fThe Laplace functional relation is expressed as:

Figure BDA0002232833140000102

and e is-tisThe approximate calculation of (a) is expressed as:

Figure BDA0002232833140000103

then, the second-order Pade approximate calculation is taken, and when the preview time is not 0, the numerator and the denominator of the formula (F14) are simultaneously divided byObtaining:

Figure BDA0002232833140000105

wherein, ai0And ai1Are all coefficients and

Figure BDA0002232833140000106

the formula (F12) is then rewritten as:

[x0f(s)-zi(s)][ai0+ai1s+s2]=-2ai1sxi(s) (F16)

step 403, the vehicle computer defines x0f(s)-xi(s)=yi(s),b1=-2ai1Performing an inverse Laplace transform on the formula (F16) yields:

redefining the state variables as:

ηi1=yi0xt=yt

Figure BDA0002232833140000108

β0=0

β1=b1-ai1β0=b1

Figure BDA0002232833140000111

and rewriting the state equation of the ith pre-aiming point as follows:

Figure BDA0002232833140000112

wherein the content of the first and second substances,

the state equation of the n preview points is obtained as follows:

wherein the content of the first and second substances,

Figure BDA0002232833140000115

Figure BDA0002232833140000116

xp=[x1x2… xn]T,

η=[η11η12η21η22… ηn1ηn2]T

step 404, the vehicle-mounted computer defines the road surface input displacement variable at the wheel as:

x0f=Dηη+Eηxp(F20)

wherein D isηAnd EηAre all constant matrices and

Figure BDA0002232833140000117

Figure BDA0002232833140000118

Figure BDA0002232833140000119

λiis a weighting coefficient of the ith preview point and12+…+λnn is the total number of the preview points as 1;

step 405, substituting the formula (F20) into the formula (F8) to obtain a system state equation:

Figure BDA0002232833140000121

step 406, the vehicle computer will q1fDefined as the front axle spring mass displacement weighting coefficient representing ride comfort, q2fDefining a front axle suspension dynamic stroke weighting coefficient which influences the vehicle body posture and is related to the structural design and the arrangement, and q3fDefining front axle tyre dynamic displacement weighting coefficient representing steering stability and defining LQR controller in T momentThe designed performance indexes are as follows:

Figure BDA0002232833140000122

wherein z iswfFor unsprung mass displacement of the front axle, zbfFor the displacement of the front axle spring-loaded mass,

Figure BDA0002232833140000123

vertical acceleration of the front axle spring load mass;

step 407, the onboard computer brings the system state equation (F8) into the equation (F22) to obtain:

wherein the state variable gain matrix

Control variable gain matrix

Figure BDA0002232833140000126

Variable correlation matrix

Figure BDA0002232833140000127

Step 408, the vehicle computer outputs an optimal control feedback gain matrix KfWherein

Figure BDA0002232833140000128

Figure BDA0002232833140000129

PcfSatisfy the Raccati equation

Figure BDA00022328331400001210

Step 409, the vehicle-mounted computer calculates according to a formula Uf=-KfxfObtaining the optimal control force U of the front axlef

Step 410, the linear motor electromagnetic thrust expression is

Faf=BδfLfIf(F24)

Wherein B isδfIs the average air gap magnetic density, L, at which the moving coil is locatedfIs the total length of a linear motor moving coil winding conductor, IfThe coil current is led into the front axle actuator;

by passing a coil current I to the front-axle actuatorfSo that Faf=UfThe front axle actuator can output the optimal active control force of the front axle;

step V, optimal pre-aiming optimization control of the rear axle base of the integrated self-powered active suspension actuator control system is carried out, and the specific process is as follows:

step 501, the road surface irregularity displacement sensor collects road surface irregularity impact on the front wheels of the vehicle, and at the moment, the road surface irregularity input displacement x at the front wheels is0fWill be delayed by a time of

Figure BDA0002232833140000131

Input displacement x equal to road surface unevenness of rear wheel0rWherein V is the running speed of the vehicle and L is the wheel base of the vehicle, thereby obtaining the road surface input x0rAnd obtaining a state equation of the rear suspension system:

Figure BDA0002232833140000132

step 502, let q1r、q2r、q3rThe vertical displacement of a vehicle body representing riding comfort, the suspension dynamic stroke which influences the posture of the vehicle body and is related to structural design and arrangement and the tire dynamic displacement representing operation stability are respectively controlled by 3 weighting coefficients, so that the performance index of the LQR controller design in the T moment is as follows:

Figure BDA0002232833140000133

the system state equation (F25) is brought into formula (F26)

Figure BDA0002232833140000134

Wherein the state variable gain matrix

Figure BDA0002232833140000141

Control variable gain matrix

Figure BDA0002232833140000142

Variable correlation matrix

Figure BDA0002232833140000143

Step 503, using adaptive multi-forgetting factor recursive least square method to Qr、Rr、NrOptimizing, and setting a coefficient optimization matrix theta ═ QrRrNr]TSetting a covariance matrix Pi=106E, wherein E is a third-order identity matrix, and a data matrix phi is set:

Figure BDA0002232833140000144

calculate time TSetting an expected optimal performance index JyAnd calculating an estimation error ek

Figure BDA0002232833140000146

Calculating the value of the variable forgetting factor:

where l is the step size, calculating a single parameterGain Ki,k

Figure BDA0002232833140000148

Where i is 1, 2 … n, where n is the number of identification parameters, calculating gain matrix Lk:

Figure BDA0002232833140000149

And (3) performing parameter estimation:

Figure BDA0002232833140000151

updating the covariance matrix:

Figure BDA0002232833140000152

step 504, repeating step 503 until the identification is finished, and obtaining the optimal performance index J meeting the expectationyQ of (2)r、Rr、NrFurther, 3 weighting coefficients q are obtained1r、q2r、q3r

505, outputting the optimal control feedback gain matrix K by the vehicle-mounted computerrWherein

Figure BDA0002232833140000153

Figure BDA0002232833140000154

PcrSatisfy the Raccati equation

Figure BDA0002232833140000155

Step 506, the vehicle-mounted computer performs calculation according to a formula Ur=-KrxrCalculating the optimal control force U of the rear axler

Step 507, the linear motor electromagnetic thrust expression is

Far=BδrLrIr(F35)

Wherein B isδrIs the average air gap magnetic density, L, at which the moving coil is locatedrIs the total length of a linear motor moving coil winding conductor, IrThe coil current is led into the front axle actuator;

by passing a coil current I to the rear-axle actuatorrSo that Far=UrThe rear axle actuator can output the optimal active control force of the rear axle;

step VI, energy recovery of the integrated self-powered active suspension actuator control system is carried out, and the specific process is as follows: when the integrated self-powered active suspension actuator control system is in MeIn the energy recovery mode, the first MOS switch module is disconnected, the second MOS switch module is connected, at the moment, wheels are excited by road surface unevenness, and when vertical acting force is generated, the lower end of the secondary iron core is connected with an axle through a second actuator mounting hole, the acting force is transmitted to the actuator, and at the moment, the S pole permanent magnet array and the N pole permanent magnet array which are fixed on the secondary iron core and the air spring piston fixedly connected with the secondary iron core vibrate vertically; the actuator mounting bracket is fixedly connected with the actuator outer shell, the actuator mounting bracket is fixedly connected with the frame through a first actuator mounting hole, the primary iron core is fixedly connected with the actuator outer shell, and a winding coil is wound in the primary iron core; the S pole permanent magnet array and the N pole permanent magnet array can generate relative displacement with the winding coil to cut the magnetic induction line to generate electric energy, and the speed of cutting the magnetic induction line by the winding coil is changed due to the fact that the excitation of the road surface unevenness is changed and the movement speed of the actuator is also changed, so that the energy feedback voltage is also changed, the energy feedback voltage needs to be processed by the rectifier filter module, and the energy feedback voltage is boosted to a voltage range which can be received by the energy storage module through the boosting module to complete vibration energy recovery.

Compared with the prior art, the invention has the following advantages:

1. the integrated self-powered active suspension actuator is suitable for heavy commercial vehicles with large full load and mass, and has the advantages of high structural integration degree, strong bearing capacity and small self rigidity, so the offset frequency is small.

2. The integrated self-powered active suspension actuator control system can perform vibration energy recovery in an energy recovery mode and vibration control in an active control mode, can control two important parameters of rigidity and damping of a suspension system, can keep energy balance, does not need to consume external energy, and achieves self-power supply of the system.

3. The invention provides a novel whole vehicle control algorithm, which is used for respectively controlling a front axle and a rear axle of a vehicle, and firstly, pre-aiming control is carried out on the front axle of the vehicle to obtain the road surface unevenness

Figure BDA0002232833140000161

Figure BDA0002232833140000162

Defining a front axis q1f、q2f、q3fThree weighting coefficients and performance index of LQR controller design in T time

Figure BDA0002232833140000163

Figure BDA0002232833140000164

Front axle gain matrix K obtained by LQR controllerfAccording to the formula Uf=-KfxfObtaining the optimal control force U of the front axlefBy passing a coil current I to the front axle actuatorfSo that the main power F of the front axle actuatoraf=UfAnd the front axle actuator is enabled to output the optimal active control force of the front axle. Inputting displacement x by collecting road surface unevenness at front wheel0fAnd is delayed

Figure BDA0002232833140000165

Then, the input displacement x of the road surface unevenness of the rear wheel is obtained0rAnd obtaining a state equation of the rear suspension system:defining a rear axis q1r、q2r、q3rThree weighting coefficients and performance index of LQR controller design in T time

Figure BDA0002232833140000167

Figure BDA0002232833140000168

Applying adaptive multi-forgetting factor recursive least squares to state variable gain matrix Qr、Rr、NrOptimizing, and setting a coefficient optimization matrix theta ═ QrRrNr]TSetting a covariance matrix Pi=106E, setting a data matrix phi and the like, and setting an expected optimal performance index JyAnd calculating an estimation error ekAnd performing parameter estimation:

Figure BDA0002232833140000169

updating the covariance matrix:

Figure BDA00022328331400001610

obtaining the optimal performance index J meeting the expectationyQ of (2)r、Rr、NrFurther obtain the optimal 3 weighting coefficients q of the rear axle1r、q2r、q3rThe LQR controller obtains a rear axle gain matrix KrAccording to the formula Ur=-KrxrCalculating the optimal control force U of the rear axlerBy passing a coil current I to the rear axle actuatorrSo that the main power F of the rear axle actuatorar=UrAnd the rear axle actuator outputs the optimal active control force of the rear axle.

4. The invention provides a mode switching rule that an actuator can realize self-powered under the condition of meeting energy balance, and a dynamic deflection threshold value f is setdyAnd duration threshold tyAnd setting the residual electric quantity of the energy storage module to be divided into a state of lacking electric quantity SOC1Normal state of charge SOC2State of sufficient charge SOC3And three levels are used for comprehensive judgment, and a mode switching rule is formulated.

5. The integrated self-powered active suspension actuator has high working stability and reliability, is not easy to break down, and does not need frequent maintenance and repair.

6. When the controller of the actuator breaks down, the system works in a passive mode, the function of the system is the same as that realized by a traditional common suspension, the damping force of the system is the electromagnetic resistance of the energy recovery mode of the linear motor, and the rigidity of the system is the combined rigidity of the air suspension and the spiral spring. The vibration damping system is prevented from being paralyzed due to system failure, and the vehicle running smoothness and the operation stability are prevented from being deteriorated.

Drawings

The technical solution of the present invention is further described in detail by the accompanying drawings and embodiments.

Fig. 1 is a cross-sectional view of an integrated self-powered active suspension actuator according to the present invention.

Fig. 2 is a plan view of an integrated self-powered active suspension actuator according to the present invention.

FIG. 3 is a block diagram of an integrated self-powered active suspension actuator control system according to the present invention.

Fig. 4 is a winding layout diagram of an integrated self-powered active suspension actuator winding coil according to the present invention.

Description of reference numerals:

Figure BDA0002232833140000171

Figure BDA0002232833140000181

Detailed Description

As shown in figure 1, an integrated self-powered active suspension actuator control system comprises an actuator 23, wherein the actuator 23 comprises an actuator mounting support 2, a first actuator mounting hole 1 is formed in the actuator mounting support 2, the actuator mounting support 2 is fixed at the top end of an actuator outer shell 10, a top fixed connection limiting block fixing seat 4 in the actuator outer shell 10 is arranged on the limiting block fixing seat 4, a rubber limiting block 3 used for limiting the maximum compression amount of the actuator 23 is arranged on the limiting block fixing seat 4, an air spring bag skin 5 arranged on the actuator outer shell 10 is arranged on the periphery of the rubber limiting block 3, the lower end of the air spring bag skin 5 is connected with an air spring piston 9 in an interference fit manner, an annular roll strip is arranged at the top in the actuator outer shell 10, the top end of the air spring bag skin 5 is arranged in a, the air spring piston 9 is internally provided with a limiting buffer seat support 8, the bottom of the limiting buffer seat support 8 is fixedly connected with the bottom of the air spring piston 9, the top of the limiting buffer seat support 8 is fixedly connected with a limiting buffer seat 6, the limiting buffer seat 6 is positioned outside the air spring piston 9 and inside an air spring capsule 5, the top of the limiting buffer seat support 8 is provided with a limiting buffer seat support inner cavity vent hole 7, the bottom of the air spring piston 9 is connected with a spiral spring 11, the lower part of the spiral spring 11 is connected with an actuator outer shell 10, the lower part of the air spring piston 9 is provided with a first permanent magnet fixing seat 12, the bottom of the air spring piston 9 is provided with an air spring additional cavity vent hole 13, the lower part of an actuator 23 is a part of a driving motor and comprises a secondary iron core 14, a primary iron core 17 and a, first permanent magnet fixing base 12 bottom is connected with secondary 14, S utmost point permanent magnet array 15 and N utmost point permanent magnet array 16 are installed to the secondary 14 outside, install elementary iron 17 between actuator shell body 10 and S utmost point permanent magnet array 15, the N utmost point permanent magnet array 16, the inside winding of elementary iron 17 has winding coil 18, second permanent magnet fixing base 19 is installed to secondary 14 lower part, first permanent magnet fixing base 12 all links to each other with S utmost point permanent magnet array 15, N utmost point permanent magnet array 16 with second permanent magnet fixing base 19, secondary 14 lower parts are equipped with second actuator mounting hole 20, actuator shell body 10 top is equipped with gas pressure sensor 21 and vent 22.

The actuator mounting bracket 2 is connected with the whole vehicle frame through the first actuator mounting hole 1 in a bolt connection mode; the upper plane of the actuator outer shell 10 is attached to the lower plane of the whole vehicle frame; the secondary iron core 14 is connected with the whole vehicle axle in a bolt connection mode through a second actuator mounting hole 20.

The air spring piston 9, the actuator outer shell 10, the first permanent magnet fixing seat 12, the secondary iron core 14 and the second permanent magnet fixing seat 19 are all made of magnetism isolating materials, and the primary iron core 17 is made of magnetic conducting materials.

The air spring bag skin 5 is internally provided with a main air chamber which is communicated with a middle air chamber formed by a limit buffer seat support 8 through an air vent hole 7 in an inner cavity of the limit buffer seat support, and the middle air chamber formed by the limit buffer seat support 8 is communicated with an additional air chamber formed by a secondary iron core 14 through an air vent hole 13 of an air spring additional cavity; a gap exists between the top of the air spring piston 9 and the limiting buffer seat 6, and a main air chamber inside the air spring bag skin 5 is communicated with an inner air chamber inside the air spring piston 9. The main air chamber in the air spring capsule 5, the middle air chamber formed by the limiting buffer seat support 8 and the additional air chamber formed by the secondary iron core 14 are communicated, and the acting area of the additional air chamber is smaller than that of the main air chamber, so that under a certain gas pressure state, the integrated self-powered active suspension actuator 23 generates larger displacement under the same acting force of the air spring. Therefore, the structure is beneficial to reducing the rigidity of a suspension system, reducing the offset frequency of the vehicle and improving the smoothness of the vehicle.

The air spring piston 9 bottom is connected with coil spring 11, coil spring 11 lower part is connected with actuator shell 10, therefore air spring piston 9 and coil spring 11 constitute serial-type compound spring system, compares and adopts single air spring piston 9 or coil spring 11, compound spring system rigidity is lower, therefore this structure is favorable to reducing suspension system rigidity, reduces the vehicle and deviates frequently, improves the vehicle ride comfort.

As shown in fig. 3, the winding coil 18 is a double-layer winding, i.e. there are two coils in each slot, each phase voltage is different by 120 °, and the distribution sequence of each phase voltage is: c-, A +; a-, B +; b-, C +; c-, A + … … where A, B, C represents three phases and +, -represents current direction.

As shown in fig. 2, the actuator control system further comprises an active control system and an energy recovery system, wherein the active control system comprises a vehicle speed sensor 29, a road surface irregularity detector 30, an unsprung mass displacement sensor 31, a sprung mass displacement sensor 32, an actuator speed sensor 33, a road surface irregularity displacement sensor 34, an on-board computer 35 and a control module 36; the energy recovery system comprises a driving module 24, a first MOS switching module 25, a rectifying and filtering module 26, a boosting module 27, a second MOS switching module 28, an OCV detection module 37, a control circuit power supply module 38 and an energy storage module 39. The vehicle speed sensor 29, the road surface irregularity detector 30, the unsprung mass displacement sensor 31, the sprung mass displacement sensor 32, the actuator speed sensor 33 and the road surface irregularity displacement sensor 34 are respectively connected with a control module 36, and the control module 36 is connected with a vehicle-mounted computer 35; the actuator 23 is sequentially connected with the driving module 24, the rectifying and filtering module 26, the boosting module 27, the second MOS switch module 28 and the energy storage module 39 to form a loop; the input end and the output end of the OCV detection module 37 are respectively connected with the output end of the energy storage module 39 and the input end of the control module 36, the output end of the control module 36 is connected with a first MOS switch module, and the output end of the first MOS switch module is connected with a circuit between the driving module 24 and the rectifying and filtering module 26; the control circuit power supply module 38 is connected in parallel between the input terminal of the OCV detection module 37 and a circuit connecting the first MOS switch module 25 and the control module 36.

The control module 36 is responsible for receiving sensor signals and performing circuit control; the vehicle-mounted computer 35 is responsible for calculating a control algorithm; the vehicle speed sensor 29 is responsible for periodically collecting vehicle speed data; the road surface irregularity detector 30 is mounted at the front of the vehicle and is responsible for detecting the road surface irregularity of the road ahead by a certain distance; the unsprung mass displacement sensor 31 is responsible for periodically collecting unsprung mass displacement; the sprung mass displacement sensor 32 is responsible for periodically collecting the sprung mass displacement; the actuator speed sensor 33 is responsible for periodically acquiring the movement speed of the actuator 23; the road surface irregularity displacement sensor 34 is responsible for periodically collecting the road surface irregularity information of the ground where the front axle is located; the gas pressure sensor 21 is responsible for periodically collecting the internal gas pressure of the actuator; the driving module 24 is responsible for controlling and driving the motor according to the control signal; the first MOS switch module 25 and the second MOS switch module 28 are responsible for switching between an energy recovery mode and an active control mode in the hardware circuit, and when the first MOS switch module 25 is connected and disconnected with the second MOS switch module 28, the actuator 23 is in the active control mode; when the first MOS switch module 25 is off and the second MOS switch module 28 is on, the actuator 23 is in the energy recovery mode; the rectifying and filtering module 26 is responsible for converting the dc voltage stored in the energy storage module 39 into the ac voltage required by the motor, or converting the ac voltage generated by the motor into the dc voltage that can be stored in the energy storage module 39; the boosting module 27 is responsible for boosting the electric energy generated by the actuator 23 in the energy recovery mode so as to store the electric energy by the energy storage module 39; the OCV detection module 37 is responsible for detecting the open circuit voltage of the energy storage module 39 and transmitting an open circuit voltage signal to the control module 36; the control circuit power supply module 38 is responsible for supplying the required electric energy to the control module 36; the energy storage module 39 is responsible for storing electrical energy.

A method of controlling an integrated self-powered active suspension actuator control system, comprising the steps of:

step I, detection data acquisition, which comprises the following specific processes: when the wheel is excited by road unevenness to generate a vertical acting force, the lower end of the secondary iron core 14 is connected with the axle through the second actuator mounting hole 20, the acting force is transmitted to the actuator 23, and the S pole permanent magnet array 15 and the N pole permanent magnet array 16 which are fixed on the secondary iron core 14 and the air spring piston 9 which is fixedly connected with the secondary iron core 14 vibrate vertically; the actuator mounting bracket 2 is fixedly connected with the actuator outer shell 10, the actuator mounting bracket 2 is fixedly connected with the frame through the first actuator mounting hole 1, the primary iron core 17 is fixedly connected with the actuator outer shell 10, the winding coil 18 is wound in the primary iron core 17, the S pole permanent magnet array 15 and the N pole permanent magnet array 16 generate relative displacement with the winding coil 18, the conductor cuts the magnetic induction line to generate electric energy, and the electric energy is in an energy recovery state; after the winding coil 18 is electrified, a traveling wave magnetic field is generated and moves along the axial direction of the actuator 23, and the running speeds of the secondary iron core 14 and the primary iron core 17 are differentThe slip rate is generated during sampling, so that electromagnetic force in the linear direction is generated, the electromagnetic force is the active control force of the actuator 23, and the active control force is adjusted in real time according to different running conditions of the vehicle, so that the actuator 23 can be kept in the optimal damping force state at any time; the expression of the slip ratio is

Figure BDA0002232833140000211

When the slip ratio is between 0 and 1, the actuator 23 is in a normal working state; when the slip ratio is greater than 1, the moving direction of the secondary iron core 14 is opposite to the traveling wave magnetic field direction, and the actuator 23 is in a braking state; when the slip ratio is less than 0, the movement speed of the secondary iron core 14 is greater than the speed of the traveling wave magnetic field, and the actuator 23 is in an energy recovery state;

under the above two conditions, the vehicle speed sensor 29 detects the vehicle speed in real time, the road unevenness detector 30 installed in the front of the vehicle detects the road unevenness of the road in front of the vehicle by a certain distance, and the road input displacement x at the front wheel of the vehicle after a certain time delay is obtained0f(ii) a The road surface irregularity displacement sensor 34 detects the road surface irregularity of the front wheels of the vehicle in real time to obtain the road surface input displacement x of the rear wheels of the vehicle after a certain delay0r(ii) a The unsprung mass displacement sensor 31 detects unsprung mass displacement in real time, the sprung mass displacement sensor 32 detects sprung mass displacement in real time, the actuator speed sensor 33 detects the motion speed of the suspension actuator in real time, and the control module 36 obtains vehicle speed V and unsprung mass displacement z of the front wheel of the vehicle, which are obtained by periodic sampling respectivelywfUnsprung mass displacement z of the rear wheel of the vehiclewrVehicle front wheel sprung mass displacement zbfVehicle rear wheel sprung mass displacement zbrAnd the moving speed u of the actuator 23 mounted on the front wheel of the vehicleafAnd the moving speed u of the actuator 23 mounted on the rear wheel of the vehiclearAnd transmitted to the vehicle computer 35;

step II, establishing a system state equation, wherein the specific process is as follows:

step 201, the vehicle-mounted computer 35 establishes a half vehicle model equation of motion, taking the left side of the vehicle as an example:

Figure BDA0002232833140000221

Figure BDA0002232833140000222

Figure BDA0002232833140000223

Figure BDA0002232833140000224

wherein the content of the first and second substances,

Figure BDA0002232833140000225

vertical acceleration of the sprung mass of the front wheel of the vehicle,

Figure BDA0002232833140000226

The vertical acceleration of the sprung mass of the rear wheel of the vehicle, the displacement of the sprung mass, and mbIs a sprung mass, IbIs the moment of inertia, theta, of the vehicle bodybIs the pitch angle, a is the distance from the front axle to the center of mass of the vehicle, b is the distance from the rear axle to the center of mass of the vehicle, ksfStiffness, k, of the front active suspension actuator 23srFor the stiffness, m, of the rear active suspension actuator 23wfFor unsprung mass, m, of the front wheels of the vehiclewrIs the unsprung mass of the rear wheel of the vehicle,

Figure BDA0002232833140000227

Is the vertical acceleration of the unsprung mass of the front wheel of the vehicle,

Figure BDA0002232833140000228

For unsprung mass vertical acceleration, k, of the rear wheels of the vehicletfAs the rigidity, k, of the front wheels of the vehicletrAs the stiffness, x, of the rear wheels of the vehicle0fInputting displacement, x, to the road surface of the front wheel0rInputting displacement for the road surface of the rear wheel;

when the pitch angle thetabSmaller, approximately have

Figure BDA0002232833140000231

The equation of motion of the 4-degree-of-freedom semi-vehicle model is expressed in a matrix form:

Figure BDA0002232833140000232

in the formula, the quality matrix [ M ] is:

[M]=diag(mb,Ib,mwf,mwr)

the damping matrix [ C ] is:

the stiffness matrix [ K ] is:

Figure BDA0002232833140000234

tire stiffness matrix [ K ]t]Comprises the following steps:

Figure BDA0002232833140000235

road surface input vector x0Comprises the following steps:

x0=[x0f,x0r]T

the system input vector x is:

x=[z θbzwfzwr]T

step 202, when the value of the suspension mass distribution coefficient epsilon is close to 1, 2-degree-of-freedom models can be respectively established for the front axis and the rear axis, and at the moment, a state variable is taken as follows:

Figure BDA0002232833140000241

wherein z iswIs the unsprung mass displacement of the vehicle, and zwIs taken as zwfOr zwr,zbIs the sprung mass displacement of the vehicle, and zbIs taken as zbfOr zbr(ii) a If the road speed input is white noise, then

Figure BDA0002232833140000242

The system state equation at this time is:

Figure BDA0002232833140000243

wherein u is the main power of the actuator, and the value of u is the main power u of the front wheel actuatorfOr rear wheel actuator main power ur(ii) a A is a constant matrix relating to state variables of the vehicle model and

Figure BDA0002232833140000245

b is a constant matrix relating to the variation of the speed of movement of the actuator 23 and

Figure BDA0002232833140000246

g is a constant matrix related to the input displacement variable of the road surface at the wheel and

Figure BDA0002232833140000247

Figure BDA0002232833140000248

csis damping of vehicle suspension and csIs taken as the front wheel actuator damping csfOr rear wheel actuator damping csr、mbIs the mass of the vehicle body, ksIs the suspension stiffness and ksIs taken as the front wheel actuator stiffness ksfOr rear wheel actuator stiffness ksr、ktIs tire stiffness and ktIs taken as the front wheel tire stiffness ktfOr rear wheel tire stiffness ktr、mwIs a vehicleUnsprung mass m of vehiclewIs taken as the unsprung mass m of the front wheelwfOr unsprung mass m of the rear wheelwr

Step III, the working mode switching rule of the integrated self-powered active suspension actuator control system is as follows:

step 301, defining an energy recovery mode of the integrated self-powered active suspension actuator control system as MeDefining the active control mode as Mc

Step 302, setting a dynamic deflection threshold fdyAnd duration threshold tyWherein dynamic deflection threshold fdyThe mean square root value of the dynamic deflection of the suspension under the condition of a good road surface and a low vehicle speed is smaller than that of a specific vehicle type; when the actual road surface dynamic deflection fd<fdyWhen the actual road surface dynamic deflection f is larger than the actual road surface dynamic deflection f, the vehicle is shown to run on a better road surfaced>fdyWhen, it indicates that the vehicle is traveling on a poor road surface; duration threshold tyThe influence on the working mode switching frequency of the integrated self-powered active suspension actuator control system is great, if the time threshold t is setyToo small, which may result in frequent switching of the operating mode of the integrated self-powered active suspension actuator control system, if the time threshold t is setyToo large, which may result in slow response of the integrated self-powered active suspension actuator control system to switching modes, time threshold tyThe specific numerical value of the road is determined according to the actual vehicle type and the common road condition;

step 303, the OCV detection module 37 collects the open-circuit voltage of the energy storage module 39, the control module 36 performs SOC estimation by using AH integration method, and the rated capacity of the energy storage module 39 is set to QeReal-time remaining capacity of QcInitial state capacity of Q0η is coulombic efficiency, then at tsReal-time remaining capacity at time QcThe state, using the integral operation, can be expressed as:

Figure BDA0002232833140000251

setting the energy storage in this caseThe remaining capacity of the module 39 is divided into 3 levels, respectively the state of lacking capacity SOC1Normal state of charge SOC2State of sufficient charge SOC3

Figure BDA0002232833140000252

Step 304, the main function of the integrated self-powered active suspension actuator control system is to ensure the ride comfort and handling stability of the vehicle, and the self-powered characteristic of the integrated self-powered active suspension actuator control system is combined, so that when the energy storage module 39 is in the SOC state2、SOC3In the state or when the running road surface condition is poor and the whole vehicle has large vibration, the vehicle is in McAn active control mode; energy storage module 39 is in SOC1The state or the running road surface condition is good, and when the whole vehicle vibrates less, the vehicle is in MeAn energy recovery mode; the mode switching rule is therefore:

step IV, performing pre-aiming optimal control on the unevenness of the front axle road surface of the integrated self-powered active suspension actuator control system, wherein the specific process is as follows:

step 401, the vehicle computer 35 defines the vehicle speed v at any moment iiLower, the pre-aiming distance LiWith the speed v of the vehicleiIs the aiming time t of the ith aiming pointiIs formulated as:

Figure BDA0002232833140000262

wherein, the value of i is a natural number from 1 to n, and n is the total sampling times of the road surface unevenness detector 30;

step 402, the vehicle-mounted computer 35 inputs the road surface input displacement x at the ith pre-aiming pointiInput displacement x from road surface at front wheel of vehicle0fThe Laplace functional relation is expressed as:

and e is-tisThe approximate calculation of (a) is expressed as:

Figure BDA0002232833140000264

then, the second-order Pade approximate calculation is taken, and when the preview time is not 0, the numerator and the denominator of the formula (F14) are simultaneously divided by

Figure BDA0002232833140000265

Obtaining:

Figure BDA0002232833140000271

wherein, ai0And ai1Are all coefficients and

the formula (F12) is then rewritten as:

[x0f(s)-zi(s)][ai0+ai1s+s2]=-2ai1sxi(s) (F16)

step 403, the vehicle computer 35 defines x0f(s)-xi(s)=yi(s),b1=-2ai1Performing an inverse Laplace transform on the formula (F16) yields:

Figure BDA0002232833140000273

redefining the state variables as:

ηi1=yi0xt=yt

Figure BDA0002232833140000274

β0=0

β1=b1-ai1β0=b1

and rewriting the state equation of the ith pre-aiming point as follows:

wherein the content of the first and second substances,

Figure BDA0002232833140000277

the state equation of the n preview points is obtained as follows:

Figure BDA0002232833140000278

wherein the content of the first and second substances,

Figure BDA0002232833140000279

Figure BDA0002232833140000281

xp=[x1x2… xn]T,

η=[η11η12η21η22… ηn1ηn2]T

step 404, the vehicle-mounted computer 35 defines the road surface input displacement variable at the wheel as:

x0f=Dηη+Eηxp(F20)

wherein D isηAnd EηAre all constant matrices and

Figure BDA0002232833140000282

Figure BDA0002232833140000283

Figure BDA0002232833140000284

λiis a weighting coefficient of the ith preview point and12+…+λnn is the total number of the preview points as 1;

step 405, substituting the formula (F20) into the formula (F8) to obtain a system state equation:

Figure BDA0002232833140000285

step 406, the vehicle mount computer 35 will q1fDefined as the front axle spring mass displacement weighting coefficient representing ride comfort, q2fDefining a front axle suspension dynamic stroke weighting coefficient which influences the vehicle body posture and is related to the structural design and the arrangement, and q3fThe front axle tire dynamic displacement weighting coefficient representing the steering stability is defined, and the performance index of the LQR controller design in the T moment is defined as:

Figure BDA0002232833140000286

wherein z iswfFor unsprung mass displacement of the front axle, zbfFor the displacement of the front axle spring-loaded mass,

Figure BDA0002232833140000287

vertical acceleration of the front axle spring load mass;

step 407, the vehicle mount computer 35 substitutes the system state equation (F8) into the equation (F22) to obtain:

Figure BDA0002232833140000288

wherein the state variable gain matrix

Figure BDA0002232833140000291

Control variable gain matrixVariable correlation matrix

Figure BDA0002232833140000293

Step 408, the vehicle computer 35 outputs an optimal control feedback gain matrix KfWherein

Figure BDA0002232833140000294

Figure BDA0002232833140000295

PcfSatisfy the Raccati equation

Figure BDA0002232833140000296

409, the vehicle computer 35 calculates the formula Uf=-KfxfObtaining the optimal control force U of the front axlef

Step 410, the linear motor electromagnetic thrust expression is

Faf=BδfLfIfF24

Wherein B isδfIs the average air gap magnetic density, L, at which the moving coil is locatedfIs the total length of a linear motor moving coil winding conductor, IfThe coil current is led into the front axle actuator;

by passing a coil current I to the front-axle actuatorfSo that Faf=UfThe front axle actuator can output the optimal active control force of the front axle;

step V, optimal pre-aiming optimization control of the rear axle base of the integrated self-powered active suspension actuator control system is carried out, and the specific process is as follows:

step 501, the road surface irregularity displacement sensor 34 collects the road surface irregularity impact on the front wheels of the vehicle, and the road surface irregularity input displacement x at the front wheels is input at the moment0fWill be delayed by a time ofInput displacement x equal to road surface unevenness of rear wheel0rWherein V is the running speed of the vehicle and L is the wheel base of the vehicle, thereby obtaining the road surface input x0rAnd obtaining a state equation of the rear suspension system:

Figure BDA0002232833140000301

step 502, let q1r、q2r、q3rThe vertical displacement of a vehicle body representing riding comfort, the suspension dynamic stroke which influences the posture of the vehicle body and is related to structural design and arrangement and the tire dynamic displacement representing operation stability are respectively controlled by 3 weighting coefficients, so that the performance index of the LQR controller design in the T moment is as follows:

Figure BDA0002232833140000302

the system state equation F25 is introduced into F26

Figure BDA0002232833140000303

Wherein the state variable gain matrix

Figure BDA0002232833140000304

Control variable gain matrix

Figure BDA0002232833140000305

Variable correlation matrix

Figure BDA0002232833140000306

Step 503, using adaptive multi-forgetting factor recursive least square method to Qr、Rr、NrOptimizing, and setting a coefficient optimization matrix theta ═ QrRrNr]TSetting a covariance matrix Pi=106E, wherein E is a third-order identity matrix, and a data matrix phi is set:

Figure BDA0002232833140000307

calculate time T

Figure BDA0002232833140000308

Setting an expected optimal performance index JyAnd calculating an estimation error ek

Figure BDA0002232833140000309

Calculating the value of the variable forgetting factor:

Figure BDA0002232833140000311

where l is the step size, calculate the gain K of a single parameteri,k

Where i is 1, 2 … n, where n is the number of identification parameters, calculating gain matrix Lk:

Figure BDA0002232833140000313

And (3) performing parameter estimation:

Figure BDA0002232833140000314

updating the covariance matrix:

step 504, repeat step 503 until the identification is finished to obtainMeets the expected optimal performance index JyQ of (2)r、Rr、NrFurther, 3 weighting coefficients q are obtained1r、q2r、q3r

Step 505, the vehicle computer 35 outputs an optimal control feedback gain matrix KrWherein

Figure BDA0002232833140000316

Figure BDA0002232833140000317

PcrSatisfy the Raccati equation

Figure BDA0002232833140000318

Step 506, the vehicle computer 35 according to the formula Ur=-KrxrCalculating the optimal control force U of the rear axler

Step 507, the linear motor electromagnetic thrust expression is

Far=BδrLrIrF35

Wherein B isδrIs the average air gap magnetic density, L, at which the moving coil is locatedrIs the total length of a linear motor moving coil winding conductor, IrThe coil current is led into the front axle actuator;

by passing a coil current I to the rear-axle actuatorrSo that Far=UrThe rear axle actuator can output the optimal active control force of the rear axle;

step VI, energy recovery of the integrated self-powered active suspension actuator control system is carried out, and the specific process is as follows: when the integrated self-powered active suspension actuator control system is in MeIn the energy recovery mode, the first MOS switch module 25 is turned off, and the second MOS switch module 28 is engaged, so that the wheel is excited by the road surface unevenness to generate a vertical force, the lower end of the secondary core 14 is connected to the axle through the second actuator mounting hole 20, and the force is transmitted to the actuator 23, and the S pole permanent magnet array 15, the N pole permanent magnet array 16, and the second MOS switch module 28 fixed to the secondary core 14 are connected to the axle through the second actuator mounting hole 20The air spring piston 9, which is fixedly connected with the secondary iron core 14, will vibrate vertically; the actuator mounting bracket 2 is fixedly connected with the actuator outer shell 10, the actuator mounting bracket 2 is fixedly connected with the frame through a first actuator mounting hole 1, the primary iron core 17 is fixedly connected with the actuator outer shell 10, and a winding coil 18 is wound in the primary iron core 17; the S pole permanent magnet array 15 and the N pole permanent magnet array 16 will generate relative displacement with the winding coil 18 to cut the magnetic induction line to generate electric energy, and as the road surface unevenness excitation is changed, the movement speed of the actuator 23 is also changed, so that the speed of the winding coil 18 cutting the magnetic induction line is also changed, and the energy feeding voltage is also changed, the energy feeding voltage needs to be processed by the rectifier filter module 26, and the energy feeding voltage is boosted to a voltage range which can be received by the energy storage module 39 by the boost module 27 to complete vibration energy recovery.

The above description is only a preferred embodiment of the present invention, and is not intended to limit the present invention, and all simple modifications, changes and equivalent structural changes made to the above embodiment according to the technical spirit of the present invention still fall within the protection scope of the technical solution of the present invention.

37页详细技术资料下载
上一篇:一种医用注射器针头装配设备
下一篇:自平衡移动平台以及自平衡移动方法

网友询问留言

已有0条留言

还没有人留言评论。精彩留言会获得点赞!

精彩留言,会给你点赞!