CN102490081A - Workpiece three-dimensional surface topography simulating method based on ball head milling - Google Patents
Workpiece three-dimensional surface topography simulating method based on ball head milling Download PDFInfo
- Publication number
- CN102490081A CN102490081A CN2011103595763A CN201110359576A CN102490081A CN 102490081 A CN102490081 A CN 102490081A CN 2011103595763 A CN2011103595763 A CN 2011103595763A CN 201110359576 A CN201110359576 A CN 201110359576A CN 102490081 A CN102490081 A CN 102490081A
- Authority
- CN
- China
- Prior art keywords
- surface topography
- dimensional surface
- tool
- point
- workpiece
- 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
Links
- 238000012876 topography Methods 0.000 title claims abstract description 87
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000003801 milling Methods 0.000 title claims abstract description 14
- 238000004088 simulation Methods 0.000 claims abstract description 55
- 238000012545 processing Methods 0.000 claims abstract description 28
- 238000005520 cutting process Methods 0.000 claims abstract description 16
- 238000003754 machining Methods 0.000 claims abstract description 14
- 238000005070 sampling Methods 0.000 claims abstract description 5
- 238000010408 sweeping Methods 0.000 claims abstract 2
- 239000011159 matrix material Substances 0.000 claims description 17
- 239000000284 extract Substances 0.000 claims description 15
- 230000009466 transformation Effects 0.000 claims description 8
- 239000013598 vector Substances 0.000 claims description 6
- 238000013519 translation Methods 0.000 claims description 5
- 238000005457 optimization Methods 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 12
- 238000000605 extraction Methods 0.000 description 6
- 238000012512 characterization method Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000000227 grinding Methods 0.000 description 3
- 230000003746 surface roughness Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000007514 turning Methods 0.000 description 1
Images
Landscapes
- Numerical Control (AREA)
Abstract
Description
技术领域 technical field
本发明涉及超精密加工领域,具体为一种工件三维表面形貌的建模及仿真方法,尤其适用于采用超精密球头铣刀加工的工件。The invention relates to the field of ultra-precision machining, in particular to a method for modeling and simulating three-dimensional surface topography of workpieces, and is especially suitable for workpieces processed by ultra-precision ball end milling cutters.
背景技术 Background technique
随着国防、航空航天、能源、医疗、光学元器件等技术和相关行业的发展,越来越多的基础装备对一些关键零部件如宇航陀螺,计算机磁鼓、磁盘,多面棱镜,大直径非球面镜,以及复杂形状的立体棱镜等提出了更高的要求。这类元器件对加工精度、表面粗糙度和三维表面形貌分布要求极高,使用常规的磨削、研磨、抛光等方法进行加工,不但加工成本很高,而且难以同时满足精度和表面粗糙度的要求,普通机床更是难以满足如此高的加工要求,必须采用超精密机床才能够完成加工。目前在加工三维表面形貌的质量成为制造业关注热点的形势下,各个部门及研究机构对超精密加工零件的表面形状精度、波纹度、即表面粗糙度等三维表面形貌的要求越来越高,相关研究也在不断进行。With the development of national defense, aerospace, energy, medical, optical components and other technologies and related industries, more and more basic equipment is required for some key components such as aerospace gyroscopes, computer drums, disks, polygonal prisms, large-diameter Spherical mirrors, and stereoscopic prisms with complex shapes put forward higher requirements. Such components have extremely high requirements on processing accuracy, surface roughness and three-dimensional surface topography distribution. Using conventional grinding, grinding, polishing and other methods for processing, not only the processing cost is high, but also it is difficult to meet the accuracy and surface roughness at the same time. It is even more difficult for ordinary machine tools to meet such high processing requirements, and ultra-precision machine tools must be used to complete the processing. At present, under the situation that the quality of processed three-dimensional surface topography has become a hot spot in the manufacturing industry, various departments and research institutions have more and more requirements for three-dimensional surface topography such as surface shape accuracy, waviness, and surface roughness of ultra-precision machined parts. High, related research is also ongoing.
近半个多世纪以来,研究人员不断的尝试通过评定工件的三维表面形貌来指导加工工艺的形成过程来满足加工要求,试图找到三维表面形貌形成与加工策略及加工参数的关系,从而有效的对工件三维表面形貌的形成过程在一定程度上进行控制。目前对三维表面形貌的研究主要从两个方面进行:1)实验测量工件三维表面形貌;2)理论预测工件三维表面形貌。虽然从这两个角度出发对三维表面形貌形成机理以及三维表面形貌评定的研究已经取得了一定成果,但是缺陷和问题仍然存在,问题主要集中在测量结果分析和三维表面形貌算法上。由于工件三维表面形貌的测量结果主要是对工件表面质量的评定,所以不能清楚表征加工参数以及加工策略对三维表面形貌的影响,仍主要出于定性分析。For more than half a century, researchers have been trying to guide the formation process of the machining process by evaluating the three-dimensional surface topography of the workpiece to meet the processing requirements, trying to find the relationship between the formation of the three-dimensional surface topography and the processing strategy and processing parameters, so as to effectively The formation process of the three-dimensional surface topography of the workpiece is controlled to a certain extent. At present, the research on the three-dimensional surface topography is mainly carried out from two aspects: 1) experimentally measuring the three-dimensional surface topography of the workpiece; 2) theoretically predicting the three-dimensional surface topography of the workpiece. Although some achievements have been made in the research on the formation mechanism of 3D surface topography and the evaluation of 3D surface topography from these two perspectives, defects and problems still exist. The problems mainly focus on the analysis of measurement results and the algorithm of 3D surface topography. Since the measurement results of the three-dimensional surface topography of the workpiece are mainly for the evaluation of the surface quality of the workpiece, the inability to clearly characterize the influence of processing parameters and processing strategies on the three-dimensional surface topography is still mainly due to qualitative analysis.
在超精密车削表面形貌建模方面,国内外做了较多工作。如:香港理工大学的李荣彬教授,加拿大的P.A.Meyer等,天津大学的陈东祥对超精密磨削的表面形貌建模进行了研究。但在超精密铣削加工表面形貌建模方面,研究较少。在球头铣削普通加工表面形貌仿真方面,很多学者也进行了研究,如西北工大的谭刚,山东大学的张孝峰等。但其针对的对象均是普通铣削加工,不能满足超精密铣削加工中的高精度或者在保证高精度时计算时间过长等问题。该方法从材料去除角度描述表面形貌生成过程,并且根据该过程进行切削力分析,但在计算交点过程中刀刃每变换一次空间位置就需要计算一组交点,计算量较大。In the aspect of ultra-precision turning surface topography modeling, a lot of work has been done at home and abroad. Such as: Professor Li Rongbin of Hong Kong Polytechnic University, P.A.Meyer of Canada, etc., and Chen Dongxiang of Tianjin University have studied the surface topography modeling of ultra-precision grinding. However, there are few studies on surface topography modeling in ultra-precision milling. Many scholars have also conducted research on the simulation of surface topography in general machining of ball-end milling, such as Tan Gang from Northwestern University of Technology, Zhang Xiaofeng from Shandong University, etc. However, it is aimed at ordinary milling, which cannot meet the high precision of ultra-precision milling or the calculation time is too long when high precision is guaranteed. This method describes the surface topography generation process from the perspective of material removal, and analyzes the cutting force according to the process. However, in the process of calculating the intersection point, a set of intersection points needs to be calculated every time the cutting edge changes its spatial position, which requires a large amount of calculation.
发明内容 Contents of the invention
本发明的目的在于提供一种基于超精密球头铣刀加工的工件三维表面形貌建模及仿真方法。通过对刀刃在加工过程中的运动学描述,将三维表面形貌用离散点云数据表达。根据加工曲面信息和采样点数目对仿真区域进行划分并建立随动包容盒,对包容盒内点云数据进行数值分析和空间变换,从而获取加工三维表面形貌。本方法解决了以往三维表面形貌仿真方法中计算精度低的缺点,非常适合于高精度的超精密加工,而且在高精度的要求下,仍能保持不错的计算效率。The object of the present invention is to provide a method for modeling and simulating three-dimensional surface topography of workpieces based on ultra-precision ball end milling cutter processing. Through the kinematic description of the cutting edge during the machining process, the three-dimensional surface topography is expressed by discrete point cloud data. According to the processing surface information and the number of sampling points, the simulation area is divided and a dynamic containment box is established, and the point cloud data in the containment box is numerically analyzed and spatially transformed to obtain the processed 3D surface topography. This method solves the shortcomings of low calculation accuracy in the previous three-dimensional surface topography simulation methods, and is very suitable for high-precision ultra-precision machining, and can still maintain good calculation efficiency under the high-precision requirements.
一种工件三维表面形貌的建模及仿真方法,具体过程如下:A method for modeling and simulating a three-dimensional surface topography of a workpiece, the specific process is as follows:
(一)在加工表面上选取任意区域Ssim(u,v)作为仿真区域,其中,u、v为参数域,u∈[0,1],v∈[0,1]。(1) Select any area S sim (u, v) on the machined surface as the simulation area, where u, v are parameter domains, u∈[0,1], v∈[0,1].
(二)根据仿真区域边界信息,选取已有加工刀具路径中影响仿真区域表面形貌的刀位点作为三维表面形貌仿真的刀位点,并记录每行刀路的行数以及刀路的起始刀位点;(2) According to the boundary information of the simulation area, select the tool position points that affect the surface topography of the simulation area in the existing machining tool path as the tool position points for the three-dimensional surface topography simulation, and record the number of lines of each line of tool paths and the number of tool paths starting tool position;
(三)以步骤(二)所选取的刀位点作为运动节点,根据刀刃运动形成三维表面形貌扫掠点云模型。刀位轨迹中任意位置的刀刃扫掠点云模型表达式为:(3) Using the knife position selected in step (2) as the motion node, a three-dimensional surface topography sweep point cloud model is formed according to the blade motion. The expression of the blade sweep point cloud model at any position in the tool trajectory is:
其中,B0为刀具坐标系{TCS}下刀具静态刃部离散点坐标矩阵,Bi为工件坐标系{WCS}下刀具运动后的刀刃离散点坐标矩阵,T1(t,N,θ)为刀具坐标系下绕自身Z轴旋转的4×4变换矩阵,T2(t,p0,p1,f)为刀具坐标系相对于工件坐标系平移的4×4变换矩阵,ti为从p0点运动到p0、p1点之间任意点pi所经历时间,N为刀具转速,θ为初始切入相位角,f为刀具进给速度,p0为三维表面形貌仿真起始刀位点,p1为三维表面形貌仿真的终止刀位点;Among them, B 0 is the coordinate matrix of the discrete points of the static edge of the tool in the tool coordinate system {TCS}, B i is the coordinate matrix of the discrete points of the blade after the tool moves in the workpiece coordinate system {WCS}, T 1 (t, N, θ) is the 4×4 transformation matrix that rotates around its own Z axis in the tool coordinate system, T 2 (t, p 0 , p 1 , f) is the 4×4 transformation matrix that translates the tool coordinate system relative to the workpiece coordinate system, and t i is The time elapsed from point p 0 to any point p i between points p 0 and p 1 , N is the tool rotation speed, θ is the initial cut-in phase angle, f is the tool feed rate, p 0 is the start of three-dimensional surface topography simulation The starting knife position, p 1 is the ending knife position of the three-dimensional surface topography simulation;
将影响仿真区域所有刀位点产生的刀刃离散点云Bi叠加在一起就构成了仿真区域刀刃扫掠点云模型。The point cloud model of the blade sweep in the simulation area is formed by superimposing the blade discrete point cloud B i generated by all the knife position points in the simulation area.
(四)三维表面形貌控制点提取:(4) Extraction of three-dimensional surface topography control points:
设定仿真区域三维表面形貌采样点数pu、pv,其中pu为u向采样数目,pv为v向采样数目,进而对仿真区域进行网格划分,生成网格节点Pknot(m,n),其中m、n为整数且m∈[1,pu-1],n∈[1,pv-1];Set the number of sampling points p u and p v of the three-dimensional surface topography in the simulation area, where p u is the number of samples in the u direction, and p v is the number of samples in the v direction, and then mesh the simulation area to generate grid nodes P knot (m , n), where m, n are integers and m∈[1, p u -1], n∈[1, p v -1];
三维表面形貌的表征是由Pknot(m,n)节点处的表征点Ps(m,n)来控制。在仿真区域网格上建立包容盒,包容盒边界方向与节点Pknot(m,n)在加工曲面处法失方向平行,包容盒边界与对应网格边界相同,根据此包容盒对步骤(三)中仿真区域刀刃扫掠点云模型进行筛选,并找出筛选点中离网格距离最近点作为该网格区域三维表面形貌表征点Ps(m,n)。The characterization of the three-dimensional surface topography is controlled by the characteristic point P s (m, n) at the node of P knot (m, n). Establish a containment box on the grid of the simulation area, the boundary direction of the containment box is parallel to the normal direction of the node P knot (m, n) at the processing surface, and the boundary of the containment box is the same as the corresponding grid boundary. According to the containment box, the steps (three ) to filter the point cloud model of the knife edge sweep in the simulation area, and find out the point closest to the grid among the screened points as the 3D surface topography characterization point P s (m, n) of the grid area.
最后根据上述方法遍历仿真区域网格实现仿真区域上三维表面形貌控制点的提取,并最终得到工件表面的三维表面形貌;Finally, according to the above method, the grid of the simulation area is traversed to realize the extraction of the three-dimensional surface topography control points on the simulation area, and finally the three-dimensional surface topography of the workpiece surface is obtained;
本发明以三维表面形貌生成机理及多轴数控加工理论为基础,建立了一种结合加工工艺参数、刀具运动策略以及刀具模型的工件三维表面形貌仿真模型,该模型充分考虑到了机床运动、刀具运动以及三维表面形貌生成的关系;并在此基础之上建立了一种独立于机床结构的三维表面形貌仿真算法,能够很好的表征超精密球头铣刀加工条件下工件三维表面形貌特征。Based on the three-dimensional surface topography generation mechanism and multi-axis numerical control machining theory, the present invention establishes a workpiece three-dimensional surface topography simulation model that combines processing parameters, tool movement strategies and tool models. The model fully takes into account machine tool movement, The relationship between tool movement and 3D surface topography generation; and on this basis, a 3D surface topography simulation algorithm independent of the machine tool structure is established, which can well characterize the three-dimensional surface of the workpiece under the processing conditions of ultra-precision ball end milling cutters Morphological features.
本发明根据已设计好的刀具路径及工艺参数加工出的工件,通过表面轮廓测量仪对其表面进行测量获得三维表面形貌测量值,然后利用本发明中的仿真方法对其加工表面进行仿真,并最终与实验测量值进行对比,能够很好的表征工件三维表面形貌的各个特征。本发明充分考虑了机床加工参数以及刀具运动策略对三维表面形貌的影响,可以根据已有的刀具路径和加工参数对工件三维表面形貌进行建模及仿真,为工艺人员提供优化加工参数和刀具运动策略的方法和手段。In the present invention, according to the workpiece processed by the designed tool path and process parameters, the surface is measured by the surface profile measuring instrument to obtain the three-dimensional surface topography measurement value, and then the simulation method in the present invention is used to simulate the processed surface. And finally compared with the experimental measurement value, it can well characterize the various characteristics of the three-dimensional surface topography of the workpiece. The present invention fully considers the influence of machine tool processing parameters and tool motion strategies on the three-dimensional surface topography, and can model and simulate the three-dimensional surface topography of the workpiece according to the existing tool paths and processing parameters, and provide technicians with optimized processing parameters and Methods and means for tool motion strategies.
附图说明 Description of drawings
图1为提取参与三维表面形貌生成刀位点示意图Figure 1 is a schematic diagram of the extraction of knife positions involved in the generation of 3D surface topography
图2为刀具在机床运动过程示意图Figure 2 is a schematic diagram of the tool movement process in the machine tool
图3为提取三维表面形貌刀刃扫掠点示意图Figure 3 is a schematic diagram of the blade sweep points for extracting 3D surface topography
图4为三维表面形貌特征点提取示意图Figure 4 is a schematic diagram of the extraction of three-dimensional surface topography feature points
具体实施方式 Detailed ways
下面结合附图和具体实施例对本发明作进一步说明。The present invention will be further described below in conjunction with the accompanying drawings and specific embodiments.
本发明具体实施步骤如下:The specific implementation steps of the present invention are as follows:
(1)设定仿真区域并提取影响仿真区域三维表面形貌的刀位点。为了减少计算量,只需要将影响选定仿真区域Ssim(u,v)的刀位点提取出来进行建模即可,这样既可以减少无谓的计算也可以满足加工表面三维表面形貌的特征。(1) Set the simulation area and extract the tool points that affect the three-dimensional surface topography of the simulation area. In order to reduce the amount of calculation, it is only necessary to extract the tool position points that affect the selected simulation area S sim (u, v) for modeling, which can reduce unnecessary calculations and meet the characteristics of the three-dimensional surface topography of the machined surface .
如图1所示,在三轴加工条件下,为不失一般性,选取工件表面的任意区域作为仿真区域Ssim(u,v),并记录其边界:As shown in Fig. 1, under the condition of three-axis machining, without loss of generality, any area on the workpiece surface is selected as the simulation area S sim (u, v), and its boundary is recorded:
E1,E2,E3和E4分别为仿真区域Ssim(u,v)的四条边界曲线。E 1 , E 2 , E 3 and E 4 are four boundary curves of the simulation area S sim (u, v) respectively.
由于刀位点在刀位点偏置面上,而目前的边界是在加工曲面上选取的,所以需要对该边界进行相应偏置计算才可以作为选题三维表面形貌建模刀位点的条件。偏置方向为边界上各点在加工曲面上该点的曲面法失方向,偏置距离为刀具半径R,故有:Since the tool point is on the offset surface of the tool point, and the current boundary is selected on the processing surface, it is necessary to perform corresponding offset calculations on the boundary before it can be used as the tool point of the selected topic for 3D surface topography modeling. condition. The offset direction is the surface normal direction of each point on the processing surface at the point on the boundary, and the offset distance is the tool radius R, so:
E1off,E2off,E3off和E4off分别为仿真区域Ssim(u,v)的四条边界曲线经偏置后的边界曲线。E 1off , E 2off , E 3off and E 4off are respectively the offset boundary curves of the four boundary curves of the simulation region S sim (u, v).
根据偏置后的边界曲线所形成的仿真区域进行刀位点提取,找出在区域范围内的刀位点作为三维表面形貌建模的刀位点,并且按照走刀顺序进行保存,同时记录各条刀具路径在仿真区域内的起始点和终止点,用作后续切入以及切出相位角的计算。According to the simulation area formed by the offset boundary curve, the tool position points are extracted, and the tool position points within the area are found as the tool position points for 3D surface topography modeling, and are saved according to the order of the cutting tool, and recorded at the same time The start and end points of each tool path in the simulation area are used for the calculation of the subsequent cut-in and cut-out phase angles.
(2)根据(1)中提取出的刀位点,建立刀刃扫掠点云模型,如图2所示。具体步骤如下:(2) According to the knife position points extracted in (1), establish the point cloud model of the blade sweep, as shown in Figure 2. Specific steps are as follows:
(2.1)刀具旋转矩阵计算。设定刀轴与刀具坐标系Z轴重合,则刀具坐标系下旋转矩阵T1(ti,N,θ)的计算公式如下:(2.1) Calculation of tool rotation matrix. Assuming that the tool axis coincides with the Z axis of the tool coordinate system, the calculation formula of the rotation matrix T 1 (t i , N, θ) in the tool coordinate system is as follows:
(2.2)刀具平移矩阵计算。加工过程中的刀位文件坐标是相对于工件坐标系生成的,为了使计算方便,在建模初始状态下将刀具坐标系与工件坐标系重合,则刀具在任一两刀位点之间的移动过程可根据已生成的刀位文件计算得到。三轴加工中刀具在两刀位点之间的移动过程为直线插补,即两刀位点之间ti时刻刀具中心坐标满足Pi=P0+(P1-P0)/|(P1-P0)|×f×ti,其中P0和P1为相邻两刀位点坐标,ti为从P0开始到Pi点所经历时间。综上所述,刀具平移矩阵为:(2.2) Calculation of tool translation matrix. The tool position file coordinates in the machining process are generated relative to the workpiece coordinate system. In order to facilitate the calculation, the tool coordinate system and the workpiece coordinate system are coincident in the initial modeling state, and the tool moves between any two tool positions The process can be calculated according to the generated tool location file. In three-axis machining, the moving process of the tool between two tool positions is linear interpolation, that is, the coordinates of the tool center at time t i between the two tool positions satisfy P i =P 0 +(P 1 -P 0 )/|( P 1 -P 0 )|×f×t i , where P 0 and P 1 are the coordinates of two adjacent knife points, and t i is the time elapsed from P 0 to P i . In summary, the tool translation matrix is:
其中e1 T、e2 T、e3 T分别为矢量P0P1的单位矢量的X、Y、Z方向的分量。Where e 1 T , e 2 T , and e 3 T are components in the X, Y, and Z directions of the unit vector of the vector P 0 P 1, respectively.
(2.3)计算两刀位点之间刀刃空间位置。工件三维表面形貌是刀刃切削工件形成的,最终工件表面的形状取决于刀刃在切削过程中形成的包络面,而刀具在加工过程中进给速度、刀具转速和刀具路径决定了刀刃上任意一点的运动为空间螺旋运动。为了表示刀刃在加工过程中不同时刻的位置,定义Δω作为刀刃位置的离散参数对仿真时间进行分割,以保证在两刀位点之间有足够满足仿真精度的刀刃扫掠点点云生成。(2.3) Calculate the spatial position of the blade between the two knife positions. The three-dimensional surface morphology of the workpiece is formed by the cutting edge of the workpiece. The final shape of the workpiece surface depends on the envelope surface formed by the cutting edge during the cutting process, and the feed rate, tool speed and tool path of the tool determine any shape on the cutting edge. The movement of one point is space spiral movement. In order to represent the position of the blade at different moments in the machining process, Δω is defined as the discrete parameter of the blade position to divide the simulation time, so as to ensure that there are enough blade sweep point cloud generation between the two blade positions to meet the simulation accuracy.
其中ti为从P0开始到Pi点所经历时间,根据Δω对时间t=(p1-p0)/f进行离散,得到计算两刀位点之间不同离散时刻下的刀刃离散点位置Bi。Among them, t i is the time elapsed from P 0 to P i point, according to Δω, the time t=(p 1 -p 0 )/f is discretized, and the discretized point of the blade at different discretized moments between the two knife positions can be obtained position B i .
根据步骤(1)中提取出到仿真刀位点,对其中相邻两刀位点依次执行步骤(2)直到得到仿真区域所有刀位点的刀刃扫掠点云模型,将这些所有的刀刃离散点Bi叠加在一起就构成了刀刃扫掠点云模型。According to the simulated knife position extracted in step (1), step (2) is performed in turn for two adjacent knife positions until the blade sweep point cloud model of all knife positions in the simulation area is obtained, and all these knife edges are discretized The point B i is superimposed to form the blade sweep point cloud model.
(3)根据建立的刀刃扫掠点云模型提取三维表面形貌点。虽然工件三维表面形貌的形状有许多因素决定,例如机床振动、刀具路径、刀刃模型、切削参数和初始切入相位角等,并且许多影响因素具有不确定性,例如机床振动和刀具变形,但这些因素最终都会表现在刀刃切削工件材料后的残留形状,即最终的表面形貌由工件最终的残留形状决定。因而,表面形貌的提取问题变成了刀刃扫掠点靠近加工曲面一侧的包络面提取,具体步骤如下:(3) Extract the three-dimensional surface topography points according to the established blade sweep point cloud model. Although the shape of the three-dimensional surface topography of the workpiece is determined by many factors, such as machine vibration, tool path, cutting edge model, cutting parameters and initial cut-in phase angle, etc., and many influencing factors have uncertainties, such as machine vibration and tool deformation, but these Factors will eventually be reflected in the residual shape of the workpiece after the cutting edge cuts the material, that is, the final surface morphology is determined by the final residual shape of the workpiece. Therefore, the problem of surface topography extraction becomes the envelope surface extraction of the edge sweep point close to the processing surface. The specific steps are as follows:
(3.1)根据步骤(2)得到了仿真区域刀刃扫掠点云模型,提取出影响最终表面形貌的点云集合Bextract。(3.1) According to the step (2), the blade-swept point cloud model of the simulation area is obtained, and the point cloud set B extract that affects the final surface morphology is extracted.
如图3所示,对仿真区域Ssim(u,v)进行偏置,偏置量为切削深度,形成偏置面Ssimoff(u,v),根据Ssim(u,v)和Ssimoff(u,v)提取出位于这两个平面之间的刀刃扫掠点云Bextract。As shown in Figure 3, the simulation area S sim (u, v) is offset, and the offset amount is the cutting depth to form an offset surface S simoff (u, v), according to S sim (u, v) and S simoff (u, v) extracts the blade-swept point cloud B extract located between these two planes.
(3.2)对仿真区域Ssim(u,v)按照三维表面形貌采样点数pu、pv划分网格,记录每个网格节点信息Pknot(m,n)和仿真区域边界信息。每一个网格区域可以对应一部分刀刃扫掠点,如图4所示,以Pknot(m,n)处曲面法失为边界方向、网格边界为边界建立包容盒,位于包容盒内的点云作为该区域对应的刀刃扫掠点云Bbox(m,n)。(3.2) The simulation area S sim (u, v) is divided into grids according to the three-dimensional surface topography sampling points p u and p v , and the information of each grid node P kno t(m, n) and the boundary information of the simulation area are recorded. Each grid area can correspond to a part of the blade sweep points, as shown in Figure 4, the surface method at P knot (m, n) is used as the boundary direction, and the grid boundary is used as the boundary to establish a containment box, and the points located in the containment box The cloud is used as the blade-swept point cloud B box (m, n) corresponding to this area.
(3.3)提取每个网格包容盒内对应刀刃扫掠点云Bbox(m,n)中离加工曲面最近的点,作为该区域内三维表面形貌的表征点Ps(m,n)。(3.3) Extract the nearest point to the processing surface in the corresponding blade sweep point cloud B box (m, n) in each grid containment box, as the representative point P s (m, n) of the three-dimensional surface topography in the area .
如图4所示,在任意一个网格包容盒内的刀刃点云中,选取一个点作为该区域的表征点。为了使该区域内选取的表征点对该网格内区域具有代表性,需要选取其中到加工曲面距离最近的点作为该区域的三维表面形貌表征点。As shown in Figure 4, in any point cloud of the blade in the enclosing box of a grid, a point is selected as the representative point of the area. In order to make the selected characterization points in this area representative of the area in the grid, it is necessary to select the point with the closest distance to the processed surface as the 3D surface topography characterization point of this area.
以Pknot(m,n)结点处加工曲面法失、u向切失、v向切失和点Pknot(m,n)建立原点在Pknot(m,n)的动态坐标系{R(m,n)}。{R(m,n)}下的一组基可表示为: The dynamic coordinate system { R (m,n)}. A set of basis under {R(m,n)} can be expressed as:
为找到离加工曲面距离最近的点,对Pknot(m,n)网格对应刀刃扫掠点云Bbox(m,n)进行位姿变换,获得Bbox(m,n)变换后在坐标系{R(m,n)}下的坐标R(m,n)Bbox(m,n)In order to find the point closest to the processing surface, the P knot (m, n) grid corresponds to the blade sweep point cloud B box (m, n) for pose transformation, and after the B box (m, n) transformation is obtained, the coordinates Coordinates R(m, n) B box ( m, n) under the system {R(m , n)}
其中,T为从工件坐标系{WCS}到动态坐标系{R(m,n)}的齐次变换矩阵;提取{R(m,n)}坐标系下R(m,n)Bbox(m,n)中Z值最小的点作为该区域网格中的三维表面形貌的表征点,并进行逆变换到工件坐标系下,得到对应形貌特征点Ps(m,n)。(3.4)按照(3.3)步骤依次遍历所有网格点直到每个网格区域中的三维表面形貌表征点在in, T is the homogeneous transformation matrix from the workpiece coordinate system {WCS} to the dynamic coordinate system {R(m, n)}; extract R(m, n) B box (m, The point with the smallest Z value in n) is used as the representative point of the three-dimensional surface topography in the grid of the area, and is inversely transformed into the workpiece coordinate system to obtain the corresponding topography feature point P s (m, n). (3.4) According to the step (3.3), traverse all the grid points in turn until the three-dimensional surface topography characterization points in each grid area are at
工件坐标系{WCS}下坐标全部算出,即仿真区域中全部三维表面形貌特征点Ps。All coordinates in the workpiece coordinate system {WCS} are calculated, that is, all three-dimensional surface topography feature points P s in the simulation area.
Claims (5)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110359576 CN102490081B (en) | 2011-11-14 | 2011-11-14 | Workpiece three-dimensional surface topography simulating method based on ball head milling |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110359576 CN102490081B (en) | 2011-11-14 | 2011-11-14 | Workpiece three-dimensional surface topography simulating method based on ball head milling |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102490081A true CN102490081A (en) | 2012-06-13 |
CN102490081B CN102490081B (en) | 2013-07-24 |
Family
ID=46181964
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110359576 Active CN102490081B (en) | 2011-11-14 | 2011-11-14 | Workpiece three-dimensional surface topography simulating method based on ball head milling |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102490081B (en) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103042434A (en) * | 2012-12-08 | 2013-04-17 | 华中科技大学 | Control method of surface topography and texture in ultra-precision milling |
CN103949701A (en) * | 2014-04-08 | 2014-07-30 | 西安理工大学 | Special-shaped part milling processing method based on rotating vectors |
CN105740497A (en) * | 2014-12-31 | 2016-07-06 | 达索系统公司 | Simulating the machining of a workpiece |
CN106599389A (en) * | 2016-11-11 | 2017-04-26 | 大连理工大学 | Cutting edge contact distinguishing method for flat-end end mill strip maximization machining |
CN106873522A (en) * | 2017-02-13 | 2017-06-20 | 武汉理工大学 | A kind of numerical control turning cutter path planing method of non-axis symmetry sweeping surface |
CN107066654A (en) * | 2016-12-15 | 2017-08-18 | 广东工业大学 | Towards the TRAJECTORY CONTROL dot picking method and apparatus in the enveloped box face of point cloud model |
CN108127481A (en) * | 2017-12-15 | 2018-06-08 | 北京理工大学 | A kind of Forecasting Methodology of the workpiece surface appearance based on Flank machining |
CN108515217A (en) * | 2018-04-09 | 2018-09-11 | 吉林大学 | A kind of ball-end milling free form surface surface topography emulation mode |
CN109716248A (en) * | 2016-07-05 | 2019-05-03 | 重力.瑞士股份公司 | Method for being machined embossment |
CN109767486A (en) * | 2018-12-27 | 2019-05-17 | 广州纬纶信息科技有限公司 | A kind of special-shaped workpiece cutting modeling method, electronic equipment, storage medium and system |
CN110385601A (en) * | 2019-07-10 | 2019-10-29 | 湖南大学 | A kind of generation method of the ball cutter knife rail in car panel die bedding-in type face |
CN110497088A (en) * | 2019-08-31 | 2019-11-26 | 大连理工大学 | Laser processing error control method for flexible conformal antenna based on surface mapping |
CN111400830A (en) * | 2019-01-03 | 2020-07-10 | 上海狮迈科技有限公司 | Machining and calibrating method and device for three-dimensional blank workpiece |
CN112417537A (en) * | 2020-12-16 | 2021-02-26 | 华中科技大学 | Surface morphology simulation method for visualizing geometric errors based on turning |
CN107577882B (en) * | 2017-09-12 | 2021-04-13 | 电子科技大学 | A Surface Topography Modeling and Forming Simulation Method for Side Milled Ruled Surfaces |
WO2021174518A1 (en) * | 2020-03-06 | 2021-09-10 | 大连理工大学 | Flutter-free milling surface topography simulation method |
USRE48940E1 (en) | 2011-11-03 | 2022-02-22 | Dassault Systèmes | Simulation of the machining of a workpiece |
CN114812486A (en) * | 2022-05-13 | 2022-07-29 | 武汉理工大学 | Method and device for acquiring surface roughness of machined workpiece and electronic equipment |
CN116977652A (en) * | 2023-09-22 | 2023-10-31 | 之江实验室 | Workpiece surface morphology generation method and device based on multi-mode image generation |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS61274845A (en) * | 1985-05-29 | 1986-12-05 | Okuma Mach Works Ltd | Method of simulating track of tool in nc device |
CN1752876A (en) * | 2005-11-01 | 2006-03-29 | 西北工业大学 | Space Scanning Volume Construction Method Based on Compressed Voxel Model |
CN101881599A (en) * | 2010-07-12 | 2010-11-10 | 华中科技大学 | Method and device for measuring three-dimensional shape of nanostructures |
-
2011
- 2011-11-14 CN CN 201110359576 patent/CN102490081B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS61274845A (en) * | 1985-05-29 | 1986-12-05 | Okuma Mach Works Ltd | Method of simulating track of tool in nc device |
CN1752876A (en) * | 2005-11-01 | 2006-03-29 | 西北工业大学 | Space Scanning Volume Construction Method Based on Compressed Voxel Model |
CN101881599A (en) * | 2010-07-12 | 2010-11-10 | 华中科技大学 | Method and device for measuring three-dimensional shape of nanostructures |
Non-Patent Citations (2)
Title |
---|
K.SAWADA 等: "超精密铣削的三维微加工工艺", 《国防科技参考》 * |
彭芳瑜: "光学零件超精密铣削加工表面形貌预测与仿真", 《中国优秀硕士论文电子期刊网》 * |
Cited By (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
USRE48940E1 (en) | 2011-11-03 | 2022-02-22 | Dassault Systèmes | Simulation of the machining of a workpiece |
CN103042434B (en) * | 2012-12-08 | 2014-12-10 | 华中科技大学 | Ultraprecise milling surface topography texture control method |
CN103042434A (en) * | 2012-12-08 | 2013-04-17 | 华中科技大学 | Control method of surface topography and texture in ultra-precision milling |
CN103949701A (en) * | 2014-04-08 | 2014-07-30 | 西安理工大学 | Special-shaped part milling processing method based on rotating vectors |
CN103949701B (en) * | 2014-04-08 | 2016-04-20 | 西安理工大学 | Based on the irregular part milling method of rotating vector |
CN105740497A (en) * | 2014-12-31 | 2016-07-06 | 达索系统公司 | Simulating the machining of a workpiece |
CN105740497B (en) * | 2014-12-31 | 2021-07-09 | 达索系统公司 | Simulating the machining of a workpiece |
CN109716248A (en) * | 2016-07-05 | 2019-05-03 | 重力.瑞士股份公司 | Method for being machined embossment |
CN106599389A (en) * | 2016-11-11 | 2017-04-26 | 大连理工大学 | Cutting edge contact distinguishing method for flat-end end mill strip maximization machining |
CN106599389B (en) * | 2016-11-11 | 2019-08-13 | 大连理工大学 | A kind of blade contact method of discrimination of square end mill wide-line processing |
CN107066654A (en) * | 2016-12-15 | 2017-08-18 | 广东工业大学 | Towards the TRAJECTORY CONTROL dot picking method and apparatus in the enveloped box face of point cloud model |
CN107066654B (en) * | 2016-12-15 | 2020-03-10 | 广东工业大学 | Box surface containing track control point picking method and device facing point cloud model |
CN106873522A (en) * | 2017-02-13 | 2017-06-20 | 武汉理工大学 | A kind of numerical control turning cutter path planing method of non-axis symmetry sweeping surface |
CN107577882B (en) * | 2017-09-12 | 2021-04-13 | 电子科技大学 | A Surface Topography Modeling and Forming Simulation Method for Side Milled Ruled Surfaces |
CN108127481A (en) * | 2017-12-15 | 2018-06-08 | 北京理工大学 | A kind of Forecasting Methodology of the workpiece surface appearance based on Flank machining |
CN108515217A (en) * | 2018-04-09 | 2018-09-11 | 吉林大学 | A kind of ball-end milling free form surface surface topography emulation mode |
CN108515217B (en) * | 2018-04-09 | 2019-05-31 | 吉林大学 | A kind of ball-end milling free form surface surface topography emulation mode |
CN109767486A (en) * | 2018-12-27 | 2019-05-17 | 广州纬纶信息科技有限公司 | A kind of special-shaped workpiece cutting modeling method, electronic equipment, storage medium and system |
CN109767486B (en) * | 2018-12-27 | 2022-01-04 | 广州纬纶信息科技有限公司 | Special-shaped workpiece cutting modeling method, electronic equipment, storage medium and system |
CN111400830B (en) * | 2019-01-03 | 2023-07-11 | 西北工业大学太仓长三角研究院 | Machining calibration method and device for three-dimensional blank workpiece |
CN111400830A (en) * | 2019-01-03 | 2020-07-10 | 上海狮迈科技有限公司 | Machining and calibrating method and device for three-dimensional blank workpiece |
CN110385601A (en) * | 2019-07-10 | 2019-10-29 | 湖南大学 | A kind of generation method of the ball cutter knife rail in car panel die bedding-in type face |
CN110497088A (en) * | 2019-08-31 | 2019-11-26 | 大连理工大学 | Laser processing error control method for flexible conformal antenna based on surface mapping |
WO2021174518A1 (en) * | 2020-03-06 | 2021-09-10 | 大连理工大学 | Flutter-free milling surface topography simulation method |
CN112417537A (en) * | 2020-12-16 | 2021-02-26 | 华中科技大学 | Surface morphology simulation method for visualizing geometric errors based on turning |
CN112417537B (en) * | 2020-12-16 | 2024-02-02 | 华中科技大学 | Surface morphology simulation method for visualizing geometric errors based on turning |
CN114812486A (en) * | 2022-05-13 | 2022-07-29 | 武汉理工大学 | Method and device for acquiring surface roughness of machined workpiece and electronic equipment |
CN116977652A (en) * | 2023-09-22 | 2023-10-31 | 之江实验室 | Workpiece surface morphology generation method and device based on multi-mode image generation |
CN116977652B (en) * | 2023-09-22 | 2023-12-22 | 之江实验室 | Workpiece surface morphology generation method and device based on multi-mode image generation |
Also Published As
Publication number | Publication date |
---|---|
CN102490081B (en) | 2013-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102490081B (en) | Workpiece three-dimensional surface topography simulating method based on ball head milling | |
JP7015041B2 (en) | Cutting simulation method and equipment | |
CN104759942B (en) | A kind of milling deformation on-line measurement of thin-walled parts and compensation processing method | |
Cho et al. | Inspection planning strategy for the on-machine measurement process based on CAD/CAM/CAI integration | |
CN102528555A (en) | Geometry and mechanics integrated optimization information processing method of interference-free tool path on complex curved surface | |
CN107577882B (en) | A Surface Topography Modeling and Forming Simulation Method for Side Milled Ruled Surfaces | |
Qin et al. | CWE identification and cutting force prediction in ball-end milling process | |
CN102866672A (en) | Online detecting method of numerical control machining middle state of plane structural member | |
US20230204344A1 (en) | Method for inspecting a workpiece in the process of being manufactured | |
CN102621928A (en) | Method for generating pose angle feasible regions of five-axis tools quickly | |
CN111459097A (en) | A Calculation Method for Surface Machining Contact Area of Ball End Milling Cutter | |
CN113618488A (en) | A method of centering the center of rotation of the B-axis and the center of the circular arc of the cutting edge | |
CN104462775A (en) | Five-axis machining center milling force calculating method | |
Zeroudi et al. | Prediction of machined surface geometry based on analytical modelling of ball-end milling | |
Bo et al. | Geometry simulation and evaluation of the surface topography in five-axis ball-end milling | |
Baburaj et al. | Study of micro ball end mill geometry and measurement of cutting edge radius | |
Cui et al. | Sub-regional polishing and machining trajectory selection of complex surface based on K9 optical glass | |
CN114417526A (en) | An accurate prediction method for different grinding surface topography | |
JP7446055B2 (en) | Method and machinery for manufacturing cutting tools | |
Li et al. | A novel path generation method of onsite 5-axis surface inspection using the dual-cubic NURBS representation | |
Sun et al. | A review on theories/methods to obtain surface topography and analysis of corresponding affecting factors in the milling process | |
Xing et al. | On-machine measurement method and geometrical error analysis in a multi-step processing system of an ultra-precision complex spherical surface | |
JP2019107763A5 (en) | ||
Chen et al. | Modeling and simulation of the surface topography in ball-end milling based on biharmonic spline interpolation | |
CN106681278A (en) | Cutter and workpiece instantaneous contact contour extraction method during five-axis side-milling machining |
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 |