CN103970066B - 基于机床不同结构状态的数控机床频响函数获取方法 - Google Patents
基于机床不同结构状态的数控机床频响函数获取方法 Download PDFInfo
- Publication number
- CN103970066B CN103970066B CN201410181574.3A CN201410181574A CN103970066B CN 103970066 B CN103970066 B CN 103970066B CN 201410181574 A CN201410181574 A CN 201410181574A CN 103970066 B CN103970066 B CN 103970066B
- Authority
- CN
- China
- Prior art keywords
- lathe
- movable member
- machine tool
- state
- mode
- 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.)
- Active
Links
Landscapes
- Automatic Control Of Machine Tools (AREA)
Abstract
一种基于机床不同结构状态的数控机床频响函数获取方法,该方法包括,步骤(1)将整个机床划分为两部件,包括固定部件和可动部件;步骤(2)获得可动部件在第一位置时机床结构状态的模态参数;步骤(3)获得可动部件在第二位置时机床结构状态的模态参数;步骤(4)利用不同结构状态的模态参数先计算得出模态标定因子,后结合频响函数的模态表达式合成机床在该结构状态的频响函数,本发明方法不需要测量输入力即可获得频响函数,应用本发明方法可实时监控机床的动态特性。
Description
技术领域
本发明属于数控装备结构状态动力学分析领域,尤其涉及一种基于机床不同结构状态的数控机床频响函数获取方法。
背景技术
频响函数即频率响应函数,是互功率谱函数除以自功率谱函数得到的商。频响函数是复函数,它是被测***的动力学特征在频域范围的描述,也就是被测***本身对输入信号在频域中传递特性的描述。在数控机床中,用频响函数来表征机床结构状态的动力学特性,而机床的动力学特性对于改善其加工性能和结构状态的设计至关重要。
目前获得机床结构状态的动力学特性方法主要是基于强迫振动的实验模态分析法,该方法的主要特点是对机床结构状态施加特定的激励力,并同时测量该激励力和机床结构状态的振动响应以获得频响函数。该方法只能获得机床结构状态处于静止时的频响函数,然而,在实际应用当中,人们更关注的是机床在工作状态下的动力学特性,大量研究证明,机床结构状态在工作状态下的动力学特性会发生变化,与静止状态下机床结构状态的动力学特性完全不同。
近些年发展起来的工作模态分析法通过测量机床不同结构状态在工作时的振动响应来辨识机床结构状态的模态参数,其辨识结果较强迫振动的实验模态分析法更加精确可靠,更加适合机床运行状态下动力学特性的研究。然而,在工作模态分析中,获得固有频率、阻尼比以及未标定的模态振型后,还必须测定输入力,才能获得机床结构状态的频响函数,而测量输入力需要占用机床大量的工作时间,费时费力。
发明内容
针对以上现有技术的缺陷,本发明目的是提供一种基于机床不同结构状态的频响函数获取方法,通过改变机床可动部件的位置使机床结构状态发生变化,利用不同机床结构状态的工作模态参数计算各阶模态振型的标定因子,然后合成机床在该结构状态下的频响函数,其特征在于,该方法具体包括:
(1)将整个机床划分为两大部件,包括固定部件和可动部件,其中可动部件是机床中可通过程序控制位置的可移动部件,固定部件是机床中不可移动的部件;
(2)获得机床可动部件在第一位置时结构状态的工作模态参数;
(3)获得机床可动部件在第二位置时结构状态的工作模态参数;
(4)利用步骤(2)和步骤(3)获得的工作模态参数先计算得出模态标定因子,后结合频响函数的模态表达式分别合成可动部件在第一位置和第二位置时机床结构状态的频响函数。
进一步的,所述步骤(2)包括以下子步骤:
(2-1)控制机床可动部件移动到第一位置,然后对机床实施激励,通过加速度传感器测量机床的振动响应;
(2-2)通过机床的振动响应信号辨识出可动部件在第一位置时机床的工作模态参数,包括固有频率ωr、阻尼比ζr和未标定的模态振型向量ηr。
进一步的,所述步骤(3)包括以下子步骤:
(3-1)控制机床可动部件移动到第二位置,然后对机床实施激励,通过加速度传感器测量机床的振动响应;
(3-2)通过机床的振动响应信号辨识出可动部件在第二位置时机床的工作模态参数,包括固有频率ωr′、阻尼比ζr’和未标定的模态振型向量ηr′。
进一步的,所述的激励为空运行自激励和切削自激励,优选为空运行自激励。
进一步的,所述步骤(4)中,所述可动部件在第一位置时机床结构状态的频响函数[H1(ω)]采用以下公式计算:
式中,[H1(ω)]表示所述可动部件在第一位置时机床结构状态的频响函数,ar为所述机床结构状态的模态标定因子,r为第r阶模态,n为共有n阶模态,ωr表示可动部件在第一位置时机床固有频率,ζr表示可动部件在第一位置时机床阻尼比,ηr表示可动部件在第一位置时未标定的模态振型向量,T表示对向量进行的转置,i为复数。
进一步的,所述步骤(4)中,所述可动部件在第二位置时机床结构状态的频响函数[H2(ω)]采用以下公式计算:
式中,[H2(ω)]表示所述可动部件在第二位置时机床结构状态的频响函数,ar为所述机床结构状态的模态标定因子,r为第r阶模态,n为共有n阶模态,ωr′表示可动部件在第二位置时机床固有频率,ζr’表示可动部件在第二位置时机床阻尼比,ηr′表示可动部件在第二位置时未标定的模态振型向量,T表示对向量进行的转置,i为复数。
进一步的,所述步骤(4)中,所述的模态标定因子ar采用以下公式计算:
式中,模态标定因子ar为可动部件在第一位置和第二位置时机床的模态标定因子,m是可动部件质量组成的质量矩阵,ωr表示可动部件在第一位置时机床固有频率,ηr表示可动部件在第一位置时未标定的机床模态振型向量,ωr′表示可动部件在第二位置时机床固有频率,ηr′表示可动部件在第二位置时未标定的模态振型向量,T表示对向量进行的转置。
本发明中可动部件的选择是根据频响函数需要而决定的,根据不同可动部件各自位置的改变分别获得机床在这两个位置对应的频响函数,本发明具有以下的优点和技术效果:
(1)本发明对频响函数的获取不需要测量输入力,不占用机床的大量的工作时间,节约了成本和时间。
(2)利用本发明方法,可以快速获取所有可能的结构状态下机床的频响函数,即整个机床的完整频响函数。
(3)本发明方法可以实时测量机床的频响函数,对于监测机床的动态特性有重要的意义。
附图说明
图1为本发明实施例中机床结构状态划分和可动部件位置示意图;
图2为本发明实施例中机床滑块随机加减速运动示意图;
图3为本发明中合成的机床结构状态频响函数曲线与实测的频响函数曲
线对比图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明,应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
本实施例中以XHK5140型立式加工中心为例对本发明的方法进行说明,基于机床不同结构状态的数控机床频响函数获取方法包括以下步骤:
(1)图1为XHK5140型立式加工中心的固定部件和可动部件位置示意图。将整个机床划分为固定部件和可动部件,包括由床身与立柱组成的固定部件和由机头与工作台组成的可动部件。
(2)如图1(a)所示,移动可动部件到位置a,获得可动部件在位置a时机床结构状态的工作模态参数,具体包括以下来两个子步骤:
(2-1)可动部件在位置a时,通过数控机床程序控制滚珠丝杠上滑块按照设计在一个小的区域内进行随机加减速运动,即使机床滑块空运行,滑块随机运动产生的惯性力通过滚珠丝杠驱动***对机床实施激励,即对可动部件在位置a时的机床结构状态进行有效激振,通过加速度传感器测量机床的振动响应;
(2-2)通过机床的振动响应信号辨识出机床可动部件在位置a时结构状态的工作模态参数,包括固有频率ωr、阻尼比ζr和未标定的模态振型向量ηr。
(3)如图1(b)所示,移动可动部件到位置b,获得可动部件在位置b时机床结构状态的工作模态参数,具体包括以下来两个子步骤:
(3-1)可动部件在位置b时,通过数控机床程序控制滚珠丝杠上滑块按照设计在一个小的区域内进行随机加减速运动,即使机床滑块空运行,滑块随机运动产生的惯性力通过滚珠丝杠驱动***对机床实施激励,即对可动部件在位置b时的机床结构状态进行有效激振,通过加速度传感器测量机床的振动响应;
(3-2)通过机床的振动响应信号辨识出机床可动部件在位置b时结构状态的工作模态参数,包括固有频率ωr′,阻尼比ζr’和未标定的模态振型向量ηr′。
(4)利用不同结构状态的模态参数计算得出模态标定因子,并结合频响函数的模态表达式合成机床结构状态的频响函数。
其中,计算模态标定因子的计算过程如下,
机床可动部件在位置a时的结构状态对应的运动方程为:
式中,MA、CA和KA分别表示可动部件在位置a时机床结构状态对应的质量矩阵、阻尼矩阵和刚度矩阵,X是固定部件的位移向量,x是可动部件的位移向量。
在不考虑阻尼的情况下,式(1)的特征值方程为:
式中,表示机床固定部件已标定的模态振型向量,φr表示可动部件已标定的模态振型向量,ωr为可动部件在位置a时机床固有频率。
同理,机床可动部件在位置b时结构状态对应的特征值方程为:
式中,MB和KB分别表示可动部件在位置b时机床结构状态对应的质量矩阵和刚度矩阵,表示固定部件已标定的模态振型向量,φr′表示可动部件已标定的模态振型向量,ωr′为可动部件在位置b时机床固有频率。
式(2)减去式(3)得到如下式(4):
通过将式(4)化简,并两边同时乘以可得式(5):
进一步计算得式(6):
其中M表示固定部件的质量矩阵,m表示可动部件的质量矩阵。因为选择的数控机床的固定部件和可动部件是确定的,所以可动部件无论在位置a或位置b,其固定部件的质量矩阵M和可动部件的质量矩阵m是一定的。
已标定的模态振型向量和未标定的模态振型向量可以表示为:
其中αr是模态标定因子,表示已标定的模态振型向量和未标定的模态振型向量之间的线性关系。
将式(7)代入式(6)中,计算出可动部件在位置a和位置b时机床结构状态的模态标定因子αr如下式:
结合频响函数的模态表达式合成机床结构状态的频响函数过程如下,
频响函数的模态表达式为:
式中,表示已标定的模态振型向量,ωr表示固有频率,ζr表示阻尼比。
将模态标定因子αr代入频响函数的模态表达式中,合成可动部件在位置a时机床结构状态对应的频响函数,如下式:
式中,[H1(ω)]表示所述可动部件在第一位置时机床结构状态的频响函数,ar为所述机床结构状态的模态标定因子,r为第r阶模态,n为共有n阶模态,ωr表示可动部件在第一位置时机床固有频率,ζr表示可动部件在第一位置时机床阻尼比,ηr表示可动部件在第一位置时未标定的模态振型向量,T表示对向量进行的转置,i为复数。
合成可动部件在位置b时机床结构状态的频响函数,如下式:
式中,[H2(ω)]表示所述可动部件在第二位置时机床结构状态的频响函数,ar为所述机床结构状态的模态标定因子,r为第r阶模态,n为共有n阶模态,ωr′表示可动部件在第二位置时机床固有频率,ζr’表示可动部件在第二位置时机床阻尼比,ηr′表示可动部件在第二位置时未标定的模态振型向量,T表示对向量进行的转置,i为复数。
如图3为合成的可动部件在位置a时机床的频响函数曲线与实测的频响函数曲线对比图,两者具有很好的吻合度,说明本发明中频响函数的获取方法准确合理。
本发明中,对机床的激励不限于空运行激励,还可采用切削自激励等其它激励方式,采用切削自激励获得频响函数的方法与采用空运行激励相同。
本发明中,当需要获得可动部件在任何其他位置时机床的频响函数时,移动可动部件到指定位置,按照本发明提供的方法,计算得到模态标定因子后,合成相应的频响函数即可,因此采用本发明方法可以快速获取所有可能的结构状态下机床频响函数,即整个机床的完整频响函数。
以上所述仅是本发明的较佳实施方式,故凡依本发明申请范围所述的特征及原理所做的等效变化或修饰,均包括于本发明保护范围内。
Claims (6)
1.一种基于机床不同结构状态的数控机床频响函数获取方法,通过改变机床可移动部件的位置使机床结构状态发生变化,利用机床不同结构状态的工作模态参数计算各阶模态振型的标定因子,然后合成机床在不同结构状态下的频响函数,其特征在于,该方法具体包括以下步骤:
(1)将整个机床划分为两大部件,包括固定部件和可动部件,其中可动部件是机床中可通过程序控制位置的可移动部件,固定部件是机床中不可移动的部件;
(2)获得机床可动部件在第一位置时结构状态的工作模态参数;
(3)获得机床可动部件在第二位置时结构状态的工作模态参数;
(4)利用步骤(2)和步骤(3)获得的工作模态参数先计算得出模态标定因子,后结合频响函数的模态表达式分别合成可动部件在第一位置和第二位置时机床结构状态的频响函数,
所述的模态标定因子αr采用以下公式计算:
式中,模态标定因子αr为可动部件在第一位置和第二位置时机床的模态标定因子,m是可动部件质量组成的质量矩阵,ωr表示可动部件在第一位置时机床固有频率,ηr表示可动部件在第一位置时未标定的机床模态振型向量,ωr′表示可动部件在第二位置时机床固有频率,ηr′表示可动部件在第二位置时未标定的模态振型向量,T表示对向量进行的转置。
2.根据权利要求1所述的方法,其特征在于,所述步骤(2)包括以下子步骤:
(2-1)控制机床可动部件移动到第一位置,然后对机床实施激励,通过 加速度传感器测量机床的振动响应;
(2-2)通过机床的振动响应信号辨识出可动部件在第一位置时机床的工作模态参数,包括固有频率ωr、阻尼比ζr和未标定的模态振型向量ηr。
3.根据权利要求1所述的方法,其特征在于,所述步骤(3)包括以下子步骤:
(3-1)控制机床可动部件移动到第二位置,然后对机床实施激励,通过加速度传感器测量机床的振动响应;
(3-2)通过机床的振动响应信号辨识出可动部件在第二位置时机床的工作模态数,包括固有频率ωr′、阻尼比ζr′和未标定的模态振型向量ηr′。
4.根据权利要求2或3所述基于机床不同结构状态的数控机床频响函数获取方法,所述的激励为空运行自激励和切削自激励。
5.根据权利要求4所述的方法,所述步骤(4)中,所述可动部件在第一位置时机床结构状态的频响函数[H1(ω)]采用以下公式计算:
式中,[H1(ω)]表示所述可动部件在第一位置时机床结构状态的频响函数,αr为所述机床结构状态的模态标定因子,r为第r阶模态,n为共有n阶模态,ωr表示可动部件在第一位置时机床固有频率,ζr表示可动部件在第一位置时机床阻尼比,ηr表示可动部件在第一位置时未标定的模态振型向量,T表示对向量进行的转置,i为复数。
6.根据权利要求5所述的方法,所述步骤(4)中,所述可动部件在第二位置时机床结构状态的频响函数[H2(ω)]采用以下公式计算:
式中,[H2(ω)]表示所述可动部件在第二位置时机床结构状态的频响函 数,αr为所述机床结构状态的模态标定因子,r为第r阶模态,n为共有n阶模态,ωr′表示可动部件在第二位置时机床固有频率,ζr′表示可动部件在第二位置时机床阻尼比,ηr′表示可动部件在第二位置时未标定的模态振型向量,T表示对向量进行的转置,i为复数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410181574.3A CN103970066B (zh) | 2014-04-30 | 2014-04-30 | 基于机床不同结构状态的数控机床频响函数获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410181574.3A CN103970066B (zh) | 2014-04-30 | 2014-04-30 | 基于机床不同结构状态的数控机床频响函数获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103970066A CN103970066A (zh) | 2014-08-06 |
CN103970066B true CN103970066B (zh) | 2017-01-11 |
Family
ID=51239700
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410181574.3A Active CN103970066B (zh) | 2014-04-30 | 2014-04-30 | 基于机床不同结构状态的数控机床频响函数获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103970066B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109884985A (zh) * | 2019-03-11 | 2019-06-14 | 上海理工大学 | 数控机床整机加工状态动态特性的测量方法 |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106052743B (zh) * | 2016-05-17 | 2018-05-11 | 湖北工业大学 | 一种评估传感器质量对频响函数影响大小的方法 |
CN106502199B (zh) * | 2016-10-20 | 2019-01-04 | 华中科技大学 | 一种机械装备结构运行状态下的频响函数识别方法 |
CN108469784A (zh) * | 2018-03-07 | 2018-08-31 | 上海理工大学 | 适用于数控机床加工状态的模态参数的测量装置及方法 |
CN111680368B (zh) * | 2019-02-25 | 2023-07-25 | 中国石油天然气集团有限公司 | 底部框架支撑式塔器结构的获取方法及装置 |
CN112880813B (zh) * | 2021-01-20 | 2023-06-02 | 长沙理工大学 | 一种混凝土泵车臂架固有频率计算方法及装置 |
CN114789472B (zh) * | 2022-05-23 | 2024-02-02 | 华中科技大学 | 一种机器人频响函数的获取方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101718613A (zh) * | 2009-11-12 | 2010-06-02 | 东莞华中科技大学制造工程研究院 | 一种数控装备实验模态分析方法 |
-
2014
- 2014-04-30 CN CN201410181574.3A patent/CN103970066B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101718613A (zh) * | 2009-11-12 | 2010-06-02 | 东莞华中科技大学制造工程研究院 | 一种数控装备实验模态分析方法 |
Non-Patent Citations (3)
Title |
---|
On exact and approximated formulations for scaling-mode shapes in operational modal analysis by mass and stiffness change;M.Lo´pez-Aenlle;《Journal of Soundand Vibration》;20120130;第331卷(第3期);全文 * |
基于响应信号的数控机床模态参数提取方法研究;李长玉;《中国优秀博硕士学位论文全文数据库(硕士)工程科技I辑》;20090615(第6期);正文第6页第1.4节第1段、第3段,第11页第2.1.1节第3、5段,第13页第2.2.2节第1段,第17页第2.3节第3-6行,第26页第第3.33节(2) * |
基于无线传感器考虑施工阶段的大跨度刚构桥的模态分析;陈国瑞;《中国优秀博硕士学位论文全文数据库(硕士)工程科技II辑》;20130115(第1期);正文第14-15页 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109884985A (zh) * | 2019-03-11 | 2019-06-14 | 上海理工大学 | 数控机床整机加工状态动态特性的测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103970066A (zh) | 2014-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103970066B (zh) | 基于机床不同结构状态的数控机床频响函数获取方法 | |
CN104786101A (zh) | 立铣切削颤振的监测方法 | |
CN106960068A (zh) | 一种基于脉冲激励响应频谱的模态阻尼比快速计算方法 | |
CN107870074B (zh) | 一种纤维复合板非线性内共振表征测试方法及测试*** | |
CN106502199B (zh) | 一种机械装备结构运行状态下的频响函数识别方法 | |
CN105108583B (zh) | 基于隔振式电涡流传感器测量动态铣削力装置及测量方法 | |
CN104915498A (zh) | 基于模型识别与等效简化的高速平台运动参数自整定方法 | |
CN103175602B (zh) | 基于单点激光连续平面扫描测振的模态测试***及方法 | |
CN103543012B (zh) | 一种滚动轴承刚度测试实验装置及其实验方法 | |
CN104777054A (zh) | 一种基于软测量技术的谐振式疲劳裂纹扩展试验振动***的参数识别方法 | |
CN108520125A (zh) | 一种预测刀具磨损状态的方法及*** | |
CN104897268B (zh) | 基于激光扫描的高档数控机床模态振型测试装置及方法 | |
CN106052743B (zh) | 一种评估传感器质量对频响函数影响大小的方法 | |
CN106556501B (zh) | 基于疲劳损伤等效分析的多轴振动试验条件裁剪方法 | |
CN106289947A (zh) | 一种轻质高强梁结构损伤识别方法 | |
CN101975633A (zh) | 应用连续细化分析傅立叶变换法的发动机激励力测量方法 | |
Ghafoori et al. | Dynamic responses of a rectangular plate under motion of an oscillator using a semi-analytical method | |
Bouguerriou et al. | Novel cyclostationarity-based blind source separation algorithm using second order statistical properties: Theory and application to the bearing defect diagnosis | |
CN102744649A (zh) | 一种数控车床主传动***非切削能耗获取方法 | |
CN109884985A (zh) | 数控机床整机加工状态动态特性的测量方法 | |
He et al. | A new approach based on modal mass distribution matrix to identify weak components of machine tool structure | |
CN109190328A (zh) | 一种混合有限-复变差分的多参数结构动响应灵敏度分析方法 | |
CN106096033A (zh) | 足底压力信息分析方法、装置及*** | |
Wang et al. | Extraction of real modes and physical matrices from modal testing | |
Wang et al. | A method to identify the main mode of machine tool under operating conditions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |