CN111665785A - 一种六轴五联动激光加工开放式数控系统及其工作方法 - Google Patents

一种六轴五联动激光加工开放式数控系统及其工作方法 Download PDF

Info

Publication number
CN111665785A
CN111665785A CN202010544778.4A CN202010544778A CN111665785A CN 111665785 A CN111665785 A CN 111665785A CN 202010544778 A CN202010544778 A CN 202010544778A CN 111665785 A CN111665785 A CN 111665785A
Authority
CN
China
Prior art keywords
module
curve
speed
laser
axis
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.)
Granted
Application number
CN202010544778.4A
Other languages
English (en)
Other versions
CN111665785B (zh
Inventor
刘强
殷振朔
孙鹏鹏
王健
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangxi Research Institute Of Beijing University Of Aeronautics And Astronautics
Beihang University
Original Assignee
Jiangxi Research Institute Of Beijing University Of Aeronautics And Astronautics
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jiangxi Research Institute Of Beijing University Of Aeronautics And Astronautics, Beihang University filed Critical Jiangxi Research Institute Of Beijing University Of Aeronautics And Astronautics
Priority to CN202010544778.4A priority Critical patent/CN111665785B/zh
Publication of CN111665785A publication Critical patent/CN111665785A/zh
Application granted granted Critical
Publication of CN111665785B publication Critical patent/CN111665785B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical 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/404Numerical 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
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B23MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23KSOLDERING OR UNSOLDERING; WELDING; CLADDING OR PLATING BY SOLDERING OR WELDING; CUTTING BY APPLYING HEAT LOCALLY, e.g. FLAME CUTTING; WORKING BY LASER BEAM
    • B23K26/00Working by laser beam, e.g. welding, cutting or boring
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical 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/41Numerical 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 interpolation, e.g. the computation of intermediate points between programmed end points to define the path to be followed and the rate of travel along that path
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P70/00Climate change mitigation technologies in the production process for final industrial or consumer products
    • Y02P70/10Greenhouse gas [GHG] capture, material saving, heat recovery or other energy efficient measures, e.g. motor control, characterised by manufacturing processes, e.g. for rolling metal or metal working

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Optics & Photonics (AREA)
  • Plasma & Fusion (AREA)
  • Mechanical Engineering (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Numerical Control (AREA)

Abstract

本发明公开一种六轴五联动激光加工开放式数控系统及其工作方法,属于激光加工数控技术领域;其技术要点在于:所述的数控系统包括上位机和下位机;上位机包括四个模块:人机界面模块、译码器、全局速度规划模块、激光功率优化模块;下位机包括五个模块:插补计算模块、坐标变换模块、位置补偿模块、逻辑处理模块、伺服和激光控制模块。采用本申请的一种六轴五联动激光加工开放式数控系统及其工作方法,能够有效的提高激光加工质量。

Description

一种六轴五联动激光加工开放式数控系统及其工作方法
技术领域
本发明涉及激光加工数控技术领域,更具体地说,尤其涉及一种六轴五联动激光加工开放式数控系统及其工作方法。
背景技术
激光加工以其效率高,质量高,灵活性高,操作简便,节能环保等突出优点,被广泛应用于汽车,电子,航空航天,冶金,铁路和轮船等行业,已成为先进制造业的发展方向。随着激光加工技术的飞速发展,其对数控系统的要求越来越高。仅生产激光加工设备而不配备高质量数控系统不可能充分利用激光加工的优势。因此,开发一种高效高精、开放智能的激光加工数控系统是有必要的。
发明内容
本发明的目的在于针对上述现有技术的不足,提供一种六轴五联动激光加工开放式数控系统及其工作方法。
本发明解决其技术问题所采用的技术方案是:
一种六轴五联动激光加工开放式数控系统,包括上位机(1)和下位机(2);其特征在于:
上位机包括四个模块:人机界面模块(3)、译码器(4)、全局速度规划模块(5)、激光功率优化模块(6);
其中,人机界面模块(3)用于输入加工程序和显示各类信息;
其中,译码器(4)用于将输入的NC代码解释成指令队列;
其中,全局速度规划模块(5)根据机床动力学限制修改指令队列中的超约束值;
其中,激光功率优化模块(6)用于调整激光功率与进给速度相匹配;
下位机(2)包括五个模块:插补计算模块(7)、坐标变换模块(8)、位置补偿模块(9)、逻辑处理模块(10)、伺服和激光控制模块(11);
其中,插补计算模块(7)用来将轨迹曲线离散成位置坐标;
其中,坐标变换模块(8)用来进行六轴五联动正逆坐标变换;
其中,补偿模块(9)用来对位置误差进行补偿以提高定位精度;
其中,逻辑处理模块(10)用来处理外部输入IO信号;
其中,伺服和激光控制模块(11)用来驱动轴和激光器动作;
上位机各组成模块信息流程如下:
1)人机界面接收加工代码传至译码器;
2)译码器将加工代码转换成指令队列方便机器内部传送和处理,并将指令队列传至全局速度规划模块;
3)全局速度规划模块在接收指令队列后根据机床的动力学限制修改队列中的超约束值,并将修改后的指令队列传至激光功率优化模块;
4)激光功率优化模块接收指令队列后修改激光的功率使之与进给速度相匹配;至此上位机非实时任务完成,得到了适合加工的指令队列;
上位机处理好的指令队列通过ADS标准通信协议传至下位机;
下位机各组成模块信息流程如下:
1)上位机处理好的指令队列首先进入插补计算模块将指令队列中的曲线轨迹离散成位置坐标,并传至坐标变换模块;
2)坐标变换模块将位置坐标分解成各轴运动量传给补偿模块;
3)补偿模块根据激光干涉仪检测的机床补偿值将调整各轴运动量的值,以提高定位精度;
4)伺服和激光控制模块接收到各轴运动量和激光指令功率后利用TwinCAT PTP功能驱动伺服轴和激光器动作;
5)逻辑处理模块的任务是处理外部输入的IO信号,实时在线工作。
下位机接收到来自ADS协议的指令队列后按照指令进行插补计算,得到位置值后经过坐标变换计算和补偿计算即可驱动各轴运动,逻辑处理模块和伺服和激光控制模块可直接改变PLC状态。
进一步,人机界面模块(3)基于Qt平台开发,利用QML和C++语言开发系统前端和后端。
进一步,译码器(4)利用正则表达式进行文本匹配,提取指令信息,将指令信息存入指令队列供速度规划模块(5)使用。
全局速度规划模块(5),包括以下步骤:
步骤一:接收到来自译码器(4)的指令队列后首先根据机床联动动力学约束计算前后曲线段衔接点的过渡速度(即第i段曲线的末速度或者第i+1段曲线的起速度)vi为:
Figure BDA0002540236790000021
其中Amax为机床最大加速度,Tc为插补周期,F为指令速度,θ=arccos(τ12),τ1和τ2分别是第i曲线段的末端单位切矢和第i+1曲线段的起始单位切矢;
步骤二:计算第i段曲线段加(减)速所需最小距离Smin
其中当加(减)速过程中有匀加(减)速阶段的计算方式为:
Figure BDA0002540236790000031
反之,当加(减)速过程中没有匀加(减)速阶段的计算方式:
Figure BDA0002540236790000032
其中vmax=max(vi-1,vi),vmin=min(vi-1,vi),vi-1和vi分别为第i段曲线的起末速度,Am和Jm分别为机床所能达到的最大加速度和最大捷度;
步骤三:如果曲线段长度L>Smin,则起末速度可达,反之,则需根据L的大小调整起末速度;
步骤四:判断指令速度是否可达,如果不可达根据以下公式降低指令速度:
Figure BDA0002540236790000033
其中F为指令速度,Jm为机床所能达到的最大捷度;
完成上述步骤后,指令队列中的速度超约束值已经修改完毕,接下来修改激光功率使之与进给速度相匹配。
进一步,激光功率优化模块(6)具有以下流程:
步骤一:计算两相邻曲线段最大转角限制θcr
Figure BDA0002540236790000034
其中F为指令速度,ε为预设最大转角误差,Kf和Kp分别为伺服前馈增益和位置增益;
步骤二:计算前后两段曲线的弓高误差δ1、δ2
δ1=R(1-cosα1)
δ2=R(1-cos(π-θ-α1))
其中R为两曲线段外接圆半径,α1=arctan(l1sin(π-θ)/(l2+l1cos(π-θ))),l1和l2分别为前后两段曲线段的长度,θ为前后两曲线段的切矢夹角;
步骤三:如果前后两曲线段的切矢夹角θ≤θcr或者预设最大弓高误差δ<min(δ12),则根据衔接点速度修改拐点处的激光功率Pi为:
Pi=viPcmd/vcmd
其中vi为衔接点速度,vcmd和Pcmd分别为指令速度和指令激光功率;
完成上述步骤后,激光功率已经和进给速度相匹配,接下来通过ADS协议将指令队列传至下位机。
进一步,插补计算模块(7)具有以下流程:
步骤一:根据“S”型速度曲线计算插补步长ΔL为:
ΔL=vtTc
其中vt为速度曲线上t时刻的指令速度,Tc为插补周期;
步骤二:根据步长计算下一个周期节点矢量ui+1,其中直线、圆弧、NURBS曲线插补节点矢量计算分别如下:
直线的插补节点矢量递推计算:ui+1=ui+Tc
圆弧的插补节点矢量递推计算:ui+1=ui+ΔL/R
NURBS曲线的插补节点矢量递推计算:
Figure BDA0002540236790000041
其中,R为圆的半径值,C′(ui)和C″(ui)分别为NURBS曲线的一阶和二阶导数,||·||为向量的模,<·>为向量内积;ui是当前周期曲线在参数域上的的节点矢量,ui+1是下一周期曲线在参数域上的的节点矢量(节点矢量是参数曲线的重要参数,一个参数域的节点矢量对应曲线上一个点);
步骤三:将下一周期的节点矢量代入曲线表达式计算插补点坐标,其中直线、圆弧、NURBS曲线在ui+1处的坐标Pi+1计算分别如下:
直线的坐标计算:
Figure BDA0002540236790000042
其中,Pi+1是节点参数为ui+1时曲线上的点,Ps为起点坐标,Pe为终点坐标,us为起点的节点矢量,ue为终点的节点矢量(给定插补曲线后us、ue为已知量,默认:us=0,ue=1);
圆弧的坐标计算:Pi+1=Pc+R(ui+1-us,N)(Ps-Pc)
其中,R(ui+1-us,N)是表示向量绕N旋转ui+1-us角度,N为空间圆的法向量,Pc为圆心坐标,Ps起点坐标,us起点的节点矢量;
NURBS曲线的基函数递推计算如下:
Figure BDA0002540236790000043
其中i为节点序号,k为曲线次数,Ni,k(·)为基函数,u表示节点矢量(已知量,给定一条NURBS加工轨迹就会提供相应的节点矢量);
根据NURBS曲线的基函数计算在参数ui+1处的插补坐标Pi+1(ui+1)如下:
Figure BDA0002540236790000051
其中di为控制顶点,wi为权因子,
Figure BDA0002540236790000052
为基函数的累加;
完成上述步骤后,指令曲线已经离散成位置坐标值,接下来传至坐标变换模块分解成各轴运动量。
进一步,在测量出七个机床固有参数后,坐标变换模块(9)能够适应AB轴结构形式的正逆坐标变换;完成坐标变换后,可得到各轴的运动量,接下来传至补偿模块进行调整。
进一步,经过坐标变换后的各轴位置在补偿模块(10)中进行双向补偿,有效提高机床定位精度。
进一步,伺服和激光控制模块(12)利用TwinCAT提供的Tc2-MC2库,利用外部使能功能实现点对点运动控制。
一种六轴五联动激光加工开放式数控系统的工作方法,其特征在于,采用前述的一种六轴五联动激光加工开放式数控系统,包括以下步骤:
系统部署到机床的工控机内,系统中各轴变量、激光变量以及IO等变量分别和机床的伺服驱动器、激光器IO模块等硬件建立连接;
机床开机自动运行系统程序,人机交互模块接收加工程序和显示系统状态信息,接收的加工程序经过系统各个模块的处理优化后,驱动硬件完成加工。
进一步,译码器利用正则表达式进行文本匹配,提取指令信息,将指令信息存入指令队列供速度规划模块使用。
进一步,速度规划模块具有以下步骤:接收指令队列后首先根据机床联动动力学约束计算衔接点的过渡速度,使速度曲线更加光顺;之后判断过渡速度是否可达,如果不可达则调整过渡速度;最后判断指令速度是否可达,如果不可达则降低指令速度。
进一步,激光功率优化模块具有以下流程:首先根据最大转角限制和弓高误差约束判断路径拐点,然后根据拐点速度修改拐点处的激光功率,以此减少拐角处的烧蚀现象。
进一步,插补计算模块根据“S”型速度曲线计算插补步长,根据步长和曲线表达式计算下一个周期的位置点,插补计算模块能够完成直线、圆弧、NURBS三种类型曲线插补。
进一步,坐标变换模块在测量出七个机床固有参数后,能够适应AB轴结构形式的正逆五轴坐标变换。
进一步,位置补偿模块能够进行双向补偿,有效提高机床定位精度。
进一步,伺服和激光控制模块利用TwinCAT提供的Tc2-MC2库,利用外部使能功能实现点对点运动控制。
本发明的有益效果是:
第一,一种六轴五联动激光加工开放式数控系统能够准确的将NC代码翻译成指令队列,能够在机床动力学约束下合理的进行全局速度规划,能够准确的进行插补计算,能够实现六轴五联动的坐标变换,并且系统具有良好的开放性和通用性,可方便的扩展和移植。
第二,本发明的六轴五联动激光加工开放式数控系统包括上位机和下位机两部分。上位机是基于Qt平台设计的,执行非实时任务,运行于Windows系统。Qt是一个用于UI设计的高效强大的软件,它支持不同的开发语言,例如C++,QML,Python等。下位机在TwinCAT(The Windows Control and Automation Technology)内核中运行,处理实时任务。TwinCAT是德国倍福开发的基于PC的自动化软件,编程遵循IEC 61131-3标准,同时可方便的兼容C/C++,Matlab/Simulink等语言开发的模块。因此,TwinCAT本身具有开放性、灵活性和模块编程的特点,便于开发复杂的开放式数控系统。上位机和下位机的通讯是利用ADS协议实现的,ADS协议是倍福公司基于TCP/IP协议开发的用于各个模块之间的信息交换和PLC与外部环境通讯的接口协议,利用协议提供的端口和库函数可方便获取下位机PLC信息并传至上位机。
附图说明
下面结合附图中的实施例对本发明作进一步的详细说明,但并不构成对本发明的任何限制。
图1为本发明系统的总体架构图。
图2为系统人机界面。
图3为NC代码主要指令图。
图4为“S”型加减速曲线示意图。
图5为全局速度规划流程图。
图6为曲线分段点示意图。
图7为机床实物图。
附图标记说明如下:
上位机1,下位机2,人机界面模块3,译码器4,全局速度规划模块5,激光功率优化模块6,插补模块7,坐标变换模块8,位置补偿模块9,逻辑处理模块10,伺服和激光控制模块11。
具体实施方式
实施例1:以用于航空发动机机匣表面刻型的六轴五联动激光刻型机为例对本发明进行具体说明。
实施例1的六轴五联动激光加工开放式数控系统的人机界面是在Qt集成开发环境上利用QML语言设计的,它包含八个子窗口:主加工,参数设置,系统设置,伺服设置,文档管理,三维测量,螺距补偿,状态监控,如图2所示。
译码器的功能是将NC代码从文本转换为系统可以处理的命令队列。NC代码的主要命令如图3所示。特别的,激光的打开和关闭分别由M03和M05控制。D、L和S分别表示激光器的占空比,功率和频率。为了方便地匹配命令文本,译码器利用了正则表达式,例如:
(\\s*(([DGMNPS]\\d+)|([ABCFIJKXYZ][+]?\\d+\\.?\\d*)|([ABCFIJKXYZ][+-]?\\d*\\.?\\d+)|([AZ]{2,})))*(\\s*[;%].*)?
([A-Z])([+-]?\\d*\\.?\\d*)
第一个表达式可以用作语法检查语句,无法与其匹配的任何单行文本都被视为无效命令。第二个表达式可以用作单段命令检查语句,以提取此命令的类型和数据值。
全局速度规划模块是整个系统的重要组成部分。合理的速度曲线既可以提高精度,又可以节省加工时间。速度规划首先计算前后曲线段衔接点的过渡速度(即第i段曲线的末速度或者第i+1段曲线的起速度)vi
Figure BDA0002540236790000071
其中Amax为机床最大加速度,Tc为插补周期,F为指令速度,θ=arccos(τ12),τ1和τ2分别是第i曲线段的末端单位切矢和第i+1曲线段的起始单位切矢;
在获得每段曲线的末速度之后,判断每段曲线速度是否可以达到。如果使用如图4所示的“S”型速度曲线,假设vi-1和vi是第i段的起末速度,如果加(减)速过程中有匀加(减)速阶段,则从vi-1加速或减速到vi的最短距离Smin为:
Figure BDA0002540236790000072
如果加(减)速过程中没有匀加(减)速阶段,则最短距离Smin为:
Figure BDA0002540236790000081
其中vmax=max(vi-1,vi),vmin=min(vi-1,vi),vi-1和vi分别为第i段曲线的起末速度,Am和Jm分别为机床所能达到的最大加速度和最大捷度;
假设第i段曲线段的长度为L。如果L>Smin则vi是可达的。否则,需要调整vi的值。调整步骤如下:1.如果vi-1≤vi,此时是加速过程,则直接根据L减小vi。2.如果vi-1>vi,此时是减速过程,我们可以根据L减小vi-1。但是由于前一段的末速度已更改,则,需要回溯至前一段再次以上述方式判断是否可达。
在确定各段衔接速度都可达后,需要继续判断指令速度F是否可达。首先,将vi-1和F带入公式(2)或(3),计算最小距离Smin1。同样,将vi和F带入(2)或(3),计算最小距离Smin2。如果L>Smin1+Smin2,则指令速度可达,否则根据公式(4)修改指令速度。
Figure BDA0002540236790000082
全局速度规划的算法流程图如图5所示。
在全局速度规划后进行激光功率优化,使其与进给速度匹配,避免出现拐角烧蚀现象。调整激光功率应首先确定段间分割点。判断条件有两个:最小转角限制和最大弓高误差限制。如图6所示为曲线分段点示意图,图中实线为指令轨迹,虚线为实际轨迹,则最小转角限制为:
θcr=2×arccos(Kpε/F/(1-Kf)) (5)
其中F为指令速度,ε为预设最大转角误差,Kf和Kp分别为伺服前馈增益和位置增益。
然后计算两个最大弓高误差限制为:
δ1=R(1-cosα1) (6)
δ2=R(1-cos(π-θ-α1)) (7)
其中R=2×l1/sinα1,α1=arctan(l1sin(π-θ)/(l2+l1cos(π-θ))),l1和l2分别为前后曲线段的长度,θ为前后两曲线段的切矢夹角。
确认分割点后,可通过公式(1)计算分割点的速度(vi)。假设指令速度和指令激光功率分别为vcmd和Pcmd,则更改后的激光功率为:
Pi=viPcmd/vcmd (8)
插补器计算的任务是根据所规划的速度曲线计算位置点,驱动各轴运动。图3所示的速度曲线的表达式为:
Figure BDA0002540236790000091
其中:A和J分别为机床最大加速度和最大捷度,
Figure BDA0002540236790000092
Figure BDA0002540236790000093
v3=F-J(t5-t4)2/2,v4=F-J(t5-t4)2/2-A(t6-t5)。
根据分段点前后速度相等,可计算分段时间的表达式为:
t1=t3-t2=t5-t4=t7-t6=A/J
t2-t1=(F-A2/J-vs)/A
t6-t5=(F-A2/J-ve)/A
Figure BDA0002540236790000094
根据公式(9)和公式(10),计算插补步长为:
ΔL=vtTc (11)
其中Tc为插补周期。
根据步长计算下一个周期节点矢量ui+1,其中直线、圆弧、NURBS曲线插补节点矢量计算分别如下:
直线的插补节点矢量递推计算:ui+1=ui+Tc
圆弧的插补节点矢量递推计算:ui+1=ui+ΔL/R
NURBS曲线的插补节点矢量递推计算:
Figure BDA0002540236790000095
其中,R为圆的半径值,C′(ui)和C″(ui)分别为NURBS曲线的一阶和二阶导数,||·||为向量的模,<·>为向量内积;ui是当前周期曲线在参数域上的的节点矢量,ui+1是下一周期曲线在参数域上的的节点矢量,节点矢量是参数曲线的重要参数,一个参数域的节点矢量对应曲线上一个点;
将下一周期的节点矢量代入曲线表达式计算插补点坐标,其中直线、圆弧、NURBS曲线在ui+1处的坐标Pi+1计算分别如下:
直线的坐标计算:
Figure BDA0002540236790000101
其中,Pi+1是节点参数为ui+1时曲线上的点;Ps为起点坐标,Pe为终点坐标,us为起点的节点矢量,ue为终点的节点矢量(给定插补曲线后us、ue为已知量,默认:us=0,ue=1);
圆弧的坐标计算:Pi+1=Pc+R(ui+1-us,N)(Ps-Pc)
其中,R(ui+1-us,N)是表示向量绕N旋转ui+1-us角度,N为空间圆的法向量,Pc为圆心坐标,Ps起点坐标,us起点的节点矢量;
NURBS曲线的基函数递推计算如下:
Figure BDA0002540236790000102
其中i为节点序号,k为曲线次数,Ni,k(·)为基函数;u表示节点矢量,是已知量,给定一条NURBS加工轨迹就会提供相应的节点矢量;
根据NURBS曲线的基函数计算在参数ui+1处的插补坐标Pi+1(ui+1)如下:
Figure BDA0002540236790000103
其中di为控制顶点,wi为权因子,
Figure BDA0002540236790000104
为基函数的累加。
得到位置坐标后,需要将坐标变换后的位置发送到伺服轴。六轴五联动激光刻型机从激光焦点到工件的逆坐标变换为:
px=X+x2+x1 cos B-y1 sin B-ty cos A sin B+tz sin A sin B
py=Y+y2+x1 sin B+y1 cos B+ty cos A cos B-tz sin A cos B
pz=Z+z2+tzcos A+ty sin A
i=sin A sin B
j=-sin A cos B
k=cos A (16)
反解得到正坐标变换为:
A=-arccos(k)
B=arctan(-i/j)
X=px-x2-x1 cos B+y1 sin B+ty cos A sin B-tz sin A sin B
Y=py-y2-x1 sin B-y1 cos B-ty cos A cos B+tz sin A cos B
Z=pz-z2-tzcos A-tysin A (17)
其中ty,tz为激光焦点与A轴坐标系之间的距离,x1,y1为A轴坐标系到B轴坐标系的距离,A、B为两个旋转轴的运动角度,X,Y,Z为三个线性轴的运动值,x3,y3和z3为测得的校准常数。
坐标转换后,每个轴的位置值加上补偿值可以驱动轴移动。补偿值由激光干涉仪测量,手动输入到HMI。运动控制是使用TwinCAT PTP提供的Tc2-MC2库实现的。
如图7所示为机床实物图,配有倍福C6640-0040工控机,采用CoreTM i7 2.3GHzCPU,8G运行内存,系统采用win7专业版;机床旋转轴采用Kollmorgen AKD驱动器和支持EtherCAT协议的DDR D06X系列电机,线性轴采用SIEMENS S120驱动器和1FK7系列电机,支持PROFIdrive协议。
以上所举实施例为本发明的较佳实施方式,仅用来方便说明本发明,并非对本发明作任何形式上的限制,任何所属技术领域中具有通常知识者,若在不脱离本发明所提技术特征的范围内,利用本发明所揭示技术内容所作出局部更动或修饰的等效实施例,并且未脱离本发明的技术特征内容,均仍属于本发明技术特征的范围内。

