发明名称 |
一种机床固定结合部动态特性的检测方法 |
摘要 |
本发明公开了一种机床固定结合部动态特性的检测方法。该方法将固定结合部看成一种等截面的虚拟材料,虚拟材料与固定结合部两侧的零件皆为刚性连接。通过检测构成结合部零件的弹性模量、泊松比、密度、屈服强度、硬度和几何尺寸参数,得到该虚拟材料的弹性模量、泊松比、密度、厚度,虚拟材料与固定结合部两侧的二个零件刚性连接,根据计算出的虚拟材料弹性模量、泊松比、密度总共3个参数,输入到有限元软件中,检测出含结合部的复杂部件的一些动态特性(如振型、固有频率、位移等)。理论与实验前6阶固有频率的相对误差在(-10~10)%之间,而传统较有影响的Yoshimura模型最大相对误差约是虚拟材料模型最大相对误差的4倍。对一些精密数控机床,在不允许直接对其进行损伤性实验时,使用本理论的预测方法,实验成本可以大大降低。 |
申请公布号 |
CN101832881A |
申请公布日期 |
2010.09.15 |
申请号 |
CN201010172821.5 |
申请日期 |
2010.05.14 |
申请人 |
华中科技大学 |
发明人 |
李斌;田红亮;毛新勇;毛宽民;刘红奇 |
分类号 |
G01M19/00(2006.01)I;G01N3/00(2006.01)I;G01N9/00(2006.01)I |
主分类号 |
G01M19/00(2006.01)I |
代理机构 |
华中科技大学专利中心 42201 |
代理人 |
曹葆青 |
主权项 |
一种机床固定结合部动态特性的检测方法,其特征在于,该方法包括下述步骤:第1步将机床上的各固定结合部分别视为一种等截面的虚拟材料,获取各虚拟材料的弹性模量、泊松比和密度;设构成固定结合部的二个零件分别为第一零件与第二零件,则该固定结合部所对应的虚拟材料的弹性模量、泊松比和密度的计算过程为:步骤(1.1):获取虚拟材料的弹性模量E和切变模量Gx:利用式①、式②计算得到虚拟材料的弹性模量:虚拟材料的弹性模量E函数表示为: <mrow> <mi>E</mi> <mrow> <mo>(</mo> <msub> <mi>A</mi> <mi>r</mi> </msub> <mo>></mo> <msub> <mi>A</mi> <mi>rc</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <mi>D</mi> </mrow> <mrow> <mn>3</mn> <msup> <mi>π</mi> <mn>2</mn> </msup> </mrow> </mfrac> <msup> <mi>E</mi> <mo>′</mo> </msup> <msup> <mi>G</mi> <mrow> <mn>1</mn> <mo>-</mo> <mi>D</mi> </mrow> </msup> <msubsup> <mi>a</mi> <mi>L</mi> <mrow> <mn>0.5</mn> <mi>D</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>a</mi> <mi>c</mi> <mrow> <mo>-</mo> <mn>0.5</mn> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>a</mi> <mi>L</mi> <mrow> <mo>-</mo> <mn>0.5</mn> </mrow> </msubsup> <mo>)</mo> </mrow> </mrow>①虚拟材料的法向载荷P函数表示为: <mrow> <mi>P</mi> <mrow> <mo>(</mo> <msub> <mi>A</mi> <mi>r</mi> </msub> <mo>></mo> <msub> <mi>A</mi> <mi>rc</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mfrac> <mrow> <mn>4</mn> <msqrt> <mi>π</mi> </msqrt> <mi>D</mi> </mrow> <mrow> <mn>3</mn> <mrow> <mo>(</mo> <mn>3</mn> <mo>-</mo> <mn>2</mn> <mi>D</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <msup> <mi>E</mi> <mo>′</mo> </msup> <msup> <mi>G</mi> <mrow> <mi>D</mi> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>a</mi> <mi>L</mi> <mrow> <mn>0.5</mn> <mi>D</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>a</mi> <mi>L</mi> <mrow> <mn>1.5</mn> <mo>-</mo> <mi>D</mi> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>a</mi> <mi>c</mi> <mrow> <mn>1.5</mn> <mo>-</mo> <mi>D</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mi>D</mi> <mrow> <mn>2</mn> <mo>-</mo> <mi>D</mi> </mrow> </mfrac> <mi>K</mi> <msub> <mi>σ</mi> <mi>y</mi> </msub> <msubsup> <mi>a</mi> <mi>L</mi> <mrow> <mn>0.5</mn> <mi>D</mi> </mrow> </msubsup> <msubsup> <mi>a</mi> <mi>c</mi> <mrow> <mn>1</mn> <mo>-</mo> <mn>1.5</mn> <mi>D</mi> </mrow> </msubsup> <mi>forD</mi> <mo>≠</mo> <mn>1.5</mn> </mtd> </mtr> <mtr> <mtd> <msqrt> <mi>πG</mi> </msqrt> <msup> <mi>E</mi> <mo>′</mo> </msup> <msubsup> <mi>a</mi> <mi>L</mi> <mn>0.75</mn> </msubsup> <mi>ln</mi> <mfrac> <msub> <mi>a</mi> <mi>L</mi> </msub> <msub> <mi>a</mi> <mi>c</mi> </msub> </mfrac> <mo>+</mo> <mn>3</mn> <mi>K</mi> <msub> <mi>σ</mi> <mi>y</mi> </msub> <msubsup> <mi>a</mi> <mi>L</mi> <mn>0.75</mn> </msubsup> <msubsup> <mi>a</mi> <mi>c</mi> <mn>0.25</mn> </msubsup> <mi>forD</mi> <mo>=</mo> <mn>1.5</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>②其中,Ar表示第一、第二零件实际接触面积,Are表示第一、第二零件实际临界接触面积,D表示第一零件或第二零件的表面轮廓分形维数,E′表示第一、第二零件相接触的等效弹性模量,G表示第一零件或第二零件的分形粗糙度参数,aL表示第一、第二零件接触面上最大弹性微接触点的面积,ae表示第一、第二零件接触面划分弹、塑性区域的临界面积;K表示相关系数,σy表示第一零件和第二零件中硬度较低材料的屈服强度;虚拟材料的切变模量Gx为 <mrow> <msub> <mi>G</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>A</mi> <mi>r</mi> </msub> <mo>></mo> <msub> <mi>A</mi> <mi>rc</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>16</mn> <mi>π</mi> </mfrac> <msup> <mrow> <mroot> <mrow> <mn>1</mn> <mo>-</mo> <mfrac> <mi>β</mi> <mi>f</mi> </mfrac> </mrow> <mn>3</mn> </mroot> <mi>G</mi> </mrow> <mo>′</mo> </msup> <mo>[</mo> <msup> <mrow> <mo>(</mo> <mfrac> <msub> <mi>a</mi> <mi>L</mi> </msub> <msub> <mi>a</mi> <mi>c</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mrow> <mn>0.5</mn> <mi>D</mi> </mrow> </msup> <mo>-</mo> <mn>1</mn> <mo>]</mo> </mrow>③其中,β表示切向载荷与法向载荷之比,f表示摩擦因数,G′表示两接触粗糙表面的等效切变模量步骤(1.2):根据式④,计算得到虚拟材料的泊松比υ: <mrow> <mi>υ</mi> <mo>=</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <msup> <mi>μ</mi> <mo>′</mo> </msup> <mo>)</mo> </mrow> <msup> <mi>E</mi> <mo>*</mo> </msup> </mrow> <msubsup> <mi>G</mi> <mi>x</mi> <mo>*</mo> </msubsup> </mfrac> <mo>-</mo> <mn>1</mn> </mrow>④其中,E*表示相应无量纲的虚拟材料的弹性模量,采用式⑤计算: <mrow> <mi>E</mi> <mo>*</mo> <mrow> <mo>(</mo> <msubsup> <mi>A</mi> <mi>r</mi> <mo>*</mo> </msubsup> <mo>></mo> <msubsup> <mi>A</mi> <mi>rc</mi> <mo>*</mo> </msubsup> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <mi>D</mi> </mrow> <mrow> <mn>3</mn> <msup> <mi>π</mi> <mn>2</mn> </msup> </mrow> </mfrac> <msup> <mrow> <msup> <mrow> <mi>G</mi> <mo>*</mo> </mrow> <mrow> <mn>1</mn> <mo>-</mo> <mi>D</mi> </mrow> </msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>2</mn> <mo>-</mo> <mi>D</mi> </mrow> <mi>D</mi> </mfrac> <msubsup> <mi>A</mi> <mi>r</mi> <mo>*</mo> </msubsup> <mo>)</mo> </mrow> </mrow> <mrow> <mn>0.5</mn> <mi>D</mi> </mrow> </msup> <mo>[</mo> <msubsup> <mi>a</mi> <mi>c</mi> <mrow> <mo>*</mo> <mo>-</mo> <mn>0.5</mn> </mrow> </msubsup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>2</mn> <mo>-</mo> <mi>D</mi> </mrow> <mi>D</mi> </mfrac> <msubsup> <mi>A</mi> <mi>r</mi> <mo>*</mo> </msubsup> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>0.5</mn> </mrow> </msup> <mo>]</mo> </mrow>⑤Gx*表示相应的无量纲虚拟材料的切变模量,采用式⑥计算: <mrow> <msubsup> <mi>G</mi> <mi>x</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>A</mi> <mi>r</mi> <mo>*</mo> </msubsup> <mo>></mo> <msubsup> <mi>A</mi> <mi>rc</mi> <mo>*</mo> </msubsup> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>16</mn> <mi>π</mi> </mfrac> <mroot> <mrow> <mn>1</mn> <mo>-</mo> <mfrac> <mi>β</mi> <mi>f</mi> </mfrac> </mrow> <mn>3</mn> </mroot> <mo>[</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>2</mn> <mo>-</mo> <mi>D</mi> </mrow> <mi>D</mi> </mfrac> <mfrac> <msubsup> <mi>A</mi> <mi>r</mi> <mo>*</mo> </msubsup> <msubsup> <mi>a</mi> <mi>c</mi> <mo>*</mo> </msubsup> </mfrac> <mo>)</mo> </mrow> <mrow> <mn>0.5</mn> <mi>D</mi> </mrow> </msup> <mo>-</mo> <mn>1</mn> <mo>]</mo> </mrow>⑥其中,Ar*表示无量纲实际接触面积,Are*表示无量纲实际临界接触面积,ac*表示无量纲划分弹、塑性区域的临界面积;步骤(1.3):根据公式⑦计算得到虚拟材料的密度。虚拟材料的密度ρ为 <mrow> <mi>ρ</mi> <mo>=</mo> <mfrac> <mrow> <msub> <mi>ρ</mi> <mn>1</mn> </msub> <msub> <mi>h</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>ρ</mi> <mn>2</mn> </msub> <msub> <mi>h</mi> <mn>2</mn> </msub> </mrow> <mrow> <msub> <mi>h</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>h</mi> <mn>2</mn> </msub> </mrow> </mfrac> </mrow>⑦ρ1、ρ2分别表示第一、第二零件的密度,h1、h2分别表示第一、第二零件的法向高度;第2步将各虚拟材料的弹性模量、泊松比和密度分别输入到有限元软件中,得到各虚拟材料所对应的固定结合部的有限元模型;第3步利用各个固定结合部的有限元模型,按照各零件间刚性连接的方法计算得到包含机床结合部影响的机床整机动态特性。 |
地址 |
430074 湖北省武汉市洪山区珞喻路1037号 |