标题

Similar documents
g 100mv /g 0. 5 ~ 5kHz 1 YSV8116 DASP 1 N 2. 2 [ M] { x } + [ C] { x } + [ K]{ x } = { f t } 1 M C K 3 M C K f t x t 1 [ H( ω )] = - ω 2

PCA+LDA 14 1 PEN mL mL mL 16 DJX-AB DJ X AB DJ2 -YS % PEN

untitled

第 05 期 董房等 : 一种卫星遥测在线状态监测及分析系统的设计 WEB 1 2 总体功能及组成 2.1 总体功能 1 2 3Web 2.2 结构组成 Web WEB WEB 2.3 系统各模块接口关系

SVM OA 1 SVM MLP Tab 1 1 Drug feature data quantization table

[1-3] (Smile) [4] 808 nm (CW) W 1 50% 1 W 1 W Fig.1 Thermal design of semiconductor laser vertical stack ; Ansys 20 bar ; bar 2 25 Fig

标题

Fig. 1 1 The sketch for forced lead shear damper mm 45 mm 4 mm 200 mm 25 mm 2 mm mm Table 2 The energy dissip

mm ~

18-陈亚莉.FIT)

* CUSUM EWMA PCA TS79 A DOI /j. issn X Incipient Fault Detection in Papermaking Wa

Mixtions Pin Yin Homepage

T K mm mm Q345B 600 mm 200 mm 50 mm 600 mm 300 mm 50 mm 2 K ~ 0. 3 mm 13 ~ 15 mm Q345B 25

<D2BDC1C6BDA1BFB5CDB6C8DAD7CAB8DFB7E5C2DBCCB3B2CEBBE1C3FBB5A52E786C7378>

~ 4 mm h 8 60 min 1 10 min N min 8. 7% min 2 9 Tab. 1 1 Test result of modified

#4 ~ #5 12 m m m 1. 5 m # m mm m Z4 Z5

/MPa / kg m - 3 /MPa /MPa 2. 1E ~ 56 ANSYS 6 Hz (a) 一阶垂向弯曲 (b) 一阶侧向弯曲 (c) 一阶扭转 (d) 二阶侧向弯曲 (e) 二阶垂向弯曲 (f) 弯扭组合 2 6 Hz

SWAN min TITAN Thunder Identification Tracking Analysis SWAN TITAN and Nowcasting 19 TREC Tracking Radar Echo by Correlaction T

10 中 草 药 Chinese Traditional and Herbal Drugs 第 43 卷 第 1 期 2012 年 1 月 生 药 打 粉 入 药 的 基 本 特 点, 借 鉴 材 料 学 粉 体 学 等 学 科 的 研 究 成 果, 在 中 药 传 统 制 药 理 念 的 启 发

2 3. 1,,,.,., CAD,,,. : 1) :, 1,,. ; 2) :,, ; 3) :,; 4) : Fig. 1 Flowchart of generation and application of 3D2digital2building 2 :.. 3 : 1) :,

Microsoft Word 定版

KJWJ01 Fig. 1 Geological map of the Kajiwa landslide PDH01 60 ~ 70m 80 ~ 90m 80m 2 930m 556m 641m 598m m m 3 2 Ⅰ m 2

θ 1 = φ n -n 2 2 n AR n φ i = 0 1 = a t - θ θ m a t-m 3 3 m MA m 1. 2 ρ k = R k /R 0 5 Akaike ρ k 1 AIC = n ln δ 2

Technical Acoustics Vol.27, No.4 Aug., 2008,,, (, ) :,,,,,, : ; ; : TB535;U : A : (2008) Noise and vibr

Microsoft Word - 11-秦华伟.doc

1602 CHEMICAL INDUSTRY AND ENGINEERING PROGRESS TQ A Prospect of the structure of oxygen generating


1.0 % 0.25 % 85μm % U416 Sulfate expansion deformation law and mechanism of cement stabilized macadam base of saline areas in Xinjiang Song

标题

f 2 f 2 f q 1 q 1 q 1 q 2 q 1 q n 2 f 2 f 2 f H = q 2 q 1 q 2 q 2 q 2 q n f 2 f 2 f q n q 1 q n q 2 q n q n H R n n n Hessian

(Pattern Recognition) 1 1. CCD

西部人居环境学刊2016年第1期.indd

Microsoft Word - 刘 慧 板.doc

自然科学版 预处理 视盘粗定位 视盘垂直坐标的粗定位 视盘水平坐标的粗定位

标题

Fig. 1 Frame calculation model 1 mm Table 1 Joints displacement mm

~ 10 2 P Y i t = my i t W Y i t 1000 PY i t Y t i W Y i t t i m Y i t t i 15 ~ 49 1 Y Y Y 15 ~ j j t j t = j P i t i = 15 P n i t n Y

1

[1] [4] Chetverikov Lerch[8,12] LeaVis CAD Limas-Serafim[6,7] (multi-resolution pyramids) 2 n 2 n 2 2 (texture) (calf leather) (veins)

标题

1 引言

JOURNAL OF EARTHQUAKE ENGINEERING AND ENGINEERING VIBRATION Vol. 31 No. 5 Oct /35 TU3521 P315.

Stock Transfer Service Inc. Page No. 1 CENTURY PEAK METALS HOLDINGS CORPORATION (CPM) List of Top 100 Stockholders As of 12/31/2015 Rank Sth. No. Name

Microsoft Word - KSAE06-S0262.doc

5 551 [3-].. [5]. [6]. [7].. API API. 1 [8-9]. [1]. W = W 1) y). x [11-12] D 2 2πR = 2z E + 2R arcsin D δ R z E = πr 1 + πr ) 2 arcsin

第 2 期 王 向 东 等 : 一 种 运 动 轨 迹 引 导 下 的 举 重 视 频 关 键 姿 态 提 取 方 法 257 竞 技 体 育 比 赛 越 来 越 激 烈, 为 了 提 高 体 育 训 练 的 效 率, 有 必 要 在 体 育 训 练 中 引 入 科 学 定 量 的 方 法 许 多

报 告 1: 郑 斌 教 授, 美 国 俄 克 拉 荷 马 大 学 医 学 图 像 特 征 分 析 与 癌 症 风 险 评 估 方 法 摘 要 : 准 确 的 评 估 癌 症 近 期 发 病 风 险 和 预 后 或 者 治 疗 效 果 是 发 展 和 建 立 精 准 医 学 的 一 个 重 要 前

标题

10期( )

基 础 实 室 4 计 算 机 网 络 唐 爱 红 专 业 机 房 PROTEL 联 想 同 方 电 脑 180 台 唐 爱 红 MATLAB 计 算 机 网 络 电 工 电 子 技 能 训 练 室 电 子 基 本 技 能 示 波 器 毫 伏 表 雕 刻 机 图 示 仪 电 子 实 训 台 电 工

KUKA W. Polini L. Sorrentino Aized Shirinzadeh 6 7 MF Tech Pitbull Fox Taniq Scorpo Scorpo Compositum Windows KUKA 1 P 1 P 2 KU

Microsoft Word - 专论综述1.doc

% GIS / / Fig. 1 Characteristics of flood disaster variation in suburbs of Shang

专 技 能 1. 精 通 Matlab/Simulink 平 台 下 的 海 洋 运 载 器 运 动 控 制 系 统 与 仿 真 建 模 设 计 ; 2. 精 通 51 单 片 机 AVR 单 片 机 Arduino 开 源 板 的 开 发 和 设 计 ; 3. 精 通 基 于 Arduino 板

4 115,,. : p { ( x ( t), y ( t) ) x R m, y R n, t = 1,2,, p} (1),, x ( t), y ( t),,: F : R m R n.,m, n, u.,, Sigmoid. :,f Sigmoid,f ( x) = ^y k ( t) =

[9] R Ã : (1) x 0 R A(x 0 ) = 1; (2) α [0 1] Ã α = {x A(x) α} = [A α A α ]. A(x) Ã. R R. Ã 1 m x m α x m α > 0; α A(x) = 1 x m m x m +

27 :OPC 45 [4] (Automation Interface Standard), (Costom Interface Standard), OPC 2,,, VB Delphi OPC, OPC C++, OPC OPC OPC, [1] 1 OPC 1.1 OPC OPC(OLE f

132 包 装 工 程 2016 年 5 月 网 产 品 生 命 周 期 是 否 有 与 传 统 产 品 生 命 周 期 曲 线 相 关 的 类 似 趋 势 旨 在 抛 砖 引 玉, 引 起 大 家 对 相 关 问 题 的 重 视, 并 为 进 一 步 研 究 处 于 不 同 阶 段 的 互 联 网

11 25 stable state. These conclusions were basically consistent with the analysis results of the multi - stage landslide in loess area with the Monte

km km mm km m /s hpa 500 hpa E N 41 N 37 N 121

[1] Nielsen [2]. Richardson [3] Baldock [4] 0.22 mm 0.32 mm Richardson Zaki. [5-6] mm [7] 1 mm. [8] [9] 5 mm 50 mm [10] [11] [12] -- 40% 50%

标题

successful and it testified the validity of the designing and construction of the excavation engineering in soft soil. Key words subway tunnel

35期

CFDesign 2 1 CFDesign CFDesign CAD ~ r /min mm 1

<4D F736F F D20C8EDCDC1B5D8BBF9CDB2BBF9CAD4B2C9C6BDCCA8B5C4CACAD3C3D0D4B7D6CEF6>

Microsoft Word - 19王建华.doc

m m m ~ mm

Revit Revit Revit BIM BIM 7-9 3D 1 BIM BIM 6 Revit 0 4D 1 2 Revit Revit 2. 1 Revit Revit Revit Revit 2 2 Autodesk Revit Aut

31 17 www. watergasheat. com km 2 17 km 15 km hm % mm Fig. 1 Technical route of p

定稿

2 ( 自 然 科 学 版 ) 第 20 卷 波 ). 这 种 压 缩 波 空 气 必 然 有 一 部 分 要 绕 流 到 车 身 两 端 的 环 状 空 间 中, 形 成 与 列 车 运 行 方 向 相 反 的 空 气 流 动. 在 列 车 尾 部, 会 产 生 低 于 大 气 压 的 空 气 流

