CN103323278A - 一种滑动导轨系统虚拟材料层参数的识别方法 - Google Patents
一种滑动导轨系统虚拟材料层参数的识别方法 Download PDFInfo
- Publication number
- CN103323278A CN103323278A CN2013102276666A CN201310227666A CN103323278A CN 103323278 A CN103323278 A CN 103323278A CN 2013102276666 A CN2013102276666 A CN 2013102276666A CN 201310227666 A CN201310227666 A CN 201310227666A CN 103323278 A CN103323278 A CN 103323278A
- Authority
- CN
- China
- Prior art keywords
- material layer
- virtual material
- guide rail
- rail system
- sliding guide
- 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.)
- Pending
Links
Images
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种滑动导轨系统虚拟材料层参数的识别方法。该方法将滑动导轨的法向和侧向结合部看作一种具有一定厚度的虚拟材料层,并采用模态试验和有限元目标优化技术相结合的方法,识别得到虚拟材料层参数(弹性模量、泊松比、剪切模量、虚拟材料层厚度和面积)。通过参数识别结果计算的固有频率与模态试验结果误差低于8%,表明该方法简单实用、准确可靠,可在实际工程中得到广泛应用。
Description
技术领域
本发明属于机床整机动态特性分析的技术领域,特别是一种滑动导轨系统虚拟材料层参数的识别方法。
背景技术
现代数控机床正朝着高精度、高速度、高智能的方向快速发展,这就迫切要求机床设计在图样阶段就能预测整机性能,动态特性则是评价数控机床整机性能的关键。滑动导轨以其优良的承载能力和抗振性能仍在现代机床,特别是重型数控机床中得到广泛的应用。显然,滑动结合部的动态特性必将影响着机床整机的动态特性。中国专利:滑动导轨结合部动态特性参数测试装置及其测试方法,申请号:201010622646.5和中国专利:滑动结合部动态特性测试装置及测试方法,申请号:CN201010617228.7均将滑动导轨结合部看作一组弹簧阻尼器,虽然弹簧阻尼器可以定性地表征结合部的接触刚度和阻尼,但在机床整机建模时,弹簧阻尼器的数量和分布并不能明确给定,这就会造成机床整机动力学分析中出现漏阶甚至是虚假频率,影响机床整机性能预测的准确性。中国专利:一种机床固定结合部动态特性的检测方法,公开号:CN101832881A提供了一种计算固定结合部虚拟材料层参数的理论方法,该方法依赖的接触表面分形维数和分形粗糙度系数均难以测量,难以提供准确的虚拟材料层参数。
由上可知,现有技术均无法准确描述滑动结合部的动态特性。
发明内容
本发明所解决的技术问题在于提供一种滑动导轨系统虚拟材料层参数的识别方法。
实现本发明目的的技术解决方案为:一种滑动导轨系统虚拟材料层参数的识别方法,包括以下步骤:
步骤1、建立系统坐标系,将数控机床进给系统中的滑动承导件的承载面的外法线方向定义为z向,沿承导件长度方向定义为x向,与z向、x向均垂直的方向定义为y向;
步骤2、构建独立的滑动导轨系统,将数控机床进给系统中的移动部件和滑动承导件分离出来,将移动部件和滑动承导件之间的z向结合部a和y向结合部b均视为具有厚度的虚拟材料层,所述的虚拟材料层参数包括弹性模量E、泊松比υ、剪切模量G以及虚拟材料层厚度t、面积A;所述虚拟材料层为正交各向异性材料,弹性模量E包括三个方向上的弹性模量Ex,Ey,Ez;泊松比υ包括引起三个方向应变的泊松比υxy,υyz,υxz;剪切模量G包括三个平面的剪切模量Gxy,Gyz,Gxz。
步骤3、对上述的滑动导轨系统进行锤击模态试验,得到滑动导轨系统的各阶固有频率及振型;
步骤4、建立含有虚拟材料层的滑动导轨系统的有限元模型,任意选取一组虚拟材料层参数值输入到有限元模型中,进行模态分析,得到与锤击模态试验振型相同的六阶固有频率f1、f2、f3、f4、f5、f6,上述锤击模态试验振型与有限元模态分析对应的六阶频率为F1、F2、F3、F4、F5、F6;
步骤5、采用目标优化方法识别滑动导轨系统的虚拟材料层参数,将虚拟材料层参数中的弹性模量E、泊松比υ、剪切模量G、虚拟材料层厚度t、面积A作为优化设计变量,将作为目标函数,其中fi为与锤击模态试验振型相同的固有频率值,Fi为与有限元模态分析对应的六阶频率值,对滑动导轨系统有限元模型进行参数优化设计,从而得到滑动导轨系统虚拟材料层参数。
本发明与现有技术相比,其显著优点是:(1)采用虚拟材料层参数表征滑动导轨结合部的动力学性能突破了传统弹簧阻尼器理论模型的束缚,使得滑动导轨结合部模型更加符合实际;(2)提出基于目标优化的虚拟材料层参数的识别方法,该方法能够准确、迅速地获取滑动导轨结合部的虚拟材料层参数;(3)该方法不仅适用于滑动导轨结合部,也适用于机床整机结构中各种复杂结合部的虚拟材料层参数识别。
附图说明
图1是本发明的含虚拟材料层滑动导轨系统结构图。
图2是本发明的法向和侧向虚拟材料层示意图。
图3是本发明锤击模态试验系统图。
图4是本发明锤击模态试验中测点布置图。
图中标号所代表的含义为:
1.移动部件2.镶条3.滑动承导件4.加速度传感器5.力锤6.电荷放大器7.数据采集器8.计算机a.法向虚拟材料层b.侧向虚拟材料层。
具体实施方式
结合图1至图4,本发明的一种滑动导轨系统的虚拟材料层参数的识别方法,包括下述4个步骤:
步骤1、建立系统坐标系,将数控机床进给系统中的滑动承导件3的承载面的外法线方向定义为z向,沿承导件长度方向定义为x向,与z向、x向均垂直的方向定义为y向;构建独立的滑动导轨系统;将数控机床进给系统中的移动部件1和滑动承导件3分离出来,其中移动部件1与滑动承导件3的材料均为HT200,两滑动承导件2镜像布置,并固定连接在一个质量很大的基础底座上。将移动部件1和滑动承导件3之间的z向结合部a或y向结合部b均视为具有一定厚度的虚拟材料层,所述的虚拟材料层参数包括弹性模量E、泊松比υ、剪切模量G以及虚拟材料层厚度t、面积A。所述虚拟材料层为正交各向异性材料,弹性模量E包括三个方向上的弹性模量Ex,Ey,Ez;泊松比υ包括引起三个方向应变的泊松比υxy,υyz,υxz;剪切模量G包括三个平面的剪切模量Gxy,Gyz,Gxz。
步骤2、对上述的滑动导轨系统进行锤击模态试验;建立锤击模态试验系统,该系统包括加速度传感器4、力锤5、电荷放大器6、数据采集器7以及安装有机械结构模态分析软件Macras的计算机8,进行锤击模态试验得到滑动导轨系统的各阶固有频率及振型。
步骤3、建立含有虚拟材料层的滑动导轨系统的有限元模型;将虚拟材料层参数中弹性模量、泊松比、剪切模量以及虚拟材料层厚度、面积,任意选取一组值输入到有限元模型中,进行模态分析,得到与锤击模态试验振型相同的六阶固有频率f1、f2、f3、f4、f5、f6,上述锤击模态试验振型对应的六阶频率为F1、F2、F3、F4、F5、F6。
步骤4、采用目标优化技术识别滑动导轨系统的虚拟材料层参数;将虚拟材料层参数中的弹性模量、泊松比、剪切模量以及虚拟材料层厚度、面积作为优化设计变量,将作为目标函数,对滑动导轨系统有限元模型进行参数优化设计,从而识别得到滑动导轨系统虚拟材料层参数Eax=0.01MPa,Eay=0.01MPa,Eaz=5.54MPa;υax=0.25,υay=0.25,υaz=0.241;Gax=0.01MPa,Gay=0.01MPa,Gaz=2.23MPa;Ebx=0.01MPa,Eby=2.34MPa,Ebz=0.01MPa;υax=0.25,υay=0.237,υaz=0.25;Gax=0.01MPa,Gay=0.945MPa,Gaz=0.01MPa;ta=tb=0.9mm;Aa=0.018m2,Ab=0.006m2。
目标优化识别后的有限元模态分析的六阶固有频率和锤击模态试验六阶频率的比较如表1所示:
表1
1阶 | 2阶 | 3阶 | 4阶 | 5阶 | 6阶 | |
f(Hz) | 221.7 | 231.5 | 288.7 | 377.6 | 384.6 | 521.5 |
F(Hz) | 205.6 | 216.7 | 272.1 | 368.2 | 390.1 | 534.3 |
误差(%) | 7.83 | 6.83 | 6.1 | 2.55 | -1.41 | -2.4 |
通过上面的具体试验例子,采用本发明中的方法实现了滑动导轨系统虚拟材料层参数的识别。
以上对本发明的具体实施实例进行了详细阐述,但本发明并不限制于以上描述的具体实施实例,其只是作为范例。任何等同修改和替代也都在本发明的范畴之中。因此,在不脱离本发明的精神和范围下所作出的均等变换和修改,都应涵盖在本发明的范围内。
Claims (2)
1.一种滑动导轨系统虚拟材料层参数的识别方法,其特征在于,包括以下步骤:
步骤1、建立系统坐标系,将数控机床进给系统中的滑动承导件[3]的承载面的外法线方向定义为z向,沿承导件长度方向定义为x向,与z向、x向均垂直的方向定义为y向;
步骤2、构建独立的滑动导轨系统,将数控机床进给系统中的移动部件[1]和滑动承导件[3]分离出来,将移动部件[1]和滑动承导件[3]之间的z向结合部[a]和y向结合部[b]均视为具有厚度的虚拟材料层,所述的虚拟材料层参数包括弹性模量E、泊松比υ、剪切模量G以及虚拟材料层厚度t、面积A;
步骤3、对上述的滑动导轨系统进行锤击模态试验,得到滑动导轨系统的各阶固有频率及振型;
步骤4、建立含有虚拟材料层的滑动导轨系统的有限元模型,任意选取一组虚拟材料层参数值输入到有限元模型中,进行模态分析,得到与锤击模态试验振型相同的六阶固有频率f1、f2、f3、f4、f5、f6,上述锤击模态试验振型与有限元模态分析对应的六阶频率为F1、F2、F3、F4、F5、F6;
2.根据权利要求1所述的滑动导轨系统虚拟材料层参数的识别方法,其特征在于,步骤2中的虚拟材料层为正交各向异性材料,弹性模量E包括三个方向上的弹性模量Ex,Ey,Ez;泊松比υ包括引起三个方向应变的泊松比υxy,υyz,υxz;剪切模量G包括三个平面的剪切模量Gxy,Gyz,Gxz。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013102276666A CN103323278A (zh) | 2013-06-07 | 2013-06-07 | 一种滑动导轨系统虚拟材料层参数的识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013102276666A CN103323278A (zh) | 2013-06-07 | 2013-06-07 | 一种滑动导轨系统虚拟材料层参数的识别方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103323278A true CN103323278A (zh) | 2013-09-25 |
Family
ID=49192166
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2013102276666A Pending CN103323278A (zh) | 2013-06-07 | 2013-06-07 | 一种滑动导轨系统虚拟材料层参数的识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103323278A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103900778A (zh) * | 2014-04-16 | 2014-07-02 | 南京理工大学 | 双驱滑块式滚滑复合导轨动态特性测试实验装置与方法 |
CN107844656A (zh) * | 2017-11-15 | 2018-03-27 | 江苏亚威机床股份有限公司 | 一种滚动直线导轨有限元模型建立方法 |
CN110348166A (zh) * | 2019-07-19 | 2019-10-18 | 辽宁工程技术大学 | 一种玄武岩纤维树脂混凝土结合面的虚拟材料参数可视化识别方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101832881A (zh) * | 2010-05-14 | 2010-09-15 | 华中科技大学 | 一种机床固定结合部动态特性的检测方法 |
CN102095562A (zh) * | 2010-12-31 | 2011-06-15 | 北京工业大学 | 滑动结合面动态特性测试装置及测试方法 |
CN102183363A (zh) * | 2010-12-30 | 2011-09-14 | 南京理工大学 | 滑动导轨结合面动态特性参数测试装置及其测试方法 |
-
2013
- 2013-06-07 CN CN2013102276666A patent/CN103323278A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101832881A (zh) * | 2010-05-14 | 2010-09-15 | 华中科技大学 | 一种机床固定结合部动态特性的检测方法 |
CN102183363A (zh) * | 2010-12-30 | 2011-09-14 | 南京理工大学 | 滑动导轨结合面动态特性参数测试装置及其测试方法 |
CN102095562A (zh) * | 2010-12-31 | 2011-06-15 | 北京工业大学 | 滑动结合面动态特性测试装置及测试方法 |
Non-Patent Citations (2)
Title |
---|
LINGLING XIA 等: "Analysis of Holistic Dynamic Characteristics of A Vertial Machining Center KVC1050N", 《IEEE论文》, 31 December 2011 (2011-12-31), pages 67 - 71 * |
杜圣雨: "直线滚动导轨结合部动力学模型参数识别研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》, no. 05, 15 May 2009 (2009-05-15) * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103900778A (zh) * | 2014-04-16 | 2014-07-02 | 南京理工大学 | 双驱滑块式滚滑复合导轨动态特性测试实验装置与方法 |
CN107844656A (zh) * | 2017-11-15 | 2018-03-27 | 江苏亚威机床股份有限公司 | 一种滚动直线导轨有限元模型建立方法 |
CN110348166A (zh) * | 2019-07-19 | 2019-10-18 | 辽宁工程技术大学 | 一种玄武岩纤维树脂混凝土结合面的虚拟材料参数可视化识别方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103792143B (zh) | 一种单轴拉伸全程真应力应变曲线的快速获取方法 | |
CN107742018B (zh) | 基于路面材料模量应力和应变依赖模型的沥青路面结构分析增量方法 | |
Chatterjee et al. | Numerical simulations of pipe–soil interaction during large lateral movements on clay | |
CN103630450A (zh) | 考虑疲劳-蠕变交互损伤作用的沥青混合料寿命预测方法 | |
CN103047959B (zh) | 一种面向精密装配的基于熵理论的平面形状误差检测方法 | |
CN103412994B (zh) | 一种确定高速风洞大型飞机模型缩比的方法 | |
CN103955604A (zh) | 一种含裂纹金属梯度材料剩余强度预测方法 | |
CN107016173B (zh) | 一种高速压力机底座动态特性的可靠性设计方法 | |
CN106991219A (zh) | 一种考虑三维分形的法向界面刚度预测方法 | |
CN111339703A (zh) | 大应变条件下材料硬化行为的虚拟预测方法 | |
CN102567578B (zh) | 航天器振动试验夹具评价系统 | |
CN103323278A (zh) | 一种滑动导轨系统虚拟材料层参数的识别方法 | |
CN109839493A (zh) | 地下工程岩体质量评价方法、装置、存储介质及电子设备 | |
CN106354955A (zh) | 一种基于磨机振动模态参数的滑动轴承刚度识别方法 | |
Zheng et al. | Arch-dam crack deformation monitoring hybrid model based on XFEM | |
CN102279126B (zh) | 测试与cae仿真相结合确定材料性能参数的方法 | |
CN114519284A (zh) | 一种基于数值模拟的台阶爆破岩石块度预测方法 | |
CN105864335A (zh) | 非端部接触式少片斜线型副簧根部厚度的设计方法 | |
CN113836625A (zh) | 一种基于能力谱法的输电塔抗震性能等级划分方法及系统 | |
CN103267621B (zh) | 一种基于滚动导轨系统虚拟材料层参数的识别方法 | |
CN105825008A (zh) | 非端部接触式少片变截面主副簧的副簧起作用载荷验算方法 | |
CN109518574B (zh) | 湿化作用下路基顶面当量回弹模量的确定方法 | |
Yao et al. | Calculation and Control Methods for Equivalent Resilient Modulus of Subgrade Based on Nonuniform Distribution of Stress | |
CN113587880A (zh) | 一种测定扩孔率的试验与计算方法 | |
CN107192447A (zh) | 一种搜索结构物振动频率的窄带搜峰法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20130925 |