熱壓罐成型工藝是目前廣泛應用于先進復合材料結構、蜂窩夾層結構及復合材料膠接結構的主要成型方法。在成型過程中，復合材料制件是在高溫高壓下與工裝一起放進熱壓罐中固化成型，制件固化成型后幾乎不再做任何加工，外表尺寸應滿足裝配協調要求，不允許強迫裝配。然而，工裝在整個成型過程中因承受高溫熱載荷、自身重力及成型輔助件壓力的共同作用而發生變形，工裝的變形直接影響到制件固化變形而最終影響到復合材料制件的實際形狀和尺寸精度。
針對目前較為常用的熱壓罐成型工裝，國內外對熱壓罐成型工藝參數和工裝結構因素對工裝溫度場均勻性的影響進行了研究[1-3]，并研究了固化工裝變形對復合材料內部殘余應力的影響[4-6]。岳廣全等[7] 針對當前應用較多的框架式工裝在工藝過程中的變形進行了研究。
本文以柱支撐形式建立大型復合材料固化工裝模型，在保證支撐形式滿足工裝靜力剛度的前提下，研究其在熱壓罐成型過程中熱- 結構耦合作用下的結構變形。 123,123

柱支撐形式固化工裝

對大型復合材料固化工裝，工裝模具的型面需要根據固化成型后的復合材料曲面進行調整，復合材料成型的多樣性勢必增加工裝制作的經濟成本和制作周期。針對此問題，對以柱支撐為支撐形式的工裝，支撐部分與工裝模具的接觸面積小，工裝模具型面改變后，支撐柱的生產加工簡單易行；此外，工裝內部結構簡單，利于熱空氣的流通，可以提高工裝和熱空氣之間的傳熱效率。
為方便后續研究，對柱支撐形式的工裝，以結構相對簡單的平板型面工裝為例（見圖1），分析工裝在復合材料固化成型過程的變形情況。圖1中工裝的型面厚度為20mm，模型尺寸為15000mm×5000mm×1000mm，支撐柱截面設為正方形，尺寸為100mm×100mm。模具材料為殷鋼，密度ρ=7900kg/m3，彈性模量E=214GPa, 泊松比ν=0.3，線膨脹系數α=1.7×10-6/K。 copyright 123456



對工裝進行靜力分析，可通過調節支撐柱數量來控制工裝型面變形量大小。為了比較靜力作用下的變形量和熱- 結構耦合作用下的變形量，使用70 個支撐柱。對圖1中工裝進行靜力分析（見圖2），工裝總變形量為0~0.00635mm，垂直型面方向變形量為0~0.00627mm ；靜力作用下，70個支撐柱的工裝結構變形量微小，說明以支撐柱為支撐形式的工裝完全可以滿足靜力作用下的剛度要求。



成型工裝熱- 結構耦合分析

由經驗和工藝試驗可知工裝的變形均在毫米量級，工裝的尺寸在幾米至十幾米，工裝結構變形不會對熱壓罐內的氣體流動產生明顯影響。由于流體計算對網格要求較高，文中使用專業流體仿真軟件FLUENT進行工裝及熱壓罐內整個流場的網格劃分，如圖3所示。對工裝建立的有限元模型既可以適用于溫度場的計算，也可以用于工裝的結構變形計算。熱壓罐工藝規范簡化如圖4所示，溫度場計算所用相關材料參數如表1所示。

對復合材料固化工裝而言，工裝型面的變形將直接影響與其接觸的復合材料制件的成型質量，垂直型面方向的翹曲變形將直接影響到大尺寸復合材料成型的曲率精度，所以型面變形場的表征方法應該與型面的變形直接相關。為對成型工裝的型面變形量進行表征，所用支撐柱如圖5所示。



圖1所示工裝中，其型面存在眾多可用于表征型面變形量的位置，為方便解釋問題，以支撐柱坐標所對應的型面位置的變形量來表征型面變形。鑒于模型所用支撐柱數量較多且考慮到模型結構的對稱性，選取代表性的支撐柱坐標：U1V1、U7V1、U14V1、U1V3、U7V3、U14V3、U1V5、U7V5、U14V5、U3V1、U3V3 和U3V5，共12 個坐標位置。對流體仿真在整個熱歷程中計算得到的工裝溫度場，每隔100s存儲一個溫度場數據文件；將包含節點溫度值的有限元模型導入有限元分析軟件ANSYS，同時施加重力載荷，計算得到對應時刻熱-結構耦合作用下工裝的結構變形量。 copyright 123456

計算結果分析

復合材料在熱壓罐中固化分為3個熱歷程：升溫、恒溫和降溫，此過程中工裝溫度場分布的不均勻性導致各部分熱應力分布不均勻，工裝在熱應力和自重載荷耦合作用下產生變形。熱歷程中工裝變形量最大值分布如圖6所示，從圖中可得出：
（1）整個熱歷程中，工裝總變形量最大值達4.662mm。
（2）工裝沿自身長度方向變形量最大值為4.596mm，與總變形量最大值較為接近。
（3）工裝沿自身寬度方向和高度方向的變形量最大值分別為0.805mm 和0.339mm，其變形量相對于工裝長度方向變形量較小。



熱應力不均勻導致的變形量最大值為靜力載荷導致工裝變形量的53~734 倍，這說明在復合材料固化成型過程中，熱應力載荷對工裝結構變形起著決定性作用。
圖6反映出工裝總變形量最大值隨時間的變化，不能得出熱應力不均勻導致的型面翹曲對型面的影響。對工裝型面翹曲的表征從沿工裝長度方向（圖5 中U 方向）和沿工裝寬度方向（圖5 中V 方向）這兩個方面來闡述：沿工裝長度方向的翹曲，采用V1 上的U1V1、U7V1、U14V1，V3上的U1V3、U7V3、U14V3 和V5上的U1V5、U7V5、U14V5 位置作為監測點；沿工裝寬度方向的翹曲，采用U3上的U3V1、U3V3、U3V5，U7上的U7V1、U7V3、U7V5 和U14 上的U14V1、U14V3、U14V5 位置作為監測點。由圖7得知：
（1）沿V3方向，其中間位置、迎風端及背風端位置由于均受到加熱而使Z向變形量在熱歷程升溫階段逐漸增大；進入恒溫階段后，中間位置處Z向變形量基本趨于平穩，迎風端和背風端由于與環境熱空氣的溫差不一致導致升溫幅度不同，各自的Z向變形量不同；降溫階段，其變形量趨勢基本與升溫階段相反。
（2）沿V1和V5方向，迎風端位置Z 向變形量在熱歷程中先是低于后又高于中間位置和背風端位置的Z向變形量。 123456
（3）迎風端位置相對于中間位置的Z向變形量先是低于后又高于背風端位置相對于中間位置的Z向變形量。這是因為迎風端雖然首先與環境空氣相接處，但工裝內部的支撐柱在空氣流經工裝內部時對其形成擾流作用，使工裝內局部空氣流動速度明顯大于工裝外部，中間位置的支撐柱與空氣換熱速率增大。
（4）V3方向，內部支撐柱的擾流作用使背風端位置相對于中間位置的Z向變形量先急劇增大后趨于平穩，然后迅速下降；V1和V5方向，其處于工裝外圍與工裝外部流動空氣直接接觸，背風端位置相對于中間位置的Z向變形量在整個熱歷程中變化較平緩。



另外，研究表明：
（1）沿工裝寬度方向，翹曲變形基本處于對稱狀態。
（2）在U7方向上，兩端位置相對中間位置的Z向變形量在升溫階段急劇增大，在恒溫階段達到最大值并緩慢下降，進入降溫階段后急劇減小。 123,123
（3）在U3和U14方向上，監測位置處的Z向變形量曲線在熱歷程中基本保持一致；因為工裝寬度方向尺寸相對較小，在熱歷程不同階段，空氣流態在寬度方向上對溫度場差值影響較小。
（4）在U3和U14方向上，兩端位置相對中間位置的Z向變形量在升溫階段急劇增大，進入恒溫階段后迅速下降；在恒溫階段后期出現緩慢變化，進入降溫階段又迅速下降。
在復合材料的熱壓罐固化過程中，對柱支撐形式的固化工裝型面，溫度不均勻性導致的型面翹曲變形較為復雜。在工裝結構內部，支撐柱存在位置所對應型面的局部變形量大于沒有支撐柱存在的位置，這是因為支撐柱對空氣的擾流作用加快了支撐柱與周圍空氣之間的熱傳遞。圖8為3200s時工裝溫度場及Z向結構變形量分布云圖。



結束語 123456

本文以柱支撐形式建立大型復合材料固化工裝模型，在保證工裝結構靜力剛度的前提下，計算得到工裝在熱壓罐固化工藝過程中的溫度場分布，然后把重力載荷，施加到帶有溫度值的有限元模型上進行工裝的熱- 結構耦合分析。得到結論如下：
（1）成型工裝在熱- 結構耦合作用下，相對于重力載荷熱應力載荷對工裝結構變形起決定性作用。
（2）對工裝型面在垂直型面方向的變形量，有支撐柱存在的位置大于無支撐柱存在的位置。（責編　良辰）