周晨，王志瑾，候天驕

?

V型皺褶芯材夾層結(jié)構(gòu)強迫對流傳熱與熱應(yīng)力分析

周晨1，王志瑾1，候天驕2

（1. 南京航空航天大學(xué)飛行器先進(jìn)設(shè)計技術(shù)國防重點學(xué)科實驗室；2. 南京航空航天大學(xué)航天學(xué)院：南京210016）

V型皺褶芯材夾層結(jié)構(gòu)與強迫對流冷卻相結(jié)合，可應(yīng)用于高超聲速飛行器或大功率電子器件承載和散熱的多功能設(shè)計。文章采用數(shù)值仿真方法對強迫對流條件下V型皺褶芯材夾層結(jié)構(gòu)的換熱及其產(chǎn)生的熱應(yīng)力進(jìn)行研究，獲得了相應(yīng)的速度場、溫度場及結(jié)構(gòu)熱應(yīng)力分布。結(jié)果表明：強迫對流條件下，V型皺褶芯材夾層結(jié)構(gòu)的換熱性能明顯提高，同時導(dǎo)致較大的壓力損失。流體速度場隨著幾何構(gòu)型的改變而發(fā)生周期性的變化，在結(jié)構(gòu)偏折處波峰一側(cè)流速達(dá)到最大，對流換熱得到增強；沿著冷卻液流動方向，結(jié)構(gòu)與冷卻液的溫度逐漸升高，并伴隨著周期性的波動；皺褶芯材夾層結(jié)構(gòu)的整體換熱能力隨著入口流速的增大而增強。熱流輸入側(cè)的面板熱應(yīng)力和變形均較大，結(jié)構(gòu)皺褶部位存在應(yīng)力集中。

皺褶芯材；強迫對流；夾層結(jié)構(gòu)；換熱性能；熱應(yīng)力

0 引言

皺褶芯材是將平板基材按照一定規(guī)律的線系網(wǎng)格進(jìn)行局部折疊而得到的一種具有周期性胞元構(gòu)型的三維結(jié)構(gòu)[1]。其設(shè)計思想源于手工折紙，從夾層結(jié)構(gòu)出現(xiàn)之初就被提出用于芯材設(shè)計，但由于欠缺高效的生產(chǎn)制備工藝，對皺褶芯材的研究在相當(dāng)長的一段時間內(nèi)都處于停滯狀態(tài)[2]。近年來，新型制備工藝[3-5]的誕生令皺褶芯材重新受到國內(nèi)外研究人員的關(guān)注。王志瑾等對皺褶芯材幾何設(shè)計方法[1]進(jìn)行了研究，設(shè)計了平面皺褶芯材和單曲面皺褶芯材；針對皺褶芯材夾層板在隔聲[6]、防熱[7-8]及吸能[9]等方面進(jìn)行了數(shù)值計算和實驗研究。張延昌等[10]利用非線性有限元分析了皺褶芯材幾何參數(shù)對結(jié)構(gòu)耐撞性能的影響。Heimbs[11]和Fischer[12]等對皺褶芯材夾層板進(jìn)行了一系列壓縮、剪切和沖擊的實驗研究和數(shù)值仿真。皺褶芯材成型簡單，可設(shè)計性強，在隔熱、隔聲等領(lǐng)域有著廣闊的應(yīng)用前景；同時，其開環(huán)的連續(xù)空腔構(gòu)型使其具有優(yōu)良的散熱性能和流通性能[12-13]，更有利于實現(xiàn)結(jié)構(gòu)的多功能設(shè)計。

以多孔材料作為芯材的夾層結(jié)構(gòu)通常具有至少一個方向的流動通道，采用強迫對流形式通入冷卻液，可用于高超聲速飛行器、發(fā)動機燃燒室或高功率電子器件等的主動冷卻。高亮、孫雨果等[14-17]對點陣夾層結(jié)構(gòu)的流動特性和換熱特性進(jìn)行了理論分析和數(shù)值模擬。Rakow等[18]分別從實驗、數(shù)值分析和微觀角度對強迫對流主動冷卻下的金屬泡沫夾層板進(jìn)行了傳熱和結(jié)構(gòu)熱響應(yīng)研究。Feng等[19-20]采用數(shù)值仿真和實驗的方法研究了加裝肋片的金屬泡沫在強迫對流下的換熱性能。王儲等[21]建立了層板內(nèi)部流動的數(shù)值模型，研究了層板幾何參數(shù)及雷諾數(shù)對流動特性的影響。

