板壳有限单元法的变分.PDFVIP

第 16卷第 5期 水 利 水 电 科 技 进 展 1996年 10月 板壳有限单元法的变分原理与单元类别综述 ! 邵国建 卓家寿 一 ’ (河海 大学土木工程 学院 南京 210098) {摘要 基于不同种类的泛函 式 ，对 目前有限元法分析板、壳 问题时所采用的各种单元进行归类 ， 介绍非协调元 、协调元 、混台元 、杂交元 、广义协调元及拟协调元等在 实施过程 中的处理方法 ，井 分 析各种单元 的适用场台及优缺 点。 关键词 兰 墨 挛茎垦銎 墼 褪 皇元堇 有 限元法中首要 问题是场变量选取及其 ① = (OIl 1) 插值 函数 的构成 ，对于要求 c0连续性 的插值 函散 的构成 比较 简单 ，而 以板 、壳弯曲为代表 ② 1= ，M l= M (OIl2) 结构分析时要求场变量具有 c 连续性 的插 ③f )f．： (0n 值函数的构成 ，除了一维 问题 (如果 梁弯曲、 轴对 称壳)中易于实施外 ，比 C 连续 问题 中 基 于式 (1)泛函以 作 为场 变量构造单 要复杂得多 ，这 导致 用有 限元法分析板 、壳问 元 ，每个单元结点有挠度 (1个)、转动 (2个) 题 时遇到麻烦 ，这 也是人们长期从事这方 面 共 3个参数 ，根据 单元 结点数及 结 点参数选 研究 的原 因。板 壳单 元按构造 单元种 类来分 取 ，可建立适用于规则边界 的薄板 四结点非 有非协调元…、协调元 J、位移和转动分别插 协调矩形元 ，见图 1。适用 于任 意边界 的三结 值 的板单 元 、离散 Kirehhoff单 元 【、混合 点非协调三角形单 元 ，见 图2。需指 出采用三 元 j【J、杂 交 元 J、广 义 协 调 元 及 拟 协调 0— ～ ．| 元_日J。事实上 ，构造上述各种单元基于不 同类 型泛 函 下面 以薄板结构为主要分析对 象 ，就 构造不 同种类单元所采用 的各种泛 函进行分 析 ，其理论也可推广应用于一般板壳结构 。 1 薄板非协调元及协调元 势能泛函表示 为 = 』『(1 一)y一 丽 ㈤ 式 中 ：L∞为广义应变 ；D为弹性矩阵； q为横 向分布载荷 ； 为挠度 。 图 2 三结 点三 角形板单元 边界条件 ： 角形板 单元时宜用面积坐 标。这类 单元 可保 个边界中结点的a ／an值不再作 为独立 的 证单元边界面 和 a ／as协调 ，而无法保证 参数 ，而是作为限制各个边界上a ／an的三 a n协调 ，但 计 算 结果表 明效果 良好… 。 次变化 的附加 条件 。多结 点参 数 的四边形单 这是 因为利用最小势能原理求得 的近似解 一 元使用较少 ，尽管矩形单元可采用 Hermite多 般使 结构呈现过 于刚硬 ，而非协调 元虽未精 项式作为位移 函数 ，但这种单元难 以适应一 确满 足最 小势能原理 的要求 ，但 在单元交 界 般几何形状 的结构 ，而一般 四边形板 单 元除 面上有较多 的 自由度 ，使结构趋于柔软 ，正好 取三角形组合 的方 法外 ，直接建立 单 元的位 部分 地抵 消 了上述 过 于 刚硬 带来 的误 差 。 移函数 比较复杂 。 Irons[9_提 出了分片检验方法 ，为非协调元