文章快速检索     高级检索
  浙江大学学报(理学版)  2017, Vol. 44 Issue (3): 307-313  DOI:10.3785/j.issn.1008-9497.2017.03.011
0

引用本文 [复制中英文]

李军成, 宋来忠. 利用带形状参数的有理势函数构造基于Metaball的过渡曲线[J]. 浙江大学学报(理学版), 2017, 44(3): 307-313. DOI: 10.3785/j.issn.1008-9497.2017.03.011.
[复制中文]
LI Juncheng, SONG Laizhong. Construction of transition curves based on Metaball technique using rational potential function with a shape parameter[J]. Journal of Zhejiang University(Science Edition), 2017, 44(3): 307-313. DOI: 10.3785/j.issn.1008-9497.2017.03.011.
[复制英文]

基金项目

国家自然科学基金资助项目(11171181);湖南省教育厅资助科研项目(14B099)

作者简介

李军成(1982-), ORCID:http://orcid.org/0000-0002-1904-4068,男, 博士, 副教授, 主要从事计算机辅助几何设计及其应用研究,E-mail:lijuncheng82@126.com

文章历史

收稿日期:2016-01-26
利用带形状参数的有理势函数构造基于Metaball的过渡曲线
李军成1 , 宋来忠2     
1. 湖南人文科技学院 数学与金融学院, 湖南 娄底 417000;
2. 三峡大学 理学院, 湖北 宜昌 443002
摘要: 利用现有势函数构造基于Metaball的过渡曲线, 此过渡曲线无法兼具拟高阶连续性与形状可调性.针对这一问题, 巧妙地从一种带形状参数的曲线模型出发, 构造一类带形状参数的有理势函数, 并研究该势函数的性质.所构造的有理势函数具有统一的数学模型, 不仅能使过渡曲线在端点处达到拟Ck连续, 而且还可通过修改形状参数的值调整过渡曲线的形状.实例表明, 通过调整有理势函数的次数及形状参数的取值可构造出满足不同拟连续性且形状不同的过渡曲线, 以满足实际应用需要.
关键词: 势函数    形状参数    Metaball技术    过渡曲线    形状调整    
Construction of transition curves based on Metaball technique using rational potential function with a shape parameter
LI Juncheng1 , SONG Laizhong2     
1. College of Mathematics and Finances, Hunan University of Humanities, Science and Technology, Loudi 417000, Hunan Province, China;
2. College of Science, China Three Gorges University, Yichang 443002, Hubei Province, China
Abstract: Using the existing potential functions to construct transition curve based on Metaball technique, it can not have both quasi high order continuity and shape adjustability. To solving this problem, a rational potential function with a shape parameter is constructed ingeniously from a curve with shape parameters. Some properties of the rational potential function are studied. The rational potential function is expressed as a unified mathematical model, which can not only make the transition curve achieve quasi Ck continuity at the end points, but also adjust the shape of the transition curve by modifying the value of the shape parameter. Some examples showed that the transition curves with different continuities and shapes could be constructed by changing the degree and the shape parameter of the rational potential function, which could be used to meet different needs of the practical application.
Key words: potential function    shape parameter    Metaball technique    transition curve    shape adjustment    

在计算机辅助几何设计(CAGD)中, 曲线造型一直以来都是重要的研究课题.随着几何造型工业的发展, 人们往往需要改变曲线的形状以满足实际工程的需要.因此, 带形状参数的曲线造型方法逐渐成为研究热点.这些方法的主要目的是在曲线模型中引入形状参数, 并通过修改形状参数的取值实现对曲线形状的调整.例如, 带形状参数的Bézier曲线[1-2], 带形状参数的B样条曲线[3-4], 带形状参数的三角曲线[5-6]等.

