徐志強(qiáng),任毅斌,劉云剛,胡國(guó)強(qiáng),宋小雨
(中鋁材料應(yīng)用研究院有限公司,北京102209)
薄壁梁結(jié)構(gòu)作為客車(chē)結(jié)構(gòu)的主要承力部件,在車(chē)身結(jié)構(gòu)中有著十分重要的作用,會(huì)影響整車(chē)的力學(xué)性能,因此研究薄壁梁的剛度對(duì)提高整車(chē)的性能有著非常重要的意義[1-3]。
一般在客車(chē)中采用的薄壁梁多為矩形梁[4],因此,本文將從理論及有限元兩方面分析矩形薄壁梁結(jié)構(gòu)的截面幾何參數(shù)和材料特性對(duì)其彎曲和扭轉(zhuǎn)剛度的影響,并通過(guò)算例對(duì)不同壁厚鋁合金矩形薄壁梁的抗彎及抗扭剛度進(jìn)行理論及數(shù)值計(jì)算,得出理論與數(shù)值相同的計(jì)算結(jié)果。
考慮到薄壁梁在客車(chē)上的設(shè)計(jì),梁的兩端均會(huì)與其他的部件有機(jī)械連接或者焊接,因此,在單獨(dú)分析薄壁梁時(shí),可假定薄壁梁兩端固支,研究在兩端固支的約束條件下梁上任意一點(diǎn)的剛度值。
一長(zhǎng)度為L(zhǎng)的兩端固支梁如圖1所示,在距左端x處受一垂向載荷P作用。
由平衡方程知
式中,F(xiàn)AX、FAY、FCX、FCY分別為A、C點(diǎn)的X、Y向的支座反力;MA、MC分別為A、C點(diǎn)的彎矩。
公式中(1)含有6個(gè)未知數(shù),而平衡方程只有3個(gè),因此多出3個(gè)約束,屬于三次超靜定問(wèn)題。要求解上述問(wèn)題,需要變形協(xié)調(diào)方程。由于該梁為兩端固支,固支梁的兩端撓度與轉(zhuǎn)角均為0,即:
由公式(1)、(2)可得
其中,ωB為B點(diǎn)的撓度。
B點(diǎn)的撓度可以通過(guò)對(duì)上述固支梁進(jìn)行轉(zhuǎn)化解得,即上述固支梁轉(zhuǎn)變?yōu)閷點(diǎn)固支,C點(diǎn)的約束采用MC和FCY等效,則B點(diǎn)處的撓度可以看作分別受到P、MC、FCY以后疊加在一起,即
其中,ωP為載荷P在B點(diǎn)產(chǎn)生的撓度,ωMC為彎矩MC在B點(diǎn)產(chǎn)生的撓度,ωFCY為載荷FCY在B點(diǎn)產(chǎn)生的撓度,其值由材料力學(xué)可知[5-6]:
式中,E、I分別為材料的彈性模量及截面對(duì)x軸的慣性矩。
將公式(6)、(7)、(8)同時(shí)代入到公式(5)中,可得
將式(9)代入公式(4),可得B點(diǎn)的抗彎剛度為:
對(duì)于矩形薄壁梁,如圖2所示,其抗彎剛度同樣可以用公式(12)來(lái)表示。
圖2 矩形薄壁斷面
其中矩形梁的慣性矩為:
其中,b,h分別為薄壁梁截面厚度t中線(xiàn)處的寬和高。
將公式(12)代入公式(11),可得矩形薄壁梁的抗彎剛度為:
考慮AC梁兩端固支,如下圖所示,假定AC梁的剪切模量為G,對(duì)軸心的慣性矩為Iρ。
由靜力學(xué)平衡方程可知:
公式(15)中含有2個(gè)未知數(shù),而平衡方程只有1個(gè),多出1個(gè)約束,屬于一次超靜定問(wèn)題,因此要求解上述問(wèn)題,需要變形協(xié)調(diào)方程。由于該梁為兩端固支,固支梁AC的扭轉(zhuǎn)角為0,即:
將AC梁看成AB段和BC段,AB段扭矩為T(mén)A,BC段扭矩為T(mén)C,AC的扭轉(zhuǎn)角等于AB與BC扭轉(zhuǎn)角之和,即
由于AB段的扭轉(zhuǎn)方向和BC段的扭轉(zhuǎn)方向相反,對(duì)于任意閉口薄截面桿件,如圖4所示。
圖4 任意閉口薄壁斷面
根據(jù)扭轉(zhuǎn)角的計(jì)算公式[6],則有AB段與BC段的扭轉(zhuǎn)角分別為
式中,s為中線(xiàn)長(zhǎng)度,t為壁厚,A為斷面中線(xiàn)所圍成的面積。
將公式(18)代入公式(17),可得:
聯(lián)立公式(15)、(19),可得:
則梁AC上一點(diǎn)B處的扭轉(zhuǎn)角為:
因此,梁AC上受力點(diǎn)B處的扭轉(zhuǎn)剛度等于B處的扭矩除以B處的扭轉(zhuǎn)角,即:
對(duì)于圖2所示的等壁厚矩形薄壁梁,則:
將公式(23)、(24)代入到公式(21)、(22)中,可得受力點(diǎn)B處的扭轉(zhuǎn)角及扭轉(zhuǎn)剛度,即:
設(shè)置一算例,假定一長(zhǎng)度L=1 000 mm的矩形薄壁梁,中線(xiàn)寬度b=60 mm,中線(xiàn)高度h=100 mm,壁厚為t,t取1~9 mm中的整數(shù),矩形薄壁梁材料為鋁合金,其彈性模量為70 GPa,密度為2 700 kg/m3,泊松比為0.33。
為了所得到的抗彎剛度具有普遍性,矩形薄壁梁所受載荷不施加在梁的中間位置,設(shè)定矩形薄壁梁在x=400 mm處受到P=100 kN的載荷,兩端固定。
2.1.1 理論計(jì)算
將上述已知條件代入到公式(13)中,可得出在受力點(diǎn)處的彎曲撓度及抗彎剛度值KC,具體值見(jiàn)于表1。
2.1.2 有限元數(shù)值計(jì)算
采用一維梁?jiǎn)卧M(jìn)行建模,單元尺寸設(shè)置為10 mm,按照上述條件設(shè)置梁的截面參數(shù),在兩端施加全約束,并在x=400 mm處施加P=100 kN的集中載荷,如圖5所示。
圖5 彎曲有限元模型
通過(guò)計(jì)算可得出在受力點(diǎn)處的位移UY,如圖6所示。而抗彎剛度值可由下列方程式求得。
圖6 不同壁厚受力點(diǎn)處的位移云圖
式中,KBI為數(shù)值計(jì)算抗彎剛度。將不同厚度的矩形薄壁梁計(jì)算出的位移及抗彎剛度列于表1。
表1 矩形薄壁梁彎曲的理論及數(shù)值解
從表1可以看出,通過(guò)理論計(jì)算得出的抗彎剛度與通過(guò)有限元數(shù)值計(jì)算得出的抗彎剛度完全一致,從而驗(yàn)證了抗彎剛度的理論計(jì)算公式的有效性。
同矩形薄壁梁的抗彎剛度:設(shè)定矩形薄壁梁在x=400 mm處受到T=100 kN·mm的扭矩,梁的兩端固定。
2.2.1 理論計(jì)算
將上述已知條件代入到公式(26)中,可得出在受力點(diǎn)處的扭轉(zhuǎn)剛度值KT,具體值見(jiàn)表2。
2.1.2 有限元數(shù)值計(jì)算
采用三維實(shí)體單元進(jìn)行建模,按照上述條件設(shè)置梁的截面參數(shù),在兩端施加全約束,并在x=400 mm處施加T=100 kN·mm的扭矩,如圖7所示。
圖7 扭轉(zhuǎn)有限元模型
通過(guò)數(shù)值計(jì)算,得出各厚度的位移云圖,如圖8所示。而扭轉(zhuǎn)剛度為:
圖8 不同壁厚的位移云圖
式中,θTI為在受力截面處的扭轉(zhuǎn)角,KTI為數(shù)值計(jì)算出的扭轉(zhuǎn)剛度。
θTI可通過(guò)選取中間截面上沿軸心對(duì)稱(chēng)的節(jié)點(diǎn)相對(duì)位移除以?xún)烧呔嚯x得到[7],然后將計(jì)算出的扭轉(zhuǎn)角代入到公式(28),可得到不同厚度鋁合金矩形薄壁梁的扭轉(zhuǎn)剛度,見(jiàn)表2。
表2 矩形薄壁梁扭轉(zhuǎn)的理論及數(shù)值解
從表2可以看出,通過(guò)理論計(jì)算得出的抗扭剛度與通過(guò)有限元數(shù)值計(jì)算得出的抗扭剛度完全一致,從而驗(yàn)證了抗扭剛度的理論計(jì)算公式。
對(duì)兩端固支的任意截面梁進(jìn)行了彎曲及扭轉(zhuǎn)的理論研究,對(duì)加載處的抗彎剛度及抗扭剛度進(jìn)行了理論推導(dǎo),并以鋁合金矩形薄壁梁作為算例,對(duì)其抗彎剛度及抗扭剛度進(jìn)行了理論與數(shù)值計(jì)算,結(jié)論如下:
(1)通過(guò)理論研究,得出了兩端固支的任意截面梁在受力處的抗彎剛度計(jì)算公式,并通過(guò)該理論公式推出了鋁合金薄壁矩形截面梁的抗彎剛度計(jì)算公式。
(2)通過(guò)理論推導(dǎo)得出了任意閉口截面梁在兩端固支的情況下其受力處的抗扭剛度的計(jì)算公式,該理論計(jì)算公式適用于任意閉口截面梁。然后通過(guò)該理論計(jì)算公式推出了薄壁矩形截面梁在加載處的抗扭剛度。
(3)通過(guò)算例,在兩端固支的鋁合金矩形薄壁梁400 mm處位置進(jìn)行加載,對(duì)不同厚度的薄壁矩形梁的抗彎剛度及扭轉(zhuǎn)剛度進(jìn)行了理論及數(shù)值計(jì)算。計(jì)算結(jié)果表明,對(duì)不同厚度的薄壁矩形梁,其抗彎剛度及抗扭剛度理論值與數(shù)值計(jì)算值保持一致,從而驗(yàn)證了抗彎及抗扭剛度的理論計(jì)算公式的有效性,為客車(chē)的矩形梁設(shè)計(jì)提供了一定的基礎(chǔ)。