708 北 京 工 业 大 学 学 报 2011 年 以 往 的 试 验 结 果 进 行 对 比, 选 取 15D 20D 作 为 对 比 参 数, 试 件 参 数 见 表 1. Fig. 1 图 1 试 件 尺 寸 及 配 筋 图 ( mm) Geometry and reinforcement

H 2 SO ml ml 1. 0 ml C 4. 0 ml - 30 min 490 nm 0 ~ 100 μg /ml Zhao = VρN 100% 1 m V ml ρ g

[1] Liu Hongwei,2013, Study on Comprehensive Evaluation of Iron and Steel Enterprises Production System s Basic Capacities, International Asia Confere

Microsoft Word 聂雪梅.doc

\\Pc6\f\湖工大(自然科学)\自科

JOURNAL OF EARTHQUAKE ENGINEERING AND ENGINEERING VIBRATION Vol. 31 No. 6 Dec

% 8. 48% 3 80 Alcalase Novozymes Alcalase 2. 4 L Bacillus licheniformis 2. 4 AU /g 1. 2 Hitachi S-4700 JEOL JEM-1200EX Olympus Bu

cm /s c d 1 /40 1 /4 1 / / / /m /Hz /kn / kn m ~

LaDefense Arch Petronas Towers 2009 CCTV MOMA Newmark Hahn Liu 8 Heredia - Zavoni Barranco 9 Heredia - Zavoni Leyva

doc

(L d L q ). i d = 0,, [8].,,,,, / (maximum torque per ampere, MTPA) [9],. MTPA [9] [10],,,. MTPA, MTPA.,,,. MTPA,,,.,. 2 (Model of PMSM) d-q, :

蒋 哲等 基于 CVM 的有机农业认知程度和消费意愿研究 201 (P 0.000). On the basis of these results, enhancing publicity, pricing reasonably, focusing on high- income groups ar

Vol. 22 No. 4 JOURNAL OF HARBIN UNIVERSITY OF SCIENCE AND TECHNOLOGY Aug GPS,,, : km, 2. 51, , ; ; ; ; DOI: 10.

六 到 八 歲 兒 童, 設 計 並 發 展 一 套 以 van Hiele 幾 何 思 考 層 次 理 論 為 基 礎 的 悅 趣 化 學 習 數 位 教 材, 取 名 為 米 德 玩 形 狀, 同 時 探 討 低 年 級 學 童 在 使 用 本 數 位 教 材 之 後, 在 平 面 幾 何 的

<4D F736F F D D DBACEC0F25FD0A3B6D4B8E55F2DB6FED0A32D2D2DC8A5B5F4CDBCD6D0B5C4BBD8B3B5B7FBBAC52E646F63>

Dan Buettner / /

月32期

34 7 S R θ Z θ Z R A B C D PTP θ t 0 = θ 0 θ t 0 = 0 θ t 0 = 0 θ t = θ θ t = 0 θ t = 0 θ t V max θ t a max 3 θ t A θ t t 0 t / V max a max A = 3 4 S S

39 10 : PWM Fig.1 Flow ChartofCrowd MonitoringApproachBasedonPWM (pixel)d r (mm) 2(b) Q r A D r (1) : D r = H/ cosθ = H/ cos rε f (1) ε CCD Q r

E T 0 = γ 0 = 1 + R γ = nσ n nσ n ΔT 2 i - Σ n ΔT i T Pi - Σ n Σ n

中国媒体发展研究报告

标题

mm 5 1 Tab 1 Chemical composition of PSB830 finishing rolled rebars % C Si Mn P S V 0 38 ~ 1 50 ~ 0 80 ~ ~

m m 5m ~ PTFE Q45B 1 670MPa PMSAP MIDAS SAP2000 MIDAS MIDAS s 0.

THE APPLICATION OF ISOTOPE RATIO ANALYSIS BY INDUCTIVELY COUPLED PLASMA MASS SPECTROMETER A Dissertation Presented By Chaoyong YANG Supervisor: Prof.D

