岩土力学 ›› 2021, Vol. 42 ›› Issue (11): 3137-3146.doi: 10.16285/j.rsm.2021.0559

• 数值分析 • 上一篇    下一篇

有限单元内部场量的高精度求解

次惠岭1, 2, 3,白冰1, 2, 3,雷宏武1, 3,崔银祥1, 3   

  1. 1. 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室,湖北 武汉 430071;2. 中国科学院大学,北京 100049; 3. 中国科学院武汉岩土力学研究所 环境岩土工程湖北省重点实验室,湖北 武汉 430071
  • 收稿日期:2021-04-14 修回日期:2021-09-23 出版日期:2021-11-11 发布日期:2021-11-12
  • 通讯作者: 白冰,男,1980年生,博士,研究员,博士生导师,主要从事能源环境与生态岩土工程等方面的研究工作。E-mail: bai_bing2@126.com E-mail:huilingci005@163.com
  • 作者简介:次惠岭,女,1995年生,硕士研究生,主要从事节理岩体的流固热耦合特性等方面的研究工作。
  • 基金资助:
    国家自然科学基金项目(No. 41972316)。

A high-precision scheme for field variables in finite element method

CI Hui-ling1, 2, 3, BAI Bing1, 2, 3, LEI Hong-wu1, 3, CUI Yin-xiang1, 3   

  1. 1. State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan, Hubei 430071, China; 2. University of Chinese Academy of Sciences, Beijing 100049, China; 3. Hubei Key Laboratory of Environmental Geotechnology, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan, Hubei 430071, China
  • Received:2021-04-14 Revised:2021-09-23 Online:2021-11-11 Published:2021-11-12
  • Supported by:
    This work was supported by the National Natural Science Foundation of China (41972316).

摘要: 有限单元法是目前最为广泛使用的数值分析方法,能够有效求解各类工程和科学问题。但在一些需要求单元内部场变量的情况下,传统形函数插值方法往往不能保证其结果精度。针对此问题,提出了一种高精度求解方案,基本思路是将有限单元与Taylor展开相结合,建立并求解单元内部未知场量关于已知单元结点场量的线性方程组,各方程中系数取决于点的相对位置关系。研究发现:所提方案算法简单,可适用于多维度、多阶次单元,在非线性场和高阶场中性能卓越,是一种高效高精度的新型求解算法,可服务于非线性计算的中间步骤,及满足高精度后处理的需求。该方法并不局限于有限元方法,可推广用于包括无网格方法在内的各类离散数值方法。

关键词: 有限单元法, 场量, Taylor展开, 高精度

Abstract: Finite element method (FEM) is the most widely used numerical analysis method, which is an effective solution to all kinds of engineering and scientific problems. However, the interpolation based on shape function often does not have desired accuracy in cases where the internal field variables of the finite element need to be solved. To deal with this kind of problems, a high-precision solution is proposed, and the basic idea of the solution is the combination of finite element and Taylor expansion. The value of internal field is obtained by establishing and solving linear equations of unknown field variables in the element with respect to known field values at element nodes. The coefficients of each equation depend on the relative position of the points. The results illustrate that the proposed method, which implement a relatively simple algorithm, is applicable to multi-dimension and multi-order elements. It also shows an excellent performance in the application of nonlinear and high-order field functions. This new solution algorithm with high computational efficiency and high precision is expected to serve the intermediate steps of nonlinear computation and satisfy the requirements of high precision post-processing. The method is not limited to FEM and can be extended to all kinds of discrete numerical methods including meshless method.

Key words: finite element method, field quantity, Taylor expansion, high precision

中图分类号: 

  • O 241
[1] 赵海鹏, 黎学优, 万建宏, 郑翔之, 刘思威, . 基于高性能有限单元法的多层地基水平 受荷桩受力分析[J]. 岩土力学, 2021, 42(7): 1995-2003.
[2] 高磊, 周乐, 刘汉龙, 黄坚, 王洋, 韩川, . 灌注桩桩身变形高精度测量现场试验研究[J]. 岩土力学, 2021, 42(7): 2004-2014.
[3] 宋义敏, 凌小康, 张敬宗, 朱晨利, 任何, 苑德顺. 基于数字散斑相关方法和有限元法的岩土 材料力学参数反演[J]. 岩土力学, 2021, 42(10): 2855-2864.
[4] 周凤玺, 马强, 周志雄, . 二维地基中空沟−波阻板联合隔振屏障分析[J]. 岩土力学, 2020, 41(12): 4087-4092.
[5] 姚仰平,王俊博. 北斗卫星定位在机场高填方远程变形监测中的应用[J]. , 2018, 39(S1): 419-424.
[6] 罗 滔,E. T. Ooi,A. H. C Chan,傅少君,. 一种模拟堆石料颗粒破碎的离散元-比例边界有限元结合法[J]. , 2017, 38(5): 1463-1471.
[7] 张 锟,徐 青,王一凡,阿湖宝. 自适应差分进化算法在边坡滑面搜索中的应用[J]. , 2017, 38(5): 1503-1509.
[8] 金晓飞 ,梁书亭 ,朱筱俊 ,张玉良,. 软土地区深基坑连续墙最大变形简易计算方法[J]. , 2015, 36(S2): 583-587.
[9] 刘远锋 ,郑 刚,. 地基排水控制法的沉降和孔压数值分析[J]. , 2015, 36(S1): 592-596.
[10] 赵明华,张 锐,刘 猛. 下限分析有限单元法的非线性规划求解[J]. , 2015, 36(12): 3589-3597.
[11] 扈胜霞 ,陈育民 ,闫竹玲,. 弹黏塑性模型及其在预压处理软基沉降计算中的应用[J]. , 2014, 35(4): 1173-1180.
[12] 冯又全,杨 敏,陈俊岭. 线性分布基床系数弹性地基梁有限单元法改进[J]. , 2014, 35(10): 3027-3034.
[13] 林永生,陈胜宏. 基于有限元计算的边坡三维滑裂面搜索[J]. , 2013, 34(4): 1191-1196.
[14] 杜佐龙 ,黄茂松 . 非均质与各向异性黏土基坑抗隆起稳定分析[J]. , 2013, 34(2): 455-461.
[15] 彭芳乐 ,李福林 ,白晓宇. 针对砂土应变软化强非线性问题的动态松弛有限元法研究[J]. , 2012, 33(2): 590-596.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 魏 丽,柴寿喜,蔡宏洲,王晓燕,李 敏,石 茜. 麦秸秆加筋材料抗拉性能的实验研究[J]. , 2010, 31(1): 128 -132 .
[2] 黄庆享,张 沛,董爱菊. 浅埋煤层地表厚砂土层“拱梁”结构模型研究[J]. , 2009, 30(9): 2722 -2726 .
[3] 荆志东,刘俊新. 红层泥岩半刚性基床结构动态变形试验研究[J]. , 2010, 31(7): 2116 -2121 .
[4] 刘争宏,廖燕宏,张玉守. 罗安达砂物理力学性质初探[J]. , 2010, 31(S1): 121 -126 .
[5] 王登科,刘 建,尹光志,韦立德. 突出危险煤渗透性变化的影响因素探讨[J]. , 2010, 31(11): 3469 -3474 .
[6] 樊恒辉,高建恩,吴普特,娄宗科. 水泥基土壤固化剂固化土的物理化学作用[J]. , 2010, 31(12): 3741 -3745 .
[7] 张成平,张顶立,骆建军,王梦恕,吴介普. 地铁车站下穿既有线隧道施工中的远程监测系统[J]. , 2009, 30(6): 1861 -1866 .
[8] 王 军,曹 平,李江腾,刘业科. 降雨入渗对流变介质隧道边坡稳定性的分析[J]. , 2009, 30(7): 2158 -2162 .
[9] 张 渊,万志军,康建荣3,赵阳升. 温度、三轴应力条件下砂岩渗透率阶段特征分析[J]. , 2011, 32(3): 677 -683 .
[10] 张雪婵 ,龚晓南 ,尹序源 ,赵玉勃. 杭州庆春路过江隧道江南工作井监测分析[J]. , 2011, 32(S1): 488 -0494 .