CN112948979B - 汽车电控阻尼悬架***冲击工况最优宾汉数控制方法 - Google Patents
汽车电控阻尼悬架***冲击工况最优宾汉数控制方法 Download PDFInfo
- Publication number
- CN112948979B CN112948979B CN202110339869.9A CN202110339869A CN112948979B CN 112948979 B CN112948979 B CN 112948979B CN 202110339869 A CN202110339869 A CN 202110339869A CN 112948979 B CN112948979 B CN 112948979B
- Authority
- CN
- China
- Prior art keywords
- electric control
- formula
- impact
- damping
- control damping
- 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
Links
- 238000013016 damping Methods 0.000 title claims abstract description 160
- 238000000034 method Methods 0.000 title claims abstract description 109
- 239000000725 suspension Substances 0.000 title claims abstract description 42
- 230000003139 buffering effect Effects 0.000 claims abstract description 87
- 238000007906 compression Methods 0.000 claims abstract description 54
- 238000006073 displacement reaction Methods 0.000 claims abstract description 24
- 230000006835 compression Effects 0.000 claims abstract description 17
- 230000005484 gravity Effects 0.000 claims description 21
- 230000035939 shock Effects 0.000 claims description 12
- 230000001133 acceleration Effects 0.000 claims description 9
- 238000013178 mathematical model Methods 0.000 claims description 7
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 3
- 241000695274 Processa Species 0.000 claims 1
- 230000006378 damage Effects 0.000 abstract description 6
- 230000005284 excitation Effects 0.000 abstract description 5
- 238000011217 control strategy Methods 0.000 description 5
- 238000010521 absorption reaction Methods 0.000 description 2
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 239000006096 absorbing agent Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Vehicle Body Suspensions (AREA)
- Vibration Prevention Devices (AREA)
Abstract
本发明公开了一种汽车电控阻尼悬架***冲击工况最优宾汉数控制方法,其依附硬件包括车辆电控阻尼悬架***以及位移传感器与力传感器,步骤:1)根据车辆悬架在拉伸与压缩状态下的力学特性建立冲击工况下汽车悬架***动力学模型;2)利用动力学模型与悬架控制初始条件计算电控阻尼悬架在冲击激励下的运动状态,并计算压缩与回弹过程的最优宾汉数与相应期望力;3)基于多算子模型,建立汽车电控阻尼执行器非线性力学模型;4)利用电控阻尼执行器非线性力学模型跟踪期望力,计算与最优宾汉数匹配的最优输出功率。本发明能够快速、精确地实现冲击缓冲的自适应控制,从而能有效降低冲击激励对车身/车架结构和驾乘人员的结构性损伤。
Description
技术领域
本发明涉及汽车垂向动力学建模及控制领域,具体地说是一种汽车电控阻尼悬架冲击缓冲动力学建模与控制方法。
背景技术
实现冲击激励的自适应控制,降低冲击载荷对工程设备损坏和操作人员带来的危害,是许多工程应用中无法忽视的难题,诸如航天器的着陆装置,座椅碰撞保护装置,火炮枪械后座缓冲装置等。因此,最大化地降低冲击载荷造成的危害,不仅是冲击控制领域一个重要的研究课题,同时也具有巨大且广泛的应用价值。由于冲击激励的突发性与瞬时性,传统结构往往难以满足冲击缓冲控制的快速响应和自适应的要求。而磁流变智能材料因毫秒级的响应时间以及力学特性连续可变,在冲击缓冲控制领域具有非常好的应用前景。目前电控阻尼悬架***已广泛应用于汽车振动控制领域,但关于冲击缓冲的控制策略还不是很多。
ChulheeHan和Byung-HyukKang以电控飞机起落架为研究对象,在天棚控制器的基础上设计了一种混合控制器用于降低飞机着陆时的冲击载荷,提高了飞机的着陆效率,降低了飞机着陆时对成员造成的不舒适性。但该混合控制方法没有充分利用控制减振器的行程,不能最大程度降低冲击载荷。
恒定力冲击缓冲控制策略最大化地用了减振器行程,且在冲击缓冲过程中阻尼力波动很小。但恒定力控制策略需要实时根据上一时刻的加速度状态调整当前状态阻尼力输出,导致控制器的需要在每一个控制步长处理位移、速度、加速度信号,并调整输出功率。恒定力控制方法参数较多,对于控制器处理数据的能力具有一定的要求。
发明内容
本发明为了解决上述技术所存在的不足提出了一种汽车电控阻尼悬架***冲击工况最优宾汉数控制方法,以期能够快速、精确地实现冲击缓冲的自适应控制,从而能有效降低冲击激励对车身/车架结构和驾乘人员造成的结构性损伤和伤害。
本发明为了达到上述的发明目的,采用以下的技术方案:
本发明一种汽车电控阻尼悬架***在冲击工况下最优宾汉数控制方法的特点是按如下步骤进行:
步骤1:利用式(1)建立冲击缓冲压缩过程中的电控阻尼悬架动力学模型:
mz″(t)+cz′(t)-fyc+mg+kz(t)=0 (1)
式(1)中,m为簧载质量;z″(t)、z′(t)与z(t)是电控阻尼执行器在缓冲过程中t时刻的加速度、速度与位移;t表示时刻;c表示电控阻尼执行器的阻尼系数;fyc为冲击缓冲压缩过程中电控阻尼执行器的可控阻尼力;g为重力加速度;k为弹簧刚度系数并有:
Fyc=fyc-mg (3)
步骤2:利用冲击缓冲压缩过程的初始状态计算式(2),得到如式(4)-式(6)所示的冲击缓冲压缩过程中电控阻尼执行器的运动状态方程:
式(4)-式(6)中,wd为汽车电控阻尼悬架***中有阻尼固有频率;P1和P2为待求解的参数;并通过将冲击缓冲压缩过程的初始状态带入式(4)和式(5)中得到:
式(7)中,v0表示冲击缓冲压缩过程的初始状态中的初始速度,且初始移位为零;
步骤3:利用式(8)建立冲击缓冲压缩过程的期望终止状态:
zdesiring(tsc)=-S,z′desiring(tsc)=0 (8)
式(8)中,zdesiring(tsc)表示冲击缓冲压缩过程的终止状态的期望位移,z′desiring(tsc)表示冲击缓冲压缩过程的终止状态的期望速度,S表示冲击缓冲压缩过程的终止时,电控阻尼执行器的期望位置;tsc表示冲击缓冲压缩过程的时间,并有:
步骤4:利用式(10)建立抵消重力后的最优宾汉数Bio c+:
设定牛顿迭代法的迭代终止条件为从而利用牛顿迭代法对式(10)进行求解并得到最终迭代出的抵消重力后的最优宾汉数Bio *c+;其中,表示第n次迭代时抵消重力后的最优宾汉数,表示第n+1次迭代时抵消重力后的最优宾汉数;δ表示终止阈值;进而利用式(11)得到冲击缓冲压缩过程中的最优宾汉数
步骤5:对汽车电控阻尼悬架***中的电控阻尼执行器进行数学建模,并通过拟合实验数据进行参数辨识,从而得到如式(12)-式(14)所示的电控阻尼执行器的非线性力学特性的数学模型:
c=c1I2+c2I+c3 (13)
α=α1I2+α2I+α3 (14)
式(12)-式(14)中,Fmodel(t)为电控阻尼执行器数学模型建立的力学特性;α表示为电控阻尼执行器的磁滞系数;c1、c2、c3表示三个阻尼参数,并用于拟合阻尼系数c;α1、α2、α3表示三个磁滞参数,并用于拟合磁滞系数α;为磁滞输出系数;f为初始状态下电控阻尼执行器的偏置力;I为电控阻尼执行器的输出功率;
步骤6:利用式(15)计算冲击缓冲压缩过程中电控阻尼执行器的期望阻尼力Fdesiring_c:
Fdesiring_c=Bio ccv0 (15)
利用式(16)得到t时刻的力跟踪的误差Ferror(t):
Ferror(t)=Fmodel(t)-Fdesiring_c=AIc 2+BIc+C (16)
步骤7:利用式(12)对冲击缓冲压缩过程期望阻尼力Fdesiring_c进行力跟踪,并使得误差Ferror(t)最小,从而得到冲击缓冲压缩过程中电控阻尼执行器的期望力对应的最优输出功率Ic_op;
步骤8:若|z′(t)-z′(t-1)|<ε,则判定冲击缓冲压缩过程完成,并切换为冲击缓冲回弹状态后,执行步骤9,否则返回步骤7执行,其中,z′(t)表示电控阻尼执行器在缓冲过程中t时刻的速度,z′(t-1)表示电控阻尼执行器在缓冲过程中t-1时刻的速度,设定初始值为z′(0)=0;
步骤9:利用式(17)建立冲击缓冲回弹过程中的电控阻尼悬架动力学模型:
mz″(t)+cz′(t)+fyr+mg+kz(t)=0 (17)
式(17)中,fyr为冲击缓冲回弹过程中电控阻尼执行器的可控阻尼力,并有:
式(18)中,Fyr为冲击缓冲回弹过程中抵消重力之后的电控阻尼执行器的等效阻尼力,并有:
Fyr=fyr+mg (19)
步骤10:利用传感器获取冲击缓冲回弹过程的初始状态,包括:速度z(tsc)与位移z′(tsc),从而按照步骤2的过程,利用冲击缓冲回弹过程的初始状态计算式(18),得到冲击缓冲回弹过程中电控阻尼执行器的运动状态;
步骤11:利用式(20)建立冲击缓冲回弹过程的期望终止状态:
式(11)中,zdesiring(tsr)表示冲击缓冲回弹过程的终止状态的期望位移,z′desiring(tsr)表示冲击缓冲回弹过程的终止状态的期望速度;tsr表示冲击缓冲回弹过程的时间并按照步骤3过程计算得到;
步骤12:按照步骤4的过程计算冲击缓冲回弹过程中的最优宾汉数,从而按照步骤6的过程得到冲击缓冲回弹过程中电控阻尼执行器的期望阻尼力Fdesiring_r;并按照步骤7的过程得到冲击缓冲回弹过程中电控阻尼执行器的期望力对应的最优输出功率Ir_op。
与现有技术相比,本发明的有益效果在于:
1、本发明通过控制电控阻尼悬架***的可变阻尼力实现冲击缓冲控制,在减少了控制变量的同时,仍能够充分利用电控阻尼悬架的可用行程,最大程度降低了冲击载荷。
2、诸如混合控制器和恒定力控制器,不仅需要位移、速度、加速度等传感器实时地采集处理数据,对控制***的硬件要求比较高,同时,传感器数据实时的回馈与计算,对控制器数据处理的能力提出了很高的要求。而本发明采用最优宾汉数冲击缓冲控制策略,在***簧载质量确定后,仅需通过传感器获取冲击缓冲回弹与压缩过程的初始状态,便能够实现精确的冲击缓冲控制,极大地降低了冲击缓冲控制的硬件要求。
3、本发明保持电控阻尼悬架***在冲击缓冲压缩与回弹过程中可控阻尼力不变,受被动阻尼力影响,电控阻尼悬架输出阻尼力与速度变化呈正相关,降低了当执行器速度较小时力跟踪误差对控制效果带来的影响。
附图说明
图1为本发明控制原理图;
图2为本发明方法的流程图;
图中标号,1车身,2位移传感器,3力传感器,4电控阻尼悬架,5半轴,6轮胎。
具体实施方式
本实施例中,如图1所示,一种汽车电控阻尼悬架***冲击工况最优宾汉数缓冲控制方法,如图1所示,其依附硬件包括车辆电控阻尼悬架***以及位移传感器与力传感器,车辆悬架***包括:车身1、电控阻尼悬架4、半轴5、轮胎6;传感器***包括:位移传感器2、力传感器3;其中,车身1构成簧载质量。电控阻尼悬架4、半轴5、轮胎6共同构成了非簧载质量。电控阻尼悬架4由电控阻尼执行器与螺旋弹簧构成。其控制方法包括建立电控阻尼悬架***非线性力学模型;建立冲击缓冲状态下的悬架***动力学模型;基于电控阻尼悬架的冲击缓冲控制控制阻尼力计算;利用电控阻尼执行器非线性力学模型跟踪期望力,计算与最优宾汉数匹配最优输出功率。当控制开始时,位移传感器获取实时位移信息,并微分得到实时速度信息。力传感器采集电控阻尼悬架实际输出力数据并与期望力相比较,以此判断控制策略是否成功实施。具体的说,如图2所示,按如下步骤进行:
步骤1:利用式(1)建立冲击缓冲压缩过程中的电控阻尼悬架动力学模型:
mz″(t)+cz′(t)-fyc+mg+kz(t)=0 (1)
式(1)中,m为簧载质量;z″(t)、z′(t)与z(t)是电控阻尼执行器在缓冲过程中t时刻的加速度、速度与位移;t表示时刻;c表示电控阻尼执行器的阻尼系数;fyc为冲击缓冲压缩过程中电控阻尼执行器的可控阻尼力;g为重力加速度;k为弹簧刚度系数;将式(1)两边同时除以m,并用wn与ζ表示:
Fyc=fyc-mg (3)
步骤2:利用冲击缓冲压缩过程的初始状态计算式(2),得到如式(4)-式(6)所示的冲击缓冲压缩过程中电控阻尼执行器的运动状态方程,其中式(5)与式(6)分别为式(4)的一阶微分与二阶微分:
式(4)-式(6)中,wd为汽车电控阻尼悬架***中有阻尼固有频率;P1和P2为待求解的参数;并通过将冲击缓冲压缩过程的初始状态带入式(4)和式(5)中得到:
式(7)中,v0表示冲击缓冲压缩过程的初始状态中的初始速度,且初始移位为零;
步骤3:利用式(8)建立电控阻尼悬架***冲击缓冲压缩过程的期望终止状态:
zdesiring(tsc)=-S,z′desiring(tsc)=0 (8)
式(8)中,zdesiring(tsc)表示冲击缓冲压缩过程终止状态的期望位移,z′desiring(tsc)表示冲击缓冲压缩过程终止状态的期望速度,S表示冲击缓冲压缩过程的终止时,电控阻尼执行器的期望位置;tsc表示冲击缓冲压缩过程的时间,将z′desiring(tsc)=0带入公式(5)求得tsc:
步骤4:将tsc与zdesiring(tsc)=-S带入式(4)建立抵消重力后的最优宾汉数Bio c+:
设定牛顿迭代法的迭代终止条件为从而利用牛顿迭代法对式(10)进行求解并得到最终迭代出的抵消重力后的最优宾汉数Bio *c+;其中,表示第n次迭代时抵消重力后的最优宾汉数,表示第n+1次迭代时抵消重力后的最优宾汉数;δ表示终止阈值;进而利用式(11)得到冲击缓冲压缩过程中的最优宾汉数
步骤5:对汽车电控阻尼悬架***中的电控阻尼执行器进行数学建模,并通过拟合实验数据进行参数辨识,从而得到如式(12)-式(14)所示的电控阻尼执行器的非线性力学特性的数学模型:
c=c1I2+c2I+c3 (13)
α=α1I2+α2I+α3 (14)
式(12)-式(14)中,Fmodel(t)为电控阻尼执行器数学模型所建立的力学特性;α表示为电控阻尼执行器的磁滞系数;c1、c2、c3表示三个阻尼参数,并用于拟合阻尼系数c;α1、α2、α3表示三个磁滞参数,并用于拟合磁滞系数α;为磁滞输出系数;f为初始状态下电控阻尼执行器的偏置力;I为电控阻尼执行器的输出功率;
步骤6:利用式(15)计算冲击缓冲压缩过程中电控阻尼执行器的期望阻尼力Fdesiring_c:
Fdesiring_c=Bio ccv0 (15)
求出利用式(16)得到t时刻的力跟踪的误差Ferror(t):
Ferror(t)=Fmodel(t)-Fdesiring_c=AIc 2+BIc+C (16)
步骤7:利用电控阻尼执行器数学模型所建立的力学特性Fmodel(t)拟合期望阻尼力Fdesiring_c,利用式(12)对期望阻尼力Fdesiring_c进行力跟踪,利用二次函数求最优解方法使力跟踪误差Ferror(t)最小,从而得到冲击缓冲压缩过程中电控阻尼执行器的期望力对应的最优输出功率Ic_op;
步骤8:利用传感器获取冲击缓冲过程电控阻尼执行器的运动状态,若|z′(t)-z′(t-1)|<ε,则判定冲击缓冲压缩过程完成,并切换为冲击缓冲回弹状态后,执行步骤9,否则返回步骤7执行,其中,z′(t)表示电控阻尼执行器4在缓冲过程中t时刻的速度,z′(t-1)表示电控阻尼执行器4在缓冲过程中t-1时刻的速度,设定初始值为z′(0)=0;
步骤9:利用式(17)建立冲击缓冲回弹过程中的电控阻尼悬架动力学模型:
mz″(t)+cz′(t)+fyr+mg+kz(t)=0 (17)
式(17)中,fyr为冲击缓冲回弹过程中电控阻尼执行器的可控阻尼力,并将式(17)变换为式(18),变换过程同步骤2:
式(18)中,Fyr为冲击缓冲回弹过程中抵消重力之后的电控阻尼执行器的等效阻尼力,并有:
Fyr=fyr+mg (19)
步骤10:利用传感器获取电控阻尼执行器冲击缓冲回弹过程的初始状态及冲击缓冲压缩过程的终止状态,包括:速度z(tsc)与位移z′(tsc),从而按照步骤3的过程,利用冲击缓冲回弹过程的初始状态计算式(18),从而得到如式(20)-(22)所示的冲击缓冲回弹过程中电控阻尼执行器的运动状态;
式(20)-式(21)中,P3和P4为待求解的参数;并通过将冲击缓冲回弹过程的初始状态带入式(20)和式(21)中得到:
步骤11:利用式(25)建立电控阻尼悬架***冲击缓冲回弹过程的期望终止状态:
式(25)中,zdesiring(tsr)表示冲击缓冲回弹过程的终止状态的期望位移,z′desiring(tsr)表示冲击缓冲回弹过程的终止状态的期望速度;tsr表示冲击缓冲回弹过程的时间并按照步骤(4)过程计算得到。tsr也可表示为式(26):
步骤12:按照步骤4的过程计算冲击缓冲回弹过程中的最优宾汉数,从而按照步骤6的过程得到冲击缓冲回弹过程中电控阻尼执行器的期望阻尼力Fdesiring_r;并按照步骤7的过程得到冲击缓冲回弹过程中电控阻尼执行器的期望力对应的最优输出功率Ir_op。
Claims (1)
1.一种汽车电控阻尼悬架***在冲击工况下最优宾汉数控制方法,其特征是按如下步骤进行:
步骤1:利用式(1)建立冲击缓冲压缩过程中的电控阻尼悬架动力学模型:
mz″(t)+cz′(t)-fyc+mg+kz(t)=0 (1)
式(1)中,m为簧载质量;z″(t)、z′(t)与z(t)是电控阻尼执行器在缓冲过程中t时刻的加速度、速度与位移;t表示时刻;c表示电控阻尼执行器的阻尼系数;fyc为冲击缓冲压缩过程中电控阻尼执行器的可控阻尼力;g为重力加速度;k为弹簧刚度系数并有:
Fyc=fyc-mg (3)
步骤2:利用冲击缓冲压缩过程的初始状态计算式(2),得到如式(4)-式(6)所示的冲击缓冲压缩过程中电控阻尼执行器的运动状态方程:
式(4)-式(6)中,wd为汽车电控阻尼悬架***中有阻尼固有频率;P1和P2为待求解的参数;并通过将冲击缓冲压缩过程的初始状态带入式(4)和式(5)中得到:
式(7)中,v0表示冲击缓冲压缩过程的初始状态中的初始速度,且初始移位为零;
步骤3:利用式(8)建立冲击缓冲压缩过程的期望终止状态:
zdesiring(tsc)=-S,z′desiring(tsc)=0 (8)
式(8)中,zdesiring(tsc)表示冲击缓冲压缩过程的终止状态的期望位移,z′desiring(tsc)表示冲击缓冲压缩过程的终止状态的期望速度,S表示冲击缓冲压缩过程的终止时,电控阻尼执行器的期望位置;tsc表示冲击缓冲压缩过程的时间,并有:
步骤4:利用式(10)建立抵消重力后的最优宾汉数Bio c+:
设定牛顿迭代法的迭代终止条件为从而利用牛顿迭代法对式(10)进行求解并得到最终迭代出的抵消重力后的最优宾汉数Bio *c+;其中,Bio c+ n表示第n次迭代时抵消重力后的最优宾汉数,Bio c+ n+1表示第n+1次迭代时抵消重力后的最优宾汉数;δ表示终止阈值;进而利用式(11)得到冲击缓冲压缩过程中的最优宾汉数
步骤5:对汽车电控阻尼悬架***中的电控阻尼执行器进行数学建模,并通过拟合实验数据进行参数辨识,从而得到如式(12)-式(14)所示的电控阻尼执行器的非线性力学特性的数学模型:
c=c1I2+c2I+c3 (13)
α=α1I2+α2I+α3 (14)
式(12)-式(14)中,Fmodel(t)为电控阻尼执行器数学模型建立的力学特性;α表示为电控阻尼执行器的磁滞系数;c1、c2、c3表示三个阻尼参数,并用于拟合阻尼系数c;α1、α2、α3表示三个磁滞参数,并用于拟合磁滞系数α;为磁滞输出系数;f为初始状态下电控阻尼执行器的偏置力;I为电控阻尼执行器的输出功率;
步骤6:利用式(15)计算冲击缓冲压缩过程中电控阻尼执行器的期望阻尼力Fdesiring_c:
Fdesiring_c=Bio ccv0 (15)
利用式(16)得到t时刻的力跟踪的误差Ferror(t):
Ferror(t)=Fmodel(t)-Fdesiring_c=AIc 2+BIc+C (16)
步骤7:利用式(12)对冲击缓冲压缩过程期望阻尼力Fdesiring_c进行力跟踪,并使得误差Ferror(t)最小,从而得到冲击缓冲压缩过程中电控阻尼执行器的期望力对应的最优输出功率Ic_op;
步骤8:若|z′(t)-z′(t-1)|<ε,则判定冲击缓冲压缩过程完成,并切换为冲击缓冲回弹状态后,执行步骤9,否则返回步骤7执行,其中,z′(t)表示电控阻尼执行器在缓冲过程中t时刻的速度,z′(t-1)表示电控阻尼执行器在缓冲过程中t-1时刻的速度,设定初始值为z′(0)=0;
步骤9:利用式(17)建立冲击缓冲回弹过程中的电控阻尼悬架动力学模型:
mz″(t)+cz′(t)+fyr+mg+kz(t)=0 (17)
式(17)中,fyr为冲击缓冲回弹过程中电控阻尼执行器的可控阻尼力,并有:
式(18)中,Fyr为冲击缓冲回弹过程中抵消重力之后的电控阻尼执行器的等效阻尼力,并有:
Fyr=fyr+mg (19)
步骤10:利用传感器获取冲击缓冲回弹过程的初始状态,包括:速度z(tsc)与位移z′(tsc),从而按照步骤2的过程,利用冲击缓冲回弹过程的初始状态计算式(18),得到冲击缓冲回弹过程中电控阻尼执行器的运动状态;
步骤11:利用式(20)建立冲击缓冲回弹过程的期望终止状态:
式(11)中,zdesiring(tsr)表示冲击缓冲回弹过程的终止状态的期望位移,z′desiring(tsr)表示冲击缓冲回弹过程的终止状态的期望速度;tsr表示冲击缓冲回弹过程的时间并按照步骤3过程计算得到;
步骤12:按照步骤4的过程计算冲击缓冲回弹过程中的最优宾汉数,从而按照步骤6的过程得到冲击缓冲回弹过程中电控阻尼执行器的期望阻尼力Fdesiring_r;并按照步骤7的过程得到冲击缓冲回弹过程中电控阻尼执行器的期望力对应的最优输出功率Ir_op。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110339869.9A CN112948979B (zh) | 2021-03-30 | 2021-03-30 | 汽车电控阻尼悬架***冲击工况最优宾汉数控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110339869.9A CN112948979B (zh) | 2021-03-30 | 2021-03-30 | 汽车电控阻尼悬架***冲击工况最优宾汉数控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112948979A CN112948979A (zh) | 2021-06-11 |
CN112948979B true CN112948979B (zh) | 2022-08-30 |
Family
ID=76230527
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110339869.9A Active CN112948979B (zh) | 2021-03-30 | 2021-03-30 | 汽车电控阻尼悬架***冲击工况最优宾汉数控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112948979B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113378408B (zh) * | 2021-07-01 | 2022-09-13 | 合肥工业大学 | 一种电控悬架整车耦合的最优控制方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3208402A (en) * | 1962-10-31 | 1965-09-28 | Sidney H Bingham | Railroad and rapid transit car and undercarriage suspension system |
CN103195858A (zh) * | 2013-04-15 | 2013-07-10 | 沈阳航空航天大学 | 一种面向冲击负荷的磁流变缓冲单元结构及其控制方法 |
CN105260530A (zh) * | 2015-09-29 | 2016-01-20 | 合肥工业大学 | 一种磁流变阻尼器逆模型的建模方法及其应用 |
CN108170870A (zh) * | 2016-12-07 | 2018-06-15 | 黑龙江傲立辅龙科技开发有限公司 | 一种基于电流变材料的减振器设计及控制策略 |
CN109677305A (zh) * | 2019-02-26 | 2019-04-26 | 合肥工业大学 | 一种半主动可控座椅悬架的振动和冲击混合控制方法 |
CN111152616A (zh) * | 2020-01-08 | 2020-05-15 | 合肥工业大学 | 一种磁流变阻尼悬架及其测控方法 |
CN111959218A (zh) * | 2020-07-24 | 2020-11-20 | 中北大学 | 一种基于阻尼力闭环控制的半主动悬架天棚阻尼控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1255935A1 (en) * | 2000-01-31 | 2002-11-13 | Delphi Technologies, Inc. | Tuneable steering damper using magneto-rheological fluid |
-
2021
- 2021-03-30 CN CN202110339869.9A patent/CN112948979B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3208402A (en) * | 1962-10-31 | 1965-09-28 | Sidney H Bingham | Railroad and rapid transit car and undercarriage suspension system |
CN103195858A (zh) * | 2013-04-15 | 2013-07-10 | 沈阳航空航天大学 | 一种面向冲击负荷的磁流变缓冲单元结构及其控制方法 |
CN105260530A (zh) * | 2015-09-29 | 2016-01-20 | 合肥工业大学 | 一种磁流变阻尼器逆模型的建模方法及其应用 |
CN108170870A (zh) * | 2016-12-07 | 2018-06-15 | 黑龙江傲立辅龙科技开发有限公司 | 一种基于电流变材料的减振器设计及控制策略 |
CN109677305A (zh) * | 2019-02-26 | 2019-04-26 | 合肥工业大学 | 一种半主动可控座椅悬架的振动和冲击混合控制方法 |
CN111152616A (zh) * | 2020-01-08 | 2020-05-15 | 合肥工业大学 | 一种磁流变阻尼悬架及其测控方法 |
CN111959218A (zh) * | 2020-07-24 | 2020-11-20 | 中北大学 | 一种基于阻尼力闭环控制的半主动悬架天棚阻尼控制方法 |
Non-Patent Citations (5)
Title |
---|
A novel semi-active control strategy based on the quantitative feedback theory for a vehicle suspension system with magneto-rheological damper saturation;R.Jeyasenthil .etal;《Mechatronics》;20181031;第54卷;36-51页 * |
一种改进的磁流变阻尼器模型及其对振动主动控制的应用;王唯 等;《力学季刊》;20080630;第29卷(第02期);194-198页 * |
汽车电控半主动悬架的混合仿真;罗浩铭;《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》;20121015(第10期);C035-197页 * |
磁流变半主动横向稳定杆对汽车侧倾的影响;白先旭 等;《机械工程学报》;20191205;第55卷(第24期);145-152页 * |
车辆磁流变半主动悬架混合阻尼模糊电流控制研究;胡国良 等;《现代制造工程》;20181031(第10期);94-101页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112948979A (zh) | 2021-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110154666B (zh) | 一种可实现路况预测的车辆悬架***自适应反推控制方法 | |
Williams et al. | Active suspension control to improve vehicle ride and handling | |
CN108891220B (zh) | 一种车辆半主动悬架***的天棚控制改进方法 | |
CN112339517A (zh) | 一种半主动悬架控制方法及控制*** | |
JPH05221223A (ja) | 懸架装置のエンドストップ衝突を低減する装置 | |
Rao et al. | Modeling and simulation of quarter car semi active suspension system using LQR controller | |
CN108859648B (zh) | 一种悬架减震器阻尼控制切换加权系数确定方法 | |
CN107976904B (zh) | 磁流变半主动悬架泰勒级数-二重h2时滞补偿控制方法 | |
CN112948979B (zh) | 汽车电控阻尼悬架***冲击工况最优宾汉数控制方法 | |
Kawabe et al. | New semi-active suspension controller design using quasi-linearization and frequency shaping | |
Du et al. | Constrained H∞ control of active suspension for a half-car model with a time delay in control | |
CN113467233B (zh) | 一种基于模糊模型的时滞有限频域输出反馈控制方法 | |
Onat et al. | LPV gain-scheduling controller design for a non-linear quarter-vehicle active suspension system | |
Sulaiman et al. | Groundhook control of semi-active suspension for heavy vehicle | |
Samsuria et al. | Enhanced sliding mode control for a nonlinear active suspension full car model | |
Najm et al. | Mathematical Modelling and PID Controller Implementation to Control Linear and Nonlinear Quarter Car Active Suspension | |
Rao et al. | Modelling and vibration control of suspension system for automobiles using LQR and PID controllers | |
Vaishnav et al. | Mathematical modelling and comparison of two degree of freedom suspension system of quarter car | |
Lauwerys et al. | Model free control design for a semi-active suspension of a passenger car | |
Lauwerys et al. | A model free control design approach for a semi-active suspension of a passenger car | |
Sandage et al. | Simulation analysis of 2dof quarter car semi-active suspension system to improve ride comfort-A review | |
CN113183711B (zh) | 一种振动与冲击融合控制方法 | |
Yuan et al. | Analysis and tests of magnetorheological adaptive shock mitigation systems | |
Mesa et al. | ROBUST CONTROL H∞ OVER AN ACTIVE DAMPING SYSTEM OF A VEHICLE ROOM | |
Bakar et al. | Ride comfort performance of electric vehicle conversion with active suspension system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |