构力学II第一套王焕定7弹性薄板.pptVIP

五、弹性薄板和平板壳单元分析初步及程序 五、弹性薄板和平板壳单元分析初步及程序 弹性薄板基本知识 弹性薄板矩形（R12）单元 薄板分析程序的使用 矩形平板壳体单元 平板壳体程序的使用 计算结构力学基础的结束语 5.1 弹性薄板基本知识 5.1.1 弹性薄板基本概念 x y z u v w 中面 所谓薄板是指板厚h比板最小尺寸b在如下范围的平板 板面位移如图所示。 当w小于板厚h时，有克希霍夫(G.kirchhoff)假定成立： a) 板中面是中性面，没有变形。 b) 中面法线变形后仍为挠曲面法线，长度不变。 c) 忽略应力?z和应变? z 。 平分厚度的平面称中面。 5.1 弹性薄板基本知识 由克希霍夫假定c) 忽略应变? z可推得w与z无关，由b)可知? zx和? yz等于零，再加上中面无变形，最终可得 由此结果可得 x向曲率 y向曲率 扭率 他们完全确定板的变形，因此称他们组成的矩阵为形变矩阵，记作[?]，也即 位移只与挠度w有关 5.1 弹性薄板基本知识 由此可得薄板应变矩阵为[?]=z[?]。 5.1.2 薄板内力和总势能 1) 设平面应力弹性矩阵为[D]’，则薄板应力矩阵为 [?]=-z[D]’ [?]。 x y z 扭矩 弯矩 2) 薄板内力 微元体如图所示。 由图可得 5.1 弹性薄板基本知识 由此可得薄板单位长度内力为Mx、My、Mxy= Myx (dx=dy=1)，依此顺序排列的列阵称内力矩阵，记作[M]。 将应力应变关系代入并对z进行积分，可得 [M]=[D][?] 式中 [D]=（h3/12）[D]’ 称作薄板的弹性矩阵。 3) 薄板的应变能 5.1 弹性薄板基本知识 4) 薄板的总势能 设薄板受z方向分布荷载q(x)作用，则线弹性薄板的总势能为 上式就是下面作有限元分析的理论依据。 能写出各向同性弹性体的[D]矩阵吗? 能写出正交各向异性弹性体的[D]矩阵吗? 5.2 弹性薄板矩形（R12）单元 5.2.1 薄板单元位移模式 设局部编号1、2、3、4，x 、y方向长度分别为2a、2b的矩形板单元如图所示。 x y z w3 ?y3 ?x3 Q1 My1 Mx1 1 2 4 3 1) 结点位移和结点力矩阵 图中还给出了各结点位移和结点力的示意图。 5.2 弹性薄板矩形（R12）单元 2) 形函数的确定 薄板的形函数可以用广义坐标法，也可以用试凑法得到。由于单元自由度为12，因此可有12个广义坐标，位移模式可设为如下不完全四次多项式 x y z w3 ?y3 ?x3 Q1 My1 Mx1 1 2 4 3 利用12个结点位移条件，由广义坐标法可建立形函数，显然十分麻烦。龙驭球提出利用对称性较直接广义坐标法要容易一些，也还有很大工作量。 为此介绍试凑法,首先引入自然坐标?=x/a，?=y/b。 ？ 5.2 弹性薄板矩形（R12）单元 2-1) 试凑形函数N1 由形函数性质，对N1有: N1(1)=1；N1(j)=0，j=2,3,4 N1对x，y的偏导数在结点处均为零。 x y z w3 ?y3 ?x3 Q1 My1 Mx1 1 2 4 3 利用所有点N1的导数为零条件，P.125 经式(c)~(l)的推导，可得 考虑到挠度是非完全四此式，为使自动满足它点为零N1(j)=0 ，可设 再由本点处位移的条件，可得d=-1/8，由此 5.2 弹性薄板矩形（R12）单元 2-2) 其他形函数Nj、Nix 、Niy 记 ?0 = ??i ；? 0=? ?i 仿N1可得: x y z w3 ?y3 ?x3 Q1 My1 Mx1 1 2 4 3 对于转角?xi相关的形函数，同样思路推导可得 对于转角?yi相关的形函数，可推导得 5.2 弹性薄板矩形（R12）单元 3) 薄板的挠度场 有了每一结点的形函数，记 x y z w3 ?y3 ?x3 Q1 My1 Mx1 1 2 4 3 则薄板的挠度场可由结点位移表示为 4) 单元间位移的协调性 可以证明，上述w在边线上任意一点的挠度和转角都是三次多项式。 5.2 弹性薄板矩形（R12）单元 因此，边线的挠度和转角可由两端点的挠度和沿边线导数对应的转角唯一地确定。 但是，边界法向转角只有两端两个法向转角位移条件，当然无法唯一地确定，所以相邻单元法向转角位移不协调。 x y z w3 ?y3 ?x3 Q1 My1 Mx1 1 2