閆 棣， 蘇 祺， 李四平

(上海交通大學(xué) 船舶海洋與建筑工程學(xué)院， 上海 200240)

固體力學(xué)研究領(lǐng)域中，工程結(jié)構(gòu)的非線性屈曲問題一直是非?；钴S的話題[1].在使用通用有限元軟件進(jìn)行屈曲模擬的過程中，一些文獻(xiàn)[2-3]直接使用軟件中提供的線性屈曲分析模塊進(jìn)行屈曲臨界載荷的預(yù)測與屈曲模態(tài)的判斷.線性屈曲分析方法相當(dāng)于彈性屈曲分析，然而材料的非線性、結(jié)構(gòu)形狀的初始缺陷以及幾何大變形使得實際結(jié)構(gòu)并不在理論彈性范圍內(nèi)屈曲.因此，特征值屈曲分析所得的解為非保守解，一般不能用于實際工程的分析[4-5].除此之外，線性屈曲分析方法對材料非線性等問題無法解決，有一定的局限性.一些文獻(xiàn)[6-7]使用線性屈曲分析中所得的一階模態(tài)或前幾階模態(tài)進(jìn)行缺陷的引入，而在引入缺陷時所用的比例因子參數(shù)沒有定論.這種方法通過人為調(diào)控比例因子的大小可能得到不同的結(jié)果，計算結(jié)果仍需要通過進(jìn)一步的試驗比對來調(diào)整參數(shù)，而且在線性屈曲方法無法計算的材料非線性等問題上[8]，繼續(xù)使用此方法進(jìn)行屈曲的分析，所得結(jié)論是不可信的.

本文提出一種更為普適的模擬方法，即通過隨機(jī)的方法對理想結(jié)構(gòu)模型引入初始缺陷進(jìn)行屈曲的分析.這種方法不僅能夠得到上述文獻(xiàn)中所用方法計算得到的多階屈曲模態(tài)和其對應(yīng)臨界力，而且將特征值問題轉(zhuǎn)化為幾何非線性靜力變形問題進(jìn)行計算，克服了線性屈曲方法對材料非線性等問題的局限性，獲得了更多中間加載過程的信息.此外，以彈性地基上梁的屈曲問題和橫向纖維作用下穿透型矩形脫層屈曲問題為例，計算討論隨機(jī)缺陷法的普適性和以上文獻(xiàn)中所使用方法的局限性.同時，通過對比不同大小缺陷對計算結(jié)果的影響，提出合理的隨機(jī)缺陷引入方法.

1 隨機(jī)缺陷法

1.1 理論分析

初始缺陷是指與理論計算模型相比，實際構(gòu)件在尚未受到載荷作用前就已經(jīng)存在的各種缺陷，其中包含力學(xué)初始缺陷以及幾何初始缺陷.力學(xué)初始缺陷主要是由各種原因引起的初始?xì)堄鄳?yīng)力；而幾何初始缺陷則包含整體初彎曲，初偏心，局部初始缺陷，等等.

考慮如圖1所示的兩端鉸支的細(xì)長桿，當(dāng)它處于微彎曲平衡的臨界狀態(tài)時，若桿內(nèi)應(yīng)力不超過材料的比例極限，則壓桿任一截面m-n處的彎曲變形與截面彎矩的關(guān)系為

(1)

式中：E為桿件彈性模量；I為桿件截面慣性矩；M(x)為桿件在距離下支座x處的由外力F引起的彎矩.

圖1 兩端鉸支細(xì)長桿Fig.1 Schematic of hinged-hinged slender beam

考慮模型有初始缺陷，由于初始缺陷是一個隨機(jī)的無窮小量，忽略高階小量，做疊加處理得

(2)

式中：M0(x)為初始缺陷引起的初彎矩，是一個隨機(jī)的非連續(xù)小量.由于M(x)=-Fy，所以式(2)可寫為

(3)

(4)

解二階非線性常系數(shù)微分方程，其通解及特解分別為

Y=asinkx+bcoskx

