有限元法与伪谱法混合求解弹性波动方程   总被引：6，自引：0，他引：6  

马德堂  朱光明 《地球科学与环境学报》2004,26(1):61-64

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.  相似文献   

Effect of element size on the static finite element analysis of steep slopes     

Scott A. Ashford  Nicholas Sitar 《国际地质力学数值与分析法杂志》2001,25(14):1361-1376

The accuracy of the computed stress distribution near the free surface of vertical slopes was evaluated in this study as a function of the element size, including aspect ratio. To accomplish this objective, a parametric study was carried out comparing stresses computed using the finite element method (FEM) to those obtained from a physical model composed of photoelastic material. The results of the study indicate a reasonable agreement between a gelatin model and the FEM model for shear stresses, and an overall good agreement between the two models for the principal stresses. For stresses along the top of the slope, the height of the element tends to be more important than width or aspect ratio, at least for aspect ratios up to 4. In all cases, the greatest difference between the two models occurs in the vicinity of the slope. Specifically, if H is defined as the slope height, an element height of H/10 appears to be adequate for the study of stresses deep within the slope, such as for typical embankment analyses. However, for cases where tensile stresses in the vicinity of the slope face which are critical, such as for the stability analysis of steep slopes, element heights as small as H/32, or higher‐order elements, are necessary. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

A 3D finite element simulation model for TBM tunnelling in soft ground     

Thomas Kasper  Günther Meschke 《国际地质力学数值与分析法杂志》2004,28(14):1441-1460

  相似文献   

广西沙田柚农业地质   总被引：1，自引：0，他引：1  

冯群耀 《南方国土资源》2001,14(1):47-50

沙田柚生长好坏、产量品质高低 ,与所在地质环境有密切关系 ,种植在碎屑岩和花岗岩环境区及河流阶地、江心洲、山麓地貌环境区的沙田柚 ,长势好、产量品质高。种植在碳酸盐岩环境区及山坡型、岩溶洼地地貌区的沙田柚 ,长势差、产量品质低。影响沙田柚长势、产量品质的关键元素是Na,Na是沙田柚的特需元素。  相似文献   

成都市冬季重污染日气象条件成因分析及预报     

张波  杨柏生 《成都信息工程学院学报》1992,(3)

本文通过对成都10个重污染日进行天气学分析,将污染浓度与气象要素进行聚类、研究了重污染日形成原因及污染浓度与气象要素的关系。在此基础上,建立了SO_2、TSP日平均浓度分级预报方程。  相似文献   

A micromechanical finite element model for linear and damage‐coupled viscoelastic behaviour of asphalt mixture     

Qingli Dai  Martin H. Sadd  Zhanping You 《国际地质力学数值与分析法杂志》2006,30(11):1135-1158

This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

Fast iterative solution of large undrained soil‐structure interaction problems     

Kok‐Kwang Phoon  Swee‐Huat Chan  Kim‐Chuan Toh  Fook‐Hou Lee 《国际地质力学数值与分析法杂志》2003,27(3):159-181

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

REVIEWS OF GEOGRAPHIC SOFTWARE     

《The Professional geographer》1987,39(1):85-94

Reviews of geographic software in this article: DEMO-GRAPHICS: WORLD POPULATIONS AND PROJECTIONS. ESP GAUSS. CEMODEL S. Damus LIMDEP. William H. Greene MICROSTAT 4.1 OTIS PCIPS. (Personal Computer Image Processing System) . H.J. Meyers and R. Bernstein. REGRESSION ANALYSIS OF TIME SERIES (RATS) SPSS/PC+ URBAN DATA MANAGEMENT SOFTWARE (UDMS)  相似文献   

时空概念中的椭球改正     

夏哲仁 《测绘学报》1991,20(3):190-200

  相似文献   

岩体三维节点偶对分析及其应用          下载免费PDF全文

王思敬  张菊明 《地质科学》1991,(3):277-287

本文讨论了节点偶对分析的三维扩展。节点偶对分析是针对摩擦滑动节理单元的约束边界进行的序列矢量判定。在三维节理单元的分析中,几何约束和相应的力学判据极为复杂,采用节点偶对等效结合力方法才能迎刃而解。文中对三维转向节理模型的分析得到合理的结果,它表明本文所提出的方法可成功地应用于节理岩体的有限元分析。  相似文献