過渡曲面生成算法探究

韓松 吳麗娟 李博

摘 要：本文主要介紹了幾種過渡曲面生成算法的特點、應用范圍、優(yōu)缺點，包括滾球法、偏微分方程法、能量構造法、蒙皮構造法和基于裁剪線的過渡曲面構造法，以期為相關學者的研究提供借鑒。

關鍵詞：過渡曲面;滾球法;偏微分方程法

中圖分類號：TP391 文獻標識碼：A 文章編號：1003-5168（2019）02-0044-03

Research on Surface Smooth Mosaic Method

Absrtact： This paper mainly introduced the characteristics， application scope， advantages and disadvantages of several transition surface generation algorithms， including ball rolling method， partial differential equation method， energy construction method， skin construction method and transition surface construction method based on clipping line， in order to provide reference for relevant scholars.

Keywords：transition surface;rolling ball method;partial differential equation method

典型的過渡曲面生成算法主要有滾球法、偏微分方程法、能量構造法、蒙皮構造法和基于裁剪線的過渡曲面構造法。由此，本文主要對過渡曲面生成算法的特點、應用范圍、優(yōu)缺點進行分析[1]。

1 滾球法

1984年，Rossignac和Requicha提出滾球法。滾球過渡是兩片參數曲面間構造過渡曲面的經典方法，這樣所生成的過渡曲面和原曲面間是Gl連續(xù)的[2]。該方法可以看作是用一個球，沿著兩個過渡基曲面滾過，球滾過所形成的包絡線就是過渡曲面，滾球中心所走過的曲線被稱為脊線，而滾球與原基曲面相切的地方在滾動過程中所走過的曲線就是過渡切觸線，滾球法后來發(fā)展出了常半徑滾球過渡法和變半徑滾球過渡法[3]。

常半徑滾球過渡法的特點是球在滾動時，球的半徑不發(fā)生改變，形成的過渡曲面截線都具有相同的曲率。用這種方法來構造過渡曲面，首先要計算兩基曲面與過渡面的偏移半徑，設偏移半徑等于球的半徑，然后計算待拼接面與偏移面間的切線。切線就是滾球的中心線，也是過渡曲面的脊線。令滾球與兩基曲面相交的點為[C0和C1]。基曲面[S0（u0，v0）]在[C0]處的切平面為[Q0]，基曲面[S1（u1，v1）]在[C1]處的切平面為[Q1]，球心為O，則過[O，C0，C1]點有平面[Q2]，那么[Q0，Q1，Q2]相交于[C2]點，再以[C0，C1，C2]為控制定點來構造如下的二次截線[4]：

[Kα=ω0（1-α）2C0+2ω11-ααC2+ω2（α）2C1ω0（1-α）2+2ω11-αα+ω2（α）2] ? ? ? ? ? （1）

變半徑滾球過渡法的特點是滾球在滾動過程中，滾球的大小是變化的，但用變半徑滾球過渡法構造過渡曲面時，其脊線[P（u）]和半徑函數[F（u）]都是已知的，存在：

[Gu/Lu≤I] ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?（2）

構造過渡曲面時，首先根據已給定的初始脊線[LO（u）]和半徑[G（u）]，接下來調整初始脊線的位置使其與兩基曲面間等距，調整后的初始脊線稱為等距線[L（u）]，最后通過等距線來確定過渡曲面。

2 偏微分方程法

偏微分方程法（PDE）在幾何造型中的最初應用就是構造過渡面[5]。該方法具有以下優(yōu)點：曲面由其參數的超越函數表示而并不是簡單的多項式，因此生成的曲面自然光順;構造過渡曲面簡單易行，只需要給定邊界曲線和跨界導矢就可以生成光順的過渡曲面;除邊界曲線和跨界導矢外，也可調整方程中的一個物理參數來調整曲面的形狀。偏微分方程法受到了很多人的關注，又相繼發(fā)展出許多新的PDE求解過渡曲面方法，如Lihua.You先后提出的四階偏微分方程法和特征函數法。

3 能量構造法

能量構造法是在偏微分方程法基礎上發(fā)展出來的一種方法。1987年，加拿大學者Terzopoulos等將基于物理能量模型的可變形曲線曲面造型技術引入計算機圖形學中，使用物理能量模型來創(chuàng)建一個微分方程，引入約束處理和外載荷。此外，增加幾個補充條件，對這個方程進行一些變化，最后使用數值積分的方法來對這種物理能量模型所構造出的微分方程進行求解。之后，Celniker.G和Gossard.D又進一步優(yōu)化了這一方法，他們把能量模型設為目標函數，將特征線作為約束加到補充函數里去，這樣在加上外載荷以后就可以對所構造的曲面進行較好的調整和控制。

簡言之，能量構造法就是將曲面看成是彈性變形的薄殼，引入能量泛函，通過某些力學原理建立變形曲面的控制方程，然后引入一定的數值方法求解滿足一定約束條件的數值解[6]。把在待拼接面的過渡切觸線及其在待拼接面上的法矢和其他邊界線作為約束，采用基于物理的曲面造型技術來生成過渡曲面，再通過數值分析的方法來求解能量曲面的型值點，進而得到過渡曲面。

曲面能量泛函公式為：

[Esur=α11w2u+2α11wuwv+α11w2v+β11wuu+2β12wuv+β22wvv-2wfu，vdudv]（3）

4 蒙皮構造法

