CN108732990A - 一种三维圆及椭圆的等角度逼近方法 - Google Patents

一种三维圆及椭圆的等角度逼近方法 Download PDF

Info

Publication number
CN108732990A
CN108732990A CN201810374025.6A CN201810374025A CN108732990A CN 108732990 A CN108732990 A CN 108732990A CN 201810374025 A CN201810374025 A CN 201810374025A CN 108732990 A CN108732990 A CN 108732990A
Authority
CN
China
Prior art keywords
dimensional
circle
elliptical
point
dimensional circle
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
CN201810374025.6A
Other languages
English (en)
Other versions
CN108732990B (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.)
Huaqiao University
Original Assignee
Huaqiao 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 Huaqiao University filed Critical Huaqiao University
Priority to CN201810374025.6A priority Critical patent/CN108732990B/zh
Publication of CN108732990A publication Critical patent/CN108732990A/zh
Application granted granted Critical
Publication of CN108732990B publication Critical patent/CN108732990B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/19Numerical 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 positioning or contouring control systems, e.g. to control position from one programmed point to another or to control movement along a programmed continuous path
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/32Operator till task planning
    • G05B2219/32317Smallest ratio for imminent processing time divided by total processing time

Abstract

本发明涉及一种三维圆的等角度逼近方法和三维椭圆的等角度逼近方法,具体为根据从DXF文件中获取的三维圆和椭圆的图元信息,用向量法推导出三维圆和椭圆的参数方程;根据三维圆和椭圆的参数方程,每隔一定的角度增量△θ取值,计算三维圆和椭圆上的点坐标;将点坐标导入AutoCAD中与原图进行对比,对比结果验证了本发明逼近方法的正确性。本发明的等角度逼近方法是一种根据参数方程(以角度为参数)对三维圆和椭圆进行数据点密化的等角度逼近方法;根据DXF文件中的数据特点,通过等角度分割,方便、高效地对三维圆和椭圆进行数据点密化。

Description

