CN105528474A - 飞机样板激光切割编程图形的碎线拟合圆弧方法 - Google Patents
飞机样板激光切割编程图形的碎线拟合圆弧方法 Download PDFInfo
- Publication number
- CN105528474A CN105528474A CN201510839598.8A CN201510839598A CN105528474A CN 105528474 A CN105528474 A CN 105528474A CN 201510839598 A CN201510839598 A CN 201510839598A CN 105528474 A CN105528474 A CN 105528474A
- Authority
- CN
- China
- Prior art keywords
- circular arc
- fitting
- point
- circle
- array
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Numerical Control (AREA)
Abstract
本发明为一种飞机样板激光切割编程图形的碎线拟合圆弧方法,实现该方法的程序基于ObjecARX(VC)开发。方法为先分解选定的CAD二维图形中的所有元素,提取分解所得元素的特征点。然后对这些点进行最小二乘拟合,得到拟合的圆心、半径结果。之后根据最初选择的碎线在该圆上截得一个能拟合这些碎线的劣弧,计算拟合公差,如果公差在允许范围内则拟合成功。本发明能将碎线快速拟合成一个完整的圆弧,并能精确的控制误差。与现有技术对比,本发明所得结果的图形质量更高,适用范围更广,能从根本上杜绝各种拟合错误,还能处理现行技术无法处理的特殊情况。
Description
技术领域
本发明为一种飞机样板激光切割编程图形的碎线拟合方法,属于飞机数字化制造技术领域。
背景技术
飞机样板是一类将飞机部、组、零件的某些特征以1:1的比例制成且带有工艺信息的钢性量具,其制造简单,成本低廉,使用方便,因而在飞机制造过程中得到了广泛的应用。制造飞机样板,需要先在三维数模中将要表达的零件特征转化为二维图形,将此二维图形进行一系列处理后导入专门的软件中自动生成激光切割程序。如果图形中存在大量连续的短线(以下成为碎线,如图1)会造成激光切割机卡顿,缩短激光切割机的使用寿命。对碎线的处理中,尤以呈圆弧状排列的碎线难度最大,因为圆弧的起止点、半径等参数均需肉眼判断,误差难以控制,可靠性差。目前业内使用的图形处理程序SYLBDxfProcess中提供了类似的拟合碎线功能,试图解决这类问题,但此程序的处理结果并未改变图形含有大量连续短线的本质,有若干常见情形无法处理,并且偶尔会出现误差超过允许范围的结果,使用完毕后需要人工校对,性能与实际需求相差甚远。
发明内容
本发明提供一种飞机样板激光切割编程图形的碎线拟合圆弧方法,目的是提高拟合精度及准确性。
本发明是通过以下技术方案实施的:飞机样板激光切割编程图形的碎线拟合圆弧方法,基于ObjecARX(VC)开发,具体方法为:
(1)获得选定的图形元素,分解其中的组合元素,具体过程为:
①提示用户选择,当用户完成选择后,将选择集中的图形的ID加入实体ID数组;
②如果实体ID数组不为空,依次打开实体ID数组中的每个图形,分解该图形;
③如果分解成功,将该图形的ID从实体ID数组中移除,并将分解的结果加入到实体ID数组中;
(2)获取选择集中所有元素的特征点,具体过程为:
①遍历实体ID数组,依次打开每一个图形,如果图形是曲线,即直线或圆弧,则取其两个端点和中点放入点数组中;
②建立内外两层循环,将外循环的点与内循环的逐一比较,如果发现坐标值相等,将该点移出数组;
(3)用特征点拟合圆,调整圆心与半径,根据碎线信息截得圆弧,具体过程为:
①使用最小二乘法获得拟合特征点的圆;
②所有特征点依次与圆心构建矢量,方向沿半径向外,根据这些矢量之间的最小正角判断出圆弧的两个端点,并用这两个端点定位圆心;
③获得由②所得的圆与两端点确定的劣弧的起止角与终止角,创建圆弧;
(4)判断拟合误差,具体过程为:
①获得拟合圆弧的起点、中点与终点,构建形状相同的几何类对象;
②对每一个特征点,计算点与该几何类对象的最小距离,并取其最大值为误差值;
③如果误差值小于等于行业规定的允许值,则删除(1)中选择的所有元素,将创建的圆弧发送至模型空间并突出显示,同时发送拟合误差信息;否则关闭创建的圆弧,发送拟合失败信息。
本发明的有益效果:本发明采用上述方案,可以快速地将圆弧状排列的碎线拟合成整体的圆弧(如图2),而不是由多段直线或圆弧拼接而成;对于直线与圆弧相切的情况,不改变切点位置;拟合不会出现误差大于行业规定最大误差的结果;对于局部法向误差过大、局部断开等情形,仍能给出符合上述特性的结果。
本发明改变了现行技术拟合结果仍是多条短线的尴尬局面(如图3),拟合精度极高,并且能由程序轻而易举地控制误差,避免了人为控制误差而引起的超差现象。另外,本发明可以处理诸如局部断开之类的特殊情况,弥补了现行技术该方面的空白。本发明在适用性与准确性上有了质的提高,为飞机样板激光切割编程图形的碎线拟合圆弧提高了可靠的解决方案。
附图说明
图1为需要拟合的圆弧形碎线。
图2为本发明对图1的处理结果。
图3为现有技术对图1的处理结果。
具体实施方式
一种飞机样板激光切割编程图形的碎线拟合圆弧方法,具体方法为:
(1)获得选定的图形元素,分解其中的组合元素;
(2)获取选择集中所有元素的特征点;
(3)用特征点拟合圆,调整圆心与半径,根据碎线信息截得圆弧;
(4)判断拟合误差。
所述的步骤(1)获得所选择图形元素,分解其中的组合元素:在命令激活后提示用户选择,获得选择集,建立循环,将获得所有实体的ID直接加入ID数组。在这些实体中,难免会有块、多段线等组合元素,应对其进行分解。对获得的所有实体执行分解动作,组合元素的分解会成功。如果分解成功,则将所得的实体依次发送至模型空间,并将ID加入之前的ID数组,然后在ID数组中移除该组合元素的ID,并将该组合元素从模型空间删除。
步骤(2)获取选择集中所有元素的特征点:依次以“读”模式打开ID数组中的所有元素,打开的图形必须是曲线(直线或圆弧,下同)。由于本发明拟合采用的方法是最小二乘法,需要用点来构造曲线,而用于拟合碎线的圆弧,其波动范围应由这些碎线的最小包围盒确定。对于直线,其最小包围盒只由起点和终点确定,对于圆弧则由起点、中点和终点确定。然而如果圆弧取三个点而直线取两个点,会造成直线的权值小于圆弧,这在实际情况下是不合理的,所以不论直线还是圆弧,均取其起点、中点和终点放入点数组中;
步骤(3)用特征点拟合圆,调整圆心与半径,根据碎线信息截得圆弧:最小二乘拟合的任务是根据N个点(i=0,1,……,n-1)得出圆心(A,B)和半径R。计算过程需要定义一下几个量:
∑X1=X0+X1+X2+……+Xn-1
∑X1 2=X0 2+X1 2+X2 2+……X2 n-1
∑X1 3=X0 3+X1 3+X2 3+……X3 n-1
∑Y1=Y0+Y1+Y2+……+Yn-1
∑Y1 2=Y0 2+Y1 2+Y2 2+……Y2 n-1
∑Y1 3=Y0 3+Y1 3+Y2 3+……Y3 n-1
∑X1Y1=X0Y0+X1Y1+X2Y2+……Xn-1Yn-1
∑X2 1Y1=X2 0Y0+X2 1Y1+X2 2Y2+……X2 n-1Yn-1
∑X1Y2 1=X0Y2 0+X1Y2 1+X2Y2 2+……Xn-1Y2 n-1
由以上变量构建参数C、D、E、G和H:
C=N*∑X1 2-∑X1*∑X1
D=N*∑X1Y1-∑X1*∑Y1
E=N*∑X1 3+V*∑X1Y2 1-(∑X1 2+∑Y1 2)*∑X1
G=N*∑Y1 2-∑Y1*∑Y1
H=N*∑X2 1Y1+N*∑Y1 3-(∑X1 2+∑Y1 2)*∑Y1
由以上参数构建参数a、b、c:
a=(H*D-E*G)/(C*G-D*D)
b=(H*C-E*D)/(D*D-G*C)
c=-(a*∑X1+b*∑Y1+∑X2 1+∑X2 1)/N
则
A=a/(-2.0)
B=b/(-2.0)
R=sqrt(a*a+b*b-4*c)/2.0
根据圆心和特征点可以在上述的拟合圆上截得一个劣弧。按特征点在数组里的顺序,首先构建圆心到第一个特征点的矢量v0,然后用之后的特征点与圆心构建矢量与v0计算夹角,找到夹角最大的那个矢量v1,则点p1就是弧的一点端点。同理,将上述所有矢量与v1计算夹角并比较即可得到弧的另一个端点。这里注意圆是那和得到的,意味着该圆可能不会通过特征点。考虑到切点位置不能改变,应修正正圆心的参数。方法好似按照拟合的半径以两端点为圆心做成两个圆,如果两圆有2个交点(否则圆心不需修正),则离当前圆心最近的一个点就是新的圆心。至此已经可以建立拟合圆弧。
步骤(4)判断拟合误差:拟合误差是指碎线上各点与拟合圆弧的最小距离中的最大值。由于该最大值必定出现在特征点上,故只计算各个特征点与拟合圆弧的距离即可。如果误差小于等于行业规定的允许值,意味着拟合结果符合要求,需要用拟合图形来代替原图形,即删除所选择的碎线,并将拟合圆弧发送至模型空间;如果拟合超差,出于图形数据库的管理要求仍需要关闭拟合圆弧。
Claims (5)
1.飞机样板激光切割编程图形的碎线拟合圆弧方法,其特征在于实现该方法的程序基于ObjecARX(VC)开发,具体方法为:
(1)获得选定的图形元素,分解其中的组合元素;
(2)获取选择集中所有元素的特征点;
(3)用特征点拟合圆,调整圆心与半径,根据碎线信息截得圆弧;
(4)判断拟合误差。
2.如权利要求1所述的飞机样板激光切割编程图形的碎线拟合圆弧方法,其特征在于:所述步骤(1)具体过程为:
①提示用户选择,当用户完成选择后,将选择集中的图形的ID加入实体ID数组;
②如果实体ID数组不为空,依次打开实体ID数组中的每个图形,分解该图形;
③如果分解成功,将该图形的ID从实体ID数组中移除,并将分解的结果加入到实体ID数组中。
3.如权利要求1所述的飞机样板激光切割编程图形的碎线拟合圆弧方法,其特征在于:所述步骤(2)具体过程为:
①遍历实体ID数组,依次打开每一个图形,如果图形是曲线,即直线或圆弧,则取其两个端点和中点放入点数组中;
②建立内外两层循环,将外循环的点与内循环的逐一比较,如果发现坐标值相等,将该点移出数组。
4.如权利要求1所述的飞机样板激光切割编程图形的碎线拟合圆弧方法,其特征在于:所述步骤(3)具体过程为:
①使用最小二乘法获得拟合特征点的圆;
②所有特征点依次与圆心构建矢量,方向沿半径向外,根据这些矢量之间的最小正角判断出圆弧的两个端点,并用这两个端点定位圆心;
③获得由②所得的圆与两端点确定的劣弧的起止角与终止角,创建圆弧。
5.如权利要求1所述的飞机样板激光切割编程图形的碎线拟合圆弧方法,其特征在于:所述步骤(4)具体过程为:
①获得拟合圆弧的起点、中点与终点,构建形状相同的几何类对象;
②对每一个特征点,计算点与该几何类对象的最小距离,并取其最大值为误差值;
③如果误差值小于等于行业规定的允许值,则删除(1)中选择的所有元素,将创建的圆弧发送至模型空间并突出显示,同时发送拟合误差信息;否则关闭创建的圆弧,发送拟合失败信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510839598.8A CN105528474A (zh) | 2015-11-27 | 2015-11-27 | 飞机样板激光切割编程图形的碎线拟合圆弧方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510839598.8A CN105528474A (zh) | 2015-11-27 | 2015-11-27 | 飞机样板激光切割编程图形的碎线拟合圆弧方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105528474A true CN105528474A (zh) | 2016-04-27 |
Family
ID=55770696
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510839598.8A Pending CN105528474A (zh) | 2015-11-27 | 2015-11-27 | 飞机样板激光切割编程图形的碎线拟合圆弧方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105528474A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109741262A (zh) * | 2019-01-07 | 2019-05-10 | 凌云光技术集团有限责任公司 | 一种基于位置关系的轮廓图像拼接方法 |
CN110348163A (zh) * | 2019-07-19 | 2019-10-18 | 合肥工业大学 | 一种轮胎花纹边界连续相切圆弧的拟合方法 |
CN111788604A (zh) * | 2019-02-15 | 2020-10-16 | 深圳配天智能技术研究院有限公司 | 图像处理方法、系统、电子设备、机器人和存储装置 |
CN113857789A (zh) * | 2021-10-26 | 2021-12-31 | 杨官全 | 一种高精度异形曲面的加工方法 |
CN114819306A (zh) * | 2022-04-13 | 2022-07-29 | 南通恒康数控机械股份有限公司 | 一种应用于海绵切割文件的拟合方法 |
CN115828398A (zh) * | 2023-01-10 | 2023-03-21 | 深圳市大乐装建筑科技有限公司 | 碎线数据处理方法、装置、设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102000916A (zh) * | 2010-08-31 | 2011-04-06 | 深圳市大族激光科技股份有限公司 | 一种激光切割方法 |
US20120109358A1 (en) * | 2005-11-02 | 2012-05-03 | Goldman David A | Printer driver systems and methods for automatic generation of embroidery designs |
-
2015
- 2015-11-27 CN CN201510839598.8A patent/CN105528474A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120109358A1 (en) * | 2005-11-02 | 2012-05-03 | Goldman David A | Printer driver systems and methods for automatic generation of embroidery designs |
CN102000916A (zh) * | 2010-08-31 | 2011-04-06 | 深圳市大族激光科技股份有限公司 | 一种激光切割方法 |
Non-Patent Citations (3)
Title |
---|
傅强: "图形交互式数控激光切割自动编程系统的研发", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
姚必强等: "数控加工曲线的等弧长圆弧拟合方法", 《四川大学学报(工程科学版)》 * |
贾东海: "平板件数控下料编程前处理技术研究与实现", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109741262A (zh) * | 2019-01-07 | 2019-05-10 | 凌云光技术集团有限责任公司 | 一种基于位置关系的轮廓图像拼接方法 |
CN111788604A (zh) * | 2019-02-15 | 2020-10-16 | 深圳配天智能技术研究院有限公司 | 图像处理方法、系统、电子设备、机器人和存储装置 |
CN110348163A (zh) * | 2019-07-19 | 2019-10-18 | 合肥工业大学 | 一种轮胎花纹边界连续相切圆弧的拟合方法 |
CN110348163B (zh) * | 2019-07-19 | 2022-09-13 | 合肥工业大学 | 一种轮胎花纹边界连续相切圆弧的拟合方法 |
CN113857789A (zh) * | 2021-10-26 | 2021-12-31 | 杨官全 | 一种高精度异形曲面的加工方法 |
CN113857789B (zh) * | 2021-10-26 | 2024-03-22 | 杨官全 | 一种高精度异形曲面的加工方法 |
CN114819306A (zh) * | 2022-04-13 | 2022-07-29 | 南通恒康数控机械股份有限公司 | 一种应用于海绵切割文件的拟合方法 |
CN115828398A (zh) * | 2023-01-10 | 2023-03-21 | 深圳市大乐装建筑科技有限公司 | 碎线数据处理方法、装置、设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105528474A (zh) | 飞机样板激光切割编程图形的碎线拟合圆弧方法 | |
US9984181B2 (en) | Junction meshing for lattice structures | |
CN109976262B (zh) | 一种针对微线段加工的全局曲率连续光顺方法 | |
Han et al. | Iso-phote based tool-path generation for machining free-form surfaces | |
KR20070017070A (ko) | 요구되는 기하학적 연속성을 가지는 파라메트릭 표면 생성프로세스 | |
CN110703685B (zh) | 蒙皮镜像铣削加工刀具轨迹自适应调整方法及装置 | |
CN104898554A (zh) | 一种基于离散刀位点的复合刀具路径生成方法 | |
CN108073138B (zh) | 适用于高速高精加工的椭圆弧平滑压缩插补算法 | |
JP2017205975A (ja) | 3次元データ生成装置、3次元造形装置、造形物の製造方法及びプログラム | |
CN102794488A (zh) | 类直纹面整体叶轮曲面的侧铣加工方法 | |
CN113953685A (zh) | 一种平面板材激光切割路径规划方法 | |
CN105867310A (zh) | 基于occ实现t样条模型的数控加工方法及系统 | |
CN101515168A (zh) | 产品三角网格模型的环切粗加工nc刀轨生成方法 | |
Chatelain et al. | A balancing technique for optimal blank part machining | |
CN114943144A (zh) | 利用Phi函数进行距离控制的卫星布局优化设计方法 | |
EP0240569A1 (en) | Method of forming compposite curved surface | |
CN112613830B (zh) | 一种物资储备中心选址方法 | |
CN116339242B (zh) | 一种自由曲面刀具路径生成方法及相关设备 | |
WO2019087032A1 (en) | Method for the reconstruction of cad models through parametric adaptation | |
CN107403469B (zh) | 一种提高折角成型质量的自适应扫描速度方法 | |
CN116484635A (zh) | 基于改进圆弧拟合算法的正弦弧线拟合局部光顺方法 | |
US9152743B2 (en) | Computer process for determining best-fitting materials for constructing architectural surfaces | |
Chan et al. | Volume decomposition of CAD models for rapid prototyping technology | |
CN112507410B (zh) | 轨道梁图纸的生成方法和生成装置 | |
US8872819B2 (en) | Computational geometry design for sheet metal machinery corner |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160427 |