本文以V型皺褶芯材夾層結(jié)構(gòu)作為研究對象，采用數(shù)值仿真的方法研究了在以水作為冷卻液的強迫對流條件下，該結(jié)構(gòu)的傳熱以及相應(yīng)的熱-結(jié)構(gòu)耦合問題。通過計算先后得到了V型皺褶芯材夾層結(jié)構(gòu)的溫度場及相應(yīng)的結(jié)構(gòu)應(yīng)力場和變形場；根據(jù)結(jié)果分析了其流道中流體的流動特性，進(jìn)一步探討了對結(jié)構(gòu)換熱性能的影響機制；并對不同冷卻液流速下該結(jié)構(gòu)的總體換熱性能進(jìn)行了分析和評價。

1 模型描述

圖1為典型的強迫對流冷卻下的V型皺褶芯材夾層結(jié)構(gòu)示意圖。熱量通過上面板傳入夾層結(jié)構(gòu)內(nèi)部，沿方向通入冷卻液利用強迫對流進(jìn)行散熱，從而提高結(jié)構(gòu)承受熱載荷的能力。如圖1所示，V型皺褶芯材夾層結(jié)構(gòu)的流道有2種類型：綠線所經(jīng)過的正三角型流道（由芯材與下面板構(gòu)成）以及藍(lán)線所經(jīng)過的倒三角型流道（由芯材與上面板構(gòu)成）。

圖1 受熱V型皺褶芯材夾層結(jié)構(gòu)強迫對流冷卻示意

V型皺褶芯材胞元可由4個獨立的參數(shù)表示[1]，如圖2(a)所示。在對有厚度的平板基材進(jìn)行皺褶成型時，可能造成折痕位置出現(xiàn)局部缺陷或材料積聚[22]。假設(shè)有厚度的芯材經(jīng)過皺褶后如圖2(b)所示，本文所研究的V型皺褶芯材夾層結(jié)構(gòu)的幾何尺寸參數(shù)見表1，其中c和f分別為芯材和面板的厚度。

(a) V型皺褶芯材胞元示意

(b) 帶厚度的V型皺褶芯材胞元示意

圖2 V型皺褶芯材夾層結(jié)構(gòu)胞元及幾何尺寸示意

Fig. 2 Geometric parameters of a V-pattern folded core unit cell

表1 皺褶芯材夾層結(jié)構(gòu)幾何尺寸

本研究采用水作為冷卻液，面板和芯材采用鋁合金Al 2024-T851。根據(jù)初步分析可知，在本研究所涉及的溫度變化范圍內(nèi)，鋁合金材料屬性變化不大，可假設(shè)夾層結(jié)構(gòu)固體部分的熱物理屬性及力學(xué)性能不隨溫度發(fā)生變化，其相關(guān)性能參數(shù)[23]如表2所示。

表2 Al 2024-T851的物理性能參數(shù)

2 強迫對流共軛傳熱分析

2.1 耦合傳熱計算

強迫對流條件下，皺褶芯材夾層結(jié)構(gòu)的共軛傳熱過程包括：流體的流動與傳熱、固體結(jié)構(gòu)的熱傳導(dǎo)以及固體結(jié)構(gòu)與流體之間的對流換熱。本文采用Fluent軟件對計算區(qū)域進(jìn)行整場離散、整場求解。由于皺褶芯材與面板構(gòu)成周期性的流道，為了減少計算量，對圖1所示V型皺褶芯材夾層板沿向取1個周期（見圖3）作為分析對象。

圖3 簡化后的周期模型示意

在進(jìn)行數(shù)值模擬時，假設(shè)：1）流動和傳熱是穩(wěn)態(tài)的；2）流體為不可壓縮的牛頓型流體；3）忽略結(jié)構(gòu)與周圍環(huán)境的自然對流和輻射換熱；4）不考慮黏性耗散熱。

基于上述假設(shè)，流場的控制方程采用笛卡兒張量符號可表示為[24]：

連續(xù)方程：

式中為流體速度矢量。

動量方程：

式中：為冷卻液的密度；為壓力；τ為黏性應(yīng)力張量，在不可壓縮條件下可表示為

， (3)

其中為黏性系數(shù)。

能量方程：

式中：c、f分別為冷卻液的比熱容和熱傳導(dǎo)系數(shù)；f為流體溫度。

對于固體區(qū)域，穩(wěn)態(tài)導(dǎo)熱微分方程可直接由式(4)在不計流體本身流動的情況下推導(dǎo)得到，

式中：s為固體材料的熱傳導(dǎo)系數(shù)；s為固體的溫度。

狀態(tài)方程：

=。 (6)

在固液界面上，根據(jù)溫度和熱流的連續(xù)性條件可得：

。 (8)