一种三维圆及椭圆的等角度逼近方法
技术领域
本发明涉及数控自动化加工领域,尤其涉及一种三维圆及椭圆的等角度逼近方法,具体涉及一种根据参数方程(以角度为参数)对三维圆和椭圆进行数据点密化的等角度逼近方法。
背景技术
目前对曲线进行直线段逼近的算法有等间距法、等弦长法、等误差法等,上述逼近算法若推广到三维空间,其计算较为繁琐。数控自动化加工中往往需从DXF图形文件中获取图形信息,而DXF图形文件中存储了三维圆、椭圆的中心点坐标O′(xC,yC,zC)和起始角度θS、终止角度θE等信息。因此,根据这些数据的特点,方便、高效地将三维圆和椭圆细分,并进行数据点密化,具有显著的意义。
发明内容
本发明的目的在于克服现有技术的不足,提供一种三维圆及椭圆的等角度逼近方法,具体为根据参数方程(以角度为参数)对三维圆及椭圆进行数据点密化的等角度逼近方法;根据DXF文件中的数据特点,通过等角度分割,方便、高效地对三维圆和椭圆进行数据点密化。
为了实现上述目的,本发明的技术方案是:
一方面,一种三维圆的等角度逼近方法,根据角度参数方程对三维圆进行数据点密化,包括:
从DXF文件中获取三维圆的图元信息;
用向量法推导出三维圆的角度参数方程;
根据所述角度参数方程计算三维圆上的点坐标。
优选的,所述从DXF文件中获取三维圆的图元信息,包括:
从DXF文件中获取三维圆的中心点坐标O′(xC,yC,zC)、半径R、起始角度θS、终止角度θE及所在平面的法线
优选的,所述用向量法推导出三维圆的角度参数方程,包括:
根据所述三维圆的图元信息,将三维圆上任意一点Pi用向量表示为:
写成坐标形式的参数方程,如下:
其中,i为非负整数,i∈0,1,2…;O为三维坐标的坐标原点;θi为角度参数;为既垂直于又互相垂直的单位向量。
优选的,所述根据所述角度参数方程计算三维圆上的点坐标,包括:
根据所述参数方程,令θi从θS到θE每隔△θ取值,即θi=θS+i·△θ,计算圆上的点坐标Pi(xi,yi,zi)。
另一方面,本发明一种三维椭圆的等角度逼近方法,根据角度参数方程对三维椭圆进行数据点密化,包括:
从DXF文件中获取三维椭圆的图元信息;
用向量法推导出三维椭圆的角度参数方程;
根据所述角度参数方程计算三维椭圆上的点坐标。
优选的,所述从DXF文件中获取三维椭圆的图元信息,包括:
从DXF文件中获取三维椭圆的中心点坐标O′(xC,yC,zC)、长半轴RL、短半轴RS、起始角度θS、终止角度θE及所在平面的法线
优选的,所述用向量法推导出三维椭圆的角度参数方程,包括:
根据所述三维椭圆的图元信息,将三维椭圆上任意一点Pi用向量表示为:
写成坐标形式的参数方程,如下:
其中,i为非负整数,i∈0,1,2…;O为三维坐标的坐标原点,θi为角度参数,是既垂直于又互相垂直的单位向量。
优选的,所述根据所述角度参数方程计算三维椭圆上的点坐标,包括:
根据所述参数方程,令θi从θS到θE每隔△θ取值,即θi=θS+i·△θ,计算椭圆上的点坐标Pi(xi,yi,zi)。
本发明的有益效果如下:
(1)根据参数方程(以角度为参数),通过等角度分割,方便、高效地对三维圆和椭圆进行数据点密化。
(2)直接根据圆、椭圆三维空间任意位置的参数方程进行计算,无需坐标变换,本发明逼近方法能在一定程度上扩大数控自动化加工范围。
以下结合附图及实施例对本发明作进一步详细说明,但本发明的一种三维圆及椭圆的等角度逼近方法不局限于实施例。
附图说明
图1是本发明的三维圆等角度逼近方法流程图;
图2是本发明的三维圆空间示意图;
图3是本发明的三维圆所在平面示意图;
图4是本发明的三维圆对比示意图;
图5是本发明的三维椭圆等角度逼近方法流程图;
图6是本发明的三维椭圆空间示意图;
图7是本发明的三维椭圆所在平面示意图;
图8是本发明的三维椭圆对比示意图。
具体实施方式
下面结合附图,对本发明的技术方案进行具体说明。
实施例1
如图1所示,一方面,本发明一种三维圆的等角度逼近方法,是一种根据参数方程(以角度为参数)对三维圆进行数据点密化的等角度逼近方法,具体包括如下步骤:
步骤101,从DXF文件中获取三维圆的图元信息
具体的,从DXF文件中获取三维圆的中心点坐标O′(xC,yC,zC)、半径R、起始角度θS、终止角度θE(闭合圆θS=0°、θE=360°),及其所在平面的法线
步骤102,用向量法推导出三维圆的参数方程
如图2至3所示,根据所述三维圆的图元信息,R为三维圆的半径,O′(xC,yC,zC)为三维圆的中心点,Pi(xi,yi,zi)为三维圆上的任意一点,则
1)Pi点用向量表示为
2)向量用圆所在平面上相互垂直的两个向量表示,即
3)
其中,方向的单位向量,方向的单位向量。
综上,三维圆上任意一点Pi用向量表示为:
写成坐标形式的参数方程,如下:
式(1)-(2)中,i为非负整数,i∈0,1,2…;O为三维坐标的坐标原点;θi为角度参数;是既垂直于又互相垂直的单位向量,其作用是将三维圆限定在指定的平面内。具体的,可用下述方法求出向量。
的方法:用Z坐标轴的方向向量叉乘如果叉乘结果不为零,那么叉乘结果必然垂直于将此叉乘结果单位化并作为向量如果叉乘结果为零,再用剩下的Y坐标轴的方向向量与X坐标轴的方向向量中任意一个叉乘单位化叉乘结果作为向量
的方法:叉乘上一步得到的向量叉乘结果必然垂直于单位化此叉乘结果作为向量
步骤103,根据所述参数方程计算三维圆上的点坐标
根据所述参数方程,令θi从θS到θE每隔△θ取值,即θi=θS+i·△θ(i=0,1,2…),计算圆、椭圆上的点坐标Pi(xi,yi,zi)。
下面结合图4对本发明一种三维圆的等角度逼近方法做进一步说明:
通过AutoCAD绘制:三维圆,中心点O′(0,0,50)、半径R=100mm;令△θ=3根据步骤103计算出三维圆上的点坐标,把点坐标导入AutoCAD中与原图进行对比,对比结果如图4所示,从对比结果看出,通过式(2)计算出的点完全在原曲线上,验证了步骤102推导出的三维圆的参数方程的正确性及本发明逼近方法的正确性。
实施例2
如图5所示,另一方面,本发明一种三维椭圆的等角度逼近方法,是一种根据参数方程(以角度为参数)对三维椭圆进行数据点密化的等角度逼近方法,具体包括如下步骤:
步骤201,从DXF文件中获取三维椭圆的图元信息
具体的,从DXF文件中获取三维圆的中心点坐标O′(xC,yC,zC)、长半轴RL、短半轴RS、起始角度θS、终止角度θE(闭合椭圆θS=0°、θE=360°),及其所在平面的法线
步骤202,用向量法推导出三维椭圆的参数方程
如图6至7所示,根据所述三维椭圆的图元信息,RL、RS分别为三维椭圆的长半轴、短半轴,图中位于椭圆外的大圆以RL为半径,位于椭圆内的小圆以RS为半径,O′(xC,yC,zC)为三维椭圆的中心点,Pi(xi,yi,zi)为三维椭圆上的任意一点,则
1)Pi点用向量表示为
2)向量用椭圆所在平面上相互垂直的两个向量表示,即
3)
其中,方向的单位向量,方向的单位向量。
综上,三维椭圆上任意一点Pi用向量表示为:
写成坐标形式的参数方程,如下:
式(3)-(4)中,i为非负整数,i∈0,1,2…;O为三维坐标的坐标原点,θi为角度参数,是既垂直于又互相垂直的单位向量,其作用是将三维椭圆限定在指定的平面内,可用下述方法求出向量:
的方法:用Z坐标轴的方向向量叉乘如果叉乘结果不为零,那么叉乘结果必然垂直于将此叉乘结果单位化并作为向量如果叉乘结果为零,再用剩下的Y坐标轴的方向向量与X坐标轴的方向向量中任意一个叉乘单位化叉乘结果作为向量
的方法:叉乘上一步得到的向量叉乘结果必然垂直于单位化此叉乘结果作为向量
步骤203,根据所述参数方程计算三维椭圆上的点坐标
根据所述参数方程,令θi从θS到θE每隔△θ取值,即θi=θS+i·△θ(i=0,1,2…),计算椭圆上的点坐标Pi(xi,yi,zi)。
下面结合图8对本发明一种三维椭圆的等角度逼近方法做进一步说明:
通过AutoCAD绘制:三维椭圆,中心点O′(0,0,50)、长半轴RL=100mm、短半轴RS=85mm;令△θ=3根据步骤203计算出三维椭圆上的点坐标,把点坐标导入AutoCAD中与原图进行对比,对比结果如图8所示,从对比结果看出,通过式(4)计算出的点完全在原曲线上,验证了步骤202推导出的三维椭圆参数方程的正确性及本发明逼近方法的正确性。
以上仅为本发明实例中一个较佳的实施方案。但是,本发明并不限于上述实施方案,凡按本发明所做的任何均等变化和修饰,所产生的功能作用未超出本方案的范围时,均属于本发明的保护范围。

