D0I:10.13374/i.issn1001053x.1992.03.031 第14卷第3期 北京科技大学学报 Vol.14 No.3 1992年5月 Journal of University of Science and Technology Beijing May 1992 一种新的造型曲面 王兵团 摘要:提出一种新的造型曲面一BRB样条曲面,此曲面兼有Bzi©r曲面和B样条曲 面的基本造型优点,造型性能综合指标斋。又给出新曲面的定义和几何性质,并讨论了 BRB样条曲面的拼接向题。 关键词:曲面,造型,样条 A New Kind of Surface for Modelling Wang Bingtuan' ABSTRACT:A new kind of surface for modelling-BRB spline surface-is in- troduced here,Comparing with both Bezier surface and B-spline surface,BRB spline surface can keep basic advantages of modelling and has a high composite target of modelling capability.The definition and geometric properties are given for BRB spline surface,Some adjacent problems between BRB spline surfaces are discussed, KEY WORDS:surface,modelling,spline 曲面造型是计算几何领域中的一个重要研究方向,其应用很广,并且还在发展之中。当 前用于曲面造型的两大主要方法是Bé2ir曲面和B样条曲面,这两种曲面各有长处,都是由 某类造型曲线与自身的直积构成的1)。本文提出一种新的造型曲面一Bézir-B样条曲面(简 称BRB样条曲面)。它由不同类造型曲线的直积生成,兼有前两种造型曲面的基本长处,其 造型性能的综合指标高于这两种造型曲面。在一定意义下,可以说它是对Bézir曲面和B样 1991-12-26收稿, ·数力系(Department of Mathematics and Mechanics) 383
第 卷第 期 年 月 北 京 科 技 大 学 学 报 。 。 一种新的造型 曲面 王 兵 团 ‘ 摘 要 提出一 种新的造型 曲面一 样条曲面 , 此 曲面 兼有 曲面和 样 条 曲 面的基本造型优点 , 造型性能综合指标高 。 又 给 出新曲面 的定 义 和 几 何性质 , 并 讨 论了 样条曲面 的拼接 向题 。 关健词 曲面 , 造型 , 择 条 牙 ” “ ’ 口口尸 一 一 亏 己 以 一 , ,泣 。 , , 、 曲面造型是计算几何领域中的 一个重要研究方向 , 其应用 很 广 , 并且还在发 展 之中 。 当 前用于 曲面造型的两大 主要方法是 乙 曲面和 样 条 曲面 , 这两种 曲面各有 长处 , 都 是 由 某类造型 曲线与 自身的直积构成的 〔 ‘ ’ 。 本文提 出一种 新的造型 曲面一 吮 卜 样条 曲面 简 称 样条 曲面 。 它 由不同类造型 曲线的直积生成 , 兼有前两种 造型 曲面的 基本长处 , 其 造 型性能 的 综合 指标高于这两种造型 曲面 。 在一定意义 下 , 可 以说 它是对 “ 曲 面 和 样 一 一 收稿 , 数力系 对 , , DOI :10.13374/j .issn1001-053x.1992.03.034
条曲面的综合和推广。 1BRB样条曲面的定义和性质 定义1:设{b,}j=0,1…n;j=0,1,…m为(m+1)(n+1)个空间点列,称m×n次参数曲面 P(u,w)=∑∑B,.()Fin(w)b:i 0≤u,w≤1 (1) 00 为m×n次BRB曲面片。 式中B.()=(行)u'(1-u)-‘为n次Bernstein多项式函数,Fn(w) 1 m(-) (m十1)(四+m-j-)为m次B样条基函数,其中(日),(m1)为二项式系数。6, 称为控制点或顶点,逐次用线段连接点列{b,}相邻两点组成的空间网格称为特征网格,它控 制着曲面P(“,)的形状。 在式(1)中,令b,(w)=∑b,F.(w),i=0,1,n,则有: P(u,w)=∑b,(w)B(u) 0≤4,w≤1 (2) 1■0 说明m×n次BRB曲面片可以如下构造:先把顶点列b,按第一个角标分类得到n+1类 顶点集b,=b,}a,i=0,1,…n,然后以顶点集b,构造参数为w的m次B样条曲线b,(w), =0,1,…n。任取w。∈[0,1]用顶点列{b,(wg)}:-o构造参数为4的n次Bezier曲线,这就得到 了P(u,w)的一条u线P(u,w),当w。取遍区间[0,1]时,就得出m×n次BRB曲面片P(u,w)。 这个构造过程能容易地想象出BRB曲面片的形状并做出相应的图形。 注意到曲面片P(u,w)的w线是m次B样条曲线,4线是n次Bezier曲线的特点,易知由式 (1)确定的BRB曲面片具有几何不变、凸包、端点、对称和逼近等良好的几何性质。 把同次的BRB曲面片沿w方向按B样条曲面片的连接方式连接起来,就得到BRB样条曲 面。 定义2:设{b,}8:1:k≥m,为(n+1)(k+1)个空间点列,称m×n次参数曲面。 P(u,w)=∑∑B,.(4)N,m+1(o)b,0≤u≤1,0≤w≤k (3) 00 为m×n次BRB样条曲面。 式中B:,()同定义1,N,m+1(w)是区间C0,]上递推形式的m次B样条基函数r) 显然,m×n次BRB样条曲面是C-1连续的,它不但具有BRB曲面片所有的几何性质, 而且还有局部修改特性。从式(3)看出,改动特征网格上的一个顶点最多只影响到以这点为 中心的m+1块m×次BRB曲面片。例如当m=n=2时,改动一个顶点只影啊到邻近的3块 2×2次BRB曲面片(图1)。 由于参数曲面上参数曲线网体现着参数曲面的特性,注意到m×n次BRB样条曲面的“,w 线的特点,可以期望BRB样条曲面会反映出Bézir曲面和B样条曲面所具有的一些性质,并 呈现出它们的双重特性。这样,Bézier曲面和B样条曲面的长处在一定意义下被融为一体, 使得BRB样条曲面具有更丰富的造型功能。 384
条曲面的绘合和推广 。 样条曲面 的定义 和性质 定义 , 设 ‘ , 二 。 , ” 勺 二 , , … 二为 十 。 十 个空 间点列 , 称“ 次参数 曲面 , 二 二 为, 火 。 次 曲面片 。 ‘ 。 了, 。 ‘ 了 。 。 , 。 毛 一 ‘ 一 式 中 , 。 了梦 。 ‘ 一 。 一 为。 次 , 多项式函数 , 二 , 。 、 , 认 产 而下乙 一 亨 ‘ 二 二 一 ,一 ‘ 。 为 , 次 样条基 函数 , 其中 梦 , “ 广 ‘ 为二项式 系数 。 “ , 称为控制点或顶 点 ,逐次用线段 连 接点列 ‘ ,, 相 邻两 点组成的 空间网格称为特征 网 格 , 它 控 制着曲面 。 , 。 的 形状 。 在式 中 , 令“ 。 ,, ,, 。 , , , … 。 , 则 有 了 。 , 。 二 艺 ‘ 。 ‘ 。 “ “ , “ 说明, 次 曲面片可以如下构造 先把顶点列 ‘ , 按第一个角标分类得 到 十 类 顶点集‘ , 谧 ,, , 。 , ‘ , , … , 然后 以顶点 集 “ 构造参数为,的, 次 样条 曲 线 ‘ , , 才 , , … 。 任取, 。 〔 〔 , 〕用 顶点 丈 ‘ , 。 李一 。 构造参数为 。 的 。 次 ‘ 曲线 , 这就得到 , 。 的一条。 线 。 ,甜 。 , 当, 。 取遍区间 , 〕时 , 就得 出“ 次 曲面 片 , ‘ 这个 构造过程能容易地 想 象 出 曲面片的 形状并做 出相应 的图形 。 注意到曲面片 。 , , 的“ 线是” 次 样条 曲线 , 。 线是 。 次 能 曲线的特点 , 易知 由 式 确定的 曲面片具有 几何不变 、 凸包 、 端点 、 对 称和逼 近等 良好的 几何性质 。 把 同次的 曲面片沿。 方向按 样条 曲面片的连 接方式连 接起 来 , 就得到 样 条 曲 面 。 定义 设 铸 ‘ , 吕 二特 , 为 舟十 个 空间点 列 , 称 义 。 次参数 曲面 。 , 。 二 乙 ‘ 。 ,。 , 。 ‘ , 簇 “ 毛 , 。 儿 ‘ 一 为 。 次 样条曲面 。 式 中 ‘ 。 同定义 , ,。 十 , 是 区’ 〔 , 〕上递推形式的 协 次 样条基 函数 ‘ “ ’ 。 显 然 , 。 次 样条曲面是 ” 一 ’ 连 续的 , 它 不但具有 曲面片所有的 几何性质 , 而 且还有局部修改特性 。 从式 看出 , 改动特征网格上的 一个顶点最多只影响到 以这 点为 中心的 十 块 。 次 曲面 片 。 例如 当, 二 。 时 , 改 动一 个顶点 只影 响 到 邻 近 的 块 次 曲面片 图 。 由于参数 曲面上参数曲线 网体现着参数曲面的特性 , 注意到 , 次 样条 曲面 的 。 , 。 线的 特点 , 可以期望 样条曲面会反映 出 能 曲面和 样条 曲面所具有的一些 性 质 , 并 呈 现 出它们的双重特性 。 这样 , 能 曲面和 样 条 曲面的 长处在一定意义 下被融为一体 , 使得 样条 曲面具有 更丰富的造型功能
201J 2×2BRB patches 图1(a)2×2次BRB样条曲面 (b)改动顶点b,后的2×2次样条曲面 Fig.1 (a)2x2 BRB spline surface (b)The BRB spline surface by modifying vertex bij 为便于理解BRB样条曲面,现仅就m=n=3的BRB曲面片进行讨论。显然3×3次BRB曲 面片以两条3次B样条曲线和两条3次Bézier曲线为边界。w线是3次B样条曲线,u线是3次 Bézir曲线。特征网格的中间4个顶点控制着曲面片的内部形状。把3×3次BRB曲面片分别写 成矩阵形式和双3次Coons曲面形式,比较它们的系数便得出3×3次BRB曲面片的角点信息 和其特征网格顶点之间的关系: P(,)=合(6,t401+6) P.(,0=合g-b) P.(,)=2[bg11-bg+4(o1-b2)+(6g-6门 P.(,)=8[og18-bg)-(a1-b2] k=0,1:1=0,1 27 U12 2 b:2 011 b21 b11 b35 b2 b13 13 b3 bi0 b20 b30 010 U的 bos (a)双3次Bezicr曲面片 (b)双3次B样东曲面片 385
卜 七砖甲份吟艺 片尸刀 哪 巍 圈 月 匕“ 味城沁‘ 州禅扫“ 喊 户 代 ’卜 、 ’ 飞 、 ‘ 、 份 户 ‘ , ”即… ’考 沪 、 夕天全妞生,曰‘ ‘ 一,动翻‘ ,尸目 琴 、 ‘ , 分 尹 臼 俨 嗽‘ 之七‘ 知二书从刀日以 口 曰 巴曰益‘ 曰鞠币一台格曰哭夕 图 次 样条 曲面 改动 顶点 ‘ ,后的 次 样条 曲面 一 、 为便于理解 样 条 曲面 , 现仅就 二 月 二 的 曲面 片进 行讨论 。 显然 次 曲 面 片以两 条 次 样条 曲线和两 条 次 曲线 为边界 。 二 线 是 次 样 条 曲 线 , 线 是 次 既 曲线 。 特征网格的 中间 个 顶点 控制着 曲面 片的 内部形状 。 把 次 曲面片分别 写 成 矩阵形式和双 次 曲面 形式 , 比较 它 们的 系数 便 得 出 次 曲面片的 角点 信 息 和 其特征 网格顶点之 间的 关系 付 一 ‘ 一二一 ‘ , ‘ 吞 、 , , , ,十 , , 一 。 , 一一 一 二 , 。 , , , 一 吐 ‘ 、 , , 一 。 ,十 , 石 。 十 , 一 。 , , 。 、 , 〔 。 十 一 。 , 一 吞 。 十 , ‘ 一 , , 〕 沪尹 、﹄沂﹄陌 二 , 二 , 双 次 么 曲面 片 双 次 样条曲面 片
712 022 01z b5 00 b23 b23 11 D10 b20 0211 60 010 b00 b oo b05 (C)3×3次BRB曲面片 (d)改动顶点b12后的3×3次BRB曲面片 图?一些不同的曲面片 Fig.2 Some different patches 这个关系式指出了网格顶点对曲面片的影响。例如顶点b11,b12,b21,b22对曲面片的角 点扭矢和u线的切矢有影响,从而它们对曲面片的内部和Bézir边界线产生影响,调整这些 顶点就能构造不同形式的曲面片。图2给出在同一特征网格下的3×3BRB、Bezier和B样条曲 面片。从图中可知BRB曲面介于Bezier曲面和B样条曲面之间。为便于比较,图中还画出了 改动某控制点后的BRB曲面片。对一般的m×n次BRB曲面片也有类似的性质。 2 BRB样条曲面片的连接 在实用中,为了不减弱网格对曲面的控制并得到满意的曲面,可取BRB样条曲面的次数 m,n≤5并用若干这样的曲面片连接起来进行造型。使用复合曲面的关键是处理好各曲面片 的连接问题。B样条曲面具有整体配置顶点的特点,没有连接问题,但适应性有限制;Bézier 曲面片装配灵活,造型适应性强,但曲面片的光滑连接麻烦,处理不当就容易出问题3),这 主要是由Bezier曲面片的连接是“,w两个参数方向进行而引起的,BRB样条复合曲面也有连接 问题,但其只在一个“参数方向进行,于是 BRB样条曲面的连接比Bezier曲面片简单且较 少出问题(如图3)。 现在讨论次数分别为m×n和m×I的两张 BRB样条曲面片的连接条件,次数的这种限制 基本上不影响曲面的造型。注意到BRB样条曲 The boundaries of B3RB patches 面片的连接实际上可转化为相对应的BRB曲面 国3由BRB样条曲面片A和B构成的复合曲面 片沿u参数方向的连接,这是Bezier曲面片的单 Fig.3 The surfacc made by a mosaic of 向连接,类似Bezier曲面片的连接条件r4),易 BRB splinc patch A and B 得BRB样条曲面片下列连接条件: 设m×n和m×1次BRB曲面片S:和S,为 386
次 曲面 片 改 动 顶点 后的 次 曲面 卜 图 一 些不 同 的曲面 片 这个 关系式指 出了网格顶 点 对 曲面片的 影响 。 例如顶点 , , , ‘ , ‘ , 如 对 曲面 片 的 角 点 扭矢和“ 线的切 矢有 影响 , 从而它们对 曲面 片的 内部和 亡 “ 边界线产生影响 , 调 整 这 些 顶点就能 构造不同形式的 曲面片 。 图 给 出在同一特征 网格下的 、 亡 和 样条曲 面片 。 从 图中可知 曲面介于 扛 曲面和 样条 曲面之 间 。 为 便于 比较 , 图中还 画 出 了 改动某控制点后的 曲面片 。 对一 般的 , , 次 曲面片也有 类似的性质 。 样条曲面 片的连接 在实用 中 , 为 了不减弱 网格对 曲面的控制并得到满意的 曲面 , 可取 样条 曲面的次数 二 ,, 并用若干这样的 曲面片连接起来进 行造型 。 使用 复合 曲面的关键是处理好各 曲 面 片 的连 接问题 。 样条曲面具有整体配置顶点的特点 , 没有连 接问题 ,但适应性有限制 此 曲面片装配灵 活 , 造型适应性强 , 但曲面片的光滑连 接麻烦 , 处理不当就容易出问题 〔 “ ’ ,这 主要是 由 比 曲面片的连接是 “ ,,两个参数方 向进行而 引起的 样 条复合曲面也有连 接 上 ‘ 丈玉 阵。 由 样条曲面片 和 构成的复合曲面 问题 , 但其只在一个 “ 参数方 向 进 行 , 于 是 样条曲面的连 接比 曲面片简单且较 少 出问题 如图 。 现在讨论次数分别为。 。 和 的两 张 样 条曲面片的连接条件 , 次数的这种限 制 基本上不影响 曲面的造型 。 注意到 样 条曲 面片的连 接实际上可转化为相对应的 曲面 片沿 参数方向的连接 ,这是 比 曲面片的单 向连 接 , 类似 曲面片的连 接条件 〔 弓 ’ , 易 得 样 条曲面片下列连接条件 设 。 和, 次 曲面片 和 为 图扒
S1:P(u,w)=∑∑B.(u)F,n(w)bi 04,w1 ▣0f S2:Q(u,w)=∑∑B:(u)F,n(w)h1 0u,w1 t量0j=0 再设m×n及m×l次BRB样条曲面A和B为 A:r(u,2w)=∑∑B.()Nm1(w)b 0u1 i=0j=0 B:q(u,w)=>B(u)N(w)h:i 0wk i=0i-0 则有: 定理1:曲面片S:和Sz沿公共边界C连续的充要条件是: b.4=h。,j=0,1,2,…m (4) 龙厘2:曲而片S1和S2沿公共边界线GC'连续的允要条件是D(w)Q.(0,w)=E(w)P.(1,w) +F(w)P.(1,w)及式(4)同时成立。其中D(w),E(w)和F(w)分别是次数不超过2m-1,2m 和2m-1的多项式,P.6,P.8胎,Q.-品。 .面定理中C表示位置连续连接,GC′表示光滑连接。由BRB样条曲面片和BRB曲面 片的关系得到BRB样条曲面片的连接条件如下。 定理3:曲面片A与B沿公共边界C连续的充要条件是b,=h8,j=0,1,2,…k。 定理4:曲面片A与B沿公共边界GC'连续的充要条件是其相对应的每一对BRB曲面片满足定 理2的条件。 在定理2中,可取D(w)=1,E(w)=4。+41w,F(w)=元,其中4。,1和元为待定 常数。此取法有3个自由度,常可以满足实用中的光滑连接要求。利用这种取法,当m= n=3时,得到BRB曲面片S1和S2的光滑连接条件为: [b31=hoi j=0,1,2,3 Q.(0,w)=(4。+41w)P(1,w)+元P.(1,w) (5) 将式(5)的第二式展开并比较w的次数,得: h-h0=(b0-b20+Co(-11b3g+18b:-9b3+2b) +41(7b30-8b31+b32)门 h1-A1=a(b1-b2)+8C,(-2bs0-3b1+6bga-bg) +41(-2b3。+4b1-2b32)门 (6) A:-h=A(bga-bg:+8C4,(b-60g1+3ba+2b) +41(b30-8b31+7b32)〕 h-A=2(b:-bg)+8CA,(-2ba+96s1-186ga+11b) +4。(-2b3。+10b51-26b32+18b3s)门 387
, , 。 乙 乙 ‘ 。 “ ,, , ‘ , 气 , 留 , 。 二 乙 乙 ‘ ,。 刃 ‘ , 叹 “ , 簇 再 设。 , 及 ,。 次 样 条曲面 和 为 , 叼 艺 艺 ‘ , ,。 。 六 又 一 一 五 叮 , , 二 艺 乙 “ 十 ,。 份 , 、 二留 丈二 ‘ 一 则 有 定 理 曲面片 和 沿公 共边 界 ” 连 续的 充要 条 件是 。 , 入 。 , , , , “ … 机 定 理 曲面片 , 和 沿 公 共边 界线 ‘ 连 续的 充 要 条 件是 川 , 川 , , 劝 尸 。 尸 “ , 。 及 式 同 时成立 。 其中 , 和 川 分 别是 次数不 超 过 , 一 , 加 和 。 一 的 多 项 式 , 尸 二 哄 , , 。 二 塑 , 口 嘿 , ‘ ,一 曰 一 ’ 一 上 ,,功 少 兴 少、 ’ 一 “ ’ 一 沙 夕 ” “ ’ 上面 定理 中 。 表示 位置连 续 连 接 , ‘ 表示 光 滑连接 。 由 样条曲面片和 曲面 片的 关系得到 样 条曲面片的连 接条件如 下 。 定理 曲面片 与 沿公共边 界 。 连续的充要 条件是 共 乙 。 , , , “ · ” 。 定理 曲面 片 与 沿公 共边 界‘ ‘ 连 续 的 充要 条件是 其相对 应的每一对 曲面片满足定 理 的 条件 。 在定理 中 , 可取 。 , 二 户 。 那 二 , 。 久, 其中声。 , 声 和 几 为待定 常数 。 此取法 有 个 自由度 , 常 可 以满足实用 中的 光滑连 接要 求 。 利 用 这 种 取法 , 当川 。 时 , 得到 曲面片 , 和 的 光滑连接条件为 , “ 。 , , , , 火 。 , 切 声。 产 , 二 ,功 几 。 , 田 将式 的 第二式展开 并比较 。 的 次数 , 得 “ 。 一 。 。 “ “ 。 一 , 。 命 〔 。 ‘ 一 ,“ 。 ‘ ‘ 工 一 ” ‘ ‘ ‘ , 料 , 。 一 、 石 〕 一 一 。 久 一 , 一 〔 拜。 一 ‘ 。 一 吞 、 一 , 群 , 一 。 」 一 〕 入 一 。 几 一 〔拼。 。 一 产 石 。 一 〕 一 一 。 久 。 一 〔 声。 一 。 一 , 拼 。 一 。 一 〕
于是式(5)的第一式及式(6)一起构成了由网格顶点确定的光滑连接儿何条件,它的 适用范围较大。在式(6)中取4。=1=0时,得到更简单的光滑连接条件: b3)=h。y j=0,1,2,3 h1-h。y=1(b1-b2) (7) 式(7)要求S:和S,交于边界的特征多边形“楼边”必须共线且跨界网格交点处切线比 值相等(图4a),但式(6)没此限制(图4b)。 b25b55=ho3h11 n) b22 b21 b33-ho3 》 bos b21 32h02 b31=o1 b20 b30=0 boo h30 b00 b3o=h oo b20 hio S1 52 图4S1和S2的特征多边形 Fig.4 The characteristic of S1 and Sa 3结 论 BRB样条曲面兼有Bezier曲面和B样条曲面的基本造型优点,曲面片的光滑连接比Bezier 曲面片简单,造型性能综合指标高,具有可行性,实用性。 致谢:本文受到刘饮圣老师的热情支持和指子,陈仁宽老师对此文提出了宝贵意见,谨此深表谢意。 参考文献 1苏布青,刘鼎元。计算几何,上海科技出版社,1981 2 Willian J,Gordon and Richard F,Riesenfeld,B-Spline Curves and Surfa- ces,Computer Aided Geometric Design Academic Press,1974 3 Faux I D Pratt M J,Computational Gcometry for Dcsign and Manufacturc, John Wiley and Sons,1979 4刘鼎元。应用数学学报,1986,9:432~442 3888
子是式 的第一式及式 一起构成了由网格顶点确定的光滑连接几何 条 件 , 它 的 适用 范围较大 。 在式 中取 产。 拼 时 , 得到更简单的光滑连接条件 ‘ ,, “ 盖。 , · , , , 气 , 一 , 久 一 式 要求 和 交于边界的特征多边形 “ 棱边 ” 必须共线且跨界网格交点处 切线 比 值相等 图 , 但式 没此限 制 图 。 图 和 的 特征 多边形 ‘ 万 结 论 样条曲面兼有 既 曲面和 样条 曲面的基本造型优点 , 曲面片的 光滑连 接 比 行 汀 曲面 片简单 , 造型性能 综 合指标高 , 具有可行性 , 实用性 。 致谢 本文受到刘钦 圣老师的热情支持和指导 , 陈 仁宽老师对此文提出了宝贵意见 , 谨此深表谢意 。 参 考 文 献 苏布青 , 刘 鼎元 。 计 算几何 , 上海科技 出版社 , , , 一 盯 , , 。 一 , , 刘 鼎元 。 应用数学学 报 , , 吐 匀嘴