原创《基于动态迟滞单元的磁流变阻尼器建模和辨识》
该文是阻尼器方面有关开题报告范文跟阻尼器和动态迟滞单元和流变有关参考文献格式范文。
摘 要:磁流变阻尼器(Magnetorheologicaldamper,简称MR阻尼器)在振动控制方面具有很好的应用前景,但是由于MR阻尼器固有的非线性的迟滞特性会影响控制精度,因此需要对磁流变阻尼器进行建模分析.首先提出了改进的动态迟滞单元(Modifieddynamichysteresisoperator,MDHO),MDHO 能够描述不同的输出幅值;然后在MDHO的基础上提出了MR阻尼器的迟滞模型,在该模型中磁流变阻尼器的输出阻尼力表示为MDHO 的滞变阻尼力与粘滞力两部分;最后利用改进的粒子群算法(Modifiedparticleswarmoptimization,MPSO)进行模型参数辨识.实验结果表明,基于动态迟滞单元的磁流变阻尼器的迟滞模型能够有效模拟MR阻尼器的实际特性.
关键词:磁流变阻尼器;改进动态迟滞单元;迟滞模型;改进粒子群算法
中图分类号:TP２０ 文献标志码:A 文章编号:１６７３Ｇ３８５１(２０１８)０３Ｇ０１９５Ｇ０５
０ 引 言
磁流变阻尼器被认为是最具应用前景的振动控制智能驱动器,但是由于本身固有的非线性的迟滞特性限制了它在很多方面的应用.因此要充分利用磁流变阻尼器的低功耗、高阻尼力、响应速度快等优点,需要建立较精确的数学模型.已有研究中模型一般可以分为两大类:参数化模型和非参数化模型.
MR阻尼器的参数化模型是将线性和非线性的弹簧模块、阻尼器以及其他的物理元件集合在一起,基于该类结构的迟滞模型能够模拟磁流变阻尼器非线性的特性[１Ｇ６].已有的MR 阻尼器的参数化模型包括:宾汉模型(Bingham model)[１]以及改进的宾汉模型[２]、非线性双粘性模型[３]、非线性滞回模型、BoucＧWen模型[４Ｇ５]、LuGre模型[６].宾汉模型[１]的结构是将库伦元件和粘塑性模块结合在一起.该模型能够较好地描述MR 阻尼器的位移和阻尼力的关系,已被广泛运用于迟滞的建模中,但是在低速下并不能较好地描述速度和阻尼力的关系.改进的宾汉模型[２]的预测结果与实验结果的差异明显,且数值计算较复杂耗时,不利于工程应用.非线性滞回模型[３]是非线性双粘性模型的改进,由四条线构成的滞环模型,它的局限性在于无法描述阻尼力饱和区域速度和阻尼力的特性.BoucＧWen模型[４Ｇ５]的粘塑性模块和混合的弹簧结构能够较好反映MR阻尼器的动态特性,因此该模型运用非常广泛,但是模型参数需要进行辨识,参数辨识结果对模型影响很大.LuGre模型[６]是一种基于MR 阻尼器内部磁流变液流动粘滞状态行为的阻尼器模型,利用动静力摩擦较准确地解释了磁流变液的粘性特性,但是该模型的参数确定非常困难.参数化模型是在理想化的假设条件下提出的,可能会造成质量为负值等与实际不符的结果.
非参数化模型是对数据进行分析和训练,从而模拟实际磁流变阻尼器特性的模型.非参数化模型包括神经网络模型[７Ｇ８]、模糊模型[９Ｇ１０]等.Chang等[１１]提出了一种多层神经网络微分最优的方法,Wang等[１２]提出了一种递归的神经网络模型.以上两种基于神经网络的方法均能够高精度地描述MR阻尼器的特性,但是神经网络的架构和训练过程很复杂.Truong等[９]提出了基于模糊理论的非参数化的黑箱模型,该模型根据映射方式对MR 阻尼器的特性进行了估计;Liem 等[１３]提出了基于自调节的模糊机理的非参数化模型.非参数化模型能够较好地描述磁流变阻尼器的动态特性,但是该类模型存在实现过程复杂的缺点.
本文提出了一种基于动态迟滞单元的参数化的磁流变阻尼器迟滞模型.首先改进了动态迟滞单元,增加了调整参数,从而改变迟滞特性的高度和偏移,以适应磁流变阻尼器固有的非线性特性,使得模型具有较好的灵活性;然后利用改进的粒子群算法对模型参数进行辨识.实验结果表明,该模型能精确表达磁流变阻尼器的力Ｇ速度特性,适合磁流变阻尼器控制器的设计.
１ 改进动态迟滞单元
本文提出了改进的动态迟滞单元(Modifieddynamichysteresisoperator,MDHO),其结构如图１所示.
其中:U 是速度作为模型输入;Φ(U )是迟滞单元的输出;η 是状态参数;ki 是积分系数;k 是死区特性的斜率;x０ 是死区宽度;b 是迟滞环的偏置.实际情况中磁流变阻尼器在不同的电流下其阻尼力是不同的,表现为迟滞的幅值不同.为了表达不同幅值的迟滞现象,本文增加了死区宽度x０.
在x０ 不同的情况下,当输入U ＝１０?sin(２πt/３０),其迟滞环幅值的变化如图２所示.从图２可以看出,当x０＝１,２,３时,即随着死区宽度的增加,迟滞环的输出幅值逐渐增大.另外参数k、ki 对迟滞环的影响如图３、图４所示.
Φ(U)为迟滞单元的输出,由图３可以看出,死区的斜率可以改变迟滞环小范围的高度和宽度.由图４ 可以看出,积分系数ki 可以改变迟滞环的宽度.
文 献[１４]中的动态单元输出幅值是限定的,而实际情况下磁流变阻尼器中不同电流对应不同的输出幅值,不能直接用于磁流变阻尼器的建模.本文与文献[１４]所提出的迟滞单元不同,通过改变死区宽度x０ 来改变输出的幅值,从而更精确地表达不同电流下的阻尼力Ｇ速度的关系.
２ 模型结构
本文提出的磁流变阻尼器模型结构如图５所示.磁流变阻尼器的输出阻尼力为改进迟滞单元的滞变阻尼力与粘滞力之和,该模型的数学表达式为:
其中:Φ(U )是改进迟滞单元的滞变阻尼力输出;c０?x是粘滞力输出;x 是位移;F 是阻尼力,为模型输出.
电流和速度均会影响磁流变阻尼器的输出,为了使得模型在变电流的情况下能够适应速度的变化,需要确定模型参数和电流的关系.这样,模型参数x０、ki、k、w、b、c０ 可表示为:
其中:x０、ki、k、w、b、c０ 分别表示电场为零时的粘滞阻尼力、积分系数、屈服强度、前馈系数、偏置力系数、粘滞阻尼系数;r 是中间变量,与电流i的关系可以用一阶微分方程表示为:
?r＝－α(r－pi－q),
从而建立参数x０、ki、k、w、b、c０ 与电流i 的联系.因此,模型中有x０１、x０２、ki１、ki２、k１、k２、w１、w２、b１、b２、c０１、c０２、α、p、q 共１５个参数需要辨识.
３ 改进的粒子群算法
粒子群算法中,所有的粒子都被用于寻找最优解,每个粒子寻找最优解的速度跟它当前的位置和全局的位置是相关的,传统的粒子群算法表示为:
其中:Vikg＋１、Vkgi 、Xkgi 、N 跟传统粒子群算法中一样,ωmin和ωmax分别表示为权值的最大值和最小值,ωα是中间变量,ω 是权值,ckg１ 、ckg２ 分别是第kg 个粒子的学习因子,c１min,c２min、c１max,c２max分别表示学习因子的上限和下限,rand１ 是０到１之间的随机数.改进的粒子群算法主要有两点优势,首先增加的ωm 作为Vkgi 和ckg１r１(pkgibest－Xkgi )共同的权值,传统粒子群算法中容易陷入局部最优是由于没有对局部最优解取权值,改进的粒子群算法在ckg１r１(pkgibest－Xkgi )中包含pkgibest,pkgibest能够反映局部搜索情况,对其取权值能克服陷入局部最优的问题;另外ckg１ 、ckg２ 用于保证收敛速度[１５],随着搜寻的进行,学习因子随着采样次数的变化而减小,所以改进粒子群算法有更快的收敛速度.
４ 实验验证
本实验选用美国Lord公司型号为RDＧ１００５Ｇ３的磁流变阻尼器.验证的电流频率为１５Hz,粒子个数N＝５００,学习因子c１ 的最大值和最小值分别为c１min＝１１,c１max＝１３;c２ 的最大值和最小值分别为c２min＝１４,c２max＝１９.权值的最大值和最小值分别为ωmin＝０４,ωmax＝０９.利用改进的粒子群算法辨识出参数分别为x０１＝０２９１０,x０２＝０００３０,ki１ ＝１４９２３７,ki２ ＝－０１４７９,b１ ＝００２６５,b２ ＝－００００３,c０１＝００５８２,c０２＝００００６,k１＝１２４８３０,k２＝０１０６４,w１＝０６７３８,w２＝０００７４,Z＝－３８００００,p＝－９８５００,q＝１２００００.实验结果与模型输出的比较如图６所示.图７分别是实际输出和模型输出的比较,图８为建模误差,可以看出,该模型能比较精确的逼近实
际输出,由实验结果可知相对误差均在７８％以下,满足控制建模的要求.
５ 结 论
本文提出了一种迟滞单元与粘塑性模块结合的磁流变阻尼器模型,将磁流变阻尼器的阻尼力分解成滞变阻尼力与粘滞力之和.改进的迟滞单元通过增加死区宽度来调整输出的幅值,从而更加适应磁流变阻尼器的迟滞特性.运用改进的粒子群算法来辨识模型参数,克服了传统算法容易陷入局部最优和速度慢的缺陷.该模型能精确表达变电流的情况下阻尼力Ｇ速度的关系,结构简单,易于分析,便于控制器设计.
参考文献:
[１]HongSR,WereleyN M,ChoiY T,etal．AnalyticalandexperimentalvalidationofanondimensionalBinghammodelformixedＧmode magnetorheologicaldampers[J]．JournalofSound& Vibration,２００８,３１２(３):３９９Ｇ４１７．
[２]HuW,WereleyN M．RateＧdependentelastoslidemodelformagnetorheologicaldamper[J]．JournalofGuidanceControl& Dynamics,２０１５,３１(３):４７９Ｇ４８９．
[３]YangM G,LiCY,ChenZQ．AnewsimplenonＧlinearhysteretic modelfor MR damperand verification ofseiicresponsereductionexperiment[J]．EngineeringStructures,２０１３,５２(９):４３４Ｇ４４５．
[４]DominguezＧgonzalezA,STIHARUI,SEDAGHATIR．Practicalhysteresiodelformagnetorheologicaldampers[J]．JournalofIntelligentMaterialSystems&Structures,２０１４,２５(８):９６７Ｇ９７９．
[５]DomiinguezA,SedaghatiR,StiharuI．Modellingthehysteresisphenomenonof magnetorheologicaldampers[J]．SmartMaterials&Structures,２００４,１３(６):１３５１．
[６]ShiraziFA,MohammadpourJ,GrigoriadisK M,etal．Identificationandcontrolofanmrdamperwithstictioneffectanditsapplicationinstructuralvibrationmitigation[J]．IEEETransactionsonControlSystemsTechnology,２０１２,２０(５):１２８５Ｇ１３０１．
[７]Witters M,SweversJ．BlackＧbox modelidentificationforacontinuouslyvariable,electroＧhydraulicsemiＧactivedamper[J]．MechanicalSystems & SignalProcessing,２０１０,２４(１):４Ｇ１８．
[８]EkkachaiK,TungpimplrutK,NilkhamhangI．ForcecontrolofamagnetorheologicaldamperusinganelementaryhysteresiodelＧbasedfeedforwardneuralnetwork[J]．SmartMaterials&Structures,２０１３,２２(１１):５０３０．
[９]TruongDQ,AhnKK．NonlinearblackＧboxmodelsandforceＧsensorlessdampingcontrolfordampingsystemsusingmagnetoＧrheologicalfluiddampers[J]．Sensors& ActuatorsAPhysical,２０１１,１６７(２):５５６Ｇ５７３．
[１０]TruongDQ,AhnK K．IdentificationandapplicationofblackＧboxmodelforaselfＧsensingdampingsystemusingamagnetoＧrheologicalfluiddamper[J]．Sensors& ActuatorsAPhysical,２０１０,１６１(１/２):３０５Ｇ３２１．
[１１]ChangCC,RoschkeP．Neuralneworkmodelingofamagnetorheologicaldamper[J]．JournalofIntelligentMaterialSystems&Structures,１９９８,９(９):３５１Ｇ３５８．[１２]Wang D H,Liao W H． Modeling and controlof magnetorheologicalfluiddampersusingneuralnetworks[J]．SmartMaterials&Structures,２００４,１４(１):１１１Ｇ１２６．
[１３]Liem D T,Truong D Q,Ahn K K． HysteresiodelingofmagnetoＧrheologicaldamperusingselfＧtuningLyapunovＧbasedfuzzyapproach[J]．InternationalJournalofPrecision Engineering & Manufacturing,２０１５,１６(１):３１Ｇ４１．
[１４]FriedlandB,GurfilP．Newsimplemodelofdynamichysteresis[J]．TransactionsoftheInstituteofMeasurement& Control,２００３,２５(３):２０２Ｇ２１６．
[１５]KhadhraouiT,KtataS,BenzartiF,etal．FeaturesselectionbasedonmodifiedPSOalgorithmfor２dfacerecognition[C]//InternationalConferenceonComputerGraphics,ImagingandVisualization．IEEE,２０１６:９９Ｇ１０４．
阻尼器论文参考资料：
本文结束语：上述文章是一篇关于阻尼器和动态迟滞单元和流变方面的阻尼器论文题目、论文提纲、阻尼器论文开题报告、文献综述、参考文献的相关大学硕士和本科毕业论文。