Claims (8)

1.一种三维圆的等角度逼近方法,其特征在于,根据角度参数方程对三维圆进行数据点密化,包括:
从DXF文件中获取三维圆的图元信息;
用向量法推导出三维圆的角度参数方程;
根据所述角度参数方程计算三维圆上的点坐标。
2.根据权利要求1所述的三维圆的等角度逼近方法,其特征在于,所述从DXF文件中获取三维圆的图元信息,包括:
从DXF文件中获取三维圆的中心点坐标O′(xC,yC,zC)、半径R、起始角度θS、终止角度θE及所在平面的法线
3.根据权利要求2所述的三维圆的等角度逼近方法,其特征在于,所述用向量法推导出三维圆的角度参数方程,包括:
根据所述三维圆的图元信息,将三维圆上任意一点Pi用向量表示为:
写成坐标形式的参数方程,如下:
其中,i为非负整数,i∈0,1,2…;O为三维坐标的坐标原点;θi为角度参数;为既垂直于又互相垂直的单位向量。
4.根据权利要求3所述的三维圆的等角度逼近方法,其特征在于,所述根据所述角度参数方程计算三维圆上的点坐标,包括:
根据所述参数方程,令θi从θS到θE每隔△θ取值,即θi=θS+i·△θ,计算圆上的点坐标Pi(xi,yi,zi)。
5.一种三维椭圆的等角度逼近方法,其特征在于,根据角度参数方程对三维椭圆进行数据点密化,包括:
从DXF文件中获取三维椭圆的图元信息;
用向量法推导出三维椭圆的角度参数方程;
根据所述角度参数方程计算三维椭圆上的点坐标。
6.根据权利要求5所述的三维椭圆的等角度逼近方法,其特征在于,所述从DXF文件中获取三维椭圆的图元信息,包括:
从DXF文件中获取三维椭圆的中心点坐标O′(xC,yC,zC)、长半轴RL、短半轴RS、起始角度θS、终止角度θE及所在平面的法线
7.根据权利要求6所述的三维圆的等角度逼近方法,其特征在于,所述用向量法推导出三维椭圆的角度参数方程,包括:
根据所述三维椭圆的图元信息,将三维椭圆上任意一点Pi用向量表示为:
写成坐标形式的参数方程,如下:
其中,i为非负整数,i∈0,1,2…;O为三维坐标的坐标原点,θi为角度参数,是既垂直于又互相垂直的单位向量。
8.根据权利要求7所述的三维椭圆的等角度逼近方法,其特征在于,所述根据所述角度参数方程计算三维椭圆上的点坐标,包括:
根据所述参数方程,令θi从θS到θE每隔△θ取值,即θi=θS+i·△θ,计算椭圆上的点坐标Pi(xi,yi,zi)。
CN201810374025.6A 2018-04-24 2018-04-24 一种三维圆及椭圆的等角度逼近方法 Active CN108732990B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810374025.6A CN108732990B (zh) 2018-04-24 2018-04-24 一种三维圆及椭圆的等角度逼近方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810374025.6A CN108732990B (zh) 2018-04-24 2018-04-24 一种三维圆及椭圆的等角度逼近方法

