有限元第七章弹性板壳有限元分析.pptVIP

第7章弹性板壳有限元分析初步 板壳结构式土木工程中广为采用的一种结构形式，但是只有少量几何形状、荷载及边界条件较简单的情况才能获得解析结果。因此，用有限元方法进行板壳结构受力分析研究，自有限元产生开始，直至今天仍在发展。本书仅就一种矩形单元简单介绍有限元板壳分析的初步知识，为读者进一步深入学习打下基础。 7-2-2非完全协调元的收敛性准则 对拟协调元艾恩斯指出： 对于不能完全满足协调元准则的单元，若对任意下片的几个单元集合施加与任一常应变状态相应的结点位移，如果能在外部荷载作用下满足结点平衡条件并且获得常应力状态，则称此单元能通过小片检验。而且实践证明，小片检验所要求的条件是保证非协调任意单元收敛性的充分条件。 对于某种薄板单元位移模式，进行小片检验的具体步骤如下： （1）取某一单元小片，并在小片的边界上给出对应 于完全二次多项式的边界条件。 （2）按此位移模式进行（在无荷载作用下的）单元、 整体分析，并在上述位移边界条件下求解。 （3）若所求得的结点位移构造的小片上的挠度为一 完全二次多项式，则单元的位移模式通过小片检验。 当薄板程序不能解已知边界支撑位移时，也可按如下步骤进行小片检验： （1）取某一单元小片，对小片的每一结点给以对应于完全二次多项式的结点位移。 （2）每一单元按 求单元结点力，式中 为对应所考察位移模式的单元刚度矩阵。 （3）检验小片内部结点处是否均满足 若条件均成立，则单元位移模式能通过小片检验。 例如图7-4所示单元小片，取 则 按上述式子给结点以位移后，可求得 7-2-3 单元分析 因为单元挠度为 * 7-1 弹性薄板基本理论 所谓薄板（thin plate）系指板厚 与板面最小尺寸b的 比值在下列范围内的平板( ~ ) ( ~ ) 1.板的中面在自身平面内没有变形。也即在弯曲时中面是中性面。 2.弯曲前板内垂直于中面的直线段在弯曲后仍保持为直线，并垂直于中 性曲面，这线段的长度不变。 平分板厚度 的平面为中面（中平面）。在板的挠度 小于厚度 时有克希霍夫（G.Kirchoff）假设： 3.忽略应力 及应变 。 根据以上假设可得 由 可推知 由 可推知 由 可推知 图7-1 弹性薄板示意 因为假设中面无变形，也即 时 所以最终可得 综上可知薄板分析主要是求解挠度 ，它是x,y的函数。由此出发则薄板的其余 为应变分量为 在小变形情况下 、 、 分别为弹性曲面在x方向和y方 向的曲率以及x、y方向的扭率，它们完全确定了板内各点的应变，故称如下 矩阵的形变矩阵（deformation matrix）,记作 （7-1） 由此可见薄板的应变矩阵(strain matrix)为 与此应变矩阵相应的各向同性体的弹性矩阵（elasticity matrix）为 薄板应力矩阵（stress matrix）为 根据图7-2示意可得内力为 （7-2） （7-3） 图7-2 内力与应力关系示意 考虑到式（7-3）则内力与形 变间的关系为 式中：D为薄板的弹性矩阵。 对各向同性体 （7-4） 从式（7-2式）和式（7-3）可得薄板应变能为 对于受垂直薄板中面、集度为q(x,y)的分布荷载作用的 板弯曲问题，外力势能 为 由此可得线弹性小变形的薄板势能为 （7-5） （7-6） （7-7） 7-2 矩阵（12自由度）薄板单元分析 设薄板被离散成若干矩形单元的集合，单元的结点位移与结点力（正向）如 图7-3所示。 图7-3 矩形单元结点位移和结点力 若引入如下结点位移和结点力矩阵符号 （7-8） 则单元结点位移矩阵和单元结点力矩阵分别为 因为单元结点位移参数（每结点的挠度和绕两坐标轴的转角） 总计有12个（故称12自由度），所以从广义坐标法角度来说 位移模式可取作 （7-9） （7-10） 式中四次项所以取 是为了 保证坐标互换时的不变性和曲率、扭率具 有相同方次（后面还将说明）。对式（7-10）求导可得 （7-11） （7-12） 利用式（7-10）、式（7-11）和式（7-12）及4个结点的位移条件即可确定