蔡 強(qiáng)，趙曉寧，李新田，淡林鵬，蒲曉航，黃慧慧

(中國(guó)運(yùn)載火箭技術(shù)研究院，北京 100076)

0 引言

纖維纏繞復(fù)合材料壓力容器(以下簡(jiǎn)稱復(fù)材容器)由于具有輕質(zhì)高強(qiáng)的優(yōu)點(diǎn)，在航天技術(shù)領(lǐng)域得到了廣泛的應(yīng)用[1-7]。復(fù)材容器一般由前后封頭和筒段組成，三者一體纏繞成型，在前后封頭纏繞過程中，纏繞角、纖維厚度和纖維應(yīng)力等參數(shù)是不斷變化的，不同的纏繞線型和子午線形狀直接影響復(fù)材容器的承壓能力，因此封頭設(shè)計(jì)及分析工作是復(fù)材容器設(shè)計(jì)的核心內(nèi)容。

復(fù)材容器封頭的設(shè)計(jì)分析多采用網(wǎng)格理論計(jì)算，內(nèi)壓完全由纖維承擔(dān)，樹脂基體僅起到黏合作用，纖維纏繞通常采用測(cè)地線纏繞理論，不用考慮摩擦力的因素，纖維能夠在芯模上穩(wěn)定鋪敷不打滑[8-10]。近年來，有關(guān)學(xué)者引入經(jīng)典層板理論和非測(cè)地線纖維纏繞軌跡，考慮了樹脂基體的承力并擴(kuò)大了纖維纏繞軌跡設(shè)計(jì)空間[11-19]。

在方案論證階段，通常需要快速開展封頭方案的迭代設(shè)計(jì)、優(yōu)化分析，因此亟需建立不同封頭的快速方案評(píng)估模型，并通過對(duì)比分析優(yōu)選出滿足設(shè)計(jì)要求的方案。為此，本文系統(tǒng)整理出螺旋纏繞封頭、螺旋纏繞橢球封頭、平面纏繞封頭和平面纏繞橢球封頭4種封頭的控制方程，設(shè)計(jì)上采用成熟可靠的網(wǎng)格理論和測(cè)地線理論，通過實(shí)例計(jì)算對(duì)比分析纏繞線型和封頭子午線類型對(duì)復(fù)材容器性能的影響。

1 復(fù)材容器參數(shù)模型

1.1 幾何參數(shù)模型

本文主要研究對(duì)象是復(fù)材容器封頭，但是筒段的長(zhǎng)度、外徑和壁厚直接影響封頭的參數(shù)，因此采用圖1所示的封頭和筒段一體化的幾何參數(shù)模型。封頭的主要參數(shù)有子午線形狀、纏繞角、封頭壁厚和纖維應(yīng)力。

圖1 復(fù)材容器參數(shù)模型Fig.1 Simplified parametric model of filament-wound composite pressure vessels

復(fù)材容器成型工藝是在芯模上按照一定的軌跡連續(xù)纏繞浸膠的復(fù)材纖維帶，常用的纖維纏繞線型軌跡有螺旋纏繞和平面纏繞2種方式。

螺旋纏繞線型如圖2所示，其中r1為給定的筒段徑向位置；α1為筒段纖維纏繞角；r2為給定的封頭徑向位置；α2為封頭纖維纏繞角。纖維軌跡嚴(yán)格滿足穩(wěn)定的測(cè)地線，即r1sinα1=r2sinα2，不存在滑紗的風(fēng)險(xiǎn)，同時(shí)封頭各點(diǎn)纖維應(yīng)力相等，但是要求前后極孔直徑相等。

圖2 螺旋纏繞線型Fig.2 Helical filament winding

平面纏繞線型如圖3所示，其中φFront pole為前極孔直徑，φBack pole為后極孔直徑，前后極孔直徑可以不等，α為纖維纏繞角。平面纏繞線型軌跡是從一個(gè)極孔直接纏繞至另一個(gè)極孔，這種線型只能用于長(zhǎng)度較小的殼體，要求殼體長(zhǎng)徑比L/D≤2。

圖3 平面纏繞線型Fig.3 Polar filament winding

以上2種線型，封頭和筒段一體纏繞，纖維應(yīng)力連續(xù)變化，承內(nèi)壓能力較好，但是封頭子午線形狀受限于纏繞線型無法自主設(shè)計(jì)，工程上為了設(shè)計(jì)和工藝方便，并提高體積裝填率，經(jīng)常采用橢球封頭，即封頭子午線方程滿足橢圓弧曲線方程

r2/R2+z2/B2=1

(1)

令橢球比

m=R/B

將不同纏繞線型和封頭子午線形狀組合，常用的封頭形式可以分為螺旋纏繞封頭、螺旋纏繞橢球封頭、平面纏繞封頭和平面纏繞橢球封頭4類。

1.2 網(wǎng)格理論假設(shè)

封頭設(shè)計(jì)及分析時(shí)采用如下的網(wǎng)格理論，有如下假設(shè)條件：

1)復(fù)材容器的載荷全部由纖維的拉伸強(qiáng)度承擔(dān)，樹脂基體僅對(duì)纖維起黏結(jié)、定位的作用，使纖維強(qiáng)度充分發(fā)揮；

2)在整個(gè)復(fù)材容器上纖維連續(xù)纏繞，分布均勻且對(duì)稱；

3)復(fù)材容器為薄壁結(jié)構(gòu)件，應(yīng)力沿壁厚均勻分布，無彎曲應(yīng)力，內(nèi)部壓力處處相等，均勻一致；

4)復(fù)材容器前后封頭和筒段螺旋向一體成型，然后筒段補(bǔ)充環(huán)向纏繞，即前后封頭為“螺旋向”或“平面”纏繞線型，筒段為“螺旋向+環(huán)向”或“平面+環(huán)向”纏繞線型。

1.3 筒段壁厚計(jì)算方法

筒段長(zhǎng)度L和外徑R一般根據(jù)約束條件給定，筒段壁厚hfc與封頭形式無關(guān)，本文采用參考文獻(xiàn)[20]中方法計(jì)算，為了方便下文引用，列出筒段壁厚

(2)

式中：pb為容器最大內(nèi)壓；hfα為纖維螺旋向或平面纏繞壁厚；hfθ為環(huán)向壁厚；σfb為纖維的實(shí)際發(fā)揮強(qiáng)度。

2 復(fù)合容器封頭設(shè)計(jì)模型

2.1 設(shè)計(jì)規(guī)則及控制方程

一般封頭的幾何構(gòu)型如圖4所示，封頭上采用螺旋纏繞或平面纏繞線型，纖維分布滿足以下3個(gè)規(guī)則：

規(guī)則1：纖維排列關(guān)于子午線對(duì)稱；

規(guī)則2：纖維與子午線的夾角即纏繞角α是關(guān)于平行圓半徑r的函數(shù)，α=α(r)，赤道圓上的纏繞角等于圓筒上的縱向纏繞角α0，在極孔上，α=π/2；

規(guī)則3：根據(jù)物質(zhì)守恒定律，通過平行圓法截面上的纖維總量等于通過赤道圓法截面上的纖維總量，且等于通過圓筒橫截面上的纖維總量。因平行圓半徑是變化的，因而封頭厚度hf也是r的函數(shù)，即hf=hf(r)，在赤道圓處，封頭壁厚等于筒段壁厚hfα，即hf=hfα。

圖4 一般封頭幾何構(gòu)型Fig.4 Geometry of general dome

極坐標(biāo)條件下，一般回轉(zhuǎn)體封頭的曲面方程

S(θ,z)={rcosθ,rsinθ,z}

(3)

式中：r=r(z)為封頭的母線方程，r是平行圓半徑坐標(biāo)；z為回轉(zhuǎn)體軸線坐標(biāo)；r和z均大于0，θ為沿緯度圓的包角。

引入無量綱參數(shù)

(4)

根據(jù)規(guī)則1～規(guī)則3，可以推導(dǎo)出一般封頭的無量綱化控制方程[20]

(5)

2.2 螺旋纏繞封頭

(6)

圖5 螺旋纏繞封頭子午線拐點(diǎn)及處理方法Fig.5 Inflection point and modified meridian of helical winding dome

2.3 螺旋纏繞橢球封頭

橢球封頭的子午線方程是給定的橢圓弧，將式(4)代入式(1)得出無量綱化的子午線方程

ρ2+m2ε2=1

(7)

對(duì)ρ求一階和二階導(dǎo)數(shù)，可得

(8)

(9)

值得注意的是該封頭不再滿足等應(yīng)力條件。

2.4 平面纏繞封頭

從圖3中可知，平面纏繞筒段纖維纏繞角α0滿足

(10)

根據(jù)圖6中的平面纏繞封頭的幾何關(guān)系，得出

ρsinθ=ρe+εtanα0

(11)

(12)

式中ρe為赤道處纖維軌跡起點(diǎn)坐標(biāo)在子午線平面的投影。

圖6 平面纏繞封頭幾何關(guān)系Fig.6 Geometry of polar winding dome

聯(lián)合式(5)、式(11)和式(12)，得出平面纏繞封頭的控制方程

(13)

與橢球封頭相似，平面封頭曲率在纏繞角α=54.7°時(shí)也會(huì)出現(xiàn)拐點(diǎn)，采用相同的處理方法，在等曲率點(diǎn)α=45°以后用圓弧替代原子午線。

2.5 平面纏繞橢球封頭

將式(8)中橢圓弧二階導(dǎo)數(shù)代入式(13)，得出平面纏繞橢球封頭的控制方程

(14)

3 實(shí)例及結(jié)果分析

3.1 模型驗(yàn)證

為了驗(yàn)證本文計(jì)算模型的正確性，選取極孔半徑ρ0=0.1、0.2、0.3、0.4共4種工況，分別對(duì)比螺旋纏繞封頭曲線與參考文獻(xiàn)[21]中的等應(yīng)力封頭曲線，計(jì)算結(jié)果如圖7所示，可見吻合性非常好。

