周家來 李縣準(zhǔn)
(蘭州大學(xué)土木工程與力學(xué)學(xué)院,甘肅 蘭州 730000)
在力學(xué)分析中,研究對象常被視為理想材料,但在實(shí)際研究中不可能有符合這些條件的材料,需要運(yùn)用試驗(yàn)或建模對其理論計(jì)算值進(jìn)行校驗(yàn)。繁瑣而又千篇一律的計(jì)算,使得人們?nèi)菀缀雎粤W(xué)建模和力學(xué)原理的學(xué)習(xí)。通過ANSYS建模,人們可以了解構(gòu)件或者結(jié)構(gòu)的變形、應(yīng)力分布規(guī)律等,對加強(qiáng)人們的建模與力學(xué)原理的思想具有積極的推動作用。國內(nèi)外學(xué)者重視運(yùn)用ANSYS分析結(jié)果來檢驗(yàn)理論方法的準(zhǔn)確性。2005年,楊笑冬[1]運(yùn)用ANSYS對懸臂梁進(jìn)行模態(tài)分析,結(jié)果表明ANSYS對模態(tài)分析具有可靠性。2014年,Yaylac1等人[2]運(yùn)用ANSYS、彈性理論和積分變換法分析退讓接觸問題,理論計(jì)算值與ANSYS分析結(jié)果符合良好。本文詳盡地闡述ANSYS對懸臂工字梁不同位置集中荷載作用下的靜力響應(yīng)分析;利用材料力學(xué)撓度公式計(jì)算L/4,L/2,3L/4和懸臂端的撓度,通過對比分析,得到材料力學(xué)撓度公式的準(zhǔn)確性。通過完整的ANSYS建模和分析不同工況下懸臂工字梁的力學(xué)性能,進(jìn)一步加深人們對力學(xué)基本概念的理解。
圖1為一工字懸臂梁。根據(jù)《鋼結(jié)構(gòu)設(shè)計(jì)規(guī)范》[3],梁的材料可選取:14號熱軋普通工字鋼,其截面參數(shù)如表1所示。梁跨L=1 800 mm,Q235鋼的彈性模量E=2.0×105N/mm2,泊松比μ=0.33。
表1 工字型截面參數(shù)
型號尺寸/mm截面特性hbtwtIx/cm4Wx/cm3ix/cmⅠ14140805.59.1712101.75.75
Q235鋼的許用應(yīng)力[σ]=160 MPa[4],懸臂梁正應(yīng)力強(qiáng)度條件為:
由上式得:Mmax=16.272 kN·m。由此可得集中力在不同位置時,所能施加的最大集中力,如表2所示。
表2 最大集中力
定義圖形界面過濾參數(shù)是為了得到一個相對簡潔的分析菜單。拾取菜單Main Menu中Preferences,在彈出過濾對話框中選擇“Structural”項(xiàng),單擊“OK”按鈕。
本文普通熱軋工字鋼采用六面體Solid45單元,而且只設(shè)置一種材料模型,為統(tǒng)一單位,彈性模量單位為N/mm2。設(shè)置材料物理參量之后,根據(jù)所選材料,確定14個關(guān)鍵點(diǎn):1(-40,0,0),2(0,0,0),3(40,0,0),4(-40,9.1,0),5(2.75,9.1,0),6(2.75,130.9,0),7(40,130.9,0),8(40,140,0),9(0,140,0),10(-40,140,0),11(-40,130.9,0),12(-2.75,130.9,0),13(-2.75,9.1,0),14(-40,9.1,0)。采用“自底而上”的方法,通過確定關(guān)鍵點(diǎn),依次連成線生成面,最后建立跨度是1 800 mm的懸臂梁三維實(shí)體(見圖2)。
映射網(wǎng)格劃分是對規(guī)整模型的一種規(guī)整網(wǎng)格劃分方法,要求體的形狀遵循一定的規(guī)則。工字鋼梁并非磚形、楔形,故在建立懸臂梁三維實(shí)體過程中,工字梁截面通過拉伸生成三維復(fù)雜的實(shí)體。在某個方向上,拓?fù)湫问绞冀K保持一致。對于處理復(fù)雜幾何實(shí)體,掃掠網(wǎng)格劃分是一個非常好的方法。它經(jīng)過簡單的切分處理,就自動形成規(guī)整的六面體網(wǎng)格(Hex/Wedge),比映射網(wǎng)格劃分具有更大的優(yōu)勢和靈活性。因此,本文采用掃掠網(wǎng)格劃分(見圖3)。
劃分網(wǎng)格之前要設(shè)定工字懸臂梁單元屬性,包括單元類型、實(shí)常數(shù)、材料性質(zhì)等。其次,設(shè)定單元大小。點(diǎn)擊“Mesh Tool”,彈出的對話框,在“Size Controls”的Global一欄點(diǎn)擊Set,然后在彈出的對話框Size一欄中填10,表示每10 mm劃分一個單元。
工字懸臂梁一端固定,一端自由。GUI操作過程為:Main Menu→Preprocessor→Define loads→Loads→Apply→Structural→Displacement→On Areas。彈出對話框,然后點(diǎn)取一端截面,點(diǎn)擊“OK”,彈出另一對話框(見圖4)。選擇“All DOF”,確定該截面XYZ方向全部約束。通過上述操作可得到如圖5所示的邊界約束條件。
集中荷載施加的位置有:L/4,L/2,3L/4,懸臂端。在操作的過程中,根據(jù)集中荷載施加的位置建立兩個不同長度的實(shí)體。例如當(dāng)集中荷載施加于L/2時,截面Z軸負(fù)方向拉伸900 mm,Z軸正方向拉伸900 mm,生成兩個緊靠又相互獨(dú)立的實(shí)體,運(yùn)用“Add”的Booleans操作,實(shí)現(xiàn)“合二為一”。圖6列出了集中荷載P作用于L/2和懸臂端的情況。
該工字懸臂梁為實(shí)腹式壓彎構(gòu)件,故需要進(jìn)行集中荷載作用下的局部穩(wěn)定驗(yàn)算(所選的材料滿足整體穩(wěn)定驗(yàn)算、強(qiáng)度驗(yàn)算和剛度驗(yàn)算)。該工字鋼外伸翼緣板寬厚比:
經(jīng)計(jì)算,其局部穩(wěn)定驗(yàn)算滿足要求,故本算例可以不用設(shè)置加勁肋和墊塊,集中載荷可直接施加在相應(yīng)位置一點(diǎn)處,這不會引起局部失穩(wěn)。
由于各工況的模擬結(jié)果表現(xiàn)出來的規(guī)律大同小異,故只列出了懸臂端在集中載荷P=5 kN作用下的位移、應(yīng)變(見圖7,圖8)。
從圖7,圖8可以看出,腹板截面應(yīng)力應(yīng)變小,中部應(yīng)變甚至為0;上下翼緣應(yīng)力應(yīng)變大,而且大小相近,在整個截面的受力基本
符合平截面假定。
在建模過程中,整體坐標(biāo)系分別設(shè)置在L/4,L/2,3L/4,懸臂端,實(shí)現(xiàn)對荷載的精確施加,方便在云圖中找到相應(yīng)位置的撓度,結(jié)合撓度計(jì)算公式[4],相應(yīng)位置的-Y方向位移如表3所示。
表3 撓度對比
表3數(shù)據(jù)顯示,集中載荷接近Pmax時,模擬結(jié)果和解析解均沒有超出撓度限值[L/200]=9 mm。當(dāng)載荷作用于L/4,L/2,3L/4和懸臂端時,解析解與模擬結(jié)果相對誤差的平均值分別為0.15/0.11,0.11/0.07,0.07/0.12,0.11/0.13,說明解析解準(zhǔn)確性良好。
通過ANSYS建模和材料力學(xué)撓度計(jì)算公式,得到了八個工況的撓度模擬值與計(jì)算值,兩者的相對誤差較小,解析解有較好的準(zhǔn)確性。通過將完整的ANSYS建模引入《材料力學(xué)》的學(xué)習(xí)中,對加強(qiáng)人們的建模與力學(xué)原理的思想具有積極的推動作用,具體表現(xiàn)為:1)通過ANSYS模擬平臺,引導(dǎo)人們樹立建模思想,提高人們的工程意識。2)加強(qiáng)人們對基本力學(xué)概念的理解,提高空間想象力,為人們創(chuàng)新思維的拓展提供空間。
[1] 楊笑冬.基于ANSYS的懸臂梁模態(tài)分析[J].機(jī)電一體化,2008,14(6):58-60.
[2] Yaylac1 M,?ner E,Birinci A.Comparison between Analytical and ANSYS Calculations for a Receding Contact Problem[J].Journal of Engineering Mechanics,2014,140(9):70-71.
[3] GB 50017—2003,鋼結(jié)構(gòu)設(shè)計(jì)規(guī)范[S].
[4] 孫訓(xùn)方.材料力學(xué)(Ⅰ)[M].第5版.北京:高等教育出版社,2009:20-125.