CN116339238B - A beam motion control method for five-axis laser processing equipment with independent beam scanning - Google Patents
A beam motion control method for five-axis laser processing equipment with independent beam scanning Download PDFInfo
- Publication number
- CN116339238B CN116339238B CN202310336954.9A CN202310336954A CN116339238B CN 116339238 B CN116339238 B CN 116339238B CN 202310336954 A CN202310336954 A CN 202310336954A CN 116339238 B CN116339238 B CN 116339238B
- Authority
- CN
- China
- Prior art keywords
- axis
- machine tool
- motion
- coordinate system
- plane
- 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
- 238000012545 processing Methods 0.000 title claims abstract description 87
- 238000000034 method Methods 0.000 title claims abstract description 43
- 239000013598 vector Substances 0.000 claims description 43
- 230000009466 transformation Effects 0.000 claims description 29
- 239000011159 matrix material Substances 0.000 claims description 24
- 238000003754 machining Methods 0.000 claims description 23
- 238000013519 translation Methods 0.000 claims description 11
- 230000001131 transforming effect Effects 0.000 claims 2
- 239000002131 composite material Substances 0.000 abstract description 3
- 230000009471 action Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 9
- 230000036544 posture Effects 0.000 description 7
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 101100012902 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) FIG2 gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
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/408—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 data handling or data format, e.g. reading, buffering or conversion of data
- G05B19/4086—Coordinate conversions; Other special calculations
-
- 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/33—Director till display
- G05B2219/33258—Common coordinate conversion for multiple heads, spindles
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Engineering & Computer Science (AREA)
- Human Computer Interaction (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Numerical Control (AREA)
Abstract
本发明属于一种光束运动控制方法,为解决多轴激光加工装备加工时,若激光光束的运动图形不是圆形,会导致错位加工,造成零件报废的技术问题,提供一种独立光束扫描的五轴激光加工设备光束运动控制方法,通过分析五轴机床运动结构及连接关系,将光束扫描装置的运动轨迹分解到各运动轴上,获取各轴复合作用下的光束扫描装置相对工件的位恣,继而建立当前位恣下的局部扫描坐标系,实现带独立光束扫描装置的五轴激光加工设备运动控制。
The present invention relates to a method for controlling the motion of a light beam. In order to solve the technical problem that, when a multi-axis laser processing equipment is processed, if the motion pattern of the laser beam is not circular, misaligned processing will occur, resulting in the scrapping of parts, a method for controlling the motion of a light beam of a five-axis laser processing equipment with independent light beam scanning is provided. By analyzing the motion structure and connection relationship of a five-axis machine tool, the motion trajectory of a light beam scanning device is decomposed into each motion axis, the position of the light beam scanning device relative to a workpiece under the composite action of each axis is obtained, and then a local scanning coordinate system under the current position is established to realize the motion control of a five-axis laser processing equipment with an independent light beam scanning device.
Description
技术领域Technical Field
本发明属于一种光束运动控制方法,具体涉及一种独立光束扫描的五轴激光加工设备光束运动控制方法。The invention relates to a light beam motion control method, and in particular to a light beam motion control method for five-axis laser processing equipment with independent light beam scanning.
背景技术Background technique
多轴机床装备由于具有更加灵活的加工位置控制方法,被广泛应用于航空航天等高科技领域,用以实现叶片、机匣、肋板等复杂零件的加工。为了提高多轴机床的适用性,降低机床控制代码的编写难度,一般采用带刀尖跟随(RTCP)功能的数控代码编写方式,即工艺人员只需要关注于如何控制刀具与工件的相对运动,而不需要关注机床各个轴的控制方法。Multi-axis machine tools are widely used in high-tech fields such as aerospace due to their more flexible processing position control methods, to achieve the processing of complex parts such as blades, casings, ribs, etc. In order to improve the applicability of multi-axis machine tools and reduce the difficulty of writing machine tool control codes, the CNC code writing method with tool tip following (RTCP) function is generally adopted, that is, the process personnel only need to focus on how to control the relative movement of the tool and the workpiece, and do not need to pay attention to the control method of each axis of the machine tool.
随着激光加工技术的发展,激光加工领域也出现了大量的多轴激光加工装备。这些装备通常包含了多个机械轴和多个光束轴。为了减轻代码的编写难度,在激光加工装备中也采用了RTCP功能。但是,该功能来自于机械加工装备,只适用于机械轴的控制代码,无法控制光束轴,使得加工过程中机械轴坐标系与光束轴坐标系之间始终存在一定的角度偏差。一旦激光光束的运动图形不是圆形,则待加工图形的角度很可能发生偏转,导致错位加工,造成零件报废。With the development of laser processing technology, a large number of multi-axis laser processing equipment have also appeared in the field of laser processing. These equipment usually contain multiple mechanical axes and multiple beam axes. In order to reduce the difficulty of code writing, the RTCP function is also used in laser processing equipment. However, this function comes from mechanical processing equipment and is only applicable to the control code of the mechanical axis. It cannot control the beam axis, so that there is always a certain angle deviation between the mechanical axis coordinate system and the beam axis coordinate system during the processing. Once the motion pattern of the laser beam is not circular, the angle of the pattern to be processed is likely to deflect, resulting in misaligned processing and scrapped parts.
发明内容Summary of the invention
本发明为解决多轴激光加工装备加工时,若激光光束的运动图形不是圆形,会导致错位加工,造成零件报废的技术问题,提供一种独立光束扫描的五轴激光加工设备光束运动控制方法。The present invention provides a beam motion control method for five-axis laser processing equipment with independent beam scanning to solve the technical problem that if the motion pattern of the laser beam is not circular during multi-axis laser processing, it will cause misaligned processing and scrapped parts.
为达到上述目的,本发明采用以下技术方案予以实现:In order to achieve the above object, the present invention adopts the following technical solutions:
一种独立光束扫描的五轴激光加工设备光束运动控制方法,其特殊之处在于,包括以下步骤:A beam motion control method for a five-axis laser processing device with independent beam scanning is characterized in that it includes the following steps:
S1,使机床A轴正向旋转,获取机床C轴运动分量;所述机床为五轴机床,五轴分别为X轴、Y轴、Z轴、A轴和C轴;S1, making the A axis of the machine tool rotate in the positive direction to obtain the motion component of the C axis of the machine tool; the machine tool is a five-axis machine tool, and the five axes are respectively X-axis, Y-axis, Z-axis, A-axis and C-axis;
S2,根据振镜加工平面法向TM到加工坐标系的-YMOMZM平面的旋转变换矩阵R-C1,对振镜加工平面法向TM进行变换,得到变换后的矢量TC,M1;其中,加工坐标系为OM-XMYMZM;S2, transform the galvanometer processing plane normal TM according to the rotation transformation matrix R -C1 from the galvanometer processing plane normal TM to the -Y M O M Z M plane of the processing coordinate system, and obtain the transformed vector TC , M1 ; wherein the processing coordinate system is O M -X M Y M Z M ;
S3,结合变换后的矢量TC,M1,获取机床A轴运动分量;S3, combining the transformed vector TC , M1 , to obtain the motion component of the machine tool A axis;
S4,结合机床C轴运动分量和机床A轴运动分量,计算机床在其X轴、Y轴、Z轴的平移分量;S4, combining the C-axis motion component of the machine tool and the A-axis motion component of the machine tool, calculates the translation component of the machine tool on its X-axis, Y-axis, and Z-axis;
S5,判断机床A轴正向旋转下,机床X轴、Y轴、Z轴、A轴是否超过机床行程限制,若其中任一轴超过机床行程限制,则执行步骤S6,否则,执行步骤S8;S5, determining whether the X-axis, Y-axis, Z-axis, and A-axis of the machine tool exceed the machine tool travel limit when the machine tool A axis rotates in the forward direction. If any of the axes exceeds the machine tool travel limit, execute step S6; otherwise, execute step S8;
S6,使机床A轴向负向旋转,获取机床C轴运动分量;根据振镜加工平面法向TM到加工坐标系的+YMOMZM平面的旋转变换矩阵R-C2,对振镜加工平面法向TM进行变换,得到变换后的矢量TC,M2;结合变换后的矢量TC,M2,获取机床A轴运动分量;结合机床C轴运动分量和机床A轴运动分量,计算机床在其X轴、Y轴、Z轴、A轴的平移分量;S6, rotate the A axis of the machine tool in the negative direction to obtain the motion component of the C axis of the machine tool; transform the normal direction TM of the galvanometer processing plane to the + YMOMZM plane of the processing coordinate system according to the rotation transformation matrix R -C2 , to obtain the transformed vector TC ,M2 ; obtain the motion component of the A axis of the machine tool by combining the transformed vector TC ,M2 ; calculate the translation components of the machine tool on its X axis, Y axis, Z axis and A axis by combining the motion component of the C axis of the machine tool and the motion component of the A axis of the machine tool;
S7,判断机床A轴正向旋转下,机床X轴、Y轴、Z轴、A轴是否超过机床行程限制,若其中任一轴超过机床行程限制,更换零件装夹位置或更换机床,否则,执行步骤S8;S7, determining whether the X-axis, Y-axis, Z-axis, and A-axis of the machine tool exceed the machine tool travel limit when the machine tool A axis rotates forward. If any of the axes exceeds the machine tool travel limit, change the part clamping position or replace the machine tool. Otherwise, execute step S8;
S8,以机床A轴旋转,并按照在其X轴、Y轴、Z轴、A轴的平移分量运动后的位姿为参考,计算XG轴、YG轴在加工坐标系下的轴向;其中,OG-XGYG为2D振镜的扫描平面坐标系;S8, with the A axis of the machine tool rotating, and taking the position after the translational components of the X axis, Y axis, Z axis, and A axis move as a reference, calculate the axial directions of the XG axis and the YG axis in the machining coordinate system; wherein OG - XGYG is the scanning plane coordinate system of the 2D galvanometer;
S9,根据XG轴、YG轴在加工坐标系下的轴向,建立振镜扫描坐标系,并在该坐标系中生成振镜扫描轨迹,结合机床各轴运动分量,得到光束运动的完整参数。S9, according to the axial directions of the XG axis and the YG axis in the machining coordinate system, a galvanometer scanning coordinate system is established, and a galvanometer scanning trajectory is generated in the coordinate system, and the complete parameters of the beam motion are obtained by combining the motion components of each axis of the machine tool.
进一步地,步骤S1具体为:Furthermore, step S1 is specifically as follows:
S1.1,将振镜加工平面法向TM绕加工坐标系的ZM轴旋转到加工坐标系的-YMOMZM平面上;S1.1, rotate the normal T M of the galvanometer processing plane around the Z M axis of the processing coordinate system to the -Y M O M Z M plane of the processing coordinate system;
S1.2,通过下式,得到机床A轴向正向旋转时机床C轴运动分量CS1:S1.2, the motion component C S1 of the C axis of the machine tool when the A axis of the machine tool rotates in the positive direction is obtained by the following formula:
当振镜加工平面法向TM在YM轴上的分量jM=0时,机床C轴运动分量CS1为:When the component j M of the normal direction TM of the galvanometer processing plane on the Y M axis = 0, the motion component CS1 of the C axis of the machine tool is:
其中,iM为振镜加工平面法向TM在XM轴上的分量;Wherein, i M is the component of the normal T M of the galvanometer processing plane on the X M axis;
当振镜加工平面法向TM在YM轴上的分量jM≠0时,机床C轴运动分量CS1为:When the component j M of the normal direction TM of the galvanometer processing plane on the Y M axis ≠ 0, the motion component CS1 of the C axis of the machine tool is:
其中, in,
进一步地,步骤S2中,所述旋转变换矩阵R-C1通过下式得到:Further, in step S2, the rotation transformation matrix R -C1 is obtained by the following formula:
变换后的矢量TC,M1通过下式得到:The transformed vector T C,M1 is obtained by the following formula:
TC,M1=TM·R-C1。 TC,M1 = TM ·R -C1 .
进一步地,步骤S3具体为,通过下式获取机床A轴运动分量:Furthermore, step S3 is specifically to obtain the motion component of the machine tool A axis by the following formula:
当kC,M1=0时,机床A轴运动分量AS1=π/2;其中,kC,M1为变换后的矢量TC,M1在ZM轴上的分量;When k C, M1 = 0, the motion component of the machine tool A axis A S1 = π/2; where k C, M1 is the component of the transformed vector T C, M1 on the Z M axis;
当kC,M1≠0时,机床A轴运动分量AS1为:When k C,M1 ≠ 0, the motion component AS1 of the machine tool A axis is:
其中,jC,M1为变换后的矢量TC,M1在YM轴上的分量。in, j C,M1 is the component of the transformed vector T C,M1 on the Y M axis.
进一步地,步骤S4具体为,通过下式计算机床在其X轴、Y轴、Z轴的平移分量XS1、YS1和ZS1:Furthermore, step S4 is specifically to calculate the translation components X S1 , Y S1 and Z S1 of the tool machine in the X-axis, Y-axis and Z-axis by the following formula:
其中,PM为振镜加工平面原点,(Δx1,Δy1,Δz1)为机床A轴回转中心到C轴回转中心的偏移量,(Δx2,Δy2,Δz2)为机床原点到A轴回转中心的偏移量,R-A1为变换后的矢量TC,M1到+ZM轴的旋转变换矩阵:Among them, PM is the origin of the galvanometer processing plane, ( Δx1 , Δy1 , Δz1 ) is the offset from the A-axis rotation center of the machine tool to the C-axis rotation center, ( Δx2 , Δy2 , Δz2 ) is the offset from the machine tool origin to the A-axis rotation center, R -A1 is the transformed vector TC, and the rotation transformation matrix from M1 to the + ZM axis is:
进一步地,步骤S5中,所述机床行程限制通过下式确定:Further, in step S5, the machine tool travel limit is determined by the following formula:
其中,[x-,x+]为机床的X轴运动范围,[y-,y+]为机床的Y轴运动范围,[z-,z+]为机床的Z轴运动范围,[α-,α+]为机床的A轴运动范围,(x0,y0,z0,α0,γ0)为机床对刀后的原点位置。Among them, [x - , x + ] is the X-axis motion range of the machine tool, [y - , y + ] is the Y-axis motion range of the machine tool, [z - , z + ] is the Z-axis motion range of the machine tool, [α - , α + ] is the A-axis motion range of the machine tool, and (x 0 , y 0 , z 0 , α 0 , γ 0 ) is the origin position of the machine tool after tool setting.
进一步地,步骤S6具体为:Further, step S6 is specifically as follows:
S6.1,将振镜加工平面法向TM绕加工坐标系的ZM轴旋转到加工坐标系的+YMOMZM平面上;S6.1, rotate the normal T M of the galvanometer processing plane around the Z M axis of the processing coordinate system to the +Y M O M Z M plane of the processing coordinate system;
S6.2,通过下式,得到机床A轴向负向旋转时机床C轴运动分量CS2:S6.2, the motion component CS2 of the C-axis of the machine tool when the A-axis of the machine tool rotates in the negative direction is obtained by the following formula:
当振镜加工平面法向TM在YM轴上的分量jM=0时,机床C轴运动分量CS2为:When the component j M of the normal direction TM of the galvanometer processing plane on the Y M axis = 0, the motion component CS2 of the C axis of the machine tool is:
当振镜加工平面法向TM在YM轴上的分量jM≠0时,机床C轴运动分量CS2为:When the component j M of the normal direction TM of the galvanometer processing plane on the Y M axis ≠ 0, the motion component CS2 of the C axis of the machine tool is:
S6.3,通过下式得到变换后的矢量TC,M2:S6.3, the transformed vector T C,M2 is obtained by the following formula:
TC,M2=TM·R-C2 TC,M2 = TM ·R -C2
其中:in:
S6.4,通过下式获取机床A轴运动分量:S6.4, obtain the motion component of the machine tool A axis through the following formula:
当kC,M2=0时,机床A轴运动分量AS2=-π/2;When k C,M2 = 0, the motion component of the machine tool A axis AS2 = -π/2;
当kC,M2≠0时,机床A轴运动分量AS2为:When k C,M2 ≠0, the motion component AS2 of the machine tool A axis is:
其中,kC,M2为变换后的矢量TC,M2在ZM轴上的分量,jC,M2为变换后的矢量TC,M2在YM轴上的分量;in, k C,M2 is the component of the transformed vector T C,M2 on the Z M axis, j C,M2 is the component of the transformed vector T C,M2 on the Y M axis;
S6.5,通过下式计算机床在其X轴、Y轴、Z轴的平移分量XS2、YS2和ZS2:S6.5, calculate the translation components X S2 , Y S2 and Z S2 of the tool in its X-axis, Y-axis and Z-axis by the following formula:
其中,R-A2为变换后的矢量TC,M2到+ZM轴的旋转变换矩阵:Among them, R -A2 is the rotation transformation matrix of the transformed vector T C,M2 to the +Z M axis:
进一步地,步骤S7中,所述机床行程限制通过下式确定:Further, in step S7, the machine tool travel limit is determined by the following formula:
进一步地,步骤S8具体为,通过下式得到XG轴、YG轴在加工坐标系下的轴向:Furthermore, step S8 is specifically to obtain the axial directions of the XG axis and the YG axis in the machining coordinate system by the following formula:
其中,XG,M为XG轴在加工坐标系下的轴向,YG,M为YG轴在加工坐标系下的轴向,R-A为变换后的矢量TC,M1到+ZM轴的旋转变换矩阵R-A1或变换后的矢量TC,M2到+ZM轴的旋转变换矩阵R-A2,R-C为旋转变换矩阵R-C1或旋转变换矩阵R-C2。Among them, XG ,M is the axial direction of the XG axis in the machining coordinate system, YG ,M is the axial direction of the YG axis in the machining coordinate system, R -A is the rotation transformation matrix R -A1 of the transformed vector TC ,M1 to the + ZM axis or the rotation transformation matrix R -A2 of the transformed vector TC ,M2 to the + ZM axis, and R -C is the rotation transformation matrix R -C1 or the rotation transformation matrix R -C2 .
进一步地,步骤S9中,所述振镜扫描坐标系,是以当前加工特征对应的振镜加工平面原点PM为原心,XG,M轴为对应X轴,YG,M轴为对应Y轴建立的坐标系。Furthermore, in step S9, the galvanometer scanning coordinate system is a coordinate system established with the origin PM of the galvanometer processing plane corresponding to the current processing feature as the origin, the XG ,M axis as the corresponding X axis, and the YG ,M axis as the corresponding Y axis.
与现有技术相比,本发明具有以下有益效果:Compared with the prior art, the present invention has the following beneficial effects:
本发明提出一种独立光束扫描的五轴激光加工设备光束运动控制方法,通过控制带独立光束扫描装置的五轴激光加工设备各轴运动分量,保证光束的复合运动轨迹遵循预设路径,从而实现任意零件上指定形状特征的正确加工。通过分析五轴机床运动结构及连接关系,将光束扫描装置的运动轨迹分解到各运动轴上,获取各轴复合作用下的光束扫描装置相对工件的位姿,继而建立当前位姿下的局部扫描坐标系,实现带独立光束扫描装置的五轴激光加工设备运动控制。The present invention proposes a beam motion control method for a five-axis laser processing equipment with independent beam scanning, which ensures that the composite motion trajectory of the beam follows a preset path by controlling the motion components of each axis of the five-axis laser processing equipment with an independent beam scanning device, thereby realizing the correct processing of specified shape features on any part. By analyzing the motion structure and connection relationship of the five-axis machine tool, the motion trajectory of the beam scanning device is decomposed into each motion axis, the position and posture of the beam scanning device relative to the workpiece under the composite action of each axis is obtained, and then a local scanning coordinate system under the current position and posture is established to realize the motion control of the five-axis laser processing equipment with an independent beam scanning device.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
图1为XYZAC机械结构+2D振镜形式五轴激光加工设备示意图;Figure 1 is a schematic diagram of a five-axis laser processing equipment in the form of an XYZAC mechanical structure + 2D galvanometer;
图2为本发明一种独立光束扫描的五轴激光加工设备光束运动控制方法实施例中工件与振镜加工位姿示意图;FIG2 is a schematic diagram of the processing postures of a workpiece and a galvanometer mirror in an embodiment of a beam motion control method for a five-axis laser processing device with independent beam scanning according to the present invention;
图3为本发明一种独立光束扫描的五轴激光加工设备光束运动控制方法实施例中,矢量TM从加工坐标系的-YMOMXM平面旋转到-YMOMZM平面的变换关系示意图;3 is a schematic diagram of the transformation relationship of the vector TM from the -Y M O M X M plane of the processing coordinate system to the -Y M O M Z M plane in an embodiment of a beam motion control method for a five-axis laser processing equipment with independent beam scanning of the present invention;
图4为本发明一种独立光束扫描的五轴激光加工设备光束运动控制方法实施例中,jC,M1≤0时矢量TC,M1旋转到ZM轴的变换关系示意图;4 is a schematic diagram of the transformation relationship of the vector T C,M1 rotating to the Z M axis when j C,M1 ≤0 in an embodiment of a beam motion control method for a five-axis laser processing device with independent beam scanning according to the present invention;
图5为本发明一种独立光束扫描的五轴激光加工设备光束运动控制方法实施例中,矢量TM旋转到坐标系+YMOMZM平面的变换关系示意图;5 is a schematic diagram of the transformation relationship of the vector TM rotating to the +Y M O M Z M plane of the coordinate system in an embodiment of a beam motion control method for a five-axis laser processing device with independent beam scanning according to the present invention;
图6为本发明一种独立光束扫描的五轴激光加工设备光束运动控制方法实施例中,jC,M2≥0时矢量TC,M2旋转到ZM轴的变换关系示意图;6 is a schematic diagram of the transformation relationship of the vector T C, M2 rotating to the Z M axis when j C, M2 ≥ 0 in an embodiment of a beam motion control method for a five-axis laser processing device with independent beam scanning according to the present invention;
图7为本发明一种独立光束扫描的五轴激光加工设备光束运动控制方法实施例中,机床各轴运动后的位姿及原始图形下的振镜加工结果示意图;7 is a schematic diagram of the position and posture of each axis of the machine tool after movement and the galvanometer processing result under the original graphics in an embodiment of a beam motion control method for a five-axis laser processing equipment with independent beam scanning of the present invention;
图8为本发明一种独立光束扫描的五轴激光加工设备光束运动控制方法实施例中,机床中的振镜扫描坐标系示意图;FIG8 is a schematic diagram of a galvanometer scanning coordinate system in a machine tool in an embodiment of a beam motion control method for a five-axis laser processing device with independent beam scanning according to the present invention;
图9为本发明一种独立光束扫描的五轴激光加工设备光束运动控制方法实施例中,加工坐标OM-XMYMZM中的振镜坐标系PM-XG,MYG,M示意图。Fig. 9 is a schematic diagram of a galvanometer coordinate system PM - XG , MYG, M in a processing coordinate system OM- XMYMZM in an embodiment of a beam motion control method for a five-axis laser processing equipment with independent beam scanning of the present invention.
具体实施方式Detailed ways
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Generally, the components of the embodiments of the present invention described and shown in the drawings here can be arranged and designed in various different configurations.
本发明为解决多轴激光加工装备在加工时,若激光光束的运动图形不是圆形,会导致错位加工,造成零件报废的问题,提出一种独立光束扫描的五轴激光加工设备光束运动控制方法,如下是本发明的一个具体实施例,包括以下步骤:In order to solve the problem that when a multi-axis laser processing equipment is processing, if the motion pattern of the laser beam is not circular, it will cause misaligned processing and cause the parts to be scrapped, the present invention proposes a beam motion control method for a five-axis laser processing equipment with independent beam scanning. The following is a specific embodiment of the present invention, which includes the following steps:
如图1所示,以XYZAC机械运动结构+2D振镜构成的五轴激光加工设备为例,其中,A轴回转中心到C轴回转中心的偏移量为(Δx1,Δy1,Δz1),机床原点到A轴回转中心的偏移量为(Δx2,Δy2,Δz2),2D振镜的扫描平面坐标系OG-XGYG中XG、YG轴分别与机床的X轴、Y轴平行。机床X轴运动范围[x-,x+],Y轴运动范围[y-,y+],Z轴运动范围[z-,z+],A轴运动范围[α-,α+],C轴运动范围[0°,360°](支持连续回转),机床对刀后的原点位置为(x0,y0,z0,α0,γ0),其具体运动控制实现步骤如下:As shown in Figure 1, taking the five-axis laser processing equipment composed of XYZAC mechanical motion structure + 2D galvanometer as an example, the offset from the A-axis rotation center to the C-axis rotation center is ( Δx1 , Δy1 , Δz1 ), the offset from the machine tool origin to the A-axis rotation center is ( Δx2 , Δy2 , Δz2 ), and the XG and YG axes in the scanning plane coordinate system OG - XGYG of the 2D galvanometer are parallel to the X-axis and Y-axis of the machine tool , respectively. The machine tool X-axis motion range is [x - , x + ], the Y-axis motion range is [y - , y + ], the Z-axis motion range is [z - , z + ], the A-axis motion range is [α - , α + ], and the C-axis motion range is [0°, 360°] (supports continuous rotation). The origin position of the machine tool after tool setting is (x 0 , y 0 , z 0 , α 0 , γ 0 ). The specific motion control implementation steps are as follows:
步骤1:如图2,设加工坐标系OM-XMYMZM下,当前加工特征对应的振镜加工位姿SM=(PM,TM),其中PM=(xM,yM,zM)为振镜加工平面原点,TM=(iM,jM,kM)为振镜加工平面法向。Step 1: As shown in Figure 2, assume that in the processing coordinate system O M -X M Y M Z M , the galvanometer processing posture S M = (P M , T M ) corresponding to the current processing feature, where P M = (x M , y M , z M ) is the origin of the galvanometer processing plane, and T M = (i M , j M , k M ) is the normal of the galvanometer processing plane.
步骤2:如图3所示,为使机床A轴向正向旋转,应取A≥0,故需将矢量TM从加工坐标系的-YMOMXM平面绕ZM轴旋转到加工坐标系的-YMOMZM平面上,从而获取机床C轴运动分量。其中,当TM在YM轴上的分量jM=0时,机床C轴运动分量CS1可由下式计算:Step 2: As shown in Figure 3, in order to make the machine tool A axis rotate in the positive direction, A ≥ 0 should be taken, so the vector TM needs to be rotated from the -Y M O M X M plane of the machining coordinate system around the Z M axis to the -Y M O M Z M plane of the machining coordinate system, so as to obtain the machine tool C axis motion component. Among them, when the component j M of TM on the Y M axis = 0, the machine tool C axis motion component CS1 can be calculated by the following formula:
其中,iM为TM在XM轴上的分量;Where i M is the component of T M on the X M axis;
当jM≠0时,机床C轴运动分量CS1可由下式计算:When j M ≠ 0, the C-axis motion component of the machine tool, C S1 , can be calculated by the following formula:
其中, in,
步骤3:矢量TM到-YMOMZM平面的旋转变换矩阵R-C1可表示为:Step 3: The rotation transformation matrix R -C1 of the vector TM to the -YMOMZM plane can be expressed as:
则利用R-C1可将矢量TM变换为矢量TC,M1=(iC,M1,jC,M1,kC,M1),即:Then, using R -C1, the vector TM can be transformed into a vector TC ,M1 = (i C,M1 , j C,M1 , k C,M1 ), that is:
TC,M1=TM·R-C1 TC,M1 = TM ·R -C1
其中,必然存在jC,M1≤0。Among them, there must exist j C,M1 ≤0.
步骤4:如图4所示,在YMOMZM平面上将矢量TC,M1绕XM轴旋转到+ZM轴上,即可得到机床的A轴运动分量。其中,当kC,M1=0时,机床A轴运动分量AS1=π/2;当kC,M1≠0时,机床A轴运动分量AS1可由下式计算:Step 4: As shown in Figure 4, rotate the vector T C, M1 around the X M axis to the + Z M axis on the Y M O M Z M plane to obtain the A-axis motion component of the machine tool. When k C, M1 = 0, the A-axis motion component of the machine tool AS1 = π/2; when k C, M1 ≠ 0, the A-axis motion component of the machine tool AS1 can be calculated by the following formula:
其中/> Where/>
步骤5:矢量TC,M1到+ZM轴的旋转变换矩阵R-A1可表示为:Step 5: The rotation transformation matrix R -A1 of the vector TC,M1 to the + ZM axis can be expressed as:
步骤6:在得到机床两个旋转轴运动分量CS1和AS1后,可通过下式计算机床在X轴、Y轴、Z轴的三个平移分量XS1、YS1和ZS1:Step 6: After obtaining the motion components CS1 and AS1 of the two rotation axes of the machine tool, the three translation components XS1 , YS1 and ZS1 of the machine tool in the X-axis, Y-axis and Z-axis can be calculated by the following formula:
步骤7:由于机床C轴可连续旋转,因此C轴不会超出机床行程。机床X、Y、Z、A轴受行程限制,需要利用以下条件判定当前各轴是否超过机床行程限制,即:Step 7: Since the C axis of the machine tool can rotate continuously, the C axis will not exceed the machine tool travel. The X, Y, Z, and A axes of the machine tool are limited by travel. The following conditions need to be used to determine whether the current axes exceed the machine tool travel limit, namely:
当各轴运动分量满足机床行程条件时,转到步骤11;否则,转到步骤8。When the motion components of each axis meet the machine tool travel conditions, go to step 11; otherwise, go to step 8.
步骤8:如图5所示,使机床A轴向负向旋转,即取A<0,故需将矢量TM绕ZM轴旋转到坐标系的+YMOMZM平面上,从而获取机床的C轴运动分量。Step 8: As shown in Figure 5, make the machine tool A-axis rotate in the negative direction, that is, take A<0, so it is necessary to rotate the vector TM around the ZM axis to the + YMOMZM plane of the coordinate system to obtain the C-axis motion component of the machine tool.
其中,当jM=0时,机床C轴运动分量CS2可由下式计算When j M = 0, the C-axis motion component CS2 of the machine tool can be calculated by the following formula:
当jM≠0时,机床C轴运动分量CS2可由下式计算When j M ≠ 0, the C-axis motion component CS2 of the machine tool can be calculated by the following formula:
其中/> Where/>
则矢量TM到+YMOMZM平面的旋转变换矩阵Then the rotation transformation matrix of vector TM to + YMOMZM plane is
矢量TM经过R-C2变换后的得到TC,M2(iC,M2,jC,M2,kC,M2)=TM·R-C2,则必有jC,M2≥0。After R -C2 transformation, the vector TM obtains TC,M2 (iC ,M2 ,jC ,M2 ,kC ,M2 )= TM ·R -C2 , then jC , M2≥0.
步骤9:如图6所示,在YMOMZM平面上将矢量TC,M2绕XM轴旋转到+ZM轴上,即可得到机床的A轴运动分量。其中,当kC,M2=0时,机床A轴运动分量AS2=-π/2;当kC,M2≠0时,机床A轴运动分量AS2可由下式计算:Step 9: As shown in Figure 6, rotate the vector T C, M2 around the X M axis to the + Z M axis on the Y M O M Z M plane to obtain the A-axis motion component of the machine tool. Among them, when k C, M2 = 0, the A-axis motion component of the machine tool AS2 = -π/2; when k C, M2 ≠ 0, the A-axis motion component of the machine tool AS2 can be calculated by the following formula:
其中/> Where/>
则矢量TC,M2到+ZM轴的旋转变换矩阵R-A2可表示为:Then the rotation transformation matrix R -A2 of the vector T C,M2 to the +Z M axis can be expressed as:
故机床的三个平移分量XS2、YS2和ZS2可通过下式计算:Therefore, the three translation components of the machine tool X S2 , Y S2 and Z S2 can be calculated by the following formula:
步骤10:判断机床X、Y、Z、A轴运动分量是否满足以下限制Step 10: Determine whether the motion components of the machine tool X, Y, Z, and A axes meet the following restrictions
当满足时,转到步骤11;否则,可判定当前设备不能用于加工该零件,需要更换零件装夹位置或更换机床。If the conditions are met, go to step 11; otherwise, it can be determined that the current equipment cannot be used to process the part, and the part clamping position needs to be changed or the machine tool needs to be replaced.
步骤11:如图7所示,机床各轴经过上述旋转和平移运动后,如果直接利用振镜进行原始图形的扫描,则会导致加工出的特征发生旋转而得到错误的形状。为处理这一问题,令(R-A,R-C)为机床A轴向正向旋转过程中的变换矩阵(R-A1,R-C1)或A轴向负向旋转过程中的变换矩阵(R-A2,R-C2),(XS,YS,ZS,AS,CS)为机床A轴取正时各轴的当前运动分量(XS1,YS1,ZS1,AS1,CS1)或A轴取负时各轴的当前运动分量(XS2,YS2,ZS2,AS2,CS2),并转到步骤12。Step 11: As shown in FIG7 , after the axes of the machine tool undergo the above rotation and translation motion, if the galvanometer is directly used to scan the original graphics, the processed features will rotate and the wrong shape will be obtained. To deal with this problem, let (R -A , R -C ) be the transformation matrix (R -A1 , R -C1 ) during the positive rotation of the A axis of the machine tool or the transformation matrix (R -A2 , R -C2 ) during the negative rotation of the A axis, ( XS , YS , ZS , AS , CS ) be the current motion components of each axis when the A axis of the machine tool takes the positive direction ( XS1 , YS1 , ZS1 , AS1 , CS1 ) or the current motion components of each axis when the A axis takes the negative direction ( XS2 , YS2 , ZS2 , AS2 , CS2 ), and go to step 12.
步骤12:如图8所示,以机床各轴运动后的位姿为参考,将当前位姿下的X轴运动方向(1,0,0)作为振镜扫描图形的XG轴方向,当前Y轴运动方向(0,1,0)作为振镜扫描图形的YG轴方向,利用以下公式计算XG轴、YG轴在加工坐标系OM-XMYMZM下的轴向XG,M、YG,M:Step 12: As shown in Figure 8, taking the position and posture of each axis of the machine tool after movement as a reference, the X-axis movement direction (1, 0, 0) in the current position and posture is used as the XG- axis direction of the galvanometer scanning pattern, and the current Y-axis movement direction (0, 1, 0) is used as the YG- axis direction of the galvanometer scanning pattern. The axial directions XG , M and YG , M of the XG and YG axes in the machining coordinate system OM - XMYMZM are calculated using the following formula:
步骤13:以XG,M方向为振镜扫描平面X轴,YG,M方向为振镜扫描平面Y轴,当前加工特征对应的振镜加工平面原点PM为原心,建立如图9所示的振镜扫描坐标系PM-XG,MYG,M,并在该坐标系中生成振镜扫描轨迹,结合机床各轴运动分量(XS,YS,ZS,AS,CS),即可得到带独立光束扫描装置的五轴激光加工设备完整运动参数,从而实现零件特征的准确加工。Step 13: With the XG ,M direction as the X-axis of the galvanometer scanning plane, the YG ,M direction as the Y-axis of the galvanometer scanning plane, and the origin PM of the galvanometer processing plane corresponding to the current processing feature as the origin, establish the galvanometer scanning coordinate system PM - XG,MYG ,M as shown in Figure 9, and generate the galvanometer scanning trajectory in the coordinate system. Combined with the motion components of each axis of the machine tool ( XS , YS , ZS , AS , CS ), the complete motion parameters of the five-axis laser processing equipment with an independent beam scanning device can be obtained, thereby realizing accurate processing of part features.
以上仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The above are only preferred embodiments of the present invention and are not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.
Claims (9)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310336954.9A CN116339238B (en) | 2023-03-31 | 2023-03-31 | A beam motion control method for five-axis laser processing equipment with independent beam scanning |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310336954.9A CN116339238B (en) | 2023-03-31 | 2023-03-31 | A beam motion control method for five-axis laser processing equipment with independent beam scanning |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116339238A CN116339238A (en) | 2023-06-27 |
CN116339238B true CN116339238B (en) | 2024-04-05 |
Family
ID=86894738
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310336954.9A Active CN116339238B (en) | 2023-03-31 | 2023-03-31 | A beam motion control method for five-axis laser processing equipment with independent beam scanning |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116339238B (en) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008185607A (en) * | 2007-01-26 | 2008-08-14 | Mitsubishi Electric Corp | Laser processing apparatus and laser processing method |
CN101424512A (en) * | 2008-12-18 | 2009-05-06 | 浙江大学 | Method for detecting radial motion error of high speed principal axis by multi-ring coincidence three-point method |
CN106354094A (en) * | 2016-09-14 | 2017-01-25 | 大连理工大学 | Method for calibrating machine tool follow-up laser scanning coordinates on basis of space standard balls |
CN107121062A (en) * | 2016-12-07 | 2017-09-01 | 苏州笛卡测试技术有限公司 | A kind of robot three-dimensional scanning means and method |
CN109765848A (en) * | 2019-03-18 | 2019-05-17 | 西南交通大学 | Calculation and impact assessment method of geometric error contribution value of motion axis of five-axis machine tool |
CN110385522A (en) * | 2019-08-07 | 2019-10-29 | 西安中科微精光子制造科技有限公司 | Laser machine control method and control system |
CN110500978A (en) * | 2019-09-25 | 2019-11-26 | 上海交通大学 | On-line Calibration Method of Beam Direction Vector and Zero Point Position of Point Laser Sensor |
CN112325773A (en) * | 2020-10-30 | 2021-02-05 | 上海交通大学 | Method for calibrating beam direction vector and origin position of laser displacement sensor |
WO2021179460A1 (en) * | 2020-03-08 | 2021-09-16 | 华中科技大学 | Laser light exit direction calibration method employing standard ball |
CN113427133A (en) * | 2021-06-16 | 2021-09-24 | 西安交通大学 | Laser equipment and method for guiding automatic processing based on three-dimensional vision online measurement |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9063533B2 (en) * | 2009-03-27 | 2015-06-23 | Mitsubishi Electric Corporation | Multi-spindle translation control for multiple coordinate systems |
WO2019090487A1 (en) * | 2017-11-07 | 2019-05-16 | 大连理工大学 | Highly dynamic wide-range any-contour-error monocular six-dimensional measurement method for numerical control machine tool |
US10684609B2 (en) * | 2018-10-30 | 2020-06-16 | Hurco Automation, Ltd. | Program code generating method for tilted plane machining by multi-axis machine tool and device thereof |
-
2023
- 2023-03-31 CN CN202310336954.9A patent/CN116339238B/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008185607A (en) * | 2007-01-26 | 2008-08-14 | Mitsubishi Electric Corp | Laser processing apparatus and laser processing method |
CN101424512A (en) * | 2008-12-18 | 2009-05-06 | 浙江大学 | Method for detecting radial motion error of high speed principal axis by multi-ring coincidence three-point method |
CN106354094A (en) * | 2016-09-14 | 2017-01-25 | 大连理工大学 | Method for calibrating machine tool follow-up laser scanning coordinates on basis of space standard balls |
CN107121062A (en) * | 2016-12-07 | 2017-09-01 | 苏州笛卡测试技术有限公司 | A kind of robot three-dimensional scanning means and method |
CN109765848A (en) * | 2019-03-18 | 2019-05-17 | 西南交通大学 | Calculation and impact assessment method of geometric error contribution value of motion axis of five-axis machine tool |
CN110385522A (en) * | 2019-08-07 | 2019-10-29 | 西安中科微精光子制造科技有限公司 | Laser machine control method and control system |
CN110500978A (en) * | 2019-09-25 | 2019-11-26 | 上海交通大学 | On-line Calibration Method of Beam Direction Vector and Zero Point Position of Point Laser Sensor |
WO2021179460A1 (en) * | 2020-03-08 | 2021-09-16 | 华中科技大学 | Laser light exit direction calibration method employing standard ball |
CN112325773A (en) * | 2020-10-30 | 2021-02-05 | 上海交通大学 | Method for calibrating beam direction vector and origin position of laser displacement sensor |
CN113427133A (en) * | 2021-06-16 | 2021-09-24 | 西安交通大学 | Laser equipment and method for guiding automatic processing based on three-dimensional vision online measurement |
Non-Patent Citations (2)
Title |
---|
"十轴六联动双光束激光焊接及其 跟踪控制技术";龚时华 等;《航空制造技术》;第61卷(第11期);34-39 * |
"激光微加工系统中运动控制的实现";裴红星 等;《组合机床与自动化加工技术》(第4期);53-57 * |
Also Published As
Publication number | Publication date |
---|---|
CN116339238A (en) | 2023-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9002502B2 (en) | Method and program for calculating correction value for machine tool | |
JP4083554B2 (en) | 3D model data generator | |
CN107507177B (en) | Processing of robots object localization method based on 3-D scanning and device | |
WO2022170841A1 (en) | Benchmark selection method for minimizing hole position errors in large-component hole group machining | |
CN110716497B (en) | Registration method based on plane reference constraint and margin constraint | |
CN107999967A (en) | A kind of large-scale three dimensional piece surface parallel laser lithography method and apparatus | |
Zhang et al. | Singularity-aware motion planning for multi-axis additive manufacturing | |
CN107775065A (en) | The coordinated movement of various economic factors synchronisation control means of the wall thickness such as dual robot mirror image milling processing | |
US20130204426A1 (en) | Tool Path Generation Method and Apparatus | |
CN113119450A (en) | AC corner optimization method for five-axis 3D printing line cutting path planning | |
CN107097122A (en) | A kind of robot for independently grinding large-scale free form surface | |
KR20030055250A (en) | Three-dimensional laser beam machine | |
CN116339238B (en) | A beam motion control method for five-axis laser processing equipment with independent beam scanning | |
CN118192418A (en) | A nonlinear error compensation method for non-standard five-axis linkage machine tool machining | |
Garmendia et al. | Geometrical control of DED processes based on 3D scanning applied to the manufacture of complex parts | |
CN107942942B (en) | Inclined coordinate system establishing method applied to intersected inclined planes of machine tool equipment | |
CN111992909B (en) | Three-dimensional laser drilling positioning method | |
CN118295321A (en) | Method for generating finish machining tool path of open-pore groove on surface of thin-wall part based on milling process | |
Yun et al. | A geometric postprocessing method for 5-axis machine tools using locations of joint points | |
CN114001683B (en) | Complex component machining residual wall thickness self-adaptive compensation method based on coordinate origin translation | |
JP2005349520A (en) | Wire-cut electric-discharge machining method, machining program producing method in wire-cut electric-discharge machining method, and wire cut electric-discharge machining program producing device | |
CN114739290A (en) | Path planning method and system for line laser scanning chemical milling of rubber reticle | |
CN112419415B (en) | Ultrasonic scanning method for realizing pose planning based on CRS curve fitting | |
CN112846323B (en) | Three-dimensional vibration-assisted milling processing system and three-dimensional vibration-assisted milling method for structural surface | |
Cai et al. | Trajectory planning and simulation for intersecting line cutting of the industry robot |
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 |