›› 2016, Vol. 37 ›› Issue (10): 2984-2992.doi: 10.16285/j.rsm.2016.10.033

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

自由度具有物理含义的线性无关高阶数值流形法

徐栋栋1,杨永涛2,郑 宏2,邬爱清1   

  1. 1. 长江科学院 水利部岩土力学与工程重点实验室,湖北 武汉 430010; 2. 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室,湖北 武汉 430071
  • 收稿日期:2014-12-23 出版日期:2016-10-11 发布日期:2018-06-09
  • 作者简介:徐栋栋,男,1986年生,博士,工程师,主要从事计算岩土力学方面的研究。
  • 基金资助:

    国家自然科学基金(No.11502033, No.51579016 , No. 41272350, No. 41672320);国家重点基础研究发展计划(973)项目(No. 2011CB710603)。

A linearly independent high-order numerical manifold method with physically meaningful degrees of freedom

XU Dong-dong1, YANG Yong-tao2, ZHENG Hong2, WU Ai-qing1   

  1. 1. Key Laboratory of Geotechnical Mechanics and Engineering of Ministry of Water Resources, Yangtze River Scientific Research Institute, Wuhan, Hubei 430010, China; 2. State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan, Hubei 430071, China
  • Received:2014-12-23 Online:2016-10-11 Published:2018-06-09
  • Supported by:

    This work was supported by the National Natural Science Foundation of China (11502033, 51579016, 41272350 and 41672320) and the National Program of Key Basic Research Project of China (973 Program) (2011CB710603).

摘要: 数值流形方法(NMM)基于两套覆盖(数学和物理覆盖)和接触环路而建立,能够统一地处理岩土工程中的连续和非连续变形分析问题。与其他基于单位分解理论的数值方法一样,NMM可以自由地提高物理片上局部位移函数(多项式)的阶次,从而在不加密网格的情况下显著地提高计算精度,但有可能会使总体刚度矩阵奇异,产生线性相关问题。针对这种情况,引入了一种新的高次多项式形式的局部位移函数,在此基础上,建立了新的NMM求解体系,并应用于求解一般的弹性力学问题。结果表明:它有效地消除了线性相关问题;较之传统局部位移函数取一次多项式的NMM,达到了更高的精度;节点应力是连续的;定义在物理片上的所有自由度都具有明确的物理含义,其中第3~5个刚好是物理片所对应插值点处的应变分量,因此,直接获得此处的应力状态。该方法可以很容易地推广到其他基于单位分解的数值方法中。

关键词: 数值流形法, 物理片, 线性相关, 亏秩, Cook斜梁, 带孔无限板, 均质边坡

Abstract: Numerical manifold method (NMM) is established based on the two cover systems (including mathematical cover and physical cover) and the contact loop, and can be used to solve the continuous and discontinuous deformation problems in the geotechnical engineering in a unified way. Similar to other numerical methods based on the partition of unity (PU) theory, the NMM can also improve the computational accuracy by increasing the orders of local displacement functions freely without mesh refinement, though this may cause the global stiffness matrix singular, leading to the linear dependence issue. In this study, a new local displacement function of high-order polynomials is proposed. The new function is applied to solve the general elastic problems. The results show that the linear dependence is eliminated. Compared with the traditional NMM based on the first order polynomials, higher precision is reached. Stresses at nodes are continuous. All the degrees of freedom defined on a physical patch are physically meaningful, with the third to fifth simply being the strain components at the interpolation point of the patch. As a result, the stresses at the interpolation points can be directly obtained. The proposed procedure can be easily extended to other PU-based methods.

Key words: numerical manifold method, physical patch, linear dependence, rank deficiency, Cook skew beam, infinite plate with a circular hole, homogeneous slope

中图分类号: 

  • O 241

