地下水非稳定流的LTFAM算法

渗流域内应用拉普拉斯变换(LT)建立相应的有限分析(FAM)方程,顾及渗流域内地下水流的初始条件和边界条件,可在LT空间构成一个封闭的以水头像函数为变量的线性方程组。将此方程组所得的解,通过Stehfest数值反演公式,可归化为时间域的解(水头)。由于时间t被隐含在数值方程内,从而克服了传统数值法按时段(△t)逐步迭代的缺陷,提高了计算效率,也为用嵌入法建立地下水流管理模型提供了一条捷径。     

“三角形法”在区域国土资源综合评价中的应用*——以湖北大别山区为例

本文建立了国土资源评价的指标系统,应用“三角形法”进行了综合评价,高度概括了区域国土资源的宏观特征。     

Les failles récentes des Grands-Moulins (Savoie) : un Sackung (tassement gravitaire) majeur dans les Alpes françaises

In Savoy, the Grands-Moulins recent fault scarps, previously interpreted as seismic fault ruptures, are in fact part of a major Sackung (deep seated gravitational spreading) of the French Alps (9 km long). We mapped more than 60 sackung scarps, some of them reaching 1330 m long and 30 m high. These antislope scarps stop the active screes and offset relict Dryassic rock glaciers by 16 m. We present geomorphologic observations attesting for their gravitational origin. This Sackung is primarily due to glacial debuttressing, while seismic shaking could be a triggering mechanism. To cite this article: J.-C. Hippolyte et al., C. R. Geoscience 338 (2006).     

Analysis of a rigid footing lying on three-layered soil using the finite difference method

A numerical procedure is described for the analysis of the vertical deformation and the stress distribution of the strip footings on layered soil media. Three layers of soil with different stiffness are considered with the middle soil layer the thinnest and most stiff layer. The soil media is discretized and using the theory of elasticity, the governing differential equations are obtained in terms of vertical and horizontal displacements. These equations along with appropriate boundary and continuity conditions are solved by using the finite difference method. The vertical and horizontal displacements, strains and stresses are found at various nodes in the soil media. Parametric studies are carried out to study the effect of the placement depth of the middle soil layer, the relative ratios of the moduli of deformation of the soil layers on the vertical displacement of the footing and the vertical stress distribution. These studies reveal that the middle thin but very stiff layer acts like a plate and redistributes the stresses on the lower soft soil layer uniformly. The displacement on the top and bottom of the middle soil layer is almost the same showing that the compression of the middle layer is negligible as it is very stiff.     

Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.     

A traveltime reciprocity discrepancy in the Podvin & Lecomte time3d finite difference algorithm

We have found that the extensively used finite difference scheme time3d produces time fields which are 'asymmetric' in the sense that traveltimes computed to the right of the source are faster than traveltimes computed to the left. All finite difference schemes will, as they are approximations to the wave equation, to some extent fail to obey reciprocity perfectly. We show, however, that the errors in time3d may be significant—and unnecessarily large. An asymmetry in the point source initialization has been identified, and after correction time3d produces time fields with an improved reciprocity.     

基于有限差分法的抗滑桩计算机辅助设计

基于地基系数“m-m”法、“m-k”法、“k-k”法的原理,考虑桩顶和桩底边界条件以及桩在滑动面处位移、转角、弯矩和剪力的连续条件,可解得桩身各节点的位移和内力,提出了进行抗滑桩全桩内力计算的有限差分法。根据差分方程并用VB6.0编制了实用的计算程序,既可避免繁琐的查表计算,提高计算速度,又可提高计算精度,直观生动,真正实现了人机交互,在界面的引导下,设计人员可完成全部计算,并绘出内力图形和抗滑桩截面配筋图,使设计更方便快捷,该软件可以极大地提高生产效率,降低工程造价,从而实现抗滑桩的优化设计。最后采用上述方法对某滑坡的悬臂抗滑桩进行了设计与计算。     

51463迁站后与原站址的气候要素差异分析

迁站造成气象要素的不连续性，对资料使用带来困难，利用迁站后的对比观测资料，分析气温、风、降水资源的差异发现，季平均气温、最高、最低气温老站址均高于新站址；季平均风速、降水量老站址小于新站性，同时提出了一些差异量，供资料使用者参考。