[高等教育]有限元板壳——王勖成.pptVIP

第十章 平板弯曲问题 主要内容： 薄板弯曲理论基本假设和基本方程 基于薄板理论的非协调板单元 基于薄板理论的协调板单元 1、薄板弯曲理论基本假设和基本方程 基本假设 ——Kirchhoff假设 （1）直法线假设：薄板中面法线变形后仍保持为法线。由此，板中面内剪应变为零。 （2）忽略板中面的法线应力分量，且不计其引起的应变。 （3）薄板中面内的各点没有平行于中面的位移，即中面不 变形。 利用上述假设将平板弯曲问题转化为二维问题，且全部应力和应变可以用板中面挠度w表示。 基本方程 (1)位移：由假设（1）、（3），有 (2)应变 由假设（1）、（2），薄板弯曲问题只需要考虑三个分量。根据几何方程，应变可表示为 形变分量：中面x和y方向的曲率与x，y方向的扭率。 应力 内力：板单位宽度上弯矩Mx 、 My和 Mxy ，为应力分量在板截面上的合力矩 ： 薄板弯曲问题中的弹性矩阵[D] 平衡方程 由广义应力应变关系及几何关系代入平衡方程得由W的微分方程： 边界条件 （1）位移边界条件 （2）混合边界条件 其中 （3）力边界条件 泛函表达式： 2、基于薄板理论的非协调板单元 结点位移 将结点坐标和结点位移代入上式，可解出a1~a12，再代入该式并整理得位移函数 单元刚度矩阵基本公式 等效结点力： 2.2 薄板三角形单元 1、位移模式 三角形单元能较好地适应斜边界，实际中广泛应用。单元的结点位移仍然为结点处的挠度wi和绕x，y轴的转角θxi、θyi，独立变量为wi。 对于三角形单元，面积坐标的一、二、三次齐次分别有以下项： 将三个结点的位移和面积坐标代入上式，可得：a1=wi ， a2=wj， a3=wm。代入上式对Li，Lj求导，注意Lm=1-Li-Lj，可得 将结点的面积坐标代入上述两式，可得6个关于a4~a9的方程，求解后可得a4~a9： 最后，待定常数a1~a9代入位移模式，整理后得： 3、基于薄板理论的协调板单元 1）增加结点参数 结点参数中包含w的二次导数 项 2）保持结点3个参数前提下采取其他措施，如附加校正函数法、在分割法等。 1、3结点的协调元 引入校正函数，eg在2-3边上引入φ23，其性质为： 1）在全部边界上φ23=0； 2）在边界1-2和1-3上 3）在2-3上 ,在2-3的中点4取单位值. 类似可有φ12、 φ31.则位移函数为 其中， 是非协调元的位移函数， 、 、 为待定常数，可以通过调整其大小使相邻单元交界面上的协调性得到实现。 实现方法： 通过调整的大小，使得 在各个边界中点的值等于各个边界两端结点 值的平均值，即w在各个边界法向导数 成线性变化 则 3、多结点参数的协调元 问题： 1、为什么结点位移参数 要取 ？ 2、在矩形单元中四结点12自由度薄板单元位移四次项为什么选x3y和xy3,而没有选其它四次项？选用x4和y4将导致什么结果？ 3、四结点12自由度薄板单元挠度和转角协调，法向转角不协调，所以称其为非协调元，如何来证明？ 4、对于非协调的矩形单元除广义坐标法外还有什么方法可更容易的确定插值函数？ 广义应变 弹性矩阵 内力矩表示薄板应力的公式 薄板应变能： 单元结点位移列阵 1)、结点位移 2.1、矩形板元 2、位移模式 矩形薄板单元有4个结点，12个结点位移分量，1个挠度独立变量，根据选取位移函数的原则，取： 1-3项刚体位移 4-5项常应变 非协调元 单元间法线导数可能不连续 式中形函数 广义应变 广义内力 板单元受横向均布载荷p作用，则 等效结点力为 非协调板单元可以通过分片试验，当单元划分不断缩小时，计算结果可以收敛于精确解答，但是收敛并非一定是单调的，即不一定是精确解的上界或下界。 三角形单元位移模式应包含9个参数。若考虑完全三次多项式，则有10个参数。若以此为基础构造位移函数，则必须去掉一项。无法保证对称。经过许多研究，问题最后在面积坐标下得以解决。 Li Lj Lm及其一阶导数在三个结点为零， 对于确定待定参数无用；考虑到用结点位移表示待定参数时计算方便，不考虑二次和三次的前三项。因此，只能在剩下的6个三次项中选择三个或利用某种线性组合。考虑对称性，假设位移模式为： 将w,Lii和w,Lji变换成θxi、 θyi，从而得到相应于θxi、 θyi的形函数Nxi、 Nyi 利用： 非协调板单元的收敛性是以分片试验为条件的即使收敛也非单调，不能对解的上下界作出估计。 经典薄板理论范围内使板单元满足协调性要