变厚度曲壳元在蛋形消化池计算中的应用研究.pdfVIP

第十届全国工程设计计算机应用学术会议论文集 广东．广州2000 变厚度曲壳元在蛋形消化池计算中的应用 昊玉华 楼文娟 浙江交通设计事务所 浙江大学土木系 、{I自，^f{■，lI{4‘●o{l=●《}}屯#r、{l球j《 【摘要】对于一个变厚度双曲壳体结构，如果采用平面壳体单元．将会引起几何E近似所带来的讽差． 而采用变犀度曲面壳单元可以很好地描述壳体的真实几何形状，从而褥到具有更高精度的*簪．因此， 本文采用一种变犀度八绪点四十自由度的曲面单元．井用c增言将其辘橱成实用程序，对杭州四量污水 处理厂的蛋形澹化池进行了计算．计算结果表明，该有瞑元形式是一种黻更高．收敛更快的有限元模 式． 1引言 随着国民经济的发展和人民生活质量的提高，环境工程已越来越受到人们韵重视，而 污承处理工程也艘视为现代城市环境质量的一种标志和文明的象征。近十余年来，为适应 改革开放豹鲁要，城市污承她理厂豹建设日趋增多，污水处理构筑物的容量也日趋增大。 困此，采用结构设计台理．容量丈．经侪可靠舳消化池显得十分重要。预应力蛋形消化池 以其在结构外形，容量以及技术经济效益等方面的优越牲．而将逐步成为污泥消化抛舶发 展方向。 蛋形消化池是—种双曲变厚度薄壳结构，常用的平板壳单元虽然计算簿单，但平壳单 元除了由于选取的拉移模式所弓I入的误差井，还引进了且何形状近自i所带来舶误差。这种 误差虽然随着单元教蚋增加而减小．但对于大型结j町，增加单元数将大大增加计算量和计 算时间，甚至超出计算软件对结点自白瘦总量的限制．此外．平壳单元只能处理等厚度壳 单元，因此，在划分有限元罔格和输入单元信息时，通常需要计算出每～单元豹平均厚度。 而采用曲面单元可以描述壳体的真实几何形状，使之用不太多的单元来代替复杂形状的壳 体，将得到具有更高精度的解答。同时对于变厚度壳体，曲面单元艟自动计算每一结点处 的厚度．大大减少了这部分的工作量。 2曲面单元公式 图l所示为考虑横向剪切变形的八结点四十自由度舶 曲面壳单元，在壳单元中引入一个自然坐标系。翻‘．令岛q 为壳体中面上的曲线坐标，对应于(．1的表面称为顶面(或 上表面)。对应于‘一l的表面称为底面(或下表面)。在单 元的中面上选取八个点称为结点，过备结点i(i=l2，…8) 作中面的法线，交顶面和底面的称为结点i的对点，则1 中面上结点t的整体坐标值是； 田I八结点四十十自度的曲壳单元 苎±鱼全璺兰堡堡苎生墨!生旦茎苎全垫堕查苎 ￡查：￡|州!!!1 2．1 = + l一2 ，耋=．㈨旧、 、●，●●●●●●L，●●，●●J 啦、●●●●●t，●●●J 璜 ，●●●●●●●●，，●●●●●●L t只‘ 显然．结点t处中面法线方向可以由下列单位矢量所确定： 。={是。}=砉({耄}礓一{蓼}矗] 2—2 式中．矗，，mm和码，是结点i处中面法线方向对于整体坐标轴∞9疆的方向余弦·而槐是 结点i处的壳体厚度。 于是．单元内任意点的坐标值可以通过形函数以睡’1)的插值得到，即 {i}=喜Ⅳ，g，叩)[售}+鲁f以』] 2—3 懈稚讥姗 胛：璧{重}蛐