在Fluent共軛傳熱計算過程中，流體域和固體域采用通用的控制方程，即雷諾平均N-S方程（RANS），區(qū)別在于固體域只需求解能量方程即可。固-液耦合界面的對流換熱系數(shù)由Fluent計算直接得到。選用RNG-雙方程湍流模型，在近壁區(qū)域采用增強壁面函數(shù)進(jìn)行處理。固體域采用結(jié)構(gòu)化網(wǎng)格劃分；流體域采用非結(jié)構(gòu)網(wǎng)格劃分，并在靠近壁面位置采用邊界層網(wǎng)格進(jìn)行局部加密。

2.2 熱載荷和邊界條件

在V型皺褶芯材夾層結(jié)構(gòu)的上面板施加恒定的熱流密度[25]=100000W/m2。沿流道方向（向）通入冷卻液（水）對結(jié)構(gòu)進(jìn)行對流冷卻，采用速度入口邊界條件，入口處冷卻液溫度為0= 300K，流速0=0.1m/s；采用壓力出口邊界條件，出口處設(shè)置參考壓力為0Pa。由于沿向只取了1個周期作為研究對象，所以向側(cè)面均設(shè)置為周期性邊界條件，固-液界面為無滑移、熱耦合邊界，其他壁面均設(shè)置為絕熱壁面。

2.3 網(wǎng)格獨立性驗證

CFD計算結(jié)果是否準(zhǔn)確很大程度上依賴于計算區(qū)域網(wǎng)格的劃分。但是若網(wǎng)格密度過大，會造成計算量急劇增大。因此，需要在保證計算精度的前提下，盡可能地減少網(wǎng)格數(shù)量。本文對網(wǎng)格劃分進(jìn)行了獨立性驗證，針對多種網(wǎng)格劃分方案分別進(jìn)行了計算，圖4給出了整體結(jié)構(gòu)的平均努塞爾數(shù)隨網(wǎng)格數(shù)量的變化曲線?？梢钥吹?，隨著網(wǎng)格數(shù)量的增大，努塞爾數(shù)逐漸變大，當(dāng)網(wǎng)格數(shù)量達(dá)到200000左右時逐漸趨于穩(wěn)定。最終，本文采用的網(wǎng)格單元總數(shù)為185300。

圖4 網(wǎng)格獨立性驗證

3 結(jié)構(gòu)熱應(yīng)力分析

利用共軛傳熱分析得到的結(jié)構(gòu)溫度場作為載荷，給定邊界條件，求解熱彈性方程，得到相應(yīng)的結(jié)構(gòu)應(yīng)力場和變形場，在不受外加機械載荷情況下的熱彈性基本方程如下[26]：

平衡方程：

幾何方程：

。 (10)

在溫度發(fā)生變化的情況下，彈性體的應(yīng)變由2部分組成：1）由于自由膨脹而產(chǎn)生的應(yīng)變分量，即=Δ，對應(yīng)的剪應(yīng)變分量為0；2）在熱膨脹時彈性體內(nèi)各部分之間由相互約束產(chǎn)生的應(yīng)變分量，它們與熱應(yīng)力之間滿足胡克定律。因此，變溫情況下的物理方程為

式(9)～式(11)中：σ為應(yīng)力張量；ε為應(yīng)變張量；為位移矢量；δ為克羅內(nèi)克符號；s為柔度系數(shù)；為熱膨脹系數(shù)；Δ為溫度變化。

在熱應(yīng)力計算時，將上述流場與傳熱分析時得到的皺褶結(jié)構(gòu)溫度場作為載荷，入口端約束u=u=u=0，出口端約束u=0，上、下面板的向側(cè)面則采用對稱約束u=0。

4 仿真結(jié)果及分析

4.1 V型皺褶芯材夾層結(jié)構(gòu)速度場

為了使仿真結(jié)果顯示更加直觀，圖5～圖9均沿向顯示了5個周期。冷卻液沿方向流入，由于面板和芯材的干擾和阻礙，其速度將不斷發(fā)生變化；同時，由于V型皺褶芯材夾層結(jié)構(gòu)在方向上也是周期排列，所以在穩(wěn)態(tài)條件下，速度場將呈現(xiàn)周期性變化規(guī)律。

V型皺褶芯材夾層結(jié)構(gòu)在幾何構(gòu)型上以胞元為單位，周期排列，在每個流道內(nèi)，峰、谷依次交替。圖5為流體區(qū)域緊靠上面板位置（=0.0109m）的截面速度分布云圖。整體來說，速度分布隨胞元發(fā)生周期性變化，進(jìn)口和出口效應(yīng)的主要影響區(qū)域均為約1個胞元的長度。圖6為中間胞元在同一位置截面的平面速度分布矢量圖，從圖中可以看出，冷卻液在峰前得到加速，隨后一方面在峰后產(chǎn)生回流，另一方面繼續(xù)加速直接沖擊下游另一側(cè)的峰前壁面。

