李征宇,劉 沛

(1.中國(guó)船舶集團(tuán)有限公司第七二三研究所,江蘇 揚(yáng)州 225101;2.海軍青島雷達(dá)聲納修理廠,山東 青島 266100)

0 引 言

高超聲速飛行器在服役過程中會(huì)面臨高溫、噪聲和振動(dòng)等復(fù)雜動(dòng)力學(xué)環(huán)境,研究表明飛行器大面積區(qū)域溫度在750～1 450 ℃之間[1]。高溫對(duì)飛行器動(dòng)力學(xué)特性的影響主要表現(xiàn)在以下2個(gè)方面:一是降低了材料的力學(xué)性能,使結(jié)構(gòu)剛度明顯下降;二是在結(jié)構(gòu)中產(chǎn)生不均勻的熱應(yīng)力,導(dǎo)致結(jié)構(gòu)“剛化”或“軟化”。因此,高溫可能會(huì)使飛行器結(jié)構(gòu)承載能力降低,進(jìn)而導(dǎo)致結(jié)構(gòu)更容易發(fā)生破壞。復(fù)雜的動(dòng)力學(xué)環(huán)境對(duì)高超聲速飛行器的結(jié)構(gòu)設(shè)計(jì)提出了更高要求。近年來,復(fù)合材料以其輕質(zhì)、高比剛度、高比強(qiáng)度和熱穩(wěn)定性好等特點(diǎn)被廣泛應(yīng)用于高超聲速飛行器蒙皮等結(jié)構(gòu)中[2],其中復(fù)合材料加筋板是航空航天領(lǐng)域結(jié)構(gòu)設(shè)計(jì)中常用的典型構(gòu)件之一。

國(guó)內(nèi)外針對(duì)高溫環(huán)境下結(jié)構(gòu)的振動(dòng)特性已有較多相關(guān)研究。Behnke等[3]研究了熱環(huán)境對(duì)金屬薄壁結(jié)構(gòu)低頻和中頻振動(dòng)特性的影響;Malekzadeh等[4]分析了溫變、材料屬性等參數(shù)對(duì)結(jié)構(gòu)固有頻率的影響;王宏宏等[5]研究了加熱狀態(tài)下變厚度導(dǎo)彈翼面模型的振動(dòng)特性;賀旭東等[6]研究了熱應(yīng)力對(duì)飛行器機(jī)翼結(jié)構(gòu)固有頻率的影響;李躍明等[7-8]以X43飛行器為例,分別考慮了氣動(dòng)加熱環(huán)境下材料力學(xué)性能變化和熱應(yīng)力對(duì)結(jié)構(gòu)振動(dòng)特性的影響。目前研究大多針對(duì)傳統(tǒng)的金屬材料,針對(duì)復(fù)合材料的相關(guān)研究較少。此外,由于氣動(dòng)加熱及太陽輻射的影響,高速飛行器在服役過程中通常面臨非均勻的熱載荷[9],而載荷的非均勻程度對(duì)結(jié)構(gòu)的穩(wěn)定性及振動(dòng)特性有較大影響。因此,有必要針對(duì)非均勻溫度載荷下復(fù)合材料結(jié)構(gòu)的振動(dòng)特性開展研究。

本文以復(fù)合材料加筋板為研究對(duì)象,同時(shí)考慮溫度對(duì)材料力學(xué)性能和結(jié)構(gòu)熱應(yīng)力的影響,基于計(jì)及熱效應(yīng)有限元方程研究非均勻溫度載荷對(duì)結(jié)構(gòu)振動(dòng)特性的影響,計(jì)算均勻溫度場(chǎng)及不同溫度梯度的非均勻溫度載荷下結(jié)構(gòu)的熱應(yīng)力分布和振動(dòng)模態(tài),分析非均勻溫度載荷對(duì)結(jié)構(gòu)固有頻率和模態(tài)振型的影響。

1 理論基礎(chǔ)

1.1 溫度對(duì)材料力學(xué)性能的影響

溫度對(duì)材料力學(xué)性能的影響主要表現(xiàn)在對(duì)結(jié)構(gòu)彈性模量的改變。根據(jù)若干溫度下給定的彈性模量值,對(duì)每一個(gè)溫度區(qū)間線性插值,可得溫度t下的彈性模量Ee(t)為:

Ee(t)=Ee(t0)+ΔEe

(1)

對(duì)于線彈性材料,單元彈性矩陣De由彈性模量E及泊松比μ決定,因此溫度t下單元彈性矩陣De可以表示為初始溫度下的彈性矩陣與溫差引起的附加矩陣之和,即:

De(t)=De(t0)+ΔDe

(2)

單元?jiǎng)偠染仃嚍?

(3)

式中:Be為單元幾何矩陣。

將式(2)代入式(3)中,可得溫度t下單元?jiǎng)偠染仃?

(4)

Ke(t)=Ke(t0)+KT

(5)

1.2 熱應(yīng)力對(duì)結(jié)構(gòu)剛度的影響

