数值分析

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

展开
  • 清华大学 土木工程系,北京 100084
袁帅,男,1984年生,博士研究生,主要从事土力学数值方法方面的研究工作。

收稿日期: 2014-07-14

  网络出版日期: 2018-06-09

基金资助

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

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

Expand
  • Department of Civil Engineering, Tsinghua University, Beijing 100084, China

Received date: 2014-07-14

  Online published: 2018-06-09

Supported by

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

摘要

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

本文引用格式

袁 帅,钟宏志 . 无穷域问题的弱形式求积元分析[J]. 岩土力学, 2016 , 37(4) : 1187 -1194 . DOI: 10.16285/j.rsm.2016.04.034

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.
文章导航

/