圖5 截面z=0.0109m上的速度分布云圖

圖6 中間胞元在截面z=0.0109m上的速度分布矢量圖

圖7為沿方向上各截面的速度分布云圖，其中圖7(a)為沿流動方向前4個胞元的中央截面，圖7(b)為后3個胞元的出口截面。從圖中可以清晰地看出，除了第1個和最后1個胞元外，中間各胞元在相對應(yīng)的截面上速度分布基本一致。流體在波峰處流速最大；在波谷貼近固體壁面的位置，由于受到壁面黏性和下游芯材的阻礙，流速最低。

(a) 截面=0.015, 0.045, 0.075, 0.105m

(b) 截面=0.09, 0.12, 0.15m

圖7 沿向各截面速度分布云圖

Fig. 7 Velocity distribution at cross sections alongdirection

4.2 V型皺褶芯材夾層結(jié)構(gòu)溫度場

對于單相流體的強迫對流換熱，強化傳熱最根本的措施就是減薄或破壞邊界層，增強流體的徑向混合[27]。圖8為穩(wěn)態(tài)下V型皺褶芯材夾層板的結(jié)構(gòu)溫度與冷卻液溫度分布云圖。從圖中可以看出，熱量由上面板向下傳遞，當(dāng)冷卻液沿方向流過腔體時，與上、下面板及芯材進(jìn)行換熱，通過自身溫度升高的方式存儲能量，帶走一部分原本將沿夾層板厚度方向傳遞的熱量，從而有效降低了下面板的溫度。沿冷卻液流動方向，結(jié)構(gòu)和流體的溫度均不斷升高，且溫度載荷在夾層板厚度方向上的滲透不斷加深，但滲透加深的速率逐漸減緩，這是由于參與下游換熱的冷卻液在上游已經(jīng)得到了加熱。

(a) 結(jié)構(gòu)溫度分布

(b) 冷卻液溫度分布

圖8 結(jié)構(gòu)溫度及冷卻液溫度分布云圖

Fig. 8 Temperature contours in structure and coolant

圖9為沿方向各胞元截面上的結(jié)構(gòu)溫度與冷卻液溫度分布云圖。從圖中可以看出，由于上面板為熱源輸入側(cè)，且直接構(gòu)成倒三角流道的底邊，所以在同一截面中，倒三角流道中冷卻液的溫度明顯高于正三角流道中的冷卻液溫度。熱量自上而下傳遞的過程中，在倒三角流道底邊中部區(qū)域，通過與冷卻液之間的對流進(jìn)行換熱；而在倒三角流道的兩個底角區(qū)域，上面板與芯材相互接觸，熱量由芯材通過熱傳導(dǎo)的方式向下傳遞，同時面板、芯材均與冷卻液進(jìn)行對流換熱，因此各截面的最大溫度值均出現(xiàn)在倒三角流道的底邊中部。此外，在出口截面上還有相當(dāng)一部分冷卻液溫度較低，說明其利用率還可進(jìn)一步提高。

圖9 沿x向各截面溫度分布云圖

為了便于觀察和比較結(jié)構(gòu)溫度分布，分別沿上、下面板外壁面向邊界（參見圖3）提取節(jié)點溫度。圖10為2種結(jié)構(gòu)沿上、下面板2條路徑（Polyline 1和Polyline 2）的溫度變化曲線。其中橫坐標(biāo)為各節(jié)點對應(yīng)的方向坐標(biāo)。

(a) 沿上面板Polyline 1溫度分布

(b) 沿下面板Polyline 2溫度分布

圖10 面板溫度沿流動方向分布

Fig. 10 Temperature distribution along the flow direction

從圖10可看出，上面板在入口端溫度迅速上升，在最初半個胞元位置處出現(xiàn)峰值，而后溫度先下降再緩慢上升，其中上面板的溫度呈現(xiàn)出小幅度的周期性波動。冷卻液平行于方向勻速流入，在前半個胞元內(nèi)，以一種相對平穩(wěn)的狀態(tài)流動，流速較低（參見圖5），因此該處對流換熱較弱；隨后由于流道發(fā)生偏折，皺褶芯材對來流形成阻礙和拖曳作用，冷卻液徑向運動加?。▍⒁妶D7），使得對流換熱得到增強。圖10(a)所對應(yīng)的路徑為圖3中的Polyline 1，該路徑位于倒三角流道的中線在上面板外壁面的投影處，因此，冷卻液與上面板之間的換熱為影響沿該路徑溫度變化的主要因素。從圖10(a)可看出，除去入口端和出口端，局部溫度最小值出現(xiàn)的位置均緊隨皺褶偏折位置，而局部溫度最大值均出現(xiàn)在每半個胞元的中后方位置。在每半個胞元中，溫度呈現(xiàn)出先上升、后下降的趨勢。這是由于上游冷卻液越過波峰達(dá)到流速最大，該處對流換熱最強，所以對應(yīng)出現(xiàn)溫度局部最小值，而后在一側(cè)產(chǎn)生回流，另一側(cè)遇到壁面先減速再不斷加速，從而導(dǎo)致溫度先升后降。

