共查询到20条相似文献，搜索用时 109 毫秒
1.
本给出了用边界元法积分方程的算法。通过处例验算表明，这种算法综合了边界元法精度高和用矩最法解积分方程速度快的优点，是一种求解积分方程十分有效的算法。 相似文献
2.
本文给出了用边界元法求积分方程的算法。通过算例验算表明，这种算法综合了边界元法精度高和用矩量法解积分方程速度快的优点，是一种求解积分方程十分有效的算法。 相似文献
3.
4.
一种新的二维电阻率模拟方法 总被引：2，自引：0，他引：2
目前，对地下介质作电阻率成像需要将研究区域网格化。本文以格林函数为数学工具导出了网格化地电模型二维问题的边界积分方程，该方程适合于电阻率成像。用边界单元法求解该方程可以得出空间电位分布，从而进行视电阻率模拟。计算与实验结果的对比表明，这种方法是正确、可行的。 相似文献
5.
结合边界元法和离散元法,提出一种可以进行计算颗粒内部应力和破碎路径的方法。该方法利用离散元法求解颗粒的相互作用和每个颗粒上的荷载。然后利用边界元法计算颗粒的应力分布,为了实现动态平衡,将颗粒的加速度视为恒定大小的体力。但体力导致边界积分方程中出现域积分,故采用直线积分法将域积分转化为边界积分,以保证边界元法降维的优势。为了提高边界元的计算效率,对于几何形状相似的颗粒,以其中一个颗粒作为模板颗粒,只需要计算模板颗粒在局部坐标系中的系数矩阵,其他相似颗粒可以通过局部和全局坐标系之间的映射获得。在得到应力后,基于Hoek-Brown准则来判断颗粒是否破碎。此外,将破坏路径简化为直线,并采用最小二乘法拟合得到破坏路径。 相似文献
6.
针对传统弹性理论法过高地估计了桩-桩相互作用效应,基于杆件有限单元法建立了半无限土体中群桩基础的桩侧剪应力求解方程,通过简化桩-桩的相互作用效应,将群桩基础桩侧摩阻力的求解方程近似解耦,实现了群桩基础桩身剪应力和位移的快速求解目标。通过两桩相互作用系数以及柔性承台下群桩基础差异沉降的对比分析,验证了该简化方法的合理性。参数分析结果表明,该简化方法计算得到的相互作用系数与严格的边界积分方程法解答较为接近,稍小于Poulos弹性理论法的计算结果;柔性承台下群桩的差异沉降在桩间距较小时与经典解答较为接近,而在桩间距较大时则存在一定的差别。该简化方法大幅减少了群桩计算工作量,适用于大规模群桩基础的快速计算要求。 相似文献
7.
本文讨论了在均匀谐变场激发下导电导磁体电磁响应的边界单元法计算问题。文中首先建立矢量位的边值问题,然后导出其边界积分方程,最后用边界元法进行求解。本方法的计算结果具有精度高、速度快的优点。 相似文献
8.
9.
第二类非线性Fredholm型积分方程数值解 总被引：1，自引：0，他引：1
配置法研究了地球物理中常见的第二类非线性 Fredholm 型积分方程的数值解法,将第二类非线性 Fredholm 型积分方程转化为非线性代数方程组进行求解,采用高斯数值积分公式,给出了数值计算的具体实例.利用Matlab软件的符号运算功能编程计算,克服了非线性方程难于变成求解的困难,数值例子表明该方法编程简便有效.对非线性积分方程和非线性代数方程组的求解都有重要价值. 相似文献
10.
扶永铭 《成都理工学院学报》1995,22(3):109-114
该文研究了低频人工时变场（似稳场）中任意导电导磁矿体的边界元素法。通过对矢量位的边界积分方程，再以样条边界元素法进行离散，得到实分量的常系数方程组和虚分量的二阶微分方程组。此方程可利用时间积分的样条函数方法求解，解出适量位后，利用磁场强度与矢量位Ａ的相应关系而求得二次磁场的水平分量和垂直分量。 相似文献
11.
Abstract The calculations of unsteady flow to a multiple well system with the application of boundary element method (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem of parabolic differential equation. It is changed into an elliptic one by Laplace transform to eliminate time variable. The image function of water head H can be solved by BEM. We derived the boundary integral equation of the transformed variable H and the discretization form of it, so that there is no need to discretize the boundaries of well walls and it becomes easier to solve the groundwater head H by numerical inversion. 相似文献
12.
Boundary integral equations for traction boundary-value problems of two-dimensional elastostatics are derived by the indirect boundary element method. Quadratic variation functions for the representation of geometry, fictitious forces and displacements over each boundary element are described. A system of equations approximating to the boundary integral equations is obtained by a Galerkin formulation in which the integral equation is written at Gauss integration points of elements. The method of computation of the Cauchy principal value is described. Examples of application to the analysis of stress and displacement around underground excavations demonstrate the accuracy and efficiency of the formulation. 相似文献
13.
C. Fidelibus 《国际地质力学数值与分析法杂志》2007,31(11):1329-1348
This paper describes the essential features of a numerical technique for the simulation of the coupled fluid flow and deformation in a 2D assembly of poroelastic blocks and transmissive fractures. The boundary element method (BEM) is applied to each block to reduce Navier and diffusion equations to a set of integral equations involving block boundary terms, whereas a Galerkin weighted‐residuals finite element method (FEM) is applied to the fracture diffusion equations. In addition, fracture local equilibrium is rendered through spring‐like equations relating the stresses to the relative displacements of the fracture walls. A time‐marching process is implemented leading to an algebraic system where the right‐hand side vector is built based on the collected solutions of the previous time steps. The technique requires the meshing of the fracture network only. The accuracy of the results is adequate even with relatively coarse meshes without the resort to small time steps at the beginning of the simulation. It furnishes outputs that focus only on the salient features of the response. The efficiency of the technique is demonstrated through the illustration of the results of three examples. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
14.
A time domain boundary element method (BEM) for evaluating stresses in an axisymmetric soil mass undergoing consolidation has been developed. Previous BEM work on axisymmetric poroelasticity for boundary displacements and pore pressures is extended to permit the computation of stresses at both boundary and interior points. The stress formulation preserves the surface-only discretization. The boundary displacement integral equation is progressively differentiated to obtain the related stress and strain integral equations. Explicit expressions for the steady-state axisymmetric fundamental solutions are derived in this process. The transient components of the integrands are obtained directly from the transformation of the three-dimensional kernels into a cylindrical system. Numerical implementation of these integral equations is carried out within a general purpose BEM computer code and several illustrative examples are presented to validate the method. 相似文献
15.
Boundary setting method for the seismic dynamic response analysis of engineering rock mass structures using the discontinuous deformation analysis method 
下载免费PDF全文
下载免费PDF全文 Large deformations and discontinuous problems can be calculated using the discontinuous deformation analysis (DDA) method by solving time steps, and this method is suitable for simulating the seismic dynamic response of engineering rock mass structures. However, the boundary setting must be carefully analyzed. In this paper, four boundary settings for the DDA method are investigated. First, the contributions to the DDA equations for nonreflecting boundaries (including the viscous boundary and the viscoelastic boundary) are deduced based on the Newmark method. Second, a free‐field boundary is introduced in the DDA method with boundary grid generation and coupling calculation algorithms to accurately simulate external source wave motion, such as earthquakes. Third, seismic input boundary treatments are intensively examined, and the force input method is introduced based on nonreflecting boundaries. Finally, the static‐dynamic unified boundary is implemented to ensure consistent boundary transformation. The boundary setting method in the DDA method is discussed, and the suggested treatments are used to analyze the seismic dynamic response of underground caverns. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
16.
17.
An analytical solution is presented in this paper to study the time‐dependent settlement behaviour of a rigid foundation resting on a transversely isotropic saturated soil layer. The governing equations for a transversely isotropic saturated soil, within Biot's poroelasticity framework, are solved by means of Laplace and Hankel transforms. The problem is subsequently formulated in the Laplace transform domain in terms of a set of dual integral equations that are further reduced to a Fredholm integral equation of the second kind and solved numerically. The developed analytical solution is validated via comparison with the existing analytical solution for an isotropic saturated soil case, and adopted as a benchmark to examine the sensitivities of the mesh refinement and the locations of truncation boundaries in the finite element simulations using ABAQUS. Particular attention is paid to the influences of the degree of soil anisotropy, boundary drainage condition, and the soil layer thickness on the consolidation settlement and contact stress of the rigid foundation. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
18.
19.
介绍了联合运用积分方程法与有限元法(简称混合法)来计算三维非均匀介质中电磁场分布的理论方法,并进行了数值模拟实验。混合法的原理是引入一个包围非均匀目标体的虚构边界,在边界内部的场用有限元法模拟,在边界外部(包括边界)的场用积分方程表达,二者在边界上通过场的连续性耦合起来。数值实验结果表明,混合法既能显著地减小网格规模,又能灵活地模拟复杂的介质情况,且计算精度较高。 相似文献
20.
三维地形大地电磁场的边界元模拟方法 总被引：1，自引：0，他引：1
提出了一种用边界元法计算大地电磁场三维地形影响的数值模拟方法.首先用矢量积分理论和电磁场边界条件, 将上半空间(空气)和下半空间(地下介质)两个区域电磁场边值问题变为仅对地形界面的两个矢量面积分方程, 其中一个计算磁场, 称磁场方程; 另一个计算电场, 称电场方程.然后将对地形界面的积分剖分为一系列的三角单元积分.在三角单元积分中, 假设单元中电磁场为水平均匀大地空间电磁场与地形影响的迭加, 并假设地形影响为常项, 这样既保证了计算精度又使得计算方法简便.通过分解和计算, 每一个矢量面积分方程分解为对应3个坐标方向的3个常量线性方程, 这些线性方程组成了对角占优的线性方程组, 可用SSOR方法求解.文中给出了2个三维地形上大地电磁视电阻率曲线的计算结果. 相似文献