Claims (10)

1.一种六轴五联动激光加工开放式数控系统,包括上位机(1)和下位机(2);其特征在于:
上位机包括四个模块:人机界面模块(3)、译码器(4)、全局速度规划模块(5)、激光功率优化模块(6);
其中,人机界面模块(3)用于输入加工程序和显示各类信息;
其中,译码器(4)用于将输入的NC代码解释成指令队列;
其中,全局速度规划模块(5)根据机床动力学限制修改指令队列中的超约束值;
其中,激光功率优化模块(6)用于调整激光功率与进给速度相匹配;
下位机(2)包括五个模块:插补计算模块(7)、坐标变换模块(8)、位置补偿模块(9)、逻辑处理模块(10)、伺服和激光控制模块(11);
其中,插补计算模块(7)用来将轨迹曲线离散成位置坐标;
其中,坐标变换模块(8)用来进行六轴五联动正逆坐标变换;
其中,补偿模块(9)用来对位置误差进行补偿以提高定位精度;
其中,逻辑处理模块(10)用来处理外部输入IO信号;
其中,伺服和激光控制模块(11)用来驱动轴和激光器动作;
上位机各组成模块信息流程如下:
1)人机界面接收加工代码传至译码器;
2)译码器将加工代码转换成指令队列方便机器内部传送和处理,并将指令队列传至全局速度规划模块;
3)全局速度规划模块在接收指令队列后根据机床的动力学限制修改队列中的超约束值,并将修改后的指令队列传至激光功率优化模块;
4)激光功率优化模块接收指令队列后修改激光的功率使之与进给速度相匹配;至此上位机非实时任务完成,得到了适合加工的指令队列;
上位机处理好的指令队列通过ADS标准通信协议传至下位机;
下位机各组成模块信息流程如下:
1)上位机处理好的指令队列首先进入插补计算模块将指令队列中的曲线轨迹离散成位置坐标,并传至坐标变换模块;
2)坐标变换模块将位置坐标分解成各轴运动量传给补偿模块;
3)补偿模块根据激光干涉仪检测的机床补偿值将调整各轴运动量的值,以提高定位精度;
4)伺服和激光控制模块接收到各轴运动量和激光指令功率后利用TwinCAT PTP功能驱动伺服轴和激光器动作;
5)逻辑处理模块的任务是处理外部输入的IO信号,实时在线工作。
2.根据权利要求1所述的六轴五联动激光加工开放式数控系统,其特征在于,人机界面模块(3)基于Qt平台开发,利用QML和C++语言开发系统前端和后端。
3.根据权利要求1所述的六轴五联动激光加工开放式数控系统,其特征在于,译码器(4)利用正则表达式进行文本匹配,提取指令信息,将指令信息存入指令队列供速度规划模块(5)使用。
4.根据权利要求1所述的六轴五联动激光加工开放式数控系统,其特征在于,全局速度规划模块(5),包括以下步骤:
步骤一:接收到来自译码器(4)的指令队列后首先根据机床联动动力学约束计算前后曲线段衔接点的过渡速度(即第i段曲线的末速度或者第i+1段曲线的起速度)vi为:
Figure FDA0002540236780000021
其中Amax为机床最大加速度,Tc为插补周期,F为指令速度,θ=arccos(τ12),τ1和τ2分别是第i曲线段的末端单位切矢和第i+1曲线段的起始单位切矢;
步骤二:计算第i段曲线段加(减)速所需最小距离Smin
其中当加(减)速过程中有匀加(减)速阶段的计算方式为:
Figure FDA0002540236780000022
反之,当加(减)速过程中没有匀加(减)速阶段的计算方式:
Figure FDA0002540236780000023
其中vmax=max(vi-1,vi),vmin=min(vi-1,vi),vi-1和vi分别为第i段曲线的起末速度,Am和Jm分别为机床所能达到的最大加速度和最大捷度;
步骤三:如果曲线段长度L>Smin,则起末速度可达,反之,则需根据L的大小调整起末速度;
步骤四:判断指令速度是否可达,如果不可达根据以下公式降低指令速度:
Figure FDA0002540236780000024
其中F为指令速度,Jm为机床所能达到的最大捷度;
完成上述步骤后,指令队列中的速度超约束值已经修改完毕,接下来修改激光功率使之与进给速度相匹配。
5.根据权利要求1所述的六轴五联动激光加工开放式数控系统,其特征在于,激光功率优化模块(6)具有以下流程:
步骤一:计算两相邻曲线段最大转角限制θcr
Figure FDA0002540236780000031
其中F为指令速度,ε为预设最大转角误差,Kf和Kp分别为伺服前馈增益和位置增益;
步骤二:计算前后两段曲线的弓高误差δ1、δ2
δ1=R(1-cosα1)
δ2=R(1-cos(π-θ-α1))
其中R为两曲线段外接圆半径,α1=arctan(l1sin(π-θ)/(l2+l1cos(π-θ))),l1和l2分别为前后两段曲线段的长度,θ为前后两曲线段的切矢夹角;
步骤三:如果前后两曲线段的切矢夹角θ≤θcr或者预设最大弓高误差δ<min(δ12),则根据衔接点速度修改拐点处的激光功率Pi为:
Pi=viPcmd/vcmd
其中vi为衔接点速度,vcmd和Pcmd分别为指令速度和指令激光功率;
完成上述步骤后,激光功率已经和进给速度相匹配,接下来通过ADS协议将指令队列传至下位机。
6.根据权利要求1所述的六轴五联动激光加工开放式数控系统,其特征在于,插补计算模块(7)具有以下流程:
步骤一:根据“S”型速度曲线计算插补步长ΔL为:
ΔL=vtTc
其中vt为速度曲线上t时刻的指令速度,Tc为插补周期;
步骤二:根据步长计算下一个周期节点矢量ui+1,其中直线、圆弧、NURBS曲线插补节点矢量计算分别如下:
直线的插补节点矢量递推计算:ui+1=ui+Tc
圆弧的插补节点矢量递推计算:ui+1=ui+ΔL/R
NURBS曲线的插补节点矢量递推计算:
Figure FDA0002540236780000032
其中,R为圆的半径值,C′(ui)和C″(ui)分别为NURBS曲线的一阶和二阶导数,||·||为向量的模,<·>为向量内积;ui是当前周期曲线在参数域上的的节点矢量,ui+1是下一周期曲线在参数域上的的节点矢量(节点矢量是参数曲线的重要参数,一个参数域的节点矢量对应曲线上一个点);
步骤三:将下一周期的节点矢量代入曲线表达式计算插补点坐标,其中直线、圆弧、NURBS曲线在ui+1处的坐标Pi+1计算分别如下:
直线的坐标计算:
Figure FDA0002540236780000041
其中,Pi+1是节点参数为ui+1时曲线上的点,Ps为起点坐标,Pe为终点坐标,us为起点的节点矢量,ue为终点的节点矢量(给定插补曲线后us、ue为已知量,默认:us=0,ue=1);
圆弧的坐标计算:Pi+1=Pc+R(ui+1-us,N)(Ps-Pc)
其中,R(ui+1-us,N)是表示向量绕N旋转ui+1-us角度,N为空间圆的法向量,Pc为圆心坐标,Ps起点坐标,us起点的节点矢量;
NURBS曲线的基函数递推计算如下:
Figure FDA0002540236780000042
其中i为节点序号,k为曲线次数,Ni,k(·)为基函数,u表示节点矢量(已知量,给定一条NURBS加工轨迹就会提供相应的节点矢量);
根据NURBS曲线的基函数计算在参数ui+1处的插补坐标Pi+1(ui+1)如下:
Figure FDA0002540236780000043
其中di为控制顶点,wi为权因子,
Figure FDA0002540236780000044
为基函数的累加;
完成上述步骤后,指令曲线已经离散成位置坐标值,接下来传至坐标变换模块分解成各轴运动量。
7.根据权利要求1所述的六轴五联动激光加工开放式数控系统,其特征在于,在测量出七个机床固有参数后,坐标变换模块(9)能够适应AB轴结构形式的正逆坐标变换;完成坐标变换后,可得到各轴的运动量,接下来传至补偿模块进行调整。
8.根据权利要求1所述的六轴五联动激光加工开放式数控系统,其特征在于,经过坐标变换后的各轴位置在补偿模块(10)中进行双向补偿,有效提高机床定位精度。
9.根据权利要求1所述的六轴五联动激光加工开放式数控系统,其特征在于,伺服和激光控制模块(12)利用TwinCAT提供的Tc2-MC2库,利用外部使能功能实现点对点运动控制。
10.一种六轴五联动激光加工开放式数控系统的工作方法,其特征在于,采用如权利要求1至9任意一项所述的一种六轴五联动激光加工开放式数控系统,包括以下步骤:
系统部署到机床的工控机内,系统中各轴变量、激光变量以及IO等变量分别和机床的伺服驱动器、激光器以及IO模块等硬件进行扫描连接;
机床开机自动运行系统程序,人机交互模块接收加工程序和显示系统状态信息,接收的加工程序经过系统各个模块的处理优化后,驱动硬件完成加工。
CN202010544778.4A 2020-06-15 2020-06-15 一种六轴五联动激光加工开放式数控系统及其工作方法 Active CN111665785B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010544778.4A CN111665785B (zh) 2020-06-15 2020-06-15 一种六轴五联动激光加工开放式数控系统及其工作方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010544778.4A CN111665785B (zh) 2020-06-15 2020-06-15 一种六轴五联动激光加工开放式数控系统及其工作方法