圖11為在冷卻液流動方向不同位置沿向的溫度分布情況。其中各路徑所在位置參見圖3。由圖11可知，冷卻液溫度沿向有著顯著的變化，靠近上面板（熱流輸入側(cè)）的位置溫度很高，溫度梯度很大，而中心區(qū)域溫度較低，這是由傳熱過程決定的，同時也表明處在流道中央?yún)^(qū)域的冷卻液利用率較低。另外，可以發(fā)現(xiàn)在靠近下面板區(qū)域的冷卻液溫度要高于中心區(qū)域的溫度，這是由于固體結(jié)構(gòu)的導(dǎo)熱能力更強，一部分熱量通過固體芯材側(cè)壁向下傳導(dǎo)至下面板，再由下面板與冷卻液之間進(jìn)行換熱造成的。同時，從圖11也可看出，冷卻液溫度沿流動方向不斷升高。

圖11 流體溫度沿橫向分布

4.3 不同雷諾數(shù)對結(jié)構(gòu)換熱性能的影響

針對V型皺褶芯材夾層結(jié)構(gòu)，以芯材高度為特征長度定義以下量綱為1的參數(shù)：

。 (13)

式中：為冷卻液入口流速；和f分別為冷卻液的密度、黏度和熱傳導(dǎo)系數(shù)；為結(jié)構(gòu)的平均換熱系數(shù)。

設(shè)置一系列不同的入口速度值，并保持其他參數(shù)設(shè)置不變，分別進(jìn)行求解。對于整個夾層板結(jié)構(gòu)，不同流速下對應(yīng)的雷諾數(shù)與平均努塞爾數(shù)之間的關(guān)系如圖12所示，隨著雷諾數(shù)的增加，結(jié)構(gòu)平均努塞爾數(shù)也不斷增大。強迫對流換熱中，結(jié)構(gòu)的平均努塞爾數(shù)可表示為雷諾數(shù)和普朗特數(shù)的函數(shù)，這里假設(shè)流體的普朗特數(shù)為常數(shù)，則平均努塞爾數(shù)與雷諾數(shù)之間的關(guān)系可擬合為

4.4 V型皺褶芯材與波紋芯材夾層結(jié)構(gòu)熱效率對比

通過上述分析可知，V型皺褶芯材夾層結(jié)構(gòu)由于芯材的皺褶使換熱得到增強。然而，曲折流道必然會導(dǎo)致比平直流道更高的壓力損失，因此有必要對結(jié)構(gòu)的綜合熱效率進(jìn)行評估。為了便于比較，取與上述V型皺褶芯材夾層結(jié)構(gòu)具有相同相對密度的波紋芯材夾層板的1個周期（如圖13所示）進(jìn)行研究。

圖13 波紋芯材夾層板周期模型示意

針對本文研究對象，泵功率可表示為

≈Δ·20。 (15)

將平均對流換熱系數(shù)表示為

式中：Δ為進(jìn)口和出口間的壓力損失；為總熱流；為熱流輸入側(cè)面板的面積，=20′（為了減少表達(dá)所涉及的參數(shù)，對于本文所討論的對象，寬即2，長即5×2），根據(jù)式(13)、(15)和(16)，采用文獻(xiàn)[16]中定義的指標(biāo)來表征結(jié)構(gòu)的綜合熱效率，可得

。 (17)

其中，f和10/2分別體現(xiàn)了流體屬性和外結(jié)構(gòu)幾何尺寸的影響，令

具有熱流密度倒數(shù)的量綱，其綜合考慮了結(jié)構(gòu)的換熱性能及對驅(qū)動冷卻液流動的循環(huán)動力的要求。