(5)

(6)

式中：D為微分算子.代入桿端邊界條件(x=0時，y=0，M0(0)=0；x=l時，y=0，M0(l)=0)，可以得到式(4)的解為

(7)

(8)

即當(dāng)M0(x)趨近于0時，式(7)同式(1)的解等價，表示當(dāng)加入的隨機(jī)缺陷足夠小時并不會影響模型整體的屈曲結(jié)果.

然而理想不含初始缺陷的構(gòu)件在外力作用下并不會發(fā)生分支點屈曲，分支點屈曲發(fā)生的原因也正是由于實際構(gòu)件中包含一定的初始缺陷.在通用有限元軟件線性屈曲分析中，將屈曲問題作為特征值問題進(jìn)行求解.而隨機(jī)缺陷法，通過隨機(jī)的方法在理想構(gòu)件中引入缺陷將屈曲作為幾何非線性靜力計算問題進(jìn)行求解，這更符合實際情況.

在對模型施加缺陷時，隨機(jī)缺陷法強(qiáng)調(diào)隨機(jī)性，充分模擬實際構(gòu)件中缺陷發(fā)生的隨機(jī)性.缺陷可以是構(gòu)件的幾何缺陷，也可以是力學(xué)初始缺陷，引入這些小缺陷的目的就是觸發(fā)構(gòu)件屈曲.缺陷大小和隨機(jī)性是計算結(jié)果精度的主要影響因素.

1.2 有限元實現(xiàn)方法的探討

本文算例中采用隨機(jī)選取模型網(wǎng)格節(jié)點進(jìn)行初始偏移的方法隨機(jī)地添加缺陷，根據(jù)計算結(jié)果可以得知偏移的大小取梁截面徑向尺寸的千分之一或更小為宜.這是一種簡單方便的幾何缺陷引入方法.除此之外，本文也嘗試通過程序隨機(jī)選取結(jié)構(gòu)中某一微小單元并賦予較小的模量來模擬材料的空心幾何缺陷.力學(xué)缺陷上，可以使用虛擬的觸發(fā)載荷來啟動不穩(wěn)定性，也可以通過不均勻的溫度場致使屈曲的發(fā)生.

在保證隨機(jī)性的前提下，引入初始缺陷的方法有很多，其根本的目的是為了觸發(fā)屈曲.初始缺陷越小，有限元模型越接近理想模型，計算結(jié)果與理論結(jié)果也越接近，但此時為觸發(fā)屈曲，模型的計算量也會相應(yīng)增加.所以有限元計算中觸發(fā)屈曲時，初始缺陷的大小、模型的計算量以及計算結(jié)果的精度需要有一個平衡的過程.

2 兩端鉸支細(xì)長桿多階屈曲模擬

設(shè)有長l=100 cm，截面面積A=4 cm2，慣性矩Ix=Iy=1.33 cm4的兩端鉸支細(xì)長桿件，長細(xì)比λ=173.2，如圖2所示.材料的彈性模量E=200 GPa.

圖2 兩端鉸支壓桿Fig.2 Schematic of hinged-hinged compression bar

使用通用有限元軟件進(jìn)行建模，采用梁單元進(jìn)行模擬，并將梁沿軸向劃分為100個單元.通過線性屈曲方法得到前5階屈曲模態(tài)及臨界載荷，結(jié)果以第1階和第2階模態(tài)為例，如圖3所示，其中u2表示桿件徑向位移.參照Euler公式

(9)

式中：Fcr為臨界力；n為屈曲模態(tài)數(shù).線性屈曲方法模擬結(jié)果與Euler公式理論結(jié)果(Fcr1=26.32 kN，F(xiàn)cr2=105.28 kN)較為吻合.

圖3 線性屈曲方法所得壓桿屈曲前兩階模態(tài)Fig.3 The first two order modes of compression bar by linear buckling method

圖4 隨機(jī)缺陷法計算壓桿屈曲第1階模態(tài)Fig.4 The first mode of compression bar by stochastic defect method

