求解具有規(guī)定測地曲率曲線的高效幾何方法

緱建杰，張 鵬，楊 揚(yáng)，陳代鑫

（1.成都飛機(jī)工業(yè)（集團(tuán)）有限責(zé)任公司檢驗(yàn)檢測部，四川 成都 610091 2.長沙理工大學(xué)汽車與機(jī)械工程學(xué)院，湖南 長沙 410114）

1 引言

曲面上的曲線在曲面裁剪[2]、曲面過渡[2]、數(shù)控加工刀具路徑規(guī)劃[3]等眾多應(yīng)用中發(fā)揮著重要的作用。在一些應(yīng)用場景中，需要設(shè)計出的曲線上測地曲率滿足一定的要求，這種曲線被稱為具有規(guī)定測地曲率的曲線。

例如，在自動鋪帶工藝中，鋪放路徑上測地曲率的大小直接決定了復(fù)合材料預(yù)浸帶的面內(nèi)變形[4]。由于預(yù)浸帶面內(nèi)變形能力很小，因此帶料通常沿著模具曲面上的測地線鋪放[5]。

測地線是具有規(guī)定測地曲率曲線中的一個特例，其上各處測地曲率均為零。自動鋪帶軌跡規(guī)劃中經(jīng)常要處理測地線求解的初值問題，即給出鋪放起始點(diǎn)位置以及初始鋪放方向，求解曲面上的唯一一條測地線。

到目前為止，該問題已經(jīng)有很多經(jīng)典的求解方法，具體可分為數(shù)值解法[6-7]、離散解法[8-9]以及幾何解法[10-11]。

相較于測地線，本研究主要關(guān)注參數(shù)曲面上測地曲率不為零的曲線，這類曲線在圖案設(shè)計、齒輪設(shè)計、機(jī)械結(jié)構(gòu)設(shè)計等場合均有所應(yīng)用。

文獻(xiàn)[12]提出了一種保測地曲率的曲面曲線設(shè)計方法，其目的是在光滑曲面上進(jìn)行圖案設(shè)計。

在該項(xiàng)研究中，所需的曲面圖案首先被構(gòu)造為一平面圖案。在曲面上指定起始點(diǎn)和方向后，該平面圖案被逐段映射到目標(biāo)曲面上。

映射過程中，為保證視覺效果，要求曲面曲線的測地曲率與平面曲線的曲率相等。

文獻(xiàn)[13]基于保測地曲率的曲線設(shè)計方法對非圓錐齒輪傳動進(jìn)行了設(shè)計與分析，利用球面上具有規(guī)定測地曲率的曲線生成了非圓錐齒輪的空間三維齒廓。

文獻(xiàn)[14]則基于球面上具有規(guī)定測地曲率的曲線設(shè)計了一種不等速行星輪系分插機(jī)構(gòu)。

具有規(guī)定測地曲率的曲線也常應(yīng)用于復(fù)合材料自動化加工領(lǐng)域中。

文獻(xiàn)[15]提出了一種針對圓錐面模具的自動鋪絲軌跡規(guī)劃方法，利用等測地曲率曲線實(shí)現(xiàn)了纖維變形與纖維方向之間的平衡。

該方法所依據(jù)的幾何學(xué)原理是可展曲面上的等測地曲率曲線實(shí)際上是其對應(yīng)展開面上的一段圓弧。

文獻(xiàn)[16]給出了一種針對自由曲面模具的自動鋪帶軌跡規(guī)劃方法，利用具有微小測地曲率的曲線實(shí)現(xiàn)了預(yù)浸帶在模具曲面上的轉(zhuǎn)向，消除了相鄰帶料間過大的鋪放間隙。

本研究提供一種參數(shù)曲面上具有規(guī)定測地曲率曲線的求解方法，當(dāng)給出目標(biāo)曲面上的起始點(diǎn)、方向以及測地曲率的分布時，可追蹤出符合設(shè)計要求的空間曲線。

其主要由以下幾個部分組成：第二節(jié)對現(xiàn)有的計算方法進(jìn)行了分析；第三節(jié)詳細(xì)介紹了提出方法的原理；第四節(jié)將所提出的方法與現(xiàn)有的方法在計算效率上進(jìn)行了比較；第五節(jié)展示了本方法的兩種應(yīng)用場景；第六節(jié)則進(jìn)行總結(jié)并給出結(jié)論。

2 研究現(xiàn)狀與分析

在設(shè)計制造的一些應(yīng)用場景中，需要設(shè)計出的曲面曲線上測地曲率滿足一定的要求，并且測地曲率大小不為零。然而，現(xiàn)有的研究主要處理測地線的計算問題，其中多數(shù)方法無法直接用于構(gòu)造具有非零測地曲率的曲線，例如數(shù)值解法以及離散解法這兩類。

到目前為止，具有非零測地曲率曲線的計算方法主要可分為數(shù)值解法、分析解法以及幾何解法這三類。

文獻(xiàn)[16]基于微分幾何建立了具有規(guī)定測地曲率曲線的數(shù)學(xué)模型，該數(shù)學(xué)模型由一個微分方程組表達(dá)，需要利用數(shù)值積分器進(jìn)行求解。

其中曲線的設(shè)計精度與數(shù)值積分器的階數(shù)直接相關(guān)。文獻(xiàn)[15]提出的分析解法可以構(gòu)造出準(zhǔn)確的目標(biāo)曲線，然而該方法只適用于可展曲面，應(yīng)用范圍較窄。

文獻(xiàn)[12]提出的幾何解法在每一步計算中需要將點(diǎn)投影回曲面，降低了計算效率。

文獻(xiàn)[11]提出了一種構(gòu)造測地線的幾何方法，經(jīng)過調(diào)整，該方法可用于求解具有非零測地曲率的曲線。該方法則直接在參數(shù)空間中進(jìn)行計算，無需額外的投影操作。

現(xiàn)有的幾何方法主要基于伏雷內(nèi)公式與泰勒展開設(shè)計目標(biāo)曲線。設(shè)Q(s)—目標(biāo)曲線上一點(diǎn)，Q'(s)—曲線在該點(diǎn)的延伸方向，Q''(s)—曲線在該點(diǎn)的曲率。根據(jù)現(xiàn)有的幾何方法可知，Q(s)在目標(biāo)曲線上臨近點(diǎn)的位置以及方向可表達(dá)為：

式中：h—計算步長。

在得到臨近點(diǎn)位置及延伸方向后，可采用公式追蹤其它離散點(diǎn)從而得到整段曲線。

需要注意的是，文獻(xiàn)[12]中提出幾何方法是在三維空間中追蹤曲線，而文獻(xiàn)[11]中的幾何方法則在參數(shù)空間中追蹤曲線。

由方程可知，現(xiàn)有的幾何方法可以保證二階位置計算精度與一階方向計算精度。

設(shè)延伸方向Q'(s)計算誤差為Ah2，其中A—常數(shù)。若忽略曲率Q''(s)的計算誤差，則位置計算誤差主要來源于以下兩個部分：

（1）由于方向計算誤差引起的誤差：

（2）泰勒展開的余項(xiàng)：

顯然，在每一步計算中，由于方向計算誤差引起的誤差ε1相對于泰勒展開的余項(xiàng)ε2不能忽略。

因此，現(xiàn)有的幾何方法雖然具有二階位置計算精度，但計算過程中由于方向計算誤差引起的誤差ε1不斷累積，降低了最終精度。

由前述分析可知，若想獲得較好的曲線設(shè)計精度，需保證方向計算精度為高階，而這正是現(xiàn)有的幾何方法難以保證的，其主要原因在于這些方法難以得到比曲率Q''(s)更高階的幾何信息。

對現(xiàn)有的規(guī)定測地曲率曲線計算方法的主要特點(diǎn)進(jìn)行了總結(jié)與展示，如表1所示。

表1 現(xiàn)有計算方法的主要特點(diǎn)Tab.1 Main Features of Existing Methods

3 求解方法

本小節(jié)介紹文中所提出方法的數(shù)學(xué)原理。設(shè)r(u，v)=[x(u，v)，y(u，v)，z(u，v)]為模具曲面S的參數(shù)方程，其中x、y、z是參數(shù)u和v的可微函數(shù)。

