共查询到20条相似文献，搜索用时 107 毫秒
1.
本文提出的有限单元和边界单元交连耦合法(CFEBEM)可作为地下水流场和水文地质参数正反演计算的有效工具。该方法改进了Liggett(1983)的代入耦合法,通过内边界节点组和交连方程,实现有限单元和边界单元的多区耦合。CFEBEM方法发挥了有限单元和边界单元各自的特长,具有单元区划分灵活、方程组装简便和计算精度高等特点,可有效地解决结构复杂含水层及大型密集抽水条件下的流场计算和水文地质参数反演。在太原岩溶地下水资源评价和合理开发研究中的应用,证明了这种方法的有效性。 相似文献
2.
简述了复合单元法的基本原理和算法,该方法的主要优点在于可以明显地减少前处理的工作量和难度,且可以与有限单元很好地耦合在一个系统中。介绍了宝珠寺重力坝的物理模型超载试验。分别采用三维弹黏塑性复合单元法和三维弹黏塑性有限单元法对宝珠寺重力坝的超载工况进行了数值试验,并与物理试验进行比较。结果表明,三种方法得到的特征点位移吻合得较好,两种数值试验方法的整体位移和应力也基本一致。研究验证了复合单元法的正确性,也说明了该算法的主要优点：前处理简单,计算精度与有限单元法相当,同时还证实了复合单元单元法可以较好地解决复杂岩基上重力坝的变形与破坏问题。 相似文献
3.
三维复杂域非结构化六面体网格全自动生成方法 总被引：3，自引：0，他引：3
使用自适应有限单元法在三维复杂域内求解水工及岩土工程问题,讨论了应用行波法在多介质三维复杂块体内六面体单元网格的自动生成过程。其特点是对不规则区域具有较强的适应性,并考虑了节理单元的自动生成。 相似文献
4.
《岩土力学》2017,(5):1463-1471
试验研究和数值模拟都表明,颗粒破碎对颗粒岩土材料的宏观力学响应具有重要的影响。结合离散元与比例边界有限元法,提出一种用于研究堆石料颗粒破碎行为的新的数值计算方法,可以将两种方法各自的优势同时发挥出来。离散单元法将被用于求解颗粒的运动以及相互作用。在每一个时间步尾,采用比例边界有限单元法计算颗粒内部的应力分布。比例边界有限单元法可以仅用一个任意边数的多边形来描述一个颗粒,大大节约应力计算的成本。当应力状态确定之后,Hoek-Brown准则将被引入用于确定颗粒内部的破坏点,当破坏点达到一定的比例时,即视为这个颗粒将破坏。为简化起见,假设破坏路径为直线,如果破坏发生,则该颗粒一分为二,各个新的颗粒将直接参与离散元和比例边界有限元的计算。该方法无需预定义任何子颗粒或网格重划分。最后以一个简单的双轴试验的模拟展现方法的可行性。 相似文献
5.
岩土工程数值计算中的无网格方法及其全自动布点技术 总被引：14，自引：3，他引：11
自然单元法采用无网格的思想全域构造插值函数,它的求解精度高,计算时间少,可准确地施加边界条件,兼具有无网格法和有限单元法的优点和特点,是一种理想的用于岩土及地下工程分析计算的数值方法。文中简要地介绍了自然单元法的基本理论,并针对岩土及地下工程问题特点,给出了一种无网格离散点的全自动布置方法。 相似文献
6.
Du Fengwen 《吉林大学学报(地球科学版)》1989,(2)
本文根据灰色系统理论,通过建立渗流域边界上水位灰色预测模型来模拟有限单元法预报水位的边值条件,从而提高有限单元法预报水位的可靠性。 相似文献
7.
通过对复合地基加固路堤稳定性的极限平衡法和有限单元法的计算结果进行对比分析,揭示出当路堤的稳定性由复合地基决定时,极限平衡法和有限单元法的计算结果存在较大差异。应用岩土有限元软件Plaxis中的强度折减法可得到合理的路堤稳定性分析结果。对复合地基中桩体破坏模式的分析认为,桩体发生非剪断破坏之外的弯曲、转动、拉伸等破坏模式是极限平衡法与有限单元法计算结果产生差异的根本原因。由此进一步指出,当处理存在土与结构物相互作用的边坡稳定问题时,极限平衡法的分析结果可能会高估了边坡的稳定性,应慎重判别其合理性。 相似文献
8.
9.
基于三维颗粒离散单元法,赋予颗粒相应的细观参数,并采用黏结发生在接触颗粒间有限范围内的模型来考虑冻土颗粒中冰的胶结作用,建立了冻结黏土三维离散元数值模型.在相同围压、不同温度和相同温度、不同围压下对冻结黏土的室内三轴试验进行数值模拟,对比了数值试验与室内测试的应力-应变曲线,两者吻合较好.数值模拟结果表明:围压增大会使得接触黏结逐渐失效,在剪切带中胶结冰的破坏区域将增大,而温度的降低则会产生相反结果,这些微观变化都将对冻结黏土的宏观力学变形产生较大影响,同时,细观参数对温度的依赖性也很明显.冻结黏土三轴试验微观变形离散元模拟思路及方法可为今后运用离散单元法研究冻土力学行为提供一定的参考. 相似文献
10.
滑坡稳定性分析与安全系数取值研究 总被引：5，自引：0，他引：5
在地质灾害中,滑坡灾害分布最广、发生频率最高、危害最大,是我国地质灾害防治的主要对象。滑坡研究的主要任务之一是稳定性分析。论文针对三板溪水电站进水口滑坡和水布垭水电站大岩淌滑坡,分别采用刚体极限平衡法、块体单元法、有限单元法对稳定安全系数和条间推力分布曲线进行分析和比较,研究各种方法的差异和计算精度,以及各种方法稳定安全系数的取值标准,希望为水利水电工程边坡设计规范有关有限单元法、块体单元法条目的编写和相应的允许安全系数取值标准的制定提供一定的参考。研究结果表明:①块体单元法计算的安全系数与刚体极限平衡法计算的安全系数可以采用相同的取值标准,而有限单元法计算的安全系数取值标准可适当降低;②块体单元法计算的条间水平推力、条间剪切力与刚体极限平衡法(尤其是Sarma法)计算的条间水平推力、条间剪切力非常相近,有限单元法计算的条间水平推力、条间剪切力较小;③分析方法、力学模型与参数取值应与安全系数取值标准相匹配。 相似文献
11.
强度折减有限元法中的单元阶次影响分析 总被引：1，自引：0，他引：1
强度折减有限元法是当前较为有效的边坡稳定性评价方法,且应用越来越广泛。但影响强度折减有限元法的因素有很多,单元阶次是其中比较重要的一个。通过3个经典算例,这些算例分别是二维地基承载力问题、二维边坡和三维边坡问题,分析了单元阶次的选择对强度折减法的影响。计算结果表明：随着单元的增多,线性单元和二次单元都从大于真实解的一侧来逼近真解;相对于二次单元,由于线性单元过“刚”,因此,会过高地估计安全系数,对于实际工程会偏于危险,且误差大,二次单元的误差是线性单元误差的1/8左右。在采用系统最大位移收敛与否的评判标准的基础上,利用二次单元来进行强度折减分析,则可以弥补这种线性单元的不足,得到更加合理的安全系数。二次单元比线性单元更适合于强度折减有限元法。 相似文献
12.
Robert Kirby 《Computational Geosciences》2003,7(3):197-214
In this paper, residual-based a posteriori error bounds are derived for the mixed finite element method applied to a model second order elliptic problem. A global upper bound for the error in the scalar variable is established, as well as a local lower bound. In addition, due to the fact that the scalar and vector variables are approximated to equal order accuracy, the dual problem may be modified to give an upper bound for the vector variable. Some comments on estimating more general error quantities are also made. The estimate effectively guides adaptive refinement for a smooth problem with a boundary layer, as well as detects the need to refine near a singularity. 相似文献
13.
用区间有限元计算边坡荷载组合效应 总被引：2，自引：2，他引：0
边坡工程的稳定性受多种因素的影响,其荷载具有不确定性,可以用具有上、下限的区间数来表示,提出了用区间数学和有限元结合构成的区间有限元方法来求解边坡变形。以古树包边坡的变形为例,研究了该方法的可行性,并与常规有限元计算结果进行了比较,实例计算结果表明,该方法用于计算具有不确定性参数和荷载的工程问题的变形是有效的,与常规有限元的要考虑多种工况相比,计算工作量大大减少,使用简便,而计算结果却能满足工程需要,适合计算边坡工程。 相似文献
14.
针对分量式钻孔仪记录同震应变阶跃时常出现的ΔεⅠ+ΔεⅢ与ΔεⅡ+ΔεⅣ观测值不等, 无法自检, 且观测到的应变/应力阶跃大于理论计算值等现象, 利用有限元方法计算分析可能的原因, 建立不同围压和地震波应力扰动条件下的多种模型, 运用各向同性“杀死单元”和横向各向同性"杀伤单元"法, 模拟分量式钻孔仪的金属外壁和固结水泥环之间产生张裂隙或剪切滑动的情况。计算结果显示, 若钻孔内的金属圆筒和水泥环之间发生剪切滑动, 基本不会影响到观测; 然而两者之间发生张裂会使得ΔεⅠ+ΔεⅢ与ΔεⅡ+ΔεⅣ出现较大差异, 严重影响到同震阶变的观测。今后在井下安置探头时应该尝试新的改进方法。 相似文献
15.
David Trujillo 《Computational Geosciences》2004,8(2):173-185
A mixed finite volume method is applied in order to simulate a simplified Far Field model. We show how this method can be adapted to both the elliptic and the convection–diffusion equations describing the water flow and the variations of the concentration of the radioactive element respectively. For the water flow problem, we justify an a posteriori error estimator. Finally, we present some numerical results. 相似文献
16.
研究了三维渗流的自适应有限元法,建立了集成网格自动剖分、渗流有限元分析、网格调整、后处理等四大模块的四面体单元h型自适应系统。给出的算例证实了系统的有效性和优越性。 相似文献
17.
基于边界元法的边坡矢量和稳定分析 总被引：4，自引：0，他引：4
矢量和法物理力学意义明确,计算简单,且能根据边坡当前的应力分布状态合理地评价其整体稳定性状态。其中边坡的应力状态通常是采用有限元法来求解。由于边界元法具有研究问题降阶、离散化带来的误差值仅产生在边界以及计算量小等优点,在工程中得到了广泛应用;对于平面问题,以源点作为原点,以所积分单元的切向和法向为坐标轴建立局部坐标系,对于线性单元可以得到所有积分的解析解。因此,可以得到计算区域内部任意点的场变量的解析解,这就保证了位于边界附近区域场变量的精度。利用边界元法得到二维边坡体内连续的应力分布状态,使用矢量和法对该边坡进行稳定性分析,并且与基于有限元的矢量和法、极限平衡法进行对比分析。边坡圆弧滑面和折线滑面的计算结果表明,基于边界元法得到的矢量和安全系数和基于有限元的矢量和法、极限平衡法基本一致;边界元法对应的矢量和安全系数对边界单元尺寸不敏感。 相似文献
18.
弹簧元法是一种将单元离散为一系列弹簧的数值计算方法。不同的单元具有不同的离散方式,确定相应的离散弹簧的刚度系数表达式是弹簧元法的关键。将四节点矩形单元离散为6个基本弹簧,每个基本弹簧包括法向弹簧和切向弹簧两个派生弹簧,并用泊松弹簧和纯剪弹簧描述单元的泊松效应和剪切效应,用有限元的单元刚度矩阵标定各弹簧的刚度系数,实现了一种四节点矩形弹簧元的构造形式。该单元的同类弹簧具有相同的表达形式。法向与切向弹簧的刚度表达式中分别含有法向和切向弹簧刚度待定系数。通过改变待定系数的值可使该单元分别对应于有限元的常应变、双线性及Wilson非协调单元。将上述弹簧元方法进行理论推导,并应用于基于连续介质的离散单元法（CDEM）的核心计算进行简单算例验证,证明了提出方法的正确性。通过以上研究发现,四节点矩形弹簧元有以下特点：对于相同问题,不同单元有不同的计算精度;对于梁弯曲问题,应用该单元可显著提高离散单元法的求解精度;改变待定系数的值,可得到更高或者更低精度的单元。 相似文献
19.
岩土介质三维快速拉格朗日数值分析方法研究 总被引：12，自引：0，他引：12
研究了FLAC3D（fast lagrangian analysis of continua in 3 dimensions）的特点,并与有限单元法作了比较。FLAC3D方法以结点运动方程为支配方程,追踪了介质从受荷到达到平衡状态的过程,而有限元法是根据介质力学平衡方程直接求解,这是二者主要区别。FLAC3D没有采用介质真实的阻尼特性和结点质量,给出的不是介质所经历的真实过程,不能正确反映过程的影响,因此给出的介质应力和变形计算结果的物理意义是不甚明确的。求解过程中的介质振动,是一种噪音,可引起弹塑性介质计算结果误差,而弹性介质的计算结果几乎不受影响。研究指出了FLAC3D方法的优缺点。还通过算例作了四个方面的研究：(1) Drucker-Prager屈服准则与Mohr-Coulomb屈服准则比较;(2) 膨胀角取值对计算结果的影响;(3) 大变形与小变形对计算结果的影响;(4) 精度设置对计算结果的影响。研究表明,Drucker-Prager准则与Mohr-Coulomb准则结果差异颇大;膨胀角取值对结果的影响是敏感和显著的;一般情况下,取小变形模式是合适的,计算精度取10-5是足够的。 相似文献
20.
针对移动荷载引起地基振动问题的特点,在随荷载移动坐标系下建立有限元方程进行分析,能够与常规有限元方法较好地符合。建立了可用于移动坐标有限元法的黏弹性边界,从而有效地模拟波在边界的透射。以Euler-Bernoulli梁模拟轨道和路基的作用,以三维实体单元模拟黏弹性分层地基,得到的数值模拟结果能够与瑞典Ledsgård地区高速铁路实测数据较好地吻合,表明移动坐标有限元法可以有效地分析此类问题。 相似文献