使用隨機(jī)缺陷法進(jìn)行模擬.通過編寫程序操作計算文件，隨機(jī)選取桿件上節(jié)點施加10-5R0(R0為橫截面徑向尺寸)的橫向偏移，加載位移載荷，使其從0開始緩慢增長，其第1階屈曲模態(tài)如圖4所示.

為保證隨機(jī)性不影響最終結(jié)果，進(jìn)行了多次計算，每次選取的偏移節(jié)點均由程序隨機(jī)生成，所得結(jié)果相同.在使用節(jié)點偏移法的同時，嘗試使用空心幾何缺陷，通過實體單元進(jìn)行建模，隨機(jī)選擇模型中某一位置劃分出微小的單元塊，并賦予較低的材料模量來模擬空心，載荷施加不變，所得到的結(jié)果與前述方法相同.

在第1階模態(tài)基礎(chǔ)上，第2次計算時不改變隨機(jī)缺陷位置，使位移載荷在加載開始時跳過第1階的臨界位移載荷，而后緩慢增長，尋找第2階模態(tài)的臨界載荷，計算結(jié)果如圖5所示.

圖5 隨機(jī)缺陷法計算壓桿屈曲第2階模態(tài)Fig.5 The second mode of compression bar by stochastic defect method

通過以上算例可以看出，隨機(jī)缺陷法不僅可以計算模型屈曲的最低階模態(tài)，而且可以通過加載參數(shù)的改變來計算高階模態(tài)，計算結(jié)果與理論結(jié)果較為吻合.

3 彈性地基上梁的屈曲

3.1 嵌入于地基中的梁

將上文中的桿件嵌入于彈性地基中.地基使用Winkler地基模型進(jìn)行模擬，彈簧勁度系數(shù)β=20 MN/m，模型如圖6所示.

圖6 Winkler彈性地基梁Fig.6 Schematic of Winkler foundation model

使用通用有限元軟件進(jìn)行建模，采用梁單元進(jìn)行模擬，并將梁沿軸向劃分為100個單元，在除兩端外的每個網(wǎng)格節(jié)點處添加線性彈簧單元[9-10].

使用線性屈曲方法進(jìn)行計算，在梁的一端施加單位位移載荷，所得結(jié)果如圖7所示.計算特征值為7.655 86×10-2，所以臨界應(yīng)變εcr1=7.655 86×10-2.

3.2 擱置于地基表面的梁

在以上模型的基礎(chǔ)之上，加入非線性彈簧，即模擬梁未嵌入地基，而是放置在地基表面上的情況.此時，梁在向上彎曲時，與地基表面脫離；而向下彎曲時，則受到來自地基的壓力.模型相當(dāng)于兩端鉸支的壓桿穩(wěn)定問題.在建模過程中彈簧壓縮時的勁度系數(shù)β=20 MN/m，而拉伸時勁度系數(shù)趨近于0.

使用線性屈曲方法進(jìn)行計算，發(fā)現(xiàn)計算結(jié)果與使用線性彈簧時的結(jié)果相同.線性屈曲方法會選取彈簧初始剛度或者上一分析步結(jié)束時的剛度作為線性剛度進(jìn)行計算，而忽略后續(xù)彈簧參數(shù)的變化[11]，即不能計算非線性彈簧或者非線性材料等問題.而在此基礎(chǔ)之上通過影響因子的方式將線性屈曲方法的計算結(jié)果作為缺陷引入到模型中再進(jìn)行靜力計算，顯然也有問題，其結(jié)果的有效性有待驗證.同時可以認(rèn)為，對于材料非線性等問題，使用此方法仍然會失效.

使用隨機(jī)缺陷法進(jìn)行計算，同樣地隨機(jī)選取梁上節(jié)點引入10-5R0的橫向偏移，其計算結(jié)果如圖9與10所示.隨機(jī)缺陷法計算所得的臨界載荷Fcr=26.31 kN，與Euler公式的計算結(jié)果吻合， 說明隨機(jī)缺陷法對非線性彈簧問題同樣適用，相比于線性屈曲方法有較強(qiáng)的普適性.

