CN107102620B - 一种数控机床双驱进给系统模态建模方法 - Google Patents
一种数控机床双驱进给系统模态建模方法 Download PDFInfo
- Publication number
- CN107102620B CN107102620B CN201710283520.1A CN201710283520A CN107102620B CN 107102620 B CN107102620 B CN 107102620B CN 201710283520 A CN201710283520 A CN 201710283520A CN 107102620 B CN107102620 B CN 107102620B
- Authority
- CN
- China
- Prior art keywords
- lead screw
- equation
- energy
- nut
- feed system
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
- G05B19/408—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by data handling or data format, e.g. reading, buffering or conversion of data
- G05B19/4086—Coordinate conversions; Other special calculations
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/35—Nc in input of data, input till input file format
- G05B2219/35356—Data handling
Abstract
本发明提出了一种数控机床双驱进给系统模态建模方法:步骤1、确定系统的广义坐标;步骤2、建立机床进给系统的动能表达式;步骤3、建立进给系统的势能表达式;步骤4、建立进给系统耗能表达式;步骤5、建立滑块坐标与工作台坐标的的位置转化关系式;步骤6、将能量表达式带入拉格朗日方程,得到进给系统的运动微分方程;步骤7、求解模态参数:求解模态固有频率;本发明的建模方法更加直观,变量可以控制,根据变量的改变可以观察固有频率的改变,变量与固有频率之间的关系更加直观。
Description
技术领域
本发明属于机械振动分析技术领域,特别是一种数控机床双驱进给系统模态建模方法。
背景技术
作为机床重要部件之一的进给系统,其动态性能是决定机床性能的重要因素,在设计阶段就需要考虑机床的动态性能,机床模态要远离共振区域,减小共振对加工的影响。目前主要有建模和有限元软件模拟两种方法。
许向荣(许向荣.滚珠丝杠副直线导轨进给单元动态性能研究[D].济南:山东大学,2011)考虑滚动结合面的影响,对滚珠丝杠副直线导轨单元进行动力学建模,根据拉格朗日方程推导了进给单元的微分振动方程,获得了进给单元的质量矩阵、刚度矩阵及阻尼矩阵。张会端(张会端.机床进给系统的动力学分析[D].吉林:吉林大学,2009)根据机床进给系统的特点,建立其动力学模型,列出系统的动能、势能、耗能函数和外载荷所作的虚功,利用拉格朗日方程推导系统的运动学方程。上述研究由于均把系统当做连续系统研究,导致无法求解出具体的振型函数,所以无法求解出模态参数。王永强(王永强.滚珠丝杠进给系统自适应建模理论与方法研究[D].济南:山东大学,2013)运用拉个朗日方程建立了滚珠丝杠进给系统的二维动力学模型,求解出系统的一阶固有频率,由于建立的模型是在平面上的单自由度模型,所以不能求解二阶及以上的固有频率。
在模态分析中,利用有限元软件进行分析也是常见的手段之一,在边界条件设置准确的情况下,有限元软件分析可以做到比较精准的分析。但是也存在明显的局限性:软件的通用性不强,每次分析都需要准确设置边界情况,操作比较繁琐;软件后处理只输出结果,对于每个不同的模型会有不同的输出结果,但是不能显示出对模态影响大的因素。
发明内容
本发明所解决的技术问题在于提供一种数控机床双驱进给系统模态建模方法,以解决现有技术中连续系统振型方程无法求解和有限元软件模拟操作繁琐等问题;本发明提供一种能够具有可操作性的且求解效率更高的模态分析方法。
实现本发明目的的技术解决方案为:
一种数控机床双驱进给系统模态建模方法,包括以下步骤:
步骤1、确定系统的广义坐标:判断系统的自由度数,并选取广义坐标qj来描述系统的运动状态,其坐标数目和自由度数相同;
2、步骤2、建立机床进给系统的动能表达式:即分别建立丝杠的动能和工作台的动能表达式,进一步建立机床进给系统总的动能表达式;
步骤3、建立进给系统的势能表达式:即分别建立导轨和滑块接触面的势能、轴承内外圈接触面的势能、丝杠和螺母接触面的势能、丝杠的轴向势能和扭转势能;进一步建立进给系统总势能表达式;
步骤4、建立进给系统耗能表达式:即分别建立导轨和滑块接触面的耗能、轴承内外圈接触面的耗能、丝杠和螺母接触面的耗能、丝杠的耗能;进一步建立进给系统总耗能表达式;
步骤5、建立滑块坐标与工作台坐标的的位置转化关系式:通过坐标转换得到滑块坐标与工作台坐标的数学关系式;
步骤6、将能量表达式带入拉格朗日方程,得到进给系统的运动微分方程;将运动微分方程写成矩阵形式:
式中M为质量矩阵,C为阻尼矩阵,K为刚度矩阵,F为激振力矩阵;
步骤7、求解模态固有频率;简化方程(7)为判断矩阵行列式的值是否恒为零,分别采用矩阵迭代法和根据来求解模态固有频率。
本发明与现有技术相比,其显著优点:
(1)本发明的建模方法更加直观:变量可以控制,根据变量的改变可以观察模态的改变,变量与模态之间的关系更加直观。
(2)本发明的建模方法改进更具有针对性:从质量矩阵M和刚度矩阵K的数值可以分析推导出对模态影响较大的参数,在改进的时候可以通过分析薄弱环节,有针对性的改进。
(3)本发明的建模方法具有更好的可操作性:由能量法求解进给系统固有频率的方法,具有良好的可操作性,求解的效率比较高,借助于正交试验法可以实现机床进给系统的固有频率反向设计,对比有限元仿真具有一定的优势。
下面结合附图对本发明作进一步详细描述。
附图说明
图1为本发明方法流程示意图。
图2为进给系统位置关系和系统坐标图。
图3(a-c)分别为滑块及工作台坐标的空间示意图、主视图、俯视图。
具体实施方式
为了说明本发明的技术方案及技术目的,下面结合附图及具体实施例对本发明做进一步的介绍。
结合图1的流程示意图,本发明的一种数控机床双驱进给系统模态建模方法,包括以下步骤:
步骤1、确定系统的广义坐标:判断系统的自由度数,并选取广义坐标来描述系统的运动状态,其坐标数目和自由度数相同。
结合图2,以处于静止状态时的工作台重心位置为坐标原点,建立O-XYZ主坐标系,X轴的正方向垂直于导轨,Y轴的正方向垂直于工作台上表面向上,Z轴的正方向为工作台远离电机的方向,Z轴方向为工作台进给方向。
建立广义坐标:
qj={xw,yw,zw,α,β,γ,z1,z2,z3,θ1,θ2,θ3,x,y,…}j=(1,2,3,....n);
其中j为广义坐标系qj的坐标编号,如q1={xw};“…”为省略的其他自由度,n为系统的自由度数目。
其中,xw、yw、zw分别为工作台在X、Y、Z方向上的横向振动位移;α、β、γ分别为工作台绕X、Y、Z方向的扭转角度;z1为丝杠左端(远离电机一端)的轴向位移,θ1为丝杠绕Z方向的转动角度;z2为螺母中心的轴向位移,θ2为螺母绕Z方向的转动角度;z3为丝杠右端(靠近电机一端)的轴向位移,θ3为丝杠绕Z方向的转动角度;x、y分别为丝杠的径向平行与X轴和Y轴的位移;
下面均以自由度数n=14为例,即
qj={xw,yw,zw,α,β,γ,z1,z2,z3,θ1,θ2,θ3,x,y}j=(1,2,3,....n)。
步骤2、建立机床进给系统的动能模型:即分别建立丝杠的动能和工作台的动能表达式,进一步建立机床进给系统总的动能表达式;
步骤2.1、建立工作台的动能方程:
工作台的动能包括在X、Y、Z方向上的横向振动动能和绕X、Y、Z方向的扭转动能,工作台的动能Tw表达式为:
式中mw为工作台的质量,Jx为工作台绕X轴的转动惯量,Jy工作台绕Y轴的转动惯量,Jz工作台绕Z轴的转动惯量。
步骤2.2建立丝杠的动能方程:
丝杠的动能包括丝杠的扭转动能、轴向动能以及横向动能,
(1)建立丝杠轴向动能:
对于螺母左端丝杠,设距离丝杠左端为z的位置有长度为dz的极小一段,则该位置在Z向的位移方程为:
其中b为螺母左端丝杠长度;
对上述公式进行积分得到螺母左端丝杠的轴向动能Tl表达式:
同理,对螺母右端丝杠进行积分运算,右端丝杠动能Tr为
mbs为丝杠质量,a为丝杆的有效长度(两个轴承之间丝杠的长度);
(2)建立丝杠横向动能表达式为:
(3)建立丝杠扭转动能表达式为:
Jbsl为螺母左端丝杠的转动惯量,Jbsr为螺母右端丝杠的转动惯量;
(4)建立丝杠总动能Tbs表达式为:
2.3、机床进给系统总动能为:
步骤3、建立进给系统的势能表达式:即分别建立导轨和滑块接触面的势能、轴承内外圈接触面的势能、丝杠和螺母接触面的势能、丝杠的轴向势能和扭转势能;进一步建立进给系统总势能表达式;
步骤3.1、建立导轨和滑块接触面的势能方程:
导轨和滑块接触面的势能Us表达式为:
其中,Ksx为滑块与导轨侧向接触刚度,Ksy为滑块与导轨法向接触刚度。
步骤3.2、建立轴承内外圈结合面势能方程:
轴承内外圈结合面势能Ub表达式为:
其中,Kbx、Kby为轴承径向刚度,Kbz为轴承轴向刚度。
步骤3.3、建立丝杠螺母结合面势能势能方程:
杠螺母结合面势能势能Un表达式为:
其中,Knx、Kny为丝杠螺母副径向接触刚度,Knz为丝杠螺母副轴向接触刚度。
步骤3.4、建立丝杠的势能方程:
丝杠的势能包括丝杠的轴向势能和扭转势能:Ubs=Uz+Uθ
(1)丝杠轴向势能Uz表达式:Uz=Kbszl(z2-z1)2+Kbszr(z3-z2)2
其中,Kbszl为螺母左侧丝杠轴向刚度,Kbszr为螺母右侧丝杠轴向刚度。
(2)丝杠扭转势能Uθ表达式:Uθ=Kθl(θ2-θ1)2+Kθr(θ3-θ2)2
其中,Kθl为螺母左侧丝杠扭转刚度,Kθr为螺母右侧丝杠扭转刚度。
步骤3.5、建立进给系统总势能的表达式:
步骤4、建立进给系统耗能表达式:即分别建立导轨和滑块接触面的耗能、轴承内外圈接触面的耗能、丝杠和螺母接触面的耗能、丝杠的耗能;进一步建立进给系统总耗能表达式;
步骤4.1、建立导轨和滑块接触面耗能方程:
轨和滑块接触面耗能Ds表达式:其中,Csx为滑块与导轨侧向接触阻尼,Csy为滑块与导轨法向接触阻尼;xsi为第i个滑块在X方向的位移,ysi为第i个滑块在Y方向的位移;
步骤4.2、建立轴承内外圈结合面耗能方程:
轴承内外圈结合面耗能Db表达式:
其中,Cbx、Cby为轴承径向阻尼,Cbz为轴承轴向阻尼。
步骤4.3、建立丝杠螺母结合面耗能方程:
丝杠螺母结合面耗能Dn表达式:
其中,Cnx、Cny为丝杠螺母副径向接触阻尼,Cnz为丝杠螺母副轴向接触阻尼。
步骤4.4、建立丝杠的耗能方程:
丝杠的耗能Dbs表达式:
其中,Cbsl为螺母左侧丝杠轴向阻尼,Cbsr为螺母右侧丝杠轴向阻尼。
步骤4.5、建立进给系统总耗能的表达式:
步骤5、建立滑块坐标与工作台坐标的的位置转化关系式:建立动态坐标系,通过坐标转换得到滑块坐标与工作台坐标在动态坐标系的数学关系式;
结合图3滑块与工作台坐标示意图,以工作台任意时刻的重心为坐标原点,建立O1-X1Y1Z1动态坐标系,X1、Y1、Z1的方向分别与X、Y、Z方向一致。四个滑块中心点分别为s1,s2,s3,s4,第i个滑块中心与原点O在X方向上的距离为dsi,在Y方向上的距离为hsi,在Z方向上的距离为lsi。工作台在X、Y、Z三个方向上的横向振动位移分别为xw、yw、zw,绕X、Y、Z的扭转位移分别为α、β、γ。
步骤5.1通过坐标转换得到滑块坐标与工作台坐标在动态坐标系下的关系为:
5.2由于α、β、γ的值很小,所以tanα=α,tanβ=β,tanγ=γ,所以上式可以化简为:
步骤6:将能量表达式带入拉格朗日方程,
得到进给系统的运动微分方程;将运动微分方程写成矩阵形式。
其中qj为系统的坐标,为系统的速度,T为系统的动能,U为系统的势能,t为时间,Q'j为对应的力,n为系统的自由度数。
下面仍以自由度n=14为例,分别计算广义坐标
qj={xw,yw,zw,α,β,γ,z1,z2,z3,θ1,θ2,θ3,x,y}中每个坐标系下拉格朗日方程中各项导数;
(1)以q1={xw}为例,联立方程式(1)-(6)对xw计算拉格朗日方程中各项导数:
将上述各式整理得运动微分方程为:
(2)对广义坐标qj={xw,yw,zw,α,β,γ,z1,z2,z3,θ1,θ2,θ3,x,y}每个坐标分别求导,微分方程为:
式中F1-F14分别为每组方程对应的激振力,由于求解过程中令方程右端为零,所以F1-F14值未知对求解过程没有影响。
(3)将运动微分方程写成矩阵形式
式中M为质量矩阵,C为阻尼矩阵,K为刚度矩阵,F为激振力矩阵。其中:
步骤7、求解模态固有频率:
步骤7.1、矩阵行列式的值不恒为零时解法:
由于有阻尼的固有频率ωd和无阻尼的固有频率ωn存在的关系,当进给系统阻尼比ξ=10%时,ωd=0.99499ωn,且一般情况下进给系统阻尼比ξ≤10%,所以求解进给系统固有频率时不考虑阻尼矩阵C,同时也不考虑作用力F。所以方程(7)化简为
式(8)有的解为
X为振幅,ωn为固有频率,为相位角。
将式(9)解代入式(8),得到如下振型方程
上式存在非零解的条件是:当且仅当系数行列式为零,即:
式中Δ(ω2)为特征值行列式,将其展开后可以得到ωn的n次代数方程式:
ωn 2n+a1ωn 2(n-1)+a2ωn 2(n-2)+…+an-1ωn 2+an=0
式中a1,a2…an为常系数。
这一方程有n个根ωn 2(r=1,2,…n),这些根称为特征值,特征值的平方根ωn称为进给系统的固有频率。
步骤7.2、矩阵行列式的值恒为零时的解法:
特别地:在机床进给系统中,当工作台位于丝杠中点时矩阵为半正定矩阵,所以矩阵行列式的值恒为零,此时用上述方法无法求解方程,需采用矩阵迭代法求解。同时由于刚度矩阵是半正定矩阵,所以用Matlab中pinv语言求解刚度矩阵的伪逆矩阵,通过迭代求取工作台位于中点时进给系统的固有频率。
本方法具有一定的通用性,适用于双导轨双丝杠四滑块进给系统,符合上述要求的进给系统都可以用此建模方法求解。
例1.以MCH63数控机床进给系统为例,整理得到矩阵中所需参数(具体参见表3),通过MATLAB编程求解得到系统固有频率参数,整理之后得到进给系统前六阶固有频率计算值,如表1所示。
表1.MCH63数控机床进给系统参数
K<sub>sx</sub>(N/m) | K<sub>sy</sub>(N/m) | K<sub>bx</sub>(N/m) | K<sub>by</sub>(N/m) | K<sub>bz</sub>(N/m) | K<sub>bszr</sub>(N/m) | K<sub>bs</sub>(N/m) |
7.7*10<sup>8</sup> | 3.45*10<sup>9</sup> | 2.4*10<sup>8</sup> | 2.4*10<sup>8</sup> | 3.2*10<sup>8</sup> | 3.2*10<sup>8</sup> | 1.6*10<sup>8</sup> |
K<sub>nz</sub>(N/m) | K<sub>bszl</sub>(N/m) | K<sub>θl</sub>(N/m) | K<sub>θr</sub>(N/m) | K<sub>θ</sub>(N/m) | K<sub>nx</sub>(N/m) | K<sub>ny</sub>(N/m) |
1*10<sup>9</sup> | 3.2*10<sup>8</sup> | 1.734*10<sup>4</sup> | 1.734*10<sup>4</sup> | 8.67*10<sup>3</sup> | 8*10<sup>8</sup> | 8*10<sup>8</sup> |
d<sub>s</sub>(m) | h<sub>s</sub>(m) | J<sub>x</sub>(kgm<sup>2</sup>) | J<sub>y</sub>(kgm<sup>2</sup>) | J<sub>z</sub>(kgm<sup>2</sup>) | a(m) | b(m) |
0.327m | 0.170 | 45.732 | 85.229 | 58.875 | 1.2 | 0~1.2 |
m<sub>w</sub>(kg) | m<sub>bs</sub>(kg) | J<sub>bs</sub>(kgm<sup>2</sup>) | J<sub>bsl</sub>(kgm<sup>2</sup>) | J<sub>bs</sub>(kgm<sup>2</sup>) | E | G |
727.64 | 13.267 | 2.540*10<sup>-3</sup> | J<sub>bs</sub> b/a | 2.540*10<sup>-3</sup> | 2.06*10<sup>11</sup> | 7.938*10<sup>10</sup> |
l<sub>s1</sub>=l<sub>s2</sub>(m) | l<sub>s3</sub>=l<sub>s4</sub>(m) | ρ<sub>bs</sub>(kg/m<sup>3</sup>) | ρ<sub>w</sub>(kg/m<sup>3</sup>) | J<sub>bsr</sub>(kgm<sup>2</sup>) | ||
0.311 | 0.271 | 7850 | 7300 | J<sub>bs</sub>(a-b)/a |
表中ds为滑块与原点O在X方向上的距离,hs为滑块与原点O在Y方向上的距离,E为丝杠材料弹性模量,G为丝杠材料剪切弹性模量(可由E、G计算出刚度值K,上表已算出),ls1为第一个滑块与原点在Z方向的距离,ls2为第二个滑块与原点在Z方向的距离,ls3为第三个滑块与原点在Z方向的距离,ls4为第四个滑块与原点在Z方向的距离,ρbs为丝杠材料的密度,ρw为工作台材料的密度。
表2固有频率计算值
阶数 | 建模理论计算频率/HZ | 仿真计算频率/HZ | 误差 |
1 | 45.2 | 48.6 | 7.5% |
2 | 91.66 | 94.1 | 2.6% |
3 | 105.5 | 106.6 | 1.0% |
4 | 227.4 | 227.7 | 0.1% |
5 | 269.5 | 306.9 | 12.1% |
6 | 615.1 | 615.22 | 0.01% |
7 | 615.1 | 616.04 | 0.15% |
表2中建模理论计算结果与ANSYSY软件仿真结果误差最大为第五阶频率:12%,其他阶数的误差值都比较小,说明了建模方法具有一定的准确性。通过矩阵可以解耦出结构和固有频率的关系,从矩阵中的参数可以比较出对固有频率影响大的因素,在进行优化时,可以首先优化这些参数,使优化的目标更加明确。
虽然本发明已以较佳实施例揭露如上,然其并非用以限定本发明。本发明所属技术领域中具有通常知识者,在不脱离本发明的精神和范围内,当可作各种的更动与润饰。因此,本发明的保护范围当视权利要求书所界定者为准。
Claims (8)
1.一种数控机床双驱进给系统模态建模方法,其特征在于,包括以下步骤:
步骤1、确定系统的广义坐标:判断系统的自由度数,并选取广义坐标qj来描述系统的运动状态,其坐标数目和自由度数相同;广义坐标为:
qj={xw,yw,zw,α,β,γ,z1,z2,z3,θ1,θ2,θ3,x,y,…}j=(1,2,3,....n);
其中j为广义坐标系qj的坐标编号,如q1={xw};n为系统的自由度数目;
xw、yw、zw分别为工作台在X、Y、Z方向上的横向振动位移;α、β、γ分别为工作台绕X、Y、Z方向的扭转位移;z1为丝杠左端的轴向位移,θ1为丝杠绕Z方向的转动角度;z2为螺母中心的轴向位移,θ2为螺母绕Z方向的转动角度;z3为丝杠右端的轴向位移,θ3为丝杠绕Z方向的转动角度;x、y分别为丝杠的径向平行与X轴和Y轴的位移;“…”为省略的其他自由度;
步骤2、建立机床进给系统的动能表达式:即分别建立丝杠的动能和工作台的动能表达式,进一步建立机床进给系统总的动能表达式;具体步骤为:
步骤2.1、建立工作台的动能方程:
式中mw为工作台的质量,Jx为工作台绕X轴的转动惯量,Jy工作台绕Y轴的转动惯量,Jz工作台绕Z轴的转动惯量;
步骤2.2、建立丝杠的动能方程:
其中,b为螺母左端丝杠长度;mbs为丝杠质量,a为丝杆的有效长度;Jbsl为螺母左端丝杠的转动惯量,Jbsr为螺母右端丝杠的转动惯量;z1为丝杠左端的轴向位移,θ1为丝杠绕Z方向的转动角度;z2为螺母中心的轴向位移,θ2为螺母绕Z方向的转动角度;z3为丝杠右端的轴向位移,θ3为丝杠绕Z方向的转动角度;
步骤2.3、机床进给系统总动能为:
步骤3、建立进给系统的势能表达式:即分别建立导轨和滑块接触面的势能、轴承内外圈接触面的势能、丝杠和螺母接触面的势能、丝杠的轴向势能和扭转势能;进一步建立进给系统总势能表达式;
步骤4、建立进给系统耗能表达式:即分别建立导轨和滑块接触面的耗能、轴承内外圈接触面的耗能、丝杠和螺母接触面的耗能、丝杠的耗能;进一步建立进给系统总耗能表达式;
步骤5、建立滑块坐标与工作台坐标的的位置转化关系式:通过坐标转换得到滑块坐标与工作台坐标的数学关系式;
步骤6、将能量表达式带入拉格朗日方程,得到进给系统的运动微分方程;将运动微分方程写成矩阵形式:
式中M为质量矩阵,C为阻尼矩阵,K为刚度矩阵,F为激振力矩阵;qj为系统的坐标,为系统的速度;
步骤7、求解模态固有频率;简化方程(7)为判断矩阵行列式的值是否恒为零,分别采用矩阵迭代法和根据来求解模态固有频率。
2.如权利要求1所述的一种数控机床双驱进给系统模态建模方法,其特征在于,建立广义坐标的的过程为:以处于静止状态时的工作台重心位置为坐标原点,建立O-XYZ主坐标系,X轴的正方向垂直于导轨,Y轴的正方向垂直于工作台上表面向上,Z轴的正方向为工作台远离电机的方向,Z轴方向为工作台进给方向。
3.如权利要求1所述的一种数控机床双驱进给系统模态建模方法,其特征在于,步骤2.2建立丝杠的动能方程,具体步骤为:
(1)建立丝杠轴向动能:对于螺母左端丝杠,设距离丝杠左端为z的位置有长度为dz的极小一段,则该位置在Z向的位移方程为:
对上述公式进行积分得到螺母左端丝杠的轴向动能Tl表达式:
对螺母右端丝杠进行积分运算,右端丝杠动能Tr为:
(2)建立丝杠横向动能表达式为:
(3)建立丝杠扭转动能表达式为:
(4)建立丝杠总动能Tbs表达式。
4.如权利要求1所述的一种数控机床双驱进给系统模态建模方法,其特征在于,步骤3建立进给系统的势能表达式,具体步骤为:
步骤3.1、建立导轨和滑块接触面的势能方程:
其中,Ksx为滑块与导轨侧向接触刚度,Ksy为滑块与导轨法向接触刚度;xsi为第i个滑块在X方向的位移,ysi为第i个滑块在Y方向的位移;
步骤3.2、建立轴承内外圈结合面势能方程:
其中,Kbx、Kby为轴承径向刚度,Kbz为轴承轴向刚度;
步骤3.3、建立丝杠螺母结合面势能势能方程:
其中,Knx、Kny为丝杠螺母副径向接触刚度,Knz为丝杠螺母副轴向接触刚度;
步骤3.4、建立丝杠的势能方程:
(1)丝杠轴向势能Uz表达式:Uz=Kbszl(z2-z1)2+Kbszr(z3-z2)2
其中,Kbszl为螺母左侧丝杠轴向刚度,Kbszr为螺母右侧丝杠轴向刚度;
(2)丝杠扭转势能Uθ表达式:Uθ=Kθl(θ2-θ1)2+Kθr(θ3-θ2)2
其中,Kθl为螺母左侧丝杠扭转刚度,Kθr为螺母右侧丝杠扭转刚度;
步骤3.5、建立进给系统总势能的表达式。
5.如权利要求1所述的一种数控机床双驱进给系统模态建模方法,其特征在于,步骤4建立进给系统耗能表达式,具体步骤为:
步骤4.1、建立导轨和滑块接触面耗能方程:其中,Csx为滑块与导轨侧向接触阻尼,Csy为滑块与导轨法向接触阻尼;xsi为第i个滑块在X方向的位移,ysi为第i个滑块在Y方向的位移;
步骤4.2、建立轴承内外圈结合面耗能方程:
其中,Cbx、Cby为轴承径向阻尼,Cbz为轴承轴向阻尼;
步骤4.3、建立丝杠螺母结合面耗能方程:
其中,Cnx、Cny为丝杠螺母副径向接触阻尼,Cnz为丝杠螺母副轴向接触阻尼;
步骤4.4、建立丝杠的耗能方程:
其中,Cbsl为螺母左侧丝杠轴向阻尼,Cbsr为螺母右侧丝杠轴向阻尼;
步骤4.5、建立进给系统总耗能的表达式。
6.如权利要求2所述的一种数控机床双驱进给系统模态建模方法,其特征在于,步骤5建立滑块坐标与工作台坐标的的位置转化关系式,具体步骤为:
以工作台任意时刻的重心为坐标原点,建立O1-X1Y1Z1动态坐标系,X1、Y1、Z1的方向分别与X、Y、Z方向一致;第i个滑块中心与原点O在X方向上的距离为dsi,在Y方向上的距离为hsi,在Z方向上的距离为lsi;
步骤5.1通过坐标转换得到滑块坐标与工作台坐标在动态坐标系下的关系为:
5.2将上式可以化简为:
其中,xw、yw、zw分别为工作台在X、Y、Z三个方向上的横向振动位移,α、β、γ分别为绕X、Y、Z的扭转角度。
7.如权利要求1所述的一种数控机床双驱进给系统模态建模方法,其特征在于,步骤6将能量表达式带入拉格朗日方程,得到进给系统的运动微分方程,具体过程为:将能量表达式和广义坐标中每个坐标带入拉格朗日方程分别求导,得到系统的微分方程;
将运动微分方程写成矩阵形式
其中为系统的速度,T为系统的动能,U为系统的势能,t为时间,Q'j为对应的力,n为系统的自由度数。
8.如权利要求1所述的一种数控机床双驱进给系统模态建模方法,其特征在于,步骤7求解模态固有频率,具体过程为:将方程(7)化简为判断矩阵行列式的值是否恒为零;
(1)若非恒为零,根据展开后可以得到ωn的n次代数方程式:ωn 2n+a1ωn 2(n-1)+a2ωn 2(n-2)+…+an-1ωn 2+an=0
式中Δ(ω2)为特征值行列式,a1,a2…an为常系数,ωn 2(r=1,2,…n)为方程的根,即特征值,特征值的平方根ωn即为进给系统的固有频率;
(2)若恒为零,采用矩阵迭代法求解工作台位于中点时进给系统的固有频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710283520.1A CN107102620B (zh) | 2017-04-26 | 2017-04-26 | 一种数控机床双驱进给系统模态建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710283520.1A CN107102620B (zh) | 2017-04-26 | 2017-04-26 | 一种数控机床双驱进给系统模态建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107102620A CN107102620A (zh) | 2017-08-29 |
CN107102620B true CN107102620B (zh) | 2019-06-21 |
Family
ID=59657777
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710283520.1A Expired - Fee Related CN107102620B (zh) | 2017-04-26 | 2017-04-26 | 一种数控机床双驱进给系统模态建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107102620B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109141872B (zh) * | 2018-05-14 | 2020-05-01 | 重庆大学 | 基于复模态分析法的双排行星齿轮传动均载系数试验方法 |
CN109800447B (zh) * | 2018-12-06 | 2023-04-07 | 天津大学 | 一种基于有限元思想的进给系统刚柔耦合动力学建模方法 |
CN112799304B (zh) * | 2021-01-06 | 2022-06-14 | 武汉理工大学 | 一种基于时变摩擦补偿的双电机同步控制方法及设备 |
CN113239486B (zh) * | 2021-04-30 | 2023-10-17 | 中国华能集团清洁能源技术研究院有限公司 | 一种基于动力学分析的双风轮风机传动系统共振预测方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3383282B2 (ja) * | 2000-11-10 | 2003-03-04 | 東洋ゴム工業株式会社 | 6自由度シミュレーション方法 |
CN101598632B (zh) * | 2009-07-07 | 2010-12-01 | 华中科技大学 | 丝杠驱动进给系统中支撑点轴向刚度和阻尼的测量方法 |
CN102156033B (zh) * | 2011-03-24 | 2013-06-12 | 武汉理工大学 | 一种数控机床扭转振动模态的测量装置及其测量方法 |
CN103323248B (zh) * | 2013-07-04 | 2015-07-08 | 南京理工大学 | 角接触球轴承动、静态特性参数测试装置 |
CN104992038B (zh) * | 2015-07-27 | 2018-10-16 | 天津大学 | 一种刚柔-机电耦合进给系统的动态性能优化设计方法 |
-
2017
- 2017-04-26 CN CN201710283520.1A patent/CN107102620B/zh not_active Expired - Fee Related
Non-Patent Citations (2)
Title |
---|
重心驱动进给系统机电耦合建模及性能分析;李学言;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑(月刊)》;20170315(第3期);第8-14页 |
铣车复合加工中心双驱进给系统动力学分析与仿真;唐余林;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑(月刊)》;20141015(第10期);第19-38页 |
Also Published As
Publication number | Publication date |
---|---|
CN107102620A (zh) | 2017-08-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107102620B (zh) | 一种数控机床双驱进给系统模态建模方法 | |
CN106202647B (zh) | 电主轴的多轴疲劳寿命预测方法及疲劳寿命可靠性评估方法 | |
Cao et al. | A new dynamic model of ball-bearing rotor systems based on rigid body element | |
CN105608297A (zh) | 一种角接触球轴承动刚度的数值模拟方法 | |
Yang et al. | Static and dynamic characteristics modeling for CK61125 CNC lathe bed basing on FEM | |
Shi et al. | Cutting force and chatter stability analysis for PKM-based helical milling operation | |
Guo et al. | A tolerance analysis method for rotating machinery | |
Zhao et al. | Forced vibration of two-degrees-of-freedom machine tool feed system with clearance and friction | |
Yu et al. | Oil film shape prediction of hydrostatic thrust bearing under the condition of high speed and heavy load | |
Fegade et al. | Unbalanced response of rotor using ansys parametric design for different bearings | |
Sawicki et al. | A Nonlinear Model for Prediction of Dynamic Coefficients in a Hydrodynamic Journal Bearing. | |
CN102393679B (zh) | 一种多轴加工系统的相对动刚度获取方法及其应用 | |
Sheng et al. | Dynamic model and vibration characteristics of planar 3-RRR parallel manipulator with flexible intermediate links considering exact boundary conditions | |
Zhou et al. | Torsion vibration analysis of lead-screw feed drives with changeable table position and work-piece mass | |
Chen et al. | Tooth fracture detection in spiral bevel gears system by harmonic response based on finite element method | |
Feng et al. | Analysis of static and dynamic characteristic of spindle system and its structure optimization in camshaft grinding machine | |
Vukov et al. | Mechanic-Mathematical Model for Investigations of the Natural Frequencies and Mode Shapes of the Free Spatial Vibrations of Wood Shaper and its Spindle | |
Zhu et al. | Research on 3D chatter stability of blade by high-speed turn-milling | |
Fuzhong | CNC milling machine spindle characteristics analysis by FEM | |
Wu et al. | Nonlinear analysis of axial vibration of five-axis machine tool worktable with double turntable | |
Vukov et al. | Numerical Investigations of the Natural Frequencies and Mode Shapes of the Free Spatial Vibrations of Wood Shaper and its Spindle | |
Wang et al. | Research on Optimization of Face Gear Grinding Machine Structure Design Based on Dynamic Performance of Machine Tool | |
Shen et al. | Experimental investigation and improvement of dynamic performance of high-speed grinding machine | |
Zhang et al. | Modal analysis for a new double-side blade grinding machine | |
Xu et al. | Modeling and measurement of precision reducers' efficiency curve |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190621 Termination date: 20210426 |
|
CF01 | Termination of patent right due to non-payment of annual fee |