基于有限元分析板時(shí)厚跨比的選擇

板和殼在生產(chǎn)和建設(shè)中得到廣泛使用。隨著科學(xué)技術(shù)的發(fā)展，除傳統(tǒng)的薄膜（外殼）外，多個(gè)特殊的板殼結(jié)構(gòu)也被用于。例如，氣冷核反應(yīng)堆的預(yù)測(cè)混凝土壓力容器、各種中間板（外殼）和環(huán)保板（外殼）。在工程中，廣泛應(yīng)用了kirchhoff-hele-hele-hele-hele-hele-hele-hele-hele的薄板理論和薄板的薄板理論。在實(shí)際工程中，由于結(jié)構(gòu)幾何形狀、邊界條件和負(fù)荷的復(fù)雜性，使用分析方法幾乎是不可能的。其中大多數(shù)使用各種數(shù)值方法進(jìn)行分析。然而，由于有限的元法只需劃分區(qū)域中的單元，所需分解的精度和計(jì)算成本較大，因此大多數(shù)研究人員開始研究邊界積分方程（biem）。在邊界分析方法的介紹下，邊界分析單元被引入不同的區(qū)域，其中f.vadweec.a.brebas，j.c.f.teles在邊界分析理論中研究了厚板的邊界方法。本文采用Reissner中厚板理論邊界元程序、Mindlin中厚板理論有限元程序和三維彈性力學(xué)有限元程序計(jì)算各種厚跨比下四邊固支和四邊簡(jiǎn)支矩形板,并進(jìn)行對(duì)比,來說明三維理論、中厚板理論和薄板理論分析板時(shí)厚跨比的適用范圍以及精確程度,為采用何種理論來分析板提供理論依據(jù)和參考;同時(shí)也驗(yàn)證了采用縮減積分的Mindlin中厚板有限元法消除了剪切鎖死現(xiàn)象.1厚板控制方程Reissner中厚板理論是考慮了橫向剪切變形和擠壓變形的理論.考慮笛卡爾坐標(biāo)系(x1,x2,x3),其中x1,x2軸在板中面內(nèi),x3則垂直于板中面,如圖1所示.在板的上下表面切應(yīng)力為零而正應(yīng)力為分布載荷集度,即x3=±h/2,σ33=±q/2.各應(yīng)力在厚度方向上的變化關(guān)系如下:σij=12x3h3Μij?σi3=32h[1-(2x3h)2]Qi?σ33=14(2x3h)[3-(2x3h)2]q.(1)σij=12x3h3Mij?σi3=32h[1?(2x3h)2]Qi?σ33=14(2x3h)[3?(2x3h)2]q.(1)式中:i,j=1,2;Mij為單位長(zhǎng)度彎矩;Qi為剪力.為了求解Reissner中厚板控制方程的基本解,在此假定體力在厚度方向上也按應(yīng)力一樣變化,即Ρi=Bi(x1,x2)12x3h3?Ρ3=B3(x1,x2)32h[1-(2x3h)2].(2)Pi=Bi(x1,x2)12x3h3?P3=B3(x1,x2)32h[1?(2x3h)2].(2)由以上假定,Reissner中厚板的控制方程變?yōu)?{Μij,j-Qi=0,Qi,i+q=0.(3)其中的彎矩、扭矩和剪力用轉(zhuǎn)角和位移表示為:Μij=D(1-ν)2(βi,j+βj,i+2ν1-νβk,kδij)+ν1-νqλ2δij,Qi=D(1-ν2)λ2(βi+w,i).(4)式中:D=Eh3/[12(1-ν2)]為板的彎曲剛度;ν為泊松比;λ=√10/h是厚度參數(shù).式(4)中轉(zhuǎn)角和位移為沿板厚的平均值,即:βi=∫h/2-h/212x3hvidx3,w=∫h/2-h/232h[1-(2x3h)2]v3dx3.(5)式中:vi和v3為坐標(biāo)軸方向的位移分量.為了簡(jiǎn)便,將βi(i=1,2)和w用ui(i=1,2,3)表示,則對(duì)于區(qū)域Ω的邊界Γ,位移和廣義面力的邊界條件可以表示為:ui=ˉui,在Su上;ti=ˉti,在St上.其中S=Su∪St,ti為廣義面力.假定邊界上某點(diǎn)外法線方向余弦為ni=cos(n,xi),則廣義面力為ti=Mijnj,t3=Qini.(i,j=1,2)(6)2厚板邊界積分基本解可由式(3)中體力項(xiàng)定義為一塊無限大中厚板作用單位橫向集中力或單位彎矩求得,具體推導(dǎo)見文獻(xiàn).僅列出關(guān)于位移和面力的基本解:u*ij=12πD{21-ν[B(z)δij-A(z)r,ir,j]-12δij(lnz-12)-12r,ir,j}?u*i3=-u*3i=12πD12rr,i(lnz-12)?u*33=-12πD1λ2[21-νlnz-14z2(lnz-1)].(7)t*ij=-14πr[(4A(z)+2zΚ1(z)+1-ν)(δijr,nnj)+(4A(z)+1-ν)r,jni-2(8A(z)+2zΚ1(z)+1-ν)r,ir,jr,n]?t*i3=λ22π[B(z)ni-A(z)r,ir,n]?t*3i=-1-ν8π[(21+ν1-νlnz-1)ni+2r,ir,n]?t*33=-12πrr,n.(i,j=1,2)(8)上兩式中:A(z)=Κ0(z)+2z[Κ1(z)-1z];B(z)=Κ0(z)+1z[Κ1(z)-1z];z=λr;r=[(x1-ξ1)2+(x2-ξ2)2]1/2;r,i=(xi-ξi)/r;r,n=r,ini.其中K0(z)和K1(z)為修正貝塞爾函數(shù).使用Betti功互等定理、格林函數(shù)或加權(quán)殘值法,可以得到Reissner中厚板的邊界積分方程:cijuj+∫Γt*ijujdΓ=∫Γu*ijtjdΓ+∫Ω[u*i3-ν(1-ν)λ2u*ik,k]qdΩ?i,j=1,2,3；k=1,2.(9)式中:cij是一個(gè)與邊界形狀有關(guān)的系數(shù),對(duì)于光滑邊界cij=δij/2,對(duì)于有角點(diǎn)和尖點(diǎn)的邊界法線不連續(xù)的邊界,可以解析計(jì)算cij,但更方便的是間接計(jì)算cij的值.當(dāng)整個(gè)邊界無外力和載荷作用時(shí),式(9)變?yōu)?cijuj+∫Γt*ijujdΓ=0.(10)上式的非零解為板的3種獨(dú)立剛體位移的任意組合,把這些位移代入式(10)有cia(ξ)=-∫Γ[t*ia(ξ,x)+(ξa-xa)t*i3(ξ,x)]dΓ,ci3(ξ)=-∫Γt*i3(ξ,x)dΓ,(11)i=1,2,3;α=1,2.當(dāng)均布荷載情況下,可以應(yīng)用散度定理把式(9)中的域積分簡(jiǎn)化為邊界積分.數(shù)值方法求解式(9)后,未知邊界位移和邊界力可以全部求出.令cij=δij可以求出內(nèi)部點(diǎn)的位移,而內(nèi)部點(diǎn)的內(nèi)力解,需要對(duì)式(9)關(guān)于源點(diǎn)ξ求微分,有如下形式:Ρij=ν(1-ν)λ2qδij+∫Γ[Dijktk-Sijkuk+qwij]dΓ.(12)式中:i,j≠3.P的分量為Pij=Mij,P3j=Qj.其中wij,Dijk和Sijk見文獻(xiàn)中附錄Ⅰ.3標(biāo)準(zhǔn)高斯積分法把所考慮的板的邊界劃分為M個(gè)單元,以每個(gè)節(jié)點(diǎn)為源點(diǎn)并在M個(gè)單元上積分,離散式(9),并將其中域積分轉(zhuǎn)換為邊界積分,則可得線性方程組HU=GT+D.(13)式中:H,G為影響矩陣;D為對(duì)應(yīng)于載荷q的邊界積分列向量.代入已知的邊界條件后,將其中的已知項(xiàng)移到方程等號(hào)右邊與D向量合并,未知項(xiàng)移到等號(hào)左邊,得到標(biāo)準(zhǔn)代數(shù)方程組AX=B.(14)采用標(biāo)準(zhǔn)高斯消去法求解上述方程組,即可求得未知量.之后可以利用式(9)計(jì)算內(nèi)部任意一點(diǎn)的位移,利用式(12)求內(nèi)部點(diǎn)的內(nèi)力.由于式中的奇異性,在計(jì)算奇異積分時(shí),采用了三次多項(xiàng)式非線性變換.此種變換可以在奇異點(diǎn)附近集中較多的高斯點(diǎn),因此對(duì)奇異積分采用標(biāo)準(zhǔn)高斯積分能達(dá)到相當(dāng)高的精度.對(duì)于存在角點(diǎn)和尖點(diǎn)等邊界法線不連續(xù)的情況,采用雙節(jié)點(diǎn)或部分不連續(xù)單元.雙節(jié)點(diǎn)單元中2個(gè)節(jié)點(diǎn)有不能在同一方向上給定位移的限制,而不連續(xù)單元沒有這種限制.根據(jù)前述離散化的公式,編制了Reissner中厚板邊界元的FORTRAN程序.4固支矩形板.采用了2個(gè)具有代表性的算例來研究薄板和中厚板理論適用范圍和精確程度.算例為:1)受均布荷載作用的四邊固支矩形板,如圖2所示;2)受均布荷載作用的四邊簡(jiǎn)支矩形板.兩算例的長(zhǎng)、寬取a=b=1,橫向荷載集度q=1,板的彎曲剛度D=1,泊松比ν=0.3.4.1接觸模型的建立利用本文編制的Reissner中厚板邊界元程序計(jì)算了在不同厚跨比(h/a=0.001,0.005,0.01,0.05,0.1,0.2,0.5)下的板中心點(diǎn)的撓度w,板中心點(diǎn)下表面應(yīng)力σx和板固支邊中點(diǎn)的彎矩Mx.并與利用有限元軟件ABAQUS中對(duì)應(yīng)的三維有限元解和板殼有限元解,還有部分薄板解析解進(jìn)行對(duì)比.利用ABAQUS的三維8節(jié)點(diǎn)單元、二維4節(jié)點(diǎn)板殼單元和2節(jié)點(diǎn)邊界單元進(jìn)行分析計(jì)算,其中以三維8節(jié)點(diǎn)有限元分析的結(jié)果作為基準(zhǔn),以說明二維4節(jié)點(diǎn)板殼單元和2節(jié)點(diǎn)邊界單元計(jì)算結(jié)果的準(zhǔn)確性.在三維有限元、板殼有限元和中厚板邊界元計(jì)算中,取中面中心點(diǎn)的撓度和中心點(diǎn)板的下表面點(diǎn)的應(yīng)力σx作為比較對(duì)象.因?yàn)槿S單元只能得到節(jié)點(diǎn)的位移和應(yīng)力,要想得到截面上的內(nèi)力時(shí)需要進(jìn)行合成計(jì)算,為此取中心點(diǎn)板的下表面點(diǎn)的應(yīng)力σx作為比較對(duì)象.板殼有限元分析結(jié)果既可輸出內(nèi)力,也可輸出應(yīng)力,但由于邊界元只能輸出內(nèi)力,所以需要利用彎矩Mx與應(yīng)力σx的關(guān)系求得:σx=12Μxh3z.(15)式中:z是應(yīng)力計(jì)算點(diǎn)到中面的距離.為了更好、更精確地對(duì)比各方法的計(jì)算結(jié)果,對(duì)Reissner中厚板邊界元模型每邊劃分16個(gè)線性單元,整個(gè)板共64個(gè)單元;對(duì)板殼有限元模型同樣每邊劃分16個(gè)線性單元,共16×16個(gè)單元;三維有限元模型則在板面內(nèi)每邊劃分16個(gè)線性單元,厚度方向上網(wǎng)格的劃分是依照三維單元盡量要成正六面體的準(zhǔn)則來進(jìn)行的,這樣能保證作為基準(zhǔn)的三維實(shí)體有限元結(jié)果的精度.圖3畫出了各種厚跨比下板中點(diǎn)的撓度.從圖3中可以看出在所有厚跨比情況下,3種方法的計(jì)算結(jié)果十分吻合,說明了考慮剪切變形的Reissner中厚板理論是對(duì)三維彈性力學(xué)的較精確的簡(jiǎn)化.從圖中還可以看出,當(dāng)厚跨比大于0.2時(shí),與三維單元的結(jié)果誤差明顯增大.從圖4可看到,在厚跨比較小(h/a<0.2)時(shí),邊界元、板殼有限元和三維有限元計(jì)算出的σx十分接近,但是在h/a>0.2后,誤差越來越大.但是邊界元解精度比板殼有限元的解精度高,證明了邊界元法在理論上比有限元法應(yīng)力(內(nèi)力)的精度高.綜合各種厚跨比撓度和應(yīng)力的誤差,中厚板與厚板的厚跨比分界線一般取0.2為宜,厚跨比超過了0.2,無論是撓度還是應(yīng)力的誤差都會(huì)急劇增加.同時(shí),從圖3和圖4還可以看出:利用縮減積分技術(shù)的Mindlin板殼有限元,在計(jì)算非常薄的板時(shí)大大緩解了剪切鎖死現(xiàn)象.圖5還給出了固支邊中點(diǎn)處的彎矩Mx.在計(jì)算彎矩時(shí),三維實(shí)體單元選取20節(jié)點(diǎn)二次單元,因?yàn)槎螁卧奈灰剖嵌畏植嫉?但是應(yīng)力是線性分布的,所以二次單元積分出的彎矩應(yīng)該更精確.從圖5可以看出邊界元解比板殼有限元的解精確得多,而且在薄板階段和解析解十分靠近,充分說明了中厚板邊界元解對(duì)薄板有很好的精度.邊界元解和有限元三維實(shí)體單元解有差異,應(yīng)當(dāng)是來源于彎矩計(jì)算時(shí)積分的誤差,因?yàn)槿S單元只能輸出應(yīng)力,想獲得彎矩還需對(duì)應(yīng)力進(jìn)行積分,而節(jié)點(diǎn)并不是在高斯積分點(diǎn)上,只能使用梯形法近似積分,致使精度較差.4.2單元?jiǎng)澐旨罢鍎澐秩鐖D6所示四邊簡(jiǎn)支矩形板,所有數(shù)據(jù)同例4.1.板殼有限元仍然采用雙線性板殼單元,而中厚板邊界元采用2節(jié)點(diǎn)線性單元.另外單元?jiǎng)澐譁?zhǔn)則也跟算例4.1相同,邊界元在每邊上劃分16個(gè)單元,整板共64個(gè)單元;板殼有限元也在每邊劃分16個(gè)單元.從圖7中可以看出,3種單元的四邊簡(jiǎn)支板中心的撓度的解十分接近.在厚板階段(h/a>0.1),圖中清楚地表明了本文所用中厚板邊界元法比板殼有限元法更接近三維有限元的解,說明了邊界元法精度比有限元高.從圖8中可以看出,無論是邊界元還是板殼有限元的解都十分接近三維有限元的解.從圖中可以看到,在薄板階段(h/a<0.1),邊界元解和薄板理論十分接近,和三維有限元解、板殼有限元解有些許誤差,因?yàn)榘鍤び邢拊獾膬?nèi)力精度要低于邊界元解.5板殼有限元法本文將Reissner中厚板理論運(yùn)用到邊界元法中,并通過數(shù)值方法將其編制成FORTRAN程序,通過幾個(gè)典型算例的計(jì)算,證明了邊界元法分析中厚板的程序的正確性.另外,通過與三維有限元、板殼有限元的對(duì)比,證明了Reissner理論是對(duì)三維彈性力學(xué)