›› 2016, Vol. 37 ›› Issue (4): 1187-1194.doi: 10.16285/j.rsm.2016.04.034

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

无穷域问题的弱形式求积元分析

袁 帅,钟宏志   

  1. 清华大学 土木工程系,北京 100084
  • 收稿日期:2014-07-14 出版日期:2016-04-11 发布日期:2018-06-09
  • 作者简介:袁帅,男,1984年生,博士研究生,主要从事土力学数值方法方面的研究工作。
  • 基金资助:

    国家自然科学基金(No. 51178247)。

Analysis of unbounded domain problems by the weak form quadrature element method

YUAN Shuai,ZHONG Hong-zhi   

  1. Department of Civil Engineering, Tsinghua University, Beijing 100084, China
  • Received:2014-07-14 Online:2016-04-11 Published:2018-06-09
  • Supported by:

    This work was supported by the National Natural Science Foundation of China (NSFC) (51178247).

摘要: 岩土工程中经常会遇到无穷域问题,而采用无限单元可以实现对其有效地模拟。弱形式求积元法是一个有效的数值工具,它常通过提高积分阶次来提高计算精度。建立了无限弱形式求积单元并被应用于求解岩土工程中的无穷域问题,该单元基于坐标映射,将无穷域变换到标准域,在标准域上进行数值积分和数值微分,保留了传统弱形式求积元的积分点坐标和权系数。求解了瞬态渗流、固结和静力分析等数值算例,并与解析解或截断方法进行了对比。结果表明:基于坐标映射的无限弱形式求积单元使用简单,可以模拟各种类型的无穷域问题,仅需要将感兴趣的范围进行有限域划分并通过提高积分阶次来减小对极点位置的依赖,极大地节省了计算资源,提高了计算精度。

关键词: 无限弱形式求积单元, 无穷域, 静力分析, 渗流, 固结, 坐标映射

Abstract: Unbounded domain problems are frequently encountered in geotechnical engineering and infinite elements are often used effectively for simulation. The weak form quadrature element method is an effective numerical tool in which the computational accuracy is often improved through increasing the order of integration. An infinite weak form quadrature element is developed and applied to the analysis of unbounded domain problems in geotechnical engineering. Based on coordinate transformation, an unbounded domain is mapped onto a standard region where numerical integration and numerical differentiation are conducted and the conventional numerical integration points and weights in the weak form quadrature element method are retained. Numerical examples in the areas of transient seepage, consolidation and elastostatic analysis are given and the results are compared with analytical solutions or those of other numerical methods. It is shown that the infinite weak form quadrature element is simple and applicable to solution of various unbounded domain problems, while conventional elements are used during discretization of the domain of interest. In addition, the dependence on the pole of coordinate transformation can be alleviated considerably with the increase of the integration order of the element. Consequently, computational resources are reduced significantly and accuracy of results is enhanced.

Key words: infinite weak form quadrature element, unbounded domain, static analysis, seepage, consolidation, coordinate transformation

中图分类号: 

  • TU 43

[1] 王珂, 盛金昌, 郜会彩, 田晓丹, 詹美礼, 罗玉龙, . 应力−渗流侵蚀耦合作用下粗糙裂隙渗流特性研究[J]. 岩土力学, 2020, 41(S1): 30-40.
[2] 安笑, 潘华利, 欧国强, 孔玲, 李炳志, . 恒定渗流作用泥石流碎屑物质起动判别模型研究[J]. 岩土力学, 2020, 41(S1): 115-122.
[3] 庄心善, 赵汉文, 王俊翔, 黄勇杰, 胡智. 循环荷载下重塑弱膨胀土滞回曲线 形态特征定量研究[J]. 岩土力学, 2020, 41(6): 1845-1854.
[4] 胡盛斌, 杜国平, 徐国元, 周天忠, 钟有信, 石重庆, . 基于能量测量的声呐渗流矢量法及其应用[J]. 岩土力学, 2020, 41(6): 2143-2154.
[5] 孙德安, 薛垚, 汪磊, . 变荷载作用下考虑半透水边界热传导性的 一维饱和土热固结特性研究[J]. 岩土力学, 2020, 41(5): 1465-1473.
[6] 李红坡, 梅国雄, 肖涛, 陈征. 涂抹区重叠竖井地基固结特性研究[J]. 岩土力学, 2020, 41(5): 1560-1566.
[7] 江留慧, 李传勋, 杨怡青, 张锐. 变荷载下双层地基一维非线性固结近似解析解[J]. 岩土力学, 2020, 41(5): 1583-1590.
[8] 陈琼, 崔德山, 王菁莪, 刘清秉. 不同固结状态下黄土坡滑坡滑 带土的蠕变试验研究[J]. 岩土力学, 2020, 41(5): 1635-1642.
[9] 杨赫, 程卫民, 刘震, 王文玉, 赵大伟, 王文迪. 注水煤体有效渗流通道结构分形特征 核磁共振试验研究[J]. 岩土力学, 2020, 41(4): 1279-1286.
[10] 师旭超, 孙运德. 线性卸荷作用下软土超孔隙水压力 变化规律分析[J]. 岩土力学, 2020, 41(4): 1333-1338.
[11] 任宇晓, 闫玥, 付登锋. 浅层地基上管道轴向运动的阻力研究[J]. 岩土力学, 2020, 41(4): 1404-1411.
[12] 刘建民, 邱月, 郭婷婷, 宋文智, 谷川, . 饱和粉质黏土静剪强度与振动后 静剪强度对比研究[J]. 岩土力学, 2020, 41(3): 773-780.
[13] 郎瑞卿, 杨爱武, 闫澍旺, . 修正等应变假定下刚性桩复合地基固结特性分析[J]. 岩土力学, 2020, 41(3): 813-822.
[14] 米博, 项彦勇, . 砂土地层浅埋盾构隧道开挖渗流稳定性的 模型试验和计算研究[J]. 岩土力学, 2020, 41(3): 837-848.
[15] 杨高升, 白冰, 姚晓亮, . 高含冰量冻土路基融化固结规律研究[J]. 岩土力学, 2020, 41(3): 1010-1018.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!