磁流变减振器基于matlab的仿真分析

《磁流变减振器基于matlab的仿真分析》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、黑龙江工程学院本科生毕业设计磁流变减振器基于Matlab的仿真分析摘要：基于磁流变减振器在汽车悬架减振系统半主动控制中的广泛应用，根据磁流变液的特点和磁流变减振器阻尼力与结构参数的关系，设计了新型的磁流变减振器，并对影响磁流变减振器性能的参数进行了仿真。仿真表明，该磁流变减振器设计计算是一种能优化阻尼力的有效算法。关键词：磁流变减振器；半主动控制；磁流变液1.1减振器的阻尼力计算模型本文选用剪切阀式磁流变阻尼器工作模式进行结构设计，在结构设计前，必须明确该工作模式磁流变液的流变方程，继而推导出磁流变阻尼力的计算模型，这是结构设

2、计过程中的依据所在。基于剪切阀式磁流变阻尼器的阻尼通道的宽度远大于其阻尼间隙，因而可简化成磁流变液在两相对运动平板之间的运动。为了简化分析，工作于剪切阀式的磁流变阻尼力可以看成是在阀式工作模式下的阻尼力和剪切工作模式下阻尼力的叠加。在外加磁场作用下，磁流变液表现Bingham流体，其磁流变液在平板的流动和速度分布如图1.1所示，其本构关系可用下列方程描述：（1.1）（1.2）8黑龙江工程学院本科生毕业设计图1.1磁流变液在平板中的流动和速度分布在阀式工作模式下磁流变液的速度分布如图1.1所示。假设磁流变液的体积流速Q在x方向上

3、一维流动，在y方向上不流动。设两平板之间的间隙为h，长度为L，宽度为b，由流体力学可得下列微分方程：（1.3）式中u、v分别是磁流变液在x、y方向上的流动速度；是磁流变液在x方向的压力梯度，为了简化将压力梯度是为x方向线性变化=，l是阻尼通道的长度；是阻尼通道两端的压力差；是磁流变液的密度；t是时间变量；由于流动速度低，可不计惯性效应，；令沿x的剪切应力，由于磁流变流动的连续性，沿x方向的速度不变即则方程（1.3）简化为：（1.4）对其积分可得：（1.5）D是待定的积分常数。由公式（1.4）可知，磁流变液受到的剪切应力沿平板间

4、隙是按线性分布的，靠近平板的磁流变液受到的剪切力最大，而中间对称面上的磁流变液受到的剪切应力最小，根据极板两端压差产生的剪切应力与极板附近磁流变液的临界剪切屈服应力比较，当前者小于后者磁流变液静止不动；当前者大于后者将产生如图1.1所示的流体状态，即靠近平板处得磁流变液流动；而中间对称区间的磁流变液不流动。可将此时的磁流变液的流动分为屈服流动，刚性流动，屈服流动三个区域。区域：屈服流动剪切应变率，由公式（1.1）可得：（1.6）8黑龙江工程学院本科生毕业设计将公式（1-6）代入公式（1.5）中，并注意u（0）=0，求解微分方程

5、如下：（1.7）（1.8）（1.9）区域：刚性流动，剪切应变率，同理可得：（1.10）区域：屈服流动，剪切应变率（1.11）将公式（1.11）代入公式（1.5），已知u(h)=0，，求解微分方程得：（1.12）（1.13）由公式和（1.8）公式（1.13）相减可得刚性流动区得厚度为（1.14）由于存在，由公式（1.9）和公式（1.13）可得（1.15）（1.16）8黑龙江工程学院本科生毕业设计由公式（1.14）和公式（1.16）可得：；（1.17）流经平板间隙的磁流变液的体积流量Q可有下列得到：（1.18）代入化简可得（1.1

6、9）经进一步化简可得压差近似公式：（1.20）考虑到阻尼器的实际阻尼通道为环形通道，流动模式下的阻尼力可以表示为：(1.21)式中为活塞受压的有效面积。在移动平板的影响下，磁流变液发生屈服流动，剪切模式下磁流变液的速度分布如图1.2所示。剪切应变率，则由公式（1.1），剪切应力可表示为：(1.22)假如磁流变液的速度是沿y方向分布如图1.2所示，即图1.2剪切模式下磁流变液的速度分布8黑龙江工程学院本科生毕业设计剪切模式下的阻尼力：(1.23)混合工作模式的阻尼力可视为流动模式、剪切模式两种工作模式下的阻尼力的叠加。即，由于符

7、号的正负只反映活塞运动的方向，因此，整理上式得：(1.24)式中参数c变化范围2-3，本文c=2，因此剪切阀式磁流变阻尼器阻尼力为：(1.25)；公式可以改为：(1.26)(1.27)(1.28)式中粘滞阻尼力系数：；库伦阻尼力：；为磁流变阻尼器活塞运动速度；sgn为符号函数；为30-50K从上式可以看出磁流变阻尼器的阻尼力由两部分组成，一部分由液体流动时液体粘性产生的粘滞阻尼力，而另一部分由磁流变效应产生的库伦阻尼力组成。8黑龙江工程学院本科生毕业设计1.2磁流变减振器的仿真分析磁流变减振器的数学模型采用公式1.25，建立磁

8、流变减振器的仿真模型如图1.3所示。图1.3仿真模型图1.4Matlab仿真图8黑龙江工程学院本科生毕业设计由公式（1.25）作为数学模型可进行计算。F=（）+（）=1504.5+2108上式的计算结果是在阻尼间隙为0.6mm是计算而得。在不同的速度下可计算出不同的磁流变阻尼