考慮不同大小初始缺陷對最終臨界載荷Fcr計算結(jié)果的影響.對模型分別引入不同大小初始缺陷，所得計算結(jié)果如表1所示.其中:δ=R/R0，R為初始缺陷.

圖9 擱置地基表面梁受壓端軸向位移與整體支反力關(guān)系Fig.9 Relation of axial displacement of foundation surface beam compression end and integral bearing reaction

圖10 擱置地基表面梁中點處橫向位移與整體支反力關(guān)系Fig.10 Relation of transverse displacement at midpoint of foundation surface beam and integral bearing reaction

表1 桿初始缺陷與臨界載荷關(guān)系Tab.1 Relation table of stochastic defect and critical load

通過對比分析可以發(fā)現(xiàn)，當(dāng)δ=10-3或更小時，所得計算結(jié)果的精度較高.

4 橫向纖維搭橋下穿透型矩形脫層屈曲

設(shè)有長2a=200 mm，寬2b=100 mm的穿透型矩形脫層如圖11所示，受到兩側(cè)均布載荷Px的作用.脫層的厚度h=1 mm，在橫向纖維的影響下，受到均布載荷q(x)的作用抑制屈曲.材料的彈性模量E=2.3 GPa，泊松比μ=0.3.文獻(xiàn)[12]中也進(jìn)行了類似的計算，其所使用的方法是首先使用線性屈曲模塊，然后將1階屈曲模態(tài)位移乘以比例因子來進(jìn)行初始缺陷的引入，最終計算結(jié)果與理論結(jié)果沒有很好的對應(yīng)關(guān)系.

圖11 橫向纖維搭橋下矩形脫層Fig.11 Schematic of rectangular delamination bridged by fibers

根據(jù)文獻(xiàn)[13]所述理論進(jìn)行分析，

q(x)=βω(x),β>0

(10)

(11)

(12)

式中：q(x)為纖維約束力；β為纖維假設(shè)為彈簧時的勁度系數(shù)；ω(x)為屈曲撓度；D為抗彎剛度；a0為屈曲時特征長度.

根據(jù)表2中計算結(jié)果與理論值的對比可以發(fā)現(xiàn)，使用隨機(jī)缺陷法所得到的計算結(jié)果有較高的可信度.這種方法通過隨機(jī)地改變初始缺陷，可充分模擬實際構(gòu)件的缺陷形式，降低人為因素對計算結(jié)果的干擾，對工程應(yīng)用中分支點屈曲問題的模擬有較高的參考價值.

圖12 屈曲模態(tài)圖Fig.12 Schematic of buckling mode

5 結(jié)論

通過使用線性屈曲方法和隨機(jī)缺陷法分別對壓桿穩(wěn)定、彈性地基上梁的屈曲以及橫向纖維搭橋下穿透型矩形脫層屈曲問題進(jìn)行了有限元模擬.得到以下幾點結(jié)論：

(1) 隨機(jī)缺陷法將屈曲作為幾何非線性靜力計算問題進(jìn)行有限元模擬，改進(jìn)了線性屈曲計算中將屈曲作為特征值問題進(jìn)行求解的不足，與實際情況更為相符.

(2) 討論了隨機(jī)缺陷法的理論依據(jù)，并給出了通用有限元軟件使用初始隨機(jī)缺陷法的幾種思路.

(3) 隨機(jī)缺陷法具有更強(qiáng)的普適性，解決了線性屈曲分析無法計算非線性彈簧等問題.同時，隨機(jī)缺陷法同樣可以計算模型高階的屈曲模態(tài).

(4) 通過對比不同大小初始隨機(jī)缺陷的計算結(jié)果，發(fā)現(xiàn)初始隨機(jī)缺陷取構(gòu)件截面徑向尺寸千分之一或更小的值時，有限元計算結(jié)果與實際情況較為吻合.