磁流变阻尼器等效线性阻尼系数解析模型建模方法技术

本发明专利技术公开了一种磁流变阻尼器等效线性阻尼系数解析模型建模方法，其特征是：通过MTS力学测量系统实验获取磁流变阻尼器的样本数据，利用磁流变阻尼器一个周期耗散能量计算不同工况的磁流变阻尼器等效现行阻尼系数，采用线性拟合和幂函数拟合构造磁流变阻尼器的三个独立工况参数分别与等效线性阻尼系数的关系，通过耦合计算得到耦合了三个独立工况参数的磁流变阻尼器等效线性阻尼系数的解析模型，采用最小二乘伪逆方法利用试验样本数据对等效线性阻尼系数解析模型的待辨识参数进行辨识。本发明专利技术可用于建立磁流变阻尼器等效动力学模型，可用于安装有磁流变阻尼器结构系统的稳定性分析，也可用于磁流变阻尼器的性能状态的预测和判定。

【技术实现步骤摘要】
磁流变阻尼器等效线性阻尼系数解析模型建模方法
本专利技术属于机械系统动力学建模及振动控制领域。具体涉及一种磁流变阻尼器等效线性阻尼系数解析模型建模方法。
技术介绍
磁流变阻尼减振器是一种新型的结构简单、响应速度快、阻尼力连续可调的半主动减振装置，经过近几十年的研究发展，其越来越广泛地应用在车辆工程、建筑以及桥梁等结构系统上。目前，磁流变阻尼减振器的力学模型可以分为两类：参数化模型和非参数化模型。参数化模型主要有：Bingham模型、非线性双粘性模型、修正的Bouc_Wen模型、修正的Dahl模型、现象模型等，这些基本上都是由结构复杂的非线性方程构成，虽能较好地描述磁流变阻尼器的性能，但求解困难，制约了磁流变阻尼器的应用。非参数化模型主要有：多项式模型、Sigmoid模型等半几何模型以及利用模糊理论、神经网络和黑箱理论等建立的智能化模型。多项式模型在激励条件改变后，需更新模型参数，不利于实时控制；Sigmoid模型待辨识参数较多，且不能较好地描述磁流变阻尼器的性能；而智能化模型理论复杂，并需要大量的训练数据，难以在线实现。亟待开展结构简单，宜于应用的磁流变阻尼器力学模型相关研究工作。由于磁流变阻尼减振器自身力学性能具有非线性，加上磁流变阻尼减振器使用工况的复杂性，磁流变阻尼器的阻尼力是时变的，其等效线性阻尼系数也是时变的。实时在线识别磁流变阻尼减振器的等效线性阻尼系数具有重要的意义。利用等效线性阻尼系数建立磁流变阻尼器等效动力学模型，可以简化控制系统，有利于系统控制的实时性；并且对分析安装有磁流变阻尼器结构系统的稳定性提供了便利；此外，利用磁流变阻尼器等效线性阻尼系数随时间在同等工况下的改变可以对磁流变阻尼器的寿命进行预测和判定。磁流变阻尼减振器的等效线性阻尼系数与三个参数，即输入电压、位移激励振幅以及激励频率等因素有关。目前，求解磁流变阻尼减振器等效线性阻尼主要用能量法。能量法的原理为：将阻尼器在给定电流以及给定的工况下循环运行若干周期后，根据实际消耗的能量以及等效线性阻尼器消耗的能量相等，计算该工况下的等效线性阻尼系数。这种方法的局限性在于，只能获得单一工况下的等效线性阻尼系数等。当工况（即输入电压、位移激励振幅以及激励频率）改变后，需要进行重新计算。能量法计算量大，繁琐，不适合实时在线获得磁流变阻尼减振器的等效线性阻尼系数。
技术实现思路
本专利技术是为避免上述现有技术所存在的不足之处，提供一种磁流变阻尼器等效线性阻尼系数解析模型建模方法，以避免能量法的局限性，为磁流变阻尼减振器建立表征其固有特征的等效线性阻尼系数解析模型。该模型耦合输入电压、位移激励振幅以及激励频率，能随着阻尼器工况的改变而实时、在线地调节其等效线性阻尼系数。进而可以简化应用磁流变阻尼减振器的结构系统动力学模型的构造，获得系统时变的动态特性，便于振动控制；同时能精确控制系统的等效线性阻尼；还能对磁流变阻尼减振器的性能状态进行监测。本专利技术为解决技术问题采用如下技术方案：本专利技术磁流变阻尼器等效线性阻尼系数解析模型建模方法的特点是按如下过程进行：步骤一：获取样本数据，应用MTS力学测量系统对需要建立等效线性阻尼系数解析模型的磁流变阻尼器进行力学性能实验；分别改变磁流变阻尼器的输入电压V、正弦位移激励振幅X以及激励频率Ω来获得不同的工况，其中，输入电压V有L个元素，正弦位移激励振幅X有M个元素，激励频率Ω有N个元素，共获得L×M×N组不同工况的实验数据；测量获取不同工况下，在时间段T内，磁流变阻尼器杆件的位移时间曲线x(t)和阻尼力时间曲线f(t)；步骤二：计算不同工况下等效线性阻尼系数Ceq，对不同工况下实验得到的位移时间曲线x(t)进行微分获得磁流变阻尼器杆件的运动速度时间曲线则由式(1)获得磁流变阻尼器在对应工况下的等效线性阻尼系数Ceq为：式(1)中，E为磁流变阻尼器往复运动一个周期所耗散的能量，并有步骤三：构造数学关系，对于所述L×M×N组的不同工况下的输入电压V、正弦位移激励振幅X和激励频率Ω三个参数，分别固定其中两个，变化第三个参数，得到等效线性阻尼系数Ceq分别单独与每个参数的关系曲线；利用所得关系曲线的形状，根据形状相似的原则，对等效线性阻尼系数Ceq和输入电压V采用式(2)所示线性拟合，对等效线性阻尼系数Ceq和正弦位移激励振幅X采用式(3)所示幂函数拟合，对等效线性阻尼系数Ceq和激励频率Ω采用式(4)所示幂函数拟合：Ceq=a1V+b1(2)Ceq=c1X-1+d1(3)Ceq=e1Ω-1+g1(4)式(2)、式(3)和式(4)中，ɑ1、b1、c1、d1、e1和g1为模型过渡参数；步骤四：构造解析模型，根据输入电压V、正弦位移激励振幅X以及激励频率Ω三个参数的独立性，将步骤三中的式(2)、式(3)和式(4)进行耦合计算，得到等效线性阻尼系数Ceq的解析模型如式(5)所示：Ceq=(a2V+b2)(h(XΩ)-1+j)(5)式(5)中：ɑ2、b2、h、j为模型过渡参数；记为：A=a2h；B=a2j；C=b2h；D=b2j，则有：式(6)中，A、B、C、D为待辨识参数；步骤五：对待辨识参数A、B、C、D进行参数辨识：应用MTS力学测量系统对需要建立等效线性阻尼系数解析模型的磁流变阻尼器进行力学性能试验，用正弦信号作为位移激励信号，改变输入电压V、正弦位移激励振幅X以及激励频率Ω这三个参数以获得不同工况下的试验样本数据，输入电压V有L1个元素，正弦位移激励振幅X有M1个元素，激励频率Ω有N1个元素，共获得L1×M1×N1组试验样本数据，利用所述实验样本数据采用最小二乘伪逆法实现待辨识参数A、B、C、D的参数辨识。与已有技术相比，本专利技术有益效果体现在：1、本专利技术磁流变阻尼器等效线性阻尼系数的解析模型能利用一个简单的数学模型把影响阻尼器等效线性阻尼系数的三个工况参数：输入电压V、位移激励振幅X以及激励频率Ω耦合在一起，结构简单，参数容易辨识，方便应用于构造结构系统的动力学模型。而已有的能量法，没有直接将工况参数耦合在一个数学公式，而是需要多次变化，比较繁琐复杂，不利于实时获取等效线性阻尼系数。2、利用本专利技术提出的等效线性阻尼系数解析模型，能实时获取系统时变的阻尼信息，结合系统的刚度和质量信息，从而能够得到系统时变的动态特性，便于振动控制。而已有的能量法只能获得特定工况下的等效线性阻尼系数，不具有实时辨识的功能。3、利用本专利技术提出的等效线性阻尼系数解析模型的逆模型，可以实现由输入电压到阻尼系数的精确控制。在激励频率和位移振幅相对稳定的情况下，调节输入电压，能比较精确地控制结构系统的等效线性阻尼系数，满足工程实际的需求。这在一些特殊场合下，比如移动缆绳的振动，绳系卫星绳索的振动，都需要精确地的阻尼来最大限度地消除绳上的振动。目前还没有通过磁流变阻尼器提供精确阻尼控制的解析模型。4、本专利技术模型能够应用于磁流变阻尼减振器固有力学性能变化的监测。通过辨识同一个阻尼器在不同时期的等效线性阻尼系数解析模型，在同样的工况下，检测和比较两个模型的等效线性阻尼系数，可以监测磁流变阻尼减振器固有力学性能变化。附图说明图1为本专利技术方法流程图。具体实施方式本实施例中磁流变阻尼器等效线性阻尼系数解析模型建模方法是按如下过程进行：步骤一：应用MTS力本文档来自技高网...

【技术保护点】
一种磁流变阻尼器等效线性阻尼系数解析模型建模方法，其特征是按如下过程进行：步骤一：获取样本数据，应用MTS力学测量系统对需要建立等效线性阻尼系数解析模型的磁流变阻尼器进行力学性能实验；分别改变磁流变阻尼器的输入电压V、正弦位移激励振幅X以及激励频率Ω来获得不同的工况，其中，输入电压V有L个元素，正弦位移激励振幅X有M个元素，激励频率Ω有N个元素，共获得L×M×N组不同工况的实验数据；测量获取不同工况下，在时间段T内，磁流变阻尼器杆件的位移时间曲线x(t)和阻尼力时间曲线f(t)；步骤二：计算不同工况下等效线性阻尼系数Ceq，对不同工况下实验得到的位移时间曲线x(t)进行微分获得磁流变阻尼器杆件的运动速度时间曲线则由式(1)获得磁流变阻尼器在对应工况下的等效线性阻尼系数Ceq为：Ceq=EπΩX2---(1)式(1)中，E为磁流变阻尼器往复运动一个周期所耗散的能量，并有步骤三：构造数学关系，对于所述L×M×N组的不同工况下的输入电压V、正弦位移激励振幅X和激励频率Ω三个参数，分别固定其中两个，变化第三个参数，得到等效线性阻尼系数Ceq分别单独与每个参数的关系曲线；利用所得关系曲线的形状，根据形状相似的原则，对等效线性阻尼系数Ceq和输入电压V采用式(2)所示线性拟合，对等效线性阻尼系数Ceq和正弦位移激励振幅X采用式(3)所示幂函数拟合，对等效线性阻尼系数Ceq和激励频率Ω采用式(4)所示幂函数拟合：Ceq=a1V+b1???????(2)Ceq=c1X?1+d1???????(3)Ceq=e1Ω?1+g1?????????(4)式(2)、式(3)和式(4)中，ɑ1、b1、c1、d1、e1和g1为模型过渡参数；步骤四：构造解析模型，根据输入电压V、正弦位移激励振幅X以及激励频率Ω三个参数的独立性，将步骤三中的式(2)、式(3)和式(4)进行耦合计算，得到等效线性阻尼系数Ceq的解析模型如式(5)所示：Ceq=(a2V+b2)(h(XΩ)?1+j)??????(5)式(5)中：ɑ2、b2、h、j为模型过渡参数；记为：A=a2h；B=a2j；C=b2h；D=b2j，则有：Ceq=AVXΩ+BV+C1XΩ+D---(6)式(6)中，A、B、C、D为待辨识参数；步骤五：对待辨识参数A、B、C、D进行参数辨识：应用MTS力学测量系统对需要建立等效线性阻尼系数解析模型的磁流变阻尼器进行力学性能试验，用正弦信号作为位移激励信号，改变输入电压V、正弦位移激励振幅X以及激励频率Ω这三个参数以获得不同工况下的试验样本数据，输入电压V有L1个元素，正弦位移激励振幅X有M1个元素，激励频率Ω有N1个元素，共获得L1×M1×N1组试验样本数据，利用所述实验样本数据采用最小二乘伪逆法实现待辨识参数A、B、C、D的参数辨识。FDA0000400470240000014.jpg,FDA0000400470240000012.jpg...

