下载文档原格式
薄板弯曲问题的有限元法一、 薄板弯曲问题的基本方程什么是薄板？薄板就是指厚度t 远小于其长度、宽度的板。
1. 三个基本假设（克希霍夫假设）： (1) 法线假设，εz =0,γyz =γzx =0 (2) 正应力假设，σz <<σx ,σy ,τxy (3) 小挠度假设，w<t/4根据假设，可以得到位移分量()()()()()(),,,,,,,,,,, x y z u x y z z x x y z v x y z z y x y z x y ωωωω∂⎧=-⎪∂⎪∂⎪=-⎨∂⎪⎪=⎪⎩式4-1图 1 薄板弯曲后某点B 的位移2. 应变分量{}222222x y z x z y x y ωεωεεεω⎧⎫∂-⎪⎪∂⎪⎪⎧⎫⎪⎪∂⎪⎪⎪⎪==-⎨⎬⎨⎬∂⎪⎪⎪⎪⎩⎭⎪⎪∂-⎪⎪∂∂⎪⎪⎩⎭式4-23. 曲率{}222222x y z x y x y ωχωχχχω⎧⎫∂-⎪⎪∂⎪⎪⎧⎫⎪⎪∂⎪⎪⎪⎪==-⎨⎬⎨⎬∂⎪⎪⎪⎪⎩⎭⎪⎪∂-⎪⎪∂∂⎪⎪⎩⎭式4-3 22=x x ωχ∂-∂——薄板弹性曲面在x 方向的曲率22=y yωχ∂-∂——薄板弹性曲面在y 方向的曲率2=z x yωχ∂-∂∂——薄板弹性曲面在x 方向和y 方向的扭率4. 应力分量与应变分量间的关系：{}[]{}2222222222221 11D Ez xy Ez x y Ez x y σεωωμμωωμμωμ=⎧⎫⎛⎫∂∂-+⎪⎪ ⎪-∂∂⎝⎭⎪⎪⎪⎪⎛⎫∂∂⎪⎪=-+⎨⎬ ⎪-∂∂⎝⎭⎪⎪⎪⎪∂⎪⎪--∂∂⎪⎪⎩⎭式4-4 5. 线力矩{}()2222222101012110022x y z x M Et M M y M x y ωμωμμμω⎧⎫∂-⎪⎪⎡⎤∂⎪⎪⎢⎥⎧⎫⎪⎪⎢⎥∂⎪⎪⎪⎪==-⎨⎬⎨⎬⎢⎥∂-⎪⎪⎪⎪⎢⎥-⎩⎭⎪⎪⎢⎥∂⎣⎦-⎪⎪∂∂⎪⎪⎩⎭式4-5a广义应力与广义应变之间的关系式{}[]{}D M χ= 式4-5b式中：[D]—薄板弯曲问题的弹性矩阵6. 薄板弯曲问题的基本方程（双调和方程）()32222222121Et p xx y y ωωωμ⎛⎫∂∂∂++= ⎪∂∂∂∂-⎝⎭ 式4-6()32121Et μ-——薄板弯曲刚度 二、 矩形薄板单元分析 1、矩形薄板单元图 2 矩形薄板单元2、位移函数22123456322333789101112 a a x a y a x a xy a y a x a x y a xy a y a x y a xy ω=+++++++++++ 式4-73、形状函数[]{}k i i xi xi yi yi j j xj xj yj yj k kxk xk yk y l l xl xl yl yl N N N N N N N N N N N N N q ωωθθωθθωθθωθθ=+++++++++++= 式4-8式中：i,j,k,l ——节点号N i ，N xi ，N yi ，……，N yl ——形状函数()()()()()()()()()()2211128N 111 ,,,8111 8y i i i i i xi i i iyi i i i b N i i j h l N a x a b ξξηηξξηηξηηξξηηηξξξηηξξη⎧⎫++++--⎪⎪⎧⎫⎪⎪⎪⎪⎪⎪=-++-=⎨⎬⎨⎬⎪⎪⎪⎪⎩⎭⎪⎪++-⎪⎪⎩⎭==， 式4-94、单元刚阵[][][][]S K TB D B dxdy =⎰ 式4-10式中：[]22222222222222222222 2222yi yl i xiyi yl ixi yi yl i xi N N N N x x x x N N NN B y y y y N N N N x yx yx yx y ⎡⎤∂∂∂∂⎢⎥∂∂∂∂⎢⎥⎢⎥∂∂∂∂=⎢⎥∂∂∂∂⎢⎥⎢⎥∂∂∂∂⎢⎥⎢⎥∂∂∂∂∂∂∂∂⎣⎦式4-11 5、节点力与节点位移的关系式{}[]{}F K q = 式4-12三、 三角形薄板单元分析1、三角形薄板单元当薄板具有斜交边界或曲线边界时，可采用三角形单元较好地反映边界形状。
