共查询到10条相似文献，搜索用时 31 毫秒
1.
数值流形方法中一般采用有限元网格或规则网格作为其数学覆盖系统,而规则的网格突出的优点是不需要适应求解域的边界和各种不连续面。采用规则的矩形网格作为数值流形方法中的数学网格,并借助适合分析的T样条实现了数值流形方法中的局部加密。适合分析的T样条定义在一个限制的T网格上,其基函数具有线性无关、单位分解、局部加密等许多重要性质,使得其非常适合用于工程设计及分析。当对一个适合分析的T网格加密后,所产生的新的网格往往不再是适合分析的T网格。基于此,提出了一种简单的数学网格加密算法,该算法能保证局部加密后的数学网格仍然是适合分析的。算例结果表明:在应力集中区和裂纹尖端等应力梯度较大区域,该算法均具有较强的适用性。 相似文献
2.
数值流形方法（NMM）基于两套覆盖（数学和物理覆盖）和接触环路而建立,能够统一地处理岩土工程中的连续和非连续变形分析问题。与其他基于单位分解理论的数值方法一样,NMM可以自由地提高物理片上局部位移函数（多项式）的阶次,从而在不加密网格的情况下显著地提高计算精度,但有可能会使总体刚度矩阵奇异,产生线性相关问题。针对这种情况,引入了一种新的高次多项式形式的局部位移函数,在此基础上,建立了新的NMM求解体系,并应用于求解一般的弹性力学问题。结果表明：它有效地消除了线性相关问题;较之传统局部位移函数取一次多项式的NMM,达到了更高的精度;节点应力是连续的;定义在物理片上的所有自由度都具有明确的物理含义,其中第3～5个刚好是物理片所对应插值点处的应变分量,因此,直接获得此处的应力状态。该方法可以很容易地推广到其他基于单位分解的数值方法中。 相似文献
3.
为了改进计算区域离散化问题,本文利用自适应非结构化网格有限单元法求解二维地电结构下大地电磁场满足的加权余量表达式。在有限元求解电磁场的过程中,网格剖分越精细、计算精度越高,计算量也会越大。此外,结构化网格难以适应任意地形以及复杂地质构造。而自适应非结构化网格在电性变化剧烈的区域会自动加密,在电性缓变的区域则生成粗疏的网格,从而优化网格质量与数量。因此,文中引入COMSOL Multiphysics软件,以实现若干地电模型的构建及非结构化自由四边形单元网格化。将网格数据信息导入本文算法,计算大地电磁场响应,并与解析解及数值解对比。结果表明,基于非结构化网格的正演模拟精度高、适应性强,为计算区域网格化提供了新的方法。 相似文献
4.
常规的有限元强度折减法在确定边坡塑性区时,不能直观地反映出塑性区的形状和滑裂面位置,计算结果精度低。在塑性区范围内进行网格局部加密,不仅使安全系数的计算精度大幅度提高,而且可求出滑裂面的形状和位置。针对求解塑性区搜索滑裂面这一具体问题,提出塑性应变界限值网格加密准则,采用自适应性网格局部加密技术求解边坡的滑裂面和稳定安全系数。数值计算结果表明,所提出的加密准则能正确地确定滑裂面的形状和位置,计算所得稳定安全系数和理论稳定安全系数仅相差0.2 %。 相似文献
5.
6.
常规变加密网格有限差分波动方程数值模拟方法采用水平分层加密网格,该网格剖分策略在适应地形起伏和近地表速度结构变化特征方面效果较差。针对该问题,提出一种起伏多重变加密网格有限差分波动方程数值模拟方法。该方法根据地形起伏和近地表低速层到高速层的速度分布特征进行网格加密;采用不同网格中的变系数差分格式离散声波方程,在保证波场模拟精度的同时兼顾计算效率;同时,为了进一步保障地表附近波场模拟的精度,地表附近最细网格中的差分格式不做降阶处理,针对位于地表以上的虚像点的波场值,提出一种融入自由地表边界条件的法向虚像外推法。算例分析验证该算法对速度模型不同区域进行的网格多次加密显著提高了计算效率,以黄土塬实际模型为例,耗时为常规1 m×1 m网格耗时的43.3%,并可达到和细网格基本一致的模拟精度,模拟误差控制在10-12范围内,同时表现出很好的近地表散射压制和边界吸收效果,且算法能稳定地适应实际复杂地表介质。 相似文献
7.
8.
在岩土工程分析中求解精度控制常常是必需的,在数值流形法中可以通过控制数学覆盖网格的稀疏和覆盖位移的阶数来达到精度的要求。提出了基于等几何分析的数值流形方法,定义了相应的数学覆盖的构造形式,推导了基于二次B样条的9节点数值流形方法分析格式;针对基于Lagrange插值函数的4节点数值流形方法提出了基于T样条思想的数学覆盖网格的局部加密方法。算例计算结果表明,相对于4节点的数值流形方法,基于非均匀有理B样条的9节点数值流形方法具有更高的精度;基于T样条思想的加密网格在保持计算精度的前提下降低了自由度的数量,表明T样条加密是一种自然的局部加密算法。 相似文献
9.
天然水体具有形状复杂的计算区域和水下地形。本文采用无结构的网格以适合这一情况,并能方便地根据工程应用的要求局部地和适应性地加密网格。相应地,对二维水流计算问题采用了有限体积法的数学表示。跨单元边界的法向数值通量通过用Osher格式求解黎曼问题得到。文中给出了适合于二维浅水方程组的有关单元界面和各种物理边界的法向数值通量公式。分析了这一格式所具有的优点,包括:普适性、守恒性、逆风性、单调保持性、高效性、对间断的高分辨率、边界处理和内部格式相容,以及不引入数值边界条件等。最后,通过长江口南支水流计算的实例,阐明其良好性能。 相似文献
10.
《物探化探计算技术》2019,(5)
有限元求解大地电磁正演问题时,对研究区域的剖分常规做法是采用规则化的网格。但规则化网格在剖分地形、断层、褶皱等复杂模型会产生较大的几何离散误差。针对上述情况,这里采用非结构化的四边形网格对二维地电模型进行网格剖分,并与自适应有限元相结合,由剖分的粗网格出发,利用每个频点下网格单元的后验误差估计值指导网格的局部加密,优化网格质量和数量,从而提高正演模拟的精度。通过一维K型地电模型利用本文方法算出数值解与解析解进行对比分析,验证了自适应有限元法在求解大地电磁二维正演模拟中的有效性。并通过对断层、褶皱模型的正演模拟,分析了其大地电磁正演响应特征。 相似文献