非饱和带有限分析数值模拟的误差分析

由于描述包气带水分运移的Richard方程具有高度非线性,常用数值方法求解。有限分析法可以较好地保持原有问题的物理特性,近年来被广泛应用于环境及农业工程领域。基于不同的原理,目前有两种有限分析法模拟非饱和带水分运移问题,即基于局部线性化有限分析法和基于Kirchhoff变换有限分析法。为探讨两者在求解非饱和带水分运移问题时数值表现的差异性,采用已有解析解对两种方法进行初步验证;使用两个数值实验比较两种方法的计算精度;应用实测数据验证Kirchhoff变换有限分析法。相比于局部线性化有限分析法,Kirchhoff变换有限分析法能够更好地控制误差及获得更高精度的数值解。该研究对于完善非饱和带水分运移有限分析数值模拟的理论与方法有重要意义。     

复杂岩体地下水运动问题的有限分析法

在数值模拟方法基础上,提出了一种全局坐标下解复杂岩体地下水运动问题的有限分析法。算例表明,这种方法不仅可行,而且具有较高的数值求解精度。     

求解地下水非稳定流问题的局部坐标有限分析法

给出了一种求解地下水非稳定流方程的局部坐标有限分析法，证明了格式的稳定性。用该格式计算了两个典型问题，表明该方法有较高的计算精度，并能运用几何形状较复杂的渗流问题的计算。     

求解地下水非稳定流问题的局部坐标有限分析法

给出了一种求解地下水非稳定流方程的局部坐标有限分析法，证明了格式的稳定性。用该格式计算了两个典型问题，表明该方法有较高的计算精度，并能适用几何形状较复杂的渗流问题的计算。     

边坡大变形弹塑性有限元分析[Ⅰ]: 理论

本文采用以物理坐标为变量的Ｕｐｄａｔｅｄ－Ｌａｇｒａｎｇｉａｎ描述方法,在考虑边坡非线性变形破坏的基础上,建立了边坡大变形破坏的弹塑性有限元模型,详细推导了有限单元方程和数值分析公式,并考虑了边坡的初始应力和超孔隙水压力。     

有限分析法在浅水湖泊水流水质计算中的应用

根据有限分析法能准确地模拟对流效应和有效消除伪数值扩散等特点，求解了流函涡量方程，以滇池为例，模拟计算了湖泊在有风条件下的流场，浓度场，计算结果与实测情况吻合较好。     

软岩巷道围岩流变大变形有限元计算方法

本文针对软岩工程中围岩流变大变形问题，以有限变形力学理论为依据，采用更新拖带坐标法编制非线性大变形有限元计算程序，对软岩巷道流变大变形问题进行了数值模拟计算。在计算中，选用了适合描述软岩流变大变形的物性方程，并考虑了锚杆和拱两种支护方式。结合表明，本文采用的理论和计算方法合理，计算程序具有精度高、收敛速度快等优点，展示了良好的应用前景。     

移动荷载下瞬态振动的三维移动坐标有限元分析

针对移动荷载引起地基振动问题的特点,在随荷载移动坐标系下建立有限元方程进行分析,能够与常规有限元方法较好地符合。建立了可用于移动坐标有限元法的黏弹性边界,从而有效地模拟波在边界的透射。以Euler-Bernoulli梁模拟轨道和路基的作用,以三维实体单元模拟黏弹性分层地基,得到的数值模拟结果能够与瑞典Ledsgård地区高速铁路实测数据较好地吻合,表明移动坐标有限元法可以有效地分析此类问题。     

三维磁场有限元—无限元耦合数值模拟

利用传统有限单元法在有限空间范围内开展三维地球物理场正演模拟时,由于截断边界的影响,会引起局部异常的畸变,影响数值模拟的精度,对该问题通常采用扩边的办法加以解决,但需要的范围较大,从而大大增加运算成本,影响正演模拟效率。本文基于COMSOL Multiphysics软件,在求解域外部边界设置无限元以替代传统边界条件,达到减小计算区域目的。通过孤立球体和组合模型磁场正演模拟,考虑退磁、剩磁及地表起伏条件,与传统有限元方法相比,有限元—无限元耦合算法能够有效克服边界效应,提高计算精度,降低运算量,从而提高了有限单元法正演数值模拟效率。     

饱和砂土局部变形带模拟的有限元数值实现

基于有限变形理论,推导了Newton-Raphson 迭代算法在k+1步增量表达的矩阵形式,实现了饱和砂土变形局部化的有限元数值计算,得到了饱和砂土发生局部化变形的准则。基于Galerkin 方法,得到了位移场和应力场的空间离散化矩阵方程;由土体局部变形带的连续性条件,引入第1切线算子,推导出了砂土等颗粒状媒介发生局部化变形的必要条件。基于此核心算法,编制了有限元计算程序,模拟了饱和砂土在不排水条件下平面压缩过程中剪切带的形成与发展;通过比较分析,研究了有限元网格粗细对于土体局部变形带的影响,结果表明,网格粗细的病态依赖只是微小的,它只与变形条带的宽度有关,对于土体所表现出来的其他力学特性没有影响。     

土工格栅加筋砂土的变形与破坏机理解析