在曲线造型中, 过渡曲线的构造在许多实际工程问题中有着十分广泛的应用[7-9].为满足过渡曲线的设计要求, 李凌丰等[10]提出了一种基于势函数与Metaball技术构造过渡曲线的方法, 该法采用WYVILL等定义的六次多项式势函数构造能光滑连接2条曲线的过渡曲线, 虽然对被连接曲线的种类没有限制, 但所构造的过渡曲线在端点处的连续性较低,且其形状无法调整.为解决此问题, 高晖等[11]构造了2类势函数:第1类为可使过渡曲线在端点处达到拟Ck连续的多项式势函数; 第2类为可使过渡曲线在端点处达到拟C1连续且具有形状可调性的混合三角势函数.文献[11]构造的第1类势函数,虽然提高了过渡曲线在端点处的拟连续性, 但仍然无法调整其形状; 构造的第2类势函数,虽然过渡曲线的形状可调, 但在端点处的拟连续性较低.另外, 2类势函数构造过程较为烦琐,且均无统一的数学模型, 因此需要重新构造满足不同情形的势函数.与文献[10-11]类似, 刘华勇等[12-14]研究了基于调配函数的过渡曲线构造与连续性问题, 并给出了可使过渡曲线在两端点处满足拟Ck(k=0, 1, 2) 连续的调配函数.注意到为了使过渡曲线在满足一定连续性的条件下同时具有形状可调性, 文献[12-14]将2条被过渡曲线取为带形状参数的曲线模型.由于所选取的2条被过渡曲线均带有形状参数, 所以所构造的过渡曲线在不改变几何连续性的情形下,可通过其所带的形状参数进行形状调整.然而, 在大量实际过渡曲线构造问题中, 2条被过渡曲线的模型往往并不特定.为此, 李军成等[15]构造了一类带参数的多项式势函数, 并研究了该势函数在构造过渡曲线中的应用.虽然文献[15]提出的方法对2条被过渡曲线的种类没有限制, 而且还可利用势函数所带的参数对过渡曲线的形状进行调整, 但所构造的过渡曲线在两端点处仅能满足拟C2连续, 不适合连续性要求更高的场合.

为满足更高要求的过渡曲线设计, 本文在文献[15]的基础上, 巧妙地从一种带形状参数的曲线模型出发, 构造一类带形状参数的有理势函数, 该势函数不仅具有统一的数学模型, 而且可使过渡曲线同时具有拟Ck(k≥1) 连续性和形状可调性.

1 问题的描述

在解决大量实际工程问题的过程中, 常常需要将数段零散的曲线段连接成一个整体, 且要求不破坏曲线段的光滑度, 基于Metaball的过渡曲线构造就能满足这一要求[10].文献[10-11]描述了基于Metaball构造过渡曲线的问题:给定平面上相交于点C的2条参数曲线P(t)与Q(t), 2条曲线的端点分别记为AB, 如图 1所示, 希望构造1条能光滑连接AB点的过渡曲线R(t).

图 1 构造基于Metaball的过渡曲线 Fig. 1 Construction of transition curves basedon Metaball technique

针对上述问题, 文献[11]给出了过渡曲线的方程:

$ \mathit{\boldsymbol{R}}\left( t \right) = \mathit{\boldsymbol{P}}\left( t \right)f\left( t \right) + \mathit{\boldsymbol{Q}}\left( t \right)\left( {1 - f\left( t \right)} \right), $ (1)

其中,0≤t≤1, f(t)为势函数.

由式(1) 可知, 为使所构造的过渡曲线在2个端点处达到拟Ck连续, 势函数f(t)(0≤t≤1) 在端点处须满足

