有限元课件ch7板壳结构.pptVIP

第7章 板壳结构分析 2、薄板理论 基本假设 （1）直法线假设：薄板中面法线变形后仍保持为法线。由此，板中面内剪应变为零。 （2）忽略板中面的法线应力分量，且不计其引起的应变。 （3）薄板中面内的各点没有平行于中面的位移，即中面不 变形。 利用上述假设将平板弯曲问题转化为二维问题，且全部应力和应变可以用板中面挠度w表示。 基本方程 （1）位移：由假设（1）、（3），有 (2)应变 由假设（1）、（2），薄板弯曲问题只需要考虑三个分量。根据几何方程，应变可表示为 形变分量：中面x和y方向的曲率与x，y方向的扭率。 应力 内力：板单位宽度上弯矩Mx 、 My和 Mxy ，为应力分量在板截面上的合力矩 ： 薄板弯曲问题中的弹性矩阵 [Df] 7.2 薄板矩形单元[shell63] 结点位移 单元结点力 将结点坐标和结点位移代入上式），可解出a1~a12，再代入该式并整理得位移函数 3、单元应变 3、单元应力 5、单元刚度矩阵 基本公式 6、等效结点力 板单元受横向均布载荷p作用，则 等效结点力为 受中心集中力的四边支承板的计算结果(边长为1，厚度为0.01，弹模为1，波松比为0.3) 7.3 薄板三角形单元 1、位移模式 三角形单元能较好地适应斜边界，实际中广泛应用。单元的结点位移仍然为结点处的挠度wi和绕x，y轴的转角θxi、θyi，独立变量为wi。三角形单元位移模式应包含9个参数。若考虑完全三次多项式，则有10个参数: 对于三角形单元，面积坐标的一、二、三次齐次分别有以下项： 将三个结点的位移和面积坐标代入上式，可得：a1=wi ， a2=wj， a3=wm。代入上式对Li，Lj求导，注意Lm=1-Li-Lj，可得 将结点的面积坐标代入上述两式，可得6个关于a4~a9的方程，求解后可得a4~a9： 最后，待定常数a1~a9代入位移模式，整理后得： 利用求得位移函数，可以得到应变列阵和相应的应变矩阵[B]，进一步可得到形变列阵{Χ}和相应的形变矩阵[B’]。 四边将简支板的中心挠度系数 7.4 考虑横向剪切影响的平板弯曲单元 在薄板单元中，构造协调单元的困难在于单元间要求斜率的连续性。如果放弃薄板理论的直法线假设，考虑横向剪切的影响，有可能绕过这一困难。假设：中面法线变形后仍为直线，但绕x、y轴转动了θx、θy.。 代入几何方程，应变矩阵： 平板的变形由中面挠度w和法线绕x、y轴的转角θx、θy.确定。每个结点取它们作为自由度，采用8结点平板单元。可以参照8结点等参单元。 中面上任意点的挠度和转角可以表示为： 内力计算 单元刚度矩阵 等效结点荷载： 设单元表面作用有均布荷载q(x,y)，等效结点荷载为 7.5 平面壳体单元 平面壳体单元的应力状态是由平面应力和弯曲应力单元的叠加，因此，构造壳体平面单元时，只需要将前面的两种单元简单叠加组合即可。值得注意的时，壳体平面单元用于分析壳体结构时，需要对单元刚度矩阵、等效结点荷载和结点位移等进行坐标变换。 结点位移：5个位移，即{ui,vi,wi,θxi, θyi},前两个对应为平面应力问题，后三个对应平板弯曲问题。对应结点力为 在局部坐标中，节点位移不含θz，但为了将局部坐标下的刚度矩阵转换到整体整体坐标系，须将θz加入节点位移中。即： 7.6 考虑横向剪切影响的壳体单元 采用平面单元模拟曲面结构，不够理想，而采用曲面单元效果更好，单元数量也会相应减少。 采用薄板理论的直法线假设使得中面转动依赖于中面位移，给位移模式的构造带来困难。如果考虑横向剪切变形的影响，就可以认为中面转动是独立变量而不在依赖于位移的一阶导数。因此，只需要单元边界上的位移函数的连续性就可以，从而避开了要求一阶导数的连续性。 实际中常用八结点40自由度四边形单元来考虑剪切变形对壳体单元的影响。该类单元适合于薄壳和厚壳。 1、单元几何形状的确定 8结点壳单元，象类似空间等参单元一样，引入一个自然坐标系oξηζ。命ξη为壳体中面上的曲线坐标；ξ=1对应顶面， ξ=-1对应底面。在单元中面上选8个结点，过各结点i(i=1,2…8)作中面的法线，交顶面和底面的点称为结点i的对点。其坐标记为： 单元类任意一点的坐标可以通过形函数Ni(ξ,η)的插值表示： 2、位移模式 假定中面法线V3i变形后仍为直线，但是不再是变形后的中面法线，有绕V1i和V2i的转角βi、αi。则单元内任意点的位移可以表示为 写成标准形式 应变计算 应力计算（单元坐标） 由此可得位移模式 应变分量 应力分量 例：承受均布荷载q的方板和板中心处受集中荷载P的方板，四边固定和简支。4X4网格，挠度？ 0.04437 0.07557 0.07500 0.4 0.04437 0.0