共查询到20条相似文献，搜索用时 407 毫秒
1.
This paper presents a new numerical tool to model the two-dimensional contaminant transport through saturated porous media
using a meshfree method, called radial point interpolation method (RPIM) with polynomial reproduction. In RPIM, an approximate
solution is constructed entirely in terms of a set of nodes and no characterisation of the interrelationship of the nodes
is needed. The advection–dispersion equation with sorption is considered to illustrate the applicability of the RPIM. The
Galerkin weak form of the governing equation is formulated using 2D meshfree shape functions constructed using thin plate
spline radial basis functions. MATLAB code is developed to obtain the numerical solution. Three numerical examples are presented
and the results are compared with those obtained from the finite element method and analytical solutions. In order to test
the practical applicability and performance of the RPIM, two case studies of contaminant transport through landfill liners
are presented. A good agreement is obtained between the results of the RPIM and the field investigation data. 相似文献
2.
径向基函数点插值无网格法(radial point interpolation method,RPIM)是一种新型的无网格法,其形函数具有插值特性,且形式简单,易于施加本质边界条件。文中介绍了径向基函数点插值无网格法的基本原理,推导了三维情况下点插值无网格法的基本公式。从变分原理出发,结合比奥固结理论,建立了流-固耦合的三维点插值无网格法基本方程和数值积分方法,并开发了相应计算程序。通过三维悬臂梁和单向固结问题的数值试验,验证了该方法对三维弹性问题和流-固耦合问题的适用性和有效性 相似文献
3.
无网格法是一类新型数值算法,具有精度高、高阶形函数构造与物性加载便利等特点,在计算力学领域应用广泛。将无网格方法（PIM、RPIM及EFGM）用于重力异常场二维正演计算:首先从重力异常二维变分问题出发,利用Galerkin法结合高斯积分公式推导了对应的无网格离散系统矩阵表达式;其次通过数值试验得出了RPIM-MQ、RPIM-exp及EFGM-exp形状参数的建议值,最后比较分析了最优形状参数下不同无网格法的计算效果。结果表明:无网格法适用于介质物性分布变化较大的重力异常二维正演,exp函数形状参数αc最优取值区间为[1.5,1.7],β建议值为0.6,MQ函数q取值区间为-4.1~1.9;EFGM较PIM及RPIM具有更高的计算精度。 相似文献
4.
针对线弹性断裂力学问题,提出扩展径向点插值无网格法（X-RPIM）。该方法基于单位分解思想,在传统径向点插值无网格法的位移模式中加入扩展项来描述裂纹两侧的不连续位移场和裂尖奇异场。由于其形函数具有Kronecker ? 函数性质,易于施加本质边界条件。详细描述了X-RPIM不连续位移模式的建立,支配方程的离散形式以及J积分计算混合模式裂纹的应力强度因子的实现过程,讨论了不同积分区域对应力强度因子的影响。数值算例分析证明了该方法在求解断裂问题时的可行性和有效性,同时说明扩展径向点插值无网格法在模拟裂纹扩展问题时具有良好的前景。 相似文献
5.
6.
A fully coupled meshfree algorithm is proposed for numerical analysis of Biot’s formulation. Spatial discretization of the governing equations is presented using the Radial Point Interpolation Method (RPIM). Temporal discretization is achieved based on a novel three-point approximation technique with a variable time step, which has second order accuracy and avoids spurious ripple effects observed in the conventional two-point Crank Nicolson technique. Application of the model is demonstrated using several numerical examples with analytical or semi-analytical solutions. It is shown that the model proposed is effective in simulating the coupled flow deformation behaviour in fluid saturated porous media with good accuracy and stability irrespective of the magnitude of the time step adopted. 相似文献
7.
8.
A meshfree node‐based smoothed point interpolation method (NS‐PIM), which has been recently developed for solid mechanics problems, is applied to obtain certified solutions with bounds for hydraulic structure designs. In this approach, shape functions for displacements are constructed using the point interpolation method (PIM), and the shape functions possess the Kronecker delta property and permit the straightforward enforcement of essential boundary conditions. The generalized smoothed Galerkin weak form is then applied to construct discretized system equations using the node‐based smoothed strains. As a very novel and important property, the approach can obtain the upper bound solution in energy norm for hydraulic structures. A 2D gravity dam problem and a 3D arch dam problem are solved, respectively, using the NS‐PIM and the simulation results of NS‐PIM are found to be the upper bounds. Together with standard fully compatible FEM results as a lower bound, we have successfully determined the solution bounds to certify the accuracy of numerical solutions. This confirms that the NS‐PIM is very useful for producing certified solutions for the analysis of huge hydraulic structures. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
9.
目前隧道及大型地下工程往往在裂隙岩体中开挖,而裂隙与地下空间的距离及裂隙的扩展条件,制约着隧道及地下工程的稳定性。应用能考虑孔洞和裂纹问题的新型边界积分方程与无网格加辽金法结合,建立一种新型的边界无单元法。在该方法中基本的未知量是由边界上的面力和边界上位移密度函数构成的复变量边界函数 。文中应用的边界积分公式和Muskhelishvili的积分公式直接相关。将无网格构造方法引入新型的边界积分方程,建立了新型的边界无单元法。应用该方法详细分析了含隧道和裂纹间相互关系等问题,其数值结果与解析结果吻合很好,说明该方法的正确性和可行性。 相似文献
10.
无网格法是近几年来发展的一种新的基于变分原理的数值计算方法,由于在计算形函数中不需要划分网格.在力学、电磁学等领域得到了广泛的研究.基于无网格法在大地电磁勘探正演中的应用进行了研究,首先对无网格法的基本原理进行了阐述,并利用广义变分原理推导出了相应的离散方程,编制了相应的程序,最后通过两个理论模型的计算结果检验了算法的正确性. 相似文献
11.
假定各物理量沿径向不变化,建立了低应变瞬态集中荷载作用下现浇大直径管桩振动响应的计算模型和波动方程。采用Laplace变换法,求得了波动方程的频域解析解,采用Fourier逆变换求得了时域响应。将文中解计算结果与三维频域解析解进行了对比分析,文中解的入射峰-反射峰时间差与三维频域解差别很小,说明对于PCC桩这种大直径薄壁管桩,采用径向不变假定对计算结果几乎没有影响。因此,在PCC桩低应变检测波形分析时,采用基于径向不变假定的二维解是完全合理的。将文中二维频域解的计算结果与二维时域解的结果进行了对比分析,结果表明,2种解在入射波、第一个桩底反射波、第二个桩底反射波峰值大小和到达时间非常吻合,这说明考虑环向位移与否对计算结果没有太大的影响,研究PCC桩低应变动力响应的问题时忽略水平方向的位移是完全可以的。 相似文献
12.
13.
14.
An alternative approach for quasi‐static large deformation analysis of saturated porous media using meshfree method 
下载免费PDF全文
下载免费PDF全文 An alternative coupled large deformation formulation combined with a meshfree approach is proposed for flow–deformation analysis of saturated porous media. The formulation proposed is based on the Updated Lagrangian (UL) approach, except that the spatial derivatives are defined with respect to the configuration of the medium at the last time step rather than the configuration at the last iteration. In this way, the Cauchy stresses are calculated directly, rendering the second Piola–Kirchhoff stress tensor not necessary for the numerical solution of the equilibrium equations. Moreover, in contrast with the UL approach, the nodal shape function derivatives are calculated once in each time step and stored for use in subsequent iterations, which reduces the computational cost of the algorithm. Stress objectivity is satisfied using the Jaumann stress rate, and the spatial discretisation of the governing equations is achieved using the standard Galerkin method. The equations of equilibrium are satisfied directly, and the nonlinear parts of the system matrix are derived independent of the stresses of the medium resulting in a stable numerical algorithm. Temporal discretisation is effected based on a three‐point approximation technique that avoids spurious ripple effects and has second‐order accuracy. The radial point interpolation method is used to construct the shape functions. The application of the formulation and the significance of large deformation effects on the numerical results are demonstrated through several numerical examples. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
15.
16.
Modelling of contaminant transport through landfill liners and natural soil deposits is an important area of research activity in geoenvironmental engineering. Conventional mesh‐based numerical methods depend on mesh/grid size and element connectivity and possess some difficulties when dealing with advection‐dominant transport problems. In the present investigation, an attempt has been made to provide a simple but sufficiently accurate methodology for numerical simulation of the two‐dimensional contaminant transport through the saturated homogeneous porous media and landfill liners using element‐free Galerkin method (EFGM). In the EFGM, an approximate solution is constructed entirely in terms of a set of nodes and no characterization of the interrelationship of the nodes is needed. The EFGM employs moving least‐square approximants to approximate the function and uses the Lagrange multiplier method for imposing essential boundary conditions. The results of the EFGM are validated using experimental results. Analytical and finite element solutions are also used to compare the results of the EFGM. In order to test the practical applicability and performance of the EFGM, three case studies of contaminant transport through the landfill liners are presented. A good agreement is obtained between the results of the EFGM and the field investigation data. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
17.
18.
This study introduces the prediction of probabilistic settlements with the uncertainty in the spatial variability of Young’s modulus to illustrate the preliminary development of a spectral stochastic meshless local Petrov–Galerkin (SSMLPG) method. Generalized polynomial chaos expansions of Young’s moduli and a two-dimensional meshfree weak–strong formulation in elasticity are combined to derive the SSMLPG formulation. Because of the local and truly meshless nature, the SSMLPG method is more computationally efficient than available stochastic numerical methods. Two examples further show that SSMLPG-based predictions remain sufficiently accurate even in case of scattered nodes. Therefore, the SSMLPG method can be a valuable alternative for solving stochastic boundary-value problems. 相似文献
19.
无网格法在地下水水位预测中的应用 总被引：1，自引：0，他引：1
基于移动最小二乘理论的无网格法是近几年来兴起的一种新的数值计算方法,与有限元法相比,它的主要优点在于无需单元信息,只需节点信息。用无网格法构造了场函数,包括基函数和权函数的选取,形函数及其导数的计算。根据鞍山市首山区水文地质条件,建立了求解双层渗流二维平面系统的数学模型,详细推导了模型求解的无网格方程。应用已识别的参数,用无网格法对该数学模型进行了求解, 并对该区的地下水水位进行了预测,预测的水位与实际水位变化规律基本一致。 相似文献
20.
Three‐dimensional simulations of tensile cracks in geomaterials by coupling meshless and finite element method 下载免费PDF全文
下载免费PDF全文 Failure in geotechnical engineering is often related to tension‐induced cracking in geomaterials. In this paper, a coupled meshless method and FEM is developed to analyze the problem of three‐dimensional cracking. The radial point interpolation method (RPIM) is used to model cracks in the smeared crack framework with an isotropic damage model. The identification of the meshless region is based on the stress state computed by FEM, and the adaptive coupling of RPIM and FEM is achieved by a direct algorithm. Mesh‐bias dependency, which poses difficulties in FEM‐based cracking simulations, is circumvented by a crack tracking algorithm. The performance of our scheme is demonstrated by two numerical examples, that is, the four‐point bending test on concrete beam and the surface cracks caused by tunnel excavation. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献