$ \left\{ \begin{array}{l} f\left( 0 \right) = 1,\;\;\;\;f\left( 1 \right) = 0,\\ {f^{\left( i \right)}}\left( 0 \right) = {f^{\left( i \right)}}\left( 1 \right) = 0,\;\;\;i = 1,2, \cdots ,k. \end{array} \right. $ (2)

由于构造过渡曲线R(t)的目的是使曲线P(t)连续平滑地过渡到Q(t), 因此过渡曲线R(t)在靠近曲线P(t)处应与P(t)具有相似的形状, 而在靠近曲线Q(t)处则应与Q(t)具有相似的形状.由式(1) 可知, 当给定曲线P(t)与Q(t)时, 过渡曲线R(t)的形状完全由势函数f(t)决定, 故势函数的选取是构造过渡曲线的关键.本文的主要目的就是构造一类同时具有多种特性的势函数, 并将其用于构造基于Metaball的过渡曲线.

2 势函数的构造及其性质 2.1 势函数的构造

构造势函数的主要思想:首先选取一种恰当的带形状参数的曲线模型, 然后利用曲线模型在端点处满足性质来构造具有某种特定要求的势函数.

文献[16]构造了一种带5个参数k, ωi, αi(i=1, 2) 的曲线, 其表达式为

$ \mathit{\boldsymbol{b}}\left( t \right) = \sum\limits_{i = 0}^3 {{b_i}\left( t \right){\mathit{\boldsymbol{V}}_i}} , $ (3)

其中,0≤t≤1, ViRd(d=2, 3;i=0, 1, 2, 3) 为给定的控制顶点, bi(t)(i=0, 1, 2, 3) 为调配函数, 满足

$ \left\{ \begin{array}{l} {b_0}\left( t \right) = \frac{{\left( {1 - {\alpha _1}} \right){{\left( {1 - t} \right)}^{k + 1}}}}{{R\left( t \right)}},\\ {b_1}\left( t \right) = \frac{{\left( {{\alpha _1}{{\left( {1 - t} \right)}^{k + 1}} + {\omega _1}{{\left( {1 - t} \right)}^{k + 1}}t} \right)}}{{R\left( t \right)}},\\ {b_2}\left( t \right) = \frac{{\left( {{\omega _2}\left( {1 - t} \right){t^{k + 1}} + \left( {1 - {\alpha _2}} \right){t^{k + 1}}} \right)}}{{R\left( t \right)}},\\ {b_3}\left( t \right) = \frac{{{\alpha _2}{t^{k + 1}}}}{{R\left( t \right)}}. \end{array} \right. $ (4)

其中,

$ \begin{array}{*{20}{c}} {R\left( t \right) = {{\left( {1 - t} \right)}^{k + 1}} + {\omega _1}{{\left( {1 - t} \right)}^{k + 1}}t + }\\ {{\omega _2}\left( {1 - t} \right){t^{k + 1}} + {t^{k + 1}},} \end{array} $

整数k≥1, 实数ωi > 0(i=1, 2), α1∈[0, 1), α2∈(0, 1].

由式(3) 定义的曲线在端点处满足

$ \left\{ \begin{array}{l} \mathit{\boldsymbol{b}}\left( 0 \right) = {\mathit{\boldsymbol{V}}_0} + {\alpha _1}\left( {{\mathit{\boldsymbol{V}}_1} - {\mathit{\boldsymbol{V}}_0}} \right),\\ \mathit{\boldsymbol{b}}\left( 1 \right) = {\mathit{\boldsymbol{V}}_2} + {\alpha _2}\left( {{\mathit{\boldsymbol{V}}_3} - {\mathit{\boldsymbol{V}}_2}} \right),\\ {\mathit{\boldsymbol{b}}^{\left( i \right)}}\left( 0 \right) = {\left( { - 1} \right)^{i - 1}}{A_i}\left( {{\mathit{\boldsymbol{V}}_1} - {\mathit{\boldsymbol{V}}_0}} \right),\\ {\mathit{\boldsymbol{b}}^{\left( i \right)}}\left( 1 \right) = {B_i}\left( {{\mathit{\boldsymbol{V}}_3} - {\mathit{\boldsymbol{V}}_2}} \right), \end{array} \right. $ (5)

其中,1≤ik, Ai是与i, ω1, α1有关的常数, Bi是与i, ω2, α2有关的常数.

显然, 式(3) 对应的函数可表示为

$ f\left( t \right) = {b_0}\left( t \right){y_0} + {b_1}\left( t \right){y_1} + {b_2}\left( t \right){y_2} + {b_3}\left( t \right){y_3}, $ (6)

其中, bi(t)(i=0, 1, 2, 3) 为由式(4) 定义的4个调配函数, yiR(i=0, 1, 2, 3).

由式(5) 可知, 式(6) 定义的函数在端点处满足

$ \left\{ \begin{array}{l} f\left( 0 \right) = {y_0} + {\alpha _1}\left( {{y_1} - {y_0}} \right),\\ f\left( 1 \right) = {y_2} + {\alpha _2}\left( {{y_3} - {y_2}} \right),\\ {f^{\left( i \right)}}\left( 0 \right) = {\left( { - 1} \right)^{i - 1}}{A_i}\left( {{y_1} - {y_0}} \right),\\ {f^{\left( i \right)}}\left( 1 \right) = {B_i}\left( {{y_3} - {y_2}} \right). \end{array} \right. $ (7)

为了使所构造的过渡曲线在2个端点处达到拟Ck连续, 由式(2) 和(7) 则可得

$ \left\{ \begin{array}{l} {y_0} + {\alpha _1}\left( {{y_1} - {y_0}} \right) = 1,\\ {y_2} + {\alpha _2}\left( {{y_3} - {y_2}} \right) = 0,\\ {\left( { - 1} \right)^{i - 1}}{A_i}\left( {{y_1} - {y_0}} \right) = 0,\\ {B_i}\left( {{y_3} - {y_2}} \right) = 0. \end{array} \right. $ (8)

由式(8) 可得

$ \left\{ \begin{array}{l} {y_0} = {y_1} = 1,\\ {y_2} = {y_3} = 0. \end{array} \right. $ (9)

将式(9) 代入式(6), 并令ω1=ω2=ω, 可定义如下一类带形状参数的有理势函数.

定义1  对于整数k≥1, 实数ω > 0, 称关于变量t(0≤t≤1) 的函数

$ {f_k}\left( t \right) = \frac{{{{\left( {1 - t} \right)}^{k + 1}} + \omega {{\left( {1 - t} \right)}^{k + 1}}t}}{{{{\left( {1 - t} \right)}^{k + 1}} + \omega {{\left( {1 - t} \right)}^{k + 1}}t + \omega \left( {1 - t} \right){t^{k + 1}} + {t^{k + 1}}}} $ (10)

为带形状参数ω的有理势函数.

2.2 势函数的性质

带形状参数ω的有理势函数fk(t)具有如下性质:

性质1  端点性:势函数fk(t)在端点处满足

$ \begin{array}{*{20}{c}} {{f_k}\left( 0 \right) = 1,\;\;\;\;{f_k}\left( 1 \right) = 0,}\\ {f_k^{\left( i \right)}\left( 0 \right) = f_k^{\left( i \right)}\left( 1 \right) = 0,i = 1,2, \cdots k.} \end{array} $

证明  由势函数fk(t)的构造过程可知, 端点性显然成立.

注1  由性质1及式(1) 可得R(0)=P(0), R(1)=Q(1).此即表明, 在端点A处, 仅曲线P(t)对过渡曲线R(t)有影响; 而在端点B处, 仅曲线Q(t)对过渡曲线R(t)有影响.

另外, 由性质1和式(1) 并经简单推导可得R(i)(0)=P(i)(0), R(i)(1)=Q(i)(1), i=1, 2, …, k.表明势函数fk(t)可使过渡曲线R(t)在两端点处达到拟Ck(k≥1) 连续.例如, 使过渡曲线在两端点处分别达到拟C1, C2, C3连续的势函数为

$ \begin{array}{l} {f_1}\left( t \right) = \frac{{{{\left( {1 - t} \right)}^2} + \omega {{\left( {1 - t} \right)}^2}t}}{{{{\left( {1 - t} \right)}^2} + \omega {{\left( {1 - t} \right)}^2}t + {{\left( {1 - t} \right)}^2} + \omega \left( {1 - t} \right){t^2} + {t^2}}},\\ {f_2}\left( t \right) = \frac{{{{\left( {1 - t} \right)}^3} + \omega {{\left( {1 - t} \right)}^3}t}}{{{{\left( {1 - t} \right)}^3} + \omega {{\left( {1 - t} \right)}^3}t + \omega \left( {1 - t} \right){t^3} + {t^3}}},\\ {f_3}\left( t \right) = \frac{{{{\left( {1 - t} \right)}^4} + \omega {{\left( {1 - t} \right)}^4}t}}{{{{\left( {1 - t} \right)}^4} + \omega {{\left( {1 - t} \right)}^4}t + \omega \left( {1 - t} \right){t^4} + {t^4}}}. \end{array} $

性质2  中点性:fk(0.5)=0.5.

证明  由式(10) 可得

$ {f_k}\left( {0.5} \right) = \frac{{{{0.5}^{k + 1}} + \omega {{0.5}^{k + 2}}}}{{2\left( {{{0.5}^{k + 1}} + \omega {{0.5}^{k + 2}}} \right)}} = 0.5. $

注2  由性质1及式(1) 可得

$ \mathit{\boldsymbol{R}}\left( {0.5} \right) = \frac{{\mathit{\boldsymbol{P}}\left( {0.5} \right) + \mathit{\boldsymbol{Q}}\left( {0.5} \right)}}{2}. $

表明在t=0.5处, 曲线P(t)与Q(t)对过渡曲线R(t)有相同的影响.

性质3  对称性:fk(t)+fk(1-t)≡1.

证明  由式(10) 可得

$ {f_k}\left( {1 - t} \right) = \frac{{{t^{k + 1}} + \omega {t^{k + 1}}\left( {1 - t} \right)}}{{{t^{k + 1}} + \omega {t^{k + 1}}\left( {1 - t} \right) + \omega t{{\left( {1 - t} \right)}^{k + 1}} + {{\left( {1 - t} \right)}^{k + 1}}}}, $ (11)

故由式(10) 与式(11) 易得

$ {f_k}\left( t \right) + {f_k}\left( {1 - t} \right) \equiv 1. $

性质4  单调性:固定kω时, 势函数fk(t)关于变量t单调递减.

证明  令

$ \begin{array}{l} {g_k}\left( t \right) = {\left( {1 - t} \right)^{k + 1}}\left( {1 + \omega t} \right),\\ {h_k}\left( t \right) = {t^{k + 1}}\left( {1 + \omega \left( {1 - t} \right)} \right), \end{array} $

固定kω时, 由式(10) 经计算可得

$ \begin{array}{*{20}{c}} {{{f'}_k}\left( t \right) = \frac{{{{g'}_k}\left( t \right){h_k}\left( t \right) - {g_k}\left( t \right){{h'}_k}\left( t \right)}}{{{{\left( {{g_k}\left( t \right) + {h_k}\left( t \right)} \right)}^2}}} = }\\ { - \frac{{{t^k}{{\left( {1 - t} \right)}^k}s\left( t \right)}}{{{{\left( {{g_k}\left( t \right) + {h_k}\left( t \right)} \right)}^2}}},} \end{array} $ (12)

其中,

$ s\left( t \right) = \left( {2\omega - k{\omega ^2}} \right){\left( {t - 0.5} \right)^2} + \frac{{\left( {\omega + 2} \right)\left( {2k + k\omega + 2} \right)}}{4}. $ (13)

k≥1, ω > 0时, 由式(13) 可得

$ \left\{ \begin{array}{l} s\left( 0 \right) = s\left( 1 \right) = \left( {1 + k} \right)\left( {1 + \omega } \right) > 0,\\ s\left( {0.5} \right) = \frac{{\left( {\omega + 2} \right)\left( {2k + k\omega + 2} \right)}}{4} > 0. \end{array} \right. $ (14)

由式(14) 可知, 当0≤t≤1时, 有s(t) > 0, 于是由式(12) 可得fk(t)≤0, 故势函数fk(t)关于变量t单调递减.当ω=2时, 势函数f1(t), f2(t),f3(t)如图 2所示.

图 2 ω=2时的势函数fk(t)(k=1, 2, 3) Fig. 2 The potential function fk(t)(k=1, 2, 3) with ω=2

注3  由性质4及式(1) 可知, 随着t(0≤t≤1) 的增大, 曲线P(t)对过渡曲线R(t)的作用逐渐减小, 而曲线Q(t)对过渡曲线R(t)的作用则逐渐增大.

性质5  形状可调性:固定k, 当0≤t≤0.5时, 势函数fk(t)关于参数ω单调递减; 当0.5 < t≤1时, 势函数fk(t)关于参数ω单调递增.

证明  令

$ \begin{array}{l} {g_k}\left( t \right) = {\left( {1 - t} \right)^{k + 1}}\left( {1 + \omega t} \right),\\ {h_k}\left( t \right) = {t^{k + 1}}\left( {1 + \omega \left( {1 - t} \right)} \right), \end{array} $

固定kt, 由式(10),计算得

$ \frac{{{\rm{d}}{f_k}\left( t \right)}}{{{\rm{d}}\omega }} = \frac{{{t^k}{{\left( {1 - t} \right)}^k}\left( {2t - 1} \right)}}{{{{\left( {{g_k}\left( t \right) + {h_k}\left( t \right)} \right)}^2}}}. $ (15)

当0≤t≤0.5时, 有$\frac{{{\rm{d}}{f_k}(t)}}{{{\rm{d}}\omega }}$ ≤0, 此时势函数fk(t)关于参数ω单调递减; 当0.5 < t≤1时, 有$\frac{{{\rm{d}}{f_k}(t)}}{{{\rm{d}}\omega }}$ ≥0, 此时势函数fk(t)关于参数ω单调递增.参数ω取不同值时势函数f1(t)如图 3所示.

图 3 参数ω取不同值时的势函数f1(t) Fig. 3 The potential function f1(t) withdifferent values of the parameter ω

注4  由性质5知, 可通过修改参数ω(ω > 0) 的值调整势函数fk(t)的形状.于是, 由式(1) 知, 当曲线P(t)与Q(t)保持不变时, 过渡曲线在两端点处满足拟Ck(k≥1) 连续的情形下,可通过修改参数ω的值对其形状进行调整.

综上所述, 本文所构造的有理势函数fk(t)不仅与文献[10-15]中的势函数或调配函数性质完全相同, 而且还具有以下特性:

(1) 可令过渡曲线在两端点处达到拟Ck(k≥1) 连续.

(2) 具有统一的数学模型.

相对于文献[10-15]给出的势函数或调配函数, 本文构造的有理势函数fk(t)具有如下优势:

(1) 文献[10]所采用的势函数虽然结构简单, 但所构造的过渡曲线在端点处仅满足拟C1连续, 且无法调整过渡曲线的形状.而有理势函数fk(t)不仅可令过渡曲线在端点处达到拟Ck连续, 而且过渡曲线的形状还可通过所带的形状参数ω进行调整.

(2) 文献[11]所构造的2k+1次多项式势函数虽然可令过渡曲线在端点处达到拟Ck连续, 但无法调节过渡曲线的形状, 且满足不同连续性要求的势函数都需要通过求解方程组重新获得.而有理势函数fk(t)不仅具有统一的数学模型, 而且在令过渡曲线在端点处达到拟Ck连续的情形下仍可通过修改形状参数ω的值实现对过渡曲线形状的调整.

(3) 文献[11]所构造的几类三角混合势函数虽然可令过渡曲线具有形状可调性, 但过渡曲线在两端点处仅满足拟C1连续.而有理势函数fk(t)不仅可实现过渡曲线的形状可调, 而且还令过渡曲线在两端点处达到拟Ck连续.

(4) 虽然文献[12-14]所构造的调配函数可令过渡曲线同时具有形状可调性和在两端点处满足拟Ck(k=0, 1, 2) 连续性, 但都是将两被过渡曲线取为特定的带形状参数的曲线.而利用有理势函数fk(t)构造过渡曲线时, 不仅实现了形状可调和拟Ck连续, 而且两被过渡曲线可取为任意参数曲线.

(5) 文献[15]所构造的多项式势函数,不仅对两被过渡曲线的种类没有限制, 而且还可通过势函数中所带的参数调整过渡曲线的形状, 但过渡曲线在两端点处仅满足拟C2连续.而有理势函数fk(t)不仅对两被过渡曲线的种类没有限制, 而且在令过渡曲线在两端点处达到拟Ck(k≥1) 连续的情形下还可通过势函数所带的参数对其形状进行调整.

3 应用实例

给出以下几种情形下, 利用有理势函数fk(t)构造基于Metaball过渡曲线的实例.

(1) 直线与直线间的过渡曲线

设两直线P(t)与Q(t)的方程分别为

$ \left\{ \begin{array}{l} x\left( t \right) = 2\left( {1 - t} \right),\\ y\left( t \right) = 1 - t, \end{array} \right.\;\;\;\;\left\{ \begin{array}{l} x\left( t \right) = t,\\ y\left( t \right) = 2t, \end{array} \right.0 \le t \le 1. $

形状参数ω取不同值时, 分别利用有理势函数fk(t)(k=1, 2, 3) 构造满足拟Ck(k=1, 2, 3) 连续的不同形状的过渡曲线,如图 4所示.

图 4 直线与直线间的过渡曲线 Fig. 4 The transition curve between straight lines

(2) 直线与圆弧间的过渡曲线

设直线P(t)与圆弧Q(t)的方程分别为

$ \left\{ \begin{array}{l} x\left( t \right) = - \left( {1 - t} \right),\\ y\left( t \right) = 1 - t, \end{array} \right.\;\;\;\;\left\{ \begin{array}{l} x\left( t \right) = \sin \frac{\pi }{2}t,\\ y\left( t \right) = 1 - \cos \frac{\pi }{2}t, \end{array} \right.\;\;\;0 \le t \le 1. $

形状参数ω取不同值时, 分别利用有理势函数fk(t)(k=1, 2, 3) 构造满足拟Ck(k=1, 2, 3) 连续的不同形状的过渡曲线,如图 5所示.

图 5 直线与圆弧间的过渡曲线 Fig. 5 The transition curve betweenstraight line and arc

(3) 圆弧与圆弧间的过渡曲线

设两圆弧P(t)与Q(t)的方程分别为

$ \left\{ \begin{array}{l} x\left( t \right) = - \cos \frac{\pi }{2}t,\\ y\left( t \right) = 1 - \sin \frac{\pi }{2}t, \end{array} \right.\;\;\;\;\left\{ \begin{array}{l} x\left( t \right) = 1 - \cos \frac{\pi }{5}t,\\ y\left( t \right) = \sin \frac{\pi }{5}t, \end{array} \right.\;\;\;0 \le t \le 1. $

形状参数ω取不同值时, 分别利用有理势函数fk(t)(k=1, 2, 3) 构造满足拟Ck(k=1, 2, 3) 连续的不同形状的过渡曲线,如图 6所示.

图 6 圆弧与圆弧间的过渡曲线 Fig. 6 The transition curve between arcs

(4) 直线与曲线间的过渡曲线

设直线P(t)与曲线Q(t)的方程分别为

$ \left\{ \begin{array}{l} x\left( t \right) = - \left( {1 - t} \right),\\ y\left( t \right) = 1 - t, \end{array} \right.\;\;\;\;\left\{ \begin{array}{l} x\left( t \right) = \frac{\pi }{2}t - \sin \frac{\pi }{2}t,\\ y\left( t \right) = 1 - \cos \frac{\pi }{2}t, \end{array} \right.\;\;\;0 \le t \le 1. $

形状参数ω取不同值时, 分别利用有理势函数fk(t)(k=1, 2, 3) 构造满足拟Ck(k=1, 2, 3) 连续的不同形状的过渡曲线,如图 7所示.

图 7 直线与曲线间的过渡曲线 Fig. 7 The transition curve betweenstraight line and curve

(5) 圆弧与曲线间的过渡曲线

设圆弧P(t)与曲线Q(t)的方程分别为

$ \left\{ \begin{array}{l} x\left( t \right) = - \cos \frac{\pi }{2}t,\\ y\left( t \right) = 1 - \sin \frac{\pi }{2}t, \end{array} \right.\;\;\;\;\left\{ \begin{array}{l} x\left( t \right) = \pi t - \sin \frac{\pi }{2}t,\\ y\left( t \right) = 1 - \cos \frac{\pi }{2}t, \end{array} \right.\;\;\;0 \le t \le 1. $

形状参数ω取不同值时, 分别利用有理势函数fk(t)(k=1, 2, 3) 构造满足拟Ck(k=1, 2, 3) 连续的不同形状的过渡曲线,如图 8所示.

图 8 圆弧与曲线间的过渡曲线 Fig. 8 The transition curve between arcand curve

(6) 曲线与曲线间的过渡曲线

设两曲线P(t)与Q(t)的方程分别为

$ \left\{ \begin{array}{l} x\left( t \right) = \frac{1}{2}\left( {1 - t} \right),\\ y\left( t \right) = 2{\left( {1 - t} \right)^2}, \end{array} \right.\;\;\;\;\;\left\{ \begin{array}{l} x\left( t \right) = \pi t - \sin \frac{\pi }{2}t,\\ y\left( t \right) = 1 - \cos \frac{\pi }{2}t, \end{array} \right.\;\;\;0 \le t \le 1. $

形状参数ω取不同值时, 分别利用有理势函数fk(t)(k=1, 2, 3) 构造满足拟Ck(k=1, 2, 3) 连续的不同形状的过渡曲线,如图 9所示.

图 9 曲线与曲线间的过渡曲线 Fig. 9 The transition curve between curves

上述实例表明, 利用本文提出的带形状参数的有理势函数构造基于Metaball的过渡曲线时, 可通过选取势函数的次数及形状参数值获得连续的、形状不同的过渡曲线, 过渡曲线形状自然, 且均可光滑地连接2条被过渡曲线, 能满足不同的实际应用需求.

4 结语

为了满足更高要求的过渡曲线设计, 本文巧妙地利用一种带形状参数的曲线模型,构造了一类带形状参数的有理势函数.所构造的势函数不仅具有统一的数学模型, 而且在利用该势函数构造基于Metaball的过渡曲线时, 可令过渡曲线达到拟Ck(k≥1) 连续, 还可通过修改形状参数值调整过渡曲线的形状.

参考文献
[1] YAN L L, LIANG Q F. An extension of the Bézier model[J]. Applied Mathematics and Computation, 2011, 218(6): 2863–2879. DOI:10.1016/j.amc.2011.08.030
[2] BASHIR U, ABBSA M, ALI J M. The G2 and C2 rational quadratic trigonometric Bézier curve with two shape parameters with applications[J]. Applied Mathematics and Computation, 2013, 219(20): 10183–10197. DOI:10.1016/j.amc.2013.03.110
[3] 王文涛, 汪国昭. 带形状参数的均匀B样条[J]. 计算机辅助设计与图形学学报, 2004, 16(6): 783–788.
WANG W T, WANG G Z. Uniform B-spline with shape parameter[J]. Journal of Computer-Aided Design & Computer Graphics, 2004, 16(6): 783–788.
[4] 左传桂, 汪国昭. 多形状参数的四阶均匀B样条曲线设计[J]. 浙江大学学报:理学版, 2007, 34(4): 401–404.
ZUO C G, WANG G Z. Curve design of multi-parameter uniform B-spline blending function of order four[J]. Journal of Zhejiang University: Science Edition, 2007, 34(4): 401–404.
[5] HAN X L, ZHU Y P. Curve construction based on five trigonometric blending functions[J]. BIT Numerical Mathematics, 2012, 52(4): 953–979. DOI:10.1007/s10543-012-0386-0
[6] 徐迎博, 喻德生. 带形状参数的二次三角Bézier曲线形状分析[J]. 浙江大学学报:理学版, 2013, 40(1): 35–41.
XU Y B, YU D S. Shape analysis of quadratic trigonometric polynomial Bézier curves with a shape parameter[J]. Journal of Zhejiang University: Science Edition, 2013, 40(1): 35–41.
[7] 王健立, 毕合春, 刘奎丰. 滚齿加工的齿轮齿根过渡曲线[J]. 机械传动, 2013, 37(1): 82–86.
WANG J L, BI H C, LIU K F. Tooth root fillet curve of gear by hobbling[J]. Journal of Mechanical Transmission, 2013, 37(1): 82–86.
[8] 郭蓓, 赵远扬, 李连生, 等. 旋叶式压缩机的气缸型线研究[J]. 西安交通大学学报, 2003, 37(3): 256–259.
GUO B, ZHAN Y Y, LI L S, et al. Research on cylinder profile of rotary vane compressor[J]. Journal of Xi'an Jiaotong University, 2003, 37(3): 256–259.
[9] 宋立权, 赵学科, 李智成, 等. 任意缸体旋叶式压缩机的叶片型线设计理论研究及应用[J]. 机械工程学报, 2011, 47(15): 143–148.
SONG L Q, ZHAO X K, LI Z C, et al. Study and application of the designing theory of the vane profile for any rotary vane compressor[J]. Journal of Mechanical Engineering, 2011, 47(15): 143–148.
[10] 李凌丰, 谭建荣, 赵海霞. 基于Metaball的过渡曲线[J]. 中国机械工程, 2005, 16(6): 483–486.
LI L F, TAN J R, ZHAO H X. Construction method of transition curve based on Metaball technique[J]. China Mechanical Engineering, 2005, 16(6): 483–486.
[11] 高晖, 寿华好. 势函数的构造及基于Metaball的过渡曲线[J]. 计算机辅助设计与图形学学报, 2015, 27(5): 900–906.
GAO H, SHOU H H. Construction of potential function and transition curve based on Metaball technique[J]. Journal of Computer-Aided Design & Computer Graphics, 2015, 27(5): 900–906.
[12] 刘华勇, 张大明, 李璐. 基于参数连续HCBézier-like曲线的过渡曲线的构造[J]. 纯粹数学与应用数学, 2011, 27(1): 69–74.
LIU H Y, ZHANG D M, LI L. Constructing of blending HC Bézier-like curve preserving continuity[J]. Pure and Applied Mathematics, 2011, 27(1): 69–74.
[13] 刘华勇, 张大明, 李璐, 等. 几何连续的Bézier-like曲线的形状调配[J]. 山东大学学报:理学版, 2012, 47(3): 51–55.
LIU H Y, ZHANG D M, LI L, et al. Shape blending Bézier-like curve with geometric continuity[J]. Journal of Shandong University: Natural Science, 2012, 47(3): 51–55.
[14] 刘华勇, 段小娟, 张大明, 等. 基于三角Bézier-like的过渡曲线构造[J]. 浙江大学学报:理学版, 2013, 40(1): 42–46.
LIU H Y, DUAN X J, ZHANG D M, et al. Constructing of blending trigonometric Bézier-like curve[J]. Journal of Zhejiang University: Science Edition, 2013, 40(1): 42–46.
[15] 李军成, 宋来忠, 刘成志. 带参数的多项式势函数与构造基于Metaball的过渡曲线[J]. 中国图象图形学报, 2016, 21(7): 893–900.
LI J C, SONG L Z, LIU C Z. The polynomial potential function with a parameter and transition curve based on Metaball technology[J]. Journal of Image and Graphics, 2016, 21(7): 893–900. DOI:10.11834/jig.20160707
[16] 严兰兰, 韩旭里. 形状及光滑度可调的自动连续组合曲线曲面[J]. 计算机辅助设计与图形学学报, 2014, 26(10): 1654–1662.
YAN L L, HAN X L. Automatic continuous composite curve and surface with adjustable shape and smoothness[J]. Journal of Computer-Aided Design & Computer Graphics, 2014, 26(10): 1654–1662.