具有旋转自由度ψz的平板型矩形厚壳单元.pdfVIP

———— 一 !三翌垄堂!望型垫!!兰 ：≥23： 具有旋转自由度％的平板型矩形厚壳单元 邹佑学匡文起龙驭球 f搏华丈学土术幕，北京10。。拍 ● 摘要本文在文献flt2】基础上．构造了广义协调厚壳单元GZR24M，用于壳体结构静、动力分析．算 例表明：壳单元GZR24M厚薄通用．收敛速度快．可咀求出较高精度的高阶频率．而且能避免单元共面 或接近共面出现的奇异性问题，适合于工程应用。 关键词 广义协调．旋转自由度，平板型壳元，厚薄通用，振动 1引言 在壳体结构分析中，由于其几何形状与受力状态的复杂性，寻求解析解比较困难，甚 至有时是不可能的，故较多地采用数值方法分析。有限元法是～种有效的数值分析方法。 通常的壳体单元可以分为三种类型：平板型单元、曲面型单元和三维实体型单元。平板型 壳体单元具有表达式简单，自由度少，易于实施等优点．因此拨广泛地采用。 由于受板弯曲元和平面应力元的限制．一般平板型壳元分析壳体结构存在三个问题， (1)板弯曲元和平面应力元是否能合理匹配；(2)当存在几个单元共面或接近共面时， 总体刚度矩阵发生奇异：(3)忽略了剪切变形的影响。对于壳体结构的振动分析，板弯曲 元和膜元是否合理匹配对计算结果的影响显得更为突出。 广义协调元已经成功的应用于薄板、中厚板和壳体计算的许多领域。本文采用构造的 ● 广义协调平扳型厚壳元GZR24M用于柱壳的静、动力分析。 ， 2广义协调矩形厚壳单元GZR24M 局部坐标系下的矩形平板型壳元如图1所示，共有二十四个自由度：局部坐标的x 和y轴位于壳元平面，并分别与单元两边平行，：轴与壳元平面垂直。每个角点 i(i=1,2，3，4)上有六个位移分量： r 1一 {吼}=k．V．Wi纵，∥‘y。J。 f1) 这里“。，v。，w。是局部坐标方向上的线位移，虬。，矿，，∥，是绕局部坐标轴的角位移；位 移仉正方向定义为绕x轴从y轴正向转至z轴正向．同理定义位移¨，性的正方向。 邹佑学，男，1974．4生，工学硕士 清华大学结构工程与振动开放研究实验室研究基金资助项目 ·324· 《工程力学》增刊2001年 毖 (a) (b) (c) ● 图1局部坐标下的平板型矩形壳元 单元自由度由四个角点位移组成： … {口。)7 {92)7{茸，)7 J1 {g)。=L{g。)7 本文研究的广义协调平板型矩形壳元GZR24M，其薄膜与弯曲状态分别由GRl2M平 面膜元…和LFRI厚扳弯曲单元[21表示。 2．1平面膜元GRl2M 具有旋转自由度矿，的平面膜元GRl2M的角点位移和单元位