3.2 實(shí)例對(duì)比分析

以下針對(duì)4種不同類型復(fù)材容器封頭的實(shí)例進(jìn)行討論，采用T700/樹脂基體復(fù)合材料，纖維體積含量為60%，復(fù)合密度為1.56g/cm3，最大設(shè)計(jì)壓強(qiáng)為10 MPa，筒段外徑均為φ1 000 mm。4個(gè)實(shí)例三維模型如圖8所示，其他幾何參數(shù)如表1所示。

圖7 封頭設(shè)計(jì)模型驗(yàn)證Fig.7 Verification of dome design model

圖8 設(shè)計(jì)實(shí)例三維模型Fig.8 3D solid geometric models of cases

表1 實(shí)例幾何尺寸

由表1可以看出，4個(gè)后封頭尺寸相同，以下對(duì)比4個(gè)后封頭的設(shè)計(jì)結(jié)果。封頭的子午線、纏繞角和壁厚的設(shè)計(jì)結(jié)果對(duì)比如圖9所示。

1)子午線設(shè)計(jì)結(jié)果：除螺旋纏繞封頭高度較大外，其余3種封頭子午線輪廓設(shè)計(jì)結(jié)果相近，螺旋纏繞封頭高度(332 mm)比其余封頭高度大了約8.9%。另外，平面纏繞封頭和螺旋纏繞封頭在等曲率點(diǎn)前后變化趨勢(shì)連續(xù)無奇異點(diǎn)，證明2.2節(jié)處理方法的正確性。

2)纏繞角設(shè)計(jì)結(jié)果：相同纏繞線型的纏繞角度變化趨勢(shì)相似，且變化連續(xù)，z=250 mm之前螺旋纏繞線型纏繞角(初始纏繞角約15°)大于平面纏繞線型(初始纏繞角約30°)，之后變化趨勢(shì)一致，由于螺旋纏繞封頭高度較大，因此在極孔附近纏繞角變化較緩。

3)封頭厚度設(shè)計(jì)結(jié)果：封頭厚度的變化趨勢(shì)同纏繞角，相同線型相似，且螺旋纏繞線型平均封頭厚度大于平面纏繞線型。

圖9 后封頭設(shè)計(jì)結(jié)果對(duì)比Fig.9 Comparison of dome design results

4種封頭的纖維應(yīng)力對(duì)比如圖10所示，可見：

圖10 后封頭纖維應(yīng)力設(shè)計(jì)結(jié)果對(duì)比Fig.10 Comparison of fiber stress design results for rear dome

工程上為了解決極孔附近應(yīng)力較大問題，通常在此處進(jìn)行局部補(bǔ)強(qiáng)，另外在極孔處通過纏繞固定鋁合金或鈦合金接頭的設(shè)計(jì)來彌補(bǔ)應(yīng)力突變帶來的結(jié)構(gòu)強(qiáng)度不足。

為比較4種后封頭的綜合性能，引入?yún)⒖嘉墨I(xiàn)[8]中的綜合性能系數(shù)pV/W(p為容器內(nèi)壓；V為封頭體積；W為封頭質(zhì)量)。詳細(xì)對(duì)比如表2所示。

1)由于封頭高度的原因，螺旋纏繞封頭自由容積最大，其余三者相近。

2)由于螺旋纏繞封頭平均厚度大于平面纏繞，相同纏繞線型的封頭質(zhì)量相近，螺旋纏繞線型封頭的質(zhì)量大于平面纏繞線型。

3)相同纏繞線型的綜合性能系數(shù)pV/W相近，平面纏繞線型的封頭的pV/W值優(yōu)于螺旋纏繞線型。

表 2 后封頭性能參數(shù)對(duì)比Tab.2 Comparison of performance parameters for rear dome

綜上，相同纏繞線型的纏繞角、壁厚和性能參數(shù)接近，平面纏繞線型封頭質(zhì)量輕、pV/W值高，螺旋纏繞線型纖維受力較好。為了避免封頭曲面曲率的不連續(xù)性，同時(shí)考慮設(shè)計(jì)和加工的便捷性，推薦采用螺旋纏繞橢球封頭或平面纏繞橢球封頭。

4 結(jié)論

本文基于復(fù)材容器的參數(shù)模型，系統(tǒng)研究了殼體筒段和封頭的設(shè)計(jì)方法，完成了螺旋纏繞封頭、螺旋纏繞橢球封頭、平面纏繞封頭和平面纏繞橢球封頭等4類封頭控制方程的建模，通過對(duì)4種后封頭實(shí)例的子午線、纏繞角、壁厚、纖維應(yīng)力以及綜合性能參數(shù)對(duì)比，發(fā)現(xiàn)相同纏繞線型的封頭設(shè)計(jì)結(jié)果相似，建議根據(jù)實(shí)際應(yīng)用情況采用螺旋纏繞橢球封頭或平面纏繞橢球封頭，以避免子午線曲率的拐點(diǎn)，同時(shí)有利于加工制造。