基于CST參數(shù)化的導(dǎo)彈三維幾何建模方法

關(guān)通 陳鵬云 沈劍 張國兵 馬英琪

引用格式：關(guān)通，陳鵬云，沈劍，等.基于CST參數(shù)化的導(dǎo)彈三維幾何建模方法［J］.航空兵器，2023，30（1）：59-63.

GuanTong，ChenPengyun，ShenJian，etal.3DGeometricModelingBasedonCSTParameterizationforMissiles［J］.AeroWeaponry，2023，30（1）：59-63.（inChinese）

摘要：針對建模效率低、成本高的問題，本文在二維CST參數(shù)化方法的基礎(chǔ)上，引入縱向輪廓描述函數(shù)，將原二維建模向三維建模擴展，建立能夠準確描述導(dǎo)彈三維形狀特征的解析函數(shù)。在此基礎(chǔ)上，通過在VB環(huán)境中對CATIA進行二次開發(fā)，使其能夠在CATIA中進行自動生成、修補和更新。幾何模型可以被導(dǎo)出和保存，并作為多學科設(shè)計與分析的基礎(chǔ)。結(jié)果表明，本文提出的三維CST建模技術(shù)可以滿足復(fù)雜幾何建模的要求，有效簡化了導(dǎo)彈參數(shù)化過程。

關(guān)鍵詞：CST參數(shù)化；CATIA二次開發(fā)；類別函數(shù)；形狀函數(shù)；參數(shù)化建模；導(dǎo)彈

中圖分類號：TJ760

文獻標識碼：A

文章編號：1673-5048（2023）01-0059-05

DOI：10.12132/ISSN.1673-5048.2022.0101

0引言

實時快速的幾何形狀建模是實現(xiàn)導(dǎo)彈多學科設(shè)計優(yōu)化的關(guān)鍵技術(shù)之一，其建模速度和建模質(zhì)量關(guān)系著多學科主模型的設(shè)計周期和設(shè)計質(zhì)量［1-4］。適用于導(dǎo)彈外形優(yōu)化的參數(shù)化方法應(yīng)該具備生成效率高、表達直觀等特點。傳統(tǒng)的幾何外形建模沒有建立幾何參數(shù)和屬性的關(guān)系，當幾何參數(shù)改變時，模型需要重新進行繪制，這種建模方法效率低、成本高，且對幾何模型中的某些曲線沒有精確的描述，造成幾何建模誤差大。因此，如何快速準確地生成導(dǎo)彈的幾何模型是必須解決的問題，幾何形狀參數(shù)化技術(shù)為此提供了一個有效的解決方法。

導(dǎo)彈的幾何形狀參數(shù)化設(shè)計是指用一組外形參數(shù)來比較精確地描述幾何外形，通過修改這些外形參數(shù)中的一個或多個，自動完成導(dǎo)彈外形的改動。傳統(tǒng)的幾何外形參數(shù)化建模主要有CAD方法、特征參數(shù)描述法、樣條方法以及二次擬合方法。直接采用CAD方法對于復(fù)雜幾何模型很難進行準確描述。周超［5］將導(dǎo)彈的所有特征參數(shù)提取出來，通過建立特征參數(shù)與部件的關(guān)系，從而完成導(dǎo)彈的參數(shù)化幾何建模。張以良等［6］采用貝塞爾曲線擬合方法，通過改變控制點的位置，實現(xiàn)幾何形狀的改變。但是隨著控制點數(shù)量的增多，該方法對曲線形狀的控制能力也變?nèi)酢Ｍ醭龋?］采用均勻有理B樣條曲線方法對徑向進行擬合，實現(xiàn)參數(shù)化建模。劉麗娜等［8］用有限個設(shè)計函數(shù)進行線性疊加，通過設(shè)置各個函數(shù)的權(quán)重，實現(xiàn)翼型的參數(shù)化建模。雖然這些方法能夠有效描述幾何形狀，但是涉及到的參數(shù)比較多，計算量大，不利于進行多學科設(shè)計與優(yōu)化。

針對上述問題，Kulfan等［9-11］提出了一種基于形狀函數(shù)/類別函數(shù)變換的二維截面參數(shù)化方法（ClassFunction/ShapeFunctionTransformation，CST）。由CST參數(shù)化方法中的類別函數(shù)生成幾何形狀的基本外形，再通過改變形狀函數(shù)中Bernstein多項式系數(shù)的權(quán)重進行模型修正，從而得到需要的幾何形狀。該方法涉及到的參數(shù)具有明確的外形幾何意義，能夠產(chǎn)生連續(xù)光滑的幾何模型，同時參數(shù)較少、可控性高，能夠快速生成各種復(fù)雜的二維截面圖形，魯棒性好。但CST參數(shù)化方法仍然以二維幾何模型為主，對三維幾何建模的應(yīng)用較少。

