宏观三角形分区板壳单元的研究.pdfVIP

建筑结构分会2006年年会学术会议论文 2006年 宏观三角形分区板壳单元研究 曹杨 李杰 (同济大学建筑工程系，上晦200092) 摘要：基于经典平板壳理论和假设自然应变方法，研究了一种新型的宏观三角形分区(Macro Triangular Parti-tion，MTP)平板壳单元。引入内部虚拟节点后，所建单元将一个常规三角形平板壳单元细分为3个了三角形甲元， 通过特定取样点处的坐标组台以及相应的应力应变关系束不断更新各子单元区域的刚度矩阵，然后滤m虚拟节点形成 宏观三角形分区扳壳单元的最终刚度矩蛘表达式。一方面，由于运用了分区再整台的思怨，该攀元最终的表达形式只 具肓极少的自由度，从而大大地节省结构分析中烦冗的计算量；另一方面。该单元考虑，横向剪切效应的影响，对壳 厚具有较强的适应性。可以显著地克服以往针对板壳结构只进行二维分析所带柬的结果偏差。已有的计算结果显示， 基于通用有限元软件ABAQUS提供的用户定义单元(UEL)功能，宏观三角形分区平板壳单元可以实现较为精确的板壳结 构分析。 关键词t宏观三角形分区平板壳单元；虚拟节点；取样点i横向剪切效应；用户定义单元(UEL) 1 概述 随着有限元仿真技术的迅猛发展，三节点三角形平扳壳单元具有适应分析各类复杂形状空间结 构的拓扑优势。然而，由于剪切闭锁问题眭期以来困扰着三角形平扳壳单元的广泛应用，使得人们不 得不引入非连续剪切应变间断因子或剪切惩罚修正参数等人为因素以及离散Kirchhoff理论、杂交单 元等特定处理方式来回避剪切闭锁效应。 当前的研究结果表明，通过预先假设单元内部特定取样点处自然应变分量与节点位移之间多项 Natuml 式分布的方式，假设自然应变(AssumedStrain，ANS)已经成为～种避免剪切闭锁和膜闭锁效 应的有效手段，因而被广泛地应用于各类有限单元形式。本文所研究的三角形宏观平板壳具有以下特 征：①有效地避免包括剪切闭锁和膜闭锁效应；②最终单元格式仅基于简单的三节点三角形构型，每个 节点具有3个平动自由度和两个转动自由度并且满足Co连续性；③单元月4度矩阵不随节点编号顺序 的改变而改变；④不需引入任何依赖于分析的人为参数米控制剪切应变能。 2单元推导 2．1板壳单元基本构型 本文研究的板壳单元是基于考虑一阶剪切变形的经典壳体理论基础上建立起来的(图1)。 X=局+孝瓦 (1) 曾杨(1978年)．盘-博士研究生． 基金项目：国家“十五”攻关资助项目(2002BAB06B-4)：国家创新研巍群体资助项目(． ·361· 建筑结构分会2006年年会学术会议论文 2006年 式中x代表壳体内部某个材料点的位置矢量，由分别表示该材料点在壳体中面所对应的位置矢量X。和 垂直于中面的位置矢量五两部分组成，而自然坐标f(一l≤掌≤1)代表该点在厚度方向的位置。矢量 U表示该材料点变形前后的位置变化，关系如下“1： u=G+亭％=(同“一凰)+善(霹“一瓦) (2) 由于在经典板壳理论中不允许厚度发生变化，因而壳体法向矢量的模量为常数，即： I嚣”l-矧=叫2 (3) 蕾l垤^壳体嚣话的几何槽形 田2三角形^糟平援蠹单元的节点自由度爰面积坐标体幕 2，2常规ANS平板壳有限单元格式 常规三角形ANs平板壳单元包含3个节点，每个节点具有5个自由度：3个平动自由度(“。，v』，Ⅵ1 和2个转动自由度(q，尼)。其中平动自由度定义在全局坐标系中，而转动自由度则定义为沿青点处单 位表面切线矢量K。和K。的转动量。面积坐标的作用是构造平板壳