基于高阶非线性摩擦的金属橡胶材料迟滞动力学建模方法技术

本发明专利技术涉及一种基于高阶非线性摩擦的金属橡胶材料迟滞动力学建模方法，包括以下步骤：根据金属橡胶内部金属丝的接触形式以及空间角度分布构建高阶非线性摩檫力模型；构建传统动力学模型；将高阶非线性摩擦力模型代入传统动力学模型中，构建基于高阶非线性摩擦的金属橡胶迟滞动力学模型；对金属橡胶迟滞动力学模型进行参数识别；构建误差模型来表示动力学模型外影响试验结果的误差；对误差模型进行参数识别，确定应纳入模型的最佳模型项数量，完成误差模型的确定；完成误差模型参数优选，根据基于误差模型的高阶非线性金属橡胶迟滞动力学模型的模型相关系数完成模型构建。该方法有利于建立具有高精度和通用性的金属橡胶动力学模型。力学模型。力学模型。

【技术实现步骤摘要】
基于高阶非线性摩擦的金属橡胶材料迟滞动力学建模方法


[0001]本专利技术属于金属橡胶领域，具体涉及一种基于高阶非线性摩擦的金属橡胶材料迟滞动力学建模方法。

技术介绍

[0002]金属橡胶是一种通过金属螺旋卷细丝缠绕、或编织、或铺设等毛坯工艺，经过冷冲压制备而成，其内部结构形成无序、网状、多孔的高弹性纯金属减振材料，在特殊服役环境下已经得到了一定的应用。但是由于其结构的复杂导致金属橡胶具有较强的非线性迟滞特性，而无法精确构建金属橡胶的动力学模型，最终限制了其广泛的应用。当前学者尝试将金属橡胶内部结构的金属丝等效成微元模型，构建金属橡胶材料的本构模型，但这些模型忽略了金属橡胶的阻尼特性，导致其模型精度不足。为了解决精度问题，学者们采用唯象理论，又引入阻尼力对金属橡胶进行动力学建模，但模型仍然无法从内部金属丝空间随机接触形式研究其阻尼力的变化，因此所建模型无通用性，亟需通过从金属橡胶内部金属丝间的接触摩擦特性，构建一个全新的具有高精度和通用性的动力学模型。

技术实现思路