本文在二維CST參數(shù)化方法的基礎(chǔ)上，引入縱向輪廓描述函數(shù)，再將各個特征部件通過部件組合的方式組合起來，通過應(yīng)用VB編程語言和CATIA的二次開發(fā)技術(shù)，實現(xiàn)導(dǎo)彈外形的參數(shù)化設(shè)計，為導(dǎo)彈多學科設(shè)計與分析奠定基礎(chǔ)。

1二維CST參數(shù)化方法

二維CST參數(shù)化方法是通過形狀函數(shù)與類別函數(shù)相乘得到翼型的基本形狀［12-13］。類別函數(shù)組成基本形狀，

形狀函數(shù)對其進行修改，最終得到二維幾何模型。二維CST參數(shù)化的數(shù)學表達式為

ζU，L（ψ）=±ψN1·（1－ψ）N2·∑ni=0A·ψi+ψ·ζT（1）

式中：ψ=x/c，ψ∈［0，1］，ζ=z/c，ζT=zT/c，x為x軸坐標，z為z軸坐標，c為弦長，zT為后緣點的坐標，當ψ·ζT=0時，翼型為封閉的；ζ為相對厚度；ψ為歸一化的單位x軸長度；ψ為圓鼻翼型；（1－ψ）為一個尖的后緣翼型。

定義類別函數(shù)：

CN1N2（ψ）=ψN1·（1－ψ）N2（2）

式中：N1，N2為幾何形狀的類別參數(shù)。

通過改變幾何形狀的類別參數(shù)，可以得到不同的剖面形狀，如圖1所示。當N1=0.5，N2=1.0，基本形狀為NACA型的圓鼻剖面；當N1=0.5，N2=0.5，基本形狀為橢圓形；當N1=1.0，N2=1.0，基本形狀為雙凸尖頭型；當N1=1.0，N2=0.001，基本形狀為錐形。

定義形狀函數(shù)：

S（ψ）=ζ（ψ）－ψ·ζTψ·（1－ψ）（3）

該項為冪級數(shù)，可以用來表示任何滿足所需形狀的解析函數(shù)。

因此，前緣表示為

S（0）=2RLE/c（4）

后緣表示為

S（1）=tanβ+△ZTE/c（5）

對于閉合剖面，當S（ψ）=1，此時剖面形狀與類別函數(shù)生成的基本形狀一致。

引入Bernstein多項式，對于n階Bernstein多項式，有

S（ψ）=∑ni=0biKr，nψr（1－ψ）n－r（6）

式中：Kr，n=nr=n！r！（n－r）！；b=（b0，b1，b2，…，bn）為權(quán)重因子。

通過對Bernstein多項式中的各個項進行加權(quán)，可得各種各樣的幾何形狀。當bi=1時，即為單位形狀函數(shù)。

當已知幾何形狀的m個坐標點時，可以采用最小二乘法，將m個坐標點代入線性方程組：

ζ0（ψ0）ζ1（ψ0）…ζn（ψ0）ζ0（ψ1）ζ1（ψ1）…ζn（ψ1）ζ0（ψm）ζ1（ψm）…ζn（ψm）b0b1bn=ζ（ψ0）ζ（ψ1）ζ（ψm）（7）

式中：ζi（ψ）=CN2N1（ψ）Kin（ψ）i（1-ψ）n-i。

求解線性方程組（7），可得權(quán)重因子b=（b0，b1，b2，…，bn），進而可以得出符合建模要求的幾何形狀解析函數(shù)。

二維CST參數(shù)化方法只能夠描述模型二維剖面結(jié)構(gòu)，只有在剖面結(jié)構(gòu)的基礎(chǔ)上，通過CAD軟件建立展向特征參數(shù)與部件的關(guān)系，才能夠?qū)崿F(xiàn)三維幾何模型建模。但是該方法涉及到的參數(shù)較多，同時不能夠準確描述三維模型上下表面的特征。因此，在二維CST參數(shù)化方法的基礎(chǔ)上引入縱向輪廓描述函數(shù)，用于描述三維模型。

2三維CST描述方法

基于二維CST參數(shù)化得到翼型解析函數(shù)：

ζu，L（ψ）=±（ψN1·（1－ψ）N2·

∑ni=0biKr，nψr（1－ψ）n－r+ψ·ζT）（8）

對于三維幾何模型，可以將其看成一系列沿縱向的剖面的集合。通過在不同的縱向位置處使用不同的形狀函數(shù)，進而確定一個解析形狀函數(shù)面，最終得到整個外形的解析函數(shù)。因此在二維CST參數(shù)化方法的基礎(chǔ)上，引入縱向輪廓描述函數(shù)，得到三維幾何模型的解析函數(shù)。

將式（8）中的權(quán)重因子b沿縱向展開，令

bi=CM1M2（η）S（η）+η·ζM（9）

式中：η=x/L，η∈（0，1），L為縱向展長。將式（9）代入式（8），可得

ζu，L（ψ，η）=±（∑ni=0∑mj=0bi，jSyj（η）Sxi（ψ）·CN1N2（ψ）·CM1M2（η）+△ζM，N（ψ，η））（10）

式中：bi，j為三維幾何模型表面的權(quán)重因子;b=（m×n），m，n分別為相對y軸和x軸展開的Bernstein多項式階數(shù)；Syj（η）=Kj，mηj（1－η）m－j；Sxi（ψ）=Ki，nψi（1－ψ）n－i。

因此，三維幾何模型的上下表面形狀由N1，N2，M1，M2和b來控制，表面控制點數(shù)為（n+1）×（m+1）個。由于式（10）中的ψ和η都歸一化，ψ∈［0，1］，η∈［0，1］，不能夠描述y坐標與η和ψ之間的關(guān)系，故對y軸引入控制方程。

對于具有對稱結(jié)構(gòu)的三維幾何模型，y軸的控制方程為

y（ψ，η）=CT1T2（η）∑wt=0bt·ST（η）·（1－2·ψ）（11）

將式（11）轉(zhuǎn)換到全局坐標系下，三維CST參數(shù)化外形曲面可表示為

X（η）=η·L

Y（ψ，η）=CT1T2（η）∑wt=0bt·ST（η）·（1－2·ψ）·（W/2）

ZU，L（ψ，η）=±（CN1N2（ψ）·CM1M2（η）∑ni=0∑mj=0bi，jSyj（η）·

Sxi（ψ）+△ζM，N（ψ，η））·（H/2）（12）

式中：η∈［0，1］;ψ∈［0，1］。

假設(shè)L=W=H=1，令形狀函數(shù)S=0.52·N1，類別函數(shù)參數(shù)N1=N2=M1=M2=T1=T2=0.001，則可得到一個單位立方體，如圖2所示。

通過改變類別函數(shù)的參數(shù)，可得到不同的關(guān)于y=0對稱的三維幾何模型，如圖3所示。

對于機翼等非對稱的三維模型，y軸的控制方程為

y=c（η）·ψ+η·tanα（13）

式中：α為機翼前緣與x軸的夾角；c（η）為η處的弦長。通過改變α的大小可以控制機翼前緣的傾斜程度，通過改變c（η）的大小可以控制后緣的傾斜程度。

將式（13）轉(zhuǎn)換到全局坐標系下，可得

X（η）=η·LY（ψ，η）=c（η）·ψ+η·tanα·WZU，L（ψ，η）=±（CN1N2（ψ）·CM1M2（η）∑ni=0∑mj=0bi，jSyj（η）·

Sxi（ψ）+△ζM，N（ψ，η））·（H/2）（14）

機翼的三維幾何模型，如圖4所示。

3CATIA的二次開發(fā)

CATIA本質(zhì)上是一個OLE自動化對象服務(wù)器，CATIA對象可以有多種方式進行訪問［14-16］，目前主要訪問方式有兩種：一種是在進程內(nèi)進行訪問，另一種為進程外的。進程內(nèi)訪問即腳本和CATIA在同一進程內(nèi)運行，由CATIA的腳本引擎來解析執(zhí)行宏腳本命令，在進程內(nèi)一般通過VBScrip語言或VBA語言。進程外訪問是通過COM接口來訪問CATIA內(nèi)部對象，可以在VB程序中嵌入訪問。本文選用VB程序采用進程外的方式對CATIA進行二次開發(fā)。

進程外訪問CATIA對象，實現(xiàn)CATIA的自動啟動是二次開發(fā)中的最基本部分。

CATIA二次開發(fā)的流程如圖5所示。首先通過計算得出三維幾何模型的坐標。通過將這些坐標點讀取到CATIA中，再借用CATIA強大的曲面功能，可以生成三維幾何外形，在此基礎(chǔ)上進行其他學科的分析。當參數(shù)發(fā)生改變時，幾何圖形的坐標點發(fā)生變化，VB程序自動進行模型的更新與導(dǎo)出，實現(xiàn)整個流程的自動化。

4算例驗證

為了驗證本文方法的實用性，將某導(dǎo)彈結(jié)構(gòu)分為頭身組合體、舵面兩部分，采用三維CST參數(shù)化建模方法建立了導(dǎo)彈的各個部件，其主要幾何外形控制參數(shù)如表1所示。

在此基礎(chǔ)上，采用部件組合的方法將各部分相互組合，得到導(dǎo)彈的三維幾何模型，如圖6所示。

但是在組合時，舵面和彈體之間會產(chǎn)生一定的縫隙，因此不能直接用于后續(xù)學科的分析。借助CATIA強大的曲面構(gòu)造能力，將導(dǎo)彈的坐標點讀入到CATIA中，利用VB語言，可以在CATIA中自動生成導(dǎo)彈模型，并對產(chǎn)生的縫隙進行自動修補，如圖7所示。該模型可以保存并導(dǎo)出，直接進行其他學科的分析。

在此基礎(chǔ)上，對導(dǎo)彈的精度進行擬合，擬合精度用上下表面縱坐標的標準差進行衡量：

σ=∑ni=1（yi－y0i）2/（n－1）（15）

式中：σ為標準差；yi為CST參數(shù)化擬合后的縱坐標；n為坐標點的個數(shù)。

對導(dǎo)彈縱截面進行擬合，擬合結(jié)果如圖8所示。擬合后的標準差為6.1036×10-6。結(jié)果表明，三維CST參數(shù)建模擬合精度能夠滿足三維幾何建模要求。

5結(jié)論

本文在二維CST參數(shù)化的基礎(chǔ)上，引入縱向輪廓描述函數(shù)，得到用于描述導(dǎo)彈各個特征部件的三維解析函數(shù)，通過部件組合的方法將各部件進行組合，最后通過對CATIA進行二次開發(fā)，實現(xiàn)了對導(dǎo)彈的參數(shù)化建模。該幾何模型可以導(dǎo)出和保存，并作為其他學科分析的基礎(chǔ)。本文研究可以得到以下結(jié)論：

（1）本文提出的三維CST參數(shù)化建模能夠方便地得到幾何模型的解析函數(shù)，所需要的參數(shù)少，可以滿足不同形狀的建模要求。

（2）通過將各種部件進行組合，可以得到一個復(fù)雜的三維幾何模型。該方法可以快速、穩(wěn)定地構(gòu)建各類導(dǎo)彈的外形。

（3）三維CST參數(shù)化建模方法與CATIA二次開發(fā)相結(jié)合，能夠自動進行模型的修補與更新，實現(xiàn)整個流程的自動化，并進行后續(xù)的多學科分析與優(yōu)化。

參考文獻：

［1］徐智強，章麗麗，吳會菊，等.航空彈藥彈體結(jié)構(gòu)強度分析及優(yōu)化設(shè)計［J］.兵器材料科學與工程，2021，44（5）：96-102.

XuZhiqiang，ZhangLili，WuHuiju，etal.ProjectileStructureStrengthAnalysisandOptimizationDesignofAviationAmmunition［J］.OrdnanceMaterialScienceandEngineering，2021，44（5）：96-102.（inChinese）

［2］陳黎，段鵬飛，袁成.隱身空空導(dǎo)彈發(fā)展現(xiàn)狀及關(guān)鍵技術(shù)研究［J］.航空兵器，2022，29（1）：14-21.

ChenLi，DuanPengfei，YuanCheng.ResearchonDevelopmentStatusandKeyTechnologiesofStealthAir-to-AirMissiles［J］.AeroWeaponry，2022，29（1）：14-21.（inChinese）

［3］SumnuA，GuzelbeyIH.CFDSimulationsandExternalShapeOptimizationofMissilewithWingandTailfinConfigurationtoImproveAerodynamicPerformance［J］.JournalofAppliedFluidMecha-nics，2021，14（6）：1795-1807.

［4］ChenJF，F(xiàn)anXQ，XiongB，etal.ShapeOptimizationoftheCross-SectionforNoncircularHypersonicMissileForebody［J］.InternationalJournalofAerospaceEngineering，2020（9）：1-9.

［5］周超.基于CATIA二次開發(fā)的導(dǎo)彈參數(shù)化建模［J］.科學與財富，2015（13）：39-41.

ZhouChao.MissileParametricModelingBasedonCATIASecondaryDevelopment［J］.Sciences&Wealth，2015（13）：39-41.（inChinese）

［6］張以良，熊鷹，時立攀.一種螺旋槳參數(shù)化建模方法［J］.艦船科學技術(shù)，2015，37（1）：34-38.

