水工结构有限元分析的前处理.PDFVIP

第 18卷第 1期 水 利 水 电科 技 进 展 1998年 2月 水工结构有限元分析的前处理 Tv3f 祖 波 侯 百军 ● — — — —— — — — — 郑州 45oo03) (黑龙江省水利厅 喑尔滨 ~5oooi) 摘要 利用超单元概：盎自动生成有 限元网格 ．建立了五面体和四面体超单元阿格生成模式，并编 制 了超单元 自动生成罔格法的有限元前处理程序。 关键词兰 ! 型 竖 利用有限元法进行结构的数值分析时工 元有八结点六面体单元和六结点五面体单元。 作量很大，如何提 高有限元分析效率是普遍 1．2．1 网格结点的局部坐标和整体 坐标 关注的问题，而有 限元前处理是 问题的焦点 图 1所示为五面体母单元和三 角形面的 之一。本文利用超单 元概念 【1j，编制 比较先 网格生成 ，假定母单元 }， ，}三方向生成的 进的三维块体结构 自动剖分程序 ，程序中超 结点数为m，n，f，其 中 m=n，三角形各边结 单元类型有六面体实体超单元、五面体 实体 点数相等且等分 ，f方向可不等分 ，剖分 比例 超单元、四面体实体超单元 、六面体无体积接 因子 ，考虑程序编制 的结构性 ，在 e， 两 触面超单元 、五面体无体积接触面超单元、六 方向仍采用剖分 比例 因子 ，只是 = =1。 面体无体积连接止水超单元。主要工作是建 结点局部坐标表示为 立 了五面体和 四面体超单元 网格生成模式。 = (∑ 一1)／(∑ 一‘1) ，-I JI1 1 超单元网格生成 I ■ 似=(∑ 一1)／(∑^一1) 1．1 六面体超单元网格生成 对某个空间八结点实体超单元 ，将其变 一 I = 2(∑^ 一‘1)／(∑ 一‘1)一1 换到正六面体母单元中，依剖分要求求 出母 · ’ ， · 单元的网格点局部坐标 ，由坐标变换解 出网 (1) 格结点的整体坐标 ．并相应地定义结点和单 元编号以及单元信息-lJ。对于八结点六面体 无体积接触 面和连接止 水超 单元 的网格 生 成，其剖分与实体单元一致，只是沿夹层厚度 方向不再进行剖分 。 1．2 五面体超单元网格生成 单元剖分中，五面体超单元经常用到，对 三角形面采用各边 的结点数相等且等分 ，以 图 I 五面体母单元及三角形面的罔格生成 等腰直角三棱柱作为母单元，剖分生成的单 第一作者筒舟：冀春接，男，博士，占L事水工结构研究 ．曾发表(水闸抗振的一种实用分析方告)等论文。 · 52r· 结点整体坐标表达式为 分与