CN102679926B - 用于多点阵列式柔性工装的基于包围盒的薄壁曲面类零件定位方法 - Google Patents
用于多点阵列式柔性工装的基于包围盒的薄壁曲面类零件定位方法 Download PDFInfo
- Publication number
- CN102679926B CN102679926B CN201210167738.8A CN201210167738A CN102679926B CN 102679926 B CN102679926 B CN 102679926B CN 201210167738 A CN201210167738 A CN 201210167738A CN 102679926 B CN102679926 B CN 102679926B
- Authority
- CN
- China
- Prior art keywords
- bounding box
- thin
- location
- curved
- curve part
- 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
- Processing Or Creating Images (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及一种用于多点阵列式柔性工装的基于包围盒的薄壁曲面类零件定位方法。薄壁曲面件在柔性工装中的定位过程是一种将曲面定位到工装工作区域内的一个过程,其实质为坐标变换。国内外的多点阵列式柔性工装的工作区域通常可用一个正六面体来表述,则薄壁曲面类零件的定位问题可转换为求取包络薄壁曲面件的正六面体包围盒,根据包围盒和柔性工装工作区域的相似性来确定坐标变换矩阵中涉及到的未知量,即求取薄壁曲面件的包围盒问题。为了提高支撑定位的稳定性,还应使曲面零件与柔性工装中尽可能多的支撑杆相接触,即曲面在工装平面内的投影面积最大。定位过程主要涉及3个步骤:初定位、精定位和平移。然后根据坐标变化矩阵计算出曲面类零件在工装坐标系下的参数化曲面方程,最后根据柔性工装中支撑杆的位置计算出支撑杆的高度。
Description
技术领域
该发明涉及一种曲面类零件的定位方法,具体涉及到一种用于多点阵列式柔性工装的基于包围盒的薄壁曲面类零件定位方法。
背景技术
工业中一些薄壁曲面类零件如飞机薄壁曲面件等的典型制造工艺是先成型后加工,由于这类零件的刚度小,型面复杂,自身在重力下容易变形,给零件成形后的曲面测量、轮廓划线、裁边、钻孔等加工操作带来了很大的困难。薄壁曲面件的传统加工定位主要依靠具有产品形状特征的模线样板和标准工艺装备,即工装表面需要具有和薄壁曲面件一致的曲面形状。传统的定位方法协调环节多,对零件的准确度影响大,且工装笨重,一套工装只针对一种零件,适应性差,成本高,难以满足各行业日益提高的多品种、小批量的需求。目前国内外对于薄壁曲面件的加工装配主要通过多点阵列式柔性工装来实现零件的定位和支撑,但关于零件在柔性工装中的定位算法较少。
目前有关曲面定位的问题主要集中在将一个曲面(往往是测量数据点)和理想曲面进行匹配的问题。常用的方法是通过特征曲面、曲率匹配来实现曲面匹配的初定位及精调整,计算过程较复杂,需要较多次的迭代计算,这类方法主要用于曲面测量、非定位/欠定位加工等方面,在曲面匹配过程没有考虑薄壁曲面件和多点阵列式柔性工装的结构特点,无法解决薄壁曲面件在柔性工装中的定位问题。
发明内容
针对现有方法的不足之处,提出了一种用于多点阵列式柔性工装的基于包围盒的薄壁曲面类零件定位方法。
本发明解决其技术问题所采用的技术方案是:一种用于多点阵列式柔性工装的基于包围盒的薄壁曲面类零件定位方法,其特征在于:根据曲面零件的包围盒和多点阵列式柔性工装工作区域的相似性来确定坐标变换矩阵中涉及到的未知量,具体步骤如下:设多点阵列式柔性工装坐标系为 ,薄壁曲面件的坐标系为,则薄壁曲面件在柔性工装坐标系中定位过程可由转换矩阵T来描述,具体可写成为:
式中,Ti(i=1,2)和T3分别表示为:
式中,字母c表示cos函数,字母s表示sin函数,T1表示预定位旋转矩阵,T2表示精定位旋转矩阵,T3表示平移矩阵,具体涉及的矩阵求解方法和定位过程按照以下步骤进行:
(1)、构建薄壁曲面件的参数化方程;
(2)、用协方差方法求解方向包围盒;
(3)、基于包围盒的预定位;
(4)、基于包围盒迭代的精定位;
(5)包围盒的平移;
(6)、柔性工装中支撑杆的高度求解。
薄壁曲面件在柔性工装中的定位过程是一种将曲面定位到工装工作区域内的一个过程,其实质为坐标变换。国内外的多点阵列式柔性工装的工作区域通常可用一个正六面体来表述,则薄壁曲面件的定位问题可转换为求取包络薄壁曲面件的正六面体包围盒,根据包围盒和柔性工装工作区域的相似性来确定坐标变换矩阵中涉及到的6个参数值,即求取薄壁曲面件的包围盒问题。而为了提高支撑定位的稳定性,应使曲面零件与柔性工装中尽可能多的支撑杆相接触,即曲面在工装平面内的投影面积最大。定位过程主要涉及3个步骤:初定位、精定位和平移,最后确定出曲面在工装坐标系下的参数化曲面方程,然后根据柔性工装中支撑杆的位置确定支撑杆的上升高度。
本发明的有益效果是方法简单,易于理解,在保证支撑稳定性的前提下,减少支撑杆数量。
附图说明
图1为基于多点阵列式柔性工装的薄壁曲面类零件定位示意图;
图2为算法的流程图;
图3为曲面的OBB(Oriented Bounding Box)图形;
图中,1.曲面类零件,2.多点阵列式柔性工装中的支撑杆,OBB的中文意思就是方向包围盒。
具体实施方式
式中,Ti(i=1,2)和T3分别表示为:
式中,字母c表示cos函数,字母s表示sin函数,T1表示预定位旋转矩阵,T2表示精定位旋转矩阵,T3表示平移矩阵。具体涉及的矩阵求解方法和定位过程按照以下步骤进行:
1、构建薄壁曲面件的参数化方程
薄壁曲面件的外形采用B样条曲面进行描述,其表达式为:
式中:i=1,2,…,m;j=1,2,…,n;gi,j 为B样条曲面片的控制顶点,m,n分别为B样条曲面片的行数和列数。当m,n取不同的值,参变量u,v取为均匀基或非均匀基时,曲面具有很大的灵活性,用于描述曲面零件可获得良好的效果。
2、 用协方差方法求解方向包围盒
求解薄壁曲面件的曲面上的数据点的包围盒,根据曲面的参数化方程式,对参数u,v进行取值,设Pi,j为曲面上对应参数ui,vj的一个点,则曲面上点集的协方差矩阵为:
协方差矩阵C的3个特征向量是正交的,归一化可作为一个基底R用于确定OBB的方向。R的表示形式为:
式中n1,n2,n3为协方差矩阵的归一化特征向量,表示OBB的三个方向。如图3所示为某曲面的OBB包围盒,n1、n2、n3分别为包围盒的3个方向,P点为包围盒的一个顶点。
3、 基于包围盒的预定位
对包围盒进行旋转,使其n1、n2、n3的三个方向分别和图1所示的工装中Z、X、Y轴相对应,求取变换矩阵中的3个角度参数α1、β1、γ1。其表达式分别为:
其中,n2’为矢量n2经过α1、β1旋转变换后的矢量。找出各方向的极值点(xmin,ymin,zmin)(xmax,ymax,zmax),得到其OBB。各角度值、、,为初定位的旋转坐标变换矩阵T1中的角度值。
4、基于包围盒迭代的精定位
经过上述的包围盒预定位后,OBB可以满足曲面的预定位要求,但不是最优的结果。可通过迭代算法进行精定位,迭代过程的实质是在工装平面内最大投影面的面积最小化问题,具体步骤为:
1)将方向包围盒绕Y轴旋转变换,使变换后的薄壁曲面件的方向包围盒在XY平面内的投影面积最大,即:
2)将方向包围盒绕X轴旋转变换,使变换后的薄壁曲面件的方向包围盒在XY平面内的投影面积最大,即:
3)将方向包围盒绕Z轴旋转变换,使变换后的薄壁曲面件的方向包围盒在XY平面内的投影面积最小,即:
经过步骤1)和2)后包围盒在XY平面内的投影面积最大,即曲面类零件在工装系统中投影面积最大,此时与曲面类零件相接触的支撑杆的数量最多。但对于阵列式的柔性工装而言,此时参与工作的支撑杆的行列数量较多,不是一个优化的结果,而经过步骤3)处理后可以有效减少工装中参与工作的支撑杆行列数量,在支撑杆行列数量相同的情况下,能够满足更大尺寸的曲面类零件的需求。在上述步骤中包围盒绕各轴旋转获得投影面积最小值时的角度值、、,为精定位旋转坐标变换矩阵T2中的角度值。
5、包围盒的平移
经过预定位和精定位后,求取包围盒顶点P的坐标,根据顶点P的坐标和工装坐标系之间的位置关系,确定平移变量l,m,n。设包围盒顶点P的坐标值为(xP、yP、zP),使点P和工装上的参考点D(工装上的某个指定点)重合,可确定出变换矩阵中的3个平移参数,其表达式为:
式中:(xD、yD、zD)为点D的位置坐标。
B样条曲面的几何信息完全包含在控制顶点矩阵G中,所以薄壁曲面件在工装坐标系中的方程可通过零件数模的控制顶点gi,j来实现,即工装内的薄壁曲面件曲面方程可表示为:
6、柔性工装中支撑杆的高度求解
在多点阵列式柔性工装中各支撑杆在x,y方向的位置是固定的,z轴方向根据曲面的需求进行调整。在已知参数曲面方程F f(u,v)的情况下,求解非线性方程组,可确定出工装中各支撑杆的最终位置坐标:
式中,xi,yi为各支撑杆在工装上x,y方向的对应位置坐标,m和n分别为x,y方向支撑杆的数目。
根据求解的结果u,v,可以得到支撑杆的上升高度:
Claims (1)
1.一种用于多点阵列式柔性工装的基于包围盒的薄壁曲面类零件定位方法,其特征在于:根据曲面零件的包围盒和多点阵列式柔性工装工作区域的相似性来确定坐标变换矩阵中涉及到的未知量,具体步骤如下:设多点阵列式柔性工装坐标系为 ,薄壁曲面件的坐标系为,则薄壁曲面件在柔性工装坐标系中定位过程可由转换矩阵T来描述,具体可写成为:
式中,Ti(i=1,2)和T3分别表示为:
式中,字母c表示cos函数,字母s表示sin函数,T1表示预定位旋转矩阵,T2表示精定位旋转矩阵,T3表示平移矩阵,具体涉及的矩阵求解方法和定位过程按照以下步骤进行:
(1)、构建薄壁曲面件的参数化方程;
(2)、用协方差方法求解方向包围盒;
(3)、基于包围盒的预定位;
(4)、基于包围盒迭代的精定位;
(5)包围盒的平移;
(6)、柔性工装中支撑杆的高度求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210167738.8A CN102679926B (zh) | 2012-05-28 | 2012-05-28 | 用于多点阵列式柔性工装的基于包围盒的薄壁曲面类零件定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210167738.8A CN102679926B (zh) | 2012-05-28 | 2012-05-28 | 用于多点阵列式柔性工装的基于包围盒的薄壁曲面类零件定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102679926A CN102679926A (zh) | 2012-09-19 |
CN102679926B true CN102679926B (zh) | 2014-05-07 |
Family
ID=46812217
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210167738.8A Expired - Fee Related CN102679926B (zh) | 2012-05-28 | 2012-05-28 | 用于多点阵列式柔性工装的基于包围盒的薄壁曲面类零件定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102679926B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103093029A (zh) * | 2012-12-19 | 2013-05-08 | 江苏申模数字化制造技术有限公司 | 用于点阵式柔性工装的离散支承杆预调高度计算方法 |
CN104393717B (zh) * | 2014-12-15 | 2016-08-10 | 上海电气集团上海电机厂有限公司 | 一种多点成形柔性模具的调形方法 |
CN106475790B (zh) * | 2016-11-01 | 2018-06-15 | 河海大学常州校区 | 大曲率薄壁件多点支撑定位方法、柔性工装及设定方法 |
TWI707740B (zh) * | 2019-11-26 | 2020-10-21 | 財團法人工業技術研究院 | 可調式工件支撐系統及方法 |
TWI716252B (zh) * | 2020-01-03 | 2021-01-11 | 中華精測科技股份有限公司 | 旋轉及定位複數測物之方法 |
CN113857769B (zh) * | 2021-09-22 | 2023-11-03 | 国营芜湖机械厂 | 基于零件外形自适应调整的复杂曲面零件加工方法及装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE29616184U1 (de) * | 1996-09-17 | 1996-11-07 | Kaiser, Stephan, 51145 Köln | 3D-Baukastensystem |
CN1300545C (zh) * | 2005-06-16 | 2007-02-14 | 上海交通大学 | 大型复杂曲面冲压薄壁件的柔性检具 |
CN102254343B (zh) * | 2011-07-01 | 2013-01-30 | 浙江理工大学 | 一种基于凸包与obb的三维网格模型骨架提取方法 |
-
2012
- 2012-05-28 CN CN201210167738.8A patent/CN102679926B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102679926A (zh) | 2012-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102679926B (zh) | 用于多点阵列式柔性工装的基于包围盒的薄壁曲面类零件定位方法 | |
CN110516388B (zh) | 基于调和映射的曲面离散点云模型环切刀轨生成方法 | |
CN106568381B (zh) | 一种基于标准平面的线激光测量系统标定方法 | |
CN112396690B (zh) | 基于改进型向心参数化法的曲面高精重构方法 | |
CN105354396A (zh) | 一种不同软件协同建模的几何比对方法 | |
CN109767416B (zh) | 机械设备的定位系统及方法 | |
CN106875439B (zh) | 基于三维点云模型的单晶硅棒外形尺寸测量方法 | |
CN107091627B (zh) | 一种多表面系统的综合测量与评估方法 | |
CN103077546B (zh) | 二维图形的三维透视变换方法 | |
CN112508895B (zh) | 一种基于曲面配准的螺旋桨叶片质量评估方法 | |
CN107238374B (zh) | 一种不规则平面零件的分类和识别定位方法 | |
CN106199519A (zh) | 一种超短基线五基元立体空间基阵及其水声定位方法 | |
CN111400830B (zh) | 一种三维毛坯工件的加工校准方法及装置 | |
CN104570928A (zh) | 基于共形参数化的网格曲面上数控加工轨迹规划方法 | |
CN112685858B (zh) | 一种基于三维扫描辅助定位的锻件校调方法 | |
CN113251919A (zh) | 一种基于坐标变换的圆柱体点云数据拟合方法及装置 | |
CN111562563A (zh) | 激光雷达转台标定方法、装置和计算机可读存储介质 | |
CN103292769B (zh) | 一种基于最小区域的平面倾斜度误差检测方法 | |
CN110340738B (zh) | 一种基于pca的机器人拉丝高铁白车身工件精确标定方法 | |
CN103191958A (zh) | 一种板料成形回弹的光学检测方法 | |
CN110363801B (zh) | 工件实物与工件三维cad模型的对应点匹配方法 | |
CN115861547A (zh) | 一种基于投影的模型表面样条线生成方法 | |
CN102831281A (zh) | 在gocad软件中建立二维工作面的辅助建模方法 | |
Tan et al. | An unconstrained approach to blank localization with allowance assurance for machining complex parts | |
CN110060306A (zh) | 一种三维空间异型曲面旋转标定实体模型及3d相机旋转标定方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140507 Termination date: 20180528 |
|
CF01 | Termination of patent right due to non-payment of annual fee |