周春凱，董楓，羅亮，陳順懷，廖顯庭

1中國(guó)艦船研究設(shè)計(jì)中心，湖北武漢4300642武漢理工大學(xué)交通學(xué)院，湖北武漢430063

基于三維CAD模型的復(fù)雜艏、艉船型曲面快速光順?lè)椒皩?shí)現(xiàn)

周春凱1，董楓2，羅亮2，陳順懷2，廖顯庭1

1中國(guó)艦船研究設(shè)計(jì)中心，湖北武漢430064
2武漢理工大學(xué)交通學(xué)院，湖北武漢430063

利用船舶型值和復(fù)雜艏、艉特征參數(shù)等數(shù)據(jù)，結(jié)合三維CAD設(shè)計(jì)軟件進(jìn)行初步建模，提出一種基于三維CAD軟件模型的復(fù)雜艏、艉船型曲面快速光順?lè)椒?，?shí)現(xiàn)復(fù)雜艏、艉船型曲面的快速三向光順。利用生成的光順?lè)抡婺Ｐ秃颓媲蠼凰惴?，獲得艏、艉部位任意型值點(diǎn)數(shù)據(jù)。該方法將拼接和光順處理分為三個(gè)階段來(lái)實(shí)現(xiàn)，即曲面重建、G1光順拼接、曲面微調(diào)3個(gè)階段。相較于三維CAD軟件中基于單個(gè)控制點(diǎn)調(diào)整的曲面拼接光順?lè)椒ǎ摲椒軌虼蠓岣咔嫫唇优c光順的處理效率，同時(shí)曲面的內(nèi)在連續(xù)性和型值點(diǎn)數(shù)據(jù)精度也得到了保證。

三維CAD模型；復(fù)雜艏、艉；船型曲面；拼接；光順

網(wǎng)絡(luò)出版地址：http://www.cnki.net/kcms/detail/42.1755.TJ.20160729.0945.004.html期刊網(wǎng)址：www.ship-research.com

引用格式：周春凱，董楓，羅亮，等.基于三維CAD模型的復(fù)雜艏、艉船型曲面快速光順?lè)椒皩?shí)現(xiàn)［J］.中國(guó)艦船研究，2016，11（4）：7-13.

ZHOU Chunkai，DONG Feng，LUO Liang，et al.Fast fairing method for ship-hull surfaces with complex bow and stern based on 3D models［J］.Chinese Journal of Ship Research，2016，11（4）：7-13.

0　引　言

目前，在三維CAD設(shè)計(jì)軟件中進(jìn)行的船舶型線(xiàn)設(shè)計(jì)，主要是通過(guò)將船舶各分段對(duì)應(yīng)的型值點(diǎn)轉(zhuǎn)化為全船坐標(biāo)點(diǎn)，再將全船坐標(biāo)點(diǎn)生成樣條曲線(xiàn)，然后由樣條曲線(xiàn)生成NURBS曲面模型，由此進(jìn)行船舶型線(xiàn)設(shè)計(jì)和總體性能計(jì)算。然而，在型值表中的這些離散型值點(diǎn)并不能夠完全有效地表達(dá)整個(gè)船體線(xiàn)型。型值點(diǎn)缺乏數(shù)據(jù)不僅在三維設(shè)計(jì)軟件建模過(guò)程中容易引起曲面片之間的較大孔洞和縫隙，使曲面拼接工作量大為增加，而且造成復(fù)雜艏、艉形狀與船舶主體的連接往往不光順，進(jìn)一步加大了曲線(xiàn)、曲面光順的工作量。同時(shí)，在船舶CAD/CAE軟件的三維建模過(guò)程中，對(duì)于復(fù)雜艏、艉船型曲面，其曲面曲率變化劇烈，僅根據(jù)型線(xiàn)圖中提供的有限站位處的型值點(diǎn)無(wú)法建立光順和精確的復(fù)雜艏、艉三維模型。三維模型的不光順不僅造成三維CAD設(shè)計(jì)軟件自動(dòng)生成的二維型線(xiàn)圖不光順，也容易引起總體性能計(jì)算結(jié)果的偏差。

另一方面，在船舶復(fù)雜艏、艉船型多方案優(yōu)化設(shè)計(jì)中，設(shè)計(jì)人員需要不斷將光順后的復(fù)雜艏、艉船型數(shù)據(jù)導(dǎo)入到CAD/CAE三維設(shè)計(jì)軟件中進(jìn)行艏、艉選型設(shè)計(jì)和優(yōu)化。在三維設(shè)計(jì)軟件中建立光順的三維模型時(shí)，由于需要通過(guò)手工交互調(diào)節(jié)曲線(xiàn)曲面的控制點(diǎn)不斷地進(jìn)行曲線(xiàn)曲面的拼接和光順才能實(shí)現(xiàn)，調(diào)節(jié)過(guò)程異常繁雜，工作量大，計(jì)算和分析效率低，嚴(yán)重制約了船舶型線(xiàn)設(shè)計(jì)和復(fù)雜艏、艉多方案選型與優(yōu)化的效率。

有鑒于此，本文利用船舶型線(xiàn)設(shè)計(jì)所必需的型值和復(fù)雜艏、艉特征參數(shù)等相關(guān)數(shù)據(jù)，結(jié)合三維設(shè)計(jì)軟件進(jìn)行初步建模和型線(xiàn)設(shè)計(jì)。通過(guò)導(dǎo)出IGES文件提取三維模型數(shù)據(jù)，并采用曲線(xiàn)曲面拼接與光順?biāo)惴ǎ?-2］等處理技術(shù)，避免頻繁地手工交互調(diào)整控制點(diǎn)，自動(dòng)快速實(shí)現(xiàn)船舶復(fù)雜艏、艉三維模型的光順。利用生成的三向光順?lè)抡婺Ｐ?，采用曲面求交算法獲得船體肋位以及復(fù)雜艏、艉部位的任意橫剖面處型值點(diǎn)的數(shù)據(jù)，而且，通過(guò)相應(yīng)的IGES數(shù)據(jù)接口的研究開(kāi)發(fā)，可以快速將生成的船舶和復(fù)雜艏、艉外型三維光順?lè)抡婺Ｐ图靶椭禂?shù)據(jù)導(dǎo)入到其他CAD/CAE軟件中以進(jìn)一步完成模型處理和仿真分析，大幅提高復(fù)雜艏、艉船型三維建模的效率，減少大量、繁瑣的數(shù)據(jù)重復(fù)輸入，從而快速有效地進(jìn)行復(fù)雜艏、艉的設(shè)計(jì)和優(yōu)化，為復(fù)雜艏、艉的型線(xiàn)設(shè)計(jì)研究和總體性能預(yù)估等提供技術(shù)支撐。

1　曲面拼接光順?biāo)惴ɡ碚摲治?/h2>

1.1曲面拼接理論及算法原理

船體曲面特別是艏、艉曲面由輪廓復(fù)雜和不規(guī)則的三維曲面構(gòu)成，在三維CAD設(shè)計(jì)軟件系統(tǒng)中，試圖用一張曲面去構(gòu)造其幾何模型極難保證模型精度，因此，需要將船體曲面劃分為多個(gè)比較簡(jiǎn)單、容易處理的曲面片，根據(jù)實(shí)際的光滑性要求對(duì)這些曲面片進(jìn)行拼接和縫合，最終完成船體曲面設(shè)計(jì)，實(shí)現(xiàn)船體的整體幾何造型。曲面拼接和光順是曲面設(shè)計(jì)的兩項(xiàng)基本工作，兩者具有理論交叉，拼接的本質(zhì)是使鄰接的曲面片在邊界處具有滿(mǎn)足設(shè)計(jì)要求的光順性。衡量拼接質(zhì)量的一個(gè)重要因素就是光滑性，而光滑性的高低取決于拼接曲面邊界的連續(xù)性。

兩塊B樣條曲面的G1拼接算法［3-4］適用于邊界完全重疊的情況，也適用于不完全重疊的情況，如圖1所示。給定兩B樣條曲面 p（s,t）和q（u,v），其控制頂點(diǎn)網(wǎng)格分別為（i=0,1,…, m1；j=0,1,…,n1），（i=0,1,…,m2；j=0,1,…, n2）；參數(shù)s與t的次數(shù)為 k1與l1，兩個(gè)節(jié)點(diǎn)矢量與；參數(shù)u與v的次數(shù)為k2與l2，兩個(gè)節(jié)點(diǎn)矢量與。對(duì)于G0連續(xù)的兩曲面 p（s,t）與q（u,v）上的公共邊界點(diǎn)Qj，有Qj=p（s′,0）=q（u′,1）。s′和u′分別為節(jié)點(diǎn)矢量S和U中對(duì)應(yīng)Qj的參數(shù)值，pt（s′,0），qv（u′,1），qu（u′,1）分別為點(diǎn)Qj處跨界切向矢量。兩曲面在Qj處具有公共的切平面的條件如式（1）所示。

式中：a與b為常數(shù)。另外為保證兩曲面在該點(diǎn)處不形成尖棱，應(yīng)滿(mǎn)足a＞0。

圖1　B樣條曲面的G1拼接Fig.1　G1-splicing of B-spline surface

兩曲面在點(diǎn)Qj處的等參數(shù)線(xiàn)可表示為

式中，Nj,l（t）為B樣條基函數(shù)，其余類(lèi)似符號(hào)與此

1意義相同。根據(jù)B樣條曲線(xiàn)的端點(diǎn)性質(zhì)，有Qj==成立，且兩曲面在Qj處的跨界切向矢量分別為

求解式（6）可以使用拉格朗日乘子法，得到a 與b，再代入式（4），便可得到調(diào)整后的控制頂點(diǎn)dˉ1j1。對(duì)于B樣條曲面p（s,t）上等參數(shù)線(xiàn)p（s′,t）的第i0個(gè)控制頂點(diǎn)，見(jiàn)式（2），當(dāng)s′從0連續(xù)變化到1時(shí)所形成的曲線(xiàn)C（s）就是由曲面 p（s,t）第i0列控制頂點(diǎn)定義的k1次B樣條曲線(xiàn)，數(shù)學(xué)表達(dá)式為

兩個(gè)B樣條曲面間G1光滑拼接的算法描述如下：

1）以預(yù)先給定的份數(shù)N分割參數(shù)s的定義區(qū)間，獲得節(jié)點(diǎn)序列。

2）循環(huán)i=0,1,…,N，由式（7）獲得序列｛Sˉ｝中Sˉ處t參數(shù)線(xiàn)的第1個(gè)（以0為起點(diǎn)）控制頂點(diǎn)。同時(shí)判斷該位置是否位于兩邊界線(xiàn)的重合范圍，由此可由式（4）得到調(diào)整后的。

重復(fù)上述步驟，可使 p（s,t）與q（u,v）達(dá)到近似的G1連續(xù)。

如果想進(jìn)一步簡(jiǎn)單處理，即避開(kāi)式（6）定義的優(yōu)化問(wèn)題，可在步驟（2）中取a=1，b=0。

1.2曲線(xiàn)光順理論及算法原理

光順準(zhǔn)則為曲率變化均勻，即曲線(xiàn)具有光滑的曲率梳。

由于關(guān)注的是均勻的曲率變化，因此這里將曲率k和曲率導(dǎo)數(shù)k′=dk/ds作為判斷的光順依據(jù)，因此曲線(xiàn)C（t）的光順性泛函為式中：第1項(xiàng)表征應(yīng)變能，控制曲線(xiàn)的總曲率，第2項(xiàng)表征曲率導(dǎo)數(shù)，控制曲線(xiàn)的曲率變化；α和β為平衡權(quán)因子，α+β=1，控制光順的側(cè)重點(diǎn)；C′（t）和C″（t）分別為曲線(xiàn)C（t）的一階和二階導(dǎo)數(shù)。

光順?biāo)惴ǎ?-9］以式（8）為目標(biāo)函數(shù)，約束條件為設(shè)計(jì)者指定的控制頂點(diǎn)或型值點(diǎn)偏離，包括切向矢量和法向矢量等約束條件，令光順前后的曲線(xiàn)分別為C0（t）和C（t），則約束條件相應(yīng)的表達(dá)式分別為：

針對(duì)正算問(wèn)題，控制頂點(diǎn)偏離：

針對(duì)插值問(wèn)題，型值點(diǎn)偏離：

切向矢量約束：

法向矢量約束：

在構(gòu)造了光順性泛函和約束條件后，便可通過(guò)相應(yīng)的優(yōu)化問(wèn)題得到光順優(yōu)化后的曲線(xiàn)。

2　曲面快速拼接光順?biāo)惴▽?shí)現(xiàn)

2.1船型曲面拼接光順整體研究流程

利用船舶型值和復(fù)雜艏、艉特征參數(shù)等數(shù)據(jù)，結(jié)合三維CAD設(shè)計(jì)軟件進(jìn)行初步建模和型線(xiàn)設(shè)計(jì)，采用曲面拼接和曲線(xiàn)曲面光順?biāo)惴?，進(jìn)行復(fù)雜艏、艉船型的快速三向光順。利用生成的光順模型，采用曲面求交算法可以獲得復(fù)雜艏、艉曲面任意型值點(diǎn)的數(shù)據(jù)，相應(yīng)的數(shù)據(jù)接口可以快速將復(fù)雜艏、艉船型三維光順模型數(shù)據(jù)轉(zhuǎn)換成.iges文件格式，導(dǎo)入到CAD/CAE軟件中進(jìn)行進(jìn)一步的模型處理和仿真分析，從而提高復(fù)雜艏、艉船型的優(yōu)化設(shè)計(jì)效率。整體研究流程如圖2所示。

圖2　整體研究流程Fig.2　The process of whole research for fairing

2.2曲面快速拼接、光順的具體實(shí)施方案

曲面的快速拼接光順處理分為3個(gè)階段，在這3個(gè)階段之前稱(chēng)為預(yù)處理階段，即還需要經(jīng)過(guò)曲線(xiàn)的光順處理，主要是對(duì)一些邊界曲線(xiàn)進(jìn)行快速光順處理，參見(jiàn)1.2節(jié)。

第1階段：曲面重新生成。使具有共同邊界的各曲面在拼接處具有相同數(shù)目的控制頂點(diǎn)，同時(shí)在相應(yīng)的方向具有相同的次數(shù)，這一階段因?yàn)樗械那娼詾殡p三次所以可自動(dòng)滿(mǎn)足。在此階段可以通過(guò)三維設(shè)計(jì)軟件Foran中的reinterpolate命令或者用Rhino的曲面重建命令實(shí)現(xiàn)，其他如CAD/CAM軟件也基本包含類(lèi)似功能。在處理過(guò)程中曲面形狀可能會(huì)發(fā)生變化，但是只要給定適當(dāng)數(shù)量的u向和v向的控制頂點(diǎn)數(shù)，便可以較好地達(dá)到保形的要求。

第2階段：曲面G1光滑拼接。這部分內(nèi)容參見(jiàn)第1.1節(jié)的拼接算法實(shí)現(xiàn)過(guò)程。

第3階段：各曲面內(nèi)部調(diào)整。保持邊界和G1連續(xù)性條件。

船型曲面需要滿(mǎn)足船舶的設(shè)計(jì)要求，本文采用一種保持邊界和G1連續(xù)性的曲面全局調(diào)整方法來(lái)調(diào)整曲面，讓設(shè)計(jì)者交互地調(diào)整曲面，從而最終獲得較為光順的曲面。該方法雖然依舊為交互式調(diào)整，但由于每次能大范圍調(diào)整，大幅減少了設(shè)計(jì)工作量，極大提高了復(fù)雜曲面的拼接光順效率。

具體的調(diào)整方法分為兩種：一種為雙線(xiàn)性調(diào)整，另一種為雙余弦調(diào)整。兩者流程基本一致，具體步驟如下：

1）獲得曲面的u向和v向控制頂點(diǎn)數(shù)m和n，如圖3所示。圖中的方格點(diǎn)即為曲面的控制網(wǎng)格。注意這里只調(diào)節(jié)曲面控制頂點(diǎn)的y坐標(biāo)分量，圖中為控制頂點(diǎn)的拓?fù)浣Y(jié)構(gòu)圖，而非實(shí)際的控制頂點(diǎn)。

2）設(shè)計(jì)者根據(jù)曲面具體的凹凸情況選擇調(diào)整峰值點(diǎn)，即調(diào)整幅度最大的點(diǎn)，如圖3中的紅色點(diǎn)所示，調(diào)整的幅值r也由設(shè)計(jì)者給定。具體的幅值可參考該曲面和周?chē)徑忧娴钠ヅ淝闆r。以艏部曲面為例，對(duì)于拓?fù)錇?×12的曲面，由于其相對(duì)周?chē)徑忧嫱蛊鹈黠@，因此幅值可取得較大一些，例如0.1或0.2。通常幅值建議取小一些，這樣能夠更容易掌握曲面的整個(gè)細(xì)節(jié)變化過(guò)程。

圖3　單片曲面分塊示意圖Fig.3　Division of single patch

3）為了保持邊界和與其他曲面間的G1連續(xù)性條件，需要保證艏、艉兩行和艏、艉兩列的控制頂點(diǎn)不變，需調(diào)整的點(diǎn)僅包括中間的藍(lán)色方格點(diǎn)和紅色峰值點(diǎn)。

4）以峰值點(diǎn)為坐標(biāo)原點(diǎn)，將需調(diào)整的點(diǎn)分為4個(gè)分塊，每個(gè)分塊采用相應(yīng)的雙線(xiàn)性函數(shù)或雙余弦非線(xiàn)性函數(shù)進(jìn)行調(diào)整。函數(shù)的確定方式要保證在峰值點(diǎn)處函數(shù)取得極大值，同時(shí)保證函數(shù)在緊鄰最外層的那一層控制頂點(diǎn)處取得極小值，即調(diào)整值為0。以第1分塊為例，雙線(xiàn)性函數(shù)和雙余弦函數(shù)分別為：

式中：d1為峰值點(diǎn)到該分塊與邊界控制頂點(diǎn)相鄰的那一列控制頂點(diǎn)間的x向距離；d2為峰值點(diǎn)到該分塊與邊界控制頂點(diǎn)相鄰的那一行控制頂點(diǎn)間的z向距離。其他3個(gè)分塊的調(diào)整函數(shù)確定方式與第1分塊類(lèi)似。

5）輸出調(diào)整后的曲面。第3階段的流程如圖4所示。

圖4　曲面微調(diào)流程Fig.4　The process of surface fine-tuning

3　算法測(cè)試與數(shù)據(jù)對(duì)比分析

3.1第1階段前后對(duì)比

以某大型拖網(wǎng)漁船艏部曲面為例（艉部曲面的處理與此相同），分析3個(gè)階段的前后數(shù)據(jù)對(duì)比。第1階段為曲面重建，使共邊界曲面在邊界處具有相同數(shù)量的控制頂點(diǎn)。通過(guò)對(duì)比，重建后的曲面完全達(dá)到了保形要求，而且還消除了部分曲面褶皺，但獨(dú)立曲面片之間的不連續(xù)性依舊非常明顯，如圖5所示。

圖5　重建前、后曲面對(duì)比圖Fig.5　Comparison of surface before and after reinterpolation

3.2第2階段前后對(duì)比

G1拼接前的曲面即為重建后的曲面，如圖5（b）所示，圖6則為G1拼接后的曲面。從G1光滑拼接前后對(duì)比，可明顯地看到曲面在拼接處光順性的改善，且從圖7所示的等照度線(xiàn)圖對(duì)比中也可證明各曲面片確實(shí)達(dá)到了近似G1連續(xù)。

圖6　G1拼接后曲面圖Fig.6　Surface after G1-splicing

圖7　G1拼接前、后等照度線(xiàn)對(duì)比圖Fig.7　Comparison of surface's isophotes before and after G1-splicing

3.3第3階段前后對(duì)比

分析上一階段得到的艏部曲面，發(fā)現(xiàn)連接球艏和進(jìn)流段大曲面的曲面片在內(nèi)部有些不必要的凸起，因此調(diào)整該曲面，將行數(shù)為8，列數(shù)為4的控制頂點(diǎn)作為修改的峰值點(diǎn)，調(diào)整幅值為0.1 m，僅經(jīng)過(guò)幾次的雙線(xiàn)性和雙余弦調(diào)整后的曲面對(duì)比如圖8和圖9所示。從圖中可以明顯看出調(diào)整后的曲面大幅地削弱了該凸起面的凸起情況，調(diào)整后的曲面要明顯優(yōu)于調(diào)整前的曲面。由于可以實(shí)時(shí)顯示曲面的投影等值線(xiàn)，因此僅通過(guò)幾次調(diào)整便達(dá)到了所期望的較好效果，相對(duì)于之前的單點(diǎn)調(diào)整，大幅提高了光順修改效率。

3.4精度驗(yàn)證

圖8　微調(diào)前后主視圖對(duì)比Fig.8　Comparison of surface's front view before and after fine-tuning

圖9　微調(diào)前后俯視圖對(duì)比Fig.9　Comparison of surface's plan view before and after fine-tuning

以19.5站為例，分析拼接光順處理前后的精度偏差。表1列出了處理前和處理后的曲線(xiàn)與原插值曲線(xiàn)的型值點(diǎn)偏差值，型值點(diǎn)數(shù)據(jù)通過(guò)基于二維標(biāo)量場(chǎng)抽取等值線(xiàn)的平面與曲面求交算法［10-12］獲得，確保不會(huì)漏交。由于站線(xiàn)是橫向曲線(xiàn)，因此這里只列出了y方向的偏差值。從表中可以看出處理后要比處理前的偏差小得多。

表1　處理前后19.5站曲線(xiàn)與原曲線(xiàn)型值點(diǎn)偏差Tab.1　Relative errors of 19.5 station curve before and after splicing and fairing

3.5運(yùn)行時(shí)間

第1階段可以直接在三維CAD設(shè)計(jì)軟件上執(zhí)行，運(yùn)行時(shí)間較快，只要點(diǎn)擊相應(yīng)命令即可瞬間完成。第2階段通過(guò)執(zhí)行優(yōu)化函數(shù)來(lái)完成保形的G1拼接，由于需要處理大量的矩陣和優(yōu)化問(wèn)題，因此這個(gè)過(guò)程耗時(shí)最多，針對(duì)艏部曲面大概需要3～5 min。第3階段本質(zhì)上是一個(gè)人工交互過(guò)程，但是該過(guò)程通過(guò)用偏離函數(shù)來(lái)對(duì)整個(gè)控制頂點(diǎn)拓?fù)溥M(jìn)行變形，因此，相較于每次只調(diào)整一個(gè)控制點(diǎn)，會(huì)節(jié)約大量的時(shí)間，艏部實(shí)例第3階段耗時(shí)約需幾分鐘。顯然，這要比在三維CAD軟件中進(jìn)行單點(diǎn)光順調(diào)整的效率高得多。

4　復(fù)雜艏、艉船型快速光順系統(tǒng)

4.1曲線(xiàn)光順系統(tǒng)界面

復(fù)雜艏、艉船型曲面快速拼接光順系統(tǒng)的編寫(xiě)以MATLAB的圖形用戶(hù)界面GUI為平臺(tái)［13］。圖10所示為艏輪廓邊界曲線(xiàn)的光順實(shí)例界面圖，其中控制頂點(diǎn)容許偏差取為0.5。在該界面中，第1個(gè)坐標(biāo)系是用來(lái)表示導(dǎo)入和修改的曲線(xiàn)，第2個(gè)坐標(biāo)系是用來(lái)進(jìn)行人工交互的界面。該系統(tǒng)可以對(duì)曲線(xiàn)進(jìn)行快速光順，顯示曲線(xiàn)的控制頂點(diǎn)和曲率圖等功能。

圖10　艏輪廓邊界曲線(xiàn)光順界面Fig.10　The interface of fairing for bow's profile curve

4.2曲面拼接和微調(diào)系統(tǒng)界面

圖11所示為艏部曲面的光順拼接和微調(diào)界面，在該界面中，坐標(biāo)系是用來(lái)表示導(dǎo)入和修改的曲面，設(shè)計(jì)者可實(shí)時(shí)觀(guān)測(cè)修改后的曲面以及相應(yīng)的Z向投影等值線(xiàn)（即水線(xiàn)），以便更直接地判斷曲面的修改趨勢(shì)。該系統(tǒng)軟件首先能夠?qū)η孢M(jìn)行G1拼接光順，然后可以對(duì)曲面控制頂點(diǎn)進(jìn)行線(xiàn)性和非線(xiàn)性（余弦）的微調(diào)，同時(shí)能夠讀取曲面的控制網(wǎng)格，顯示在右側(cè)的3個(gè)表格中。如圖11所示，控制頂點(diǎn)偏差取0.1，調(diào)整曲面編號(hào)為5，即凸起最大的那片連接球艏和主船體之間的曲面，峰值點(diǎn)為8行、4列的控制頂點(diǎn)。

圖12所示為艉部曲面的光順微調(diào)界面，其中控制頂點(diǎn)偏差取0.1，調(diào)整曲面編號(hào)為4，即連接艉封板和主船體之間的曲面，峰值點(diǎn)為5行、30列的控制頂點(diǎn)。

圖11　艏部曲面調(diào)整實(shí)例Fig.11　The interface of fairing for bow's surface

圖12　艉部曲面調(diào)整實(shí)例Fig.12　The interface of fairing for stern's surface

5　結(jié)　語(yǔ)

本文深入研究了船體曲面的拼接和光順處理問(wèn)題，提出了一種基于三維CAD軟件模型的復(fù)雜艏、艉船型曲面快速光順?lè)椒ǎ瑢⑵唇雍凸忭樚幚矸譃?個(gè)階段來(lái)實(shí)施：曲面重建；G1光順拼接；曲面微調(diào)。相較于三維CAD軟件中傳統(tǒng)的單個(gè)控制點(diǎn)調(diào)整方法，本方法能夠極大地提高曲面光順修改效率，同時(shí)曲面的內(nèi)在連續(xù)性和型值點(diǎn)數(shù)據(jù)精度也得到了保證。

本方法的不足之處在于第3階段依舊是一個(gè)人工交互過(guò)程，盡管提高了效率，但還是比較依賴(lài)于設(shè)計(jì)者的實(shí)際經(jīng)驗(yàn)。如何進(jìn)一步實(shí)現(xiàn)復(fù)雜曲面光順處理的全自動(dòng)化，是今后研究的方向。

［1］ 程少華，吳華.兩相鄰3次B樣條曲面G1連續(xù)充分條件［J］.河南師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2006，34 （3）：48-50.

CHENG Shaohua，WU Hua.G1continuity between two patches of adjacent cubic B spline surface［J］.Journal of Henan Normal University（Natural Science Edition），2006，34（3）：48-50.

［2］CHE X J，LIANG X Z.G1continuity conditions of B-spline surfaces［J］.Northeastern Mathematical Journal，2002，18（4）：343-352.

［3］ 李旭.雙三次NURBS曲面拼接及光順的研究與實(shí)現(xiàn)［D］.西安：西安理工大學(xué)，2007：18-20.

［4］曲學(xué)軍，寧濤，席平.B樣條曲面的光滑拼接［J］.計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào)，2004，16（1）：138-141.

QU Xuejun，NING Tao，XI Ping.Smooth joining between adjacent B-spline surfaces［J］.Journal of Computer Aided Design&Computer Graphics，2004，16 （1）：138-141.

［5］朱心雄.自由曲線(xiàn)曲面造型技術(shù)［M］.北京：科學(xué)出版社，2000.

［6］穆國(guó)旺，涂侯杰，雷毅，等.參數(shù)三次B樣條曲線(xiàn)的一種整體光順?lè)椒ǎ跩］.工程圖學(xué)學(xué)報(bào)，1998（1）：28-34.

MU Guowang，TU Houjie，LEI Yi，et al.A global fairing method for cubic B-spline curves［J］.Journal of Engineering Graphics，1998（1）：28-34.

［7］ 穆國(guó)旺，宋秀琴，臧婷.一種選點(diǎn)法和能量法相結(jié)合的曲線(xiàn)光順?lè)椒ǎ跩］.工程圖學(xué)學(xué)報(bào)，2005（6）：118-121.

MU Guowang，SONG Xiuqin，ZANG Ting.A curve fairing method mixing bad points selection with energy minimization［J］.Journal of Engineering Graphics，2005（6）：118-121.

［8］KJELLANDER J A P.Smoothing of cubic parametric splines［J］.Computer-Aided Design，1983，15（3）：175-179.

［9］SAPIDIS N，F(xiàn)ARIN G.Automatic fairing algorithm for B-spline curves［J］.Computer-Aided Design，1990，22（2）：121-129.

［10］ 仵大偉，林焰，紀(jì)卓尚.船體曲面的NURBS表達(dá)與設(shè)計(jì)［J］.大連理工大學(xué)學(xué)報(bào)，2002，42（5）：569-573.

WU Dawei，LIN Yan，JI Zhuoshang.Ship hull surface's NURBS representation and design［J］.Journal of Dalian University of Technology，2002，42（5）：569-573.

［11］ 張明霞，紀(jì)卓尚，林焰.NURBS曲面與隱式曲面求交的計(jì)算機(jī)實(shí)現(xiàn)及應(yīng)用［J］.中國(guó)造船，2002，43 （3）：94-98.

ZHANG Mingxia，JI Zhuoshang，LIN Yan.Implementation and application of intersection between NURBS-based surface and implicit surface［J］.Shipbuilding of China，2002，43（3）：94-98.

［12］ 宋宏勛，韓毅，吳初娜.一種基于等值線(xiàn)法的NURBS曲面與平面的求交算法［J］.數(shù)字技術(shù)與應(yīng)用，2011（7）：103-105.

SONG Hongxun，HAN Yi，WU Chuna.An intersection algorithm of NURBS surface and plane based on the method of isoline［J］.Digital Technology&Application，2011（7）：103-105.

［13］ HANSELMAN D，LITTLEFIELD B.精通MATLAB綜合輔導(dǎo)與指南［M］.李人厚，譯.西安：西安交通大學(xué)出版社，1998.

Fast fairing method for ship-hull surfaces with complex bow and stern based on 3D models

ZHOU Chunkai1，DONG Feng2，LUO Liang2，CHEN Shunhuai2，LIAO Xianting11 China Ship Development and Design Center，Wuhan 430064，China

2 School of Transportation，Wuhan University of Technology，Wuhan 430063，China

Based on the point net of offset and complex ship-hull parameters，using preliminary surface models made by 3D CAD software，this paper presents a fast fairing method for ship-hull surface modelling with complex bow and stern based on 3D models.With fairing simulation 3D-models and the surface intersection algorithm，the fairing point net of offset for arbitrarily complex bow and stern of ship-hull surfaces is first obtained.The method then divides the splicing and faring process into three stages:surface reinterpolation，G1-splicing，and surface fine-tuning.In contrast to traditional splicing and faring method based on a single control point of surface，this method improves the efficiency of surface's splicing and faring greatly，while the interior continuity of the surface and the precision of the offset points can be guaranteed. Key words：3D CAD models；complex bow and stern；ship-hull surfaces；splicing；fairing

U622.9

A

10.3969/j.issn.1673-3185.2016.04.002

2015-10-23網(wǎng)絡(luò)出版時(shí)間：2016-7-29 9:45

中國(guó)艦船研究設(shè)計(jì)中心研發(fā)基金資助項(xiàng)目

周春凱（通信作者），男，1973年生，碩士，高級(jí)工程師。研究方向：艦船總體設(shè)計(jì)與信息化技術(shù)。E-mail：zhouoldk@126.com

董楓，男，1988年生，博士生。研究方向：船舶CAD設(shè)計(jì)