CN113090704A - 一种磁流变阻尼器在冲击缓冲体系中的抗干扰系统及方法 - Google Patents
一种磁流变阻尼器在冲击缓冲体系中的抗干扰系统及方法 Download PDFInfo
- Publication number
- CN113090704A CN113090704A CN202110381914.7A CN202110381914A CN113090704A CN 113090704 A CN113090704 A CN 113090704A CN 202110381914 A CN202110381914 A CN 202110381914A CN 113090704 A CN113090704 A CN 113090704A
- Authority
- CN
- China
- Prior art keywords
- magnetic induction
- damper
- induction intensity
- magnetorheological damper
- magnetorheological
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 239000007853 buffer solution Substances 0.000 title claims description 26
- 230000006698 induction Effects 0.000 claims abstract description 150
- 238000013016 damping Methods 0.000 claims abstract description 60
- 230000005284 excitation Effects 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 19
- 239000012530 fluid Substances 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 14
- 230000003139 buffering effect Effects 0.000 claims description 6
- 230000033228 biological regulation Effects 0.000 claims description 5
- 230000000630 rising effect Effects 0.000 claims description 5
- 125000004432 carbon atom Chemical group C* 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 3
- 230000004044 response Effects 0.000 description 10
- 239000000872 buffer Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000035939 shock Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000020169 heat generation Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000021715 photosynthesis, light harvesting Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 230000002277 temperature effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
- F16F—SPRINGS; SHOCK-ABSORBERS; MEANS FOR DAMPING VIBRATION
- F16F15/00—Suppression of vibrations in systems; Means or arrangements for avoiding or reducing out-of-balance forces, e.g. due to motion
- F16F15/02—Suppression of vibrations of non-rotating, e.g. reciprocating systems; Suppression of vibrations of rotating systems by use of members not moving with the rotating systems
- F16F15/03—Suppression of vibrations of non-rotating, e.g. reciprocating systems; Suppression of vibrations of rotating systems by use of members not moving with the rotating systems using magnetic or electromagnetic means
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
- F16F—SPRINGS; SHOCK-ABSORBERS; MEANS FOR DAMPING VIBRATION
- F16F9/00—Springs, vibration-dampers, shock-absorbers, or similarly-constructed movement-dampers using a fluid or the equivalent as damping medium
- F16F9/32—Details
- F16F9/50—Special means providing automatic damping adjustment, i.e. self-adjustment of damping by particular sliding movements of a valve element, other than flexions or displacement of valve discs; Special means providing self-adjustment of spring characteristics
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
- F16F—SPRINGS; SHOCK-ABSORBERS; MEANS FOR DAMPING VIBRATION
- F16F9/00—Springs, vibration-dampers, shock-absorbers, or similarly-constructed movement-dampers using a fluid or the equivalent as damping medium
- F16F9/32—Details
- F16F9/53—Means for adjusting damping characteristics by varying fluid viscosity, e.g. electromagnetically
- F16F9/535—Magnetorheological [MR] fluid dampers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Mechanical Engineering (AREA)
- Operations Research (AREA)
- Electromagnetism (AREA)
- Acoustics & Sound (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computing Systems (AREA)
- Fluid-Damping Devices (AREA)
- Vibration Prevention Devices (AREA)
Abstract
本发明提供的抗干扰系统及方法,以磁流变可调阻尼力为依据,根据磁流变阻尼器的设计参数计算相应的目标磁感应强度;力传感器和速度传感器实时检测阻尼力和速度,根据Bingham力学模型计算出磁流变阻尼器线圈产生的实际的磁感应强度;以目标磁感应强度与实际磁感应强度的误差作为滑模控制器的输入;将所述滑模控制器的输出传送至电流驱动器,电流驱动器产生相应的电流通入磁流变阻尼器的励磁线圈,线圈产生的磁感应强度和滑模控制器输出的控制电压作为扩展状态观测器的输入,观测器的输出为对系统扰动的估计,估计的总扰动作为前馈补偿给滑模磁感应强度控制器,构成磁感应强度闭环控制系统,使得实际磁感应强度跟踪目标磁感应强度,实现抗干扰控制。
Description
技术领域
本发明涉及一种磁流变阻尼器在冲击缓冲体系中的抗干扰系统及方法,属于冲击缓冲控制领域。
背景技术
磁流变阻尼器是利用磁流变液在外加磁场作用下可在毫秒量级内实现液体、粘滞体与半固体之间的连续且可逆变换而制成的一种新型的结构半主动控制的智能装置,由于具有结构简单、阻尼力连续顺逆可控且可调范围大、响应速度快以及能耗低等优点,正逐渐应用于车辆工程、土木工程以及武器系统等领域。然而,在冲击载荷下磁流变阻尼器会因能量耗散引起的发热而经历较大的温度变化,温度T会直接影响磁流变液动力粘度η,另一方面,输出阻尼力是动力粘度η的函数,温度T的变化会导致输出阻尼力的变化,因此温度效应对磁流变阻尼器的动力学性能影响不可忽视。在冲击缓冲系统中磁流变阻尼器不仅仅只存在温度扰动,还有其他的一些扰动,如磁流变阻尼器固有的磁滞非线性、在冲击载荷作用下外部环境中未知的扰动和磁流变阻尼力会因沉降、漏油等不确定干扰而发生变化等等,这些扰动都可能造成控制效果不理想。所以,抗扰动的控制策略和控制算法是解决磁流变阻尼在冲击缓冲系统控制效果不理想的可能解决途径之一。
发明内容
本发明为了解决现有技术中存在的上述缺陷和不足,提供了磁流变阻尼器在冲击缓冲体系中的抗干扰系统及方法。
为解决上述技术问题,本发明提供一种磁流变阻尼器在冲击缓冲体系中的抗干扰系统,包括磁感应强度闭环控制模块和干扰调节控制模块;
其中,所述磁感应强度闭环控制模块,将磁流变阻尼器安装在冲击缓冲系统中,由冲击力Fpt和冲击缓冲系统的要求的阻尼力计算出磁流变阻尼器的目标可调阻尼力Fτ_desired;再将磁流变阻尼器的参数与磁流变阻尼器的目标可调阻尼力Fτ_desired相结合,获得目标磁感应强度Bdesired;
再采用力传感器和速度传感器实时监测磁流变阻尼器的输出阻尼力FMRD(t)和磁流变阻尼器活塞杆移动的速度v(t),再根据Bingham力学模型计算出磁流变阻尼器线圈产生的实际磁感应强度Bmeasure;
将目标磁感应强度Bdesired和实际磁感应强度Bmeasured之间的误差e作为滑模磁感应强度控制器的输入,控制电压uc(t)作为所述滑模磁感应强度控制器的输出,控制电压由电流驱动器转化为电流后作为磁流变阻尼器的控制量输入磁流变阻尼器的励磁线圈;由此,构成了磁感应强度闭环控制组件;
所述干扰调节控制模块,将实际磁感应强度Bmeasured和控制电压uc(t)输入扰动观测器中,由扰动观测器实时输出影响磁流变阻尼器的阻尼力的扰动量将扰动量作为前馈补偿引入滑模磁感应强度控制器中,得到补偿后的控制电压;滑模磁感应强度控制器根据补偿后的控制电压自动调节磁流变阻尼器的实际磁感应强度Bmeasured,直至实际磁感应强度Bmeasured无限收敛于目标磁感应强度Bdesired。
另外,本发明还提供一种磁流变阻尼器在冲击缓冲体系中的抗干扰的方法,包括以下步骤:
步骤一,目标磁感应强度Bdesired的计算:首先,将磁流变阻尼器安装在冲击缓冲系统中,由冲击力Fpt和冲击缓冲系统的要求的阻尼力计算出磁流变阻尼器的目标可调阻尼力Fτ_desired;其次,将磁流变阻尼器的参数与磁流变阻尼器的目标可调阻尼力Fτ_desired相结合,计算目标磁感应强度Bdesired;
步骤二,实际磁感应强度Bmeasured的计算:在冲击缓冲系统中采用力传感器和速度传感器实时监测磁流变阻尼器的输出阻尼力FMRD(t)和磁流变阻尼器活塞杆移动的速度v(t),再根据Bingham力学模型计算出磁流变阻尼器线圈产生的实际磁感应强度Bmeasured;
步骤三,抗干扰控制:由磁流变阻尼器的目标磁感应强度Bdesired和实际磁感应强度Bmeasured之间的误差值e构造滑模面s,根据滑模面s设置滑模磁感应强度控制器;将实际磁感应强度Bmeaured和控制电压uc(t)输入扰动观测器中,由扰动观测器实时输出影响磁流变阻尼器的阻尼力的扰动量由并在滑模磁感应强度控制器中引入作为前馈补偿的扰动量得到补偿后的控制电压;磁流变阻尼器根据补偿后的控制电压自动调节补偿前的实际磁感应强度获得补偿后的实际磁感应强度,使得补偿后实际磁感应强度向目标磁感器强度Bdesired收敛;整定滑模磁感应强度控制器参数,重复上述过程,直至实际磁感应强度无限收敛于目标磁感应强度,实现磁流变阻尼器的抗干扰控制。
其中,步骤一的具体过程如下:
1.1根据磁流变液体的屈服剪切应力τy与磁感应强度B之间的τy~B特性曲线进行拟合,得到磁感应强度B和屈服剪切应力τy之间的函数关系式:
B=f(τy)
1.2由磁流变阻尼器的结构和工作原理可知,磁流变阻尼器的目标可调阻尼力Fτ_desired为:Fτ_desired=c2τysgn[v(t)],其中,v(t)为磁流变阻尼器活塞杆移动的速度;sgn[·]为符号函数;c2为磁流变阻尼器的阻尼系数,计算公式为:其中,L为活塞的长度,D为缸筒内径,d为活塞直径,h为阻尼通道宽度;
1.3结合步骤1.1和步骤1.2得到磁流变阻尼器的目标磁感应强度Bdesired的计算公式:
其中,步骤二的具体计算过程如下:
根据Bingham力学模型可知:
FMRD(t)=Fη+Fτ_desired
其中,Fη=c1v(t),Fτ_desired=c2τysgn[v(t)],c1和c2为与磁流变阻尼器结构参数相关的阻尼系数,其值分别为: L为活塞的长度,D为缸筒内径,d为活塞直径,h为阻尼通道宽度,η是流体的动力粘度为常数;
因此,磁流变阻尼器线圈产生的实际磁感应强度Bmeasured为:
上式中,为控制增益,其中是系统增益,λ是常数,U0是上升阶跃电压值,电阻R=(Rc+R0),电阻R0测量电阻,电阻Rc为磁流变阻尼器的线圈的等效为电阻;T1表示输出磁感应强度Bmeasured响应时间常数,电感Lc为磁流变阻尼器的线圈的等效为电阻;p是扩展状态观测器的带宽。
进一步,磁流变阻尼器的目标磁感应强度Bdesired与实际磁感应强度Bmeasured的误差e=Bdesired-Bmeasured,则构造的滑模面s为:
其中,c为增益系数,且大于0。
进一步,滑模磁感应强度控制器输出补偿后的控制电压的过程如下:
本发明所达到的有益技术效果:本发明提供的一种磁流变阻尼器在冲击缓冲体系中抗扰动系统及方法,抑制磁流变阻尼器内外部扰动对冲击缓冲控制系统的负面影响,满足冲击缓冲振动控制系统对阻尼力响应速度和控制精度的要求。
附图说明
图1本发明中磁流变阻尼器励磁线圈的等效测量电路;
图2本发明中磁流变阻尼器励磁线圈等效电路的响应过程示意图;
图3本发明之磁流变阻尼器在冲击缓冲体系中的抗扰动的原理框图。
具体实施方式
下面结合具体实施例对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
下面结合附图和实施例对本发明专利进一步说明。
如图3所示,本发明提供一种磁流变阻尼器在冲击缓冲体系中的抗干扰系统,包括磁感应强度闭环控制模块和干扰调节控制模块;
其中,所述磁感应强度闭环控制模块,将磁流变阻尼器安装在冲击缓冲系统中,由冲击力Fpt和冲击缓冲系统的要求的阻尼力计算出磁流变阻尼器的目标可调阻尼力Fτ_desired;再将磁流变阻尼器的参数与磁流变阻尼器的目标可调阻尼力Fτ_desired相结合,获得目标磁感应强度Bdesired;
再采用力传感器和速度传感器实时监测磁流变阻尼器的输出阻尼力FMRD(t)和磁流变阻尼器活塞杆移动的速度v(t),再根据Bingham力学模型计算出磁流变阻尼器线圈产生的实际磁感应强度Bmeasured;
将目标磁感应强度Bdesired和实际磁感应强度Bmeasured之间的误差e作为滑模磁感应强度控制器的输入,控制电压uc(t)作为所述滑模磁感应强度控制器的输出,控制电压由电流驱动器转化为电流后作为磁流变阻尼器的控制量输入磁流变阻尼器的励磁线圈;由此,构成了磁感应强度闭环控制组件;
所述干扰调节控制模块,将实际磁感应强度Bmeasured和控制电压uc(t)输入扰动观测器中,由扰动观测器实时输出影响磁流变阻尼器的阻尼力的扰动量将扰动量作为前馈补偿引入滑模磁感应强度控制器中,得到补偿后的控制电压;滑模磁感应强度控制器根据补偿后的控制电压自动调节磁流变阻尼器的实际磁感应强度Bmeasured,直至实际磁感应强度Bmeasured无限收敛于目标磁感应强度Bdesired。
另外,本发明还提供一种磁流变阻尼器在冲击缓冲体系中的抗干扰的方法,包括以下步骤:
步骤一,目标磁感应强度Bdesired的计算:首先,将磁流变阻尼器安装在冲击缓冲系统中,由冲击力Fpt和冲击缓冲系统的要求的阻尼力计算出磁流变阻尼器的目标可调阻尼力Fτ_desired;其次,将磁流变阻尼器的参数与磁流变阻尼器的目标可调阻尼力Fτ_desired相结合,计算目标磁感应强度Bdesired;具体计算过程如下:
1.1根据磁流变液体的屈服剪切应力τy与磁感应强度B之间的τy~B特性曲线进行拟合,得到磁感应强度B和屈服剪切应力τy之间的函数关系式:
B=f(τy)
1.2由磁流变阻尼器的结构和工作原理可知,磁流变阻尼器的目标可调阻尼力Fτ_desired为:Fτ_desired=c2τysgn[v(t)],其中,v(t)为磁流变阻尼器活塞杆移动的速度;sgn[·]为符号函数;c2为磁流变阻尼器的阻尼系数,计算公式为:其中,L为活塞的长度,D为缸筒内径,d为活塞直径,h为阻尼通道宽度;
1.3结合步骤1.1和步骤1.2得到磁流变阻尼器的目标磁感应强度Bdesired的计算公式:
步骤二,实际磁感应强度Bmeasured的计算:在冲击缓冲系统中采用力传感器和速度传感器实时监测磁流变阻尼器的输出阻尼力FMRD(t)和磁流变阻尼器活塞杆移动的速度v(t),再根据Bingham力学模型计算出磁流变阻尼器线圈产生的实际磁感应强度Bmeasured;具体计算过程如下:
根据Bingham力学模型可知:
FMRD(t)=Fη+Fτ_desired
其中,Fη=c1v(t),Fτ_desired=c2τysgn[v(t)],c1和c2为与磁流变阻尼器结构参数相关的阻尼系数,其值分别为: L为活塞的长度,D为缸筒内径,d为活塞直径,h为阻尼通道宽度,η是流体的动力粘度为常数;
因此,磁流变阻尼器线圈产生的实际磁感应强度Bmeasured为:
步骤三,抗干扰控制:由磁流变阻尼器的目标磁感应强度Bdesired和实际磁感应强度Bmeasured之间的误差值e构造滑模面s,根据滑模面s设置滑模磁感应强度控制器;磁流变阻尼器的目标磁感应强度Bdesired与实际磁感应强度Bmeasured的误差e=Bdesired-Bmeasured,则构造的滑模面s为:
其中,c为增益系数,且c>0,通过设计合适的增益系数c可以在滑模面内实现实际磁感应强度与目标磁感应强度之间的误差e以最快的速度收敛到零。
将实际磁感应强度Bmeasured和控制电压uc(t)输入扰动观测器中,由扰动观测器实时输出影响磁流变阻尼器的阻尼力的扰动量扰动观测器采用扩展状态观测器,该观测器的输入为磁流变阻尼器的实际磁感应强度Bmeasured和控制电压uc(t),输出为经扰动观测器观测到的扰动量和观测磁感应强度具体计算模型如下:
上式中,为控制增益,其中是系统增益,λ是常数,U0是上升阶跃电压值,电阻R=(Rc+R0),电阻R0测量电阻,电阻Rc为磁流变阻尼器的线圈的等效为电阻;T1表示输出磁感应强度Bmeasured响应时间常数,电感Lc为磁流变阻尼器的线圈的等效为电阻;p是扩展状态观测器的带宽。
磁流变阻尼器根据补偿后的控制电压自动调节补偿前的实际磁感应强度获得补偿后的实际磁感应强度,使得补偿后实际磁感应强度向目标磁感器强度Bdesired收敛;整定滑模磁感应强度控制器参数,重复上述过程,直至实际磁感应强度无限收敛于目标磁感应强度,实现磁流变阻尼器的抗干扰控制。
为了说明本发明提供的抗干扰系统的稳定性,现结合具体实施例进行说明。
一,首先说明观测器的收敛性,因为该扩展状态观测器的收敛性稳定,从而保证了整个抗干扰系统的稳定性,过程如下:
由磁流变阻尼器冲击缓冲体系的二阶非线性方程可知:
将上式替换为:
其中,x1=Bmeasured,x2=d(t);进一步换算得如下表达式:
扩张状态观测器可设计成:
进一步简化得:
取矩阵L,使矩阵(A-LC)稳定,则有e1≥0,从而通过扩张状态观测器的状态可以近似地估计出磁流变阻尼器冲击缓冲二阶非线性系统的所有状态。本发明采用重极点配置的方法使(A-LC)稳定,若要使矩阵(A-LC)稳定,则要使(A-LC)的特征多项的根全部位于复平面S的左半平面。(A-LC)的特征多项为:
上式中l1=2p,l2=p2。作为一个较佳的实施例中,可以按远离虚轴的重极点进行配置。分别是状态磁感应强度B和扰动量d(t)的估计,和会渐进收敛于实际磁感应强度B(t)和d(t),收敛速度取决于参数p。
二,证明本发明的抗干扰系统的稳定性
首先采用实验法建立磁感应强度响应特性理论模型,具体为:
冲击缓冲体系用的磁流变阻尼器的励磁线圈可以等效为电阻Rc和电感Lc串联电路,大功率电阻R0作为测量电阻用来测量两端电压和减小系统的响应时间,电路图如图1所示,其等效阻抗Z=Rc+jwLc,其中wLc是感抗的大小。j表示阻抗角等于90度,物理意义是电感两端的电压相位超前其电流90度。
图1所示电路的运动微分方程为:
对式(5)进行拉普拉斯变换的到测量电路的传递函数:
假定磁流变阻尼器励磁线圈电流产生磁场所需的时间忽略不计,则磁感应强度和电流呈非线性关系,则励磁线圈产生的磁感应强度为:
B(t)=λI(t)
则磁路环节的微分方程为:
图2为电路的整个响应过程,给电路一个输入信号uc(t),由于存在一阶惯性环节,对应的输出信号电流I在数值上要滞后于输入信号一个常量T1,假设线圈电流产生磁场的时间忽略不计,并且不存在磁滞等扰动,利用上述方法可获得有关磁感应强度响应特性的状态方程。
抗干扰系统的稳定性证明:
定义如下Lyapunov函数为V=0.5s2,将滑模磁感应强度控制器输出的控制电压uc(t)代入终端滑膜面s,得到:
对Lyapunov函数求导,可得:
本发明公开的一种在冲击载荷下磁流变阻尼器抗系统及抗干扰方法,以磁流变可调阻尼力为依据,根据磁流变阻尼器的设计参数计算相应的目标磁感应强度;力传感器和速度传感器实时检测阻尼力和速度,根据Bingham力学模型计算出磁流变阻尼器线圈产生的实际的磁感应强度;以目标磁感应强度与实际磁感应强度的误差作为滑模控制器的输入;将所述滑模控制器的输出传送至电流驱动器,电流驱动器产生相应的电流通入磁流变阻尼器的励磁线圈,线圈产生的磁感应强度和滑模控制器输出的控制电压作为扩展状态观测器的输入,观测器的输出为对系统扰动的估计,估计的总扰动作为前馈补偿给滑模磁感应强度控制器,构成磁感应强度闭环控制系统,使得实际磁感应强度跟踪目标磁感应强度,实现抗干扰控制。通过上述方式,本发明提供的一种在冲击载荷下磁流变阻尼器抗干扰方法,抑制磁流变阻尼器内外部扰动对冲击缓冲控制系统响应速度、控制精度的负面影响。
以上已以较佳实施例公布了本发明,然其并非用以限制本发明,凡采取等同替换或等效变换的方案所获得的技术方案,均落在本发明的保护范围内。
Claims (7)
1.一种磁流变阻尼器在冲击缓冲体系中的抗干扰系统,其特征在于:包括磁感应强度闭环控制模块和干扰调节控制模块;
其中,所述磁感应强度闭环控制模块,将磁流变阻尼器安装在冲击缓冲系统中,由冲击力Fpt和冲击缓冲系统的要求的阻尼力计算出磁流变阻尼器的目标可调阻尼力Fτ_desired;再将磁流变阻尼器的参数与磁流变阻尼器的目标可调阻尼力Fτ_desired相结合,获得目标磁感应强度Bdesired;
再采用力传感器和速度传感器实时监测磁流变阻尼器的输出阻尼力FMRD(t)和磁流变阻尼器活塞杆移动的速度v(t),再根据Bingham力学模型计算出磁流变阻尼器线圈产生的实际磁感应强度Bmeasured;
将目标磁感应强度Bdesired和实际磁感应强度Bmeasured之间的误差e作为滑模磁感应强度控制器的输入,控制电压uc(t)作为所述滑模磁感应强度控制器的输出,控制电压由电流驱动器转化为电流后作为磁流变阻尼器的控制量输入磁流变阻尼器的励磁线圈;由此,构成了磁感应强度闭环控制组件;
2.一种磁流变阻尼器在冲击缓冲体系中的抗干扰的方法,其特征在于,包括以下步骤:
步骤一,目标磁感应强度Bdesired的计算:首先,将磁流变阻尼器安装在冲击缓冲系统中,由冲击力Fpt和冲击缓冲系统的要求的阻尼力计算出磁流变阻尼器的目标可调阻尼力Fτ_desired;其次,将磁流变阻尼器的参数与磁流变阻尼器的目标可调阻尼力Fτ_desired相结合,计算目标磁感应强度Bdesired;
步骤二,实际磁感应强度Bmeasured的计算:在冲击缓冲系统中采用力传感器和速度传感器实时监测磁流变阻尼器的输出阻尼力FMRD(t)和磁流变阻尼器活塞杆移动的速度v(t),再根据Bingham力学模型计算出磁流变阻尼器线圈产生的实际磁感应强度Bmeasured;
步骤三,抗干扰控制:由磁流变阻尼器的目标磁感应强度Bdesired和实际磁感应强度Bmeasured之间的误差值e构造滑模面s,根据滑模面s设置滑模磁感应强度控制器;将实际磁感应强度Bmeasured和控制电压uc(t)输入扰动观测器中,由扰动观测器实时输出影响磁流变阻尼器的阻尼力的扰动量由并在滑模磁感应强度控制器中引入作为前馈补偿的扰动量得到补偿后的控制电压;磁流变阻尼器根据补偿后的控制电压自动调节补偿前的实际磁感应强度获得补偿后的实际磁感应强度,使得补偿后实际磁感应强度向目标磁感器强度Bdesired收敛;整定滑模磁感应强度控制器参数,重复上述过程,直至实际磁感应强度无限收敛于目标磁感应强度,实现磁流变阻尼器的抗干扰控制。
3.根据权利要求2所述的磁流变阻尼器在冲击缓冲体系中的抗干扰的方法,其特征在于,步骤一的具体过程如下:
1.1根据磁流变液体的屈服剪切应力τy与磁感应强度B之间的τy~B特性曲线进行拟合,得到磁感应强度B和屈服剪切应力τy之间的函数关系式:
B=f(τy)
1.2由磁流变阻尼器的结构和工作原理可知,磁流变阻尼器的目标可调阻尼力Fτ_desired为:Fτ_desired=c2τysgn[v(t)],其中,v(t)为磁流变阻尼器活塞杆移动的速度;sgn[·]为符号函数;c2为磁流变阻尼器的阻尼系数,计算公式为:其中,L为活塞的长度,D为缸筒内径,d为活塞直径,h为阻尼通道宽度;
1.3结合步骤1.1和步骤1.2得到磁流变阻尼器的目标磁感应强度Bdesired的计算公式:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110381914.7A CN113090704B (zh) | 2021-04-09 | 2021-04-09 | 一种磁流变阻尼器在冲击缓冲体系中的抗干扰系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110381914.7A CN113090704B (zh) | 2021-04-09 | 2021-04-09 | 一种磁流变阻尼器在冲击缓冲体系中的抗干扰系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113090704A true CN113090704A (zh) | 2021-07-09 |
CN113090704B CN113090704B (zh) | 2022-04-15 |
Family
ID=76675807
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110381914.7A Active CN113090704B (zh) | 2021-04-09 | 2021-04-09 | 一种磁流变阻尼器在冲击缓冲体系中的抗干扰系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113090704B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6694856B1 (en) * | 2001-02-22 | 2004-02-24 | The University Of Maryland | Magnetorheological damper and energy dissipation method |
CN1858461A (zh) * | 2005-09-15 | 2006-11-08 | 重庆大学 | 缸体感应式磁流变阻尼器的集成相对位移自传感方法及系统 |
CN101332817A (zh) * | 2008-07-18 | 2008-12-31 | 西安交通大学 | 车辆悬挂系统用磁流变阻尼器控制装置及控制方法 |
CN101576141A (zh) * | 2009-06-04 | 2009-11-11 | 重庆大学 | 基于旋转型磁流变器件的纵向直线振动抑制方法与装置 |
CN111103793A (zh) * | 2018-10-25 | 2020-05-05 | 中车株洲电力机车研究所有限公司 | 基于磁流变阻尼器的火炮后坐缓冲系统控制方法、装置 |
-
2021
- 2021-04-09 CN CN202110381914.7A patent/CN113090704B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6694856B1 (en) * | 2001-02-22 | 2004-02-24 | The University Of Maryland | Magnetorheological damper and energy dissipation method |
CN1858461A (zh) * | 2005-09-15 | 2006-11-08 | 重庆大学 | 缸体感应式磁流变阻尼器的集成相对位移自传感方法及系统 |
CN101332817A (zh) * | 2008-07-18 | 2008-12-31 | 西安交通大学 | 车辆悬挂系统用磁流变阻尼器控制装置及控制方法 |
CN101576141A (zh) * | 2009-06-04 | 2009-11-11 | 重庆大学 | 基于旋转型磁流变器件的纵向直线振动抑制方法与装置 |
CN111103793A (zh) * | 2018-10-25 | 2020-05-05 | 中车株洲电力机车研究所有限公司 | 基于磁流变阻尼器的火炮后坐缓冲系统控制方法、装置 |
Also Published As
Publication number | Publication date |
---|---|
CN113090704B (zh) | 2022-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104111607B (zh) | 一种考虑输入时滞的电机位置伺服系统的控制方法 | |
Zhao et al. | Fractional order control to the electro-hydraulic system in insulator fatigue test device | |
CN103660991B (zh) | 一种常导磁浮列车非线性悬浮控制方法及其控制系统 | |
CN104242769A (zh) | 基于连续终端滑模技术的永磁同步电机速度复合控制方法 | |
CN102360176A (zh) | 基于简化二阶自抗扰控制器的电厂主汽温度控制方法 | |
CN111708276B (zh) | 基于线性状态观测器观测误差补偿的自适应鲁棒控制方法 | |
CN108762083B (zh) | 一种基于加速度观测器的自动控制系统 | |
CN113517832B (zh) | 一种低压伺服离散线性自抗扰控制方法 | |
CN110784149A (zh) | 交流伺服系统机械谐振抑制方法以及系统 | |
Dai et al. | Adaptive force tracking control of electrohydraulic systems with low load using the modified LuGre friction model | |
CN104730922A (zh) | 基于扩张状态观测器的伺服系统线性反馈控制和极点配置确定参数方法 | |
CN108919652A (zh) | 一种自适应抗扰整形控制方法与系统 | |
CN110647111B (zh) | 考虑输出离散性的电液伺服系统非线性自抗扰控制方法 | |
CN108964563A (zh) | 一种基于非光滑控制技术的感应电机直接转矩控制方法 | |
CN113090704B (zh) | 一种磁流变阻尼器在冲击缓冲体系中的抗干扰系统及方法 | |
CN110045603B (zh) | 压电陶瓷驱动部件变载荷环境下的鲁棒自适应控制方法 | |
Piao et al. | Disturbance observer‐based robust motor control enhanced by adaptive neural network in the absence of velocity measurement | |
CN111103793B (zh) | 基于磁流变阻尼器的火炮后坐缓冲系统控制方法、装置 | |
CN112859608A (zh) | 一种基于rbf神经网络补偿的自适应动态面控制方法 | |
CN110323741B (zh) | 基于录波曲线的电力系统低频振荡分析与抑制方法 | |
CN114625006B (zh) | 一种高速大惯量电液位置伺服系统的输出反馈控制方法 | |
CN108448987B (zh) | 三相永磁同步电机扰动感知控制方法 | |
Zhao et al. | Design of MRAC and Modified MRAC for the Turntable | |
CN114077192A (zh) | 一种非线性系统输出反馈智能控制方法 | |
CN114137834B (zh) | 一种基于改进型eso的运动控制系统扰动抑制方法 |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20221201 Address after: Room 10C-A40, Building B, Chuangzhi Building, No. 17 Xinghuo Road, Jiangbei New District, Nanjing, Jiangsu Province, 210032 Patentee after: Nanjing Haosheng Precision Machinery Technology Co.,Ltd. Address before: Longpan road Xuanwu District of Nanjing city of Jiangsu Province, No. 159 210037 Patentee before: NANJING FORESTRY University |
|
TR01 | Transfer of patent right |