CN113110295B - 一种基于键合图的高速电主轴热误差建模方法 - Google Patents
一种基于键合图的高速电主轴热误差建模方法 Download PDFInfo
- Publication number
- CN113110295B CN113110295B CN202110487216.5A CN202110487216A CN113110295B CN 113110295 B CN113110295 B CN 113110295B CN 202110487216 A CN202110487216 A CN 202110487216A CN 113110295 B CN113110295 B CN 113110295B
- Authority
- CN
- China
- Prior art keywords
- bonding diagram
- model
- temperature field
- bonding
- thermal error
- 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
Images
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/404—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 control arrangements for compensation, e.g. for backlash, overshoot, tool offset, tool wear, temperature, machine construction errors, load, inertia
-
- 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/35408—Calculate new position data from actual data to compensate for contour error
Landscapes
- Engineering & Computer Science (AREA)
- Human Computer Interaction (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Automatic Control Of Machine Tools (AREA)
- Turning (AREA)
Abstract
本发明公开了一种基于键合图的高速电主轴热误差建模方法,该方法包括:1、构建高速电主轴轴承温度场键合图模型;2、构建电动机定子温度场键合图模型;3、构建电动机转子温度场键合图模型;4、根据高速电主轴的具体结构,连接步骤1至步骤3所建立的温度场键合图模型,建立完整的高速电主轴热误差系统温度场键合图模型;5、建立高速电主轴的热误差数学模型;6、根据模型进行系统仿真、分析得出高速电主轴热误差大小。本发明采用键合图建模方法,方便连接、添加各个零部件,可直接对存在多种能量形式的高速电主轴系统温度场进行建模、分析,并解决了产生无法分析的奇异点的难题,提高了高速电主轴系统建模、仿真、分析的快速性及准确性。
Description
技术领域
本发明涉及高速电主轴热误差建模技术领域,特别涉及一种基于键合图的高速电主轴热误差建模方法。
背景技术
在机械加工领域中,高速电主轴不仅是高精密数控机床的核心功能部件之一,还是数控机床的主要热源。在对电主轴热误差进行分析时,对定转子以及轴承等零部件发热的单独分析,忽略零部件之间的热对流和热传导产生的综合效应,并且目前ANSYS等仿真、分析软件存在无法对奇异点进行分析的难题,很难得到一个精确的数学模型,分析结果误差较大。
键合图法以状态方程作为数学模型,适用于机械、电气等所组成的复杂多能源系统问题的研究,它与现代控制工程里各种方法求得的状态方程有所不同。现代控制工程的状态方程一般表现是某单一变量的不同形式,如位移、速度、加速度等(由高阶方程转化而来),当研究的参数越出这一范围时,如具有力、速度、流量、压力、扭矩等同时存在的系统,现代控制工程比较难以讲他们组合在一起,而键合图法就特别适合于描述这样的系统,它将多种物理量统一归纳为4种状态变量,即势变量(e)、流变量(f)、变位变量与动量变量,键合图元通过传递功率的键相互连接构成系统的模型,方便地囊括上述各种不相干的参数,这也是这种方法的一大特点。
发明内容
本发明的目的是提供一种基于键合图的高速电主轴热误差建模方法,以解决现有技术中存在的上述不足。键合图理论为高速电主轴的研究提供了一个新的研究手段,该理论将多种能量形式集成的多体机械系统的分析研究转化为简明而统一的能量表达及传递方式,适合于高速电主轴的建模、控制、故障诊断等方面的研究。
为达到上述目的,本发明提供的种基于键合图的高速电主轴热误差建模方法,该方法包括以下步骤:
步骤a、以能量守恒定律为原则,构建标准键合图元件,建立高速电主轴轴承温度场键合图模型;
步骤b、构建标准键合图元件,建立电动机定子温度场键合图模型;
步骤c、构建标准键合图元件,建立电动机转子温度场键合图模型;
步骤d、连接高速电主轴轴承温度场键合图模型、电动机定子温度场键合图模型、电动机转子温度场键合图模型,建立完整的高速电主轴热误差系统温度场键合图模型;
步骤e、建立高速电主轴的热误差数学模型,根据模型进行系统仿真、分析得出高速电主轴热误差大小。
进一步地,所述步骤a、b、c中,零部件温度场键合图分别由键合图基本元件势源、流源、阻性元件、容性元件、“0”结、“1”结构成。
进一步地,所述步骤d中,通过采用势源、流源、惯性元件、阻性元件、容性元件、“0”结、“1”结等标准键合图元件对高速电主轴轴承温度场键合图模型、电动机定子温度场键合图模型、电动机转子温度场键合图模型进行连接,建立完整的高速电主轴热误差系统温度场键合图模型。
进一步地,所述步骤e中,建立多能量形式共存的高速电主轴系统温度场的热误差数学模型,通过20Sim仿真模拟平台,根据模型进行系统仿真、分析得出高速电主轴热误差大小。
本发明与现有技术相比,具有以下积极效果:
1、本发明采用以能量守恒原理为基础的键合图方法,可直接对整个存在多种能量形式的高速电主轴系统温度场进行建模、分析,并解决了ANSYS等仿真、分析软件存在无法对奇异点进行分析的难题,提高了高速电主轴热误差建模、仿真、分析的快速性及可靠性、准确性;
2、本发明采用的键合图建模方法能够定量地描述高速电主轴系统各个结构之间的物理关系,能够方便地添加零部件;
3、本发明采用键合图方法进行建模,运算速度极快,可自行设定节点,对高速电主轴工作状态进行实时监测、预测结果,同时便于更换初始条件进行多次类比分析。
附图说明
图1为高速电主轴轴承温度场键合图;
图2为电动机定子温度场键合图;
图3为电动机转子温度场键合图;
图4为高速电主轴热误差系统温度场键合图;
图5为高速电主轴热伸长示意图。
具体实施方式
本发明提供一种基于键合图的高速电主轴建模方法,为了使本发明更加明显易懂,以下结合附图和具体实施方式对本发明做进一步说明。
本发明的高速电主轴热误差系统温度场主要包括轴承部分、定子部分和转子部分,该系统的建模方法是首先自行划分节点,分别建立这三个部分的键合图模型,然后根据三者之间的位置关系将所有键合图模型用标准键合图元件连接,从而构建完整的高速电主轴热误差系统温度场的键合图模型,进一步建立高速电主轴的热误差数学模型,根据模型进行系统仿真、分析得出高速电主轴热误差大小。主要步骤如下:
步骤a、建立高速电主轴轴承部分键合图模型:
如图1所示,以能量守恒定律为原则的键合图方法适合于机电等多能量域共存系统的统一建模,首先构建标准键合图元件,并在此基础上建立轴承键合图模型。
轴承包含三个部分:轴承内环、轴承外环、陶瓷球。其中,轴承外侧还设有起支撑固定作用的轴承滑套、轴承座,温度场键合图分别由键合图基本元件势源Se、流源Sf、阻性元件R、容性元件C、“0”结、“1”结构成。
步骤b、建立电动机定子部分键合图模型:
如图2所示,构建标准键合图元件,建立电动机定子温度场键合图模型,定子包含三个部分:定子铁芯、定子绕组、机座,温度场键合图分别由键合图基本元件势源Se、流源Sf、阻性元件R、容性元件C、“0”结、“1”结构成。
步骤c、建立电动机转子部分键合图模型:
如图3所示,构建标准键合图元件,建立电动机转子温度场键合图模型,转子包含三个部分:转子铁芯、转子绕组、转轴,温度场键合图分别由键合图基本元件势源Se、流源Sf、阻性元件R、容性元件C、“0”结、“1”结构成。
步骤d、建立高速电主轴热误差系统温度场键合图模型:
如图4所示,本发明还提供了一种采用如上文所述的基于键合图的高速电主轴热误差系统温度场,其特征在于,主要包含高速电主轴轴承部分、电动机定子部分和电动机转子部分;其中,通过建立高速电主轴轴承键合图模型、电动机定子键合图模型和电动机转子键合图模型,再利用标准键合图元件连接所述高速电主轴轴承键合图模型、所述电动机定子键合图模型,所述电动机转子键合图模型,建立完整的高速电主轴热误差系统温度场键合图模型。
步骤e、建立高速电主轴热误差数学模型,根据模型进行系统仿真、分析得出高速电主轴热误差大小:
(1)建立高速电主轴的热误差数学模型:
从热弹性力学角度分析,高速电主轴单元热误差的产生可以归结于电主轴结构各部件的热变形的叠加。即电主轴结构温度场随时间发生变化,继而造成电主轴单元结构各部件的热变形及刚度变化,最终以电主轴末端热误差的形式表现出来。而电主轴结构任何零部件的热变形都可简化为一维杆热变形ΔL进行分析,由于主轴的径向被轴承所约束,轴的热伸长不能自由发生,于是在轴内产生了热应力Σ,温升时为压应力、温降时为拉应力,所以有
其中:αR为热膨胀系数,1/℃;σ为一维杆压力,MPa;P为轴向力,N;ER为弹性模量,N/m2;j为轴向刚度,N/m;AR为横截面积,mm2。
上式可进而简化为式:
(2)建模零点选取:
如图5所示,建立一维坐标系,对高速电主轴前部的热伸长量进行局部分析可得
ΔX2=ΔL5+ΔL6+ΔL7+ΔL8+ΔL9
ΔX=ΔX2-ΔX1
其中,ΔL0,ΔL1,ΔL2,ΔL3,ΔL4分别代表了中轴承,轴承滑套,前轴承,前轴承座,前法兰以中轴承内环右端面为零点O1的热伸长量,ΔX1代表第一部分总伸长量。由于中轴承、轴承滑套、前轴承座因为温升出线热伸长的现象,会导致中轴承内环右端面于主轴刀柄安装面反向伸长,因此需设定以主轴前法兰安装面为零点O2,此面对于电主轴整体是固定的,以有限元分析结果进行划分,将温度相同的一段划分在一起,则出现ΔX5,ΔX6,ΔX7,ΔX8,ΔX9五段伸长量,ΔX2代表第二部分总伸长量,主轴的热伸长为两部分伸长量的整合结果,ΔX代表了两部分的相对伸长长度。
(3)仿真、分析得出高速电主轴热误差大小:
根据构建的完整的高速电主轴热误差系统温度场键合图模型,得到高速电主轴内部温度场分布情况,分别以高速电主轴在运行过程中,前、中、后轴承、电动机定转子温度的时域变化数据为输入,以高速电主轴热伸长实时补偿量为输出,结合20Sim仿真模拟平台的热伸长功能模块,进行系统仿真、分析得出高速电主轴热误差大小。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (1)
1.一种基于键合图的高速电主轴热误差建模方法,其特征在于,包括轴承部分、定子部分和转子部分,该建模方法是首先自行划分节点,分别建立这三个部分的键合图模型,然后根据三者之间的位置关系将所有键合图模型用标准键合图元件连接,从而构建完整的高速电主轴热误差系统温度场的键合图模型,进一步建立高速电主轴的热误差数学模型,根据模型进行系统仿真、分析得出高速电主轴热误差大小;
其主要步骤如下:
步骤a、建立高速电主轴轴承部分键合图模型:
以能量守恒定律为原则的键合图方法适合于机电多能量域共存系统的统一建模,首先构建标准键合图元件,并在此基础上建立轴承键合图模型;轴承包含三个部分:轴承内环、轴承外环、陶瓷球;其中,轴承外侧还设有起支撑固定作用的轴承滑套、轴承座,温度场键合图分别由键合图基本元件势源Se、流源Sf、阻性元件R、容性元件C、“0”结、“1”结构成;
步骤b、建立电动机定子部分键合图模型:
构建标准键合图元件,建立电动机定子温度场键合图模型,定子包含三个部分:定子铁芯、定子绕组、机座,温度场键合图分别由键合图基本元件势源Se、流源Sf、阻性元件R、容性元件C、“0”结、“1”结构成;
步骤c、建立电动机转子部分键合图模型:
构建标准键合图元件,建立电动机转子温度场键合图模型,转子包含三个部分:转子铁芯、转子绕组、转轴,温度场键合图分别由键合图基本元件势源Se、流源Sf、阻性元件R、容性元件C、“0”结、“1”结构成;
步骤d、建立高速电主轴热误差系统温度场键合图模型:
所述步骤d提供了一种基于键合图的高速电主轴热误差系统温度场模型,包括高速电主轴轴承部分、电动机定子部分和电动机转子部分;其中,通过建立高速电主轴轴承键合图模型、电动机定子键合图模型和电动机转子键合图模型,再利用标准键合图元件连接所述高速电主轴轴承键合图模型、所述电动机定子键合图模型、所述电动机转子键合图模型,建立完整的高速电主轴热误差系统温度场键合图模型;
步骤e、建立高速电主轴热误差数学模型,根据模型进行系统仿真、分析得出高速电主轴热误差大小:
建立多能量形式共存的高速电主轴的热误差数学模型,通过20Sim仿真模拟平台,根据模型进行系统仿真、分析得出高速电主轴热误差大小。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110487216.5A CN113110295B (zh) | 2021-05-05 | 2021-05-05 | 一种基于键合图的高速电主轴热误差建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110487216.5A CN113110295B (zh) | 2021-05-05 | 2021-05-05 | 一种基于键合图的高速电主轴热误差建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113110295A CN113110295A (zh) | 2021-07-13 |
CN113110295B true CN113110295B (zh) | 2023-03-17 |
Family
ID=76720826
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110487216.5A Active CN113110295B (zh) | 2021-05-05 | 2021-05-05 | 一种基于键合图的高速电主轴热误差建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113110295B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113721548B (zh) * | 2021-09-06 | 2022-05-03 | 哈尔滨理工大学 | 一种数控铣头热误差补偿方法及系统 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102658499B (zh) * | 2012-04-20 | 2014-08-06 | 西安交通大学 | 一种精密卧式加工中心主轴热误差补偿方法 |
CN103853889A (zh) * | 2014-03-12 | 2014-06-11 | 南京航空航天大学 | 一种基于键合图的电力电子系统故障诊断方法 |
CN109176524B (zh) * | 2018-09-30 | 2020-07-31 | 上海神添实业有限公司 | 一种基于键合图的移动机械臂系统及其建模方法 |
CN111523212B (zh) * | 2020-04-15 | 2023-05-26 | 长春科技学院 | 一种电主轴热误差预测模型的建立方法 |
CN112434452B (zh) * | 2020-11-11 | 2024-05-07 | 西安交通大学 | 一种高速电主轴轴承轴径配合公差确定方法、系统设备 |
-
2021
- 2021-05-05 CN CN202110487216.5A patent/CN113110295B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113110295A (zh) | 2021-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109870134B (zh) | 一种旋转叶片非接触式动应变场测量方法及其系统 | |
CN106202647B (zh) | 电主轴的多轴疲劳寿命预测方法及疲劳寿命可靠性评估方法 | |
CN104200019B (zh) | 一种重型数控落地铣镗床热特性数值仿真方法 | |
CN110412462B (zh) | 一种航空用永磁同步电动机瞬态温度预测方法 | |
CN107111297A (zh) | 用于由至少一台计算机数控机器加工的工件的部件分析的计算机实现方法 | |
CN113110295B (zh) | 一种基于键合图的高速电主轴热误差建模方法 | |
CN108334661B (zh) | 一种进给轴热变形预测方法 | |
CN103366050A (zh) | 一种高速电主轴力-热耦合建模方法 | |
CN109684711B (zh) | 一种涡轴发动机气动连接多转子耦合振动分析方法 | |
CN109657320A (zh) | 一种风电机组传动链建模与扭振特性分析方法 | |
CN103487182B (zh) | 一种高速转子流阻扭矩的高精度测量方法及装置 | |
CN112364551B (zh) | 一种行星传动齿轮箱疲劳分析方法 | |
CN101628340B (zh) | 大型船用曲轴车床主轴箱轴承间距的优化方法 | |
CN117763760A (zh) | 一种发动机叶片及轮盘寿命计算方法和系统 | |
JP6605369B2 (ja) | シミュレーションプログラム | |
CN116629053A (zh) | 一种水轮机轴系湿模态的分析方法、系统、设备及介质 | |
CN113029620B (zh) | 轴-盘-叶片非轴对称旋转机械振动响应预测方法 | |
CN206756452U (zh) | 一种工业机器人动力单元可靠性试验台 | |
Feng et al. | Analysis of static and dynamic characteristic of spindle system and its structure optimization in camshaft grinding machine | |
TW201933009A (zh) | 設計支援裝置、設計支援方法、及設計支援程式 | |
Qu et al. | Dynamic design of grinding and polishing machine tool for blisk finishing | |
Liu et al. | Integrated design and analysis system for feed drive system of CNC machine tools | |
CN111086646A (zh) | 一种飞机襟翼操纵系统的状态监测方法及系统 | |
CN107942734B (zh) | 一种基于动态测试数据的进给系统闭环时域仿真方法 | |
CN112069611A (zh) | 一种分块建模和实验参数辨识的摆角铣头整机建模方法 |
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 |