共查询到20条相似文献，搜索用时 219 毫秒
1.
边界元法在数值解法中利用拉普拉斯方程的基本解和格林公式,简单易行,可以解决物性参数分区均匀的一类正演问题(电法和重、磁法勘探问题),也可计算向上延拓.特别是将它与有限元法结合,解决复杂边值问题,可取得较理想的效果. 边界元法的基本计算方法有直接法和间接法两种. 边界元法的直接法现以磁法勘探正问题为例来说明. 解磁法正问题可归结为在求解区域Ω上解拉普拉斯方程的边值问题,并可用边界元法求解, 相似文献
2.
边界元法是一种新的微分方程数值解法.本文讨论了用该法求解点电源场中任意起伏地形下二维地电体位场问题的方法原理和数值处理技术,并给出了不同的地电体视电阻率与视极化率的计算结果相应用实例. 结果表明:边界元法是求解电法勘探边值问题和资料处理的一种经济有效的数值计算方法. 相似文献
3.
本文用边界元法解决非稳定系统中水质预测的计算问题,建立了三维水质预测问题的边界元法计算模型,为简化计算手续、提高计算精度,讨论了Laplace变换边界元法,推广了Laplace变换边界元法中的Schapery方法。 相似文献
4.
扶永铭 《成都理工学院学报》1995,22(3):109-114
该文研究了低频人工时变场（似稳场）中任意导电导磁矿体的边界元素法。通过对矢量位的边界积分方程，再以样条边界元素法进行离散，得到实分量的常系数方程组和虚分量的二阶微分方程组。此方程可利用时间积分的样条函数方法求解，解出适量位后，利用磁场强度与矢量位Ａ的相应关系而求得二次磁场的水平分量和垂直分量。 相似文献
5.
本文讨论了在均匀谐变场激发下导电导磁体电磁响应的边界单元法计算问题。文中首先建立矢量位的边值问题,然后导出其边界积分方程,最后用边界元法进行求解。本方法的计算结果具有精度高、速度快的优点。 相似文献
6.
结合边界元法和离散元法,提出一种可以进行计算颗粒内部应力和破碎路径的方法。该方法利用离散元法求解颗粒的相互作用和每个颗粒上的荷载。然后利用边界元法计算颗粒的应力分布,为了实现动态平衡,将颗粒的加速度视为恒定大小的体力。但体力导致边界积分方程中出现域积分,故采用直线积分法将域积分转化为边界积分,以保证边界元法降维的优势。为了提高边界元的计算效率,对于几何形状相似的颗粒,以其中一个颗粒作为模板颗粒,只需要计算模板颗粒在局部坐标系中的系数矩阵,其他相似颗粒可以通过局部和全局坐标系之间的映射获得。在得到应力后,基于Hoek-Brown准则来判断颗粒是否破碎。此外,将破坏路径简化为直线,并采用最小二乘法拟合得到破坏路径。 相似文献
7.
离散元法及其在岩土工程中的应用综述 总被引：16，自引：0，他引：16
离散元法是基于不连续性假设的数值方法,它特别适合于求解节理岩体中的非连续性问题。在介绍离散元法基本原理的基础上,着重对离散元法在岩土工程领域的应用现状作了叙述和分析,并对其发展趋势进行了探讨。 相似文献
8.
9.
10.
考虑表皮效应竖向非均匀分布,建立各向异性承压层3D稳态井流数学模型,为求得问题的解,首先将表皮效应系数视为分段连续函数,再根据该分段函数将承压层划分成N层,然后通过差分法将问题转化为求解矩阵微分方程边界值问题,最后通过矩阵理论求得相应问题的解。应用数学软件对所求得解进行编程,并应用于实例计算和分析中。结果表明,在距竖井较近处,承压层水头降和流量沿竖向分布均与表皮效应系数函数变化趋势相反,即承压层某深度处的表皮效应系数越大（小）,则相应位置的水头降和流量越小（大）;而距竖井较远处,仅当竖向渗透系数相对于水平渗透系数较小时,才仍具有前述变化规律;承压层不同深度处的水头降沿径向的分布曲线的幅度值,不仅取决于该处的表皮效应系数大小,而且取决于竖向渗透系数与水平渗透系数相对值大小 相似文献
11.
基于边界元法的边坡矢量和稳定分析 总被引：4，自引：0，他引：4
矢量和法物理力学意义明确,计算简单,且能根据边坡当前的应力分布状态合理地评价其整体稳定性状态。其中边坡的应力状态通常是采用有限元法来求解。由于边界元法具有研究问题降阶、离散化带来的误差值仅产生在边界以及计算量小等优点,在工程中得到了广泛应用;对于平面问题,以源点作为原点,以所积分单元的切向和法向为坐标轴建立局部坐标系,对于线性单元可以得到所有积分的解析解。因此,可以得到计算区域内部任意点的场变量的解析解,这就保证了位于边界附近区域场变量的精度。利用边界元法得到二维边坡体内连续的应力分布状态,使用矢量和法对该边坡进行稳定性分析,并且与基于有限元的矢量和法、极限平衡法进行对比分析。边坡圆弧滑面和折线滑面的计算结果表明,基于边界元法得到的矢量和安全系数和基于有限元的矢量和法、极限平衡法基本一致;边界元法对应的矢量和安全系数对边界单元尺寸不敏感。 相似文献
12.
最小二乘配点法是用于地下水流计算的一种新型、高效的无网格方法。此方法是在径向基函数配点法的基础上,对计算域进行节点离散,并布置辅助点,近似函数仍然通过节点构造,微分方程在节点和辅助点上都严格满足,从而计算精度更高。而且此方法不需要背景网格,效率高,形式简单。利用该方法计算地下水流向河、渠中的非承压含水层稳定流和非稳定流问题,算例表明,该方法有很好的精度且计算量小,比径向基函数配点法有更精确的结果。 相似文献
13.
本文阐述了布西涅斯克方程的近似解析解,文中假定具有不变的初始条件,并且边界上的测压水头成阶梯函数增加。利用波尔兹曼变换将这祥的一维问题简化为常微分方程问题,并利用精确解的一些基本特性的方法导出这种问题多项式的近似解。这种方法也适用于一维非线性弥散问题。文中表示了与精确解的一些解析性质有关的数值过程。论文中进行了讨论和推断。 相似文献
14.
利用间接边界元法,在频域内求解了层状弹性半空间中洞室对入射平面SV波的散射问题。通过自由场反应分析,求得假想洞室边界上各点位移和各单元应力响应。在洞室边界各个单元上施加虚拟分布荷载,求得位移和应力的格林函数。根据应力边界条件确定虚拟分布荷载,将自由场位移响应和虚拟分布荷载产生的位移响应叠加起来,即得到问题的解答。比较了层状半空间和均匀半空间中洞室对入射平面SV波的放大作用。结果表明,层状半空间情况有可能导致较大的地表位移幅值,尤其是对于较低频率入射波。 相似文献
15.
边坡稳定性分析的极限平衡法是要做一定假设的。本文考虑了离散元与极限平衡法基于刚体理论的相似性,用离散元法计算出极限平衡法中的各个计算分量值,然后代入到常用极限平衡法中分析计算,根据离散元法的计算结果讨论了各极限平衡法中各种假设的合理性,为极限平衡法的计算提供了一个新思路。结果表明Bishop法与离散元法的计算结果比较接近,是本文经过离散元法计算后所推荐使用的方法。 相似文献
16.
弹丸侵彻无钢筋混凝土的数值模拟 总被引：1，自引：0，他引：1
采用二维梁-颗粒模型BPM2D（beam-particle model in two dimensions）模拟了刚性弹丸侵彻无钢筋混凝土的过程。离散元法（DEM）和有限元法（FEM）等数值计算方法各有其优势,同时也都存在不足之处：离散元法适于处理由连续介质向非连续介质转化的破坏问题,但对于连续体计算结果精度不高;有限元法适于预测材料破坏的区域,但难以直接用于计算脆性材料破坏过程,因此将两种方法结合可以形成一种较好的混合模型。梁-颗粒模型BPM2D是基于离散元法,结合有限元法开发的二维数值计算模型,采用3种类型梁单元形成混凝土数值试样,每种类型梁单元的力学性质均按韦伯（Weibull）分布随机赋值,以模拟混凝土细观结构的非均匀性,同时梁单元的强度随应变率不同而变化。利用此模型分析了弹丸侵彻下混凝土的破坏过程,并给出侵彻过程弹丸减速度-时间历程曲线。比较计算结果与试验数据表明梁-颗粒模型可有效应用于计算和模拟脆性材料动态破坏问题。 相似文献
17.
将复合地基视为横观各向同性体,利用半解析数值计算原理,又将复合地基位移竖向离散为多项式,在水平两个方向上选取解析函数模拟无限域地基进行计算,从而克服了有限元计算中人为设定计算边界的弊端,提高了计算精度,并将三维沉降问题转化为一维问题进行计算,收敛速度快,显著地减少了计算工作量,将计算结果与理论解进行对比,满足精度要求。通过改变复合地基置换率、加固深度和桩土模量比,采用半解析元法研究了复合地基的加固效果,计算结果表明,桩土模量比较大时,可忽略加固区压缩变形量;增大置换率对减小沉降作用很小,而增加加固深度能明显减小复合地基沉降,在以控制沉降量为主要目的复合地基优化设计中应首先考虑增加加固深度。 相似文献
18.
基于MatDEM的砂土侧限压缩试验离散元模拟研究 总被引：2，自引：2，他引：0
离散元法基于非连续介质力学理论,尤其适用于砂土等离散介质体的数值模拟研究。利用岩土体离散元模拟软件MatDEM的二次开发功能,研发了砂土侧限压缩试验三维离散元模拟器。对三个不同级配砂土试样进行了侧限压缩试验,并且进行了与之相对应的数值模拟,通过分析对比试验结果与数值模拟结果验证了所开发模拟器的有效性。模拟结果表明:离散元法可以很好地反应砂土压缩过程中的配位数变化;每个模拟样品中,粒径较小的单元受到较大的平均压力,导致平均位移较大;数值计算结果的主要误差是由离散元颗粒自身的泊松比引起的。研究突破了常规土力学研究方法的局限性,为今后岩土工程离散元模拟研究提供了参考。 相似文献
19.
本文给出了瞬变弹性动力问题的一个边界元法。该法是利用威尔逊——θ法的差分公式,把运动方程化为椭圆型微分方程。根据贝蒂定理和动力点荷载的特解,可获得动力问题的边界积分方程。这个解法是在真实时间域内逐步求解的,不需要使用拉氏变换。一个应力波传播的数值算例证实了该方法使用方便,且解答精度较高。 相似文献