您當前的位置:首頁 > 農業

(大四)板殼理論——薄板的小撓度彎曲

作者:由 Ready-Player 發表于 農業時間:2018-06-24

接著上一次寫的內容:

QuYln:板殼理論——Kirchhoff–Love Plate Theory

本次寫三個部分:

(1)控制方程的推導思路;

(2)邊界條件的推導思路;

(3)簡單小撓度彎曲的算例;

一、控制方程的推導思路:

1。首先給出小撓度彎曲的控制方程:

D\nabla^4\omega=q

,這個控制方程其實是由“板面的邊界條件”得來的,即:

(\sigma_z)_{z=-\frac{\delta}{2}}=-q

,下表面的邊界條件是恆等式零,故沒有意義。

2。接下來,說幾個注意事項:

(1)

\sigma_z

是非常小的應力,“擠壓應力”,就像材料力學建立“伯努利-尤拉梁”一樣,直接忽略掉了擠壓應力,板殼理論作為“應用彈性力學”,在此處

忽略掉了擠壓應力的影響

,得到假設

\varepsilon_z=0

,即不考慮板厚度方向上的應變,以中面該點的位移來代表板厚度方向上各點的位移,但是

\sigma_z

作為影響平衡的關鍵,不能忽略不計。

(2)

\tau_{yz},\tau_{xz}

也是非常小的應力,也忽略其影響,但本身由於起到維持平衡的作用,故不能忽略。由此得到的假設是:

\gamma_{yz},\gamma_{xz}

為零,體會其物理含義,切應變會影響法線的角度變化,角度不統一,得到的法向也不同,忽略了這個現象,得到

直法線假定

3。方程的推導思路:

前面都說了,控制方程是擠壓應力的邊界條件,那麼現在的問題就是

求出擠壓應力

(1)由於基爾霍夫假定中對且應變的假設,使得按位移求解的三個未知量

u,v,w

的獨立性消失,僅保留一個獨立的未知量

w

\gamma_{zx}=\frac{\partial w}{\partial x}+\frac{\partial u}{\partial z}=0

,因此:

\frac{\partial u}{\partial z}=-\frac{\partial w}{\partial x}

,由此可得

u=-\int\frac{\partial w}{\partial x}dz+f_{1}(x,y)

u=-\frac{\partial w}{\partial x}z+f_{1}(x,y)

,右由於中面假定沒有水平位移,所以:

u=-\frac{\partial w}{\partial x}z

,另外,

v=-\frac{\partial w}{\partial y}z

(2)已知,獨立的未知量為

w(x,y,z)

,再加上假定,所以唯一的位移未知量為

w(x,y)

有了位移,可以透過

幾何關係

,得到相應的應變:

\varepsilon_x=\frac{\partial u}{\partial x} =-\frac{\partial^ 2u}{\partial x^2}z =\chi_x z

,曲率一定,可見應變呈線性分佈;

\varepsilon_y=\frac{\partial v}{\partial y} =-\frac{\partial^ 2v}{\partial y^2}z =\chi_y z

,曲率一定,可見應變呈線性分佈;

\gamma_{xy}=\frac{\partial v}{\partial x}+\frac{\partial u}{\partial y} =-2\frac{\partial^ 2w}{\partial x\partial y}z =2\chi_{xy}z

,扭率一定,可見應變呈線性分佈;

(3)得到了由縱向位移表示的應變後,再由

本構關係

得到應力:

\sigma_x=\frac{E}{1-\mu^2}(\varepsilon_x+\mu\varepsilon_y)

 \sigma_x=-\frac{E}{1-\mu^2}(\frac{\partial^ 2u}{\partial x^2}+\mu\frac{\partial^ 2v}{\partial y^2}) z

,可見,x方向正應力也呈線性分佈;

由此,可以得到其餘的兩個應力:

 \sigma_y=-\frac{E}{1-\mu^2}(\frac{\partial^ 2v}{\partial y^2}+\mu\frac{\partial^ 2u}{\partial x^2}) z

,y方向正應力呈線性分佈;

 \tau_{xy}=-\frac{E}{1+\mu}\frac{\partial^ 2w}{\partial x\partial y}z

,切應力呈線性分佈;

(4)最終要轉向為擠壓應力,於是可以透過

平衡方程

來聯絡各個應力:

此處是把體力和下表面的面力等效成了上表面的面力,於是平衡方程中不包含體力項。

另外,要注意的是,只考慮橫向荷載,不考慮平行於中面的載荷,平行於中面的載荷可以由疊加原理重新新增一個“平面應力問題”。

透過平衡方程,先確定與z有關的兩個切應力,由於存在積分,可以透過邊界條件來確定待定函式,再透過與z有關的兩個切應力來確定z方向的擠壓應力,同樣的,透過邊界條件來確定待定函式,有了擠壓應力,即得到了控制方程。

得出的結論是:

\tau_{zx}=\frac{E}{2(1-\mu^2)}( z^2-\frac{\delta^2}{4})\frac{\partial}{\partial x}\nabla^2w

\tau_{zy}=\frac{E}{2(1-\mu^2)}( z^2-\frac{\delta^2}{4})\frac{\partial}{\partial y}\nabla^2w

兩個切應力沿著厚度方向呈

二次分佈

\sigma_{z}=-\frac{E\delta^3}{6(1-\mu^2)} (\frac{1}{2}-\frac{z}{\delta})^2(1+\frac{z}{\delta}) \nabla^4w

擠壓應力沿著厚度方向呈

三次分佈

標簽: 應力  方程  擠壓  應變  邊界條件