Publications (2)

Publication Number Publication Date
CN111665785A true CN111665785A (zh) 2020-09-15
CN111665785B CN111665785B (zh) 2021-06-15

Family

ID=72387832

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010544778.4A Active CN111665785B (zh) 2020-06-15 2020-06-15 一种六轴五联动激光加工开放式数控系统及其工作方法

Country Status (1)

Country Link
CN (1) CN111665785B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112051803A (zh) * 2020-09-17 2020-12-08 华南理工大学广州学院 一种基于空间平面法向量的小线段光顺方法
CN112405107A (zh) * 2020-11-26 2021-02-26 北京航空航天大学 一种数控机床通信异构进给系统及其调试方法
CN113848817A (zh) * 2021-10-11 2021-12-28 东北大学 一种智能五轴控制系统及方法
CN115616982A (zh) * 2021-07-12 2023-01-17 大族激光科技产业集团股份有限公司 一种拐弯轨迹的过度方法及激光装置和计算机设备
CN117226302A (zh) * 2023-11-10 2023-12-15 迈为技术(珠海)有限公司 激光切割控制方法、装置、设备以及存储介质
CN115616982B (zh) * 2021-07-12 2024-06-04 大族激光科技产业集团股份有限公司 一种拐弯轨迹的过度方法及激光装置和计算机设备

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104076742A (zh) * 2014-06-20 2014-10-01 北京航空航天大学 一种连续多类型曲线段混合插补三次多项式进给速度前瞻规划方法
CN104714483A (zh) * 2013-12-11 2015-06-17 中国科学院宁波材料技术与工程研究所 激光加工的控制方法及系统
CN107491042A (zh) * 2017-10-25 2017-12-19 卢俊 一种6轴嵌入式数控系统
CN109416529A (zh) * 2017-12-19 2019-03-01 深圳配天智能技术研究院有限公司 数控系统、具有存储功能的装置及nurbs曲线插补方法
CN208673132U (zh) * 2018-06-01 2019-03-29 武汉富乐瑞激光工程有限公司 一种辊类表面激光毛化数控系统
US10466676B1 (en) * 2019-04-08 2019-11-05 SmileDirectClub LLC Systems and methods for shearing adjustment of a dental aligner cutting tool
CN110919212A (zh) * 2019-12-12 2020-03-27 佛山市宏石激光技术有限公司 一种激光切割机功率控制方法及系统

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104714483A (zh) * 2013-12-11 2015-06-17 中国科学院宁波材料技术与工程研究所 激光加工的控制方法及系统
CN104076742A (zh) * 2014-06-20 2014-10-01 北京航空航天大学 一种连续多类型曲线段混合插补三次多项式进给速度前瞻规划方法
CN107491042A (zh) * 2017-10-25 2017-12-19 卢俊 一种6轴嵌入式数控系统
CN109416529A (zh) * 2017-12-19 2019-03-01 深圳配天智能技术研究院有限公司 数控系统、具有存储功能的装置及nurbs曲线插补方法
CN208673132U (zh) * 2018-06-01 2019-03-29 武汉富乐瑞激光工程有限公司 一种辊类表面激光毛化数控系统
US10466676B1 (en) * 2019-04-08 2019-11-05 SmileDirectClub LLC Systems and methods for shearing adjustment of a dental aligner cutting tool
CN110919212A (zh) * 2019-12-12 2020-03-27 佛山市宏石激光技术有限公司 一种激光切割机功率控制方法及系统

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112051803A (zh) * 2020-09-17 2020-12-08 华南理工大学广州学院 一种基于空间平面法向量的小线段光顺方法
CN112051803B (zh) * 2020-09-17 2022-02-08 华南理工大学广州学院 一种基于空间平面法向量的小线段光顺方法
CN112405107A (zh) * 2020-11-26 2021-02-26 北京航空航天大学 一种数控机床通信异构进给系统及其调试方法
CN112405107B (zh) * 2020-11-26 2021-09-28 北京航空航天大学 一种数控机床通信异构进给系统及其调试方法
CN115616982A (zh) * 2021-07-12 2023-01-17 大族激光科技产业集团股份有限公司 一种拐弯轨迹的过度方法及激光装置和计算机设备
CN115616982B (zh) * 2021-07-12 2024-06-04 大族激光科技产业集团股份有限公司 一种拐弯轨迹的过度方法及激光装置和计算机设备
CN113848817A (zh) * 2021-10-11 2021-12-28 东北大学 一种智能五轴控制系统及方法
CN117226302A (zh) * 2023-11-10 2023-12-15 迈为技术(珠海)有限公司 激光切割控制方法、装置、设备以及存储介质

