王宇東,馬保吉
(西安工業(yè)大學(xué)兵器科學(xué)與技術(shù)學(xué)院,西安 710021)
光纖制導(dǎo)技術(shù)綜合運(yùn)用了成像探測(cè)、圖像跟蹤、特種光纖傳輸?shù)汝P(guān)鍵技術(shù)[1],具有抗干擾能力強(qiáng)、命中精度高等優(yōu)勢(shì),是精確制導(dǎo)技術(shù)的重要發(fā)展方向[2]。光纖線包是指制導(dǎo)光纖在線軸上纏繞完成后形成的紡錘形實(shí)體,其半剖截面呈梯形結(jié)構(gòu)[3]。光纖釋放之前要求在戰(zhàn)場(chǎng)復(fù)雜的運(yùn)輸環(huán)境下不松動(dòng)坍塌,在導(dǎo)彈飛行過程中光纖能夠順利釋放、不纏結(jié)、不折斷?;谏鲜鲆?,在光纖纏繞過程中需要涂覆粘結(jié)劑,然后進(jìn)行高溫固化,粘結(jié)劑固化后能使光纖之間緊密結(jié)合,避免松動(dòng)、坍塌[4]。
由于線包厚度較厚且材料傳熱性能差,在固化過程中線包表面和內(nèi)部的溫度場(chǎng)和固化度場(chǎng)分布差異較大,導(dǎo)致線包發(fā)生不同程度的熱膨脹和固化收縮,從而引起熱應(yīng)力和固化收縮應(yīng)力。這些殘余應(yīng)力會(huì)使得線包內(nèi)部出現(xiàn)應(yīng)力集中、裂縫等缺陷,這會(huì)導(dǎo)致光纖釋放過程中出現(xiàn)斷纖、脫匝、纏結(jié)等現(xiàn)象,從而導(dǎo)致放線失敗,導(dǎo)彈失去控制,因此如果能夠?qū)€包內(nèi)部缺陷進(jìn)行有效檢測(cè)將會(huì)大大提高導(dǎo)彈的可靠性。徐振華等[5]曾試圖通過X射線掃描的方法對(duì)線包進(jìn)行纏繞缺陷檢測(cè),但是由于光纖由不同材料組成,不同材料對(duì)X射線的吸收率存在較大差異,分辨率很難保證,并且射線的衰減導(dǎo)致其很難穿透整個(gè)線包。此外,關(guān)于線包固化殘余應(yīng)力的仿真研究也很困難,因?yàn)榫€包結(jié)構(gòu)復(fù)雜,光纖間隙涂抹有粘結(jié)劑,無法建立完整的光纖和粘結(jié)劑幾何模型以及光纖和粘結(jié)劑之間的接觸模型。
通過分析線包結(jié)構(gòu)發(fā)現(xiàn)其與復(fù)合材料結(jié)構(gòu)類似,且其結(jié)構(gòu)特性符合復(fù)合材料細(xì)觀力學(xué)基本假設(shè)[6],文中擬采用復(fù)合材料細(xì)觀力學(xué)理論對(duì)線包固化過程進(jìn)行模擬。復(fù)合材料細(xì)觀力學(xué)理論在復(fù)合材料固化過程仿真研究的應(yīng)用已非常成熟,文獻(xiàn)[7-9]都曾借助于細(xì)觀力學(xué)理論對(duì)復(fù)合材料固化過程進(jìn)行仿真研究,但仿真多以簡(jiǎn)單的層合板為對(duì)象,而線包結(jié)構(gòu)是一個(gè)紡錘形實(shí)體且其與線軸是不可分割的整體,在固化仿真過程中需同時(shí)考慮線包和線軸不同的材料屬性、接觸屬性,這些都將成為線包固化過程仿真的難點(diǎn)所在。
基于以上分析,將基于復(fù)合材料細(xì)觀力學(xué)理論建立線包固化過程的熱傳導(dǎo)數(shù)學(xué)模型,與固化動(dòng)力學(xué)數(shù)學(xué)模型進(jìn)行耦合,利用有限元分析軟件ANSYS建立線包和線軸的三維有限元模型、等效材料模型,設(shè)置初始條件和邊界條件,對(duì)耦合模型進(jìn)行數(shù)值求解。分析固化過程中線包不同位置的溫度和固化度分布,并研究了固化溫度和升溫速率對(duì)線包固化均勻性的影響。研究結(jié)果將有助于明確固化機(jī)理,根據(jù)特定的固化均勻性要求設(shè)計(jì)最佳固化工藝參數(shù)。
線包的固化工藝是將繞線和涂膠完成的線包放置在烘箱中進(jìn)行加熱處理。線包固化工藝曲線如圖1所示,t1~t2溫度從初始溫度T0升高到T1,t2~t3以T1溫度保溫,t3~t4繼續(xù)升溫到T2,t4~t5以T2溫度保溫,t5~t6逐漸降溫到初始溫度T0,固化完成。熱量通過線包表面空氣對(duì)流的形式傳遞到線包表面,然后再以熱傳導(dǎo)的形式傳遞到線包內(nèi)部。在熱作用下粘結(jié)劑的主要成分硅橡膠發(fā)生交聯(lián)反應(yīng)同時(shí)釋放出熱量,導(dǎo)致其機(jī)械性能發(fā)生顯著變化(液態(tài)→固態(tài)),從而起到穩(wěn)定光纖,防止其松動(dòng)的作用。
圖1 線包固化工藝曲線圖
線包結(jié)構(gòu)符合細(xì)觀力學(xué)理論基本假設(shè):
1)光纖和粘結(jié)劑纖維均勻分布在整個(gè)線包之中;
2)光纖和粘結(jié)劑表面直接接觸,并且互相或者以化學(xué)的或者以物理的方法粘結(jié)在一起,直到光纖釋放前不脫開、不發(fā)生相對(duì)滑移;
3)線包中孔隙與氣泡體積的總和很小,可以忽略不計(jì)。
熱傳導(dǎo)模型采用傅里葉熱傳導(dǎo)定律[10]和能量平衡原理來建立:
(1)
式中:ρ,Cp為線包整體的等效密度和比熱容;kx,ky,kz分別為線包3個(gè)方向的熱傳導(dǎo)系數(shù);T為線包溫度;Q為內(nèi)部熱源,可表示為:
(2)
式中:Hr為單位體積硅橡膠固化完成后放出的總熱量;α和dα/dt分別為固化度和固化速率。
由于光纖纏繞方向都是沿著線軸周向,所以線包橫截面具有橫觀各向同性材料屬性[11],根據(jù)復(fù)合材料混合率公式,ρ,Cp可表示為:
ρ=Vfρf+(1-Vf)ρm
(3)
(4)
式中:Vf為光纖體積分?jǐn)?shù);ρf,ρm分別為光纖和粘結(jié)劑的密度;Cf,Cm分別為光纖和粘結(jié)劑的比熱容。ρm和Cm會(huì)根據(jù)粘結(jié)劑狀態(tài)的改變發(fā)生變化。
沿光纖纏繞方向的熱導(dǎo)率ky為:
ky=Vfkf+(1-Vf)km
(5)
垂直光纖纏繞方向的截面兩個(gè)方向的熱導(dǎo)率kx和kz為:
(6)
式中B=2(km/kf-1);km為粘結(jié)劑熱傳導(dǎo)系數(shù),隨著粘結(jié)劑狀態(tài)的改變而改變;kf為光纖熱傳導(dǎo)系數(shù)。
液體硅橡膠的自催化固化動(dòng)力學(xué)模型[12]是將固化速率和溫度聯(lián)系起來的唯象模型,其固化反應(yīng)速率為:
(7)
式中:α為固化度,用于表征硅橡膠的固化程度,未開始固化時(shí)α=0,固化完成后α=1;n為反應(yīng)級(jí)數(shù),模型中取0.822;k(T)為液體硅橡膠改進(jìn)自催化反應(yīng)動(dòng)力學(xué)模型的反應(yīng)速率常數(shù),其表達(dá)式為:
k(T)=k0e-ΔEa/RT
(8)
式中:k0為反應(yīng)動(dòng)力學(xué)模型的頻率因子;ΔEa為反應(yīng)動(dòng)力學(xué)模型的活化能;R為氣體常數(shù)。
考慮到溫度和固化度之間的耦合關(guān)系,分別建立熱傳導(dǎo)和固化動(dòng)力學(xué)計(jì)算模塊,利用仿真軟件ANSYS的參數(shù)化編程語言APDL實(shí)現(xiàn)兩個(gè)模塊并進(jìn)行耦合求解。
圖2所示為熱傳導(dǎo)模型和固化動(dòng)力學(xué)模型耦合模擬過程流程圖,在模型中輸入初始溫度、固化溫度、光纖和粘結(jié)劑初始材料屬性、光纖體積分?jǐn)?shù)、初始固化度、固化速率、輸出溫度場(chǎng)和固化度場(chǎng)。
圖2 熱傳導(dǎo)模型和固化動(dòng)力學(xué)模型耦合模擬過程流程圖
熱傳導(dǎo)模型為基礎(chǔ)模型,它為固化動(dòng)力學(xué)模型提供溫度參數(shù),將溫度和時(shí)間輸入固化動(dòng)力學(xué)模型便可輸出固化度,熱傳導(dǎo)模型再利用溫度和固化度參數(shù)計(jì)算出時(shí)變材料屬性,為新的固化周期初始化。
結(jié)合細(xì)觀力學(xué)理論對(duì)材料模型進(jìn)行等效,線包幾何模型結(jié)構(gòu)將不再區(qū)分光纖和粘結(jié)劑,幾何模型尺寸由某產(chǎn)品線包實(shí)體1∶1簡(jiǎn)化而成,線軸內(nèi)徑為113.89 mm(忽略錐度),線軸壁厚為5 mm,線軸長(zhǎng)度為216 mm,線包初始匝數(shù)522,兩邊退匝數(shù)為2,纏繞層數(shù)為70。
線包結(jié)構(gòu)材料屬性如表1所示,在硅橡膠固化過程中其熱物理特性將隨溫度和固化度發(fā)生變化,而光纖和線軸材料屬性隨溫度變化可以忽略不計(jì)。固化動(dòng)力學(xué)參數(shù)如表2所示。
表1 線包結(jié)構(gòu)材料參數(shù)
表2 固化動(dòng)力學(xué)參數(shù)
選取三維8節(jié)點(diǎn)六面體單元SOLID70單元進(jìn)行熱分析,采用掃略網(wǎng)格劃分的方法對(duì)線包進(jìn)行網(wǎng)格劃分。單元坐標(biāo)系方向與單元IJ邊平行,通過掃略劃分的網(wǎng)格得到規(guī)則的六面體單元可以很方便定義正交各向異性單元材料屬性。在保證計(jì)算精度的前提下,盡量選取相對(duì)較大的單元尺寸和較少的單元數(shù)量,可以很好的節(jié)省計(jì)算空間和時(shí)間,經(jīng)過算例分析,綜合選取的單元數(shù)為17 856個(gè)。劃分好網(wǎng)格的三維有限元模型如圖3所示。
圖3 線包三維有限元模型
線包各個(gè)位置的初始溫度和表面空氣初始溫度都設(shè)置為25 ℃,初始固化度為0。邊界條件為熱力學(xué)第三類邊界條件[9]:
(9)
式中:hf為線包與空氣接觸部分表面熱對(duì)流系數(shù),其值取20 W/(m2·K);Tf為線包表面空氣溫度即為固化溫度。試驗(yàn)中的另一個(gè)工藝參數(shù)升溫速率為開始加熱后表面空氣溫度Tf由初始溫度25 ℃上升到指定固化溫度過程中溫度升高的速率。
由于固化過程既與溫度有關(guān),也與時(shí)間有關(guān),因此需要瞬態(tài)熱求解器對(duì)模型進(jìn)行求解。同時(shí)由于線包線軸存在接觸問題,因此將采用非線性求解器求解。
求解過程:第一個(gè)載荷步根據(jù)初始條件計(jì)算出線包溫度分布,再根據(jù)溫度分布及固化時(shí)間計(jì)算出線包固化度分布以及固化釋放的熱量,并將釋放熱量作為載荷施加到下一個(gè)載荷步中。已知溫度分布和固化度分布后更新下一個(gè)載荷步線包每個(gè)單元的材料屬性,然后開始下一個(gè)載荷步,以此類推,直到完成固化工藝,最終得到線包溫度場(chǎng)和固化度場(chǎng)分布。
由于線包結(jié)構(gòu)是一個(gè)軸對(duì)稱結(jié)構(gòu),其溫度和固化度在圓周方向分布相同,因此選取任一橫截面上如圖4所示A~E五點(diǎn)表征線包不同部位固化度和溫度分布,后面所述的不同位置指橫截面的不同位置。其中A,B,C為線包截面中心同一厚度方向上的3個(gè)點(diǎn),分別代表線包厚度方向中點(diǎn)、線包表面以及線包與線軸接觸點(diǎn);D點(diǎn)是線包頂層最邊緣點(diǎn);E點(diǎn)是D點(diǎn)對(duì)應(yīng)的線包中間層點(diǎn)。
圖4 結(jié)果數(shù)據(jù)讀取點(diǎn)
圖5為線包在如圖1所示工藝條件下A,B,C三點(diǎn)的溫度、固化度在固化過程中的隨時(shí)間變化曲線。從圖中可以看出在固化初始階段,固化速率緩慢增長(zhǎng),隨著線包溫度升高,固化速率逐漸上升到最大值,隨著固化程度的提高固化速率逐漸放緩直至固化完成。由于B點(diǎn)處于線包最頂層,熱量最先傳遞到線包表面,因此3點(diǎn)中B點(diǎn)固化溫度最先升高,最早完成固化。C點(diǎn)位于線包初始層,與鋁制線筒接觸,而鋁的導(dǎo)熱能力強(qiáng),因此從圖中可以看出C點(diǎn)與B點(diǎn)的溫度和固化度接近,但沒有B點(diǎn)高。而A點(diǎn)為線包中心點(diǎn),由于線包中粘結(jié)劑體積含量較低導(dǎo)致反應(yīng)放熱量小,熱傳導(dǎo)占主導(dǎo)地位,且光纖和粘結(jié)劑導(dǎo)熱系數(shù)都比較低,線包厚度厚,所以在固化工藝升溫階段線包中心A點(diǎn)溫度、固化度低于B點(diǎn)和C點(diǎn)。
圖5 線包A,B,C三點(diǎn)固化度和溫度曲線圖
另外從圖中可以看出在相同時(shí)間點(diǎn)B點(diǎn)與A點(diǎn)固化度的最大差值出現(xiàn)在固化度曲線斜率最大的階段,即固化速率最高階段。
圖6為A,B,D,E四點(diǎn)在固化過程中溫度、固化度隨時(shí)間變化曲線。從圖中可以看出在固化過程中D和B點(diǎn)比E和A點(diǎn)的溫度和固化度相對(duì)要高一些,這是因?yàn)镈和B點(diǎn)所在的線包表面離熱源更近。另外D點(diǎn)位于外圓柱面與退匝錐面相交處,在同樣的時(shí)間內(nèi)D點(diǎn)比B點(diǎn)能獲得更多的熱量,所以D點(diǎn)比B點(diǎn)固化度高。同理E點(diǎn)要比A點(diǎn)獲得更多熱量,因此固化度高。從上面的分析可以看出隨著熱量由線包表面?zhèn)鬟f到內(nèi)部,固化過程也是由外及內(nèi)。
圖6 線包A,B,D,E四點(diǎn)固化度和溫度曲線圖
固化均勻性是通過溫度均勻性和固化度均勻性來衡量的。溫度均勻性是通過固化過程中在同一時(shí)間點(diǎn)線包中最高溫度點(diǎn)與最低溫度點(diǎn)差值的大小來表征的,后面簡(jiǎn)稱為最大溫度差值。同理固化度均勻性通過固化過程中在同一時(shí)間點(diǎn)線包中固化度最高點(diǎn)與固化度最低點(diǎn)差值的大小來表征,簡(jiǎn)稱最大固化度差值。
在升溫速率8 ℃/min時(shí)研究不同固化溫度(固化溫度分別取90 ℃,110 ℃,130 ℃,150 ℃,175 ℃)對(duì)線包固化均勻性的影響。圖7和圖8繪制了不同固化溫度時(shí)線包中心點(diǎn)A處固化度曲線和固化速率曲線,從圖7和圖8可以看出固化溫度的升高會(huì)加速線包固化,縮短固化時(shí)間。
圖7 不同固化溫度對(duì)應(yīng)線包中心點(diǎn)固化度曲線
圖8 不同固化溫度對(duì)應(yīng)線包中心點(diǎn)固化速率曲線
圖9是不同固化溫度下完成固化所用時(shí)間,從圖中可以看出隨著固化溫度的升高,固化所用時(shí)間縮短,固化效率升高。但當(dāng)固化溫度升高到一定程度之后,升高固化溫度所引起的固化時(shí)間的縮短不再明顯。
圖9 不同固化溫度下完成固化所用時(shí)間
圖10和圖11為固化過程中不同固化溫度下線包最大溫度差值和最大固化度差值曲線,從圖中可以看出隨著固化溫度升高,線包的最大溫度差值和最大固化度差值會(huì)隨之變大,固化更加不均勻。由于光纖纏繞層數(shù)多(一般為50~80),線包厚度厚,粘結(jié)劑和光纖的導(dǎo)熱性差,當(dāng)提高固化溫度后線包最大溫度差值變大,從前面固化溫度對(duì)固化速率的影響結(jié)果可知升高固化溫度會(huì)提高固化速率,這意味著線包中溫度高的位置固化速率快,溫度較低的位置固化速率相對(duì)緩慢,從而使得線包最大固化度差值變大。
圖10 固化溫度對(duì)線包最大溫度差值的影響
在固化溫度為165 ℃時(shí)研究升溫速率(升溫速率分別取4 ℃/min,6 ℃/min,8 ℃/min,10 ℃/min)對(duì)固化均勻性的影響。圖12和圖13描述了不同升溫速率下中心點(diǎn)A處固化度和固化速率的變化特性。從圖中可以看出同一時(shí)刻線包中心點(diǎn)的固化度和固化速率隨著升溫速率的提高而增大。這是因?yàn)樯郎厮俾实奶岣邥?huì)使得線包溫度與加熱介質(zhì)溫度的差值變大,兩者之間的熱交換變快,線包溫度升高,固化速率加快,但是升溫速率對(duì)固化速率的影響沒有固化溫度對(duì)其影響顯著。
圖12 不同升溫速率對(duì)應(yīng)線包中心點(diǎn)固化度曲線
圖13 不同升溫速率對(duì)應(yīng)線包中心點(diǎn)固化速率曲線
圖14和圖15為固化過程中不同升溫速率對(duì)線包最大溫度差值和最大固化度差值的影響。從圖中可以看出升溫速率越高最大溫度差值越大,但是升溫速率對(duì)最大固化度差值影響不顯著。
圖14 升溫速率對(duì)線包最大溫度差值的影響
圖15 升溫速率對(duì)線包最大固化度差值的影響
當(dāng)固化速率為4 ℃/min時(shí)固化溫度由初始溫度25 ℃上升到165 ℃需要2 100 s,而從固化度曲線可以看出在這段時(shí)間內(nèi)固化速率處于緩慢上升階段,固化速率不高,升溫速率對(duì)固化速率的影響較小。而從前面的結(jié)果可以看出固化度最大差值出現(xiàn)在固化速率最快的時(shí)間段,而這個(gè)時(shí)間段固化溫度已經(jīng)上升到165 ℃,從式(7)可以看出影響固化速率的因素是固化溫度,因此升溫速率對(duì)固化度最大差值的影響很小。
基于復(fù)合材料細(xì)觀力學(xué)理論對(duì)線包固化過程進(jìn)行了有限元仿真,分析了線包不同位置的溫度和固化度分布,研究了固化溫度和升溫速率對(duì)固化均勻性的影響,主要結(jié)論為:
1)由于線包厚度較大,且熱傳導(dǎo)性較差,靠近表面熱源的部位最先完成固化,線包中心位置最后完成固化,固化過程由外及內(nèi)。內(nèi)外固化度最大差值出現(xiàn)在固化速率最快階段。
2)升高固化溫度可以明顯提高固化速率,縮短固化時(shí)間,但也會(huì)導(dǎo)致最大溫度差值和最大固化度差值變大,從而導(dǎo)致線包均勻性變差,因此應(yīng)合理設(shè)置固化溫度。
3)升溫速率對(duì)固化均勻性的影響不顯著,因此升溫速率的選擇應(yīng)主要以考慮生產(chǎn)效率為主。