首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 609 毫秒
1.
多重网格有限元法在水质预测问题中的应用   总被引:4,自引:0,他引:4  
本文采用多重网格—有限元法进行模拟、预报。多重网格法是求解大型科学计算问题的有效且可以证明是最优阶的方法,但目前很少有人将其用于解决实际问题,且在理论上也局限于二维问题的讨论。本文首次利用多重网格法求解地表水地下水水质耦合问题,从而拓宽了多重网格法的实际应用范围;与实际观测值对比有较高的拟合精度,有效地避免了数值弥散和振荡;且具有计算速度快、收敛稳定性好等优点   相似文献   

2.
梅发国  高彦伟  郭华 《世界地质》2002,21(4):385-389
针对差分法求解过程中细网格计算所需时间和空间要求高以及在粗网格上计算又达不到所要求的精度等问题,利用粗网格上的残差校正特性消除迭代误差的低频分量,同时利用细网格上的松弛光滑特性消除迭代误差的高频分量,采用不同疏密的网格距消除不同频率范围内的误差分量,将多重网格方法成功应用于频域波动方程的求解,并给出了实际算例。算例中,分别用Gauss-Seidel迭代法和多重网格法求其数值解,并绘出效果图;比较得出多重网格方法具有精度高、收敛速度快和易于实现的特性,更适合于高维波动问题的求解。  相似文献   

3.
分析了求解变分不等式Newton方法的半局部收敛性,建立了类似于Kantorovich定理的收敛性结果。该结果不仅为判断Newton方法的收敛性提供了可计算的充分条件,也给出了Newton方法的收敛域以及问题解的存在区域。同时,文章也得到了Newton方法的若干收敛性质,包含收敛阶以及可计算的误差估计式等。  相似文献   

4.
ρ(L_(r。w))的界和迭代法收敛性刘光清,辛迪远(长春地质学院基础科学系.长春130026)(四平电视大学.四平136000)本文给出了广义AOR法(广义SSOR,广义SOR,广义G-S及广义西法)迭代矩阵谱半径较理想的上界及其收敛准则,推广了文?..  相似文献   

5.
利用Matlab函数识别沉积物中的米兰柯维奇旋回信号   总被引:1,自引:0,他引:1  
针对旋回分析中需要对数据进行的各种处理,利用Matlab中集成的滤波(filter)函数、插值(interp1)函数、求均值(mean)函数、去线性趋势(detrend)函数、求方差(variance)函数及谱分析(pwelch)函数等的组合,编制了识别沉积物中米兰柯维奇旋回信号的程序。选择窗口宽度(window)约等于有用信号中最低频率信号的波长、交叠点数(noverlap)为窗口宽度的50%、傅里叶变换的数据点数(nfft)大于全部采样数据的长度且为2的最小整数次幂时,可获得较好的谱估计结果。其不足之处是不能估计背景的红噪声谱及给出谱峰值的置信水平。  相似文献   

6.
赵明君 《探矿工程》2010,37(5):10-13
井眼轨道的软着陆设计模型的求解可以归结为一个七元非线性方程组的求解问题。前人给出了数值迭代求解算法,然而并没有证明该迭代算法的收敛性,并且该算法是否收敛严重依赖于用户给出的迭代初始值。通过一系列的消元、化简的数学技巧,将七元非线性方程组化简为一元多项式方程,并在此基础上给出了软着陆设计模型的一个新算法。理论分析和实际算例表明,新算法的主要计算工作量是求多项式方程的非负实数根,其他未知数与实数根是简单的函数关系,计算量很小。新算法克服了迭代算法的初值依赖性以及迭代过程可能发散等缺陷,并且在设计模型有多个解的情况下,可以同时求出这些解。  相似文献   

7.
大地电磁反演问题通常表述为目标函数最优化,难点是多参数、非线性和不适定性,局部和全局方法都不能实现快速全局优化[4].针对局部线性方法易使解陷入局部极值,严重依赖初始模型,而传统的遗传算法在优化应用中存在局部搜索能力弱、早熟收敛等问题.这里引进一种求解一维大地电磁测深反演问题的实数编码广义遗传算法.该算法利用拟网格法初始种群和综合交叉策略,克服了早熟收敛现象,从而提高了遗传寻优的效率.理论模型反演与其它方法比较,结果说明遗传算法具有不依赖初始模型,不容易陷入局部极小,多点多路径概率搜索,以及隐合并行性等优点.  相似文献   