[0003]本专利技术的目的在于提供一种基于高阶非线性摩擦的金属橡胶材料迟滞动力学建模方法，该方法有利于建立具有高精度和通用性的金属橡胶动力学模型。
[0004]为实现上述目的，本专利技术采用的技术方案是：一种基于高阶非线性摩擦的金属橡胶材料迟滞动力学建模方法，包括以下步骤：
[0005]步骤S1：根据金属橡胶内部金属丝的接触形式以及空间角度分布构建高阶非线性摩檫力模型；
[0006]步骤S2：依据传统动力学模型构建非线性弹性恢复力模型、非线性阻尼力模型以及迟滞阻尼力模型，采用切比雪夫多项式将迟滞阻尼力分解；
[0007]步骤S3：将高阶非线性摩擦力模型代入传统动力学模型中，构建基于高阶非线性摩擦的金属橡胶迟滞动力学模型；
[0008]步骤S4：采用基于非线性边界条件约束的内部点算法的Lsqlin求解器对金属橡胶迟滞动力学模型进行参数识别；
[0009]步骤S5：将识别的参数回代到金属橡胶迟滞动力学模型，构建误差模型NARMAX来表示动力学模型外影响试验结果的误差；
[0010]步骤S6：采用OLS算法对误差模型进行参数识别，并用APRESS统计量确定应纳入模型的最佳模型项数量，当识别出来的误差模型APRESS>0.95，则完成误差模型的确定，否则返回步骤S5，修正误差模型项数量；
[0011]步骤S7：完成误差模型参数优选，计算基于误差模型的高阶非线性金属橡胶迟滞动力学模型的模型精度指标R2，当R2>0.9完成高阶非线性金属橡胶迟滞动力学模型，否则返回步骤S2，重新构建传统动力学模型。
[0012]进一步地，步骤S1中，所述高阶非线性摩檫力模型的公式如下：
[0013][0014]式中，f
fH
为高阶非线性摩檫力，ε为金属橡胶形变位移，F
N
为金属丝受到的等效弹性力，μ为金属丝摩擦系数，k
′
m
′
为相邻金属丝点接触形式等效刚度，k
′
m
″
为相邻金属丝面接触等效刚度，d为概率分布多项式系数。
[0015]进一步地，步骤S2中，依据传统动力学模型构建非线性弹性恢复力模型、非线性阻尼力模型以及迟滞阻尼力模型，采用切比雪夫多项式将迟滞阻尼力分解，其公式如下：
[0016][0017]式中，f
k
(y(t))为非线性弹性恢复力，y(t)为金属橡胶形变位移ε，t为时间，为非线性阻尼力模型，z(t)为迟滞阻尼力，n1为非线性弹性恢复力的多项式项的数量，k
2i
‑1为刚度系数，n2为非线性阻尼力的多项式的项数量，c
2i
‑1为阻尼系数，a0与a
n
为切比雪夫多项式系数，n3为切比雪夫多项的项数量，为位移的导数即速度，Δy为金属橡胶最大变形位移。
[0018]进一步地，步骤S3中，所述基于高阶非线性摩擦的金属橡胶迟滞动力学模型的公式如下：
[0019][0020]式中，F(t)为金属橡胶受到的外力合。
[0021]进一步地，步骤S4中，对高阶非线性摩擦的金属橡胶迟滞动力学模型进行参数识别，并将识别的参数记作：
[0022][0023]式中，Θ为待识别参数，[θ1,θ2,...,θ
n
]为待识别参数组成的向量。
[0024]进一步地，所述误差模型NARMAX的公式如下：
[0025][0026]式中，Δ为误差值，ξ
i
为模型的残差项，L
ε
为残差的参数项数目，N
ξ
为误差模型的项数目，为误差模型的系数。
[0027]进一步地，步骤S1中，金属丝的接触形式表现为：
[0028][0029]式中，P
k
′
为金属橡胶中相邻金属丝无接触的概率，为金属橡胶中相邻金属丝点接触的概率，P
k
″′
为金属橡胶中相邻金属丝面接触的概率；
[0030]空间角度分布概率为：
[0031][0032]式中，为空间分布角度，randsrc为空间角度随机分布概率函数，n为金属橡胶中微弹簧单元的数目，为离散空间角度分布。
[0033]进一步地，步骤S6中，所述APRESS统计量为：
[0034]APRESS(n)＝c(n)MSE(n)
ꢀꢀꢀ
(8)
[0035]式中，MSE(q)＝1
‑
ERR
j
为误差输出方差，c(q)复杂性成本函数，α为大于1的数，N为最大模型的项数目，q为达到符合要求仅需的项数目。
[0036]进一步地，步骤S7中，所述模型相关系数R2为：
[0037][0038]式中，RSS是残差的平方之和；TSS指总平方之和；n是样本观测值的数量；y
i
观测值；表示样本是模型的预测值；代表样本观测值的平均值。
[0039]与现有技术相比，本专利技术具有以下有益效果：相较于现有的金属橡胶非线性迟滞动力学建模，本方法根据金属橡胶内部金属丝空间倾斜角与接触形式，首次提出高阶非线性摩擦力，并结合非线性弹性恢复力、非线性阻尼力以及库伦迟滞摩檫力，构建出一种精度更高的金属橡胶迟滞动力学模型，改善了传统金属橡胶动力学模型通用性不足的问题，优化了等效建模方式模型精度。
附图说明
[0040]图1是本专利技术实施例的方法实现流程图。
[0041]图2是本专利技术实施例中三种密度的环形金属橡胶样件。
[0042]图3是本专利技术实施例中测试试验工装图。
[0043]图4是本专利技术实施例中试验结果与模型曲线对比图。
[0044]图5是本专利技术实施例中误差模型时序图。
[0045]图6是本专利技术实施例中理论模型时域图与试验结果时域图的对比图。
[0046]图7是本专利技术实施例中金属本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于高阶非线性摩擦的金属橡胶材料迟滞动力学建模方法，其特征在于，包括以下步骤：步骤S1：根据金属橡胶内部金属丝的接触形式以及空间角度分布构建高阶非线性摩檫力模型；步骤S2：依据传统动力学模型构建非线性弹性恢复力模型、非线性阻尼力模型以及迟滞阻尼力模型，采用切比雪夫多项式将迟滞阻尼力分解；步骤S3：将高阶非线性摩擦力模型代入传统动力学模型中，构建基于高阶非线性摩擦的金属橡胶迟滞动力学模型；步骤S4：采用基于非线性边界条件约束的内部点算法的Lsqlin求解器对金属橡胶迟滞动力学模型进行参数识别；步骤S5：将识别的参数回代到金属橡胶迟滞动力学模型，构建误差模型NARMAX来表示动力学模型外影响试验结果的误差；步骤S6：采用OLS算法对误差模型进行参数识别，并用APRESS统计量确定应纳入模型的最佳模型项数量，当识别出来的误差模型APRESS>0.95，则完成误差模型的确定，否则返回步骤S5，修正误差模型项数量；步骤S7：完成误差模型参数优选，计算基于误差模型的高阶非线性金属橡胶迟滞动力学模型的模型精度指标R2，当R2>0.9完成高阶非线性金属橡胶迟滞动力学模型，否则返回步骤S2，重新构建传统动力学模型。2.根据权利要求1所述的基于高阶非线性摩擦的金属橡胶材料迟滞动力学建模方法，其特征在于，步骤S1中，所述高阶非线性摩檫力模型的公式如下：式中，f
fH
为高阶非线性摩檫力，ε为金属橡胶形变位移，F
N
为金属丝受到的等效弹性力，μ为金属丝摩擦系数，k
″
m
为相邻金属丝点接触形式等效刚度，k
″′
m
为相邻金属丝面接触等效刚度，d为概率分布多项式系数。3.根据权利要求2所述的基于高阶非线性摩擦的金属橡胶材料迟滞动力学建模方法，其特征在于，步骤S2中，依据传统动力学模型构建非线性弹性恢复力模型、非线性阻尼力模型以及迟滞阻尼力模型，采用切比雪夫多项式将迟滞阻尼力分解，其公式如下：型以及迟滞阻尼力模型，采用切比雪夫多项式将迟滞阻尼力分解，其公式如下：型以及迟滞阻尼力模型，采用切比雪夫多项式将迟滞阻尼力分解，其公式如下：型以及迟滞阻尼力模型，采用切比雪夫多项式将迟滞阻尼力分解，其公式如下：式中，f
k
(y(t))为非线性弹性恢复力，y(t)为金属橡胶形变位移ε，t为时间，为非
线性阻尼力模型，z(t)为迟滞阻尼力，n1为非线性弹性恢复力的多项式项的数量，k
2i
‑1为刚度系数，n2为非线性阻尼力的多项式的项数量，c
2i
‑1为阻尼系数，a0与a
n
为切比雪夫多项式系数，n3为切比雪夫多项的项数量...

【专利技术属性】
技术研发人员：任志英，周春辉，林有希，秦红玲，史林炜，马国政，
申请(专利权)人：福州大学，
类型：发明
国别省市：