板殼類結(jié)構(gòu)在面內(nèi)力作用下其彎曲剛度會(huì)發(fā)生變化,而熱應(yīng)力屬于面內(nèi)力,因此板殼類結(jié)構(gòu)在熱應(yīng)力作用下將會(huì)產(chǎn)生附加應(yīng)力剛度矩陣。

在板單元中,與平板橫向撓曲函數(shù)u=u(x,y,z)相關(guān)的薄膜應(yīng)變?yōu)?

耳鼻咽喉科ESS手術(shù)后術(shù)腔黏膜的恢復(fù)一直是臨床醫(yī)生關(guān)注的重點(diǎn)，其中瘢痕和粘連的形成是影響預(yù)后的重要因素。少數(shù)瘢痕體質(zhì)的患者術(shù)區(qū)瘢痕化嚴(yán)重影響鼻腔黏膜的正常功能，使手術(shù)不能達(dá)到預(yù)期的治療效果；又有部分患者因竇口黏膜瘢痕收縮導(dǎo)致竇腔狹小，甚至再度封閉，分泌物無法排出，反復(fù)擴(kuò)張竇口又會(huì)損傷黏膜纖毛功能并加重瘢痕的增生，最后導(dǎo)致鼻竇炎復(fù)發(fā)。粘連容易引起鼻腔狹窄、結(jié)構(gòu)異常、氣流動(dòng)力學(xué)改變、通氣異常等，這既不利于鼻腔鼻竇分泌物的排出，又會(huì)改變鼻腔的正常生理結(jié)構(gòu)，引起患者術(shù)后不適。耳鼻咽喉科常采用反復(fù)術(shù)后換藥、切除瘢痕、分解粘連及糖皮質(zhì)激素全身或術(shù)腔局部給藥等方式來解決這類問題。

(6)

當(dāng)面內(nèi)力不隨撓度變化時(shí),薄板的應(yīng)變能為:

u,xu,yNxy+u,yu,zNyz+u,zu,xNzx)dxdydz=

(7)

式中:Nx、Ny、Nz、Nxy、Nxz、Nyz分別為各方向上由熱荷載引起的面內(nèi)力。

用形函數(shù)N及板單元節(jié)點(diǎn)彎曲位移d表示板的撓度u=Nd,則有:

(8)

(9)

由各單元的應(yīng)力矩陣可組裝得到系統(tǒng)的總體應(yīng)力剛度矩陣:

(10)

式中:Kσ即為熱應(yīng)力引起的附加剛度矩陣。

2 數(shù)值仿真

2.1 研究對(duì)象

為研究非均勻溫度載荷對(duì)復(fù)合材料結(jié)構(gòu)的振動(dòng)特性影響,選取幾何尺寸為300 mm×200 mm×5 mm的典型復(fù)合材料加筋板為研究對(duì)象,采用殼單元建立其有限元模型,單元尺寸為5 mm,所建立的有限元模型如圖1所示。為方便后續(xù)描述,將加筋板的平板部分分為A板、B板和C板三部分。復(fù)合材料屬性如表1所示。在分析過程中設(shè)置結(jié)構(gòu)模態(tài)阻尼比為0.02。

圖1 復(fù)合材料加筋板有限元模型

2.2 荷載與邊界條件

在復(fù)合材料加筋板的四邊施加簡(jiǎn)支約束。為研究非均勻溫度載荷對(duì)加筋板響應(yīng)特性的影響,計(jì)算共選取了6種不同溫度梯度的熱載荷,為方便對(duì)比分析,6種溫度載荷的平均溫度均為900 ℃,溫度梯度dt分別為0 ℃(均勻溫度場(chǎng))、100 ℃、200 ℃、300 ℃、400 ℃、500 ℃。溫度梯度為300 ℃載荷作用下加筋板結(jié)構(gòu)的溫度分布如圖2所示。

圖2 加筋板溫度場(chǎng)分布

2.3 非均勻溫度場(chǎng)對(duì)結(jié)構(gòu)熱應(yīng)力的影響

常溫及不同非均勻熱載荷條件下結(jié)構(gòu)的前10階固有頻率如表2所示。由表2可知,在溫度梯度為0 ℃的熱載荷即均勻溫度載荷作用下,結(jié)構(gòu)的熱應(yīng)力呈現(xiàn)對(duì)稱分布,隨著熱載荷非均勻程度的增大,A板及C板的熱應(yīng)力差值越來越大。C板處的熱應(yīng)力遠(yuǎn)大于A板,這主要由于C板的溫度較高,高溫使結(jié)構(gòu)產(chǎn)生了較大的熱變形,在約束條件下產(chǎn)生了較大熱應(yīng)力。相比于均勻熱載荷作用下的結(jié)構(gòu),在溫度梯度為500 ℃的熱載荷作用下,結(jié)構(gòu)的熱應(yīng)力幅值增加了31.8%。

表2 結(jié)構(gòu)熱應(yīng)力