Vol. 15 No. 1 JOURNAL OF HARBIN UNIVERSITY OF SCIENCE AND TECHNOLOGY Feb O21 A

2011第1期第二部分

标题

Transcription:

2 0 1 3 年 8 月农业机械学报第 44 卷第 8 期 doi:10. 6041 / j. issn. 1000 鄄 1298. 2013. 08. 050 * 基于机器视觉的五坐标机床旋转轴误差检测方法 1 孙惠娟 2 蒋红海 1 殷国富 (1. 四川大学制造科学与工程学院, 成都 610065; 2. 昆明理工大学机电工程学院, 昆明 650093) 摘要 : 提出了一种利用机器视觉技术对五坐标机床旋转轴转角定位误差进行检测的方法 首先, 通过制作特定的 标志, 用 CCD 相机获取标志图像 ; 然后, 通过数字图像处理技术对所获得的图像进行分析处理 ; 最后, 根据标志在不 同位置处的相对转角偏差计算机床旋转轴的转角定位误差, 实现五坐标机床旋转轴转角定位误差的辨识和测量 同时, 将该方法与传统检测方法进行对比, 实验结果表明, 所提出的检测方法简单 高效, 可以实现机床旋转轴误差 的快速检测, 并为五坐标机床旋转轴误差的补偿提供了计算依据 关键词 : 五坐标机床旋转轴误差检测机器视觉图像处理 中图分类号 : TH161 + 郾 21; TG659 文献标识码 : A 文章编号 : 1000 鄄 1298(2013)08 鄄 0293 鄄 06 Error Measurement Method for Rotation Axes of Five 鄄 axis Machine Tool Based on Machine Vision Sun Huijuan 1 Jiang Honghai 2 Yin Goufu 1 (1. School of Manufacturing Science and Engineering, Sichuan University, Chengdu 610065, China 2. Faculty of Mechanical and Electrical Engineering, Kunming University of Science and Technology, Kunming 650093, China) Abstract: A measurement method was proposed for rotation axes of five 鄄 axis machine tool based on machine vision. Firstly, the vision technology was adopted to obtain sign images by using CCD camera. Then, the digital image processing techniques were used to analyze these images. At last, the positioning errors of rotary axes were identified by calculating the relative rotary errors at different positions of rotation axes of machine tool. Comparing to traditional method, the experimental results showed that the proposed method with a simply, high efficiency and application value provided a computation basis for the error compensation of a machine tool. Key words: Five 鄄 axis NC machine tool Rotary axis Errors measurement Machine vision Image processing 引言五坐标机床比三坐标机床增加了两个旋转轴, 由此机床加工的灵活性大大增强, 材料切除率和工件的表面质量也得到了极大提高 [1] 五坐标机床具有许多普通机床无法比拟的优点, 但是其加工精度却往往低于普通机床, 主要原因在于, 增加的 2 个旋转轴由于缺少精度标定及误差补偿的方法, 因此旋转轴的误差成为五坐标机床准静态误差和动态误 差的主要来源 对五坐标机床旋转轴进行精度标定及误差补偿是提高机床旋转精度的关键问题 [2] 传统三坐标机床各项误差元素的检测技术已经较为成熟, 对于五坐标机床旋转轴的检测技术尚无统一标准 目前国内外研究人员在机床旋转轴误差检测 [3 ~ 方面已作了大量研究 12] 上述检测方法存在以下缺点 : 采用的检测仪器主要有专用的激光干涉仪和球杆仪, 用激光干涉仪检测, 方便快捷, 但价格比较昂贵 ; 用价格低廉的球杆仪进行旋转轴联动检测, 收稿日期 : 2012 05 10 修回日期 : 2012 09 06 * 十二五冶国家科技重大专项资助项目 (2012ZX04011 031) 和四川省科技支撑计划资助项目 (2011GZ0065) 作者简介 : 孙惠娟, 博士生, 主要从事机床精度检测研究,E 鄄 mail: scdx2013@ 163. com

