CN1214232C - 一种自由曲面的测量方法 - Google Patents
一种自由曲面的测量方法 Download PDFInfo
- Publication number
- CN1214232C CN1214232C CN 03121461 CN03121461A CN1214232C CN 1214232 C CN1214232 C CN 1214232C CN 03121461 CN03121461 CN 03121461 CN 03121461 A CN03121461 A CN 03121461A CN 1214232 C CN1214232 C CN 1214232C
- Authority
- CN
- China
- Prior art keywords
- point
- cross
- curved surface
- measuring method
- measuring
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明涉及一种自由曲面的测量方法,该方法利用集成化柔性激光加工系统、三座标测量机或数控机床,根据自由曲面的特征将自由曲面划分成等距截面;然后,对各个截面曲线进行测量,即先在每个截面的曲线上用自动等距离步长测量方法测量该曲线上四个起始点P1、P2、P3、P4的三维空间坐标,计算出直线P2P3和直线P3P4之间的锐角θ,再根据已测得的四点P1、P2、P3、P4用三次贝塞尔方程拟合曲线,并以θ为判距自适应的调整变参数Δt,获取P5的空间三维坐标值,直至边界点;最后,测量出自由曲面。该方法在测量过程中能不断修正前面不合理的测量点,不但测量效率高,而且,以较高精度实现了自由曲面的数字化。
Description
技术领域
本发明涉及自动控制测量技术领域,特别涉及一种自由曲面的测量方法。
背景技术
由于自由曲面没有真正意义上的解析式,故其数字化一直是工程技术中的难点。目前,有些学者根据质心定理提出了对数学模型已知的自由曲面的测量方法,但对于未知数学模型的自由曲面的测量,则一般通过等距测量方法,即用等步长的方式来对自由曲面进行测量,其缺陷在于没有考虑曲面本身的特点,即使是对于变化比较剧烈的曲面也用要和平坦的曲面一样的步长,这样获得的测量数据有时就不能反映自由曲面的真实情况,并且测量效率较低。
发明内容
本发明的目的在于:克服用等距测量方法测量未知数学模型的自由曲面时不能反映曲面的真实情况,且测量效率低的缺点,从而提供一种自由曲面的测量方法。
本发明的目的是这样实现的:本发明的一种自由曲面的测量方法,利用集成化柔性激光加工系统、三座标测量机或数控机床,根据自由曲面的特征将自由曲面划分成等距截面,然后,对各个截面曲线进行测量,自适应地生成下一个探测点和对应的避障点,并能在测量过程中不断修正前面不合理的测量点,以测量出自由曲面。
本发明的一种自由曲面的测量方法,利用集成化柔性激光加工系统、三座标测量机或数控机床,其特征在于,包括以下步骤:
(1)确定待测曲面的测量范围,并将测量曲面等分为N个截面,所述的N为所需要测量的排数;设置截面序号变量m,所述的m为设置内存变量截面序号;
(2)在第m个截面的曲线上用自动等距离步长测量方法测量该曲线上四个起始点P1、P2、P3、P4的三维空间坐标,并将四点的坐标值存入计算机内的测量数据文件中;
(3)根据已测得的四点P1、P2、P3、P4,计算出直线P2P3和直线P3P4之间的锐角θ;
(4)再根据已测得的四点P1、P2、P3、P4用三次贝塞尔方程拟合曲线,再以θ为判距自适应的调整变参数Δt,获取下一个探测点P′5,并求出其对应的避障点Pb,5,利用探测点P′5和避障点Pb,5测得P5的空间三维坐标值,并将坐标值存入测量数据文件中,如果变参数Δt>0,则将P2、P3、P4、P5赋于P1、P2、P3、P4;
(5)如果变参数Δt<0,说明P5点是回测点,即P4点是在P5点的前面,则用P1、P2、P3、P5点的值测量下一个探测点和避障点;
(6)重复上述步骤(3)、(4)和(5)直至边界,并令截面序号值m=m+1;当m=N时,结束测量。
所述的判距θ≤1°,则Δt=1/3+0.10;若1°<θ≤3°,则Δt=1/3+0.05;若3°<θ≤6°,则Δt=1/3;若θ>6°,则Δt=-1/6。
所述的三次贝塞尔方程为P′5=T(1+Δt)MP,其中T(1+Δt)=[(1+Δt)3(1+Δt)2(1+Δt)1],是变参数矩阵,M为三次贝塞尔基函数的常数系数矩阵,
P为最新测量获得的四个点的空间三维坐标的矩阵,且P=[p1 p2 p3 p4]T。
所述的排数N为>2的正整数。
所述的截面序号m为>2的正整数。
本发明的优点在于:本发明的一种自由曲面的测量方法,基于特定的三次贝塞尔数学模型,自适应测量等距截面上的曲线,不但测量效率高,而且,能反映曲面的真实情况,以较高精度实现了自由曲面的数字化。
图面说明
图1是本发明自由曲面测量方法中边界点的测量示意图
图2是本发明自由曲面测量方法中的曲线测量示意图
图3是本发明自由曲面测量方法实施例的测量结果数据图
图4是本发明自由曲面测量方法的流程示意图
具体实施方式
利用本发明自由曲面测量方法的一个最佳实施例对本发明作出进一步说明。
采用中科院力学研究所已授权的中国专利号为ZL98101217.5的一种具有柔性传输和多轴联动的激光加工装置,也可采用三坐标测量机、利用接触式测头(或非接触式测头)并配以测量辅助电路与龙门式数控机床的I/O通讯口相连而改造成的测量加工机等设备,对汽车覆盖件冲压模具的一个角部的自由曲面进行了测量,如图4所示,其具体的测量过程如下:
1、先确定待测曲面的测量范围,如图1所示,沿Y轴正向看去,用测量头人工测量出待测曲面的左下角边界点Ps,其三维空间坐标值为(2072.342-1324.306-314.720);左上角边界点Pe,其三维空间坐标值为(2072.342-1225.630-286.319);右上角边界点Pc,获得其空间三维坐标值为(2145.276-1225.630-260.691)。
根据测量要求的精度确定自由曲面的每两个截面之间的距离d=4mm,计算出所需要测量的排数N=12,设置内存变量截面序号m,对于第一个要测量的截面,取m值为1。
2、从截面初始点开始,如图2所示,控制测量头用等距步长为4mm的自动等距步长测量方法测量得到待测自由曲面第m个截面上的四点P1、P2、P3、P4的三维空间坐标,将P1、P2、P3、P4为软件设置为全局内存变量,并将四点空间三维坐标值由测量控制软件按文本方式存入测量数据文件中;
3、根据测得的四点P1、P2、P3、P4,利用公式
计算出直线P2P3和直线P3P4之间的锐角θ,再以θ为判距,自适应的调整变参数变量Δt:若θ≤1°,则Δt=1/3+0.10;若1°<θ≤3°,Δt=1/3+0.05;若3°<θ≤6°,Δt=1/3;若θ>6°,Δt=-1/6,回测一点。
利用公式 P′5=T(1+Δt)MP (3)
P′b,5=P′5+(0 0 10) (4)
生成下一个探测点P′5和避障点Pb,5;
4、利用探测点P′5和避障点Pb,5驱动测头获得测量点P5的空间三维坐标值,并将坐标值存入文本格式的测量数据文件中;
5、如果变参数Δt>0,则将P2、P3、P4、P5赋于P1、P2、P3、P4;如果变参数Δt<0,说明P5点是回测点,即P4点是在P5点的前面,则用P1、P2、P3、P5点的值测量下一个探测点和避障点;
6、计算P5与Pe之间的距离DS,若DS<0.5mm则结束测量;否则重复步骤3~5直至结束;并令截面序号值m=m+1,进行下一个截面的测量。
当m=N时,结束测量。最后获得的测量结果采用origin绘图软件显示,如图3所示。
Claims (5)
1、一种自由曲面的测量方法,利用集成化柔性激光加工系统、三座标测量机或数控机床,其特征在于,包括以下步骤:
(1)确定待测曲面的测量范围,并将测量曲面等分为N个截面,所述的N为所需要测量的排数;设置截面序号变量m,所述的m为设置内存变量截面序号;
(2)在第m个截面的曲线上用自动等距离步长测量方法测量该曲线上四个起始点P1、P2、P3、P4的三维空间坐标,并将四点的坐标值存入计算机内的测量数据文件中;
(3)根据已测得的四点P1、P2、P3、P4,计算出直线P2P3和直线P3P4之间的锐角θ;
(4)再根据已测得的四点P1、P2、P3、P4用三次贝塞尔方程拟合曲线,再以θ为判距自适应的调整变参数Δt,获取下一个探测点P′5,并求出其对应的避障点Pb,5,利用探测点P′5和避障点Pb,5测得P5的空间三维坐标值,并将坐标值存入测量数据文件中,如果变参数Δt>0,则将P2、P3、P4、P5赋于P1、P2、P3、P4;
(5)如果变参数Δt<0,说明P5点是回测点,即P4点是在P5点的前面,则用P1、P2、P3、P5点的值测量下一个探测点和避障点;
(6)重复上述步骤(3)、(4)和(5)直至边界,并令截面序号值m=m+1;当m=N时,结束测量。
2.按照权利要求1所述的自由曲面的测量方法,其特征在于,以θ为判距,自适应的调整变参数变量Δt:若θ≤1°,则Δt=1/3+0.10;若1°<θ≤3°,Δt=1/3+0.05;若3°<θ≤6°,Δt=1/3;若θ>6°,Δt=-1/6,回测一点。
3、按照权利要求1所述的自由曲面的测量方法,其特征在于,所述的三次贝塞尔方程为P5=T(1+Δt)MP,其中T(1+Δt)=[(1+Δt)3(1+Δt)2(1+Δt)1],是变参数矩阵,M为三次贝塞尔基函数的常数系数矩阵,
P为最新测量获得的四个点的空间三维坐标的矩阵,且P=[p1p2p3p4]T。
4.按照权利要求1所述的自由曲面的测量方法,其特征在于,所述的排数N为>2的正整数。
5.按照权利要求1所述的自由曲面的测量方法,其特征在于,所述的截面序号m为>2的正整数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 03121461 CN1214232C (zh) | 2003-03-28 | 2003-03-28 | 一种自由曲面的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 03121461 CN1214232C (zh) | 2003-03-28 | 2003-03-28 | 一种自由曲面的测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1534276A CN1534276A (zh) | 2004-10-06 |
CN1214232C true CN1214232C (zh) | 2005-08-10 |
Family
ID=34285713
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 03121461 Expired - Fee Related CN1214232C (zh) | 2003-03-28 | 2003-03-28 | 一种自由曲面的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1214232C (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100387931C (zh) * | 2006-04-28 | 2008-05-14 | 南通大学 | 实现大量程自由曲面的高精度测量方法 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100418027C (zh) * | 2006-11-10 | 2008-09-10 | 大连理工大学 | 一种螺旋线制导的曲面数控加工方法 |
CN101738175B (zh) * | 2008-11-21 | 2014-09-10 | 红塔烟草(集团)有限责任公司 | 一种以测针球心坐标为目标点测量空间旋转曲面的方法 |
CN102062593B (zh) * | 2010-11-12 | 2012-11-07 | 北京信息科技大学 | 螺旋曲面测点获取方法 |
CN102501172B (zh) * | 2011-11-22 | 2013-10-23 | 清华大学 | 用于机器人修磨系统的面向空间曲面加工的在位测量方法 |
CN103808293B (zh) * | 2012-11-14 | 2018-08-21 | 中国广核集团有限公司 | 监测核电站热功率测量漂移中的孔板检测方法 |
CA2956518C (en) * | 2014-07-30 | 2021-03-30 | Ysystems Ltd. | Method and apparatus for measuring surface profile |
CN106092005A (zh) * | 2016-06-03 | 2016-11-09 | 哈尔滨东安发动机(集团)有限公司 | 三坐标机测量自由曲线的方法 |
CN107091627B (zh) * | 2017-04-18 | 2020-07-28 | 复旦大学 | 一种多表面系统的综合测量与评估方法 |
CN108225243B (zh) * | 2017-12-28 | 2019-09-10 | 中材科技股份有限公司 | 一种快速获取异形曲面结构厚度分布的方法 |
-
2003
- 2003-03-28 CN CN 03121461 patent/CN1214232C/zh not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100387931C (zh) * | 2006-04-28 | 2008-05-14 | 南通大学 | 实现大量程自由曲面的高精度测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1534276A (zh) | 2004-10-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1214232C (zh) | 一种自由曲面的测量方法 | |
CN100387931C (zh) | 实现大量程自由曲面的高精度测量方法 | |
CN110348125B (zh) | 一种块体稳定性的可视化判别方法 | |
CN110530296B (zh) | 一种线激光安装误差角确定方法 | |
CN110487210B (zh) | 蜂窝芯表面轮廓测量方法 | |
CN105066912B (zh) | 酸刻蚀物理模拟实验中岩板表面扫描数据的步长标定方法 | |
CN1991846A (zh) | 一种生成三维有限元网格的方法 | |
CN106875439A (zh) | 基于三维点云模型的单晶硅棒外形尺寸测量方法 | |
CN103884291A (zh) | 基于nurbs参数曲面的建筑物表面柔性变形监测方法 | |
CN105241392A (zh) | 一种复杂柱状工件的全表面三维测量设备及其测量方法 | |
CN104089599B (zh) | 用于接触式测头测量中提取二维轮廓的准形态学滤波方法 | |
CN202204481U (zh) | 齿轮倒角激光测量仪 | |
CN1217154C (zh) | 一种利用单光束激光进行准直/轴对中测量的方法 | |
CN209181760U (zh) | 一种基于双非平行面直线导轨的平行度测量装置 | |
CN101498572B (zh) | 凸轮二维曲线测量方法 | |
CN2903906Y (zh) | 一种长度和弧度测量工具 | |
CN1306460C (zh) | 一种曲面测量方法 | |
Costa | Fractal description of rough surfaces for haptic display | |
CN1758164A (zh) | 基于压缩体素模型的棒刀空间扫描体构造方法 | |
CN110411341A (zh) | 含螺纹目标物的位姿计算方法 | |
CN1214231C (zh) | 一种空间三维自由曲线的测量方法 | |
CN1831473A (zh) | 基于垂直位移扫描的接触式表面形貌测量方法及测量仪 | |
CN114707265B (zh) | 五轴球头铣削加工表面纹理形貌的评价指标的提取方法 | |
CN112884896A (zh) | 基于bim技术的异形坡道三维模型快速生成方法 | |
CN208171176U (zh) | 一种微型直线光栅传感器 |
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 | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |