劉永佼,王顯峰,肖軍*

南京航空航天大學(xué) 材料科學(xué)與技術(shù)學(xué)院,南京 210016

圓錐面等測地曲率曲線的軌跡規(guī)劃方法

劉永佼,王顯峰,肖軍*

南京航空航天大學(xué) 材料科學(xué)與技術(shù)學(xué)院,南京 210016

為解決復(fù)合材料錐殼鋪放過程中的纖維殘余應(yīng)變過大的問題,提出了等測地曲率曲線的軌跡規(guī)劃方法。通過分析發(fā)現(xiàn)軌跡的測地曲率決定了纖維的殘余應(yīng)變,于是通過計算圓錐上任意曲線的測地曲率推導(dǎo)出了等測地曲率曲線的表達(dá)式,并證明在圓錐展開面上等測地曲率曲線是圓弧。最后,通過計算等測地曲率軌跡的殘余應(yīng)變及其與圓錐母線的夾角分析該軌跡鋪放工藝性。結(jié)果表明:等測地曲率軌跡相較于固定角軌跡與測地線軌跡具有良好的可設(shè)計性和鋪放工藝性。

復(fù)合材料;自動鋪絲;軌跡規(guī)劃;等測地曲率;圓錐面

復(fù)合材料因其輕質(zhì)、高強度、高模量、可設(shè)計性強、易于實現(xiàn)自動化整體化制造等特點而應(yīng)用于航空航天等尖端技術(shù)領(lǐng)域[1-6]。復(fù)合材料錐殼在航空航天領(lǐng)域應(yīng)用廣泛,如各類運載器與載荷的連接過渡段、固體火箭發(fā)動機噴管、人造衛(wèi)星的碳纖維纏繞錐殼,以及航空領(lǐng)域上如飛機的雷達(dá)罩,寬體客機的尾機身等[7-8]。自動鋪絲技術(shù)可實現(xiàn)大尺寸、復(fù)雜形狀復(fù)合材料結(jié)構(gòu)件的精密自動化制造,是復(fù)合材料低成本自動化制造技術(shù)的重要發(fā)展方向[9-14]。

現(xiàn)有復(fù)合材料結(jié)構(gòu),常以固定角度鋪層設(shè)計以滿足纖維按復(fù)合材料強度設(shè)計方向排布。一般采用0°鋪層、±45°鋪層和90°鋪層按照特定順序排布[15-16]。對于圓錐一類的可展曲面,采用解析法求解軌跡具有更高的精確度。圓錐構(gòu)件采用固定角軌跡鋪放時,鋪放軌跡偏離測地線,鋪放過程有制造殘余壓縮應(yīng)變產(chǎn)生,隨圓錐截面變小而變大,即從大端到小端應(yīng)變逐漸增大[17]。纖維壓縮應(yīng)變過大的區(qū)域鋪放質(zhì)量較差,進(jìn)而影響制件的性能[18]。對此,本文對圓錐的幾何特征進(jìn)行解析,規(guī)劃出等測地曲率的鋪放軌跡,使整個制件的殘余應(yīng)變均勻分布,避免產(chǎn)生性能薄弱區(qū)。

1 等測地曲率曲線求解

1.1 應(yīng)變與測地曲率

由于預(yù)浸料中增強纖維彈性模量非常高,鋪放過程中纖維縱向拉伸變形能力有限,而纖維縱向受壓縮作用產(chǎn)生的微屈曲是預(yù)浸帶變形的主要機制。根據(jù)文獻(xiàn)[19]可得在平面上纖維鋪放的殘余應(yīng)變?yōu)?/p>

式中:w為帶寬;kg為測地曲率。式(1)表明,當(dāng)所使用的預(yù)浸帶寬度w一定時,平面上鋪放軌跡測地曲率越小,纖維應(yīng)變絕對值越小,軌跡的可鋪放性越好。

實際鋪放的型面多為曲面,以自由曲面為對象研究其上任意一條曲線的特征,如圖1所示,e1、e2、e3為單位向量。C為曲面上一條曲線,P為曲線上的任意一點,k為曲線在P點的曲率向量,t為曲面在P點的切向量,n為曲面在P點的法向量。根據(jù)微分幾何理論,曲線上的任意一點的曲率向量k可以分解成兩個部分:一個在曲面法線n的方向上,另一個在與n和曲線在該點切向量t都垂直的方向上:

式中:kn為曲線的法向曲率。kn是與曲面本身有關(guān)的曲率,它只會使預(yù)浸帶產(chǎn)生厚度方向的縱向彎曲變形,由于預(yù)浸料很薄(0.125～0.2 mm),不會對軌跡的可鋪放性產(chǎn)生影響。而測地曲率kg的存在則會使預(yù)浸帶產(chǎn)生側(cè)向彎曲變形進(jìn)而影響可鋪放性,所以對于任意的自由曲面的鋪放軌跡可將其測地曲率的大小作為評判軌跡纖維應(yīng)變大小的依據(jù)[20]。由此可以推出等測地曲率曲線上纖維殘余應(yīng)變處處相等。

1.2 等測地曲率曲線計算方法

如圖2所示,以圓錐的頂點為原點,圓錐的旋轉(zhuǎn)軸為Z軸建立直角坐標(biāo)系。圓錐面可以表示為XOZ平面內(nèi)一條直線l繞Z軸旋轉(zhuǎn)一周掃過的空間,該直線過O點且不與X 軸、Z軸垂直,l與Z軸的夾角θ為圓錐的半頂角。圓錐面上任意一點P的位置可以通過u、v兩個變量表示,其中u表示l繞Z軸旋轉(zhuǎn)的角度,v表示P點到O點的距離。點P的位置矢量OP記為r,則圓錐面可以表示為矢量函數(shù)r(u,v)的值域。該函數(shù)可以寫成直角坐標(biāo)系的形式r(x,y,z),兩者的轉(zhuǎn)換關(guān)系為

根據(jù)微分幾何的知識可知,空間曲線C可以寫成位置矢量函數(shù)r(s)的形式,s為該點到曲線起始點之間的弧長。當(dāng)曲線C位于圓錐面上時,其位置矢量函數(shù)還需滿足圓錐面的函數(shù)表達(dá)式,所以圓錐面上的曲線C的表達(dá)式為r(u(s),v(s))。由于該曲線具有等測地曲率的性質(zhì),從該性質(zhì)入手可以推導(dǎo)出u(s)和v(s),代入式(3)即可得到曲線在笛卡兒坐標(biāo)系下的參數(shù)方程。

微分幾何[21]中給出了測地曲率的計算公式為kg=r″(s)·(n(s)×r′(s)) (4)式中:r′(s)和r″(s)分別為r(s)的一階和二階導(dǎo)數(shù);n(s)為圓錐面的單位法向量。對于任意曲面r(u,v),其單位法向量n(u,v)的計算公式為

式中:ru(u,v)和rv(u,v)分別為r(u,v)關(guān)于u和v的一階偏導(dǎo)數(shù)。將式(3)代入圓錐面的直角坐標(biāo)表達(dá)式r(x,y,z)得到:

對式(6)分別求關(guān)于u、v的一階偏導(dǎo)數(shù),并代入式(5)和u=u(s)得到圓錐面的單位法向量:

因為曲線C位于圓錐面上,因此曲線C參數(shù)代入圓錐面參數(shù)方程依然成立。將u=u(s)、v=v(s)代入式(6)得到曲線C的參數(shù)方程:

對r(s)分別求一階導(dǎo)數(shù)和二階導(dǎo)數(shù),并和式(7)一并代入式(4)中可得:

式(8)無法直接得出u(s)和v(s)滿足何種條件時測地曲率kg為常數(shù),不妨對特定的u(s)進(jìn)行討論,找出盡可能簡單的曲線。

當(dāng)u′(s)=0時,則u(s)為常數(shù),u″(s)=0,此時kg恒等于0。由u(s)為常數(shù)可知曲線C是圓錐的母線,同時也是圓錐的測地線,測地線測地曲率恒為0,也是等測地曲率曲線。

當(dāng)u′(s)≠0且u″(s)=0時,不妨設(shè)u=as+b則:若要kg為常數(shù)則要通過多項式求和消去式(9)中v、v′、v″中的含s項,式中多項式各項v、v′、v″的總次數(shù)都是2,若v=v(s)對于s求導(dǎo)不改變含s項的次數(shù),則可通過調(diào)整各項系數(shù)來消去含s項。常見函數(shù)中滿足此條件的有兩種,v=m ens和v=m cos(cs)+n sin(cs)(為簡化計算設(shè)常數(shù)項為0,m、n、c為待定系數(shù))。分別代入式(9)可得:

首先分析式(10),若要消除e2ns中s的影響,則需要m=0或n=0(a≠0且圓錐sinθ＞0)。對于v=m ens而言,m=0時v=0,曲線C為圓錐的頂點,不符合等測地曲率曲線的要求;當(dāng)n=0時,v=m,此時曲線C為圓錐的截面圓,圓錐的截面圓是等測地曲率曲線。

接下來分析式(11),含s的項都是三角函數(shù)的二次項,所以最佳的消除方式是通過正余弦平方和為1的方式消除三角函數(shù)項。因此需要將系數(shù)配平,最直接的方法是令c=a sinθ,代入式(11)得:

式(12)中各參數(shù)均為常數(shù),此種條件下曲線C是等測地曲率曲線。為進(jìn)一步確定曲線C的幾何信息,則要盡可能確定式(12)中的未知參數(shù),或者明確他們所代表的幾何意義。由圖2的圓錐坐標(biāo)圖和式(6)的圓錐表達(dá)式,很容易得到在此坐標(biāo)系下曲線C的弧長微分表達(dá)式:

將u=as+b和v=m cos(cs)+n sin(cs)以及c=a sinθ代入式(13),并通過等式左右相等可得線C的表達(dá)式為

1.3 等測地曲率曲線的幾何特征

對式(13)作變量代換,令x=v cos(u sinθ),y=v sin(u sinθ),則式(13)變形為

式(15)是平面上弧長的微分表達(dá)式,從式(14)空間弧長表達(dá)式變換到式(15)平面弧長表達(dá)式的過程中,d s并未發(fā)生變化,即保長變換。該變換過程對應(yīng)的幾何意義是將圓錐的側(cè)面展開成平面扇形的過程。接下來討論曲線C在展開后的平面上具有什么樣的幾何特征。

為了計算簡便,設(shè)式(14)中的b=0且v≠0,則由式(14)可得v=m cos(u sinθ)+n sin(u sinθ),代入x=v cos(u sinθ),y=v sin(u sinθ)并整理可得

從式(16)可以看出曲線C在展開圖上是一個過原點(圓錐頂點)的圓,其半徑為其測地曲率等于該圓的曲率。由于過頂點的圓在軌跡規(guī)劃時有種種的局限性,所以該結(jié)論工程意義有限。為了盡可能滿足軌跡規(guī)劃的要求,要將該結(jié)論推廣到展開圖上一般圓上,即圓錐展開圖上的任意圓弧都是等測地曲率曲線。欲證明該結(jié)論,需要大量的微分幾何知識。圓錐曲面r(u,v)的全微分d r=rud u+rvd v,簡記為d r=rαdα,式中:rα=?r(u,v)/?α;α=u,v。后面出現(xiàn)的希臘字母α、β、γ、τ等若不特殊說明,作為變量時取值需要分別遍歷u、v。根據(jù)微分幾何曲面論給出r(u,v)的第一基本形式和第二基本形式:I=d r·d r=(rα·rβ)dαdβ,II=d2r·n=(rαβ·n)dαdβ,其中n,這樣,曲面圓錐曲面的兩個基本形式為

而第一類基本量gαβ和第二類基本量bαβ可以分別寫成一個2階矩陣,將矩陣記為(gαβ)。根據(jù)曲面論中第一基本形式的性質(zhì)可知,矩陣(gαβ)是正定矩陣,因此(gαβ)存在逆矩陣,記為(gαβ)。第一基本形式還具備另外一個重要性質(zhì):I與參數(shù)選擇無關(guān)。即當(dāng)x=x(u,v),y=y(u,v)時,存在gˉαβ使得I=gˉαβdαdβ。因此從圓錐曲面變換r(u,v)到平面r-(u,v)的過程中第一基本形式I是不變的,第二基本形式II是不具備該性質(zhì)的,即II與參數(shù)選擇相關(guān)。而下面要證明測地曲率kg的值與參數(shù)選擇無關(guān),即kg中不含第二類基本量。

由定義可得ru·rv=0,結(jié)合式(5)可知{r ;ru,rv,n} 是空間中的正交標(biāo)架,即ru、rv、n是線性無關(guān)的,而?rα/?β仍然是空間中的向量,因此不妨假定

即系數(shù)Dαβ恰好是圓錐曲面的第二類基本量bαβ。在式(18)兩邊用切向量rτ作內(nèi)積,則得

函數(shù)r(u,v)的兩次偏導(dǎo)數(shù)與求導(dǎo)的次序無關(guān),因此rαβ關(guān)于下指標(biāo)α、β是對稱的,于是記號Cαγβ和Cγαβ關(guān)于下指標(biāo)α、β是對稱的。對gαβ=rα·rβ求偏導(dǎo)數(shù)得到并代入式(21)得到

將下標(biāo)調(diào)換可得到

將這兩式相加,再減去式(22),并且利用Cγαβ關(guān)于后兩個下標(biāo)的對稱性,則得

因此

由此可見,式(19)的系數(shù)函數(shù)Cγαβ是由圓錐曲面的第一類基本量及其一階偏導(dǎo)數(shù)決定的,因為第一基本形式由第一類基本量決定,且第一基本形式與參數(shù)選取無關(guān),所以結(jié)合微分形式不變性可知,Cαγβ也是與參數(shù)選取無關(guān)的函數(shù)。

設(shè)圖1中的自由曲面為圓錐曲面r(u,v),其上的任意一條曲線u=u(s)、v=v(s)在空間中的參數(shù)方程為

r(s)=r(u(s),v(s))

沿曲線C建立正交標(biāo)架場{r;e1,e2,e3},所以

結(jié)合式(2)和圖1以及微分幾何中對曲線曲率的定義可得:

在該式中代入式(18)可得

由于e3(s)=n,所以n⊥e2(s),因此

但是

因此,rγ·e2(s)的結(jié)果是rγ、n、e1(s)的混合積,混合積的幾何意義是三個向量構(gòu)成的平行六面體的體積,由于n是單位向量,故

因此

前文已經(jīng)證明Cαγβ只與曲面的第一類基本量及其偏導(dǎo)數(shù)有關(guān),因此kg只依賴于曲面的第一類基本量和曲線的參數(shù)方程,與參數(shù)選擇無關(guān)。由于kg中不含有曲面第二類基本量,因此曲面作保長變換x=v cos(u sinθ)、y=v sin(u sinθ)時,測地曲率kg不發(fā)生變化。

由于圓錐曲面展開成平面扇形的過程kg不變,而平面上有kn=0,k=kg·e2(s),e2(s)為平面內(nèi)的單位向量,即平面曲線的測地曲率等于其曲率。平面上任意圓的曲率是常數(shù),因此圓錐展開圖上的任意圓弧在空間中都是等測地曲率曲線。

2 等測地曲率曲線工藝性分析

以錐角為60°的圓錐為例,圖3是小端直徑230 mm,大端直徑577 mm,長300 mm的一個圓錐面,紅色軌跡(左側(cè))是固定角度軌跡,與圓錐面母線的夾角始終為30°。藍(lán)色軌跡(右側(cè))是以固定角軌跡的兩端點為起始點和終止點,以圓錐曲面為支持面所做的測地線。綠色軌跡(中間)是介于固定角軌跡和測地線軌跡之間的一條等測地曲率曲線。等測地曲率曲線是三點確定的圓弧,分別是固定角軌跡的兩端點和固定角與測地線之間的一點,本例中第三點取固定角軌跡中點與測地線軌跡中點的連線的中點。

圖3中各軌跡為自動鋪放時的參考軌跡,即中心軌跡。以自動鋪絲為例,計算寬度為6.35 mm的單根預(yù)浸紗沿各軌跡鋪放時產(chǎn)生的壓縮應(yīng)變,其中沿固定角度軌跡鋪放的壓縮應(yīng)變記為εa,沿等測地曲率軌跡鋪放的壓縮應(yīng)變記為εc,沿測地線軌跡鋪放的壓縮應(yīng)變記為εg,圓錐的截面直徑記為d。從小端到大端各部分應(yīng)變值如表1所示。從表中可以看出,固定角軌跡在小端附近壓縮應(yīng)變很大。實際鋪放時纖維不產(chǎn)生宏觀褶皺變形的臨界應(yīng)變值受很多工藝參數(shù)影響,根據(jù)文獻(xiàn)[18]可知在殘余壓縮應(yīng)變達(dá)到0.5%時纖維的力學(xué)性能已有明顯下降,本文暫且以0.5%為標(biāo)準(zhǔn)進(jìn)行衡量,可以看出按固定角軌跡鋪放時,小端附近的鋪放質(zhì)量很差。從表1中可知,沿測地線鋪放纖維不產(chǎn)生應(yīng)變;沿等測地曲率曲線鋪放時纖維產(chǎn)生一個較小的應(yīng)變,在鋪放容許的范圍之內(nèi)。從鋪放質(zhì)量來看,測地線優(yōu)于等測地曲率曲線,固定角軌跡有時不能滿足鋪放要求。

表1 小端至大端纖維壓縮應(yīng)變Table 1 Fiber compressive strain from small end to large end of cone

雖然沿測地線鋪放的鋪放質(zhì)量很好,但其鋪放方向與固定角度是有較大偏差的。表2計算了三種軌跡上從小端到大端各點處軌跡與圓錐母線的夾角,其中固定角軌跡與圓錐母線的夾角記為φa,等測地曲率軌跡與圓錐母線的夾角記為φc,測地線軌跡與圓錐母線的夾角記為φg。

從表2中可以看出,沿測地線鋪放時,纖維方向與最初設(shè)計有較大偏差,其最大偏差已超過10°,這對制件的力學(xué)性能會有較大影響[22]。等測地曲率曲線在角度偏離上要小于測地線。

由此可以看出,等測地曲率曲線是一種在方向性和鋪放性上均介于固定角與測地線之間的一種曲線,即在具有較小的纖維應(yīng)變的情況下同時還能兼顧一定的方向性。等測地曲率曲線的位置可以通過圓弧的第三點來進(jìn)行調(diào)節(jié),使其更加靠近固定角曲線或是測地線,即更加注重方向性或是可鋪放性。因此在實際鋪放中可根據(jù)具體情況采用等測地曲率曲線來替代固定角曲線或測地線進(jìn)行鋪放。

為驗證上述分析的合理性,以航天材料及工藝研究所提供的T300/603環(huán)氧樹脂預(yù)浸料為原料,以平面?zhèn)葟濅伔艑嶒灤驽F面鋪放實驗進(jìn)行驗證。式(1)指出纖維的壓縮應(yīng)變受預(yù)浸料寬度和測地曲率影響,1.3節(jié)證明了錐面曲線的測地曲率與展開面上圓弧曲率等價,因此可以用平面?zhèn)葟潓嶒瀬硖娲F面鋪放實驗。預(yù)浸料鋪放工藝條件如表3所示。

表2 小端至大端纖維角度Table 2 Fiber angle from small end to the large end of cone

結(jié)合表1和式(1)可以計算出應(yīng)變值對應(yīng)的曲率半徑,固定角應(yīng)變值最大值和最小值對應(yīng)的曲率半徑分別是505 mm和1 055 mm,等測地曲率應(yīng)變值對應(yīng)的曲率半徑為1 545 mm。圖4給出了不同曲率半徑下的寬度為6.35 mm單絲鋪放結(jié)果,當(dāng)圓弧半徑R=6 350 mm時,預(yù)浸料表面光滑未出現(xiàn)褶皺,當(dāng)圓弧半徑R=1 270 mm時,預(yù)浸料開始出現(xiàn)少許皺褶,此時纖維的壓縮應(yīng)變?yōu)?.5%;當(dāng)圓弧半徑R=790 mm時,預(yù)浸料尚未出現(xiàn)脫粘,此時應(yīng)變值為0.8%,當(dāng)圓弧半徑R=580 mm時,預(yù)浸料已經(jīng)出現(xiàn)脫粘,此時應(yīng)變?yōu)?%。由此可以得出,沿固定角鋪放時從大端到小端,纖維褶皺會逐漸增多直至脫粘;而沿等測地曲率曲線鋪放時預(yù)浸料表面光滑無褶皺,鋪放效果較好。

表3 鋪放工藝參數(shù)Table 3 Laying process parameters

3 結(jié) 論

1)為解決圓錐面纖維鋪放的殘余應(yīng)變問題,提出了等測地曲率軌跡規(guī)劃方法,并以測地曲率作為殘余應(yīng)變的評判標(biāo)準(zhǔn)。探究了圓錐曲面上曲線的測地曲率與曲線方程之間的關(guān)系,并推導(dǎo)出展開圖上過圓錐頂點的圓具有等測地曲率曲線的性質(zhì)。

2)證明了圓錐曲面展開面上的任意圓在空間中都是等測地曲率曲線,使得在圓錐面上采用等測地曲率軌跡規(guī)劃成為可能。值得注意的是,該證明過程具有普適性,對于任意的可展曲面都可以采用相同的證明方法。即任意可展曲面的展開圖上的圓在空間中都是等測地曲率曲線。

3)通過計算比較了等測地曲率曲線和固定角曲線與測地線的鋪放工藝性優(yōu)劣。等測地曲率曲線纖維的殘余應(yīng)變與固定角軌跡相比更小,鋪放角度的偏差比測地線更小。等測地曲率曲線在保證較小的纖維應(yīng)變的同時還能一定程度滿足纖維方向的要求,適用于圓錐曲面的小角度鋪放。

[1] SHIRINZADEH B,CASSIDY G,OETOMO D,et al.Trajectory generation for open-contoured structures in robotic fiber placement[J].Robotics and Computer-Integrated Manufacturing,2007,23(4):380-394.

[2] 肖軍,李勇,李建龍.自動鋪放技術(shù)在大型飛機復(fù)合材料構(gòu)件制造中的應(yīng)用[J].航空制造技術(shù),2008(1):50-53.XIAO J,LI Y,LI J L.Application of automate placement technology in fabrication of composite components for large aircraft[J].Aeronautical Manufacturing Technology,2008(1):50-53(in Chinese).

[3] 唐見茂.航空航天復(fù)合材料發(fā)展現(xiàn)狀及前景[J].航天器環(huán)境工程,2013,30(4):352-359.TANG J M.Current status and trends of advanced composites in aerospace[J].Spacecraft Environmental Engineering,2013,30(4):352-359(in Chinese).

[4] 肖軍,李勇,文立偉,等.樹脂基復(fù)合材料自動鋪放技術(shù)進(jìn)展[J].中國材料進(jìn)展,2009,28(6):28-32.XIAO J,LI Y,WEN L W,et al.Progress of automate placement technology for polymer composite[J].Material China,2009,28(6):28-32(in Chinese).

[5] 梁憲珠,張鋮.淺談降低國內(nèi)航空復(fù)合材料結(jié)構(gòu)成本的途徑[J].航空制造技術(shù),2011(3):40-44.LIANG X Z,ZHANG C.Approaches to reducing cost of composites in domestic aviation[J].Aeronautical Manufacturing Technology,2011(3):40-44(in Chinese).

[6] 朱晉生,王卓,歐峰.先進(jìn)復(fù)合材料在航空航天領(lǐng)域的應(yīng)用[J].新技術(shù)工藝,2012(9):76-79.ZHU J S,WANG Z,OU F.Applications of advanced composite materials in aerospace[J].New Technology&New Process,2012(9):76-79(in Chinese).

[7] 林勝.自動鋪帶機/鋪絲機(ATL/AFP)—現(xiàn)代大型飛機制造的關(guān)鍵設(shè)備(上)[J].世界制造技術(shù)與裝備市場,2009(4):84-89.LIN S.ATL/AFP—The key machine for manufacturing of modern large airplane(A)[J].World Manufacturing Engineering&Market,2009(4):84-89(in Chinese).

[8] 林勝.自動鋪帶機/鋪絲機(ATL/AFP)—現(xiàn)代大型飛機制造的關(guān)鍵設(shè)備(中)[J].世界制造技術(shù)與裝備市場,2009(5):89-95.LIN S.ATL/AFP—The key machine for manufacturing of modern large airplane(B)[J].World Manufacturing Engineering&Market,2009(5):89-95(in Chinese).

[9] 林勝.自動鋪帶機/鋪絲機(ATL/AFP)—現(xiàn)代大型飛機制造的關(guān)鍵設(shè)備(下)[J].世界制造技術(shù)與裝備市場,2009(6):78-83.LIN S.ATL/AFP—The key machine for manufacturing of modern large airplane(C)[J].World Manufacturing Engineering&Market,2009(6):78-83(in Chinese).

[10] 富宏亞,韓振宇,路華.纖維纏繞/鋪帶/鋪絲成型設(shè)備的發(fā)展?fàn)顩r[J].航空制造技術(shù),2009(22):43-46.FU H Y,HAN Z Y,LU H.Development of the filament winding/tape layer/fiber placement machine[J].Aeronautical Manufacturing Technology,2009(22):43-46(in Chinese).

[11] 張振甫,肖軍,吳海橋,等.復(fù)合材料錐殼0°鋪層的自動鋪放成型方法研究[J].宇航材料工藝,2007(2):55-57.ZHANG Z F,XIAO J,WU H Q,et al.Study on laminating method for 0°plies in composite conical shell manufactured by automated tape laying[J].Aerospace Materials&Technology,2007(2):55-57(in Chinese).

[12] 還大軍.復(fù)合材料自動鋪放CAD/CAM關(guān)鍵技術(shù)研究[D].南京:南京航空航天大學(xué),2010.HUAN D J.Fundamental research on CAD/CAM technology of composite material automated placement[D].Nanjing:Nanjing University of Aeronautics and Astronautics,2010(in Chinese).

[13] 李善緣,王小平,朱麗君.復(fù)合材料鋪絲成型中的路徑規(guī)劃[J].宇航材料工藝,2009(2):25-29,41.LI S Y,WANG X P,ZHU L J.Path planning for composite fiber placement[J].Aerospace Materials&Technology,2009(2):25-29,41(in Chinese).

[14] 朱麗君.復(fù)合材料自動鋪絲技術(shù)中的路徑計算[D].南京:南京航空航天大學(xué),2010.ZHU L J.Trajectory computation on robotic fiber placement for composite structures[D].Nanjing:Nanjing University of Aeronautics and Astronautics,2010(in Chinese).

[15] 盧敏,周來水,王小平,等.圓筒狀構(gòu)件的多層鋪絲路徑生成算法[J].航空學(xué)報,2011,32(1):181-186.LU M,ZHOU L S,WANG X P,et al.Trajectory generation for cylindrical structures in robotic multi-fiber placement[J].Acta Aeronautica et Astronautica Sinica,2011,32(1):181-186(in Chinese).

[16] 韓振宇,李玥華,富宏亞,等.錐殼零件自動鋪絲變角度軌跡規(guī)劃算法[J].計算機輔助設(shè)計與圖形學(xué)學(xué)報,2012,24(3):400-405.

HAN Z Y,LI Y H,FU H Y,et al.Variable-angles trajectory planning algorithm of automated fiber placement for conical shell[J].Journal of Computer-Aided Design&Computer Graphics,2012,24(3):400-405(in Chinese).

[17] 黃威,王顯峰,肖軍.復(fù)合材料錐殼固定角軌跡的鋪放工藝性分析[J].玻璃鋼/復(fù)合材料,2016(3):27-30.HUANG W,WANG X F,XIAO J.Analyze of placement manufacturability of fixed-angle trajectory for composite conical shell[J].Fiber Reinforced Plastics/Composites,2016(3):27-30(in Chinese).

[18] ZHAO C,WANG B D,XIAO J.Macroscopic characterization of fiber micro-buckling and its influence on composites tensile performance[J].Reinforced Plastic&Composite,2017,36(3):196-205.

[19] 李俊斐.基于結(jié)構(gòu)設(shè)計的開孔曲面自調(diào)節(jié)鋪絲軌跡規(guī)劃算法研究[D].南京:南京航空航天大學(xué),2013.

LI J F.Research on automated fiber placement trajectory planning method of adjustment algorithm based on structural design for surface with holes[D].Nanjing:Nanjing University of Aeronautics and Astronautics,2013(in Chinese).

[20] 熊文磊,肖軍,王顯峰,等.基于網(wǎng)格化曲面的自適應(yīng)自動鋪放軌跡算法[J].航空學(xué)報,2013,34(2):434-441.

XIONG W L,XIAO J,WANG X F,et al.Algorithm of adaptive path planning for automated placement on meshed surface[J].Acta Aeronautica et Astronautica Sinica,2013,34(2):434-441(in Chinese).

[21] 陳維桓.微分幾何[M].北京:北京大學(xué)出版社,2006:229-230.

CHEN W H.Differential geometry[M].Beijing:Peking University Press,2006:229-230(in Chinese).

[22] 李順林,王興業(yè).復(fù)合材料結(jié)構(gòu)設(shè)計基礎(chǔ)[M].武漢:武漢理工大學(xué)出版社,1993:15-25.

LI S L,WANG X Y.Foundation of composite structure design[M].Wuhan:Wuhan University of Technology Press,1993:15-25(in Chinese).

Trajectory planning method for constant geodesic curvature curve of cone

LlU Yongjiao,WANG Xianfeng,XlAO Jun*

College of Material Science and Technology,Nanjing University of Aeronautics and Astronautics,Nanjing 210016,China

To solve the problem of excessive residual strain of fiber in conical shell laying process,a trajectory planning method with constant geodesic curvature is proposed.Analysis shows that the geometric curvature of the trajectory determines the residual strain of the fiber.The expression for the constant curvature curve is then derived by calculating the geometrical curvature of the arbitrary curve on the cone,and the constant geodesic curvature curve is proved to be a circular arc on the conical surface.The trajectory laying process is analyzed by calculating the residual strain of the constant geodesic curvature trajectory and its angle with the conical busbar.The results show that the constant geodesic curvature trajectory is better designability and laying processability than the fixed angle trajectory and the geodesic trajectory.

composite material;automated fiber placement;trajectory planning;constant geodesic curvature;conical surface

2016-11-01;Revised:2016-11-29;Accepted:2016-12-27;Published online:2017-02-17 16:14

URL:www.cnki.net/kcms/detail/11.1929.V.20170217.1614.014.html

s:Major National Science and Technology Special Projects(2016ZX04002-001);A Project Funded by the Priority Academic Program Development of Jiangsu Higher Education lnstitutions

V261;TB332

A

1000-6893(2017)07-420904-08

10.7527/S1000-6893.2016.420904

2016-11-01;退修日期:2016-11-29;錄用日期:2016-12-27;網(wǎng)絡(luò)出版時間:2017-02-17 16:14

www.cnki.net/kcms/detail/11.1929.V.20170217.1614.014.html

“高檔數(shù)控機床與基礎(chǔ)制造裝備”科技重大專項 (2016ZX04002-001);江蘇高校優(yōu)勢學(xué)科建設(shè)工程資助項目

*通訊作者.E-mail:j.xiao@nuaa.edu.cn

劉永佼,王顯峰,肖軍.圓錐面等測地曲率曲線的軌跡規(guī)劃方法[J].航空學(xué)報,2017,38(7):420904.LlU Y J,WANG X F,XlAO J.Trajectory planning method for constant geodesic curvature curve of cone[J].Acta Aeronautica et Astronautica Sinica,2017,38(7):420904.

(責(zé)任編輯:李世秋)

*Corresponding author.E-mail:j.xiao@nuaa.edu.cn