令C是曲面S上一條具有規(guī)定測地曲率的弧長參數(shù)化曲線，如圖1所示。

圖1 測地曲率與法曲率的定義Fig.1 Definition of Geodesic Curvature and Normal Curvature

圖中：P—曲線C上任意一點(diǎn)；N—曲面S在點(diǎn)P的單位法向量；t—曲線C在P點(diǎn)的單位切向量。它可以用如下方程描述：

式中：s—弧長參數(shù)。

令b=N×t，則N，t，b三者彼此正交并構(gòu)成一個右手坐標(biāo)系。根據(jù)微分幾何，曲線C在P點(diǎn)處的曲率向量k可表示為：

式中：kn—法曲率向量；

kg—測地曲率向量；

kn—曲面S在P點(diǎn)處沿切向量t的法曲率；

kg—曲線C在P點(diǎn)處的測地曲率。

由公式可知，曲線C在P點(diǎn)處的測地曲率kg可通過如下公式求解：

設(shè)P0為曲線C的初始點(diǎn)；

C'與P0'為曲線C與點(diǎn)P0在參數(shù)空間中的逆映射。

根據(jù)微分幾何，曲線C'在P0'處的延伸方向l0=與曲線C在P0處的單位切向量t0之間的關(guān)系可通過下式表達(dá)：

在方程左右兩側(cè)同時點(diǎn)乘列向量[ru(P0')rv(P0')]T得到：

3.1 構(gòu)建求解位置與方向信息的F算子

令P0為曲面上的正則點(diǎn)，則由公式可得，已知曲線C在三維空間中的單位切向量t0，可以求得C'在參數(shù)空間中的延伸方向l0為：

式中：E(P0')，F(xiàn)(P0')，G(P0')—曲面在P0'處的第一類基本量。

由點(diǎn)P0'位置以及延伸方向l0，可以通過一階線性逼近得到P0'臨近點(diǎn)P1'的位置：

在獲得P1'的位置之后，為保證追蹤過程能夠持續(xù)進(jìn)行，需計算曲線C在P1處的單位切向量t1以及C'在P1'處的延伸方向l1。由公式可知，單位切向量t0在P0點(diǎn)相對弧長的變化率可表達(dá)為：

在得到單位切向量t1后，可以利用與類似的公式計算l1。

顯然，通過以上步驟獲得的P1'與t1只具有一階精度，在多數(shù)場合中直接應(yīng)用會有較大的計算誤差。為提高計算精度，首先將前述具有一階精度的計算過程整體封裝為一個F算子：

如上式所示，F(xiàn)算子輸入為C'上某一點(diǎn)P'的位置以及C在P處的單位切向量t，輸出為C'在P'的延伸方向l以及單位切向量t在P點(diǎn)相對弧長的變化率

具體對于P0'來說，F(xiàn)算子可表示為：

3.2 通過多點(diǎn)平均提高位置、方向計算精度

為提高P1'與t1的計算精度，可采用與龍格-庫塔法相似的原理，在P0'與P1'之間設(shè)置多個中間點(diǎn)，然后對中間點(diǎn)的幾何信息進(jìn)行加權(quán)平均以精確預(yù)測曲線C'在P0'附近的幾何形狀。

若需保證P1'與t1具有二階計算精度，則依據(jù)二階龍格-庫塔法在P0'與P1'之間設(shè)置一個中間點(diǎn)

基于公式計算具有二階精度的P1'與t1，得到：

式中：2P1'與2t1的左上標(biāo)2表示這些變量具有二階計算精度。

若需保證P1'與t1具有四階計算精度，則依據(jù)經(jīng)典四階龍格-庫塔法在P0'與P1'之間設(shè)置三個中間點(diǎn)它們采用如下公式進(jìn)行計算：

式中：4P1'與4t1的左上標(biāo)4—這些變量具有四階計算精度。

需要注意的是，迭代過程中方程與方程無法保證2t1與4t1為單位向量，這可通過下式進(jìn)行調(diào)整：

在求得P1'的位置以及曲線C在P1處的單位切向量t1后，可以采用與前述相似的步驟求得后續(xù)點(diǎn)P2'、P3'…..，直到曲線抵達(dá)曲面邊界或離散點(diǎn)的數(shù)量滿足規(guī)定要求。

顯然，這些離散點(diǎn)[P0'、P1'、P2'......]定義了參數(shù)空間中具有規(guī)定測地曲率的曲線C'。

4 驗(yàn)證

為驗(yàn)證所提出方法的有效性，本節(jié)將該方法與現(xiàn)有的方法在計算精度以及計算時間方面進(jìn)行比較。

在接下來部分給出了三個算例。每個算例中采用現(xiàn)有的方法計算參數(shù)曲面上一系列具有規(guī)定測地曲率的曲線，直到曲線抵達(dá)曲面邊界。同時，記錄計算過程所消耗的時間。

另一方面，以計算曲線末端點(diǎn)與理想曲線末端點(diǎn)在三維空間中的距離Error作為評價計算精度的指標(biāo)。

顯然，Error越小表示方法越精確。由于很難通過解析方法獲得一般曲面上理想的具有規(guī)定測地曲率的曲線，因此本節(jié)中的理想曲線均采用數(shù)值解法求得。

為保證理想曲線具有足夠高的精度，數(shù)值解法采用了RK4數(shù)值積分器以及極小的計算步長。

為便于比較計算精度，計算過程中令各方法的計算步長保持恒定。

文中所有的代碼均在MATLAB 2016b環(huán)境中實(shí)現(xiàn)，然后運(yùn)行于個人電腦[Intel Core?i5-8400，16G]上來測試計算效率。

如圖2所示，算例1計算一個2×2B樣條曲面上的等測地曲率曲線。

圖2 2×2B樣條曲面上的等測地曲率曲線Fig.2 Curves with Constant Geodesic Curvature on a 2×2B Spline Surface

計算結(jié)果，如圖3、圖4以及表2所示。

圖4 算例1中計算時間對計算精度的影響Fig.4 Influence of Time Costs on Accuracy in Example 1

表2 算例1中各方法所消耗的計算時間Tab.2 Time Costs of Each Method in Example 1

圖3 算例1中計算步長對計算精度的影響Fig.3 Influence of Step Size on Accuracy in Example 1

其中，不同方法在相同計算步長下的計算精度以及計算步長對計算精度的影響，如表3所示。

不同方法在相同時間內(nèi)所能獲得的精度，如表4所示。計算過程消耗的時間，如表2所示。

在以下圖示中，RK4、RK2表示采用文獻(xiàn)[16]中的數(shù)值方法以及RK4、RK2數(shù)值積分器進(jìn)行求解，Hu’s表示采用文獻(xiàn)[12]中的幾何方法進(jìn)行求解，Zhang’s表示采用文獻(xiàn)[11]中的幾何方法進(jìn)行求解，Our4、Our2則表示采用提出的四階與二階方法進(jìn)行求解。

另外，為減小文章篇幅，表2、表3、表4 只列出了部分實(shí)驗(yàn)數(shù)據(jù)。

表4 算例3中各方法所消耗的計算時間Tab.4 Time Costs of Each Method in Example 3

算例2計算一個雙三次B樣條曲面上的等測地曲率曲線，如圖5所示。

圖5 雙三次B樣條曲面上的等測地曲率曲線Fig.5 Curves with Constant Geodesic Curvature on a 3×3B Spline Surface

各方法在相同計算步長下的精度以及計算步長對精度的影響，如圖6所示。

圖6 算例2中計算步長對計算精度的影響Fig.6 Influence of Step Size on Accuracy in Example 2

它們在相同時間內(nèi)所能獲得的精度，如圖7所示。

圖7 算例2中計算時間對計算精度的影響Fig.7 Influence of Time Costs on Accuracy in Example 2

計算過程消耗的時間，如表3所示。

表3 算例2中各方法所消耗的計算時間Tab.3 Time Costs of Each Method in Example 2

算例3計算一個6×7B樣條曲面上的等測地曲率曲線，如圖8所示。

圖8 6×7B樣條曲面上的等測地曲率曲線Fig.8 Curves with Constant Geodesic Curvature on a 6×7b Spline Surface

計算結(jié)果，如圖9、圖10以及表4所示。

圖9 算例3中計算步長對計算精度的影響Fig.9 Influence of Step Size on Accuracy in Example 3

其中，計算精度隨步長以及計算時間的變化規(guī)律，如圖9、圖10所示。計算過程消耗的時間，如表4所示。

圖10 算例3中計算時間對計算精度的影響Fig.10 Influence of Time Costs on Accuracy in Example 3

由圖2、圖5、圖8可知，測地曲率的幅值決定了曲線偏離測地線的程度，測地曲率的正負(fù)決定了曲線所處的方位。

具體地說，當(dāng)測地曲率為正時，曲線位于測地線延伸方向的左側(cè)；當(dāng)測地曲率為負(fù)時，曲線位于測地線延伸方向的右側(cè)。

首先，考慮計算精度問題。

由圖3、圖6、圖9可看出，隨著計算步長的減小，各方法的計算精度均有所提高。

同時，方法階數(shù)越高，收斂速度越快。對于四階方法（RK4與Our4），當(dāng)計算步長很小時，計算結(jié)果容易產(chǎn)生震蕩，此時繼續(xù)減小步長對改善計算精度的效果不顯著。由計算結(jié)果可看出，現(xiàn)有的幾何方法精度較差，這主要是由于這些方法無法保證較高的方向計算精度。

對于二階方法而言，文中提出的方法在精度上要優(yōu)于文獻(xiàn)[16]中的數(shù)值解法；對于四階方法而言，文中提出的方法在多數(shù)情況下可以獲得更好的計算精度。

另外，計算結(jié)果表明四階方法的計算精度要顯著高于二階方法，這主要是由于所構(gòu)造的曲線光滑連續(xù)，在解光滑的情況下，高階方法通常優(yōu)勢明顯。其次，考慮計算時間。由表2、表3、表4可看出，隨著計算步長的減小，計算過程消耗的時間顯著增加。

由上表可看出，與四階方法相比，現(xiàn)有的幾何方法消耗的計算時間更少，這主要是由于這些方法在計算過程中無需進(jìn)行多點(diǎn)平均、多次迭代。

文獻(xiàn)[12]中的幾何方法計算速度較慢，這主要是由于該方法在每一步計算中需要進(jìn)行點(diǎn)到曲面的投射操作。通過比較發(fā)現(xiàn)，本方法在時間消耗上與現(xiàn)有的數(shù)值方法大致相當(dāng)。

為綜合考慮計算精度與計算時間，圖4、圖7、圖10 展示了不同方法在相同時間內(nèi)所能獲得的計算精度。結(jié)果表明，在相同時間內(nèi)，本方法獲得的計算精度一般要優(yōu)于現(xiàn)有的數(shù)值解法，即計算效率更佳。同時，現(xiàn)有的數(shù)值方法在計算效率上要優(yōu)于幾何方法。

5 應(yīng)用

5.1 自動鋪帶

自動鋪帶技術(shù)是一種增量制造技術(shù)，適用于制造翼面、壁板等大尺寸、小曲率復(fù)合材料構(gòu)件，具有加工成本低、效率高等特點(diǎn)。

在自動鋪帶工藝中，復(fù)合材料預(yù)浸帶在鋪帶頭的作用下沿著特定軌跡鋪貼到模具表面，如圖11所示。

圖11 自動鋪帶工藝原理圖Fig.11 Sketch of Automated Tape Placement

鋪帶軌跡需保證單條預(yù)浸帶不產(chǎn)生過大變形進(jìn)而發(fā)生鋪放褶皺，同時需保證相鄰兩條預(yù)浸帶滿足一定的鋪放間隙要求。通常情況下，相鄰兩條預(yù)浸帶間不允許產(chǎn)生覆蓋，且鋪放間隙大小不能超過2.5mm。

現(xiàn)有的商用鋪帶軟件一般基于測地線規(guī)劃預(yù)浸帶的鋪放路徑以保證帶料變形最小[19]。

由于模具曲面高斯曲率的影響，相鄰兩條預(yù)浸帶間難以保持完全平行，導(dǎo)致鋪放間隙問題的產(chǎn)生。

一自由曲面模具，如圖12所示。其尺寸為（1000×5000）mm。L1、L2…L7是自由曲面模具0°鋪層中所鋪放的預(yù)浸帶。預(yù)浸帶型號為5228/T700，寬度為150mm。

圖12 自由曲面模具上相鄰帶料間的鋪放間隙分布Fig.12 Distribution of Gap on Freeform Mould Surface

為保證帶料變形最小，L1、L2…L7的初始鋪放路徑沿著測地線。由于該自由曲面不可展，相鄰測地線間無法保證完全平行，導(dǎo)致鋪放間隙問題的產(chǎn)生。如圖所示，相鄰兩條帶料L6與L7之間末端的鋪放間隙大小接近10mm，這將嚴(yán)重影響最終產(chǎn)品的質(zhì)量。鋪放間隙的大小采用了文獻(xiàn)[20]中的方法進(jìn)行計算。

接下來部分將采用提出的四階方法構(gòu)造自由曲面上的等測地曲率曲線，以實(shí)現(xiàn)預(yù)浸帶在模具曲面上的微小轉(zhuǎn)向，消除帶料間過大的鋪放間隙。

采用等測地曲率曲線的目的在于將整條預(yù)浸帶需要的轉(zhuǎn)向均勻賦予到軌跡上的每一點(diǎn)，可以防止預(yù)浸帶在軌跡某一點(diǎn)處出現(xiàn)變形集中進(jìn)而產(chǎn)生鋪放褶皺，從而保證鋪帶產(chǎn)品的最終質(zhì)量。

顯然，調(diào)整后Li的鋪放路徑應(yīng)該在調(diào)整前Li鋪放路徑的右側(cè)才能消除Li與Li-1之間過大的鋪放間隙，由第四節(jié)內(nèi)容可知，Li的鋪放路徑上測地曲率符號為負(fù)。為獲得鋪放路徑上測地曲率的幅值，可采用折半搜索法，具體，如圖13所示。

圖13 折半查找流程圖Fig.13 Flowchart of the Bi-Section Search

在圖13 中，Upper_bound 為折半查找法的上邊界，Lower_bound 為折半查找法的下邊界，滿足|Upper_bound|＞|Lower_bound|，Δ1、Δ2—初始上邊界與下邊界，它們的數(shù)值，如表5所示。

表5 折半查找的初始上邊界與下邊界Tab.5 The Initial Upper Bound and Lower Bound for the Bisection Search

折半搜索法的停止條件為1.5mm ≤End_Gap≤2.5mm，其中End_Gap表示Li與Li-1之間末端的鋪放間隙大小。當(dāng)折半搜索完成時，Li的鋪放路徑上測地曲率記為kg（Li）。

折半查找完成后測地曲率及鋪放間隙的大小，如表6所示。

表6 折半查找完成后測地曲率及鋪放間隙的大小Tab.6 Value of Geodesic Curvature and Gap After the Search Procedure is Finished

調(diào)整鋪放路徑前、后相鄰預(yù)浸帶間鋪放間隙的分布情況，如圖14所示。

圖14 調(diào)整鋪放路徑前、后鋪放間隙的分布情況Fig.14 The Distribution of Gap on the Mould Surface

其中一條曲線表示優(yōu)化前鋪放間隙的分布，另一條曲線表示優(yōu)化后鋪放間隙的分布。

由該圖可看出，采用本方法所構(gòu)造出的等測地曲率曲線成功實(shí)現(xiàn)了預(yù)浸帶在模具曲面上的轉(zhuǎn)向，消除了帶料間過大的鋪放間隙。

優(yōu)化前后模具曲面上的鋪放路徑，如圖15所示。

圖15 優(yōu)化前后自由曲面模具上的鋪放路徑Fig.15 The Lay-Up Trajectories on the Freeform Surface

其中，一條曲線表示原始的鋪放路徑；另一條曲線表示調(diào)整后的鋪放路徑。

采用第三節(jié)中的四階方法進(jìn)行軌跡規(guī)劃時鋪放路徑L1、L2…L7的計算精度以及7條曲線總的計算時間消耗（未包含折半搜索的時間消耗），如圖16所示。

圖16 本方法與文獻(xiàn)[16]中方法的對比Fig.16 Comparison of Our Method with Reference[16]

由該圖可看出，與現(xiàn)有的數(shù)值方法相比，本方法在相當(dāng)?shù)挠嬎銜r間內(nèi)可以得到更好的計算精度，即計算效率更高。

具體的說，當(dāng)計算步長為10與1時，本方法平均計算誤差為2.3E-7 與3.47E-11，而文獻(xiàn)[16]中的方法平均計算誤差為7.3E-7與7.92E-11。

5.2 纖維纏繞

纖維纏繞工藝同樣是一種復(fù)合材料自動化加工工藝，一般應(yīng)用于壓力管道、壓力容器等回轉(zhuǎn)體的制造。

近年來，該工藝也逐漸被應(yīng)用于制造葉片、飛機(jī)S型進(jìn)氣道、三通管等非回轉(zhuǎn)體異型面[21]。

本方法可應(yīng)用于纖維纏繞工藝中非測地線纏繞軌跡的設(shè)計。非測地線纏繞對測地曲率的要求為：纏繞軌跡上各點(diǎn)的測地曲率與軌跡在該點(diǎn)處法曲率的比值不能超過纖維與模具之間的最大靜摩擦系數(shù)[22]。

現(xiàn)有的研究在處理異型面纏繞時通?；诰W(wǎng)格化方法，即先將原始曲面離散為由許多微小面片組成的網(wǎng)格曲面，然后在網(wǎng)格曲面上設(shè)計纏繞軌跡。

網(wǎng)格化方法的計算精度與離散網(wǎng)格的精細(xì)程度直接相關(guān)，當(dāng)網(wǎng)格劃分很精細(xì)時，網(wǎng)格曲面的生成需要消耗較多的計算時間。另外，該方法也難以保證設(shè)計出的軌跡完全位于原始曲面之上。

本方法不需要生成離散曲面，只要目標(biāo)曲面能夠通過參數(shù)方程進(jìn)行表達(dá)，就可以設(shè)計出高精度的非測地線纏繞軌跡。

采用本方法設(shè)計非測地線纏繞軌跡時，需對F算子中公式進(jìn)行調(diào)整，此時單位切向量t0在P0點(diǎn)相對弧長的變化率需通過下式計算：

式中：λ—纖維纏繞的滑移系數(shù)。

圖17顯示了一個基于自由曲面表達(dá)的異型管道，以及采用第三節(jié)中四階方法所設(shè)計出的不同滑移系數(shù)的非測地線纏繞軌跡，這些軌跡具有很高的計算精度且完全位于目標(biāo)曲面之上。

圖17 異型管道上不同滑移系數(shù)的非測地線纏繞軌跡Fig.17 Non-Geodesic Winding Paths on a Non-Revolution Surface

6 結(jié)束語

提出了一種參數(shù)曲面上具有規(guī)定測地曲率曲線的幾何方法，當(dāng)給出目標(biāo)曲面上的起始點(diǎn)、方向以及測地曲率的分布時，可追蹤出符合設(shè)計要求的空間曲線。

與現(xiàn)有數(shù)值方法相比，本方法無需建立以及求解反映具有規(guī)定測地曲率曲線的復(fù)雜數(shù)學(xué)模型，且對于不同的參數(shù)曲面具有統(tǒng)一的求解格式；與現(xiàn)有的幾何方法相比，本方法利用多點(diǎn)平均可以得到很高的計算精度。

實(shí)驗(yàn)結(jié)果表明，本方法在計算效率上要顯著高于現(xiàn)有的幾何法，且一般情況下要優(yōu)于數(shù)值解法。本方法可應(yīng)用于自動鋪帶軌跡規(guī)劃以及纖維纏繞工藝中非測地線纏繞路徑的設(shè)計。

為便于比較各方法的計算效率，文中未考慮在計算過程中動態(tài)調(diào)整計算步長。

如何合理利用每一步計算所得到的幾何信息設(shè)計變步長方法將是下一步的研究重點(diǎn)。