2.4 非均勻溫度場(chǎng)對(duì)結(jié)構(gòu)固有頻率的影響

常溫及不同非均勻熱載荷條件下結(jié)構(gòu)的前10階固有頻率如表3所示。由表3可知,相比于常溫環(huán)境,熱載荷作用下的結(jié)構(gòu)其前10階固有頻率均有不同程度的上升;相比于常溫環(huán)境下的結(jié)構(gòu),高溫均勻熱載荷作用下的結(jié)構(gòu),前10階固有頻率的平均變化值為7.1%。這主要是由于高溫引起的材料力學(xué)性能的變化,溫度的升高增加了復(fù)合材料的彈性模量,從而增大了結(jié)構(gòu)的剛度,使得結(jié)構(gòu)的固有頻率升高。隨著熱載荷非均勻程度的增加,結(jié)構(gòu)的固有頻率變化整體呈下降趨勢(shì),但總體而言,非均勻溫度載荷對(duì)結(jié)構(gòu)固有頻率的影響較小。相比于均勻熱載荷作用下的結(jié)構(gòu),在溫度梯度為500 ℃的熱載荷作用下,結(jié)構(gòu)的前10階固有頻率的平均變化值為4.3%。這主要是由于溫度的非均勻分布對(duì)結(jié)構(gòu)的熱應(yīng)力影響造成的,結(jié)構(gòu)的熱應(yīng)力隨著溫度的非均勻程度的增大而增加,從而使得結(jié)構(gòu)“剛化”,提高了結(jié)構(gòu)的固有頻率。

表3 不同非均勻熱載荷下結(jié)構(gòu)的固有頻率 Hz

2.5 非均勻溫度場(chǎng)對(duì)結(jié)構(gòu)模態(tài)振型的影響

不同非均勻熱載荷條件下結(jié)構(gòu)的前10階模態(tài)振型如表4所示。由表4可知,結(jié)構(gòu)的前10階模態(tài)振型中,第7階是筋條振動(dòng)為主的模態(tài)振型,第8階是筋條及平板共同振動(dòng)為主的模態(tài)振型,其他階次為平板振動(dòng)為主的模態(tài)振型。對(duì)于前10階模態(tài)振型而言,熱載荷非均勻程度的變化并未使得結(jié)構(gòu)的模態(tài)振型發(fā)生順序互換或階次跳躍。非均勻溫度載荷對(duì)結(jié)構(gòu)模態(tài)中平板振動(dòng)為主的模態(tài)振型影響較大。以結(jié)構(gòu)第一階振型為例進(jìn)行說明,均勻熱載荷下結(jié)構(gòu)第一階呈現(xiàn)對(duì)稱以A和C 2個(gè)面板的中心振動(dòng)為主,隨著熱載荷非均勻程度的增加,C面板的模態(tài)位移逐漸減小,這主要是由于其溫度高于A面板,高溫使得材料的彈性模量增加從而增大了面板的剛度,類似結(jié)果也存在于其他階次的模態(tài)振型中。此外,不同溫度梯度的熱載荷作用下,結(jié)構(gòu)的第7階和第8階模態(tài)振型并未發(fā)生明顯變化,這說明對(duì)結(jié)構(gòu)模態(tài)中筋條振動(dòng)為主的模態(tài)振型的影響較小。

表4 不同非均勻熱載荷下結(jié)構(gòu)的模態(tài)振型

3 結(jié)束語

本文以復(fù)合材料加筋板為研究對(duì)象,同時(shí)考慮了溫度對(duì)材料力學(xué)性能和結(jié)構(gòu)熱應(yīng)力的影響,基于計(jì)及熱效應(yīng)有限元方程研究了非均勻溫度載荷對(duì)結(jié)構(gòu)振動(dòng)特性的影響,計(jì)算了均勻溫度場(chǎng)及不同溫度梯度的非均勻溫度場(chǎng)下的結(jié)構(gòu)熱應(yīng)力分布和振動(dòng)模態(tài),分析了非均勻溫度載荷對(duì)結(jié)構(gòu)固有頻率和模態(tài)振型的影響。研究結(jié)果表明:

(1) 非均勻溫度載荷的溫度梯度對(duì)結(jié)構(gòu)的熱應(yīng)力分布有著重要影響,結(jié)構(gòu)熱應(yīng)力的幅值隨熱載荷的溫度梯度的增加而增加;

(2) 熱環(huán)境下復(fù)合材料加筋板固有頻率隨溫度載荷的非均勻程度的增加整體呈下降趨勢(shì);

(3) 非均勻溫度載荷對(duì)結(jié)構(gòu)模態(tài)中平板振動(dòng)為主的模態(tài)振型影響較大,對(duì)結(jié)構(gòu)模態(tài)中筋條振動(dòng)為主的模態(tài)振型的影響較小。

上述研究結(jié)果對(duì)于高超聲速飛行器的設(shè)計(jì)與研制有一定的參考意義。