2種結(jié)構(gòu)的對比結(jié)果如表 3所示。從中可以發(fā)現(xiàn)，與波紋芯材夾層板相比，在相同入口速度下，V型皺褶芯材具有較高的換熱性能，但同時需要更大的壓力差來維持冷卻液的流動；在相同壓力差情況下，波紋芯材夾層板的努塞爾數(shù)更高，但表征結(jié)構(gòu)綜合熱效率的值要比V型皺褶夾層板小。可見，在動力允許的條件下，V型皺褶芯材夾層結(jié)構(gòu)換熱性能更優(yōu)，若對其皺褶幾何參數(shù)進(jìn)行優(yōu)化，綜合熱效率有望進(jìn)一步得到提高。

表 3 V型皺褶芯材與波紋芯材夾層結(jié)構(gòu)熱效率對比

4.5 V型皺褶芯材夾層結(jié)構(gòu)熱應(yīng)力

V型皺褶芯材夾層板在上述溫度場下結(jié)構(gòu)應(yīng)力與位移云圖如圖14所示。由應(yīng)力云圖（圖14(a)）可知，由于上面板溫度較高，所以總體來說上面板應(yīng)力水平較高，且在結(jié)構(gòu)偏折位置均存在應(yīng)力集中現(xiàn)象，其中最嚴(yán)重的部位發(fā)生在端部約束處，但尚未超出材料應(yīng)力的許用值。對于高超聲速飛行器而言，尤為關(guān)注在熱載荷作用下的結(jié)構(gòu)變形，因為較大的法向位移可能會破壞氣動外形，從而導(dǎo)致更嚴(yán)重的氣動加熱。由位移云圖（圖14(b)）可知，由于出口端約束較弱，且沿方向結(jié)構(gòu)溫度不斷升高，所以沿方向變形不斷變大。在夾層板厚度方向上，靠近熱源輸入側(cè)的結(jié)構(gòu)溫度變化較大，因此上面板變形也較大，導(dǎo)致結(jié)構(gòu)發(fā)生沿厚度方向的翹曲。向最大位移發(fā)生在結(jié)構(gòu)中后部位置（見圖14(c)）。

(a) 應(yīng)力云圖

(b) 總位移云圖

(c) z向位移云圖

5 結(jié)論

本文采用數(shù)值計算方法，對強迫對流冷卻條件下V型皺褶芯材夾層結(jié)構(gòu)的傳熱性能及熱-結(jié)構(gòu)耦合問題進(jìn)行了研究，得到以下結(jié)論：

1）在強迫對流條件下，冷卻液與固體結(jié)構(gòu)的對流換熱降低了結(jié)構(gòu)的溫度，有效提高了皺褶結(jié)構(gòu)承受熱載荷的能力，有利于實現(xiàn)其承載和散熱的多功能設(shè)計。

2）結(jié)構(gòu)的皺褶構(gòu)型導(dǎo)致冷卻液速度場發(fā)生周期性的變化，在皺褶位置對流換熱得到增強；沿流動方向，結(jié)構(gòu)與冷卻液的溫度在逐漸升高的同時伴隨有周期性的波動；結(jié)構(gòu)的整體換熱能力隨著冷卻液流速的增大而增強；與波紋芯材夾層結(jié)構(gòu)相比，V型皺褶結(jié)構(gòu)具有更優(yōu)異的換熱性能，但同時會帶來較大的壓力損失，在設(shè)計過程中對皺褶結(jié)構(gòu)幾何參數(shù)進(jìn)行優(yōu)化，可進(jìn)一步提高其效率。

3）結(jié)構(gòu)中的上面板熱應(yīng)力較大，由于幾何形狀的突變，在結(jié)構(gòu)皺褶部位出現(xiàn)了熱應(yīng)力集中，最大應(yīng)力發(fā)生在端部約束位置，最大法向位移發(fā)生在中后部位置。

作為潛在的承載/散熱多功能結(jié)構(gòu)，后續(xù)可綜合考慮驅(qū)動冷卻液所需的泵功率，結(jié)構(gòu)的溫度、應(yīng)力、變形及屈曲等因素，對V型皺褶芯材夾層結(jié)構(gòu)進(jìn)行優(yōu)化設(shè)計。

[1] 王志瑾, KHALIULIN V, SKRIPKIN E. 皺褶結(jié)構(gòu)芯格構(gòu)型的幾何設(shè)計方法[J]. 南京航空航天大學(xué)學(xué)報, 2002, 34(1): 6-11

WANG Z J, KHALIULIN V, SKRIPKIN E. Geometry design method of folded structure[J]. Journal of Nanjing University of Aeronautics and Astronautics, 2002, 34(1): 6-11

[2] LEBéE A, SAB K. Homogenization of thick periodic plates: application of the bending-gradient plate theory to a folded core sandwich panel[J]. International Journal of Solids and Structures, 2012, 49(19/20): 2778-2792