【技术特征摘要】
1.一种磁流变阻尼器等效线性阻尼系数解析模型建模方法，其特征是按如下过程进行：步骤一：获取样本数据，应用MTS力学测量系统对需要建立等效线性阻尼系数解析模型的磁流变阻尼器进行力学性能实验；分别改变磁流变阻尼器的输入电压V、正弦位移激励振幅X以及激励频率Ω来获得不同的工况，其中，输入电压V有L个元素，正弦位移激励振幅X有M个元素，激励频率Ω有N个元素，共获得L×M×N组不同工况的实验数据；测量获取不同工况下，在时间段T内，磁流变阻尼器杆件的位移时间曲线x(t)和阻尼力时间曲线f(t)；步骤二：计算不同工况下等效线性阻尼系数Ceq，对不同工况下实验得到的位移时间曲线x(t)进行微分获得磁流变阻尼器杆件的运动速度时间曲线则由式(1)获得磁流变阻尼器在对应工况下的等效线性阻尼系数Ceq为：式(1)中，E为磁流变阻尼器往复运动一个周期所耗散的能量，并有步骤三：构造数学关系，对于所述L×M×N组的不同工况下的输入电压V、正弦位移激励振幅X和激励频率Ω三个参数，分别固定其中两个，变化第三个参数，得到等效线性阻尼系数Ceq分别单独与每个参数的关系曲线；利用所得关系曲线的形状，根据形状相似的原则，对等效线性阻尼系数Ceq和输入电压V采用式(2)所示线性拟合，对等效线性阻尼系数Ceq和正弦位移激励振幅X采用式(3)所示幂...

【专利技术属性】
技术研发人员：陈恩伟，郭耀辉，夏增强，陆益民，刘正士，王勇，包猛，
申请(专利权)人：合肥工业大学，
类型：发明
国别省市：