Publications (2)

Publication Number Publication Date
CN108732990A true CN108732990A (zh) 2018-11-02
CN108732990B CN108732990B (zh) 2020-10-09

Family

ID=63939247

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810374025.6A Active CN108732990B (zh) 2018-04-24 2018-04-24 一种三维圆及椭圆的等角度逼近方法

Country Status (1)

Country Link
CN (1) CN108732990B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03245207A (ja) * 1990-02-23 1991-10-31 Okuma Mach Works Ltd 非真円形状データ作成装置
JPH06282659A (ja) * 1993-03-25 1994-10-07 Daikin Ind Ltd 半球面上への多角形投影方法、多角形投影面積算出方法、フォームファクタ算出方法およびこれらの装置
CN1540469A (zh) * 2003-10-28 2004-10-27 华北电力大学 一种用于数控机床的三坐标圆弧插补方法及装置
CN103761376A (zh) * 2014-01-10 2014-04-30 沈阳航空航天大学 基于二维dxf文件格式的零件三维真实感显示方法
CN104898556A (zh) * 2015-05-12 2015-09-09 重庆大学 由成形刀具确定被加工螺旋曲面廓形的点矢量包络法
CN111338286A (zh) * 2020-03-06 2020-06-26 吉林大学 基于加工精度控制的伞叶面超精密车削刀具轨迹生成方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03245207A (ja) * 1990-02-23 1991-10-31 Okuma Mach Works Ltd 非真円形状データ作成装置
JPH06282659A (ja) * 1993-03-25 1994-10-07 Daikin Ind Ltd 半球面上への多角形投影方法、多角形投影面積算出方法、フォームファクタ算出方法およびこれらの装置
CN1540469A (zh) * 2003-10-28 2004-10-27 华北电力大学 一种用于数控机床的三坐标圆弧插补方法及装置
CN103761376A (zh) * 2014-01-10 2014-04-30 沈阳航空航天大学 基于二维dxf文件格式的零件三维真实感显示方法
CN104898556A (zh) * 2015-05-12 2015-09-09 重庆大学 由成形刀具确定被加工螺旋曲面廓形的点矢量包络法
CN111338286A (zh) * 2020-03-06 2020-06-26 吉林大学 基于加工精度控制的伞叶面超精密车削刀具轨迹生成方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张志强: "椭圆数控编程参数的选择及R参数编程", 《机械工程与自动化》 *
柳青: "宏程序在椭圆面铣削加工中的应用", 《精密制造与自动化》 *