蒙皮構造法是在高級制造業(yè)中應用相對較為廣泛的拼接方法，其是通過一個拓撲點陣，或者離散的點云，或者截面曲線來構造一個光滑的曲面，這種方法在數學上就是通過在這些點陣和截面曲線間插值來構造過渡曲面，其可以簡單地想象成在一組曲面上蒙上一層光滑的曲面。這一方法最早是由Woodward提出的，他指出，蒙面就是在通過一組有序的截面曲線間插值來擬合一張曲面，這個過程就如在建造機翼時在骨架上加上一層蒙皮，這樣機翼蒙皮就稱為掃掠曲面[7-9]。

5 基于裁剪線的過渡曲面構造法

這種基于裁剪線的過渡曲面的構造法，生成的過渡曲面是不沿著基曲面的過渡切觸線與基曲面連接在一起，而是通過基曲面上剪裁線來進行連接的。這種方法較為靈活，但算法卻相當復雜。

在使用這種方法時，先在兩個基曲面確定兩條剪裁線，隨后依據這兩條剪裁線來構造脊線，最后計算出過渡面的截面曲線，從而得到想要的過渡面。在使用時，也可以跳過計算脊線的步驟直接構造過渡面。

基于剪裁線的過渡曲面構造法主要有Koparkar.P提出的方法和Filip.D.J提出的方法。由Koparkar.P提出的方法中，脊線是輔助曲面的交線，而這個輔助曲面可以是由基曲面的法線是沿裁剪線活動的軌跡生成的，隨后依據脊線和剪裁線可以計算截面曲線，最后用這條截面曲線掃掠即可得到過渡曲面[10-12]。

而由Filip.D.J提出的方法中，所生成的過渡面的構造函數為：

[Bs，t=H1sC1t+H2sC2t+H3sT1t+H4sT2t]（4）

式中，[Hjs，j=1，…，4]是三次埃爾米特插值多項式，[Cit，i=1，2]是裁剪線。方向矢量[Tit]為[C2t-C1t]在[Cit]基曲面的切平面上的投影。

6 曲面的光順

由于曲面拼接的廣泛應用，拼接之后對曲面的光順處理也就被重視起來，如何對其形狀進行調整就變得尤為重要。Hagen通過能量最小法對曲面進行光順處理。曲面的形狀由控制點決定，隨意對曲面進行光順處理時，只需要考慮控制頂點的位置[7]。由于測量方式和造型手段的局限性，所以，基于逆向工程重構出來的曲面的光順程度可能會不滿足需求，需要通過曲面光順來提高逆向工程曲面的品質。

Farin Sapidis提出的節(jié)點去除算法是一種實用性很強的光順算法，通過節(jié)點去除可以減少冗余的數據，并用于三次樣條曲線的局部光順，通過移動節(jié)點，提高曲線在節(jié)點處的連續(xù)率，從而達到光順。Pigounakis約束算法是全局算法中較有代表性的，給定某個約束條件，限定控制頂點的擾動量，結合此約束條件給出能量函數，通過能量函數最小化解得新的控制頂點。此算法能一次對整條曲線進行光順，但其缺點是需要求解優(yōu)化方法。

在實際應用中，存在只有個別控制點不光順的情況，如果使用全局光順，將會產生大量不必要的計算，而且會修改不需要被修改的控制點。由此，局部光順法應用而生，在每次的光順處理中只修改少數控制頂點，采用最優(yōu)化方法對其進行調整，而其他控制點不發(fā)生改變。

局部光順法有判別和修改兩部分。壞點的判別有兩種方式：①壞點由用戶決定，這稱為交互方式;②壞點由程序來決定，通過制定的光順標準，程序將自動決定壞點，這稱為自動方式。壞點的修改：當前主要的修改方法是Kjellander法。曲線不光順性的主要原因是三階導數的不連續(xù)性，可將控制點修改至一個更好的位置，使其三階導數差為0，進而達到光順。

7 結論

目前，曲面拼接在實際生產和生活中的應用越來越廣泛，曲面拼接的研究也如火如荼地進行著。因此，不同種類的拼接方式和拼接算法的需求不斷提高，選擇一個較為可行的曲面構造模型尤為關鍵。本文主要分析過渡曲面生成算法的特點、應用范圍、優(yōu)缺點，以期為學者選擇合適的曲面拼接方式提供借鑒。

參考文獻：

[1]陳浩乾.過渡曲面生成算法及其在隱形牙套模型中應用的研究[D].西安：西安科技大學，2008.

[2]陳亞娟.過渡面生成算法及曲面光順技術研究[D].北京：首都師范大學，2009.

[3]宋遒志，陳立平，鐘毅芳，等.基于物理的能量最小法構造融合過渡曲面[J].工程圖學學報，2002（3）：122-132.

[4]楊云志.過渡曲面構造方法[D].南京：南京航空航天大學，2013.

[5]陳亞娟，劉旭敏，關永.基于動態(tài)偏微分方程構造過渡面[J].計算機工程與設計，2009（14）：3352-3355.

[6]Celniker G，Gossard D. Deformable curve and surface finite-element for free-form shape design[J].Computer Graphics，1991（4）：257-266.

[7]賀晶.B樣條曲線曲面的光順[D].廣州：暨南大學，2008.

[8]華回春.曲線曲面分形逼近模型及應用[D].西安：西北工業(yè)大學，2004.

[9]王紀平.隱形牙套三位模型生成方法研究[D].西安：西安科技大學，2008.

[10]張曉軍.過渡曲面的B樣條曲線構造方法研究[D].大連：遼寧師范大學，2011.

[11]邵健萍，陳少良.CAD/CAM關鍵技術——曲面造型[J].機械工程師，2003（1）：23-25.

[12]朱永強，魯聰達.自由曲線曲面造型技術的綜述[J].中國制造業(yè)信息化，2003（5）：110-113.