8.
谱元法(SEM)是基于有限元发展起来的一种数值模拟方法,具备有限元法处理复杂结构的几何灵活性和谱方法的高精度和指数收敛性,将其用于数值模拟计算比传统高阶有限元减少了内存需求和CPU时间。为了提高大地电磁数值模拟的精度和效率,基于Galerkin加权余量法,采用Gauss-Lobatto-Legendre(GLL)正交多项式作为插值基函数,用Intel的Pardiso求解器直接求解线性方程组得到电磁场分布,实现了大地电磁(MT)二维谱元法数值模拟。给出的层状模型和COMMEMI 2D-1模型的数值算例验证了谱元算法和所编写程序的可靠性和稳定性。  相似文献   

9.
浸润线的计算是尾矿坝安全状况评估的重要环节,对尾矿库安全起到至关重要的作用。为了能够用解析的方法得出浸润线,首先略掉无关紧要的部分,将尾矿坝的渗流分析区域简化为一个矩形,然后在矩形区域中选定两种符合工程实际的边界条件,最后求解渗流微分方程即求出这两种边界条件下的浸润线。求解中涉及到的偏微分方程无论是齐次还是非齐次的,都通过傅里叶级数法得到收敛的无穷级数解。该方法的突出特点就是得到的最终结果是解析解。实际应用时通过与实测浸润线数据拟合确定浸润线求解公式中的常数。通过对某尾矿坝浸润线的计算发现前面得到的公式可以用于实际浸润线的计算。  相似文献   

10.
Laplace方程Cauchy问题的一种数值解法   总被引:2,自引:1,他引:2  
研究了Laplace方程Cauchy问题的数值求解,该问题是一个典型的病态问题。利用格林(Green)公式将Laplace方程的Cauchy问题转化为Hausdorff矩问题。本文利用一种新方法,即矩问题的积分方程方法,求解矩问题,设计了二维Laplace方程Cauchy问题稳定的算法,给出了近似解的误差估计,并对二维Cauchy问题进行了数值模拟。  相似文献   

11.
基于三维非线性有限元的边坡稳定分析方法   总被引:3,自引:0,他引:3  
刘耀儒  杨强  薛利军  周维垣 《岩土力学》2007,28(9):1894-1898
刚体极限平衡法不能反映岩体中实际的应力分布,而基于有限元的强度折减系数法在判断收敛性方面存在一些问题。为了解决这些问题,采用多重网格法,分别建立用于有限元计算的结构网格和用于计算滑面稳定安全系数的滑面网格,可以方便地获得任意滑面或滑块的稳定安全系数,从而将非线性有限元和极限平衡分析结合起来。为了提高计算规模和计算精度,采用有限元并行计算程序TFINE.Pfem进行计算,分析了网格密度对计算结果精度的影响,并应用于锦屏高边坡的稳定分析中。与刚体极限平衡法结果的对比分析表明,由于考虑了计算过程中的非线性应力调整,该方法的计算结果比刚体极限平衡法偏大,而且更符合实际情况。  相似文献   

12.
传递函数分析在边坡振动台模型试验的应用探讨   总被引:4,自引:1,他引:3  
蒋良潍  姚令侃  吴伟  徐光兴 《岩土力学》2010,31(5):1368-1374
传递函数分析对获取振动系统的动力特性参数及估算动力响应具有重要意义。边坡地震模拟振动台模型试验可求算两种传函:反映边坡各点的场地地震动特征的绝对加速度传函(基于测点实测加速度时程)和反映坡体本身振动性质的相对加速度传函(基于测点与坡脚即台面间的实测加速度相对时程)。以弹性体振动方程分析和试验实测,探讨了边坡振动台模型试验中动力特性参数计算的传函选择和以传递函数对动力响应的估算,得出结论:(1)不同于机械模态分析(以绝对传函的模和虚部计算动力特性参数具有等价性),对于边坡模型试验,绝对传函实部和模混入了台面加速度影响,以绝对传函虚部或相对传函虚部才能正确计算动力特性参数;(2)传递函数和台面输入波形Fourier谱的逐频相乘得到动力响应估算谱,与测点实测响应Fourier谱吻合较良好,这种估算方法有助于对不同地震波作用下的边坡动力响应作出总体评价,并指导模型试验中地震波形的合理选取。  相似文献   

