基于proe的变形曲面网格优化及几何重建技术研究

基于Pro/E的变形曲面网格优化及几何重建技术研究,答辩人: 潘 永 信 导 师: 张 树 有 教授 专 业: 机械设计及理论,浙江大学硕士学位论文答辩,Research on the Mesh Optimization and Geometry Reconstruction of Deformation Surface based on Pro/E,简介,随着计算机在工程设计领域应用的进一步发展，以虚拟交互技术和计算机辅助设计技术为基础的虚拟CAD日益受到重视。曲面重建技术在CAD工程与可视化领域有着广泛的应用，本文基于Pro/E系统三维造型信息，提取三角剖分网格；面向曲面变形过程，研究网格模型建模、网格模型优化及几何重建中的相关技术方法。,主要内容,绪论曲面网格层次模型及数据转换 曲面网格优化 拓扑网格的几何重建 曲面重建模块的实现及应用 总结与展望,绪论,课题背景虚拟设计研究现状曲面重建研究现状 主要的研究内容,课题背景,CAD技术成为现代工业设计的核心技术之一传统工业造型设计软件难以对三维形体曲面进行直观、方便地设计和编辑。外型创新设计只需将已有的设计曲面作局部的修改。计算机技术的成熟，产品设计向虚拟化发展基于虚拟现实的VR-CAD技术成为研究热点。虚拟环境下的曲面变形设计系统基于外部CAD模型的虚拟变形设计系统 。变形曲面网格的优化及重建。,虚拟设计技术研究现状,虚拟造型虚拟实体造型虚拟曲面造型自由曲面变形设计虚拟特征造型,虚拟样机虚拟样机可视化虚拟装配虚拟环境下人机分析等,虚拟曲面设计加拿大Alberta大学的THRED自由曲面设计系统美国爱荷华州立大学的C2系统 美国SGI公司 浙江大学CAD&CG国家重点实验室,曲面重建研究现状,拓扑重建高度点元数据截面线元数据散乱点元数据Hoppe、Amenta、Guo网格优化基于顶点聚合的算法Rossignac、 Luebke基于区域合并的算法 Kalvin、Reddy基于几何元素删除的算法 Schroeder、Turk能量函数优化算法Hoppe几何重建单值网格任意拓扑网格构造网格细化曲面Chaikin、Catmull、Clark构造分段连续拼接曲面Bezier、Farin推广B样条曲面Grimm、Hughes,主要研究内容,基于外部CAD的变形曲面模型信息组织曲面网格的简化及优化曲面网格的几何重建,主要研究内容,基于外部CAD的变形曲面模型信息组织曲面网格的简化及优化曲面网格的几何重建,主要研究内容,基于外部CAD的变形曲面模型信息组织曲面网格的简化及优化曲面网格的几何重建,第二部分,曲面网格层次模型及数据转换,零件信息的层次网格模型,层次模型定义零件层特征层 几何层 显示层模型数据结构,式中： 数据对象标志ID； VP状态属性和改变属性的操作集； VP的消息处理接口； VP的底层数据结构；,层次模型与CAD系统的数据交换,CAD系统标准接口 VRMLIGESSTL自定义数据接口中性文件接口,曲面零件信息的提取,三角形生成规则取向规则共顶点规则三角剖分网格的生成 网格模型信息的边界表示 曲面信息的NURBS化,曲面零件信息的提取,三角形生成规则取向规则共顶点规则三角剖分网格的生成 网格模型信息的边界表示 曲面信息的NURBS化,typedef struct ProSurface surface; /曲面标志符 int n_vertices; /网格顶点数目 ProVector* vertices; /顶点向量头指针 ProVector* normals; /顶点法向量指针 int n_facets; /三角片数目 ProTriangle* facets; /三角形数组头指针 ProSurfaceTessellationData;,曲面零件信息的提取,三角形生成规则三角剖分网格的生成数据结构 离散精度控制网格模型信息的边界表示 曲面信息的NURBS化,ProPartTessellate（ ProPart part， double chord_ht;/最大弦高( 0 ) ， double angle_cntrl; / 角度控制 /(between 0.0 and 1.0) ProSurfaceTessellationData *output; /离散数据指针数组 ）,曲面零件信息的提取,三角形生成规则三角剖分网格的生成数据结构 离散精度控制网格模型信息的边界表示 曲面信息的NURBS化,取向规则,取向规则,取向规则,取向规则,取向规则,曲面零件信息的提取,三角形生成规则 三角剖分网格的生成 网格模型信息的边界(B-rep)表示法翼边结构对称结构曲面信息的NURBS化,/曲面类：class SurfTriangle : public Triangle public:/属性 int m_type;/面的类型 CString m_objName; /面所属零件的名称 CArray m_ownerFeatures; /面所属的特征 CArray m_ triList; /曲面对应的三角形链表 CArray m_loops; /环及边界的索引 int m_outLoopId; /存在的话，外环索引 CArray m_triFaces;/该面的面片id集合;,曲面零件信息的提取,三角形生成规则 三角剖分网格的生成 网格模型信息的边界表示 曲面信息的NURBS化 数据结构,typedef struct ptc_b_splsrf int deg2; /* 基函数阶数 (U 、V向). */ double *u_par_arr; /* U向结点数组. */ double *v_par_arr; /* V向结点数组.*/ double *wghts; /*控制点权重. */ Pro3dPnt *c_point_arr; /* 控制点数组 */ ProBsplinesrfdata;,曲面模型信息组织,曲面模型组织过程曲面离散化重合点、孤立点、退化三角形处理显示层数据点、三角片几何及拓扑信息 几何层数据点、线、面、环层次信息边界信息及边界角点环信息特征层信息（存储曲面特征信息） 几何特征信息字符特征信息,层次模型的建立步骤,模型信息转换文件 NFF文件：显示层数据 VDS文件：几何、特征层数据建模过程STEP1：建立零件结点P，依次读取零件层属性值，进行P的属性赋值；STEP2：遍历零件P的特征树，建立特征节点集合U=ui, i=1n,根据特征间的约束关系与父子关系在特征间添加约束关系与从属关系；STEP3：建立零件P与特征U间的数据映射关系；STEP4：遍历零件边界表示，建立几何层节点集合S=Si, i=1n,根据几何面间的约束关系与邻接关系在节点间添加约束关系边集合和邻接关系边集合；STEP5：建立特征集合U中每一元素与几何层节点集合S中元素的数据映射关系；STEP6：依次读入三角形信息，建立三角面片节点集合T=Ti, i=1n；STEP7：建立几何层节点中元素与显示层三角面片的数据映射关系；STEP8：调用显示层数据进行模型显示。模型信息扩充再组织,第三部分曲面网格简化及优化,曲面网格简化及优化,曲面网格的质量指标 平面度（Planarity） 弦高（Chord_ht）、角度（Angles）偏差（Deviation） 边界曲折程度的评定标准 曲面网格的局部优化网格细化 变形网格插补 网格简化狭长三角形 退化三角形过度细化三角形边界优化,基于几何元素删除的网格简化,简化元素点线面简化元操作顶点删除 边删除三角形删除简化空洞的剖分基于弦内角最大准则的优化剖分,弦内角定义及性质,弦内角的基本定义弦内角的等同性,弦内角最大化准则,弦内角的优化性判别 弦内角 映射关系 对角线交换任意弦内角最大化准则 最大弦内角 等价 优化三角剖分最大内角弦 等价 最小内接圆半径,任意弦内角最大化准则算法,算法步骤：步骤1 从散乱数据中按逆时针方向顺序读入简单多边形顶点，并建立双向循环链表；步骤2 计算出双向循环链表中每个结点的凹凸性；步骤3 对双向循环链表中每个凸结点Q，设与其前后结点P，R组成的三角形为 PQR，PQR内不包含多边形上其它的顶点；步骤4 若链表中还存在3个以上的结点，则转步骤2，否则转步骤5；步骤5 由链表中最后3个结点所对应的多边形顶点构成一个三角形，删除链表中最后3个结点。并形成一个三角剖分的三角形表TML； 步骤6 以任意弦内角最大化准则进行局部优化变换，得到优化的三角剖分(基于任意弦内角最大化准则判定的简单多边形三角剖分)； 步骤7 程序结束。,算法实例,原始网格,优化剖分网格,曲面网格的边界优化,曲面边界缺陷狭长三角形 退化三角形锯齿形边界平滑优化边界曲折程度标准平滑剖分算法,边界曲折程度的评定参数,边界轮廓的算术平均偏差（ ）在取样长度内，轮廓偏距绝对值的算术平均值,边界网格的平滑剖分,平滑剖分算法步骤：Step1 边界退化、狭长三角形优化剖分预处理。Step2 边界轮廓中线构造。从边界链表中隔点取样，重新构造边界链表，形成初始的中线链表，其余边界点构成。Step3 边界曲折度评定。计算边界曲线的轮廓算术平均偏差（Ra），根据预先设定的平滑限值Rmin检验曲线平滑程度。在边界链表中重复隔点取样35次并检验曲线平滑程度，如果每次取样得到的中线链表和不平点链表都满足平滑检验要求，可认为边界平滑程度较好，算法结束；否则，进行下一步。Step4 平滑剖分。遍历不平点链表。对T型三角点，从当前网格面中删除该边界三角形，并将其添加到边界对应的邻接曲面中；对S型三角点，以点为中心，按逆时针顺序依次对四边形、进行对角线交换，得到剖分后的三角形、和，将剖分前的相应三角形删除，同时将、加入当前网格面，加入邻接网格面。Step5 以中线链表作边界链表，重复步骤25，直至完成平滑剖分。,平滑剖分实例,曲折边界,平滑剖分前,平滑剖分后,优化边界,第四部分拓扑网格的几何重建,变形网格的几何重建,曲面造型Bezier曲面B样条（NURBS）曲面几何重建Bezier曲面重建：复杂边界形状曲面的拟合 三角域Bezier重建矩形域Bezier重建B样条（NURBS）曲面重建：CAD/CAM曲面造型Floater的双三次B样条几何重建,矩形域上双三次Bezier曲面重建,基本步骤：估算各网格顶点的法矢； 计算型值点周围的控制顶点，确定相邻曲面公共边界的控制顶点； 计算公共边界处的连续拼接条件； 实现各子三角形曲面片间的连续拼接。数据结构 class cTVectors4public:static double u, v, w, x, y, z, aa, bb;void CalABC1(SurfTriangle *pST);/三角域void CalABCD(SurfTriangle *pST);/矩形域void controlPointsA(AcGePoint3d *pTs, SurfTriangle *pST);void controlPointsB(AcGePoint3d *pTs, SurfTriangle *pST);void controlPointsC(AcGePoint3d *pTs, SurfTriangle *pST);,保形参数化双三次B样条重建,重建过程：网格边界顶点的参数化(弦长法)；网格内点的参数化(保形参数化)；构造整体C1连续分段三次B-B曲面(Farin的C-T算法)；从构造曲面上取得张量形式样点；构造插值于这些样点的C2连续双三次B样条曲面。,边界点链表的建立,边界点空点（NULL）面边界点步骤：获取一个边界点，做为边界链表出发点（三角形方法）；逆序遍历当前边界点所在三角形，找到下一边界三角形；遍历每一个三角形，并将逆序的边界点加入边界链表；重复上述3步骤，直到找到第一边界点，边界链表构建完毕。,保形参数化,边界点的弦长参数化公式：内点的保形参数化面积坐标：内点与边界点的映射：建立矩阵方程：追赶法求解三对角矩阵求得内点的参数化,插值曲面标量样点的求取,U、V向样条曲线数目的估算样条数 反求拟合曲面样点方法面积坐标：邻接面边界切矢、法矢约束的计算边界上插值样点的反求内部标量样点的反求,四边拓扑曲面的构造与角点的求取,四边拓扑曲面的构造 保形映射算法的四边拓扑分析 基于深度搜索的四边角点分类求取算法,四边拓扑曲面的构造,拓扑曲面的分类四边拓扑：圆形拓扑和三边拓扑：补足角点来构造四边域；N边拓扑：剔除多余角点构造四边域 拓扑类型与角点关系表,特征信息 拓扑类型,特征信息 拓扑类型,保形映射算法的四边拓扑分析,四边域保形映射的局限曲面四边拓扑要求四边角点分布要求重建曲面的影响分析边界点：“弦长法”映射不受影响 内点：网格变形 ，破坏映射的保形 均布角点 非均布角点,基于深度搜索的四边角点求取算法,深度：网格面上两点间的最小“距离”深度搜索：网格面上距某点深度 为定值的所有点的计算过程。宽度优先原则： 定值横向搜索优先的原则分类优化思想 拓扑构造分类 分类深度搜索算法描述Step1 网格顶点标号初始化；Step2 求取曲面边界角点 ； Step3 边界拓扑分类；Step4 .Step5 .根据拓扑分类，采取不同的步骤深度搜索“均布”四边角点,第五部分变形曲面重建模块的实现及应用,虚拟变形设计系统方案和功能框架,针对CAD软件在曲面设计上存在的不足，我们将虚拟现实技术引入CAD环境，研究了基于外部CAD模型的虚拟变形设计系统的方案和功能框架：人机交互层场景管理层变形设计层模型构造模块 变形设计模块曲面重建模块数据处理层,系统界面,曲面重建相关模块功能分析,Pro/Engineer下的数据转换模块,曲面重建相关模块功能分析,MDT下的标准数据文件（IGES）转换接口,系统信息交互模式,中性文件接口VDS/NFF文件外部模型CAD(Pro/Engineer)变形设计模块接口FIT文件变形设计模块曲面重建模块接口IGES曲面重建模块外部模型CAD(Pro/Engineer),曲面重建模块工作流程图,网格优化及重建实例,摩托车水箱盖曲面重点分析典型四边域上的双三次Bezier曲面及B样条曲面重建。 摩托车座壳曲面 摩托车座壳子曲面划分理想，本例分析网格曲面模型优化并作为圆盘拓扑曲面B样条重建的实例。,水箱盖曲面,原始模型及层次网格模型 水箱盖原始模型 水箱盖网格模型（1484个三角形）,水箱盖曲面,重建曲面 Bezier曲面重建（4452个B-B小曲面） B样条曲面重建（40张子曲面）,水箱盖曲面,IGES文件格式零件在MDT的实现,座壳曲面,建模 座壳原始模型 座壳网格模型（23711个三角形 ）,座壳曲面,座壳网格简化模型（834个三角形） B样条曲面重建（16张子曲面）,座壳曲面,Bezier曲面重建 IGES格式的Bezier曲面重建效果图（Pro/E下显示）,总结,面向变形设计的零件信息建模方法 引入了面向曲面变形与重建的层次网格模型；根据CAD模型与虚拟环境中层次网格模型信息表达的差异，给出了CAD模型中提取层次网格信息的方法及层次模型信息建立的过程。曲面网格模型的简化与优化 以几何元素删除方法为基础，实现了层次模型网格的简化方法 ；针对网格简化过程中出现的空洞剖分问题，提出了基于弦内角最大准则的局部网格优化方法 ；以边界轮廓算术平均偏差为边界曲折程度的定量描述 ，给出了曲折边界的平滑优化方法。拓扑网格的几何重建引入单边退化概念，进行四边域上的Bezier曲面重建