CN104494387B - A vehicle inertial suspension structure and its parameter determination method - Google Patents

A vehicle inertial suspension structure and its parameter determination method Download PDF

Info

Publication number
CN104494387B
CN104494387B CN201410637469.6A CN201410637469A CN104494387B CN 104494387 B CN104494387 B CN 104494387B CN 201410637469 A CN201410637469 A CN 201410637469A CN 104494387 B CN104494387 B CN 104494387B
Authority
CN
China
Prior art keywords
suspension
vehicle
population
inertial
individual
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201410637469.6A
Other languages
Chinese (zh)
Other versions
CN104494387A (en
Inventor
沈钰杰
陈龙
杨晓峰
张孝良
杨军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangsu University
Original Assignee
Jiangsu University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jiangsu University filed Critical Jiangsu University
Priority to CN201410637469.6A priority Critical patent/CN104494387B/en
Publication of CN104494387A publication Critical patent/CN104494387A/en
Application granted granted Critical
Publication of CN104494387B publication Critical patent/CN104494387B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Vehicle Body Suspensions (AREA)

Abstract

The invention discloses a kind of vehicle inertial matter suspension frame structure and parameter determination method thereof, belong to automotive suspension vibration isolation technique field.The vehicle inertial matter suspension proposed is made up of two-stage series connection, and the first order is " spring-damper " two element in parallel, and the second level is " used container spring " two element in parallel." used container spring " two element in parallel is utilized to produce the action effect of parallel resonance under given conditions, it is achieved the suppression to body quality vibration, the road roughness that can effectively buffer and decay impacts.Utilize quantum genetic algorithm that suspension parameter is optimized calculating, simulation analysis shows, using vehicle inertial matter suspension frame structure proposed by the invention and parameter determination method, suspension property index has clear improvement compared to the passive suspension of transmission, and vehicle ride performance, safety are the most effectively promoted.

Description

一种车辆惯质悬架结构及其参数确定方法A vehicle inertial suspension structure and its parameter determination method

技术领域technical field

本发明属于车辆悬架隔振领域,涉及一种车辆惯质悬架结构及其参数确定方法。The invention belongs to the field of vehicle suspension vibration isolation, and relates to a vehicle inertia suspension structure and a parameter determination method thereof.

背景技术Background technique

“惯质”的概念自提出以来,已经受到了控制工程及机械隔振领域学者们的广泛关注。迄今为止,应用惯容器的机械隔振网络隔振潜力已经在车辆悬架、建筑物隔振及高性能摩托车转向补偿等领域得到证实。Since the concept of "inertial mass" was proposed, it has received extensive attention from scholars in the fields of control engineering and mechanical vibration isolation. So far, the vibration isolation potential of mechanical isolation networks using inerters has been demonstrated in the fields of vehicle suspension, building vibration isolation, and high-performance motorcycle steering compensation.

惯容器与质量元件有着相同的物理属性,其两端点的特点更适合于工程应用,Michael Z.Q.Chen在《Influence of inerter on natural frequencies of vibrationsystems》一文中指出,惯容器对振动系统的固有频率有着重要影响,可使振动系统的共振频率提前,并对多自由度的振动模型进行了理论推导与证明。从目前的实际研究来看,对振动系统的主振系进行振动控制主要采用的方法是利用调谐质量阻尼器来实现,许多学者应用附加在车身质量或车轮质量上的动力吸振器对主振系的振动状态进行控制,然而由于质量块的笨重与安装条件的限制,此方法仍停留在理论研究阶段。The inerter has the same physical properties as the mass element, and the characteristics of its two ends are more suitable for engineering applications. Michael Z.Q.Chen pointed out in the article "Influence of inerter on natural frequencies of vibration systems" that the inerter has an important effect on the natural frequency of the vibration system. The influence can make the resonance frequency of the vibration system advance, and the multi-degree-of-freedom vibration model is theoretically deduced and proved. Judging from the current practical research, the main method of vibration control for the main vibration system of the vibration system is to use a tuned mass damper. However, due to the heavy weight of the mass and the limitation of installation conditions, this method is still in the theoretical research stage.

中国专利201010281992.1公开了一种两自由端点的动力吸振器,利用惯容器取代原有的质量元件,可以实现对主振动系统的振动控制。然而,通过分析其结构特点可知,由于惯容器与阻尼元件需要在弹性元件的支撑保护下才能发挥作用,两端点动力吸振器无法直接作为车辆悬架结构而付诸运用,因此需要对其进行改进设计以期获得较好的实用效果。Chinese patent 201010281992.1 discloses a dynamic vibration absorber with two free ends, which uses an inerter to replace the original mass element, and can realize the vibration control of the main vibration system. However, through the analysis of its structural characteristics, it can be seen that since the inerter and the damping element can only function under the support and protection of the elastic element, the dynamic shock absorber at both ends cannot be directly used as a vehicle suspension structure, so it needs to be improved Design in order to obtain better practical effect.

发明内容Contents of the invention

为解决两自由端点动力吸振器在车辆悬架系统中的应用难题,本发明提出一种车辆惯质悬架结构,并给出确定其结构参数的方法,使得应用此悬架的车辆行驶性能得以有效提升。In order to solve the application problem of the dynamic shock absorber with two free ends in the vehicle suspension system, this invention proposes a vehicle inertial suspension structure, and provides a method for determining its structural parameters, so that the driving performance of the vehicle using this suspension can be improved. Effective promotion.

为实现以上的发明目的,本发明采用的技术方案是:提出一种两级串联式车辆惯质悬架,第一级结构由第一级弹簧和惯容器并联组成,第二级结构由第二级弹簧和阻尼器并联组成。第一级弹簧与惯容器的上铰接点连接于车身,第二级弹簧与阻尼器的下铰接点连接于车轮。In order to achieve the above object of the invention, the technical solution adopted by the present invention is: to propose a two-stage tandem vehicle inertial suspension, the first-stage structure is composed of the first-stage spring and the inerter connected in parallel, and the second-stage structure is composed of the second The stage spring and damper are connected in parallel. The upper hinge point of the first-stage spring and the inerter is connected to the vehicle body, and the lower hinge point of the second-stage spring and the damper is connected to the wheel.

优选的,依据悬架运动动力学与量子遗传算法适应度计算法则,所述的车辆惯质悬架参数的确定方法,包括如下步骤:Preferably, according to the suspension motion dynamics and the quantum genetic algorithm fitness calculation rule, the method for determining the parameters of the inertial suspension of the vehicle includes the following steps:

(1)种群初始化:对量子遗传算法的种群Q(t)(t=0)进行初始化,设定种群大小、进化代数、量子比特编码、悬架参数优化范围;(1) Population initialization: initialize the population Q(t)(t=0) of the quantum genetic algorithm, and set the population size, evolution algebra, qubit coding, and suspension parameter optimization range;

(2)对初始种群Q(t)的每一个个体进行一次测量,计算各个个体的量子遗传算法适应度值,并记录最优个体及对应的适应度值;(2) Perform a measurement on each individual of the initial population Q(t), calculate the quantum genetic algorithm fitness value of each individual, and record the optimal individual and the corresponding fitness value;

(3)判断是否满足进化代数条件,若满足条件则算法结束,否则进行下一步;(4)种群代数t=t+1;(3) Judging whether the evolution algebra condition is met, if the condition is met, the algorithm ends, otherwise proceed to the next step; (4) Population algebra t=t+1;

(5)利用量子旋转门U(t)对各个个体实施进化更新,得到新的种群Q(t);(5) Use the quantum revolving door U(t) to perform evolutionary update on each individual to obtain a new population Q(t);

(6)对种群Q(t)的每个个体进行测量,并计算适应度值;(6) Measure each individual of the population Q(t), and calculate the fitness value;

(7)记录最优个体和对应的适应度值;(7) Record the optimal individual and the corresponding fitness value;

(8)判断是否满足进化代数条件,若满足进化代数则算法结束,否则返回步骤(4)。(8) Judging whether the evolutionary algebraic condition is satisfied, if the evolutionary algebraic condition is satisfied, the algorithm ends, otherwise return to step (4).

步骤(2)中,所述量子遗传算法适应度计算法则步骤如下:In step (2), the steps of the quantum genetic algorithm fitness calculation rule are as follows:

(1)根据牛顿第二定律,建立包含车身质量、车轮质量的两自由度四分之一悬架振动模型;(1) According to Newton's second law, a two-degree-of-freedom quarter suspension vibration model including body mass and wheel mass is established;

(2)建立包含传统被动悬架“弹簧-阻尼器”二元件并联的1/4车辆悬架振动模型,采用积分白噪声进行输入,通过时域仿真,获取该悬架在车速为20m/s的随机路面输入下车身加速度响应均方根值BApas,悬架动行程响应均方根值SWSpas,轮胎动载荷响应均方根值DTLpas(2) Establish a 1/4 vehicle suspension vibration model that includes the parallel connection of the traditional passive suspension "spring-damper" two components, use integral white noise as input, and obtain the suspension at a vehicle speed of 20m/s through time domain simulation The root mean square value of body acceleration response BA pas , the root mean square value of suspension dynamic stroke response SWS pas , and the root mean square value of tire dynamic load response DTL pas under random road surface input;

(3)建立包含本发明提出的车辆惯质悬架的1/4车辆悬架振动模型,同样采用积分白噪声进行输入,通过时域仿真,获取该悬架在车速为20m/s的随机路面输入下车身加速度响应均方根值BA,悬架动行程响应均方根值SWS,轮胎动载荷响应均方根值DTL。(3) Establish a 1/4 vehicle suspension vibration model comprising the vehicle inertial suspension proposed by the present invention, also adopt integral white noise to input, and obtain the suspension at a speed of 20m/s on a random road surface by time domain simulation Input the root mean square value BA of the acceleration response of the lower body, the root mean square value SWS of the dynamic stroke response of the suspension, and the root mean square value DTL of the tire dynamic load response.

步骤(2)中,所述量子遗传算法的适应度计算公式为:In step (2), the fitness calculation formula of the quantum genetic algorithm is:

ff == -- (( BABA BABA paspas ++ SWSSWS SWSSWS paspas ++ DTLDTL DTLDTL paspas )) ..

步骤(5)中,所述量子旋转门的更新进化规则为:In step (5), the update evolution rule of the quantum revolving door is:

αα ii ′′ ββ ii ′′ == Uu (( θθ ii )) αα ii ββ ii

步骤(5)中,所述量子旋转门的调整操作为:In step (5), the adjustment operation of the quantum revolving door is:

U ( θ i ) = cos ( θ i ) - sin ( θ i ) sin ( θ i ) cos ( θ i ) , θi由表1设定的旋转角选择策略而定。 u ( θ i ) = cos ( θ i ) - sin ( θ i ) sin ( θ i ) cos ( θ i ) , θi is determined by the rotation angle selection strategy set in Table 1.

表1旋转角选择策略图Table 1 Rotation angle selection strategy map

表1中,xi为当前染色体的第i位,besti为当前的最优染色体的第i位,f(x)为适应度函数,s(αii)为旋转角方向。In Table 1, xi is the i -th position of the current chromosome, best i is the i-th position of the current optimal chromosome, f(x) is the fitness function, and s(α i , β i ) is the rotation angle direction.

步骤(1)中,所述种群初始化为:设定种群大小为40,进化代数为200,第一级弹簧(1)、第二级弹簧(2)刚度优化范围为:[10000,30000]N·m-1,惯容器(3)惯质系数优化范围为:[100,1000]kg,阻尼器(4)阻尼系数优化范围:[1000,3000]N·s·m-1;各个体的量子比特编码均初始化为:使得染色体所表达的是其全部可能状态的等概率叠加。In step (1), the population initialization is as follows: set the population size to 40, the evolution algebra to 200, and the stiffness optimization range of the first-stage spring (1) and the second-stage spring (2) to be: [10000, 30000]N m -1 , the optimization range of the inertial coefficient of the inerter (3): [100,1000]kg, the optimization range of the damper (4) damping coefficient: [1000,3000]N·s·m -1 ; The qubit codes are all initialized as: What the chromosome expresses is the equal probability superposition of all possible states.

采用本发明的有益效果是:对照于电路系统中并联的“电容-电感”二元件,当并联的“惯容器-弹簧”二元件发生并联谐振时,可有效阻断振动传递。基于此提出了两级串联式车辆惯质悬架结构。应用此悬架结构可有效抑制和衰减路面不平度的冲击,提升悬架性能。所述悬架参数的确定方法是基于量子计算的遗传优化算法,将量子的态矢量表达引入遗传算法,有效解决了遗传算法易陷入局部极值的现象,且在较小的种群规模下就可获得满意的优化参数。The beneficial effect of adopting the present invention is: compared with the parallel connection of the "capacitor-inductance" two components in the circuit system, when the parallel connection of the "inverter-spring" two components resonates in parallel, the vibration transmission can be effectively blocked. Based on this, a two-stage tandem vehicle inertial suspension structure is proposed. The application of this suspension structure can effectively suppress and attenuate the impact of road unevenness and improve the suspension performance. The determination method of the suspension parameters is based on the genetic optimization algorithm of quantum computing, and the quantum state vector expression is introduced into the genetic algorithm, which effectively solves the phenomenon that the genetic algorithm is easy to fall into a local extremum, and can be obtained under a small population size. Obtain satisfactory optimization parameters.

附图说明Description of drawings

图1是一种车辆惯质悬架结构示意图;Fig. 1 is a schematic diagram of the structure of a vehicle inertial suspension;

图2是基于量子遗传算法的车辆惯质悬架参数确定方法流程图;Fig. 2 is a flow chart of a method for determining parameters of a vehicle inertial suspension based on a quantum genetic algorithm;

图3是悬架参数确定方法中建立的1/4悬架振动模型示意图;Fig. 3 is a schematic diagram of a 1/4 suspension vibration model established in the suspension parameter determination method;

图4是车身加速度频域增益对照图;Fig. 4 is a comparison diagram of body acceleration frequency domain gain;

图5是悬架动行程频域增益对照图;Fig. 5 is a comparison diagram of the frequency domain gain of the suspension stroke;

图6是轮胎动载荷频域增益对照图。Figure 6 is a comparison diagram of tire dynamic load frequency domain gain.

附图标记说明如下:1.第一级弹簧,2.第二级弹簧,3.惯容器,4.阻尼器,5.车身质量,6.车辆悬架,7.车轮质量,8.轮胎等效弹簧。Reference signs are explained as follows: 1. First-stage spring, 2. Second-stage spring, 3. Inerter, 4. Damper, 5. Body mass, 6. Vehicle suspension, 7. Wheel mass, 8. Tires, etc. effective spring.

具体实施方式detailed description

下面结合附图以及具体实施例对本发明作进一步的说明,但本发明的保护范围并不限于此。The present invention will be further described below in conjunction with the accompanying drawings and specific embodiments, but the protection scope of the present invention is not limited thereto.

在图1中,一种车辆惯质悬架由两级串联组成。第一级结构由第一级弹簧1和惯容器3并联组成;第二级由第二级弹簧2和阻尼器4并联组成。第一级弹簧1与惯容器3的上铰接点连接于车身,第二级弹簧2与阻尼器4的下铰接点连接于车轮,由此完成一种车辆惯质悬架结构的安装。In Figure 1, a vehicle inertial suspension consists of two stages connected in series. The first stage structure is composed of the first stage spring 1 and the inerter 3 in parallel; the second stage is composed of the second stage spring 2 and the damper 4 in parallel. The upper hinge point of the first-stage spring 1 and the inerter container 3 is connected to the vehicle body, and the lower hinge point of the second-stage spring 2 and the damper 4 is connected to the wheel, thereby completing the installation of a vehicle inertial suspension structure.

为确定车辆惯质悬架各元件的具体参数,本发明提供了一种第一级弹簧1,第二级弹簧2,惯容器3及阻尼器4的参数确定方法。In order to determine the specific parameters of each component of the inertial suspension of the vehicle, the present invention provides a method for determining the parameters of the first stage spring 1, the second stage spring 2, the inerter container 3 and the damper 4.

基于量子遗传算法的车辆惯质悬架参数确定方法具体流程图如图3所示,主要包括以下步骤:The specific flow chart of the method for determining the parameters of the vehicle inertial suspension based on the quantum genetic algorithm is shown in Figure 3, which mainly includes the following steps:

(1)种群初始化。对量子遗传算法的种群初始化,设定种群大小为40、进化代数为200,各个体的量子比特编码均初始化为:使得染色体所表达的是其全部可能状态的等概率叠加。(1) Population initialization. For the population initialization of the quantum genetic algorithm, the population size is set to 40, the evolution algebra is 200, and the qubit codes of each individual are initialized as: What the chromosome expresses is the equal probability superposition of all possible states.

此处悬架各参数的优化范围设定如下:Here, the optimization range of each parameter of the suspension is set as follows:

1000010000 NN ·· mm -- 11 ≤≤ kk 11 ,, kk 22 ≤≤ 3000030000 NN ·· mm -- 11 100100 kgkg ≤≤ bb ≤≤ 10001000 kgkg 10001000 NN ·· sthe s ·· mm -- 11 ≤≤ cc ≤≤ 30003000 NN ·· sthe s ·· mm -- 11

量子遗传算法适应度值计算法则设定如下:The calculation rules for the fitness value of the quantum genetic algorithm are set as follows:

根据牛顿第二定律,建立双质量悬架振动模型如图3所示,图中,车身质量5为320kg,车轮质量7为40kg,轮胎等效刚度8为190kN·m-1。路面输入Zg(t)采用如下路面输入模型:According to Newton's second law, the dual-mass suspension vibration model is established as shown in Figure 3. In the figure, the mass of the vehicle body 5 is 320kg, the mass of the wheel 7 is 40kg, and the equivalent stiffness of the tire 8 is 190kN·m -1 . The road surface input Z g (t) adopts the following road surface input model:

zz ·· gg (( tt )) == -- 22 ππ ff 00 zz gg (( tt )) ++ 22 ππ GG 00 uu ww (( tt ))

式中,下截止频率f0为0.01Hz,路面不平度系数G0为5×10-6m3cycle-1,w(t)为积分白噪声,车速为20m/s。In the formula, the lower cut-off frequency f 0 is 0.01Hz, the road surface roughness coefficient G 0 is 5×10 -6 m 3 cycle -1 , w(t) is integral white noise, and the vehicle speed is 20m/s.

首先建立包含传统被动悬架“弹簧-阻尼器”二元件并联的1/4车辆悬架振动模型,被动悬架弹簧刚度为22kN·m-1,阻尼系数为1000N·s·m-1。通过时域仿真,获取传统被动悬架在随机路面输入下车身加速度响应均方根值BApas,悬架动行程响应均方根值SWSpas,轮胎动载荷响应均方根值DTLpasFirstly, a 1/4 vehicle suspension vibration model including the parallel connection of two components of the traditional passive suspension "spring-damper" is established. The spring stiffness of the passive suspension is 22kN·m -1 and the damping coefficient is 1000N·s·m -1 . Through time-domain simulation, the root mean square value BA pas of body acceleration response, the root mean square value SWS pas of suspension dynamic travel response and the root mean square value DTL pas of tire dynamic load response are obtained under random road input of traditional passive suspension.

建立包含本发明提出的车辆惯质悬架的1/4车辆悬架振动模型,在同等路面输入下获取车身加速度响应均方根值BA,悬架动行程响应均方根值SWS,轮胎动载荷响应均方根值DTL。Establish a 1/4 vehicle suspension vibration model including the vehicle inertial suspension proposed by the present invention, obtain the vehicle body acceleration response root mean square value BA, suspension dynamic stroke response root mean square value SWS, and tire dynamic load under the same road surface input Response RMS DTL.

量子遗传算法的适应度计算公式为:The fitness calculation formula of quantum genetic algorithm is:

ff == -- (( BABA BABA paspas ++ SWSSWS SWSSWS paspas ++ DTLDTL DTLDTL paspas )) ..

(2)对初始种群Q(t)的每一个个体进行一次测量,计算各个个体的量子遗传算法适应度值,并记录最优个体及对应的适应度值。(2) Measure each individual of the initial population Q(t), calculate the fitness value of the quantum genetic algorithm of each individual, and record the optimal individual and the corresponding fitness value.

(3)判断是否满足进化代数条件,若满足条件则算法结束,否则进行下一步。(3) Judging whether the condition of evolution algebra is satisfied, if the condition is satisfied, the algorithm ends, otherwise, proceed to the next step.

(4)种群代数t=t+1。(4) Population algebra t=t+1.

(5)利用量子旋转门U(t)对各个个体实施进化更新,得到新的种群Q(t)。(5) Use the quantum revolving door U(t) to perform evolutionary update on each individual to obtain a new population Q(t).

量子旋转门的调整操作为: U ( θ i ) = cos ( θ i ) - sin ( θ i ) sin ( θ i ) cos ( θ i ) The adjustment operation of the quantum revolving door is: u ( θ i ) = cos ( θ i ) - sin ( θ i ) sin ( θ i ) cos ( θ i )

其更新过程如下: α i ′ β i ′ = U ( θ i ) α i β i Its update process is as follows: α i ′ β i ′ = u ( θ i ) α i β i

其中,(αii)T与(αi′,βi′)T代表染色体第i个量子比特旋转门更新前后的概率幅。θi为旋转角,由表1设定的旋转角策略而定。Among them, (α i , β i ) T and (α i ′, β i ′) T represent the probability amplitudes before and after the revolving door update of the ith qubit of the chromosome. θi is the rotation angle, which is determined by the rotation angle strategy set in Table 1.

表1中,xi为当前染色体的第i位,besti为当前的最优染色体的第i位,f(x)为适应度函数,s(αii)为旋转角方向。、In Table 1, xi is the i -th position of the current chromosome, best i is the i-th position of the current optimal chromosome, f(x) is the fitness function, and s(α i , β i ) is the rotation angle direction. ,

(6)对种群Q(t)的每个个体进行测量,并计算适应度值。(6) Measure each individual of the population Q(t), and calculate the fitness value.

(7)记录最优个体和对应的适应度值。(7) Record the optimal individual and the corresponding fitness value.

(8)判断是否满足进化代数条件,若满足进化代数则算法结束,否则返回步骤4。(8) Judging whether the evolutionary algebraic condition is satisfied, if the evolutionary algebraic condition is satisfied, the algorithm ends, otherwise return to step 4.

在Matlab环境下仿真优化得到车辆惯质悬架的元件参数为:The component parameters of the inertial suspension of the vehicle obtained by simulation optimization under the Matlab environment are:

第一级弹簧1的刚度为:30000kN·m-1 The stiffness of the first stage spring 1 is: 30000kN·m -1

第二级弹簧2的刚度为:10000kN·m-1 The stiffness of the second stage spring 2 is: 10000kN·m -1

惯容器3的惯质系数为:312kgThe inertial coefficient of the inertial vessel 3 is: 312kg

阻尼器4的阻尼系数为:2080N·s·m-1 The damping coefficient of the damper 4 is: 2080N·s·m -1

将通过本发明提出的车辆惯质悬架参数确定方法获取的元件参数代入1/4悬架振动模型中,得到其车身加速度、悬架动行程、轮胎动载荷频域响应分析图如图4-6所示,从图中可以看出,相较于传统被动悬架,车辆惯质悬架可有效缓冲和衰减路面不平度的冲击,尤其是在车身共振偏频处,来自路面不平度的振动冲击得到了明显的衰减,频域分析表明,本发明提出的车辆惯质悬架体现出较好的性能优势。Substituting the component parameters obtained by the vehicle inertial suspension parameter determination method proposed by the present invention into the 1/4 suspension vibration model, the frequency domain response analysis diagram of its vehicle body acceleration, suspension dynamic stroke, and tire dynamic load is obtained as shown in Figure 4- As shown in Figure 6, it can be seen from the figure that compared with the traditional passive suspension, the inertial suspension of the vehicle can effectively buffer and attenuate the impact of the unevenness of the road surface, especially at the resonance bias frequency of the vehicle body, the vibration from the unevenness of the road surface The impact is obviously attenuated, and the frequency domain analysis shows that the vehicle inertial suspension proposed by the invention shows better performance advantages.

所述实施例为本发明的优选的实施方式,但本发明并不限于上述实施方式,在不背离本发明的实质内容的情况下,本领域技术人员能够做出的任何显而易见的改进、替换或变型均属于本发明的保护范围。The described embodiment is a preferred implementation of the present invention, but the present invention is not limited to the above-mentioned implementation, without departing from the essence of the present invention, any obvious improvement, replacement or modification that those skilled in the art can make Modifications all belong to the protection scope of the present invention.

Claims (5)

1.一种车辆惯质悬架结构参数的确定方法,所述车辆惯质悬架结构,包括第一级结构和第二级结构,所述第一级结构由第一级弹簧(1)和惯容器(3)并联组成,所述第二级结构由第二级弹簧(2)和阻尼器(4)并联组成;所述第一级结构和第二级结构为串联连接;所述第一级弹簧(1)与惯容器(3)的上铰接点连接于车身,所述第二级弹簧(2)与阻尼器(4)的下铰接点连接于车轮,其特征在于,所述车辆惯质悬架结构参数的确定方法包括如下步骤:1. A method for determining a vehicle inertial suspension structure parameter, the vehicle inertial suspension structure includes a first-level structure and a second-level structure, and the first-level structure is composed of a first-level spring (1) and a second-level structure Inerters (3) are connected in parallel, and the second-level structure is composed of second-level springs (2) and dampers (4) in parallel; the first-level structure and the second-level structure are connected in series; the first The upper hinge point of the first-stage spring (1) and the inertia container (3) is connected to the vehicle body, and the lower hinge point of the second-stage spring (2) and the damper (4) is connected to the wheel. It is characterized in that the vehicle inertia The method for determining the structural parameters of the mass suspension includes the following steps: (1)种群初始化:对量子遗传算法的种群Q(t)(t=0)进行初始化,设定种群大小、进化代数、量子比特编码、悬架参数优化范围;(1) Population initialization: initialize the population Q(t)(t=0) of the quantum genetic algorithm, and set the population size, evolution algebra, qubit coding, and suspension parameter optimization range; (2)对初始种群Q(t)的每一个个体进行一次测量,计算各个个体的量子遗传算法适应度值,并记录最优个体及对应的适应度值;(2) Perform a measurement on each individual of the initial population Q(t), calculate the quantum genetic algorithm fitness value of each individual, and record the optimal individual and the corresponding fitness value; (3)判断是否满足进化代数条件,若满足条件则算法结束,否则进行下一步;(3) Judging whether the evolutionary algebra condition is met, if the condition is met, the algorithm ends, otherwise proceed to the next step; (4)种群代数t=t+1;(4) population algebra t=t+1; (5)利用量子旋转门U(t)对各个个体实施进化更新,得到新的种群Q(t);(5) Use the quantum revolving door U(t) to perform evolutionary update on each individual to obtain a new population Q(t); (6)对种群Q(t)的每个个体进行测量,并计算适应度值;(6) Measure each individual of the population Q(t), and calculate the fitness value; (7)记录最优个体和对应的适应度值;(7) Record the optimal individual and the corresponding fitness value; (8)判断是否满足进化代数条件,若满足进化代数则算法结束,否则返回步骤(4)。(8) Judging whether the evolutionary algebraic condition is satisfied, if the evolutionary algebraic condition is satisfied, the algorithm ends, otherwise return to step (4). 2.如权利要求1所述的一种车辆惯质悬架结构参数的确定方法,其特征在于,步骤(2)中,所述量子遗传算法适应度计算法则步骤如下:2. the determining method of a kind of vehicle inertial suspension structure parameter as claimed in claim 1, is characterized in that, in step (2), described quantum genetic algorithm fitness calculation rule step is as follows: (1)根据牛顿第二定律,建立包含车身质量、车轮质量的两自由度四分之一悬架振动模型;(1) According to Newton's second law, a two-degree-of-freedom quarter suspension vibration model including body mass and wheel mass is established; (2)建立包含传统被动悬架“弹簧-阻尼器”二元件并联的1/4车辆悬架振动模型,采用积分白噪声进行输入,通过时域仿真,获取该悬架在车速为20m/s的随机路面输入下车身加速度响应均方根值BApas,悬架动行程响应均方根值SWSpas,轮胎动载荷响应均方根值DTLpas(2) Establish a 1/4 vehicle suspension vibration model that includes the parallel connection of the traditional passive suspension "spring-damper" two components, use integral white noise as input, and obtain the suspension at a vehicle speed of 20m/s through time domain simulation The root mean square value of body acceleration response BA pas , the root mean square value of suspension dynamic stroke response SWS pas , and the root mean square value of tire dynamic load response DTL pas under random road surface input; (3)建立包含车辆惯质悬架的1/4车辆悬架振动模型,同样采用积分白噪声进行输入,通过时域仿真,获取该悬架在车速为20m/s的随机路面输入下车身加速度响应均方根值BA,悬架动行程响应均方根值SWS,轮胎动载荷响应均方根值DTL。(3) Establish a 1/4 vehicle suspension vibration model including the vehicle inertial suspension, and also use integral white noise for input, and through time domain simulation, obtain the vehicle body acceleration of the suspension on a random road with a vehicle speed of 20m/s Response root mean square value BA, suspension dynamic stroke response root mean square value SWS, tire dynamic load response root mean square value DTL. 3.如权利要求2所述的一种车辆惯质悬架结构参数的确定方法,其特征在于,步骤(2)中,所述量子遗传算法的适应度计算公式为:3. the determining method of a kind of vehicle inertial suspension structure parameter as claimed in claim 2, is characterized in that, in step (2), the fitness computing formula of described quantum genetic algorithm is: ff == -- (( BB AA BABA pp aa sthe s ++ SS WW SS SWSSWS pp aa sthe s ++ DD. TT LL DTLDTL pp aa sthe s )) .. 4.如权利要求1所述的一种车辆惯质悬架结构参数的确定方法,其特征在于,步骤(5)中,4. the determining method of a kind of vehicle inertial suspension structure parameter as claimed in claim 1, is characterized in that, in step (5), 所述量子旋转门的更新进化规则为:The update evolution rule of the quantum revolving door is: αα ii ′′ ββ ii ′′ == Uu (( θθ ii )) αα ii ββ ii 所述量子旋转门的调整操作为:The adjustment operation of the quantum revolving door is: θi由表1设定的旋转角选择策略而定; θi is determined by the rotation angle selection strategy set in Table 1; 表1旋转角选择策略图Table 1 Rotation angle selection strategy map 表1中,xi为当前染色体的第i位,besti为当前的最优染色体的第i位,f(x)为适应度函数,s(αii)为旋转角方向;In Table 1, x i is the i-th position of the current chromosome, best i is the i-th position of the current optimal chromosome, f(x) is the fitness function, and s(α i , β i ) is the rotation angle direction; ii)T与(αi′,βi′)T分别代表染色体第i个量子比特旋转门更新前后的概率幅。i , β i ) T and (α i ′, β i ′) T represent the probability amplitudes before and after the revolving door update of the i-th qubit of the chromosome, respectively. 5.如权利要求1所述的一种车辆惯质悬架结构参数的确定方法,其特征在于,步骤(1)中,所述种群初始化为:设定种群大小为40,进化代数为200,第一级弹簧(1)、第二级弹簧(2)刚度优化范围为:[10000,30000]N·m-1,惯容器(3)惯质系数优化范围为:[100,1000]kg,阻尼器(4)阻尼系数优化范围:[1000,3000]N·s·m-1;各个体的量子比特编码均初始化为:使得染色体所表达的是其全部可能状态的等概率叠加。5. the determining method of a kind of vehicle inertial suspension structure parameter as claimed in claim 1, is characterized in that, in step (1), described population initialization is: setting population size is 40, and evolutionary number is 200, The optimization range of the stiffness of the first stage spring (1) and the second stage spring (2) is: [10000,30000]N m -1 , the optimization range of the inertial coefficient of the inerter (3) is: [100,1000]kg, Damper (4) damping coefficient optimization range: [1000,3000]N s m -1 ; the qubit codes of each individual are initialized as: What the chromosome expresses is the equal probability superposition of all possible states.
CN201410637469.6A 2014-11-12 2014-11-12 A vehicle inertial suspension structure and its parameter determination method Active CN104494387B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410637469.6A CN104494387B (en) 2014-11-12 2014-11-12 A vehicle inertial suspension structure and its parameter determination method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410637469.6A CN104494387B (en) 2014-11-12 2014-11-12 A vehicle inertial suspension structure and its parameter determination method

Publications (2)

Publication Number Publication Date
CN104494387A CN104494387A (en) 2015-04-08
CN104494387B true CN104494387B (en) 2016-09-14

Family

ID=52935876

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410637469.6A Active CN104494387B (en) 2014-11-12 2014-11-12 A vehicle inertial suspension structure and its parameter determination method

Country Status (1)

Country Link
CN (1) CN104494387B (en)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106051022B (en) 2016-05-09 2018-06-26 江苏大学 A kind of fluid power recalls used case and its application
CN106407507B (en) * 2016-08-26 2019-08-27 吉林大学 A method for analyzing ride comfort of in-wheel electric drive vehicles and application of motor selection
CN107563001B (en) * 2017-08-01 2020-02-21 江苏大学 A positive-real synthesis and passive implementation method of a vehicle ISD suspension electromechanical network system
CN108240286B (en) * 2018-01-29 2019-10-15 河海大学 Parameter optimization method for passive structure control device of floating offshore wind turbine based on inertia
CN108896326B (en) * 2018-06-15 2020-05-01 江苏大学 Vehicle ISD suspension parameter selection and test condition design method
CN108944323B (en) * 2018-07-16 2020-05-05 东北大学 A two-stage vibration damping automobile suspension structure and locking control method
CN109334379B (en) * 2018-09-19 2021-11-23 江苏大学 Method for identifying nonlinear inertia force of hydro-electric inertial container for heavy-duty vehicle
CN109773727A (en) * 2018-12-28 2019-05-21 南京理工大学 A Six-DOF Platform with Inertial Capacity
CN109820697A (en) * 2019-02-26 2019-05-31 南京理工大学 Multifunctional assisted walking exoskeleton with inertia
CN111086363A (en) * 2020-01-20 2020-05-01 中国北方车辆研究所 Two-stage serial-type suspension structure with double inertial containers
CN113496299A (en) * 2020-03-19 2021-10-12 上海达牛信息技术有限公司 Vehicle path planning method based on quantum genetic algorithm
CN112855824A (en) * 2021-03-11 2021-05-28 中国北方车辆研究所 Inertial volume type vibration isolation device aiming at simple harmonic displacement excitation and design method
CN113606278B (en) * 2021-08-15 2023-03-14 深圳大学 Multi freedom vibration isolation system
CN114312513B (en) * 2021-12-06 2024-04-09 江苏大学 Mechanical memory element and application thereof

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011089373A1 (en) * 2010-01-25 2011-07-28 Lotus Renault Gp Ltd Fluid inerter
CN101994776B (en) * 2010-09-14 2013-07-17 江苏大学 Two-free endpoint dynamic vibration absorber
CN102937159A (en) * 2012-10-31 2013-02-20 江苏大学 Integrated hydro-pneumatic spring device
CN203293845U (en) * 2013-06-14 2013-11-20 江苏大学 Vehicle passive suspension structure using inerter

Also Published As

Publication number Publication date
CN104494387A (en) 2015-04-08

Similar Documents

Publication Publication Date Title
CN104494387B (en) A vehicle inertial suspension structure and its parameter determination method
CN103593506B (en) A Two-Stage Series ISD Suspension Parameter Optimization Method
CN103072440B (en) Control method for automotive active suspension system
CN101327722B (en) Vehicle Suspension Using Inertial Accumulator
CN105644288B (en) A kind of vehicle ISD suspension parameter matching process
CN105172510B (en) A kind of Vehicle Suspension Control System and control method based on virtual used appearance
CN106250592A (en) A kind of vehicle inertial matter optimization of suspension parameters method for designing
CN105159094A (en) Design method of optimal control force of LQG controller of automobile active suspension bracket
CN111086363A (en) Two-stage serial-type suspension structure with double inertial containers
CN113103836B (en) Asymmetric reciprocating damping-based vehicle ISD suspension structure and optimal design method
CN103076148A (en) Drop test dual-four-degree-of-freedom half vehicle model-based vehicle parameter identification method
Sun et al. Research of simulation on the effect of suspension damping on vehicle ride
Nie et al. Passive skyhook suspension reduction for improvement of ride comfort in an off-road vehicle
CN203293845U (en) Vehicle passive suspension structure using inerter
CN103303085A (en) Vehicle passive suspension structure employing inertial container
Shi et al. Influence of equipment excitation on flexible carbody vibration of EMU
CN104691264A (en) Inerter suspension of vehicle
CN203974458U (en) A kind of passive suspension frame structure of vehicle of being used to container that contains
CN106951588B (en) A kind of vehicle ISD suspension research method based on coupling effect
Stepanov et al. Prospects of Using Hyperelastic Elastomer Materials in Vehicle Suspension Systems
Yang et al. Sliding mode control of chaos vibrations in 7-DOF nonlinear active vehicle suspension
CN105302944A (en) Generalized dissipated energy based calculation method for optimal damping ratio of vehicle passive suspension
DUMITRIU Car vertical dynamics 3D simulator using a 7 DOF model
Zhou et al. Optimization design and simulation of an air suspension with auxiliary chamber and throttle valve
CN103522862B (en) A kind of method determining semi-active suspension equivalent damping maxim

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant