改进Catmull-Clark细分算法及其在船用螺旋桨设计中的应用-王艳华苏洲

改进Catmull-Clark细分算法及其

在船用螺旋桨设计中的应用

2

王艳华1，

苏

洲2

2河海大学机电工程学院，江苏常州213022

摘

要：基于舰艇作战或巡航时隐身性能的需要，针对舰艇关键零件的设计和加工制造现状，结合应用日益广

1海军潜艇学院，山东青岛266042

泛的细分曲面造型方法，以Catmull-Clark细分算法为基础，提出一种基于网格边光顺程度计算的自适应细分算法。将该算法应用于船用螺旋桨的设计数值实验，由建模软件3DSMAX及其内部编程语言MAXScript来实现，并将实验结果与原始的Catmull-Clark算法细分结果、传统的顶点或边曲率法自适应细分结果进行了比较。数值实验结果表明，在满足零件设计现实需求的前提下，该自适应细分算法能较好地减少网格数量：当边的光顺度阈值取为0.825时，网格数降低了约28.83%，可提高计算速度，减小存储空间。该算法能克服以往基于均值计算的自适应算法存在的区分能力不足的影响。

关键词：舰艇隐身；关键零件；Catmull-Clark细分算法；自适应中图分类号：U664.33

文献标志码：A

文章编号：1673－3185（2012）05－50－05

AnImprovedCatmull-ClarkSegmentationAlgorithmandItsApplication

intheDesignofMarinePropellers

2

WANGYan-hua1，

SUZhou2

2CollegeofMechanicalandElectricalEngineering，HohaiUniversity，Changzhou213022，ChinaAbstract：Toimprovethestealthperformanceofnavalshipsduringcombatsorcruising，anadaptivesub divisionalgorithmbasedontheEdgeSmoothValueofmesheswasproposedinthispapertorefinepresenttechniquesofcriticalshippartsdesignandmanufacturing.ThisalgorithmisbasedontheCatmull-Clarksubdivisionalgorithm，andincorporatestheincreasinglypopularapplicationofsubdivisionsurfacemodel ing.Totestthealgorithm，severalnumericalexperimentsofmarinepropellerdesignwereconductedresort

1NavySubmarineAcademy，Qingdao266042，China

ingtothemodelingsoftware3DSMAXanditsinternalprogramminglanguageMAXScript.Theexperimen tionalvertexoredgecurvature-basedadaptivesubdivisionalgorithm.Ithasbeenfoundoutthatthenewadaptivesubdivisionalgorithmgreatlyreducesthetotalgridnumber：whenthethresholdofEdgeSmoothValueissetto0.825，adeductionof28.83%occursingridnumber，whichsignificantlyimprovesthecal foundinotheradaptivealgorithmsfordistinguishing.

talresultswerethencomparedwiththeoriginalCatmull-Clarksubdivisionalgorithmaswellastheconven

culationspeedandreducestherequireddatastoragespace.Also，thisalgorithmcansolvetheinadequacyKeywords：stealthofvessels；keyparts；Catmull-Clarksubdivisionalgorithm；self-adaption

收稿日期：2012－02－27

作者简介：王艳华（1976－），男，硕士研究生，讲师。研究方向：机械设计。E mail：wyh7622@

苏

洲（1967－），男，博士，副教授。研究方向：细分曲面造型、特征识别、机械设计。

通信作者：王艳华。

0引言

舰艇作为各国海军的重要国防力量，其优良的生存能力是衡量技战术性能的重要指标，其中就包括提高其自身隐身性能、减小被探测到的几率的能力。与军事强国相比，目前我国海军舰艇的隐身性能依然存在着一定的差距。提高舰艇隐身性能的方法很多，其中，在涉及舰艇隐身性能的关键零件的设计制造过程中，改进设计造型方法，主动提高舰艇零件的设计制造精度，降低可能产生的噪声辐射的等级或雷达反射波等就是一项行之有效的措施［1］。学术界认为，舰艇噪声辐射已成为影响舰艇隐蔽性的主要因素［2］，并且公认螺旋桨噪声是最主要的舰船噪声源［3］。但在舰船及其相关零部件的优化设计制造方面，我国目前还涉及较少。因此，在进行舰艇设计、制造时，对一些涉及舰艇隐身性能的关键零件，如螺旋桨等，有必要进行深入的声隐身性能研究，尽量使之符合降噪消音的原理，然后采用先进的造型方法，如自适应细分曲面技术来实现造型，用以在工作中就能达到降低噪声辐射的目的。在此，本文将只对其中的造型实现方法进行详细研究。

目前，国内外学者已对自适应细分造型方法进行了一定的研究，并也提出了一些行之有效的自适应细分算法和相关技术［4-8］。但这些算法在工程实践中仍然存在一定的局限性，特别是在舰艇零件设计方面的应用更少，这主要是因为现有的这些自适应方法在实现的过程中，其数值的计算受均值思想影响较重，对于控制网格光顺区域和非光顺区域的区分能力不足，从而使得部分已光顺区域网格过密。同时，对于模型中的一些尖锐特征和过渡区域，不能做出准确的判断，导致模型细分结果与极限曲面间的误差较大。本文将结合舰艇隐身性能的需要，以Catmull-Clark细分算法［9］为基础，提出一种能解决上述问题的自适应细分策略。

1Catmull-Clark细分算法

该细分模式由Catmull和Clark于1978年提

出，是一种针对以四边形网格为主的多边形控制网格的细分方法，细分曲面的极限是双三次B样条曲面，其规则包括两部分：几何规则和拓扑规则。

1.1几何规则

控制网格的新面点、C-C细分模式的几何规则主要是计算并生成

新边点和新顶点，具体细分面

如图1所示。

1/16

1/16γn/n

1/41/4βn

βn/n

1/4

1/4

3/83/8n/αβn/n

n（a）F-顶点

（1/16

b）E-顶点

1/16

（c）V-顶点

（1/2

d）边界E-顶点1/2

（1/8

e）边界3/4

V-顶点

1/8

Fig.1图1C-CC-Csubdivision细分模式几何规则

geometricrules

新面点，即F-顶点，如图1（a）所示。设原顶点分别为v1v2v3v4：

vF=4

(v1+v2+v3+v4)（1）

新边点，即E-顶点，如图1（b）所示。设原顶

点分别为v0v1v2v3v4v5：

vE=38(v0+v1)+116

(v2+v3+v4+v5)

（2）

新顶点，即V-顶点，如图1（c）所示：

βn2n-1γn2n-1

vV=αn+nåv2i+v2i+1（3）

i=0nåi=0

式中，βn=32n；γn=14n

；αn=1-βn-γn；n为顶

点的价。

边界E-顶点如图1（d）所示，设原顶点分别为v0v1：

vE=2

(v0+v1)（4）

边界V-顶点如图1（e）所示，设原顶点分别为v0v1v：

vV=(v0+v1)+384

v

（5）

1.2拓扑规则

面点和新顶点分别与新边点相连接，C-C细分模式的拓扑规则比较简单，由此生成一就是新

张新的控制网格。

1.3Catmull-Clark细分实验结果

利用三维建模软件3DSMAX及其内部编程

语言MAXScript，即可实现原始C-C细分算法（图2数增长。

），并且可以知道，模型控制网格的数量呈几何级2

基于Catmull-Clark模式的自适

应细分算法

仔细观察图2可发现，模型控制网格的数量

（a）初始控制网格（b）细分1次后的控制网格

（网格面片数15）

（网格面片数60）

（c）细分3次后的控制网格

（d）细分4次后的光照图（网格面片数960）

（网格面片数3860）

图Fig.22管路三通模型原始ModelsoriginalofC-CpipesubdivisionthreeC-Clinks算法细分结果

algorithm

basedon

呈几何级数急剧增长，若模型较复杂，就会给计算机的计算、显示和存储带来极大的困难，甚至会超出计算机的能力。通过对比模型细分3次和4次的结果可以看到，经过一定次数的细分后，网格数巨增，但模型的显示效果却没有明显变化，特别是在管口处尤为明显。鉴于此，为避免均值思想对网格光顺程度带来的不利影响，提出了一种针对模型中局部网格边的基于边光顺程度（EdgeSmooth因四边形网格被广泛应用于各类造型系统软

Value，ESV）的自适应策略。

件中，且其相对于其它形状的网格更适于加工，能够得到高精度的加工表面［10］，因而本文的算法将主要针对四边形网格。

2.1定义

给出了与规则相关的定义［11-12］，如图3所示。

Pk

iPk

2

fi+1

kVk

P1i

fi

P

k

n-1

Pk0

Pkn-2

图Fig.33基于边光顺度的自适应细分算法几何规则

GeometricalgorithmbasedrulesofonadaptiveESV

subdivision

定义1面边三角形（FaceEdgeTriangle）：以

某个四边形网格的面片中心为一个顶点，再以该面片任意一条边的两个端点为另两个顶点连接而

成的三角形，即

D(Vkk

ifiPi)

（6）

定义2面边三角形的法向量（Vector）：逆时

针环绕面边三角形三个顶点所得到的法向量，n即

0=(fiPi´fiPi+1)（7）

定义3边光顺度（EdgeSmoothValue，ESV）：若两个面边三角形共享同一条网格边，那么其法向量的夹角就称为共享边的光顺度，即

æçöESV=arccosn×n1ç0÷

èn0

´n1÷

（8）

ø定义4

光顺边（SmoothEdge）：给定一个误

差阈值εT，将计算得到的某一网格边的边光顺度与之比较，计算值小于或等于给定阈值的边就称为光顺边，计算值大于给定阈值的边就称为非光顺边。

定义5

不动顶点（NoMoveVertex）：对于与

同一条光顺边相邻的两个面片，可以指定其上的所有顶点均为不动顶点，这些顶点在下一级细分中的几何位置保持不变。

定义6

不动面片（NoMoveFace）：对于控制

网格中的某一面片，若其所有顶点都是不动顶点，则称该面片为不动面片。此定义也包括在网格中，与某一条光顺边相邻的面片即为不动面片。

2.2判据

针对局部网格，对于任意的网格边，给定一个阈值εT，然后对网格中各条边的光顺度进行计算，并将计算值与给定阈值进行比较。当某网格边的光顺度小于或等于给定阈值时，便将此网格边标定为光顺边，也称死边（DeadLine），其顶点可以指定为不动顶点，且在下一级细分中保持几何位置不变。否则，便将该网格边标定为非光顺边，也称活边（LiveLine），需要参加下一级的细分，直

到每个面片均为不动面片为止。

2.3算法的拓扑规则设计

自适应细分算法拓扑规则的建立过程如下：Step1

对模型控制网格中所有的边进行光

顺度计算，将计算结果与给定阈值进行比较并对相关几何元素的光顺性和可动性进行标记：当网格边的光顺度小于等于给定阈值时，将此边标定为光顺边，亦即死边；与此光顺边相邻的两个面片标定为光顺面，亦即死面（DeadFace）；若网格中环绕某一顶点的所有面均为光顺面，则标定该顶点为死点（DeadVertex）

。

Step2遍历控制网格的所有面片，按原始

C-C为不动面片细分模式生成新面点，（亦即死面），则不生成新面点。

若所遍历的某一面片Step3

遍历控制网格的所有边，按原始C-C

细分模式生成新边点，若所遍历的某一条边为光顺边，则用该边的中点代替刚生成的新边点。Step4遍历控制网格的所有顶点，按原始

C-C动顶点细分模式生成新顶点，（亦即死点），则直接用该顶点代替刚生成

若所遍历的顶点为不的新顶点。

Step5将新顶点和新面点分别与新边点相连接，生成新的自适应细分网格。

Step6

消除网格裂缝，生成最终模型。

该算法的实现流程如图4所示。

Fig.4图4Algorithm算法实现流程图flowchart

为了实现本文的算法，在大型三维建模软件3DS机交互界面如图MAX中，采用5所示。

MAXScript语言进行了编程，人Fig.5

Screenshots图5人机交互界面截图

ofhuman-computerinteraction

3数值实验结果

本文以船用螺旋桨为例（图6），通过建模软

件3DSMAX及其内部编程语言MAXScript实现了所提出的自适应细分算法，并将结果与原始C-C算法的细分结果进行了对比。

（a）初始控制网格（b）原始C-C方法细分结果

（c）本文自适应（阈值0.5）（d）本文自适应（阈值0.825）

Fig.6图6Comparison不同细分算法细分结果对比

differentsubdivisionoftheresultsalgorithms

basedon

螺旋桨模型的初始控制网格为1464个四边形面片，采用原始C-C算法进行细分，一次细分后的面片数就达到了5856个，并且由图6（b）可以看到，此种细分方式的每一次细分过程都是全局均匀细分，网格数增长非常快，并且已经光顺的区域仍参加下一次的细分，致使计算机的负担很重。

而采用本文提出的自适应细分算法对模型进行自适应细分，当边的光顺度阈值取为0.5时，便可以较准确地将已光顺区域从整个网格中区分出来，这样就可使细分过程只针对未光顺区域进行，其细分针对性更强，最终得到了4506个面片（图6当边的光顺度阈值取为（c）），相比原始算法，网格数降低了约0.825时，可得到更少的控23.05%；制网格数，为4168（图6（d）），相比原始算法，网格数降低了约28.83%，大大提高了计算机的运算、存储和显示速度，减小了存储空间。

对于采用本文算法设计的模型能否满足舰艇在降噪消音、推进功率、零件强度和刚度等方面的需求，已专门利用相关的分析软件ANSYS，FLU ENT值分析，等对模型进行了力学、验证了算法的正确性和可行性，流体力学等相应的数本文在此

不再赘述。4不同自适应算法对比分析

表2所示为采用传统顶点曲率自适应细分一

次和采用本文自适应细分一次的结果对比。

表2不同自适应细分算法细分数据对比

Tab.2Comparisonofdatabasedondifferent

adaptivesubdivisionalgorithms

顶点曲率自适应

本文自适应阈值初始控制网格数自适应网格数11与原始细分结果相比的网

34649573464962格数下降率/%

32.4332.34如图7所示，模型中圆圈所标示的范围即为采用不同的算法得到的细分结果。通过对比可以发现：采用传统的顶点曲率判别方法对模型进行自适应细分，由于顶点曲率的计算过程受均值思想的影响，在桨叶较光顺的区域网格密度仍较大，而在桨叶边缘等高曲率区域（或过渡区域），其网格密度与叶片表面较光顺区域的加密密度几乎没有差别，模型过渡特征不突出，被均匀化了；而与传统的曲率判别方法相比，在网格数量降低程度相当的情况下（数据对比见表2），采用本文的方法则能较准确地区分已光顺区域和未光顺区域（光顺区域的网格密度明显减小），能准确地捕捉到桨叶的过渡特征（例如，桨叶边缘等高曲率处），然后便可有针对性地适当加大该处网格的细分密度，满足造型的光顺性需求，最终使模型造型更合适。

（a）曲率自适应细分结果（b）本文自适应细分结果

Fig.7图Comparison7曲率判别法与本文方法细分结果对比图

oftheresultsadaptivecurvaturesubdivisionadaptivesubdivisionrespectivelyalgorithmproposedalgorithmsbasedonvertex

inandthis

thepaper

new结合表2和图7可得出结论：在网格数下降相当的情况下，采用本文的自适应算法对模型进行细分，细分过程更具针对性，模型高曲率区域的细分能得到加强，网格密度进一步加密，模型造型更加合理。

有关如何对该自适应细分算法进行误差控制

和精度分析，以及算法能否应用于其它细分模式和网格等问题，可作为下一步研究的重点。

5结论

在认真研究现有自适应细分技术的基础上，本文结合舰艇的隐身性能需求，提出了将自适应细分技术应用于舰艇关键零件设计和制造的工程实际，并以Catmull-Clark细分算法为基础，针对四边形网格提出了基于网格边光顺度计算的自适应细分方法，同时还以船用螺旋桨为例进行了数值实验。从实验结果可得到如下结论：

12）算法能较好地降低网格数量；算法具有较强的区分能力；

）对于网格中的光顺区域和非光顺区域，该的过渡区域，3）对于控制网格中的尖锐特征或曲率较高

该算法具有较强的捕捉能力，能进一

步加大这些区域的控制网格的密度，以使整个细分过程更具针对性，从而有效避免以顶点曲率等为判断阈值的均值思想所带来的对网格几何特征区分能力不足的不利影响。

参考文献：

1］龚锦伟，王振民，王立新，等.舰船隐身设计探讨［J］.

船舶，2003（3）：29-31，62.

GONGonnavalJW，WANGZM，WANGShipandLX，Boatetal.，2003Discussion29-31，62.

shipstealth［J］.（3）：2］缪旭弘，王振全.舰艇水下噪声控制技术现状及发展

对策［C］//第十届船舶水下噪声学术讨论会论文集，3］杨琼方，2005.

王永生，张明敏.舰艇螺旋桨水下噪声预测

［J］.船舶力学，2011，15（4）：435-442.

YANGQF，WANGYS，ZHANGMM.UnderwaterJournalnoisepredictionofShipMechanicsofshipand，2011submarine，15（4Adaptive）：435-442.propeller［J］.4］AMRESHA，FARING，RAZDANA.subdivi

sionandGeometricschemesforMethodstriangularinScientificmeshes［C］//Hierarchical5］linHU：Springer-VerlagHP，ZHANQ，，JIANG2003：319-327.

Visualization.Ber HB，etal.Adaptiverefine

mentsnalofHuazhonginCatmull-ClarkUniversitysubdivisionofScienceschemeandTechnology

［J］.Jour （NatureScienceEdition），2002，30（10）：56-58.6］

KobbeltL.

3-subdivision［C］

//ProceedingsoftheACMNewOrleansSIGGRAPH，2000：conference103-112.

oncomputergraphics.（下转第78页）

［［［［［［

2006，1（5/6）：83-89.

率汽轮发电机组系统仿真研究［J］.中国舰船研究，FANGWM，XUJ，HUANGGH.Simulationon

［7］孙中宁.核动力设备［M］.哈尔滨：哈尔滨工程大学

出版社，2004.

［8］GUIMARÃESLNF，OLIVEIRAND，BORGESEM，

etal.DerivationofaninevariablemodelofaU-tubesteamgeneratorcoupledwithathree-elementcontrol 1027-1043.

ler［J］.AppliedMathematicalModelling，2008，32（6）：［9］邱志强，邹海，孙建华.船用蒸汽发生器给水控制系

统半物理仿真技术研究［J］.舰船科学技术，2008，31（2）：141-145.

QIUZQ，ZOUH，SUNJH.Researchonsemi-physi calsimulationofmarinesteamgeneratorfeedwatercon 141-145.

low-powerturbo-generatorsetbasedonMatlab/Simu link［J］.ChineseJournalofShipResearch，2006，1（5/6）：83-89.

［3］SHIXP，WANGZC.Simulationstudyofanewmeth

odtosynchronouslycontrolrotatespeedandpowerofa15（6）：823-825，840.

steamturbine［J］.JournalofSystemSimulation，2003，［4］田兆斐，张志俭.船舶饱和蒸汽轮机动态特性仿真

［J］.汽轮机技术，2005，47（2）：90-92.

TIANZF，ZHANGZJ.SimulationonthedynamicTurbineTechnology，2005，47（2）：90-92.华中科技大学，1998.

XUJ.Simulationonthepowersystem［D］.Wuhan：［6］薛若军，田兆斐，赵强.核动力装置冷凝器实时仿真

［J］.核动力工程，2008，29（5）：119-123.

XUERJ，TIANZF，ZHAOQ.Real-timesimulationofnuclearpowercondensers［J］.NuclearPowerEngineer ing，2008，29（5）：119-123.

HuazhongUniversityofScienceandTechnology，1998.characteristicofmarinesaturatedstreamturbine［J］.［5］许建.常规潜艇动力系统综合仿真研究［D］.武汉：

trolsystem［J］.ShipScienceTechnology，2008，31（2）：［10］臧希年.核电厂系统及设备［M］.第2版.北京：清华

大学出版社，2010.

［11］ZHOUG，CHENX，YEWC，etal.Identificationof

dynamicsfornuclearsteamgeneratorwaterlevelpro cessusingRBFneuralnetworks［C］//The8thInterna struments，2007：379-383.

tionalConferenceonElectronicMeasurementandIn

［责任编辑：卢圣芳］

（上接第54页）

［7］CHENGKSD，WANGWP，QINH，etal.Fittingsub

［C］//Proceedingsofthe12thPacificConferenceonComputerGraphicsandApplication（PG’04）.Seoul，［8］WUJH，LIUWJ，WANGTR.Adaptiverefinement

schemeforsubdivisionsurfacesbasedontriangularConferenceonComputerAidedDesignandComputer［9］

Graphics（CAD-CG’05），2005：119-124.

CATMULLE，CLARKJ.Recursivelygeneratedmeshes［C］//ProceedingsoftheNinthInternational2004：16-24.

divisionsurfacestounorganizedpointdatausingSDM

［10］赵建民.基于Catmull-Clark模式的细分曲面NC刀

具轨迹生成技术［D］.大连：大连理工大学，2008.

［11］郑津津，朱巍，周洪军，等.非平均化自适应cat

mull-clark细分算法［J］.计算机应用与软件，2010，27（1）：6-8，22.

ZHENGJJ，ZHUW，ZHOUHJ，etal.Non-uniformComputer-AidedDesign，1978，10（6）：350-355.

adaptivecatmull-clarksubdivisionalgorithm［J］.Com 22.

puterApplicationsandSoftware，2010，27（1）：6-8，

［12］朱巍.细分曲面理论及其应用问题的研究［D］.合

肥：中国科学技术大学，2011.

B-splinesurfacesonarbitrarytopologicalmeshes［J］.

［责任编辑：卢圣芳］

本文来源：https://www.bwwdw.com/article/wrre.html