OpenCASCADE曲面交線分類

eryar@163.com

Abstract. OpenCascade classify the intersection line between two surfaces. A intersection line may be either geometric: line, circle, ellipse, parabola, hyperbola as defined in the class GLine, or analytic as defined in the class ALine, or defined by a set of points(coming from a walking algorithm) as defined in the class WLine. Or described by a restriction line on one of the surfaces as RLine.

Key Words. Surface Intersection, Intersection Line

1 Introduction

OpenCASCADE中對兩個曲面求交得到的交線進行了分類，如下類圖所示：

交線總共分為四類：

• ALine：Analytic解析曲線，主要為兩個二次曲面求交所得，如圓柱面、球面、圓錐面等之間的交線；
• GLine：Geometric幾何曲線，即交線可以表示成簡單的二次曲線，如直線、圓、拋物線等；
• WLine：Walking追蹤法得到的交線，保存了追蹤路線上的兩個曲面的交點；
• RLine：Restriction受限交線，這條交線可能只在一個面上；

下面我們使用Tcl腳本在DRAW中驗證一下這四類交線的來源，加深對曲面求交算法的理解。分類帶來了麻煩，可以帶著問題：為什么要分這幾種類型？有什么好處？來看這篇文章。要用好開源的東西，其實要求還是很高的，需要對源碼有相對深入的理解。

2 ALine

Analytic交線是二次曲面求交所得，二次曲面是因為可以統一使用二次型來表示的解析曲面，也是《解析幾何》中研究的主要內容。我們可以在DRAW中構造圓柱面與圓錐面求交驗證一下。TCL腳本如下：

# Test for IntPatch_ALine.
# Geometry surfaces.
cylinder s1 0 0 0 1 1 1 2
cone s2 0 0 0 0 0 1 1 0 0 30 3
# Topology faces
mkface f1 s1 0 2*pi -8 8
mkface f2 s2 0 2*pi -5 5
# Intersection.
bop f1 f2
bopsection r
# Display result.
vdisplay f1 f2 r

生成結果如下圖所示，其中紅色為交線：

DEBUG源碼可以看到是使用類IntPatch_ImpImpIntersection 計算求交，即兩個解析曲面求交算法類。

3 GLine

Geometric幾何曲線形式簡單，如果交線用幾何曲線來表示，對于后續算法有好處。如平面與圓錐面求交線，圓柱面與圓柱面求交等，都會得到幾何曲線。將上面的圓柱面換成平面與圓錐面求交我們可以在DRAW驗證經典的圓錐與平面交線：根據平面位置不同，可以得到圓、橢圓、雙曲線等幾何曲線。TCL腳本如下所示：

# Test for IntPatch_GLine.
cone s1 0 0 0 0 0 1 30 3
plane s2 0 0 0 0 0 1 
mkface f1 s1 0 2*pi -5 5
mkface f2 s2 -8 8 -8 8
bop f1 f2
bopsection r
vdisplay f1 f2 r

計算交線結果如下圖所示：

我們改變平面的法向，使其斜著與圓錐面求交，會得到橢圓：

還可以得到雙曲線、拋物線等，同學們可以自己嘗試一下。

4 WLine

對于NURBS曲面求交，一般會使用Marching方法，國內教材翻譯為追蹤法。在看《地球脈動》時，注意到對于大草原上的水牛、大象等動物成群結隊的遷徙使用了這個詞，這個詞的字面意思有行進、行軍，列隊行進之意，如果結合opencascade中的walking感覺翻譯為行進法更貼切，因為在opencascade中對于求交專門有個package名為IntWalk，其中類IntWalk_PWalking來使用marching method對兩個參數曲面進行求交。Walk有行走之意，所以對于使用Walk方法得到的交線命名為WLine。對于NURBS曲面求交及二次曲面與NURBS曲面求交，使用了Marching方法，WLine的來源是清晰的。繼續使用上面的腳本，只需要將上述兩個面轉換成NURBS曲面即可觸發Marching法進行求交。TCL腳本如下：

# Test for IntPatch_WLine.
cone s1 0 0 0 0 0 1 30 3
plane s2 0 0 0 1 1 2 
mkface f1 s1 0 2*pi -5 5
mkface f2 s2 -8 8 -8 8
nurbsconvert f1 f1
nurbsconvert f2 f2
bop f1 f2
bopsection r
vdisplay f1 f2 r

通過IntWalk_PWalking行進法配合三參數迭代法，將行進過程中的交點都保存在WLine中。

雖然結果與上面看上去一樣，內部交線已經不是簡單的幾何曲線了。

5 RLine

Restriction交線是受限交線，這種類型的交線只會位于一個面上。這里我們構造一個平面及與平面重疊的一個NURBS曲面來求交進行解釋。TCL腳本如下：

# Test for IntPatch_RLine.
plane s1 0 0 0 0 0 1
plane s2 0 0 0 0 0 1
mkface f1 s1 -5 5 -5 5
mkface f2 s2 -8 8 -8 8
nurbsconvert f1 f1
bop f1 f2
bopsection r
vdisplay f1 f2 r

DEBUG會發現這兩個曲面的交線為RLine，并在生成RLine時指定交線屬于哪個曲面，是在S1曲面SetArcOnS1還是在S2曲面SetArcOnS2：

生成交線如下圖所示：

當然可以使用RLine來判斷兩個曲面是否有重疊，但是在opencascade中兩個曲面重疊叫Tangent Face，可以將上述NURBS面不轉換，還是使用兩個重疊平面來驗證：

求交最后對交線數據歸并時的代碼有點不敢恭維：

6 Conclusion

綜上所述，對兩個曲面求交得到的交線進行分類，避免交線都是NURBS曲線，可以是簡單的二次曲線，提高后續算法性能。在理解源碼的基礎上，可以根據實際應用場景選擇高效的算法。如若只是求兩個模型之間的交線，可以直接使用曲面求交算法，一般情況下性能還是不錯的。當然理解源碼后，可以結合實際應用場景可以對求交算法做進一步優化。

要深入理解opencascade源碼，熟練使用DRAW是一個相對容易的路線。因為在DRAW中可以使用Tcl腳本快速驗證各種想法，甚至直接DEBUG源碼，從表向深入到與源碼作者直接對話。