13.
杨强  朱玲  薛利军 《岩土力学》2008,29(1):94-100
在以有限元法来分析评价坝基、边坡的稳定安全度时,多重网格法将有限元分析网格和滑移块体的滑面网格作为独立的两套网格,滑面网格可随意切割有限元单元。这样一套有限元网格和一套应力成果就可以进行任意的滑块稳定分析。着重分析了多重网格法的精度。针对三维柱体和锦屏高边坡,研究两套网格如何匹配以提高计算精度进行处理,初步得出了一些有关多重网格法计算精度的知识。最后,应用基于多重网格法的第2类方法分析了锦屏左坝肩高边坡的稳定性问题。  相似文献   

14.
The real spectral analysis of SP anomalies due to 2-D horizontal circular cylinder is carried out using Hartley transform. The Hartley transform is an alternate means of realizing spectral analysis in real domain unlike the Fourier spectral analysis. The Hartley transform yields a straightforward interpretation of SP anomalies caused by horizontal circular cylinder wherein all the parameters are derived independently as a function of frequency. A theoretical example illustrates the procedure. The effect of random noise, up to 10% of white Gaussian noise, on the interpretation scheme was also studied and found to be of negligible. The field example of the ??Sulleymonkey?? anomaly in the Ergoni copper district, Turkey exemplifies the applicability of the proposed method.  相似文献   

15.
A method is derived to estimate the degree of consolidation of a saturated, non-homogeneous clay layer. In this method, the diffusion equation govering the dissipation of pore water pressure is transformed so that the pore water pressure can be expressed by a Fourier series. The Fourier coefficients are then evaluated along the trajectory of the diffusion process. The Fourier series does not converge point-wise to the pore water pressure; instead, it converges to the pore water pressure in the mean., The Fourier series can be used to estimate the degree of consolidation, because, by Schwartz's inequality, convergence in the mean implies convergence in the degree of consolidation.  相似文献   

16.
本文给出了瞬变弹性动力问题的一个边界元法。该法是利用威尔逊——θ法的差分公式,把运动方程化为椭圆型微分方程。根据贝蒂定理和动力点荷载的特解,可获得动力问题的边界积分方程。这个解法是在真实时间域内逐步求解的,不需要使用拉氏变换。一个应力波传播的数值算例证实了该方法使用方便,且解答精度较高。  相似文献   

17.
A boundary element method formulation in the real time domain in the transient elastodynamic problems is presented. The equation motion of the linear isotropic elastiic bodare transformed into the elliptic partial differential eqnations by means of the differential formulations of the Wilson-θ method. By utilizing the particular solution corresponding to a dynamic point force in an infinite medium and the Betti’s theorem the bundary integral equations may be obtain. As this approach formulates and solves the problem in real time domain in conjunction with a time-stepping algorithm, it does not require a nimerical inversion.Which transformes the results from the Laplace or Fourier ransformed domain to the time domain. A numerical example was shown in order to demonstrate the applicability of the present method to the practical problems.  相似文献   

18.
针对向下延拓的不适定问题,提出了一种新的位场向下延拓方法-导数迭代法。从位场垂向一阶导数的定义出发,将观测面和向上延拓、向下延拓同等高度平面上的位场值近似联系起来,采用迭代法进行逐次逼近,推导出空间域位场向下延拓的导数迭代法的递推公式。考虑到空间域向上延拓积分方程实现的复杂性,对递推公式进行快速傅里叶变换,整理得到波数域中的导数迭代公式,同时证明了该方法的收敛性。模型检验和实例分析均表明迭代法相对直接FFT法具有稳定性强和下延深度大的优点。  相似文献   

19.
波动方程角度域共成像道集   总被引:1,自引:0,他引:1  
陈生昌  马在田 《地球科学》2007,32(4):569-574
针对基于波场深度递推的波动方程叠前深度偏移成像方法难以准确地给出如Kirchhoff积分叠前深度偏移成像方法的偏移距域共成像道集的不足,在分析研究现有的各种波动方程偏移成像共成像道集方法的基础上,利用波场外推的单平方根算子和波场的窗口Fourier框架展开与重构方法,提出一种局部角度域共成像道集方法——成像点处反射角共成像道集方法.把这种角度域共成像道集方法应用于国际标准的Marmousi模型数据和一条实际二维地震数据都取得了理想的结果.提出的波动方程反射角度域共成像道集方法可为进一步的叠前偏移数据振幅随角度变化分析和偏移速度分析工作提供基础.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号