CN111460586A - 一种同时辨识滑动副间隙的间隙值和间隙接触刚度的方法 - Google Patents
一种同时辨识滑动副间隙的间隙值和间隙接触刚度的方法 Download PDFInfo
- Publication number
- CN111460586A CN111460586A CN202010178822.4A CN202010178822A CN111460586A CN 111460586 A CN111460586 A CN 111460586A CN 202010178822 A CN202010178822 A CN 202010178822A CN 111460586 A CN111460586 A CN 111460586A
- Authority
- CN
- China
- Prior art keywords
- gap
- nonlinear
- sliding pair
- force
- clearance
- 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 38
- 239000011159 matrix material Substances 0.000 claims abstract description 37
- 238000006073 displacement reaction Methods 0.000 claims abstract description 21
- 230000004044 response Effects 0.000 claims abstract description 21
- 238000004364 calculation method Methods 0.000 claims abstract description 19
- 238000005457 optimization Methods 0.000 claims abstract description 18
- 238000012546 transfer Methods 0.000 claims abstract description 18
- 230000005284 excitation Effects 0.000 claims abstract description 17
- 238000012360 testing method Methods 0.000 claims abstract description 9
- 230000005540 biological transmission Effects 0.000 claims abstract description 6
- 238000013016 damping Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 description 3
- 239000000470 constituent Substances 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000005299 abrasion Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002301 combined effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005461 lubrication Methods 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种同时辨识滑动副间隙的间隙值和间隙接触刚度的方法,步骤:在有限元软件中构建去除滑动副间隙后的线性结构,通过模态分析得到基础线性结构的特性矩阵;形成含滑动副的非线性结构扩展力的传递矩阵;对含滑动副间隙的非线性机械结构进行振动测试,测量该结构激励位置、滑动副间隙位置和参考位置处的位移响应信号以及激励力信号;根据上述信号,构建非线性结构时域非线性传递方程,得到结构中包含的滑动副间隙非线性力;根据滑动副的间隙非线性力的分段非线性表达式结合非线性最小二乘算法,对构建的目标函数进行优化计算,输出结果为间隙值和间隙接触刚度。本发明通过目标函数的优化计算实现了滑动副间隙值和间隙接触刚度同时辨识。
Description
技术领域
本发明涉及一种非线性系统参数辨识领域,特别是涉及一种同时辨识滑动副间隙的间隙值和间隙接触刚度的方法。
背景技术
对于含滑动副的机械装备,由于零部件装配、润滑以及相对滑动的需要,滑动副结合面之间不可避免地存在间隙。一方面,由于加工误差以及装配误差等因素使得滑动副间隙参数具有不确定性。另一方面,由于机械装备服役过程中的磨损、热效应以及局部变形等因素又使得间隙参数具有可变性。/在这些因素的综合作用下,间隙参数往往会偏离设计的合理范围。滑动副间隙过小,会引起滑动机构卡死,从而导致机械装备工作失效。通过合适的方法辨识出含滑动副间隙的非线性结构的间隙值和间隙接触刚度,有助于了解该非线性机械结构的动态特性,并为控制或减弱滑动副间隙对机械装备性能的不良影响奠定基础。
目前,针对滑动副间隙的非线性结构,现有的辨识间隙值和间隙接触的刚度的方法为两步间隙参数辨识法,主要是首先通过改进的恢复力曲面法实现间隙值的识别,然后将其识别结果用以计算间隙非线性描述函数,并由此利用经典的非线性系数识别方法开展间隙接触刚度的识别。但是,这个方法需要将间隙值识别结果代入到间隙接触刚度的识别中,即间隙值的识别误差会降低间隙接触刚度的识别精度。
因此,现在亟需一种针对滑动副间隙非线性结构,同时辨识滑动副间隙的间隙值和间隙接触刚度的方法。
发明内容
本发明的一个目的是解决至少上述问题和缺陷,并提供至少后面将说明的优点。
为了实现根据本发明的这些目的和其它优点,提供了一种同时辨识滑动副间隙的间隙值和间隙接触刚度的方法,其特征在于,包括以下步骤:
步骤一:在有限元软件中构建去除滑动副间隙后的基础线性结构,通过计算模态分析得到所述基础线性结构的特性矩阵;
步骤二:根据所述基础线性结构的特性矩阵,计算得到含滑动副间隙的非线性结构的扩展力传递矩阵;
步骤三:对所述含滑动副间隙的非线性结构进行振动测试,测量该结构激励位置、滑动副间隙位置和参考位置处的位移响应信号和激励力信号;
步骤四:根据所述测量得到的激励力信号和位移响应信号,采用正则化方法求解构建所述非线性结构的时域非线性传递方程,计算得到所述非线性结构的滑动副间隙非线性力;
步骤五:根据滑动副间隙非线性力的分段表达式与所述步骤四中计算得到的滑动副间隙非线性力结合非线性最小二乘算法进行优化计算,优化输出结果即为所述非线性结构的滑动副间隙的间隙值和间隙接触刚度。
优选的是,在步骤一中,所述线性结构的特性矩阵为去除滑动副间隙后的基础线性结构的质量矩阵M、刚度矩阵K和阻尼矩阵Cα。
优选的是,在步骤二中,非线性结构的扩展力fex=[ffnl]T是指由激励力f以及间隙非线性力fnl组成的非线性结构的外力。
优选的是,在步骤三中,在步骤三中,所述振动测试中由激励力所引起的最大位移响应幅值大于间隙值。
优选的是,在步骤二中,所述含滑动副间隙的非线性结构扩展力传递矩阵为:
其中,Tex表示所述含滑动副的非线性结构扩展力的传递矩阵,N表示采集的位移响应的数据点数,Ad、Bd和Cd表示离散时间非线性状态空间模型的系数矩阵。
优选的是,在步骤五中,所述目标函数的优化计算为
Θ(dc,kc)=rT(dc,kc)r(dc,kc)
其中,表示利用正则化方法构建非线性结构的时域非线性传递方程计算得到的滑动副间隙非线性力,表示通过优化迭代计算中预测的迭代参数结合滑动副间隙非线性力的分段表达式估计的滑动副间隙非线性力。 r(dc,kc)表示计算得到的滑动副间隙非线性力与估计的滑动副间隙非线性力的残差。
本发明至少包括以下有益效果:通过滑动副间隙非线性机械结构的传递方程结合正则化方法计算出结构中包含的滑动副间隙非线性力,还通过优化迭代计算中预测的迭代参数结合滑动副间隙非线性力表达式估计滑动副间隙非线性力,由此构建优化目标函数,并利用非线性最小二乘算法进行优化求解,从而实现滑动副间隙值和间隙接触刚度的辨识,实现了滑动副间隙值和间隙接触刚度的同时辨识,避免了传统两步法中间隙值识别误差对间隙接触刚度识别结果的影响;且这种计算方法不受滑动副间隙非线性机械结构的几何构型、间隙数目以及间隙位置等的限制。
本发明的其它优点、目标和特征将部分通过下面的说明体现,部分还将通过对本发明的研究和实践而为本领域的技术人员所理解。
附图说明
图1为本发明的流程图;
图2为本发明一个实施例中含双间隙悬臂梁第一个间隙的间隙非线性参数辨识结果曲线图;
图3为本发明一个实施例中含双间隙悬臂梁第二个间隙的间隙非线性参数辨识结果曲线图。
其中,曲线上有圆圈的为间隙接触刚度—迭代步数曲线,曲线上有叉号的为间隙值—迭代步数曲线。
具体实施方式
下面结合附图对本发明做进一步的详细说明,以令本领域技术人员参照说明书文字能够据以实施。
应当理解,本文所使用的诸如“具有”、“包含”以及“包括”术语并不配出一个或多个其它元件或其组合的存在或添加。
图1示出了根据本发明的一种实现形式,其中包括以下步骤:
步骤一:在有限元软件中构建去除滑动副间隙后的基础线性结构,通过计算模态分析得到所述基础线性结构的特性矩阵;
步骤二:根据所述基础线性结构的特性矩阵,计算得到含滑动副间隙的非线性结构的扩展力传递矩阵;
步骤三:对所述含滑动副间隙的非线性结构进行振动测试,测量该结构激励位置、滑动副间隙位置和参考位置处的位移响应信号和激励力信号;
步骤四:根据所述测量得到的激励力信号和位移响应信号,采用正则化方法求解构建所述非线性结构的时域非线性传递方程,计算得到所述非线性结构的滑动副间隙非线性力;
步骤五:根据滑动副间隙非线性力的分段表达式与所述步骤四中计算得到的滑动副间隙非线性力结合非线性最小二乘算法进行优化计算,优化输出结果即为所述非线性结构的滑动副间隙的间隙值和间隙接触刚度。
在这种技术方案中,首先在在有限元软件中构建去除滑动副间隙后的线性结构,通过模态分析得到对应基础线性结构的特性矩阵,并由此组装含滑动副的非线性结构扩展力的传递矩阵。对含滑动副间隙的非线性机械结构进行振动测试,测量该结构激励位置、滑动副间隙位置和参考位置处的位移响应信号以及激励力信号,进而利用正则化方法求解构建的非线性结构时域非线性传递方程,得到结构包含的滑动副间隙非线性力。根据滑动副间隙的非线性力的分段非线性表达式结合非线性最小二乘算法,对构建的目标函数进行优化计算,从而实现滑动副间隙值和间隙接触刚度的同时辨识。
上述方案中的所述线性结构的特性矩阵的一种实现方式为:所述线性结构的特性矩阵为去除滑动副间隙后的基础线性结构的质量矩阵M、刚度矩阵 K和阻尼矩阵Cα。
上述方案中的所述非线性结构扩展力为:所述非线性结构的扩展力 fex=[ffnl]T是由激励力f以及间隙非线性力fnl组成的非线性结构的外力。
上述方案中的所述非线性机械结构的振动测试的一种实现方式为:所述非线性机械结构的振动测试需要以较大的激励力激起机械结构中包含的间隙非线性,即滑动副间隙位置的最大位移响应幅值需大于间隙值。
上述方案中所述含滑动副的非线性结构扩展力的传递矩阵的一种实现方式为,所述含滑动副非线性结构扩展力的传递矩阵Tex可以表示为:
其中,N表示采集的位移响应的数据点数,Ad、Bd和Cd表示离散时间非线性状态空间模型的系数矩阵,其与连续时间非线性状态空间模型的系数矩阵Ac、Bc和Cc的关系可以表示为
其中,ΔTs表示数据采集的时间间隔,连续时间非线性状态空间模型的系数矩阵Ac、Bc和Cc的具体表达式为
上述方案中的所述滑动副的间隙非线性力的分段非线性表达式为:
其中,x(lc,t)表示滑动副间隙值位置处测得的位移响应,dc表示滑动副间隙值,kc表示滑动副间隙接触刚度。
上述方案中的所述目标函数的优化计算为:
Θ(dc,kc)=rT(dc,kc)r(dc,kc)
其中,表示利用正则化方法对非线性结构时域非线性传递方程求解得到的滑动副间隙非线性力,表示通过优化迭代计算中预测的迭代参数结合滑动副间隙非线性力的表达式估计的滑动副间隙非线性力。r(dc,kc)表示计算的滑动副间隙非线性力与估计的滑动副间隙非线性力的残差。
实施例
在激励力为幅值20N、频率10Hz的正弦载荷,采样频率为2048Hz的条件下辨识含双间隙的悬臂梁第一个间隙和第二个间隙的间隙值和间隙接触刚度。
根据图2的曲线图,第一个间隙位置处的实验结果如表1所示。
表1
真实值 | 实验值 | |
间隙值/m | 0.00005 | 4.994×10<sup>-5</sup> |
间隙接触刚度/N/m | 5×10<sup>5</sup> | 4.9694×10<sup>5</sup> |
根据表1所示的第一个间隙位置处的真实值与实验值可发现,间隙值的识别误差为0.12%。间隙接触刚度的识别误差为0.61%,两者的识别误差均很小,这种同时辨识滑动副间隙的间隙值和间隙接触刚度的方法较为精确。
根据图3的曲线图,第二个间隙位置处的实验结果如表2所示。
表2
真实值 | 实验值 | |
间隙值/m | 0.0001 | 1.008×10<sup>-4</sup> |
间隙接触刚度/N/m | 1×10<sup>6</sup> | 1.006×10<sup>6</sup> |
根据表1所示的第二个间隙位置处的真实值与实验值可发现,间隙值的识别误差为0.80%。间隙接触刚度的识别误差为0.60%,两者的识别误差均很小,不超过1%,则这种同时辨识滑动副间隙的间隙值和间隙接触刚度的方法很精确。且通过两个间隙的识别误差可知,这种计算方法不受间隙数目的限制。
这里说明的设备数量和处理规模是用来简化本发明的说明的。对本发明的应用、修改和变化对本领域的技术人员来说是显而易见的。
如上所述,根据本发明,本发明至少包括以下有益效果:通过滑动副间隙非线性机械结构的传递方程结合正则化方法计算出结构中包含的滑动副间隙非线性力,还通过优化迭代计算中预测的迭代参数结合滑动副间隙非线性力表达式估计滑动副间隙非线性力,由此构建优化目标函数,并利用非线性最小二乘算法进行优化求解,从而实现滑动副间隙值和间隙接触刚度的辨识,实现了滑动副间隙值和间隙接触刚度的同时辨识,避免了传统两步法中间隙值识别误差对间隙接触刚度识别结果的影响;且这种计算方法不受滑动副间隙非线性机械结构的几何构型、间隙数目以及间隙位置等的限制。
尽管本发明的实施方案已公开如上,但其并不仅仅限于说明书和实施方式中所列运用。它完全可以被适用于各种适合本发明的领域。对于熟悉本领域的人员而言,可容易地实现另外的修改。因此在不背离权利要求及等同范围所限定的一般概念下,本发明并不限于特定的细节和这里示出与描述的图例。
Claims (8)
1.一种同时辨识滑动副间隙的间隙值和间隙接触刚度的方法,其特征在于,包括以下步骤:
步骤一:在有限元软件中构建去除滑动副间隙后的基础线性结构,通过计算模态分析得到所述基础线性结构的特性矩阵;
步骤二:根据所述基础线性结构的特性矩阵,计算得到含滑动副间隙的非线性结构的扩展力传递矩阵;
步骤三:对所述含滑动副间隙的非线性结构进行振动测试,测量该结构激励位置、滑动副间隙位置和参考位置处的位移响应信号和激励力信号;
步骤四:根据所述测量得到的激励力信号和位移响应信号,采用正则化方法求解构建所述非线性结构的时域非线性传递方程,计算得到所述非线性结构的滑动副间隙非线性力;
步骤五:根据滑动副间隙非线性力的分段表达式与所述步骤四中计算得到的滑动副间隙非线性力结合非线性最小二乘算法进行优化计算,优化输出结果即为所述非线性结构的滑动副间隙的间隙值和间隙接触刚度。
2.根据权利要求1所述的同时辨识滑动副间隙的间隙值和间隙接触刚度的方法,其特征在于,在步骤一中,所述基础线性结构的特性矩阵为去除滑动副间隙后的基础线性结构的质量矩阵、刚度矩阵和阻尼矩阵。
3.根据权利要求1所述的同时辨识滑动副间隙的间隙值和间隙接触刚度的方法,其特征在于,在步骤二中,所述非线性结构的扩展力是指由激励力以及间隙非线性力组成的非线性结构的外力。
4.根据权利要求1所述的同时辨识滑动副间隙的间隙值和间隙接触刚度的方法,其特征在于,在步骤三中,所述振动测试中由激励力所引起的最大位移响应幅值大于间隙值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010178822.4A CN111460586B (zh) | 2020-03-13 | 2020-03-13 | 一种同时辨识滑动副间隙的间隙值和间隙接触刚度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010178822.4A CN111460586B (zh) | 2020-03-13 | 2020-03-13 | 一种同时辨识滑动副间隙的间隙值和间隙接触刚度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111460586A true CN111460586A (zh) | 2020-07-28 |
CN111460586B CN111460586B (zh) | 2023-05-12 |
Family
ID=71678444
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010178822.4A Active CN111460586B (zh) | 2020-03-13 | 2020-03-13 | 一种同时辨识滑动副间隙的间隙值和间隙接触刚度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111460586B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106773678A (zh) * | 2016-11-30 | 2017-05-31 | 西安交通大学 | 用于多非线性参数耦合系统的参数辨识方法及其辨识设备 |
CN107480097A (zh) * | 2017-06-27 | 2017-12-15 | 西安交通大学 | 非线性机械结构的载荷识别方法及载荷识别设备 |
CN107563005A (zh) * | 2017-08-04 | 2018-01-09 | 大连理工大学 | 一种拉压不同刚度结构振动瞬时最优控制方法 |
-
2020
- 2020-03-13 CN CN202010178822.4A patent/CN111460586B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106773678A (zh) * | 2016-11-30 | 2017-05-31 | 西安交通大学 | 用于多非线性参数耦合系统的参数辨识方法及其辨识设备 |
CN107480097A (zh) * | 2017-06-27 | 2017-12-15 | 西安交通大学 | 非线性机械结构的载荷识别方法及载荷识别设备 |
CN107563005A (zh) * | 2017-08-04 | 2018-01-09 | 大连理工大学 | 一种拉压不同刚度结构振动瞬时最优控制方法 |
Non-Patent Citations (3)
Title |
---|
刘杰;李兵;金炜;韩罗峰;权双璐;: "用于连续体的间隙非线性参数辨识方法" * |
彭利坤;吕帮俊;熊先锋;邢继峰;: "液压Stewart平台谐振频率优化设计研究" * |
王博;马志赛;丁千;张欣;杨宁;: "基础激励下含间隙折叠舵面非线性系统辨识" * |
Also Published As
Publication number | Publication date |
---|---|
CN111460586B (zh) | 2023-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cao et al. | Damage identification for beams in noisy conditions based on Teager energy operator-wavelet transform modal curvature | |
Wang et al. | Identification of piezoelectric hysteresis by a novel Duhem model based neural network | |
CN101821582B (zh) | 用于确定测量机中的动态误差的方法 | |
Maes et al. | A smoothing algorithm for joint input-state estimation in structural dynamics | |
KR101634176B1 (ko) | 재료의 푸아송비 및 잔류응력의 측정 방법 | |
CN109902408B (zh) | 一种基于数值运算和改进的正则化算法的载荷识别方法 | |
CN112580236B (zh) | 一种热防护连接结构非线性动态响应的快速分析方法 | |
CN107480097B (zh) | 非线性机械结构的载荷识别方法及载荷识别设备 | |
Pelayo et al. | Modal analysis based stress estimation for structural elements subjected to operational dynamic loadings | |
Greek et al. | Deflection of surface-micromachined devices due to internal, homogeneous or gradient stresses | |
Yan et al. | A new probabilistic frequency-domain approach for influence line extraction from static transmissibility measurements under unknown moving loads | |
CN112182697A (zh) | 一种有阻尼吊杆系统张力的高精度动测法 | |
CN115455353A (zh) | 一种基于非线性时域滤波的在线参数辨识方法 | |
Li et al. | Principal force pattern and impulse response mode for structural equivalent force estimation and full-field response reconstruction | |
KR20090005769A (ko) | 등가정하중을 이용한 동적 비선형 응답 구조 최적해산출방법 | |
CN111460586A (zh) | 一种同时辨识滑动副间隙的间隙值和间隙接触刚度的方法 | |
Ming et al. | Free vibrations of elastic circular toroidal shells | |
CN106679911B (zh) | 一种基于多尺度数据融合理论的梁型结构损伤识别方法 | |
Haddab et al. | Improvement of strain gauges micro-forces measurement using Kalman optimal filtering | |
CN111428413B (zh) | 连续梁分段抗弯刚度识别方法 | |
He et al. | A general approach to construct full-field responses and high-resolution mode shapes of bridges using sparse measurements | |
CN113094640A (zh) | 一种频域下宽带多轴随机振动寿命预测方法 | |
CN117238410A (zh) | 一种基于物理信息网络的材料参数未知薄板挠度估计方法 | |
CN107220678B (zh) | 多自由度梁式结构非线性类型确定方法 | |
Onchis | A deep learning approach to condition monitoring of cantilever beams via time-frequency extended signatures |
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 |