利用可考虑局部破坏的非线性弹塑性有限元，对无加筋和加筋砂土的平面应变压缩试验结果进行了从小变形到破坏的全过程数值解析。加筋砂土试验体用土工格栅分6层和11层进行加筋加固。将等价二维有限元解析所得到的解析结果与试验得到的实测值进行了较为伞面的比较，结果表明：合理的二维非线性弹塑性有限元解析，不仅可以较为精确的模拟加筋砂土的平均心力-应变特性，而且还可以全面地调查试验体的局部应力-应变分布以及剪切破坏发生状况，从而，对加筋砂土的变彤破坏以及加筋材的加州机理有一个更加全面合理的认识。     

复合地基承载及变形特性的有限元—有限层法分析

本文采用有限元－有限层数值法分析了复合地基了载及变形特性的影响因素，并用工程实例并进行了计算对比。     

一种简单的网格图形局部放大技术

本文利用视口技术，根据坐标平移的原理，设计出一个用于有限元网格图局部放大的程序，并在微机上调试运行通过，本方法亦可应用于其他类型的曲线图放大。     

桩基水平位移的有限元模拟

有限元数值计算方法在岩土工程中已有较多的应用,但在水平受荷桩分析计算中的应用尚不普遍。有关水平受荷桩土体系水平承载力及变形的数值模拟方法少有文献报道。采用有限元数值模拟方法来研究水平受荷桩土体系的水平位移。首先简要分析了水平受荷桩的主要计算方法及原理。然后,基于ABAQUS软件建立了桩土系统二维及三维的有限单元模型,分别通过按现行规范推荐的m法和按Poulos的弹性理论法计算的两个算例,将数值计算方法与解析方法进行了对比计算分析。有限元数值计算结果与解析解基本吻合,说明了所采用有限元数值模拟方法的可靠性及有效性。     

贮量集中有限元方法在渗流计算中的优点

普通有限元方法和有限差分方法是渗流计算中常用的数值方法,而贮量集中有限元方法则很少被实际应用。为了使渗流计算者全面了解贮量集中有限元方法,并能够较多地应用这种方法,本文将渗流计算中的贮量集中有限元方法、普通有限元方法和有限差分方法进行了对比,具体比较和分析了有限元方程组的系数矩阵特征、计算工作量,通过具体的算例分析了它们的计算精度,由此指出贮量集中有限元方法的计算精度与有限差分法基本相同,而远好于普通有限元法,它的计算工作量小于其他二种方法。     

电阻率各向异性介质大地电磁二维非结构有限元数值模拟

介质的各向异性会引起大地电磁观测曲线的畸变,对大地电磁数据解释影响很大。这里从麦克斯韦方程组出发,推导了基于二阶插值基函数的二维大地电磁非结构有限元公式,实现了任意地形起伏条件下任意各向异性二维大地电磁响应的数值模拟。各向异性断层模型的数值模拟结果与解析解吻合一致,复杂各向异性模型的有限元模拟与已有的有限差分计算结果也相符合,验证了算法的正确性。在此基础上探讨了各向异性参数对二维大地电磁响应的影响。     

地下人防通道对周边建筑影响分析方法浅析

为了分析地下人防通道对周边建筑的影响,结合该场地工程地质情况,采用Peck法、Plaxis有限元分析法、ABAQUS有限元分析法三种理论分析方法,综合确定了人防通道沉降影响范围,表明数值计算与现场实际调查结果相符合。研究对分析类似场地条件下地下工程对周边建筑影响具有一定的借鉴意义。     

承压含水层中地下水向井非稳定流动的LT有限分析法

作者给出了一种求解承压非稳定井流方程的ＬＴ有限分析法．其特点是：将时间ｔ直接包含在有限分析数值计算格式中，克服了传统数值法求解非稳定过程按时段迭代的缺陷，实现了非稳定过程计算的一步完成，给多时段地下水资源管理和长期预报带来了方便．用本文提出的方法分别计算了两个算例，结果表明，该方法具有较高的精度，与解析解吻合甚好，可用于实际地下水流的计算之中．     

贮量集中有限元方法在渗流计算中的优点

普通有限元方法和有限差分方法是渗流计算中常用的数值方法,而贮量集中有限元方法则很少被实际应用。为了使渗流计算者全面了解贮量集中有限元方法,并能够较多地应用这种方法,本文将渗流计算中的贮量集中有限元方法、普通有限元方法和有限差分方法进行了对比,具体比较和分析了有限元方程组的系数矩阵特征、计算工作量,通过具体的算例分析了它们的计算精度,由此指出:贮量集中有限元方法的计算精度与有限差分法基本相同,而远好于普通有限元法,它的计算工作量小于其他二种方法      

可变网格的有限单元分析在边坡稳定计算中的应用

本文在非线性有限单元法的基础上，引进八节点的狭长单元来描述边坡中的滑动面，同时对于有限元单元网格可调节点坐标进行变更形成新的网格，于是利用原有基本网格架，通过多次有限元分析求出边坡中可能滑动面的稳定系数曲线Ｋ＝ｆ（），并进而确定最危险的滑动面，在此基础上开展一系列稳定分析工作。