【摘要】 一种保持交点约束的曲线线架整体光顺的方法，属于曲线曲面的CAD领域，其特征在于， 若曲线的控制顶点个数n与其所含交点个数m之比小于1.5，则首先增加曲线的控制顶点数 目；记录曲线交点的参数值；指定每根曲线在整体线架中的光顺权值；生成线架的总体约束 优化函数；把约束优化问题转化成无约束优化问题进行求解。基于该方法可以使线架中的曲 线曲率变化更加均匀，并且保持了横纵曲线间原有的交点约束，真正达到了线架整体光顺的 效果。 【专利类型】发明授权 【申请人】清华大学 【申请人类型】学校 【申请人地址】100084北京市100084-82信箱 【申请人地区】中国 【申请人城市】北京市 【申请人区县】海淀区 【申请号】CN200810104209.7 【申请日】2008-04-17 【申请年份】2008 【公开公告号】CN100583106C 【公开公告日】2010-01-20 【公开公告年份】2010 【授权公告号】CN100583106C 【授权公告日】2010-01-20 【授权公告年份】2010.0 【IPC分类号】G06F17/50 【发明人】刘永进; 臧彧; 来煜坤; 胡事民 【主权项内容】1.一种保持交点约束的曲线线架光顺的CAD方法，其特征在于，所述方法是在计算机上按照如下步骤实现的： 步骤(1)，初始化： 输入要进行光顺的原始线架，并设定 横向线的集合C＝{ci|i＝1，2…，u}，纵向线的集合D＝{di|i＝1，2…，v}， pij为横向线ci与纵向线dj的交点， 每条曲线采用B样条的表示方法，对于第i条横向线其形式为： 其中n(ci)是横向线ci所含有的控制顶点个数减1，k(ci)是横向线ci的曲线次数， 是横向线ci中序号为l的控制顶点，l＝0，1，…，n(ci)，采用笛卡尔坐标表示， 是横向线ci中序号为l的B样条基函数，l＝0，1，…，n(ci)，B样条基函数使用到如下节点向量： 符号t表示节点向量中的一个分量，取值在(0，1]间， 表示节点向量中下标序号从k(ci)+1到n(ci)的共n(ci)-k(ci)个分量值，取值均在(0，1]间； 对于第j条纵向线其形式为 其中n(dj)是纵向线dj所含有的控制顶点个数减1，k(dj)是纵向线dj的曲线次数， 是纵向线dj中序号为l的控制顶点，l＝0，1，…，n(dj)，采用笛卡尔坐标表示， 是纵向线dj中序号为l的B样条基函数，l＝0，1，…，n(dj)，B样条基函数使用到如下节点向量： 表示节点向量中下标序号从k(dj)+1到n(dj)的共n(dj)-k(dj)个分量值，取值均在(0，1]间； 步骤(2)，按照如下步骤建立线架的总体优化目标函数： 步骤(2.1)检查各条横向线的控制顶点个数n以及其所含交点个数m，若n/m＜1.5则用B样条曲线的升阶或顶点插入算法增加其控制顶点个数使其至少为所含交点个数m的1.5倍，反之执行下一步骤； 步骤(2.2)对于纵向曲线，重复步骤(2.1)； 步骤(2.3)对于横向线ci和纵向线dj的交点Pij，记录其在两条曲线中各自的参数坐标值pij(c)和pij(d)； 步骤(2.4)建立横向线的优化目标函数： 其中： 和 为节点向量中下标分别为n(ci)+1和k(ci)-1的分量，ci是变化后的第i条横向线，其形式为 这里的 是ci中序号为l的控制顶点，l＝0，1，…，n(ci)，ci的其它内容与ci相同， 是变化后曲线ci的二阶导数平方的积分，λi(c)为用户指定的该横向线的权值，其取值在0、1之间，α是罚函数因子，用于控制曲线控制顶点的变化范围，为设定值； 步骤(2.5)建立纵向线的优化目标函数： 其中： 和 为节点向量中下标分别为n(dj)+1和k(dj)-1的分量，dj是变化后的第j条纵向线，其形式为 这里的 是dj中序号为l的的控制顶点，l＝0，1，…，n(dj)，dj的其它内容与dj相同， 是变化后曲线dj的二阶导数平方的积分，λj(d)为用户指定的该纵向线的权值，其取值在0、1之间，α是罚函数因子其含义同步骤(2.4)； 步骤(2.6)建立线架的总体约束优化目标函数： 其中 为优化后第i条横向线与第j条纵向线的交点，pij(c)为该交点在第i条横向线中的参数值，pij(d)为该交点在第j条纵向线中的参数值； 光顺优化过程是按如下步骤实现的： 步骤(2.6.1)用户向计算机输入控制顶点所允许的最大变化范围ε，并且设置罚因子α的初始值为10-5； 步骤(2.6.2)加入u·v个拉格朗日乘子δij，1≤i≤u，1≤j≤v，把上述约束优化问题转化为无约束优化问题： 其中未知变量为控制顶点及δij，此式是关于未知变量的二次函数； 步骤(2.6.3)用二次函数求极小值的方法进行求解，取得函数的极小值和对应的未知变量； 步骤(2.6.4)对于求解出的 和 若 且 则优化结束，反之另α＝2×α，返回步骤(2.6.2)重新求解直到满足优化结束条件为止。 【当前权利人】清华大学 【当前专利权人地址】北京市100084-82信箱 【专利权人类型】公立 【统一社会信用代码】12100000400000624D 【引证次数】2.0 【他引次数】2.0 【家族引证次数】2.0