294 农业机械学报 2 0 1 3 年 检测路径的计算条件及计算过程复杂 ; 旋转轴与平动轴联动测量时, 由于在测试过程中加入了平动轴的几何和运动误差, 需对检测结果进行误差分离处理, 辨识过程计算较复杂 而利用正 12 面或 24 面多棱镜及自准直仪对机床旋转轴进行误差检测, 则需要制作特定的工装, 检测过程繁琐, 检测效率较低 此外, 以上所述方法主要是针对工作台回转式五坐标机床进行检测, 对于主轴头回转式和主轴头及工作台回转式五坐标机床回转轴的精度检测还很少见, 因此目前的检测方法存在局限性 本文提出一种采用机器视觉技术非接触测量检测机床旋转轴误差的检测方法, 利用 CCD 相机获取机床旋转轴在不同位置处的图像, 通过数字图像处理技术, 计算分析机床旋转轴转角定位误差, 最后, 将该方法与传统检测方法进行对比来验证该方法的有效性 1 基本理论 1 郾 1 机器视觉技术机器视觉检测系统是采用 CCD 相机将被检测的目标转换成图像信号, 传送给专用的图像处理系统, 根据像素分布 亮度和颜色等信息, 转换成数字化信号, 图像处理系统对这些信号进行各种运算来抽取目标的特征, 如面积 数量 位置 长度, 再根据预设的允许度和其他条件输出结果, 包括尺寸 角度 个数 合格 / 不合格 有 / 无等, 实现自动识别功能 机器视觉检测系统主要用于一些不适合于人工作业的危险工作环境或人工视觉难以满足要求的场合及一些大批量的生产 检测过程中 由于机器视觉系统可以快速获取大量信息, 而且易于自动处理, 也易于同设计信息以及加工控制信息集成, 因此, 机器视觉系统在现代自动化生产和农业生产等领域得 [13 ~ 到了广泛应用 19] 1 郾 2 图像处理方法 1 郾 2 郾 1 边缘提取算法图像的边缘是图像最基本的特征之一, 边缘处的图像集中了图像的大部分信息, 对于整个图像场景的识别与理解非常重要, 同时也是图像分割所依赖的重要特征 [20] 边缘检测的效果将直接影响到图像的分析 识别和理解, 因此在对图像进行处理时常需对其边缘信息进行提取 最常用的边缘检测算法有 Roberts 算子 Prewitt 算子 Sobel 算子 Laplacian 算子和 LOG 算子等, 其中前 3 个属于一阶微分算子, 后 2 个属于三阶微分算子 1 郾 2 郾 2 最小二乘椭圆拟合定中心最小二乘法椭圆拟合是较常用的椭圆拟合方 法 [21 ~ 22], 其基本原理是先假设椭圆参数, 得到每个 待拟合点到该椭圆的距离之和, 并求出使这个和最 小的参数 假设一般形式的椭圆方程为 Ax 2 + Bxy + Cy 2 + Dx + Ey + F = 0 (1) 求目标函数 f( A,B,C,D,E,F) = n 移 (Ax 2 i + Bx i y i + Cy 2 i + Dx i + Ey i + F) 2 i = 1 的最小值来确定式 (1) 中的各个系数, 由此得出所 拟合的椭圆 2 机器视觉检测方法 采用机器视觉进行机床旋转轴定位精度检测的 方法分为图像获取 图像处理及分析 图像获取包 括检测标志制作与固定 相机和光源安装等 ; 图像处 理包括对图像边缘提取 圆心位置确定及图像在不 同旋转角度下相对角度变化计算等 具体检测及处 理流程如图 1 所示 Fig. 1 图 1 2 郾 1 标志设计 机器视觉检测方法流程图 Flow chart of machine vision measurement 由于机床旋转轴与周围部件色差太小, 如果直 接利用相机对旋转轴上特定部位进行拍摄, 则获取 图像标志点与图像背景对比度不高, 不利于图像后 期处理 因此在利用机器视觉进行检测时, 需要根 据机床旋转轴的特点制作具有代表性的标志固定于 机床旋转轴上, 随着机床旋转轴的旋转, 圆形标志在 旋转平面发生转动, 利用同一个标志在不同图像中 旋转的相对角度差来计算机床旋转轴转动的转角误 差 2 郾 2 图像获取和图像处理 机床旋转轴旋转一定的角度, 用相机拍摄机床 旋转轴上的标志在不同位置处的图像, 并对图像进 行标记以利于后续的处理和分析 2 郾 2 郾 1 图像边缘提取 在透视相机模型下, 实际获取得到的圆形图像 是椭圆, 因此需要对椭圆的边界进行检测, 经过边缘 检测和阈值处理过程后, 得到椭圆边缘的二值图像 本文利用 Canny 算子对所获得的图像进行图像边缘 像素位置信息提取 2 郾 2 郾 2 最小二乘椭圆拟合中心 提取图像边缘点后, 为了精确地定出圆心的位 置, 需对图像的边缘点进行拟合 本文根据最小二

