共查询到20条相似文献，搜索用时 203 毫秒
1.
岩土工程数值计算中的无网格方法及其全自动布点技术 总被引：14，自引：3，他引：11
自然单元法采用无网格的思想全域构造插值函数,它的求解精度高,计算时间少,可准确地施加边界条件,兼具有无网格法和有限单元法的优点和特点,是一种理想的用于岩土及地下工程分析计算的数值方法。文中简要地介绍了自然单元法的基本理论,并针对岩土及地下工程问题特点,给出了一种无网格离散点的全自动布置方法。 相似文献
2.
《物探化探计算技术》2015,(3)
无单元Galerkin法作为较成熟的一种无网格方法,已成功应用于有限元法触及的领域,还解决了如大变形、裂纹扩展及高速冲击等网格方法较难处理的问题,但其最大的缺陷在于系统方程的离散需借助背景网格,因此该方法并非真正意义上的无网格方法。无网格局部径向基点插值法采用子域法构造系统方程,加权残量只要求在局部积分域消除,大大降低了对背景网格的依赖,向真正的无网格迈进了一大步.这里将此方法用于大地电磁二维正演,介绍了该方法的基本原理;从大地电磁二维边值问题出发,利用子域法推导了与之对应的无网格局部弱式系统方程,并用高斯积分将其离散化;论述了局部径向基点插值法较无单元Galerkin法及有限元法的优缺点;最后通过二维模型的计算验证了算法的有效性。 相似文献
3.
本文提出用二层介质格林函数作为余量加权函数,归化边界积分方程的边界单元法。通过对有水平覆盖层的点源二维、点源三维正演计算表明,本算法明显地优化了以往的边界元法,计算精度高,方程组规模小,节省内存空间,仅在一般微机上即可实现,扩大了边界单元法在电法勘探中的应用范围。 相似文献
4.
5.
提出了一种模拟裂纹扩展的水平集和无网格耦合方法。由于水平集和无网格方法都是建立在离散节点上,因而可以很自然地实现耦合。在该方法中,两个在裂尖处相互正交的水平集不仅用于描述裂纹的几何形态和裂尖位置,而且用于建立无网格伽辽金法（简称EFGM）不连续近似函数中的Heaviside跳跃项和裂尖处的Westergaard扩展项。当裂纹扩展时,则由水平集更新算法确定新裂纹的位置。水平集和无网格耦合法无需使用可视法、衍射法或透明法,克服了这些方法在裂尖处人为引入的不连续且能很好地再生 奇异场;而且节点影响域不受裂纹线切割的影响,在计算中往往使用较小的影响域,保持了整体刚度矩阵的带状、稀疏性;另外,水平集简化了扩展节点的选取和附加函数的建立,其更新过程无需求解演化方程,实现简单且易于编程。数值算例表明本文方法具有较高的计算精度,其模拟的裂纹扩展路径与试验结果吻合得很好,从而验证了本文方法的正确性和可行性。 相似文献
6.
7.
利用格林公式及点源场的边界条件,将地下半空间位函数的体积积分变为地面位函教的积分方程以后,用边界元方法解此积分方程便可得到地面各节点的电位值。其计算精度和速度取决于剖分方式和节点密度。本文采用变网格的剖分方式,用较少的点实现了计算,并兼顾了精度和速度。 相似文献
8.
无单元法（element free galerkin method）是一种基于滑动最小二乘法构造插值函数的无网格数值方法。以该方法中插值函数及其导函数为基础,从变分原理出发,建立了弹性地基上正交各向异性Mindlin板的无单元法求解控制方程,其中弹性地基采用Winkler地基模型。对该理论方法进行了相应的程序实现,最后通过算例分析表明了提出方法用于水工结构中弹性地基板问题分析的可行性和有效性。 相似文献
9.
针对溃坝水流数值模拟面临的复杂地形和不规则边界等问题,基于结构网格建立了适应复杂地形和不规则边界的溃坝水流数值模拟有限体积模型（HydroM2D）。模型基于具有守恒特性的二维浅水方程,利用HLLC格式的近似Riemann解计算网格界面通量,利用MUSCL-Hancock法不断向前积分,使模型在时空上具有二阶精度;对源项进行离散处理确保模型的稳定性;模型引入有效干湿边界和不规则地形边界处理方法,准确模拟了干湿单元的动态交替和复杂边界上的水流特性。最后分别利用水槽试验、物理模型和实际算例对模型进行验证。结果表明,该模型对不同情景下的溃坝洪水模拟结果和实测资料以及现有模型模拟结果具有较高的一致性,模拟精度较高,稳定性较好,具有推广应用价值。 相似文献
10.
裂隙岩体因含有发育程度不同的裂隙、节理和断层等不连续面,致其渗透性具有各向异性、不连续性等特点,因此传统的有限元法对分布密集的裂隙岩体渗流场求解有一定的难度。本文提出了采用无单元Glaerkin法求解有自由面裂隙渗流问题,并推导了无单元法求解渗流场的基本方程和积分格式,给出了应用罚函数法处理渗流边界条件和自由面处理方法。采用IDL语言编制了二维无单元法计算软件LIDAREFM。文中以北京怀柔桥梓镇某裂隙岩体边坡渗流场计算为例,研究了复杂裂隙共同作用下渗流场特性和自由面分布,讨论了不同开度、不同连通程度的裂隙对渗流场的影响。研究结果表明:无单元法可以较好地解决有密集裂隙的岩体渗流场的求解问题,实现了裂隙处结点任意加密以及积分网格的独立布置,避免了对有自由面和裂隙穿越的子域的重新处理,简化了渗流问题的求解过程。 相似文献
11.
This paper presents a single‐domain boundary element method (BEM) for linear elastic fracture mechanics analysis in the two‐dimensional anisotropic material. In this formulation, the displacement integral equation is collocated on the un‐cracked boundary only, and the traction integral equation is collocated on one side of the crack surface only. A special crack‐tip element was introduced to capture exactly the crack‐tip behavior. A computer program with the FORTRAN language has been developed to effectively calculate the stress intensity factors of an anisotropic material. This BEM program has been verified having a good accuracy with the previous researches. Furthermore, by analyzing the different anisotropic degree cracks in a finite plate, we found that the stress intensity factors of crack tips had apparent influence by the geometry forms of cracks and media with different anisotropic degrees. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
12.
O. Damisa 《国际地质力学数值与分析法杂志》1986,10(1):73-89
A problem common to most numerical methods of analysing the vibration of a rigid body on an elastic half space is indentified especially as regards the arbitrary truncation of the semi-infinite half space or the numerically expensive computations involved in some of the methods. This paper presents a boundary integral equation (BIE) method of defining a finite domain to replace the half space and consists of evaluating the BIE over the domain. A mathematical justification is provided for defining the domain, and the error involved in such truncation of the half space in shown to be small. Yet the computations are very cheap owing to the simple trigonometric/algebraic functions involved in the BIE. The accuracy of the method is quite commendable especially as seen from the results presented for the circular foundation for which analytical and experimental solutions are available for comparison. 相似文献
13.
14.
一种新的二维电阻率模拟方法 总被引：2，自引：0，他引：2
目前，对地下介质作电阻率成像需要将研究区域网格化。本文以格林函数为数学工具导出了网格化地电模型二维问题的边界积分方程，该方程适合于电阻率成像。用边界单元法求解该方程可以得出空间电位分布，从而进行视电阻率模拟。计算与实验结果的对比表明，这种方法是正确、可行的。 相似文献
15.
This paper presents a boundary element method (BEM) procedure for a linear elastic fracture mechanics analysis in two‐dimensional anisotropic bimaterials. In this formulation, a displacement integral equation is only collocated on the uncracked boundary, and a traction integral equation is only collocated on one side of the crack surface. A fundamental solution (Green's function) for anisotropic bimaterials is also derived and implemented into the boundary integral formulation so that except for the interfacial crack part, the discretization along the interface can be avoided. A special crack‐tip element is introduced to capture the exact crack‐tip behavior. A computer program using FORTRAN has been developed to effectively calculate the stress intensity factors of an anisotropic bimaterial. This BEM program has been verified to have a good accuracy with previous studies. In addition, a central cracked bimaterial Brazilian specimen constituting cement and gypsum is prepared to conduct the Brazilian test under diametral loading. The result shows that the numerical analysis can predict relatively well the direction of crack initiation and the path of crack propagation. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
16.
将无网格局部Petrov-Galerkin算法用于大地电磁二维正演。介绍了该方法的基本原理;从大地电磁二维边值问题出发,利用子域法详细推导了与之对应的局部Petrov-Galerkin弱式方程,并用高斯积分法将其离散化。论述了无网格局部Petrov-Galerkin法较无单元Galerkin法及有限元法的优缺点,最后通过二维模型的计算验证了算法的有效性。 相似文献
17.
A fully coupled element‐free Galerkin model for hydro‐mechanical analysis of advancement of fluid‐driven fractures in porous media 
下载免费PDF全文
下载免费PDF全文 Hydraulic fracturing (HF) of underground formations has widely been used in different fields of engineering. Despite the technological advances in techniques of in situ HF, the industry uses semi‐analytical tools to design HF treatment. This is due to the complex interaction among various mechanisms involved in this process, so that for thorough simulations of HF operations a fully coupled numerical model is required. In this study, using element‐free Galerkin (EFG) mesh‐less method, a new formulation for numerical modeling of hydraulic fracture propagation in porous media is developed. This numerical approach, which is based on the simultaneous solution of equilibrium and continuity equations, considers the hydro‐mechanical coupling between the crack and its surrounding porous medium. Therefore, the developed EFG model is capable of simulating fluid leak‐off and fluid lag phenomena. To create the discrete equation system, the Galerkin technique is applied, and the essential boundary conditions are imposed via penalty method. Then, the resultant constrained integral equations are discretized in space using EFG shape functions. For temporal discretization, a fully implicit scheme is employed. The final set of algebraic equations that forms a non‐linear equation system is solved using the direct iterative procedure. Modeling of cracks is performed on the basis of linear elastic fracture mechanics, and for this purpose, the so‐called diffraction method is employed. For verification of the model, a number of problems are solved. According to the obtained results, the developed EFG computer program can successfully be applied for simulating the complex process of hydraulic fracture propagation in porous media. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
18.
针对线弹性断裂力学问题,提出扩展径向点插值无网格法（X-RPIM）。该方法基于单位分解思想,在传统径向点插值无网格法的位移模式中加入扩展项来描述裂纹两侧的不连续位移场和裂尖奇异场。由于其形函数具有Kronecker ? 函数性质,易于施加本质边界条件。详细描述了X-RPIM不连续位移模式的建立,支配方程的离散形式以及J积分计算混合模式裂纹的应力强度因子的实现过程,讨论了不同积分区域对应力强度因子的影响。数值算例分析证明了该方法在求解断裂问题时的可行性和有效性,同时说明扩展径向点插值无网格法在模拟裂纹扩展问题时具有良好的前景。 相似文献
19.
This paper presents an indirect boundary integral equation method for analysis of quasi-static, time-harmonic and transient boundary value problems related to infinite and semi-infinite poroelastic domains. The present analysis is based on Biot's theory for poroelastodynamics with fluid viscous dissipation. The solution to a given boundary value problem is reduced to the determination of intensities of forces and fluid sources applied on an auxiliary surface defined interior to the surface on which the boundary conditions are specified. A coupled set of integral equations is established to determine the intensities of forces and fluid sources applied on the auxiliary surface. The integral equations are solved numerically in the Laplace domain for quasi-static and transient problems, and in the frequency domain for time-harmonic excitations. The kernel functions of the integral equation correspond to appropriate Green's functions for a poroelastic full space or half-space. The convergence and numerical stability of the present scheme are established by considering a number of bench mark problems. The versatility of the present method is demonstrated by studying the quasi-static response of a rigid spheroidal anchor, and time-harmonic and transient response of a rigid semi-circular tunnel. 相似文献
20.
In this paper, the superimposing principle and the theory of Mellin transforms are used to reduce the controlling equations and boundary conditions of linear cracks intersecting hole into dual integral equations in an infinite elastic body. With finite Hilbert transforms, the dual integral equations can be reduced into a Fredhole integral equation of the second kind which is related to S(t2 ). The dynamic stress intensity factor can be expressed with the solution of integral equation and its numerical results in obtained. In the particular case, the numerical results are the same as those in reference[2] . the intensity factor decreases with the propagation of cracks. With experiments, the propagation velocity of these linear cracks are measured. 相似文献