ZhangYiliang，XiongYing，ShiLipan.AMethodofPropellerPara-metricModelling［J］.ShipScienceandTechnology，2015，37（1）：34-38.（inChinese）

［7］王超，林揚，胡志強，等.基于均勻有理B樣條的小水線面雙體無人艇參數(shù)化建模方法［J］.艦船科學技術(shù)，2017，39（23）：143-148.

WangChao，LinYang，HuZhiqiang，etal.AParametricModellingMethodofUnmannedSurfaceVehiclewithSmallWaterplaneAreaTwin-HullBasedonUniformRationalBSpline［J］.ShipScienceandTechnology，2017，39（23）：143-148.（inChinese）

［8］劉麗娜，吳國新.基于Hicks-Henne型函數(shù)的翼型參數(shù)化設(shè)計以及收斂特性研究［J］.科學技術(shù)與工程，2014，14（30）：151-155.

LiuLina，WuGuoxin.ResearchonApplicationofHicks-HenneFunctioninAirfoilShapeParameterization&Convergence［J］.ScienceTechnologyandEngineering，2014，14（30）：151-155.（inChinese）

［9］KulfanBM.RecentExtensionsandApplicationsofthe“CST”UniversalParametricGeometryRepresentationMethod［C］∥7thAIAAAviationTechnology，IntegrationandOperationsConference（ATIO），2007.

［10］廖炎平，劉莉，龍騰.幾種翼型參數(shù)化方法研究［J］.彈箭與制導(dǎo)學報，2011，31（3）：160-164.

LiaoYanping，LiuLi，LongTeng.TheResearchonSomePara-meterizedMethodsforAirfoil［J］.JournalofProjectiles，Rockets，MissilesandGuidance，2011，31（3）：160-164.（inChinese）

［11］KulfanBM.UniversalParametricGeometryRepresentationMe-thod［J］.JournalofAircraft，2008，45（1）：142-158.

［12］WangXG，YangJ，GuoJ，etal.AGeneralIntegratedMethodforDesignAnalysisandOptimizationofMissileStructure［J］.Algorithms，2019，12（12）：257.

［13］HeW，LiuX.ImprovedAerofoilParameterisationBasedonClass/ShapeFunctionTransformation［J］.TheAeronauticalJournal，2019，123（1261）：310-339.

［14］YangHG，WangSJ，XieYR.ParametricDesignofCrankshaftBasedonCATIASecondaryDevelopmentPlatform［J］.AppliedMechanicsandMaterials，2014，602/603/604/605：287-290.

［15］GuanP，DongJC，WuLG.SecondaryDevelopmentofBladesParametricDesignofHydraulicTorqueConverterBasedonCATIA［J］.AdvancedMaterialsResearch，2011，383/384/385/386/387/388/389/390：1055-1061.

［16］胡挺，吳立軍.CATIA二次開發(fā)技術(shù)基礎(chǔ)［M］.北京：電子工業(yè)出版社，2006.

HuTing，WuLijun.CATIASecondaryDevelopmentTechnologyFoundation［M］.Beijing：PublishingHouseofElectronicsIndustry，2006.（inChinese）

3DGeometricModelingBasedonCSTParameterizationforMissiles

GuanTong，ChenPengyun*，ShenJian，ZhangGuobing，MaYingqi

（NorthUniversityofChina，Taiyuan030051，China）

Abstract：Aimingattheproblemoflowefficiencyandhighcostofmissilemodeling，thispaperintroduceslongitudinalcontourdescriptionfunctiononthebasisof2DCSTparameterizationtechnology，extendstheoriginal2Dmodelingto3Dmodeling，andestablishesananalyticalfunctionthatcanaccuratelydescribethe3Dshapecharacteristicsofthemissile.Onthisbasis，thegeometricmodelcanbeautomaticallygenerated，repairedandupdatedinCATIAthroughsecondarydevelopmentinVBenvironment.Geometricmodelscanbederivedandpreservedasabasisformultidisciplinarydesignandanalysis.Theresultsshowthatthe3DCSTmodelingtechnologyproposedinthispapercanmeettherequirementsofcomplexgeometricmodeling，effectivelysimplifythemissileparameterizationprocess.

Keywords：CSTparameterization；CATIAsecondarydevelopment；classfunction；shapefunction；parametricmodeling;missile

收稿日期：2022-05-19

基金項目：國家自然科學基金項目（51909245;62003314）；航空科學基金項目（2019020U0002）；

山西省基礎(chǔ)研究計劃（201901D211244;202103021224187）

作者簡介：關(guān)通（1997-），男，山西臨汾人，碩士研究生。

*通信作者：陳鵬云（1987-），男，山西太原人，博士，副教授。