平面应变问题广义有限元法-计算物理.PDFVIP

第３３卷 第３期 计　 　 算　 　 物　 　 理 Ｖｏｌ．３３，Ｎｏ．３ ２０１６年５月 　 ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＣＯＭＰＵＴＡＴＩＯＮＡＬ ＰＨＹＳＩＣＳ　 Ｍａｙ，２０１６ 文章编号：１００１⁃２４６Ｘ（２０１６）０３⁃０３５８⁃０９ 平面应变问题广义有限元法 杨　 璞，　 牛红攀，　 肖世富 （中国工程物理研究院总体工程研究所，四川 绵阳　 ６２１９００） 摘　 要：结合广义有限元和理性有限元，针对平面应变问题提出一种广义四边形单元．该单元考虑泊松效应，以节 点位移自由度约束弹性力学平面应变方程的半解析解，构造单元位移模式的附加项，较准确地反映真实位移场，提 高了单元的计算精度．首先推导广义单元及其等参单元的形函数，之后设计分片试验和数值算例来验证单元的精 度，数值算例表明，在规则和非规则网格下新单元计算精度均优于传统有限元和广义有限元．新单元精度高且易于 程序实现，可应用于实际工程的结构分析． 关键词：平面四边形单元；平面应变；广义有限元；等参元；形函数 中图分类号：Ｏ２４２２１ 文献标志码：Ａ ０　 引言 有限元法基本思想是将一个连续体的求解域离散为若干子域，并通过边界上的节点连接成为组合体，用 每个单元内的近似函数来表示求解域内的未知场变量，通过和原问题数学模型等效的变分原理或加权余量 ［１］ 法，建立代数方程组以求解基本未知量 ． 传统线性有限元直接以数学为基础采用双线性多项式构造形函数，其位移模式没有考虑弹性力学控制 方程，也不能反映不同自由度间位移的相互影响，单元精度较低．为改善计算精度，上世纪９０ 年代钟万勰院 士团队研发了理性有限元，核心思想是从弹性力学出发建立位移模式，其插值函数的各项均来源于弹性力学 精确解，有明确的物理意义［２－７］．但理性有限元形函数繁复，给实际操作带来很大困难［８－９］．近年来，张洪武 教授团队在理性有限元的基础上，进一步发展了广义有限元方法，核心思想是改变传统有限元只利用单元节 点同方向自由度的位移构造方式，以单元所有节点自由度来构造单元位移［１０－１１］ ，其位移模式仍沿用了传统 有限元的多项式插值函数，并没有有效利用弹性力学解． 本文针对平面应变问题，提出一种广义四边形单元，综合考虑理性有限元和广义有限元的优势，从弹性 力学平面应变方程出发建立形函数，并在不改变单元自由度的情况下为形函数引入位移附加项，提高了单元 计算精度．由于新单元只在传统有限元基础上增加位移附加项形函数，因此程序操作简单易于实现． １　 平面应变新型广义四边形单元 １１　 一般单元 传统平面矩形单元具有４个节点８ 自由度，如图１所示．用４个节点位移表达的位移模式为 ｕ ＝ Ｎ ｕ ＋Ｎ ｕ ＋Ｎ ｕ ＋Ｎ ｕ ， { １ Ｉ ２ Ｊ ３ Ｋ ４ Ｌ （１） ｖ ＝ Ｎ ｖ ＋Ｎ ｖ ＋Ｎ ｖ ＋Ｎ ｖ ， １ Ｉ ２ Ｊ ３ Ｋ ４ Ｌ 其中插值函数为双线性函数 Ｎ ＝ （１ －ｘ）（１ －ｙ）／４，　 Ｎ ＝ （１ ＋ｘ）（１ －ｙ）／４， １ ２ {Ｎ ＝ （１ ＋ｘ）（１ ＋ｙ）／４，　 Ｎ ＝ （１ －ｘ）（１ ＋ｙ）／４， （２） ３ ４ 收稿日期：２０１５－０２－１７；修回日期：２０１５－０６－２５ 基金项目：中物院高性能科学与工程计算战略科技项目（Ｒ２０１４ －０２０３）及国家自然科学基金资助项目 作者简介：杨璞（１