第 8 期孙惠娟等 : 基于机器视觉的五坐标机床旋转轴误差检测方法 295 乘原理 ( 残差平方和最小 ), 用椭圆来拟合圆轮廓并 找出椭圆中心 2 郾 2 郾 3 计算偏转角度误差 经椭圆拟合定出各同心圆的中心点后, 利用图 像中标志在旋转前后的变化量来计算两次旋转的角 度差, 即 驻兹忆 = 兹忆 i + 1 - 兹忆 i 并与机床理论上两次转动的角度差 驻兹 = 兹 i + 1 - 兹 i 进行对比, 计算机床的转角误差为 E 兹 = 驻兹 - 驻兹忆 3 机床旋转轴转角定位误差检测 3 郾 1 实验设备及仪器 将本文提出的基于机器视觉技术检测五坐标机 床旋转轴转角定位误差的方法应用于四川省某机床 厂生产的五坐标加工中心上 ( 图 2), 以该机床 B 轴 ( 转动范围为 0 毅 ~ 110 毅 ) 转角定位精度检测为例, 按 照第 2 节提出的方法进行检测和分析 Fig. 2 图 2 五坐标机床结构 Structure of five 鄄 axis machine tool 机器视觉检测系统由圆形标志 光源和相机组 成 检测系统的硬件部分如图 3 所示 其中数码相 机 (1 000 万像素 ) 所制作的 ( 部分 ) 圆形标志及标 志尺寸 ( 单位 :mm) 如图 4 所示 光源可以根据现 场条件进行选择, 本实验中由于现场光照条件比较 好, 因此未利用光源 Fig. 4 图 4 3 郾 2 图像获取过程 检测所用的标志 Signs used for measurement (1) 将所制作的标志 ( 图 4) 固定于机床旋转轴 上, 尽量保持标志面平整 (2) 安装好相机和光源, 调整相机位置与焦距, 保证标志在相机的成像区域内 (3) 机床旋转轴顺时针和逆时针转动不同的角 度, 当机床摆头转动到规定的角度时, 摆头停止转 动, 相机拍摄此时机床旋转轴上的标志图片, 并将图 片存于计算机内用于图片的后续处理和分析 本文 在进行实验时采集了机床旋转轴从 0 毅转动到 90 毅的 图像, 每隔 10 毅拍摄一次图片, 共采集到 10 张图片 3 郾 3 图像处理过程 3 郾 3 郾 1 标志图像边缘提取 仅对机床 B 轴在 0 毅和 10 毅位置的图像进行分析 说明, 其他角度的图像处理方法与此相同 图 5 和 图 6 分别为实验对象机床 B 轴在 0 毅和 10 毅位置所拍 摄的现场图像, 在对图像进行处理前已利用图 7 所 示标定图像对相机进行标定 图 5 和图 6 中所有的 同心圆都应进行计算和处理, 为了更清晰地显示图 中的内容, 本文及后续的图像处理介绍均以图 5 和 图 6 中大圆圈内的 4 组同心圆为例进行说明 Fig. 3 图 3 机器视觉检测系统 Measurement system of machine vision 图 5 0 毅时标志 Fig. 5 Sign in 0 毅利用 Matlab 自带的边缘提取函数对获得的图 像 ( 图 5 和图 6) 进行边缘信息提取, 提取边缘后的

296 农业机械学报 2 0 1 3 年 图 6 10 毅时标志 Fig. 6 Sign in 10 毅 图 9 提取图形边缘 (10 毅 ) Fig. 9 Boundary information extraction (10 毅 ) Fig. 7 图 7 标定图像 Calibration image 图像如图 8 和图 9 所示 从图中可以看出从原始图 像中提取的边缘点组成的图形并不是光滑的圆, 而 是有一些断点和散点, 因此为了精确地求出圆心的 位置需对这些边缘点进行椭圆拟合 图 10 椭圆拟合 (0 毅 ) Fig. 10 Ellipse fitting (0 毅 ) 图 8 提取图形边缘 (0 毅 ) Fig. 8 Boundary information extraction (0 毅 ) 3 郾 3 郾 2 椭圆拟合求圆心位置 根据最小二乘椭圆拟合算法对提取边缘后的图 像 ( 图 8 和图 9) 进行椭圆拟合, 拟合后图像如图 10 和图 11 所示 图中线条为拟合后的椭圆, 从图中也 可以看到一些白色的散点, 这些散点是从原始图像 提取出的边缘点中未被拟合的点 由拟合后的椭圆分别求出标志在 0 毅和 10 毅处各 图 11 椭圆拟合 (10 毅 ) Fig. 11 Ellipse fitting (10 毅 ) 组同心圆的圆心 R 1 R 2 R 3 R 4 和 R 忆 1 R 忆 2 R 忆 3 R 忆 4 的位 置 标志在 0 毅和 10 毅处各组圆心点的变化如图 12 所示 3 郾 3 郾 3 相对转角误差计算 由图 12 可知, 理论上连接各圆心的 6 条线 R 1 R 2 R 1 R 3 R 1 R 4 R 2 R 3 R 2 R 4 R 3 R 4 在 0 毅和 10 毅位置 处的相对角度差为 10 毅, 由于装配 震动 热变形等 因素影响, 机床旋转轴在实际转动过程中与理论值