[3] ELSAYED E A, BASILY B B. Applications of folding flat sheets of materials into 3-d intricate engineering designs: Working paper 07-015[R]. Rutgers University, 2009

[4] 蔡克乾, 靳明山, 姜開宇. 新型折疊夾芯結(jié)構(gòu)材料研究進(jìn)展[J]. 材料導(dǎo)報, 2015, 29(4): 129-133

CAI K Q, JIN M S, JIANG K Y. Review of progress in materials for sandwich structure with folded core[J]. Materials Review, 2015, 29(4): 129-133

[5] ZHAO D, LI Y, WANG M, et al. Fabrication of polymer origami-based V-type folded core[C]//International Design Engineering Technical Conferences and Computers and Information in Engineering Conference: Volume 5B: 39thMechanisms and Robotics Conference. Boston, Massachusetts, USA, 2015

[6] 王志瑾, 徐慶華. V-型皺褶夾芯板與隔聲性能實驗[J]. 振動工程學(xué)報, 2006, 19(1): 65-69

WANG Z J, XU Q H. Experimental research on soundproof characteristic for the sandwich plates with folded core[J]. Journal of Vibration Engineering, 2006, 19(1): 65-69

[7] 王志瑾, 張輝, 徐慶華. 鋁合金皺褶芯材夾層板當(dāng)量導(dǎo)熱系數(shù)[J]. 南京航空航天大學(xué)學(xué)報, 2008, 40(4): 507-512

WANG Z J, ZHANG H, XU Q H. Equivalent thermal conductivity coefficient for sandwich plates with folded core[J]. Journal of Nanjing University of Aeronautics and Astronautics, 2008, 40(4): 507-512

[8] 周晨, 王志瑾, 支驕楊. 主動冷卻皺褶芯材夾層板的熱力分析[J]. 固體火箭技術(shù), 2014, 37(4): 545-550

ZHOU C, WANG Z J, ZHI J Y. Thermal-mechanical analysis of actively cooled folded core sandwich panels[J]. Journal of Solid Rocket Technology, 2014, 37(4): 545-550

[9] 周華志, 王志瑾. M-型皺褶芯材夾層板吸能性能

研究[J]. 航空學(xué)報, 2016, 37(2): 579-587

ZHOU H Z, WANG Z J. Analysis of energy absorption capability of M-type folded core sandwich structure[J]. Acta Aeronautica ET Astronautica Sinica, 2016, 37(2): 579-587

[10] 張延昌, 王自力, 張世聯(lián). 折疊式夾芯層結(jié)構(gòu)耐撞性能研究[J]. 船舶力學(xué), 2010, 14(1/2): 114-120

ZHANG Y C, WANG Z L, ZHANG S L. Simulation analysis of folded core structure under dynamic load[J]. Journal of Ship Mechanics, 2010, 14(1/2): 114-120

[11] HEIMBS S. Foldcore sandwich structures and their impact behaviour: an overview[Z]. Dordrecht: Springer Netherlands, 2013: 491-544

[12] FISCHER S. Aluminum foldcores for sandwich structure application: mechanical properties and FE-simulation[J]. Thin-Walled Structures, 2015, 90: 31-41

[13] GRZESCHIK M. Performance of foldcores mechanical properties and testing[C]// International Design Engineering

Technical Conferences and Computers and Information in Engineering Conference: Volume 6B: 37thMechanisms and Robotics Conference. Portland, Oregon, USA, 2013

[14] 孫雨果, 高亮. 復(fù)合材料點陣夾芯結(jié)構(gòu)的換熱特性[J]. 復(fù)合材料學(xué)報, 2011, 28(4): 185-195

SUN Y G, GAO L. Heat transfer characteristics of composite sandwich structure with lattice cores[J]. Acta Materiae Compositae Sinica, 2011, 28(4): 185-195

[15] GAO L, SUN Y G. Thermal control of composite sandwich structure with lattice truss cores[J]. Journal of Thermophysics and Heat Transfer, 2015, 29(1): 47-54

[16] GAO L, SUN Y G. Fluid flow and heat transfer characteristics of composite lattice core sandwich structures[J]. Journal of Thermophysics and Heat Transfer, 2014, 28(2): 258-269

[17] 高亮, 孫雨果, 叢立新. 桿件傾角對復(fù)合點陣夾芯結(jié)構(gòu)主動換熱性能的影響[J]. 復(fù)合材料學(xué)報, 2013, 30(5): 180-186

GAO L, SUN Y G, CONG L X. Effect of inclination angle on active cooling performance of composite sandwich structure with lattice truss cores[J]. Acta Materiae Compositae Sinica, 2013, 30(5): 180-186

