CN112948979B - 汽车电控阻尼悬架***冲击工况最优宾汉数控制方法 - Google Patents

汽车电控阻尼悬架***冲击工况最优宾汉数控制方法 Download PDF

Info

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
Application number
CN202110339869.9A
Other languages
English (en)
Other versions
CN112948979A (zh
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.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
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 Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN202110339869.9A priority Critical patent/CN112948979B/zh
Publication of CN112948979A publication Critical patent/CN112948979A/zh
Application granted granted Critical
Publication of CN112948979B publication Critical patent/CN112948979B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling 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为弹簧刚度系数并有:
Figure GDA0003763346450000021
式(2)中,wn表示无阻尼刚度,且
Figure GDA0003763346450000022
ζ表示阻尼比,且ζ<1,Fyc为冲击缓冲压缩过程中抵消重力之后的电控阻尼执行器的等效阻尼力,并有:
Fyc=fyc-mg (3)
步骤2:利用冲击缓冲压缩过程的初始状态计算式(2),得到如式(4)-式(6)所示的冲击缓冲压缩过程中电控阻尼执行器的运动状态方程:
Figure GDA0003763346450000023
Figure GDA0003763346450000024
Figure GDA0003763346450000025
式(4)-式(6)中,wd为汽车电控阻尼悬架***中有阻尼固有频率;P1和P2为待求解的参数;并通过将冲击缓冲压缩过程的初始状态带入式(4)和式(5)中得到:
Figure GDA0003763346450000026
式(7)中,v0表示冲击缓冲压缩过程的初始状态中的初始速度,且初始移位为零;
步骤3:利用式(8)建立冲击缓冲压缩过程的期望终止状态:
zdesiring(tsc)=-S,z′desiring(tsc)=0 (8)
式(8)中,zdesiring(tsc)表示冲击缓冲压缩过程的终止状态的期望位移,z′desiring(tsc)表示冲击缓冲压缩过程的终止状态的期望速度,S表示冲击缓冲压缩过程的终止时,电控阻尼执行器的期望位置;tsc表示冲击缓冲压缩过程的时间,并有:
Figure GDA0003763346450000027
步骤4:利用式(10)建立抵消重力后的最优宾汉数Bio c+
Figure GDA0003763346450000028
设定牛顿迭代法的迭代终止条件为
Figure GDA0003763346450000031
从而利用牛顿迭代法对式(10)进行求解并得到最终迭代出的抵消重力后的最优宾汉数Bio *c+;其中,
Figure GDA0003763346450000039
表示第n次迭代时抵消重力后的最优宾汉数,
Figure GDA0003763346450000038
表示第n+1次迭代时抵消重力后的最优宾汉数;δ表示终止阈值;进而利用式(11)得到冲击缓冲压缩过程中的最优宾汉数
Figure GDA0003763346450000032
Figure GDA0003763346450000033
步骤5:对汽车电控阻尼悬架***中的电控阻尼执行器进行数学建模,并通过拟合实验数据进行参数辨识,从而得到如式(12)-式(14)所示的电控阻尼执行器的非线性力学特性的数学模型:
Figure GDA0003763346450000034
c=c1I2+c2I+c3 (13)
α=α1I22I+α3 (14)
式(12)-式(14)中,Fmodel(t)为电控阻尼执行器数学模型建立的力学特性;α表示为电控阻尼执行器的磁滞系数;c1、c2、c3表示三个阻尼参数,并用于拟合阻尼系数c;α1、α2、α3表示三个磁滞参数,并用于拟合磁滞系数α;
Figure GDA0003763346450000035
为磁滞输出系数;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)
式(16)中,Ic表示电控阻尼执行器的期望力对应的输出功率,
Figure GDA0003763346450000036
Figure GDA0003763346450000037
步骤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为冲击缓冲回弹过程中电控阻尼执行器的可控阻尼力,并有:
Figure GDA0003763346450000041
式(18)中,Fyr为冲击缓冲回弹过程中抵消重力之后的电控阻尼执行器的等效阻尼力,并有:
Fyr=fyr+mg (19)
步骤10:利用传感器获取冲击缓冲回弹过程的初始状态,包括:速度z(tsc)与位移z′(tsc),从而按照步骤2的过程,利用冲击缓冲回弹过程的初始状态计算式(18),得到冲击缓冲回弹过程中电控阻尼执行器的运动状态;
步骤11:利用式(20)建立冲击缓冲回弹过程的期望终止状态:
Figure GDA0003763346450000042
式(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与ζ表示:
Figure GDA0003763346450000051
式(2)中,wn表示无阻尼刚度,且
Figure GDA0003763346450000052
ζ表示阻尼比,半主动冲击缓冲***在不受控状态下为欠阻尼***,因此阻尼比ζ<1,Fyc为冲击缓冲压缩过程中抵消重力之后的电控阻尼执行器4的等效阻尼力,并有:
Fyc=fyc-mg (3)
步骤2:利用冲击缓冲压缩过程的初始状态计算式(2),得到如式(4)-式(6)所示的冲击缓冲压缩过程中电控阻尼执行器的运动状态方程,其中式(5)与式(6)分别为式(4)的一阶微分与二阶微分:
Figure GDA0003763346450000061
Figure GDA0003763346450000062
Figure GDA0003763346450000063
式(4)-式(6)中,wd为汽车电控阻尼悬架***中有阻尼固有频率;P1和P2为待求解的参数;并通过将冲击缓冲压缩过程的初始状态带入式(4)和式(5)中得到:
Figure GDA0003763346450000064
式(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
Figure GDA0003763346450000065
步骤4:将tsc与zdesiring(tsc)=-S带入式(4)建立抵消重力后的最优宾汉数Bio c+
Figure GDA0003763346450000066
设定牛顿迭代法的迭代终止条件为
Figure GDA0003763346450000067
从而利用牛顿迭代法对式(10)进行求解并得到最终迭代出的抵消重力后的最优宾汉数Bio *c+;其中,
Figure GDA0003763346450000068
表示第n次迭代时抵消重力后的最优宾汉数,
Figure GDA0003763346450000069
表示第n+1次迭代时抵消重力后的最优宾汉数;δ表示终止阈值;进而利用式(11)得到冲击缓冲压缩过程中的最优宾汉数
Figure GDA0003763346450000071
Figure GDA0003763346450000072
步骤5:对汽车电控阻尼悬架***中的电控阻尼执行器进行数学建模,并通过拟合实验数据进行参数辨识,从而得到如式(12)-式(14)所示的电控阻尼执行器的非线性力学特性的数学模型:
Figure GDA0003763346450000073
c=c1I2+c2I+c3 (13)
α=α1I22I+α3 (14)
式(12)-式(14)中,Fmodel(t)为电控阻尼执行器数学模型所建立的力学特性;α表示为电控阻尼执行器的磁滞系数;c1、c2、c3表示三个阻尼参数,并用于拟合阻尼系数c;α1、α2、α3表示三个磁滞参数,并用于拟合磁滞系数α;
Figure GDA0003763346450000074
为磁滞输出系数;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)
式(16)中,Ic表示电控阻尼执行器的期望力对应的输出功率,
Figure GDA0003763346450000075
Figure GDA0003763346450000076
步骤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:
Figure GDA0003763346450000081
式(18)中,Fyr为冲击缓冲回弹过程中抵消重力之后的电控阻尼执行器的等效阻尼力,并有:
Fyr=fyr+mg (19)
步骤10:利用传感器获取电控阻尼执行器冲击缓冲回弹过程的初始状态及冲击缓冲压缩过程的终止状态,包括:速度z(tsc)与位移z′(tsc),从而按照步骤3的过程,利用冲击缓冲回弹过程的初始状态计算式(18),从而得到如式(20)-(22)所示的冲击缓冲回弹过程中电控阻尼执行器的运动状态;
Figure GDA0003763346450000082
Figure GDA0003763346450000083
Figure GDA0003763346450000084
式(20)-式(21)中,P3和P4为待求解的参数;并通过将冲击缓冲回弹过程的初始状态带入式(20)和式(21)中得到:
Figure GDA0003763346450000085
Figure GDA0003763346450000086
步骤11:利用式(25)建立电控阻尼悬架***冲击缓冲回弹过程的期望终止状态:
Figure GDA0003763346450000087
式(25)中,zdesiring(tsr)表示冲击缓冲回弹过程的终止状态的期望位移,z′desiring(tsr)表示冲击缓冲回弹过程的终止状态的期望速度;tsr表示冲击缓冲回弹过程的时间并按照步骤(4)过程计算得到。tsr也可表示为式(26):
Figure GDA0003763346450000091
步骤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为弹簧刚度系数并有:
Figure FDA0003763346440000011
式(2)中,wn表示无阻尼刚度,且
Figure FDA0003763346440000012
ζ表示阻尼比,且ζ<1,Fyc为冲击缓冲压缩过程中抵消重力之后的电控阻尼执行器的等效阻尼力,并有:
Fyc=fyc-mg (3)
步骤2:利用冲击缓冲压缩过程的初始状态计算式(2),得到如式(4)-式(6)所示的冲击缓冲压缩过程中电控阻尼执行器的运动状态方程:
Figure FDA0003763346440000013
Figure FDA0003763346440000014
Figure FDA0003763346440000015
式(4)-式(6)中,wd为汽车电控阻尼悬架***中有阻尼固有频率;P1和P2为待求解的参数;并通过将冲击缓冲压缩过程的初始状态带入式(4)和式(5)中得到:
Figure FDA0003763346440000016
式(7)中,v0表示冲击缓冲压缩过程的初始状态中的初始速度,且初始移位为零;
步骤3:利用式(8)建立冲击缓冲压缩过程的期望终止状态:
zdesiring(tsc)=-S,z′desiring(tsc)=0 (8)
式(8)中,zdesiring(tsc)表示冲击缓冲压缩过程的终止状态的期望位移,z′desiring(tsc)表示冲击缓冲压缩过程的终止状态的期望速度,S表示冲击缓冲压缩过程的终止时,电控阻尼执行器的期望位置;tsc表示冲击缓冲压缩过程的时间,并有:
Figure FDA0003763346440000021
步骤4:利用式(10)建立抵消重力后的最优宾汉数Bio c+
Figure FDA0003763346440000022
设定牛顿迭代法的迭代终止条件为
Figure FDA0003763346440000023
从而利用牛顿迭代法对式(10)进行求解并得到最终迭代出的抵消重力后的最优宾汉数Bio *c+;其中,Bio c+ n表示第n次迭代时抵消重力后的最优宾汉数,Bio c+ n+1表示第n+1次迭代时抵消重力后的最优宾汉数;δ表示终止阈值;进而利用式(11)得到冲击缓冲压缩过程中的最优宾汉数
Figure FDA0003763346440000024
Figure FDA0003763346440000025
步骤5:对汽车电控阻尼悬架***中的电控阻尼执行器进行数学建模,并通过拟合实验数据进行参数辨识,从而得到如式(12)-式(14)所示的电控阻尼执行器的非线性力学特性的数学模型:
Figure FDA0003763346440000026
c=c1I2+c2I+c3 (13)
α=α1I22I+α3 (14)
式(12)-式(14)中,Fmodel(t)为电控阻尼执行器数学模型建立的力学特性;α表示为电控阻尼执行器的磁滞系数;c1、c2、c3表示三个阻尼参数,并用于拟合阻尼系数c;α1、α2、α3表示三个磁滞参数,并用于拟合磁滞系数α;
Figure FDA0003763346440000027
为磁滞输出系数;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)
式(16)中,Ic表示电控阻尼执行器的期望力对应的输出功率,
Figure FDA0003763346440000028
Figure FDA0003763346440000029
步骤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为冲击缓冲回弹过程中电控阻尼执行器的可控阻尼力,并有:
Figure FDA0003763346440000031
式(18)中,Fyr为冲击缓冲回弹过程中抵消重力之后的电控阻尼执行器的等效阻尼力,并有:
Fyr=fyr+mg (19)
步骤10:利用传感器获取冲击缓冲回弹过程的初始状态,包括:速度z(tsc)与位移z′(tsc),从而按照步骤2的过程,利用冲击缓冲回弹过程的初始状态计算式(18),得到冲击缓冲回弹过程中电控阻尼执行器的运动状态;
步骤11:利用式(20)建立冲击缓冲回弹过程的期望终止状态:
Figure FDA0003763346440000032
式(11)中,zdesiring(tsr)表示冲击缓冲回弹过程的终止状态的期望位移,z′desiring(tsr)表示冲击缓冲回弹过程的终止状态的期望速度;tsr表示冲击缓冲回弹过程的时间并按照步骤3过程计算得到;
步骤12:按照步骤4的过程计算冲击缓冲回弹过程中的最优宾汉数,从而按照步骤6的过程得到冲击缓冲回弹过程中电控阻尼执行器的期望阻尼力Fdesiring_r;并按照步骤7的过程得到冲击缓冲回弹过程中电控阻尼执行器的期望力对应的最优输出功率Ir_op
CN202110339869.9A 2021-03-30 2021-03-30 汽车电控阻尼悬架***冲击工况最优宾汉数控制方法 Active CN112948979B (zh)

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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113378408B (zh) * 2021-07-01 2022-09-13 合肥工业大学 一种电控悬架整车耦合的最优控制方法

Citations (7)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (7)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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