共查询到19条相似文献，搜索用时 640 毫秒
1.
有限元直接迭代算法中一种基本结构的研究 总被引：4，自引：0，他引：4
陈小斌 《物探化探计算技术》2000,22(2):132-136
首先探讨了基本结构的含义,提出笛卡尔直角坐标系下具体的基本结构,讨论了该基本结构的各种剖分方式;然后,基于有限元直接迭代算法原理,具体推导了其中一种剖分方式下的迭代格式;最后在正方形网格部分情形下,以拉普拉斯方程为例,分析了各部分方式所确定的迭代格式的截断误差。 相似文献
2.
为了对二维复杂地电断面进行三角网格剖分,在已有算法的基础上,提出了一种含约束条件的二维自适应三角剖分新算法,并编制了该算法的通用程序,可对任意含约束条件的二维区域进行自动剖分.通过引入节点剖分尺度函数,算法允许在指定位置设置网格剖分尺寸,自动生成渐变网格.实际算例表明,形态因子大于0.9的网格占95%以上,能很好地对带地形、含断层或地质体等二维复杂地电断面进行有限元三角网格剖分. 相似文献
3.
为了更好地模拟地下介质连续变化及开展连续介质的反演,对二维电导率分块线性变化的线源频率域可控源电磁法进行了有限元正演模拟,在剖分单元内同时对电场及电导率参数线性插值,使电导率参数在剖分单元之间保持连续变化。首先,提出有限元正演模拟的边值问题及变分问题,并详细论述了有限元的剖分、插值、单元分析及总体合成的各个步骤;其次,采用稀疏存储及基于不完全LU分解的BICGSTAB算法求解复系数方程组,节省了内存并提高了计算速度;然后,对一个均匀半空间模型进行模拟,计算结果表明,低频及高频的有限元数值解都与解析解吻合,证明了算法的正确性;最后,对水平层状模型及垂直断层模型进行正演计算,视电阻率及相位的等值线图均较好地反映出了异常体,说明文中算法能够对电导率连续变化的线源可控源电磁法进行有效地模拟。 相似文献
4.
MT二维正演的混合剖分算法及其在线源CSAMT正演计算中的应用 总被引：1，自引：1，他引：1
在矩形有限元和三角形有限元方法大地电磁测深（ＭＴ）二维正演计算的基础上，设计了矩形单元粗剖分、程序自动进行三角形单元细剖分的混合剖分算法程序，使得ＭＴ二维模型正演的计算速度大为提高。该方法同时被应用于线电流源变频电磁测深（ＣＳＡＭＴ）ＴＥ极化方式的二维正演计算。文章给出了在不同模型上计算的ＭＴ和线源ＣＳＡＭＴ的响应，分析了计算ＣＳＡＭＴ响应时，计算点离源的距离对视电阻率和相位值计算的影响，比较了ＭＴ和ＣＳＡＭＴ网格剖分的特点，提出了计算ＣＳＡＭＴ响应的改进的边界条件设计 相似文献
5.
针对地下水有限元数值模拟中区域三角网格剖分复杂难以处理的情况,提出适合其特点的Delaunay三角网格自动剖分方法,并对含有多个参数分区的含水层进行网格剖分,最后利用遗传算法和有限元程序相耦合来反演含水层水文地质参数。结果表明此方法可大大简化地下水数值模拟的前处理工作,并能提高有限元网格剖分的有效性和准确性,从而得到令人满意的数值模拟结果。 相似文献
6.
7.
时间域激发极化法非结构化三角网格有限元正演模拟 总被引：1，自引：0，他引：1
这里分析了传统网格剖分方式的不足,提出采用非结构化渐变三角网格对求解区域进行离散,采用固定带宽LDLT分解法求解线性方程组。设计了非结构化渐变三角网格自动剖分程序,并针对固定点源测深装置的激发极化2.5维有限元正演模拟程序。通过理论模型试算,验证了程序的正确性与有效性,并对正演模拟的速度和精度进行了初步评价。 相似文献
8.
王文科 《西安地质学院学报》1996,18(3):56-62
本文给出了一种在不同规则四边形网格上求解地下水非稳定流方程的混合拉氏变换局地标有限分析格式。由于将时间ｔ直接包含在有限分析数值计算格式中，克服了传统数值法求解非稳定过程按时段迭代的缺陷，实现了非稳定过程计算的一步完成，提高了计算效率，这对多时段地下水资源管理和长期预报非常有用。 相似文献
9.
余涛朱自强鲁光银曹书锦 《物探化探计算技术》2013,(4):446-451
重力张量是重力位的二阶空间导数,对密度体的变化和细节部份反映更为灵敏。对于复杂的重力密度模型来说,张量的解析公式是很难推导的。为了模拟这种情况,将复杂模型进行有限元剖分,计算每个单元的重力异常对测点张量的影响,最后叠加得到整个复杂模型在测点处的重力张量。采用Delaunay四面体非结构化网格对密度体进行剖分,经分析表明,这种剖分方式具有较高的精度。有限元剖分的方式拟合复杂形体缺点之一是计算量比较大,而并行化能有效解决这个问题。这里基于OpenMP并行模型编写了重力张量并行正演程序,分析了不同情况下并行的执行性能,为大规模复杂模型的重力张量正演提供了一种并行策略和思路。 相似文献
10.
二维地形不平条件下均匀外电场的有限差分模拟 总被引：3，自引：0，他引：3
本文提出了一种求解二维不平地形条件下中梯装置电阻率法正演问题的有限差分算法。引入了数学上的某些新进展:采用三角形和矩形相结合的方式对求解区作剖分,以积分插值法导出有限差分格式,用多层网格修正迭代法求解线性方程组,从而得到了比较满意的模拟精度和计算速度。 相似文献
11.
混沌控制反演系统构造及算法逻辑设计 总被引：1，自引：0，他引：1
给出了适用于非线性反演的混沌控制反演系统构造方法及求解控制矩阵的算法逻辑设计。该方法在迭代反演控制参数和迭代反演输出结果之间建立耦合关系，通过时时修改控制参数保证迭代稳定收敛到预期的解空间。在求解过程中，应用数据结构方法分析了混沌控制理论中计算控制矩阵的数据结构，采用树形结构表述出控制矩阵计算过程中下标的取值逻辑，利用可以复读栈中元素的中序遍历来遍历每一棵树，寻找各元素的组合序列，给出的数值算例说明了混沌控制反演方法及其算法的有效性。 相似文献
12.
由于计算机在地质研究中的大量运用,传统地质图件的数字化成为目前许多应用的主要 数据来源。作者在文章中主要讨论利用数字化仪输入地质图件的基本原理,数据转换以及坐标变换的核心技术和在Ｗｉｎｄｏｗｓ环境下,如何用Ｖｉｓｕａｌ Ｂａｓｉｃ６．０开发数字化仪通信程序。最后给出了一个简单的数字化仪通信实验程序。 相似文献
13.
Calculations of isochrons by different methods are compared for the case when both variables are subject to error. Least-squares operations of the kind given by York (1969) and linear methods are compared. The advantage of the linear methods is demonstrated for the case of data reported by Housner and Brennan (1948), which we have modified for nonreplicated sampling. The sampling distribution of the slope is obtained and confidence intervals for it are given. The confidence intervals for the intercept on the ordinate and its distribution are obtained for fixed values of the slope. Examples of application of the method are given, based on data from the literature compiled in Appendix A. 相似文献
14.
研究了带有约束机制的次优状态反馈律，主要对终端时间tf=∞系统建立次优状态反馈的数学模型。通过简单控制结构约束的引入，消除状态变量不完全可测所造成的障碍，从而得到了带控制结构约束的次优必要条件及改进后的迭代算法。 相似文献
15.
常规DCD(dynamic canonical descent)算法具有全局优化能力且无需考虑目标函数的可微性,只要预先定义优化空间即可,但是该算法的收敛速度有限,为此提出了变参数DCD算法,并建立了其相应的算法迭代格式。在数值试验和工程应用中将该变参数DCD算法与常规的DCD算法进行比较,其结果均表明：变参数DCD算法在全局优化能力和收敛速度上找到了一个均衡点,该算法不仅具有DCD算法的全局优化能力,而且收敛时所需的目标函数评估次数少,在优化过程中该算法展示出了稳定性强且优化结果可靠度高的一面。 相似文献
16.
Seismic inverse modeling, which transforms appropriately processed geophysical data into the physical properties of the Earth, is an essential process for reservoir characterization. This paper proposes a work flow based on a Markov chain Monte Carlo method consistent with geology, well-logs, seismic data, and rock-physics information. It uses direct sampling as a multiple-point geostatistical method for generating realizations from the prior distribution, and Metropolis sampling with adaptive spatial resampling to perform an approximate sampling from the posterior distribution, conditioned to the geophysical data. Because it can assess important uncertainties, sampling is a more general approach than just finding the most likely model. However, since rejection sampling requires a large number of evaluations for generating the posterior distribution, it is inefficient and not suitable for reservoir modeling. Metropolis sampling is able to perform an equivalent sampling by forming a Markov chain. The iterative spatial resampling algorithm perturbs realizations of a spatially dependent variable, while preserving its spatial structure by conditioning to subset points. However, in most practical applications, when the subset conditioning points are selected at random, it can get stuck for a very long time in a non-optimal local minimum. In this paper it is demonstrated that adaptive subset sampling improves the efficiency of iterative spatial resampling. Depending on the acceptance/rejection criteria, it is possible to obtain a chain of geostatistical realizations aimed at characterizing the posterior distribution with Metropolis sampling. The validity and applicability of the proposed method are illustrated by results for seismic lithofacies inversion on the Stanford VI synthetic test sets. 相似文献
17.
A challenging computational problem arises when a discrete structure (e.g. foundation) interacts with an unbounded medium (e.g. deep soil deposit), particularly if general loading conditions and non‐linear material behaviour is assumed. In this paper, a novel method for dealing with such a problem is formulated by combining conventional three‐dimensional finite‐elements with the recently developed scaled boundary finite‐element method. The scaled boundary finite‐element method is a semi‐analytical technique based on finite‐elements that obtains a symmetric stiffness matrix with respect to degrees of freedom on a discretized boundary. The method is particularly well suited to modelling unbounded domains as analytical solutions are found in a radial co‐ordinate direction, but, unlike the boundary‐element method, no complex fundamental solution is required. A technique for coupling the stiffness matrix of bounded three‐dimensional finite‐element domain with the stiffness matrix of the unbounded scaled boundary finite‐element domain, which uses a Fourier series to model the variation of displacement in the circumferential direction of the cylindrical co‐ordinate system, is described. The accuracy and computational efficiency of the new formulation is demonstrated through the linear elastic analysis of rigid circular and square footings. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
18.
This paper develops a general abstract framework for a posteriori estimates for immiscible incompressible two-phase flows in porous media. We measure the error by the dual norm of the residual and, for mathematical correctness, employ the concept of global and complementary pressures in the analysis. Our estimators allow to estimate separately the different error components, namely, the spatial discretization error, the temporal discretization error, the linearization error, the iterative coupling error, and the algebraic solver error. We propose an adaptive algorithm wherein the different iterative procedures (iterative linearization, iterative coupling, iterative solution of linear systems) are stopped when the corresponding errors do not affect significantly the overall error and wherein the spatial and temporal errors are equilibrated. Consequently, important computational savings can be achieved while guaranteeing a user-given precision. The developed framework covers fully implicit, implicit pressure–explicit saturation, or iterative coupling formulations; conforming spatial discretization schemes such as the vertex-centered finite volume method or the finite element method and nonconforming spatial discretization schemes such as the cell-centered finite volume method, the mixed finite element method, or the discontinuous Galerkin method; linearizations such as the Newton or the fixed-point one; and general linear solvers. Numerical experiments for a model problem are presented to illustrate the theoretical results. Only by stopping timely the linear and nonlinear solvers, speedups by a factor between 10 and 20 in terms of the number of total linear solver iterations are achieved. 相似文献