CN103413049B - 基于机电耦合特性的并联机床结构优化参数值的获取方法 - Google Patents
基于机电耦合特性的并联机床结构优化参数值的获取方法 Download PDFInfo
- Publication number
- CN103413049B CN103413049B CN201310363227.8A CN201310363227A CN103413049B CN 103413049 B CN103413049 B CN 103413049B CN 201310363227 A CN201310363227 A CN 201310363227A CN 103413049 B CN103413049 B CN 103413049B
- Authority
- CN
- China
- Prior art keywords
- machine tool
- connecting rod
- performance
- sectional area
- cross
- 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
- 230000008878 coupling Effects 0.000 title claims abstract description 24
- 238000010168 coupling process Methods 0.000 title claims abstract description 24
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000005457 optimization Methods 0.000 claims description 24
- 238000003754 machining Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000004088 simulation Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 19
- 238000011156 evaluation Methods 0.000 abstract description 8
- 238000004519 manufacturing process Methods 0.000 abstract description 2
- 238000012545 processing Methods 0.000 abstract description 2
- 230000007246 mechanism Effects 0.000 description 19
- 230000003993 interaction Effects 0.000 description 9
- 230000008859 change Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000036544 posture Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/80—Technologies aiming to reduce greenhouse gasses emissions common to all road transportation technologies
- Y02T10/82—Elements for improving aerodynamics
Landscapes
- Automatic Control Of Machine Tools (AREA)
Abstract
本发明公开了一种基于机电耦合特性的并联机床结构优化参数值的获取方法,属于机床制造领域,该方法包括:首先选定杆长l和连杆横截面积s作为优化参数,建立机床机电耦合控制模型;将平均条件数和平均刚度作为机床的运动学性能评价指标,将一阶固有频率和轮廓误差分别作为动态性能评价指标及加工性能评价指标,得到杆长l和连杆横截面积s的最优取值范围;根据这些参数对机床性能校核,以满足机床的性能要求。本发明将机电耦合动态特性应用于并联机床的优化设计中,建立了一套基于机电耦合特性的并联机床结构优化参数值的获取方法,可以提高机床的运动控制精度。
Description
技术领域
本发明属于机床制造领域,尤其涉及于机电耦合特性的并联机床优化设计方法。
背景技术
在过去的近二十年里,并联机构由于具有高刚度、高承载能力、低运动惯量等优点在工业界得到越来越广泛的应用,尤其是在机床领域得到了广泛关注。然而,目前只有Z3主轴头以及Tricept并联机床等取得了瞩目的成功,这与人们对并联机床的期望值相距甚远。虽然造成这一现状的原因众多,但是并联机床理论上固有的高速、高加速度等动态特性在实际应用中没有体现出来是限制并联机床广泛应用的一个重要原因。
因此,研究人员分别采用结构动力学及弹性动力学方法研究并联机床动态特性。然而,无论利用结构动力学还是弹性动力学研究并联机床的动态特性,都主要是针对机械***。机械***与伺服***的特性研究往往是分开孤立研究。如孤立的研究机械***本身的动态特性,偏重静态的定位精度,对伺服驱动***的研究,只重视控制***本身的特性,对数控加工刀位轨迹的规划也只考虑几何学问题。这种研究方法对于常速和常规精度是适用的,但是难以满足数控机床的高速高精度要求。并联机床的运动特性和动态特性不仅与机床本身的结构参数有关,而且与伺服电机的电磁参数、特性参数以及动态性能有关。作为驱动单元的电机本身就是一个复杂的动力***,电机的输出是非线性的变速运动,电机输出运动的变化会影响机构动态性能,而机构运动参数的改变也会影响电机自身的运行状态,这种交互影响使电机与其拖动的机构之间产生耦合振动问题。因此,并联机床动态性能与机电耦合特性直接相关。
为了实现高动态特性,必须最终将机电耦合动态特性应用于并联机床优化设计中。然而,纵观并联机床优化设计中流行的性能评价指标,不难发现,这些研究多从运动学(以速度为主)角度出发,以雅可比矩阵或者其衍生矩阵为特征参数,定义性能度量指标。速度是衡量位移变化率的特征量,只能说明物体的运动状态而不是其它。并联机床是一个典型的复杂机电***,***性能不仅依赖机械***设计,同时和控制***性能息息相关。现有的并联机床优化设计方法完全孤立的设计其机械***和电气控制***,不考虑机电耦合效应,从而很难设计出性能优异的并联机床。
发明内容
本发明的目的是为克服已有技术的不足之处,提出一种基于机电耦合特性的并联机床结构优化参数值的获取方法,把机械***性能和控制***性能一起纳入设计中,有效提高了整体设计的效率,节约了设计成本,同时可以提高机床的运动控制精度。
本发明的基于机电耦合特性的并联机床结构优化参数值的获取方法采用如下步骤:
(1)根据机床结构以及雅可比矩阵选定杆长l和连杆横截面积s作为优化参数;
(2)将平均条件数和平均刚度作为机床的运动学性能指标,根据平均条件数和平均刚度分别与杆长l和连杆横截面积s的相互影响关系,分别画出相互影响曲线图,得到最小的平均条件数和最大的平均刚度值;
(3)将一阶固有频率作为动态性能指标,根据一阶固有频率与杆长l和连杆横截面积s的相互影响关系,画出相互影响曲线图,得到最大的一阶固有频率;
(4)根据机床结构和机床的控制***模型,建立机床机电耦合控制模型;
(5)基于机电耦合模型,将工件轮廓误差作为机床加工性能指标,根据该指标与杆长l和连杆横截面积s的相互影响关系,画出相互影响曲线图,得到最小轮廓误差;
(6)根据机床的性能要求对平均条件数、平均刚度、一阶固有频率最大值和最小轮廓误差加权,得到杆长l和连杆横截面积s的最优取值范围;
(7)从最优取值范围中选择一组杆长l和连杆横截面积s值,根据这些参数对机床性能校核,若满足机床的性能要求,则该组杆长l和连杆横截面积s值作为优化值,否则重复步骤(6)直到满足机床的性能要求。
本发明提供一种在并联机床的初始设计阶段,基于机电耦合特性的结构优化参数值的获取方法,其特点和有益效果是:
1.与现有的并联机床优化设计方法相比,在设计的初始阶段,把控制***性能对于并联机床***性能的影响考虑进去,将机械***性能和控制***性能同时考虑,有效提高了整体设计的效率,节约了设计成本。
2.在并联机床的初始设计阶段,就可以准确的仿真出机床的运动控制性能。
3.本方法可以应用在关于机构参数的计算机辅助设计软件开发中。
附图说明
图1为本发明的基于机电耦合特性的并联机床结构优化参数值的获取方法流程框图。
图2为本实施例中平均条件数与杆长和连杆横截面积相互影响曲线图。
图3为本实施例中平均刚度与杆长和连杆横截面积相互影响曲线图。
图4为本实施例中一阶固有频率与杆长和截面积相互影响曲线图。
图5为本实施例中并联机床机电耦合模式。
图6为本实施例中轮廓误差与杆长和连杆横截面积相互影响曲线图。
图7为本实施例中优化目标函数与杆长和连杆横截面积相互影响曲线图。
图8为本实施例中优化前后并联机床各向刚度在空间中的分布情况。
图9为本实施例中优化前后轮廓误差随时间变化的曲线。
具体实施方式
以下结合附图和具体实例对本发明作进一步详细说明。
将本发明所公布的一种基于机电耦合特性的并联机床结构优化参数值的获取方法应用于一台两自由度并联机床中。该并联机床由2-P(Pa)型并联机构构成,即每条运动支链为移动副-平行四边形结构,由交流伺服电动机驱动,其刀架固联于动平台。
本实施例的实现流程如图1所示,具体步骤如下:
1)根据机床结构以及雅可比矩阵选定杆长l和连杆横截面积s作为优化参数;
本实施例中综合考虑其运动学和动力学性能,选取的性能评价指标有:平均条件数、平均刚度一阶固有频率f和***轮廓误差ε。为了论述结构参数获取方法,仅选择杆长l和连杆横截面积s作为优化参数。其它参数值采用传统的优化方法获得:两个立柱之间的宽度D=4.284m,动平台高度h=0.8m,动平台宽度d=0.7m。
2)将平均条件数和平均刚度作为机床的运动学性能指标,根据平均条件数和平均刚度分别与杆长l和连杆横截面积s的相互影响关系,画出相互影响曲线图,得到最小的平均条件数和最大的平均刚度值;具体实施步骤如下:
21)平均条件数
并联机构雅可比矩阵J的条件数定义为:
其中,σmin和σmax分别是雅可比矩阵的最小和最大奇异值。
并联机构的灵巧度是以平均条件数描述并联机构整体灵活性。当并联机构处于不同的运动姿态时,并联机构的条件数是不同的,因此为了便于描述机构在整个工作空间中的整体灵活度,定义一个运动学全局的平均条件数:
式子(2)中,W代表了并联机构的工作空间。根据机床结构和不同的运动姿态下的平均条件数与杆长l和截面积s的相互影响关系,画出相互影响曲线图,如图2所示。其中,最小的平均条件数为
22)平均刚度
并联机构的刚度一定程度上反映了机构在加工过程中的变形难易程度,在并联机床的设计中,刚度是很重要的一个性能指标。当并联机构处于不同位姿时,其各个方向上的刚度是不同的,因此,为了描述整体的刚度性能,定义一个全局的平均刚度,其数学定义式为:
其中,K=[KxKyKz]T。
以平均刚度作为描述并联机构整体刚度性能的评价指标。根据机构的基本参数和各个方向上的刚度与杆长l和截面积s的相互影响关系,画出相互影响曲线图,如图3所示。其中,最优的刚度,即最大的平均刚度为:
3)将一阶固有频率作为动态性能指标,根据一阶固有频率与杆长l和连杆横截面积s的相互影响关系,画出相互影响曲线图,得到最大的一阶固有频率;
一阶固有频率是并联机床运行时的一个安全动态性能指标,当并联机床处于不同位姿时,相应的固有频率的值也是不断变化的。因此,为了考虑全局的动态性能,以工作空间中最小的一阶固有频率作为评价指标,其必须大于机床的正常工作频率。根据一阶固有频率f与杆长l和杆件横截面积s相互影响关系,画出相互影响曲线图,如图4所示。最优的固有频率,即最大的一阶固有频率为fmax=31.02Hz。
4)根据机床结构和机床的控制***模型,建立机床机电耦合控制模式;
并联机床基于模型的控制通常采用动力学前馈控制方式,对各支链的电机、电流环、速度环、位置环等分别建模,并通过仿真选取合适的PID参数。计算前馈通道的传递函数,并构建相应的前馈子***。如图5所示,该***前馈环由动力学模型和动力学前馈控制器组成;根据加工零件的外形尺寸确定虚轴空间期望轨迹,通过运动学逆解算出实轴空间期望轨迹。运动学控制子***根据实轴的期望轨迹,通过位置环、速度环和电流环,控制电机位置输出,同时,动力学控制子***根据实轴空间期望轨迹进行动力学逆解的运算,通过动力学前馈控制器补偿动态特性对机床精度影响,提高机床运动精度。
5)基于机电耦合模式,将轮廓误差作为机床加工性能指标,根据该指标与杆长l和连杆横截面积s的相互影响关系,画出相互影响曲线图,得到最大轮廓误差;
以***轮廓误差ε作为描述并联机床加工性能的评价指标,选取典型的运动轨迹,分别为直线运动和圆周运动。其中,直线运动轨迹的速度规划分别采用梯形速度规划、S形速度规划和正弦速度规划,圆周运动轨迹速度规划为正弦速度规划。画出相应的轮廓误差ε与杆长l和连杆横截面积s之间的关系,如图6所示。其中,四种运动轨迹(Vmax为最大速度)对应的最优轮廓误差,即最大轮廓误差的最小值分别为:ε1min=0.2269μm,ε2min=0.3288μm,ε3min=0.7421μm和ε4min=1.3223μm。
6)根据机床的性能要求对平均条件数、平均刚度、一阶固有频率和最小轮廓误差加权,得到杆长l和连杆横截面积s的最优取值范围;
综合考虑运动学性能指标、动态性能指标和加工性能指标,得到符合性能要求设计参数的最优区域,如图7所示;选取平均刚度平均条件数倒数固有频率f、***轮廓误差倒数作为相应的单目标优化函数。在本实施例中,选取的一种加权方案如表1所示。
机械***性能和控制***性能权重分别占整个***的50%,平均刚度、平均条件数导数和固有频率的权重分别占机械***权重的1/3。
表1
因此,***机电耦合性能优化目标函数可以表示为:
F(l,s)=∑[Fk(l,s)+Ff(l,s)+Fκ(l,s)+Fε(l,s)](4)
其中,
分别为x、y和z方向平均刚度的最大值,fmax为一阶固有频率的最大值,为平均条件数的最小值,ε1min、ε2min、ε3min和ε4min为四种运动轨迹最大轮廓误差的最小值。
如果设定当F(l,s)≥1.7时,***的综合性能较佳,则可以依据上述仿真结果得到相应杆长和连杆横截面积的选取范围:l∈[2.95m,3.2m],s∈[0.02m2,0.04m2]。
7)根据机床性能要求,从最优取值范围中选择一组杆长l和连杆横截面积s值;基于这一组参数,对机床性能进行仿真校核,若满足机床的性能要求,则该组杆长l和连杆横截面积s值作为优化值,否则重复步骤6)直到满足机床的性能要求
根据设计参数,对并联机床的性能进行校核和综合分析,选取该范围内的一组参数,分析并联机床优化前后的各项性能变化,优化前连杆长度l=3.35m,连杆截面参数s=0.08m2;优化后的参数连杆长度l=3.15m,连杆截面参数s=0.04m2。图8是优化前后并联机床各向刚度在空间中的分布情况,能够明显的看出刚度有所提升。图9是优化前后轮廓误差随时间变化的曲线,能够看出优化前的误差比优化后的误差要大,说明此优化方法能够提高机床的整体性能。
Claims (1)
1.一种基于机电耦合特性的并联机床结构优化参数值的获取方法,其特征在于,包括如下步骤:
(1)根据机床结构以及雅可比矩阵选定杆长l和连杆横截面积s作为优化参数;
(2)将平均条件数和平均刚度作为机床的运动学性能指标,根据平均条件数和平均刚度分别与杆长l和连杆横截面积s的相互影响关系,分别画出相互影响曲线图,得到最小的平均条件数和最大的平均刚度值;
(3)将一阶固有频率作为动态性能指标,根据一阶固有频率与杆长l和连杆横截面积s的相互影响关系,画出相互影响曲线图,得到最大的一阶固有频率;
(4)根据机床结构和机床的控制***模型,建立机床机电耦合控制模式;具体包括:并联机床基于模型的控制采用动力学前馈控制方式,对各支链的电机、电流环、速度环、位置环等分别建模,并通过仿真选取合适的PID参数;计算前馈通道的传递函数,并构建相应的前馈子***;该***前馈环由动力学模型和动力学前馈控制器组成;
(5)基于机电耦合模型,将轮廓误差作为机床加工性能指标,根据该指标与杆长l和连杆横截面积s的相互影响关系,画出相互影响曲线图,得到最大轮廓误差的最小值;
(6)根据机床的性能要求对最小平均条件数、最大平均刚度、一阶固有频率最大值和最大轮廓误差最小值加权,得到杆长l和连杆横截面积s的最优取值范围;
(7)从最优取值范围中选择一组杆长l和连杆横截面积s值,根据这些参数对机床性能校核,若满足机床的性能要求,则该组杆长l和连杆横截面积s值作为优化值,否则重复步骤(6)直到满足机床的性能要求。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310363227.8A CN103413049B (zh) | 2013-08-20 | 2013-08-20 | 基于机电耦合特性的并联机床结构优化参数值的获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310363227.8A CN103413049B (zh) | 2013-08-20 | 2013-08-20 | 基于机电耦合特性的并联机床结构优化参数值的获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103413049A CN103413049A (zh) | 2013-11-27 |
CN103413049B true CN103413049B (zh) | 2016-04-06 |
Family
ID=49606060
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310363227.8A Active CN103413049B (zh) | 2013-08-20 | 2013-08-20 | 基于机电耦合特性的并联机床结构优化参数值的获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103413049B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104898568B (zh) * | 2015-05-20 | 2018-01-19 | 西安交通大学 | 基于刚度辨识的数控机床进给***控制参数优化方法 |
CN105892294B (zh) * | 2016-04-14 | 2018-08-28 | 清华大学 | 一种并联构型主轴头的伺服***控制参数设计方法 |
CN108303244A (zh) * | 2018-02-28 | 2018-07-20 | 吉林大学 | 一种数控转塔刀架可靠性常速与加速试验装置及方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102615514B (zh) * | 2012-04-10 | 2013-08-14 | 清华大学 | 一种平面两自由度冗余驱动并联机构及采用该机构的机床 |
CN103217924B (zh) * | 2013-04-15 | 2015-05-20 | 清华大学 | 面向实时控制应用的过约束重型并联机床动力学建模方法 |
-
2013
- 2013-08-20 CN CN201310363227.8A patent/CN103413049B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN103413049A (zh) | 2013-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Du et al. | An accurate adaptive NURBS curve interpolator with real-time flexible acceleration/deceleration control | |
US9423788B2 (en) | Automatic method for milling complex channel-shaped cavities | |
CN104615084B (zh) | 加工进给速度优化的刀轨曲线轮廓误差补偿方法 | |
CN112698625B (zh) | 一种面向联动轨迹误差预测的五轴机床数字孪生建模方法 | |
CN102591257A (zh) | 面向参数曲线刀具轨迹的数控***轮廓误差控制方法 | |
CN105243218A (zh) | 机床热误差精度转化和模型建立方法 | |
CN101870073B (zh) | 基于工艺***刚度特性的多轴数控加工刀具运动规划方法 | |
Huang et al. | A survey of modeling and control in ball screw feed-drive system | |
CN103413049B (zh) | 基于机电耦合特性的并联机床结构优化参数值的获取方法 | |
US11263373B2 (en) | Predictive machining for forging distortions | |
CN102081376A (zh) | 一种基于指令序列优化的加工负荷控制*** | |
CN105005673A (zh) | 一种面向机床正向设计的整机静刚度匹配设计方法 | |
CN113946922A (zh) | 一种五轴联动铣削过程的动力学集成建模和加工精度预测方法 | |
CN103869748A (zh) | 非圆曲面xy直驱加工廓形误差交叉耦合控制***及方法 | |
CN110402188B (zh) | 驱动***和其评估 | |
Kim et al. | Hardware-in-the-loop simulation for estimation of position control performance of machine tool feed drive | |
Hashim et al. | Simulation of PSO-PI Co ontroller of DC Motor in Micro--EDM System for Biomedical Application | |
US9927801B2 (en) | Automatic method for milling complex channel-shaped cavities via coupling flank-milling positions | |
CN102402198A (zh) | 一种多轴数控机床通用后置处理方法 | |
CN103309280B (zh) | 一种用于重型并联机床的双前馈控制*** | |
WO2022149569A1 (ja) | 工作機械を備える加工システム、加工システムにおけるパラメータを修正するパラメータの修正方法、加工プログラムを修正するプログラム修正システム、およびプログラムの修正方法 | |
CN111086646B (zh) | 一种飞机襟翼操纵***的状态监测方法及*** | |
Fang et al. | Theory research on NC instruction interpretation for trajectory trace | |
Zhou | Research on Digital Engineering Modeling System under Computer Virtual Reality Technology | |
Park et al. | Development of a servo-based broaching machine using virtual prototyping technology |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |