共查询到14条相似文献,搜索用时 424 毫秒
1.
2.
3.
在岩土工程分析中求解精度控制常常是必需的,在数值流形法中可以通过控制数学覆盖网格的稀疏和覆盖位移的阶数来达到精度的要求。提出了基于等几何分析的数值流形方法,定义了相应的数学覆盖的构造形式,推导了基于二次B样条的9节点数值流形方法分析格式;针对基于Lagrange插值函数的4节点数值流形方法提出了基于T样条思想的数学覆盖网格的局部加密方法。算例计算结果表明,相对于4节点的数值流形方法,基于非均匀有理B样条的9节点数值流形方法具有更高的精度;基于T样条思想的加密网格在保持计算精度的前提下降低了自由度的数量,表明T样条加密是一种自然的局部加密算法。 相似文献
4.
针对岩土及地下工程结构的特点,基于修整后的位移型Reissner泛函中引入独立转动场的变分原理,转动场采用结点真实转角来插值。结合平面四结点单元讨论了有效附加非协调位移的合理形式,引入适用于任何四边形单元的非协调位移函数,建立了一种带转动自由度的平面四结点内参型非协调元模型,给出了相关数值计算列式。四结点单元能通过分片检验,并易于与带转动自由度的梁单元相容,对若干算例及地下隧洞工程实例计算的数值结果表明,采用所建的计算模型,具有较高的计算精度,为岩土及地下工程结构计算提供了一种理想的数值分析方法。 相似文献
5.
6.
7.
对三维问题的分析是数值流形方法发展的必然,在数值流形方法覆盖位移函数的基础上构造了一种六面体有限覆盖的三维流形单元,推导了相应的应变矩阵、刚度矩阵及平衡方程等表达式。同时,由于目前数值流形方法的模拟分析主要是采用线弹性模型,而对于非线性模型分析研究很少;根据数值流形方法的特点和岩土体的本构模型,给出了适用于数值流形方法进行非线性分析的算法。该方法利用中点增量法进行求解,以改变 模型和 模型中弹性模量的方式来反映非线性,其实质是用分段线性来取代非线性。通过地基沉降计算算例表明,数值流形方法在三维岩土体中进行非线性分析中是有效的。 相似文献
8.
数值流形方法中覆盖函数选用的建议 总被引:5,自引:0,他引:5
在数值流形方法中,常用的覆盖函数基并非是最佳选择。循着刚度矩阵的形成过程,分析了选用常用覆盖函数时,在非对角元上出现绝对值很大的元素之成因,且发现这会增加刚度矩阵的条件数,尤其是在用刚性弹簧约束位移的情况下,而这在数值流形方法中普遍而基本。建议采用局部化较好的覆盖函数,取代常用的关联于全局坐标的覆盖函数,可显著消除这一情况。建议方式简单明了,程序改动极小,对改善刚度矩阵性态却有很大作用。算例验证了这一建议的合理性,通过比较局部化的覆盖函数及全局性的覆盖函数所形成的刚度矩阵,表明前者形成了较小条件数的刚度矩阵。 相似文献
9.
提出了求解有自由面渗流问题的三维数值流形方法,通过构造任意形状流形单元的水头函数,推导了流形单元的渗透矩阵和无压渗流分析的总体控制方程,并给出了自由面的迭代求解策略和渗透体积力的计算方法。典型算例的数值分析表明,该方法采用数学网格覆盖整个材料区域,在自由面的迭代求解过程中数学网格保持不变,只考虑自由面以下渗流区的介质,只对自由面以下的流形单元形成总体渗透矩阵,具有精度高、收敛速度快、编程简单等优点,而且能够通过单纯形积分精确计算被自由面穿越单元的渗透作用力,因此,特别适用于有自由面渗流问题的模拟。 相似文献
10.
针对数值流形方法特有的覆盖剖分方式,提出了一种模拟岩土工程中开挖过程的算法。该算法采取某种措施,在覆盖剖分过程中将开挖面视为特殊的不连续面,这种不连续面将其所在的数学网格剖分成不同的流形单元,但却不对所在的数学覆盖作剖分。这样,开挖面两侧虽分属不同的流形单元,但开挖面两侧同一数学网格内的流形单元却具有相同的物理覆盖。采用该算法,无需对开挖面处的单元进行特殊处理,可在整个分析域采用统一的网格形式;同时,打破了原有数值流形方法的限制,将开挖面的位置完全当作连续介质来处理,避免了因将其视为不连续面而产生的误差。验证了算法的可靠性后,将其应用于某假想隧道的开挖模拟,计算结果表明该算法具有一定的应用前景。 相似文献
11.
弹簧元法是一种将单元离散为一系列弹簧的数值计算方法。不同的单元具有不同的离散方式,确定相应的离散弹簧的刚度系数表达式是弹簧元法的关键。将四节点矩形单元离散为6个基本弹簧,每个基本弹簧包括法向弹簧和切向弹簧两个派生弹簧,并用泊松弹簧和纯剪弹簧描述单元的泊松效应和剪切效应,用有限元的单元刚度矩阵标定各弹簧的刚度系数,实现了一种四节点矩形弹簧元的构造形式。该单元的同类弹簧具有相同的表达形式。法向与切向弹簧的刚度表达式中分别含有法向和切向弹簧刚度待定系数。通过改变待定系数的值可使该单元分别对应于有限元的常应变、双线性及Wilson非协调单元。将上述弹簧元方法进行理论推导,并应用于基于连续介质的离散单元法(CDEM)的核心计算进行简单算例验证,证明了提出方法的正确性。通过以上研究发现,四节点矩形弹簧元有以下特点:对于相同问题,不同单元有不同的计算精度;对于梁弯曲问题,应用该单元可显著提高离散单元法的求解精度;改变待定系数的值,可得到更高或者更低精度的单元。 相似文献
12.
数值流形法(numerical manifold method)是一种新型的数值计算方法,已成功应用于岩土工程的诸多领域,但该方法尚未应用于岩土工程蠕变分析。近年来对高阶流形法的研究表明,对复杂的岩土工程问题,使用高阶覆盖函数可明显提高流形法的计算精度。为此,开展了用高阶流形法模拟蠕变的研究,在高阶流形法中引入“时步-初应变”法计算蠕变,以广义开尔文体为基础,推导了相关的计算公式,并编制了相应的计算程序,同时还通过算例,验证了方法的可行性和合理性。结果表明,高阶流形可以方便地与“时步-初应变”法结合用于蠕变计算,可较好地模拟蠕变变形。算例分析表明,在不改变网格密度情况下,仅通过采用高阶覆盖函数,高阶流形法可大幅提高传统流形法的计算精度。 相似文献
13.
To improve the computational efficiency of the numerical manifold method for discontinuous deformation simulations, a spatial-domain coupled explicit-implicit time integration algorithm is proposed. A subdomain partition algorithm based on a super manifold element is developed for the numerical manifold method to simulate dynamic motions of blocky rock mass. In different subdomains, explicit or implicit time integration method is employed respectively based on its contact and motion status. These subdomains interact through assembling the corresponding explicit or implicit time integration-based matrices of different rock blocks. The computational efficiency of the discontinuity system under dynamic loading is improved by partially diagonalizing the global matrices. Two verification examples of a sliding block along an inclined plane under a horizontal acceleration excitation and a multiblock system acted on by dynamic forces are studied to examine the accuracy of the proposed numerical method, respectively. A highly fractured rock mass situated on an inclined slope subjected to seismic excitations is then studied to show the computational efficiency of the developed algorithm. The simulated results are in good agreement with those from the versions using purely implicit or explicit time integration algorithm for the numerical manifold method. The computational efficiency is shown to be higher using the proposed algorithm, which demonstrates its potential for application in dynamic analysis of highly fractured rock masses. 相似文献