基面力元法(BFEM)的研究与应用进展.pdfVIP

基面力元法（BFEM）的研究及应用进展1) 2) 彭一江 ，宗娜娜，董占龙 （北京工业大学工程抗震与结构诊治北京市重点实验室，北京100124） ― Base 摘要：本文从高玉臣院士提出的基面力概念出发，针对一种基于基面力概念的新型有限元法 基面力元法 （ Force Element Method, BFEM ）的研究及应用进行综述。分析结果表明：基面力元法的计算所得结果与理论解 吻合较好；基于余能原理的基面力元法较传统平面 4 节点等参元（Q4 单元）具有较高的计算性能；在高度几 4 Q4R 何非线性分析中，基于余能原理的基面力元法较平面 节点减缩积分单元（ 单元）也具有一些优势；基面 力元法列式推导简单，是有限元方法的一种新思路，具有较好的应用前景。 关键词：基面力；有限元；基面力元；势能原理；余能原理 1 基面力元法综述 近年来，随着工程与材料科学的迅速发展，对数值模型的适用性、可靠性和非线性行为提出了新的 更高的要求。但是，一些改进的有限元法由于自身理论体系的复杂性，在非线性领域的应用由于受到各 种因素的制约，远没有在线性领域的应用成熟。目前，有限元法的进一步发展仍是国际计算力学界研究 的热点问题。此外，由于上述有限元法均采用二阶张量的描述方法，因此在数学模型的推导和建立较为 复杂，寻求有限元法的进一步突破较为困难。 近二十几年来,高玉臣院士的工作给出了弹性大变形新的描述方法[1,2]，建立了两种典型本构关系[3,4]， [5-21] [1] 并解决了一系列带有奇异点的代表性问题 。利用高玉臣提出的“基面力”（Base forces ）概念 ，可 以完全替代传统的各种二阶应力张量描述一物质点在初始构型和当前构型的应力状态，得到弹性力学基 本方程（平衡方程、边界条件、本构关系）的简洁表达式，建立势能原理和余能原理。在研究物体的力 学行为，特别是有限变形的分析中，一阶张量基面力具有传统的二阶应力张量无法比拟的优越性，提供 了一个很好的分析工具。2003 年高玉臣[2]较系统地提出了基面力的理论体系，并利用基面力概念给出了 推导空间任意多面体单元刚度矩阵和柔度矩阵显式表达式的思路，为基面力概念在有限元领域的应用奠 定了理论基础。在高玉臣工作的基础上，本文第一作者从 2003 年开始进行基于基面力概念的新型有限元 [23] 方法―基面力元法（Base force element method, BFEM ）的研究 。 在线弹性问题的余能原理基面力元法研究方面，2006 年彭一江和金明[24]利用 Lagrange 乘子法推导 出线弹性问题的余能原理有限元控制方程，利用 Matlab 语言编制了基于余能原理的基面力元分析软件及 1) 国家自然科学基金资助项目 2) E-mail ：pengyijiang@bjut.edu.cn 1 前处理程序，给出了具有边中节点的四边形基面力元的研究成果；2007 年彭一江和金明[25]利用基于基面 力概念详细推导出单元柔度矩阵的具体表达式，并进一步对具有边中节点的四边形基面力元进行了对比 分析；2007 年彭一江和金明[26]研究了基于余能原理的任意多边形网格基面力元法，利用 Matlab 语言编 制了基面力元分析软件及任意多边形网格剖分的前处理程序，给出了具有边中节点的任意多边形基面力 元的研究成果；2008 年彭一江、雷文贤和彭红涛[27]从描述平面应力状态的基线力概念出发，推导出具有 边中节点的四边形基面力元的柔度矩阵展开式及矩阵表达式，通过与常规有限元