針對不可展曲面的近似展開方法淺析

羅振華+蔣芬+蔣威

摘 要：隨著市場發展和人們生活質量的提升，現代工業產品造型越來越復雜多樣。復雜曲面大多通過平面材料加工而成，對不可展曲面的近似展開已成為工業設計的重要課題。該文對近20年來常用的不可展曲面的近似展開方法進行了分析、總結和歸納。

關鍵詞：不可展曲面 ?近似展開 ?市場

中圖分類號：T-9 文獻標識碼：A 文章編號：1674-098X（2014）10（c）-0051-02

Analysis of the Approximate Development Algorithms for Undevelopable Surfaces

LUO Zhenhua1 ?JIANG Fen1 ?JIANG Wei2

（1.Wuhan University of Technology Huaxia College， Wuhan Hubei， 430223;2.Shanghai High School， Shanghai， 200231，China）

Abstract：With the development of market and the improvement of people's living quality， modern industrial product design is more and more complicated. Most complex surfaces are processed by plane materials， the approximate development of undevelopable surfaces has become the important subject of industrial design. In this paper， the approximate development algorithms of undevelopable surfaces in recent twenty years are analyzed， summarized and concluded.

Key Words：Undevelopable Surface ?Approximate development ?Market

隨著市場發展和人們生活質量的提升，現代工業產品造型越來越復雜多樣，而在工業生產過程中，許多情形下曲面是通過對平面材料的加工而來的。為了節約原材料并提高生產效率，對不可展曲面的近似展開成為CAGD（Computer Aided Geometric Design）與工業制造中的重要課題。早在公元二世紀，希臘天文學家Ptolemy就提出了利用球面上的經線與緯線將球面展開到平面的方法。然而對于復雜曲面的展開問題一直沒能很好地解決，直到近20年來，隨著計算機技術的不斷進步，許多新的方法被學者們提出來。總的說來這些方法可以分為物理方法與幾何方法。

1 不可展曲面的物理展開方法

20世紀50年代，前蘇聯物理學家楚達列夫提出了滑移線法。滑移線是曲面上每點的兩個最大剪應力方向的包絡線。滑移線理論表明，處于屈服狀態的物體，塑性形變就沿著滑移線方向進行滑移。由于適用于滑移線方法的材料是要滿足一系列較為理想的假設條件的，因此，其應用范圍受到了這些較為嚴苛的假設條件的限制。到了20世紀80年代，比擬法應用較廣。這一方法的主要思路是，通過模擬材料內部介質在受到外力時的流動模型來探尋在實際加工中產生的形變。根據材料的不同特性衍生出了不同的模擬物理模型，常見的方法有電模擬法，熱傳導法，以及流體法。

電模擬法：對于不可壓縮的各向同性剛體，其塑性流動方程如下：

?2φ?x2+?2φ?y2=0 ? ? ? ?（1）

（1）式中φ表示流場速度的勢函數，這完全和電壓場中電壓U所應滿足的方程一樣：

?2U?x2+?2U?y2=0 ? ? ? ? （2）

于是可以通過實驗的方式用金屬絲模型來測量出電流的流動，進而模擬出變形后材料的外形。與之類似的熱傳導法針對在變形過程中滿足熱傳導方程的材料，而流體法則用來模擬易于流動的一類材料。從這里可以看出不僅材料的性質制約著比擬法的應用范圍，而且需要采用實驗的方法模擬材料的變形過程也不易于操作。有限元算法出現在20世紀60年代，目前被廣泛地應用到各個領域中。用有限元來展開曲面要先研究每一小塊物質單元的變形性質，來獲得一個合理的目標函數求解出最優值，進而得到所求的展開結果，C.H.Lee，H.Hu在1998年提出了一種簡單有效的目標函數：

minJp=en（he-h0）p1p ? ? ? （3）

（3）式中n表示節點數，he為單元的厚度，h0為初始的厚度，p為罰因子。

有限元算法是一種結果較為精確，并且可以反復連續地優化多次的方法。但是計算量十分巨大，花費時間較長。

2 不可展曲面的幾何展開方法

從幾何的角度講曲面展開就是去尋求如下的一一映射：φ：S→Ω ? ? ? （4）

這里S是一個三維曲面，Ω是平面上的一個區域。評判某種展開算法效果的優劣主要是計算曲面與其近似展開的平面區域之間的一些幾何量的變化程度，通常作為判斷標準的幾何量有邊長（弧長），角度，面積等。由于算法對象是不可展曲面，因此上述幾何量勢必會發生變化，以哪種幾何量作為判斷標準，需要根據實際需要來決定。近20年來通過幾何方法研究曲面展開主要有兩種算法：參數映射法和網格細分算法。endprint

2.1 參數映射法

參數映射法就是從某一類映射中去找合適的映射φ，在這類方法中線性參數映射是形式最簡潔，研究最深入，應用最廣泛的一類方法，具體的實現步驟可以歸納成以下幾步：第一步，將三角化曲面的邊界網格頂點，映射到平面的一個凸多邊形。第二步，在展開平面上，其內部網格點可以寫成其相鄰網格點的凸組合：

φv=w∈Nvμvwφ（w） ? ? ? ?（5）

Nv為所有與v相鄰的網格點的集合。記V1表示所有曲面內部網格點的集合，Vb表示所有曲面邊界網格點的集合，將式（5）中的內部頂點與邊界頂點進行分離，（5）式可以寫成如下形式：φv-w∈NV ∩V1μVWφw=w∈NV∩VbμVWφw （6）

將（6）寫成矩陣形式AX=B，這里X=φ（v）v∈V1為列向量，B是列向量，其元素由（6）式右端給出，A=awvw，v∈V1，且滿足：awv=1，w=v-μvw，w∈Nv0， otherwise求解這個線性方程組就可以得到φ的形式。不同的μvw的取值方式，會形成不同的展開結果。

2.2 網格細分算法

網格細分算法的發展時間不長，但是這種方法對曲面本身的形狀要求很低，因此被廣泛地應用到工業設計領域。其主要思路是將曲面細分成小的三角片或四邊形片區域，然后將這些小的三角片和四邊形片拼到同一個平面上，目前主要有以下兩種拼接方式。第一種方式是以曲面上某個網格點為中心，按照某種方式先拼接中心附近的三角片或四邊形片，再依次按算法規律向外拓展。龐林，劉玉君等提出了一種保形映射的方法，文章中將曲面近似展開時，允許每個三角片的三個內角發生變化而三角片的某兩條邊長度不變。依此方法是從展開中心開始，先橫向拼接，再縱向拼接，這樣的保形映射在展開的中心位置有比較好的吻合結果，而在離中心越遠者形狀會無法保持，同時在縱向拼接時會產生縫隙，康小明等在其文章中提出了另一種保面積的拼接方法，通過在拼接時邊長的伸縮保證映射前后的每個四邊形片的面積不發生變化，但與龐林等人的方法類似的也會遇到離中心越遠者形狀的偏差會很大。網格細分算法的研究大多是將細分后的曲面在一個平面上以某種方式拼接在一起，這類算法的特點在于簡潔方便，易于計算。其主要不可克服的缺陷就是過度的保證第一塊或前幾塊被映射到平面的曲面片的性質，同時這些方法也缺乏對近似展開后的圖形與原曲面的差別比較。

3 結語

目前，網格細分算法被廣泛應用與工業設計制造中。為保證其在展開過程中形變向量達到最小，展開結果與原曲面盡可能保持一致，還需要在網格細分算法的基礎上進一步優化。

參考文獻

[1] Lee C H，Huh H.Blank design and strain estimates for sheet metal forming pro-cesses by a finite element inverse approach with initial guess of linear deforma-tion[J].Journal of materials processing technology，1998，82（1）：145-155.

[2] Floater M，Hormann K.Parameterization of triangulations and unorganized points[M]//Tutorials on Multiresolution in Geometric Modelling.Springer Berlin Heidelberg，2002：287-316.

[3] 龐林，劉玉君，紀卓尚，等.不可展曲面近似展開的保形映射法[J].中國造船，2005，46（4）：58-64.

[4] 康小明，馬澤恩.不可展曲面近似展開的四邊形網格等面積法[J].西北工業大學學報，1998，16（3）：327-332.endprint