FEA-08-梁与板壳问题_修改版.ppt

1. 薄板单元位移模式 设局部编号1、2、3、4，x 、y方向长度分别为2a、2b的矩形板单元如图所示。 x y z w3 ?y3 ?x3 Q1 My1 Mx1 1 2 4 3 1) 结点位移和结点力矩阵 图中还给出了各结点位移和结点力的示意图。 §8.7 弹性薄板矩形单元 2) 形函数的确定 薄板的形函数可以用广义坐标法，也可以用试凑法得到。由于单元自由度为12，因此可有12个广义坐标，位移模式可设为如下不完全四次多项式 利用12个结点位移条件，由广义坐标法可建立形函数，显然十分麻烦。龙驭球提出利用对称性较直接广义坐标法要容易一些，也还有很大工作量。 为此介绍试凑法,首先引入自然坐标?=x/a，?=y/b。 x y z w3 ?y3 ?x3 Q1 My1 Mx1 1 2 4 3 3) 试凑形函数N1 由形函数性质，对N1有: N1(1)=1；N1(j)=0，j=2,3,4 N1对x，y的偏导数在结点处均为零。 利用所有点N1的导数为零条件，P.125 经式(c)~(l)的推导，可得 考虑到挠度是非完全四此式，为使自动满足它点为零N1(j)=0 ，可设 再由本点处位移的条件，可得d=-1/8，由此 x y z w3 ?y3 ?x3 Q1 My1 Mx1 1 2 4 3 4) 其他形函数Nj、Nix 、Niy 记 ?0 = ??i ；? 0=? ?i 仿N1可得: 对于转角?xi相关的形函数，同样思路推导可得 对于转角?yi相关的形函数，可推导得 x y z w3 ?y3 ?x3 Q1 My1 Mx1 1 2 4 3 5) 薄板的挠度场 有了每一结点的形函数，记 则薄板的挠度场可由结点位移表示为 6) 单元间位移的协调性 可以证明，上述w在边线上任意一点的挠度和转角都是三次多项式。 x y z w3 ?y3 ?x3 Q1 My1 Mx1 1 2 4 3 因此，边线的挠度和转角可由两端点的挠度和沿边线导数对应的转角唯一地确定。 但是，边界法向转角只有两端两个法向转角位移条件，当然无法唯一地确定，所以相邻单元法向转角位移不协调。 由此可见，由形函数所建立的挠度场 是非完全协调的。（教材上有更严密的数学证明） x y z w3 ?y3 ?x3 Q1 My1 Mx1 1 2 4 3 7) 非完全协调元的收敛性 对于薄板等位移场非完全协调的位移模式,如何才能保证收敛呢? Irons 给出了分片检验准则： 用待检验的单元组成一小片，在无荷载、单元结点位移满足“常应变”状态位移条件时，如果各结点能够保持平衡且获得“常应力”受力状态，则这种位移模式对应的单元在如此网格下一定收敛。 教材给出了具体检验的方法、步骤，请自学。 四、薄板的单元列式 1) 总势能用结点位移表示 将位移模式代入可得形变矩阵为 x y z w3 ?y3 ?x3 Q1 My1 Mx1 1 2 4 3 2) 薄板单元刚度方程 式中 由总势能的一阶变分为零可得 x y z 壳体一般其中面都是曲面，对于柱面壳体，可以用一系列小的矩形平板来逼近它。此外还有一些折板结构，它们都可以用平板壳元来分析。 柱壳 折板结构 1. 平板壳体单元刚度方程 壳体结构象拱一样，中面既受面内力作用产生面内变形外，还将象板一样产生弯曲。但由于弯曲中面无变形，面内变形又不产生弯曲，两者互不藕联。因此，可以象杆件单元一样，用平面应力和平板弯曲组合来得到壳体单元刚度方程。 §8.8 矩形平板壳体单元 2. 坐标转化问题 平面问题（除等参元外）和平板弯曲问题单元局部坐标和整体坐标一致，因此没有坐标转换问题。如图示意，各平板壳元局部坐标可能不同，所以集装前要做坐标转换。 x z y 3 2 取母线为局部x 如图由1、4点的坐标，可以确定局部坐标y的单位向量，局部坐标x单位向量已知，由这两单位向量的向量积就可获得局部z坐标单位向量。 按上述思路，需要建立各单位向量的方向余弦，由此即可获得坐标转换矩阵[T]。 1 4 壳体边界条件的提法 对称边 简支边 自由边 1/4圆柱壳 x y z 固定边 简支点 ANSYS壳单元shell181 Shell 181适用于薄到中等厚度的壳结构。基本原则是每块面板的主尺寸不低于其厚度的10倍。壳单元有四个节点，单元每个节点有六个自由度，分别为沿节点X，Y，Z方向的平动及绕节点X，Y，Z轴的转动。退化的三角形选项用于网格生成的过渡单元。Shell 181单元具有应力刚化及大变形功能。该单元有强大的非线性功能，并有截面数据定义，分析，可