共查询到16条相似文献,搜索用时 129 毫秒
1.
径向基函数点插值无网格法(radial point interpolation method,RPIM)是一种新型的无网格法,其形函数具有插值特性,且形式简单,易于施加本质边界条件。文中介绍了径向基函数点插值无网格法的基本原理,推导了三维情况下点插值无网格法的基本公式。从变分原理出发,结合比奥固结理论,建立了流-固耦合的三维点插值无网格法基本方程和数值积分方法,并开发了相应计算程序。通过三维悬臂梁和单向固结问题的数值试验,验证了该方法对三维弹性问题和流-固耦合问题的适用性和有效性 相似文献
2.
岩土体的渗透破坏、地下工程的防渗设计等无不与渗流计算有关。针对渗流自由面问题,提出一种重心拉格朗日插值的配点型无网格方法。由于渗流自由面问题的求解区域是不规则区域,该方法通过将不规则求解区域嵌入一个正则矩形区域,在正则区域上采用重心拉格朗日插值近似未知函数,利用配点法离散渗流问题的控制方程,将重心拉格朗日插值的微分矩阵离散成代数方程表达的矩阵形式。将自由面上的边界条件通过重心拉格朗日插值离散,通过置换方程法和附加方程法施加边界条件,利用正则区域上的重心插值配点法,通过迭代确定最终自由面的位置。数值算例表明所提出的无网格方法对于求解渗流自由面问题的正确性和高精度。 相似文献
3.
针对线弹性断裂力学问题,提出扩展径向点插值无网格法(X-RPIM)。该方法基于单位分解思想,在传统径向点插值无网格法的位移模式中加入扩展项来描述裂纹两侧的不连续位移场和裂尖奇异场。由于其形函数具有Kronecker ? 函数性质,易于施加本质边界条件。详细描述了X-RPIM不连续位移模式的建立,支配方程的离散形式以及J积分计算混合模式裂纹的应力强度因子的实现过程,讨论了不同积分区域对应力强度因子的影响。数值算例分析证明了该方法在求解断裂问题时的可行性和有效性,同时说明扩展径向点插值无网格法在模拟裂纹扩展问题时具有良好的前景。 相似文献
4.
本文中使用的径向基函数配点法是以时空配点法为基础来解决抛物型方程的一类问题。这种方法与近似求时间导数的隐式,显式法以及其他数值法不同,它不需要对离散系统的时间稳定性进行分析。用时空径向基函数配点法求解二维地下水非稳定流动问题,通过呈现有混合边界条件及只有一类边界条件两种情况下的计算结果,说明了该方法求解该问题的精度及效率较高,结果理想。 相似文献
5.
6.
介绍径向基函数插值配点法,将其应用于非均质多孔介质中的一维地下水稳定流、非稳定流问题,算例结果表明,该方法既计算效率高又有较高的精度。 相似文献
7.
8.
9.
10.
11.
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. 相似文献
12.
Pankaj K. Mishra Sankar K. Nath Mrinal K. Sen Gregory E. Fasshauer 《Computational Geosciences》2018,22(5):1203-1218
Scattered data interpolation schemes using kriging and radial basis functions (RBFs) have the advantage of being meshless and dimensional independent; however, for the datasets having insufficient observations, RBFs have the advantage over geostatistical methods as the latter requires variogram study and statistical expertise. Moreover, RBFs can be used for scattered data interpolation with very good convergence, which makes them desirable for shape function interpolation in meshless methods for numerical solution of partial differential equations. For interpolation of large datasets, however, RBFs in their usual form, lead to solving an ill-conditioned system of equations, for which, a small error in the data can cause a significantly large error in the interpolated solution. In order to reduce this limitation, we propose a hybrid kernel by using the conventional Gaussian and a shape parameter independent cubic kernel. Global particle swarm optimization method has been used to analyze the optimal values of the shape parameter as well as the weight coefficients controlling the Gaussian and the cubic part in the hybridization. Through a series of numerical tests, we demonstrate that such hybridization stabilizes the interpolation scheme by yielding a far superior implementation compared to those obtained by using only the Gaussian or cubic kernels. The proposed kernel maintains the accuracy and stability at small shape parameter as well as relatively large degrees of freedom, which exhibit its potential for scattered data interpolation and intrigues its application in global as well as local meshless methods for numerical solution of PDEs. 相似文献
13.
14.
A new numerical tool is presented which models the two-dimensional contaminant transport through saturated porous media using a meshfree method called the 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. An 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 two-dimensional meshfree shape functions constructed using thin plate spline radial basis functions. A computer program is developed for the implementation of the RPIM procedure. Three numerical examples are presented and the results are compared with those obtained from the analytical solution and finite element method. The experimental results are also used to validate the approach. The proposed RPIM has generated results with no oscillations and they are insensitive to Peclet constraints. 相似文献
15.
A novel procedure associated with the precise integration method (PIM) and the technique of dual vector is proposed to effectively calculate the magnitude and distribution of deformations in a homogeneous multilayered transversely isotropic medium. The planes of transverse isotropy are assumed to be parallel to the horizontal surface of the soil system. The linearly elastic medium is subjected to four types of vertically acting axisymmetric loads prescribed either at the external surface or in the interior of the soil medium. There are no limits for the thicknesses and number of soil layers to be considered. By virtue of the governing equations of motion and the constitutive equations of the transversely isotropic elastic body, and based on the Hankel integral transform and a dual vector formulation in a cylindrical coordinate system, the partial differential motion equations can be converted into first‐order ordinary differential matrix equations. Applying the approach of PIM, it is convenient to obtain the solutions of ordinary differential matrix equations for the continuously homogeneous multilayered transversely isotropic elastic soil in the transformed domain. The PIM is a highly accurate algorithm to solve the sets of first‐order ordinary differential equations, which can ensure to achieve any desired accuracy of the solutions. What is more, all calculations are based on the standard method with the corresponding algebraic operations. Computational efforts can be reduced to a great extent. Finally, numerical examples are provided to illustrate the accuracy and effectiveness of the proposed approach. Some more cases are analyzed to evaluate the influences of the elastic parameters of the transversely isotropic media on the load‐displacement responses. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
16.
为了研究采样和网格化方法对地球物理数据成图精度的影响,为野外数据采集布设提供一定的依据,采用数值模拟确定重力异常场场值,通过不同采样间距和不同插值方法计算重力异常绝对误差均方根值和节点处的绝对误差值,对比不同插值方法的误差,得到了如下认识:1)对于同一插值方法而言,存在小间距绝对误差均方根值小于大间距绝对误差均方根值的关系。2)对不同的插值方法而言:当采样间距小于最小异常地质体尺度时,绝对误差均方根值由小到大的顺序是径向基函数法、改进的谢别德法、克里金插值法、自然邻点法、反距离加权插值法、最近邻点法、最小曲率法,并且线性插值三角网法与自然邻点法具有几乎相同的数值;当采样间距大于最小异常地质体尺度时,绝对误差均方根值由小到大的顺序是径向基函数法、改进的谢别德法、克里金插值法、自然邻点法、最小曲率法、最近邻点法、反距离加权插值法,并且线性插值三角网法和自然邻点法具有几乎相同的数值。3)从绝对误差均方值看,径向基函数方法、改进的谢别德方法和克里金方法数值较小,其中径向基函数值绝对误差均方根值最小。4)从节点处绝对误差值来看,径向基函数方法、克里金方法、改进的谢别德方法相对其他插值方法具有更小的误差,不存在局部误差较小或较大的情况,是相对较好的插值方法,并且径向基函数方法是最好的。 相似文献