[18] RAKOW J F, WAAS A M. Response of actively cooled metal foam sandwich panels exposed to thermal loading[J]. AIAA Journal, 2007, 45(2): 329-336

[19] FENG S S, KUANG J J, WEN T, et al. An experimental and numerical study of finned metal foam heat sinks under impinging air jet cooling[J]. International Journal of Heat and Mass Transfer, 2014, 77: 1063-1074

[20] FENG S S, KUANG J J, LU T J, et al. Heat transfer and pressure drop characteristics of finned metal foam heat sinks under uniform impinging flow[J]. Journal of Electronic Packaging, 2015, 137(2): 21014

[21] 王儲, 王建華, 杜治能, 等. 冷卻介質(zhì)在層板內(nèi)流動特性研究: 第二部分?jǐn)?shù)值模擬復(fù)雜結(jié)構(gòu)內(nèi)流場[J]. 實驗流體力學(xué), 2007, 21(4): 22-26

WANG C, WANG J H, DU Z N, et al. An investigation of coolant flow performance within laminated plate: II - numerical simulation of flow field of complex structure[J]. Journal of Experiments in Fluid Mechanics, 2007, 21(4): 22-26

[22] HEIMBS S. Virtual testing of sandwich core structures using dynamic finite element simulations[J]. Computational Materials Science, 2009, 45(2): 205-216

[23] WILLIAMS S D, CURRY D M. Thermal protection materials: thermophysical property data[R]. NASA, 1992

[24] ANDERSON J D. 計算流體力學(xué)基礎(chǔ)及其應(yīng)用[M].北京: 機械工業(yè)出版社, 2007: 24-65

[25] BAPANAPALLI S K, MARTINEZ O M, GOGU C, et al. Analysis and design of corrugated-core sandwich panels for thermal protection systems of space vehicles[C]//47thAIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference. Newport, Rhode Island, 2006

[26] 范緒箕. 高速飛行器熱結(jié)構(gòu)分析與應(yīng)用[M]. 北京: 國防工業(yè)出版社, 2009: 39-62

[27] 張靖周, 常海萍. 傳熱學(xué)[M]. 北京: 科學(xué)出版社, 2015: 258-265

（編輯：張艷艷）

Heat transfer and thermal stress analyses of V-pattern folded core sandwich structures subjected to forced convection

ZHOU Chen1, WANG Zhijin1, HOU Tianjiao2

(1. Minister Key Discipline Laboratory of Advanced Design Technology of Aircraft, Nanjing University of Aeronautics and Astronautics; 2. College of Aeronautics, Nanjing University of Aeronautics and Astronautics: Nanjing 210016, China)

The combination of the V-pattern folded cores and the forced convection may apply to the multi-functional design the simultaneous load bearing and the heat dissipation for hypersonic vehicles or high-power electronic devices. The heat transfer and thermal stress characteristics of a V-pattern folded core sandwich panel subjected to the forced convection are investigated numerically to obtain the velocity field, the temperature field, and the structural thermal stress distribution. It is shown that the heat transfer performance of the V-pattern folded panels is evidently improved under the forced convection, however, an elevated pressure drop is also observed. The coolant velocity varies periodically with the geometrical construction and reaches the maximum value close to the peak position of the buckling, where the heat convection is therefore reinforced. The temperature increases along the flow direction and sees a periodic fluctuation. The overall heat transfer performance is refined with the increase of the inlet velocity. The thermal stress as well as the deformation are comparatively large in the face sheet near the heat source and the stress concentration is found at the folding positions.

folded core; forced convection; sandwich structure; heat transfer characteristics; thermal stress

V414.9

A

1673-1379(2017)04-0355-09

10.3969/j.issn.1673-1379.2017.04.003

周晨（1989—），男，博士研究生，主要研究方向為飛行器設(shè)計和熱防護(hù)結(jié)構(gòu)設(shè)計；E-mail: zhouchen@nuaa. edu.cn。指導(dǎo)教師：王志瑾（1963—），女，教授，博士生導(dǎo)師，主要從事飛行器結(jié)構(gòu)設(shè)計、結(jié)構(gòu)優(yōu)化設(shè)計和熱防護(hù)結(jié)構(gòu)設(shè)計等研究；E-mail: zhijin@nuaa.edu.cn。

2017-02-24；

2017-07-11

中央高?；究蒲袠I(yè)務(wù)費專項資金資助項目（編號：NZ2016101）；江蘇省普通高校研究生科研創(chuàng)新計劃資助項目（編號：CXLX13_163）；江蘇高校優(yōu)勢學(xué)科建設(shè)工程資助項目