Also Published As

Publication number Publication date
CN111665785B (zh) 2021-06-15

Similar Documents

Publication Publication Date Title
CN111665785B (zh) 一种六轴五联动激光加工开放式数控系统及其工作方法
US6823234B2 (en) Curve interpolation method
CN105500354A (zh) 一种工业机器人应用的过渡轨迹规划方法
CN101825884B (zh) 用五轴四联动的数控系统实现五轴五联动轨迹控制的方法
Zhang et al. High speed interpolation for micro-line trajectory and adaptive real-time look-ahead scheme in CNC machining
CN104678894A (zh) 数控加工路径的规划方法、数控加工系统和方法
CN103294007B (zh) 一种高速高精度柔性电子齿轮箱的控制方法
CN103699056A (zh) 高速高精度数控加工的小线段实时平滑过渡插补方法
Chen et al. Contour error–bounded parametric interpolator with minimum feedrate fluctuation for five-axis CNC machine tools
Huang et al. A newly developed corner smoothing methodology based on clothoid splines for high speed machine tools
US6528963B1 (en) Robot and method for controlling motor speed of the robot
CN1102890A (zh) 可对曲面轮廓直接插补的cnc系统
CN108062071A (zh) 参数曲线轨迹伺服轮廓误差的实时测定方法
Zhang et al. A local and analytical curvature-smooth method with jerk-continuous feedrate scheduling along linear toolpath
Tajima et al. Kinematic tool-path smoothing for 6-axis industrial machining robots
CN111590570B (zh) 一种同步交叉耦合机器人轮廓控制方法
Calabrese et al. Design and realization of a STEP-NC compliant CNC embedded controller
CN112379638A (zh) 一种基于有限脉冲响应滤波的样条曲线插补方法
CN114739290B (zh) 线激光扫描化铣胶刻线的路径规划方法及系统
Yaman et al. A novel command generation paradigm for production machine systems
JPH07152417A (ja) 数値制御装置の工具径路および送り速度制御方式
JPH01177617A (ja) インボリュート補間方式
Yi et al. Research on numerical control of the mobile robotic machine tool
Yin et al. Open CNC System Design for Multi-axis Laser Manufacturing Machine Tools
US20240103481A1 (en) Numerical controller

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