第 8 期孙惠娟等 : 基于机器视觉的五坐标机床旋转轴误差检测方法 297 Fig. 12 图 12 图像处理结果 Results of image processing 之间存在偏差, 因此机床旋转轴在实际转动后所采 集的标志图像与理论上旋转后应转动的位置之间必 然存在偏差, 该偏差主要由机床旋转轴的转角定位 误差引起 根据各条直线旋转前后斜率的变化计算 旋转的角度, 并与机床旋转轴理论转动的角度相比 较得出转角定位误差, 取各条线旋转前后的误差平 均值为机床旋转轴的旋转误差, 其计算结果如图 13 所示 3 郾 4 检测结果分析 为了验证本文所提检测方法的有效性, 将分别 利用机器视觉方法和多面棱体与自准直仪相结合的 方法对机床旋转轴 B 轴进行 0 毅到 90 毅转角定位精度 检测, 两种方法的检测结果对比如图 13 所示 由图 13 可知, 两种方法的检测结果比较接近, Fig. 13 图 13 检测结果 Measurement result 其残差最小值为 4 郾 2 义, 残差平均值为 4 郾 86 义 与利用 多面棱体方法相比, 利用本文提出的旋转轴误差检 测方法检测原理与检测过程更简单 检测效率更高, 且对检测设备及检测条件的要求不高 4 结论 (1) 利用机器视觉技术对机床旋转轴检测的原 理及方法简单 高效 与传统方法相比, 仅需制作特 定的标志固定于要检测的机床旋转轴上, 利用相机 对不同位置处的标志进行拍照即可实现非接触测 量, 且检测结果与传统检测方法比较接近 (2) 易于实现模块化集成 所提检测方法的关 键是图像的处理与分析, 对于获取的图像可通过编 制相关图像处理程序并进行分析处理即可获得机床 转角的定位误差, 无需进行人工记录与分析计算 参考文献 1 Hsu Y Y, Wang S S. A new compensation method for geometry errors of five 鄄 axis machine tools [ J]. International Journal of Machine Tools & Manufacture, 2007, 47(2):352 ~ 360. 2 廖平兰. 机床加工过程综合误差实时补偿技术 [J]. 机械工程学报,1992,28(2):65 ~ 68. Liao Pinglan. The integrated error real 鄄 time compensation of machining process[ J]. Chinese Journal of Mechanical Engineering, 1992,28(2):65 ~ 68. ( in Chinese) 3 Lei W T, Paung I M, Chi Chen. Total ballbar dynamic tests for five 鄄 axis CNC machine tools[j]. International Journal of Machine Tools & Manufacture, 2009,49(6):488 ~ 499. 4 Lei W T, Sung M P, Liu W L, et al. Double ballbar test for the rotary axes of five 鄄 axis CNC machine tools[ J]. International Journal of Machine Tools & Manufacture,2007,47(2):273 ~ 285. 5 张葳, 王娟. 五轴联动数控机床旋转轴几何误差测量与分离方法 [J]. 机电工程技术,2008,37(9):16 ~ 19. Zhang Wei, Wang Juan. A novel kind of geometric errors measurement and identification methods of 5 鄄 axis machine tool rotary axes [ J]. Mechanical & Electrical Engineering Technology, 2008,37(9):16 ~ 19. ( in Chinese) 6 周玉清, 陶涛, 梅雪松, 等. 旋转轴与平移轴联动误差的快速测量及溯源 [J]. 西安交通大学学报,2010,44(5): 80 ~ 84. Zhou Yuqing, Tao Tao, Mei Xuesong, et al. Rapid error measurement and trace of simultaneous rotary 鄄 translational axes collaborative motion[ J]. Journal of Xian Jiaotong University, 2010,44(5): 80 ~ 84. ( in Chinese) 7 付璇, 田怀文, 朱绍维. 五轴数控机床旋转轴几何误差测量与建模 [J]. 机械设计与制造,2011(2):157 ~ 159. Fu Xuan, Tian Huaiwei, Zhu Shaowei. Geometric error measurement and modeling for five 鄄 axis CNC machine tool with rotary axes [ J]. Machinery Design and Manufacture,2011(2): 157 ~ 159. ( in Chinese) 8 张宏韬. 双转台五坐标数控机床误差的动态实时补偿研究 [D]. 上海 : 上海交通大学, 2011. Zhang Hongtao. Research on dynamic and real 鄄 time error compensation for two turntable five 鄄 axis NC machine tool[d]. Shanghai: Shanghai Jiaotong University, 2011. ( in Chinese) 9 张大卫, 商鹏, 田延岭, 等. 五轴数控机床转动轴误差元素的球杆仪检测方法 [ J]. 中国机械工程,2008,19(22):2 737 ~ 2 741.