Also Published As

Publication number Publication date
CN108732990B (zh) 2020-10-09

Similar Documents

Publication Publication Date Title
US10627241B2 (en) Map-centric map matching method and apparatus
Takayama et al. Sketch-based generation and editing of quad meshes.
Rukavishnikov et al. Weighted finite-element method for elasticity problems with singularity
CN108022262A (zh) 一种基于点的邻域重心向量特征的点云配准方法
CN104134236B (zh) 一种三维平面实体的布尔运算方法
CN108830787A (zh) 图像变形的方法、装置及电子设备
CN105716618A (zh) 一种用于uuv航路规划的几何环境模型膨化处理方法
Alton et al. Fast marching methods for stationary Hamilton–Jacobi equations with axis-aligned anisotropy
Shen et al. Implicitizing rational tensor product surfaces using the resultant of three moving planes
Treschev On a conjugacy problem in billiard dynamics
CN103489221B (zh) 四边形网格共形参数化方法
CN107562779B (zh) 一种二维矢量多边形自相交的空间拓扑处理方法
CN108732990A (zh) 一种三维圆及椭圆的等角度逼近方法
McGann et al. Hamilton–Jacobi theory for continuation of magnetic field across a toroidal surface supporting a plasma pressure discontinuity
Kublik et al. An extrapolative approach to integration over hypersurfaces in the level set framework
WO2024055466A1 (zh) 一种基于二维激光雷达的环境直角点坐标角度的提取方法
Nguyen et al. Isogeometric segmentation: Construction of auxiliary curves
CN115481268A (zh) 一种平行线墙体数据识别方法、装置、设备及存储介质
CN104933205A (zh) 土地利用空间数据处理中基于地理本体的属性匹配方法
Kurz et al. Methods for Deterministic Approximation of Circular Densities
Zhou et al. A NOVEL ADAPTIVE FINITE VOLUME METHOD FOR ELLIPTIC EQUATIONS.
Frego et al. Point-clothoid distance and projection computation
Jian et al. An efficient mesh simplification method in 3d graphic model rendering
CN116342849B (zh) 在三维网格上生成牙模倒凹区域的方法
Zheng et al. An algorithm for computing geodesic curve based on digital experiment of point clouds

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