[1] 杨石扣, 任旭华, 张继勋, 艾华东, . 三维数值流形法覆盖系统生成算法研究[J]. 岩土力学, 2022, 43(S1): 633-640.
[2] 王秋生, 张瑞涛, 郑 宏. Malvern Hills边坡溃曲破坏分析及数值流形法模拟[J]. 岩土力学, 2022, 43(7): 1951-1960.
[3] 屈春来, 付迪, 刘世伟, 冷先伦, 李建贺, 孙熇远. 非均质成层边坡极限承载力上限分析[J]. 岩土力学, 2022, 43(10): 2923-2932.
[4] 杨亮, 杨永涛, 郑宏, . 相场数值流形法模拟岩石裂纹扩展[J]. 岩土力学, 2021, 42(12): 3419-3427.
[5] 聂治豹, 郑宏, 万涛, 林姗. 弹性静力学问题的边界积分形式的数值流形法[J]. 岩土力学, 2020, 41(4): 1429-1436.
[6] 杨石扣, 张继勋, 任旭华, . 基于改进数值流形法的接触裂纹问题研究[J]. 岩土力学, 2019, 40(5): 2016-2021.
[7] 王 珍, 曹兰柱, 王 东, . 非均质边坡稳定性上限分析评价研究[J]. 岩土力学, 2019, 40(2): 737-742.
[8] 杨石扣,张继勋,任旭华,. 基于数值流形法的三维裂纹扩展研究[J]. , 2018, 39(S1): 488-494.
[9] 杨石扣,任旭华,张继勋,. 基于数值流形法的重力坝水力劈裂研究[J]. , 2018, 39(8): 3055-3060.
[10] 杨石扣,任旭华,张继勋,. 改进的数值流形法在水力劈裂中的应用[J]. , 2018, 39(10): 3875-3881.
[11] 陈 立,张 朋,郑 宏,. 土石混合体二维细观结构模型的建立与数值流形法模拟[J]. , 2017, 38(8): 2402-2410.
[12] 许晓亮,李建林,龚家伟,宛良朋,陈将宏,. Copula理论下基于g-line失效域的边坡可靠性分析[J]. , 2017, 38(5): 1445-1462.
[13] 刘治军,郑 宏,董 苇,葛修润,孙冠华,. 基于加密物理片的数值流形法中局部网格加密[J]. , 2017, 38(4): 1211-1217.
[14] 张友良,刘登学,刘高敏. 数值流形法的T样条局部加密算法[J]. , 2016, 37(8): 2404-2410.
[15] 徐栋栋,杨永涛,郑 宏,邬爱清,. 基于数值流形法的重力坝多裂纹扩展研究[J]. , 2016, 37(12): 3598-3607.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 魏龙海,王明年,赵东平,吉艳雷. 翔安海底公路隧道陆域段变形控制措施研究[J]. , 2010, 31(2): 577 -581 .
[2] 陈 宇,张庆贺,朱继文,姚海明. 双圆盾构穿越下立交结构的流-固耦合数值模拟[J]. , 2010, 31(6): 1950 -1955 .
[3] 高树生,钱根宝,王 彬,杨作明,刘华勋. 新疆火山岩双重介质气藏供排气机理数值模拟研究[J]. , 2011, 32(1): 276 -280 .
[4] 宋勇军,胡 伟,王德胜,周军林. 基于修正剑桥模型的挤密桩挤土效应分析[J]. , 2011, 32(3): 811 -814 .
[5] 孙德安,孟德林,孙文静,刘月妙. 两种膨润土的土-水特征曲线[J]. , 2011, 32(4): 973 -0978 .
[6] 魏明尧,王恩元,刘晓斐,王 超. 深部煤层卸压爆破防治冲击地压效果的数值模拟研究[J]. , 2011, 32(8): 2539 -2543 .
[7] 褚福永 ,朱俊高 ,贾 华 ,安淑红. 粗粒土卸载-再加载力学特性试验研究[J]. , 2012, 33(4): 1061 -1066 .
[8] 卢 强,王占江,李 进,郭志昀,门朝举. 球面波加载下黄土线黏弹性本构关系[J]. , 2012, 33(11): 3292 -3298 .
[9] 王 苏,路德春,杜修力. 地下结构地震破坏静-动力耦合模拟研究[J]. , 2012, 33(11): 3483 -3488 .
[10] 刘 钦 ,李术才 ,李利平 ,原小帅 ,赵 勇 . 软弱破碎围岩隧道炭质页岩蠕变特性试验研究[J]. , 2012, 33(S2): 21 -0029 .