298 农业机械学报 2 0 1 3 年 Zhang Dawei, Shang Peng, Tian Yanling, et al. A DBB 鄄 based alignment error measurement method for rotary axis of 5 鄄 axis CNC machine tool[ J]. China Mechanical Engineering,2008,19(22): 2 737 ~ 2 741. ( in Chinese) 10 刘飞. 五轴数控机床回转轴的误差检测技术研究 [J]. 机械工程与自动化,2009(4):133 ~ 135. Liu Fei. Error calculation for rotation axes of five 鄄 axis machine tool[j]. Mechanical Engineering and Automation, 2009(4):133 ~ 135. (in Chinese) 11 高秀峰, 刘春时, 李焱, 等. 基于激光干涉仪的 A / C 轴双摆角铣头定位误差检测与辨识 [ J]. 机械设计与制造,2011(12): 212 ~ 214. Gao Xiufeng, Liu Chunshi, Li Yan, et al. Positioning errors detection and identification for A / C axes bi 鄄 rotary milling head based on laser interferometer[ J]. Machinery Design and Manufacture,2011(12):212 ~ 214. ( in Chinese) 12 崔勇, 徐岩. 精密转台角位移精度的测量方法研究 [J]. 计量与测试技术,2012,39(2):1 ~ 2. Cui Yong, Xu Yan. Research of angular displacement accuracy measuring method for precise turntable [ J ]. Metrology & Measurement Technique, 2012,39(2):1 ~ 2. ( in Chinese) 13 伍济钢, 宾鸿赞. 薄片零件尺寸机器视觉检测系统研究与开发 [J]. 机床与液压,2010,38(17):86 ~ 88. Wu Jigang, Bin Hongzan. Research and development of dimensional inspection system for thin sheet part based on machine vision [ J]. Machine tool and Hydraulics,2010,38(17):86 ~ 88. ( in Chinese) 14 何博侠, 张志胜, 徐孙浩, 等. 大尺寸机械零件的机器视觉高精度测量方法 [J]. 中国机械工程,2009,20(1):5 ~ 10. He Boxia, Zhang Zhisheng, Xu Sunhao, et al. Research on high 鄄 precision machine vision measurement method for large scale parts[ J]. China Mechanical Engineering,2009,20(1):5 ~ 10. ( in Chinese) 15 贺雅琴, 李琳. 基于机器视觉的零件尺寸检测的研究 [J], 煤矿机械,2010,31(6):102 ~ 104. He Yaqin, Li Lin. Research for part size inspection base on machine vision[ J]. Coal Mine Machinery, 2010,31(6):102 ~ 104. ( in Chinese) 16 司永胜, 乔军, 刘刚, 等. 基于机器视觉的苹果识别和形状特征提取 [J]. 农业机械学报,2009,40(8):161 ~ 167. Si Yongsheng, Qiao Jun, Liu Gang, et al. Recognition and shape features extraction of apples based on machine vision[ J]. Transactions of the Chinese Society for Agricultural Machinery,2009,40(8): 161 ~ 167. ( in Chinese) 17 陈兵旗, 孙旭东, 韩旭, 等. 基于机器视觉的水稻种子精选技术 [J]. 农业机械学报,2010,41(7):168 ~ 174. Chen Bingqi, Sun Xudong, Han Xu, et al. Rice seeds selection based on machine vision[ J]. Transactions of the Chinese Society for Agricultural Machinery,2010,41(7):168 ~ 174. ( in Chinese) 18 王海青, 姬长英, 顾宝兴, 等. 基于机器视觉和支持向量机的温室黄瓜识别 [J]. 农业机械学报,2012,43(3):163 ~ 168. Wang Haiqing, Ji Changying, Gu Baoxing, et al. In 鄄 greenhouse cucumber recognition based on machine vision and least squares support vector machine [ J]. Transactions of the Chinese Society for Agricultural Machinery,2012, 43 ( 3):163 ~ 168. ( in Chinese) 19 Li Daoliang, Yang Wenzhu, Wang Sile. Classification of foreign fibers in cotton lint using machine vision and multi 鄄 class support vector machine[ J]. Computers and Electronics in Agriculture, 2010,74(2): 274 ~ 279. 20 黄剑玲, 邹辉. 一种精确的自适应图像边缘提取方法 [J]. 计算机工程与科学,2009,31(9):53 ~ 55. Huang Jianling, Zou Hui. An exact method for adaptive image edge detection [ J]. Computer Engineering & Science, 2009, 31(9):53 ~ 55. (in Chinese) 21 Gander W, Golub G H, Strebel R, et al. Least 鄄 squares fitting of circles and ellipses[ J]. BIT Numerical Mathematics, 1994, 34(4):558 ~ 578. 22 Rosin P L. Ellipse fitting by accumulating five 鄄 point fits[j]. Pattern Recognition Letters, 1993,14(8):661 ~ 669.