共查询到20条相似文献,搜索用时 15 毫秒
1.
2.
3.
承压含水层中地下水向井非稳定流动的LT有限分析法 总被引:1,自引:0,他引:1
作者给出了一种求解承压非稳定井流方程的LT有限分析法.其特点是:将时间t直接包含在有限分析数值计算格式中,克服了传统数值法求解非稳定过程按时段迭代的缺陷,实现了非稳定过程计算的一步完成,给多时段地下水资源管理和长期预报带来了方便.用本文提出的方法分别计算了两个算例,结果表明,该方法具有较高的精度,与解析解吻合甚好,可用于实际地下水流的计算之中. 相似文献
4.
针对岩土工程材料应变软化问题及有限元对其数值计算时切线刚度矩阵负定造成求解困难的问题进行研究。建立了基于Drucker-Prager(D-P)强度准则的岩石弹塑性应变软化本构模型,本构积分算法采用一种完全隐式返回映射算法,它具有无条件稳定和精确的特点,详细论述了如何进行本构模型的程序化求解;考虑弧长法在判断切线刚度矩阵正定性导致效率低的缺点,在弹塑性增量有限元方程的迭代计算中尝试采用Newton-Raphson法和arc-length法(NR-AL法)联合迭代求解的思路,即在结构未达到极限荷载前采用NR迭代法,而当结构接近极限荷载时转换为AL法控制迭代,从而使结构越过峰值点进入软化区直至破坏,NR-AL法汲取了2者迭代求解中具有的优势;利用C++语言对所建应变软化模型的本构求解和弹塑性增量有限元方程迭代求解过程给予程序实现,应用所编程序进行数值计算,分析了D-P理想弹塑性模型、应变软化模型、应变硬化模型计算的应力-应变曲线的区别,同时将应变软化模型计算结果与试验数据进行了对比。研究结果表明:所建应变软化本构模型可以较好地模拟岩石材料的峰后软化特性,能够揭示峰后应变软化特性和破坏机制,同时NR-AL法能够求解由于应变软化造成的负刚度问题,也克服了单独使用弧长法时判断切线刚度矩阵正定性效率低的缺点。 相似文献
5.
求解钠长石在碱长石和斜长石中的分配问题常采用里特曼相图法,但其计算量大,计算过程复杂,易受人为因素影响.利用Matlab软件编程模拟人工查图时的思维和工作步骤,实现里特曼相图法查图,包括相图数字化、相图分区、混溶区内样品点的求解、不混溶区内样品点的求解和求解最终结果等,并给出Matlab语言编程中部分计算问题的关键代码.实例计算表明,手工计算约需50 min的计算量,本方法可在2~3 s内完成,准确率为100%. 相似文献
6.
7.
8.
9.
10.
在三维电阻率的正演计算中往往涉及到快速、准确求解大型线性方程纽Ax=b的问题。通过采用有限差分法来构造出求解点电源三维地电场的大型稀疏对称线性方程组。并引入Lanczos迭代技术,构造出三对角阵方程组,然后采用正交分解法进行求解,它是Krylov子空间方法中的一种。与传统迭代算法相比,它占用内存少,收敛速度快且稳定。针对大型稀疏矩阵及MATLAB语言的特点,采用简单记录矩阵的非零元素值及其所在行、列值的方法存储大型稀疏矩阵,可大大节省机器内存,提高运算速度。理论分析和计算实例显示,此算法是地电三维正演计算的有效方法,为下一步的反演计算打好基础。 相似文献
11.
在建立地下水流模型的过程中,模型验证一直是较为复杂的步骤之一,具体难点包括寻优方法的选用,为保持总体平衡所引起的参数峰值异常以及总体寻优需要大量的计算机时等问题。本文运用区域分解法中D-N交替法的基本思想,在计算目标函数值时,将整体求解水文地质模型的过程分解为计算各参数分区内的子模型的过程,并在双核计算机上实现了并行计算。理想算例的计算结果证明了该方法用于水文地质逆问题求解的可行性,它不仅减少了求解过程中所需要的计算机时,而且提高了参数的拟合度。最后将这种方法应用到山东鹏山水源地的水文地质参数求解问题中,由结果可以看出,运用该方法反求水文地质参数是可行的,具有很好的应用前景。 相似文献
12.
根据水质模型的具体特点,对不同的方程采用不同方法,水流问题用有限元法;对流弥散方程先用算子分裂的方法分解为两个方程,即对流方程和弥散方程,前者用高精度广义迎风格式求解,对弥散方程则采用多单元均衡格式法求解,最后合成为高精度广义迎风均衡格式求出溶质浓度。通过对数值实验例子的计算和实验溶质迁移的模拟,可以看出在求解对流弥散定解问题时,广义迎风均衡格式克服了有限元数值波动和浓度出现负值的问题,与有限元相比有较大改进。 相似文献
13.
14.
三维地面断面电阻率测深有限元数值模拟 总被引:6,自引:0,他引:6
用有限单元法进行了电导率分块均匀的三维点源电场电阻率测探的数值模拟。首先给出了三维构造中点源电场的边值问题、变分问题;然后,用有限单元法求解变分问题,将区域剖分成六面体单元,在单元中进行三线性函数插值,将变分方程化为线性代数方程组;最后解方程组,得各节点的电位值,进而计算出地表的视电阻率。对几例较典型的地电模型进行试算,结果表明本方法是行之有效的。 相似文献
15.
在地震勘探中,地震正演模拟是非常重要的技术。与时间域正演相比,频率域正演速度快,计算效率高。如何高效准确地完成频率域正演计算是目前该领域的一个重要问题。数值频散问题和如何提高计算效率降低求解分解阻抗内存占用量一直是频率域正演所需要解决的问题。与传统的直接法求解阻抗矩阵的频率域正演方法不同,本文采用最小二乘共轭梯度法(LSCG法)求解阻抗矩阵进行频率域正演,并提出了一种波数补偿的表达式来压制数值频散现象。经过简单模型和复杂模型的数值测试,采用最小二乘共轭梯度法(LSCG法)求解阻抗矩阵进行频率域正演能够有效降低计算时间,且采用波数补偿的频率域正演方法能够有效压制数值频散现象,提高波场模拟精度。 相似文献
16.
本文讨论了在均匀谐变场激发下导电导磁体电磁响应的边界单元法计算问题。文中首先建立矢量位的边值问题,然后导出其边界积分方程,最后用边界元法进行求解。本方法的计算结果具有精度高、速度快的优点。 相似文献
17.
有自由面渗流分析的高斯点法 总被引:11,自引:0,他引:11
有自由面的渗流问题是比较复杂的非线性问题,这类问题的有限元求解一般采用固定网格法。本文根据饱和区和非饱和区的渗流特点,提出了高斯点法,该算法原理简单,易于实现。 相似文献
18.
19.
岩土工程数值计算中的无网格方法及其全自动布点技术 总被引:14,自引:3,他引:11
自然单元法采用无网格的思想全域构造插值函数,它的求解精度高,计算时间少,可准确地施加边界条件,兼具有无网格法和有限单元法的优点和特点,是一种理想的用于岩土及地下工程分析计算的数值方法。文中简要地介绍了自然单元法的基本理论,并针对岩土及地下工程问题特点,给出了一种无网格离散点的全自动布置方法。 相似文献