弹性波的三维有限元模拟

本文讨论了弹性波瞬态传播问题的三维有限元计算方法及当前存在的实际困难.针对要求计算机内存大和计算时间长的问题,采取了改进措施.由于采用了集中质量矩阵和修正的中心差分时间积分显格式相结合的方法,可以使计算机内存和计算时间大为减少;由于采用结点定位法,最适合用于目前发展的并行计算机系统,可使计算速度大大增快;还采用了有效激发震源法,有效激发区是随时间步进的增加而逐步增大,这不仅能节省计算时间,而且使波场的传播过程一目了然,本文计算了由两种介质组成的三维楔形问题,得到若干典型剖面的瞬时波场图及随时间变化的合成地震图.     

使用价值分析优化方法在功能—技术矩阵求解中的应用

功能—技术矩阵,经众多工作者的努力研究,已渐趋成熟与实用。但如何利用功能—技术矩阵构造在广义现代设计领域中适用的科学模式,解决多学科设计方案的组成与评价问题,还存在薄弱环节。本文利用“使用价值分析优化”的概念,从设计方案技术先进性和经济合理性的协调统一的角度出发,较好地实现了功能—技术矩阵可行优化组合链的组成与评价,弥补了以往的缺陷。     

20世纪末广东省三大城市大气环境质量综合评价——改进密切值法的扩展应用

大气环境质量综合评价是一个城市环境质量考核的重要指标，改进的密切值法是其中一种切实可行的评价方法。通过对1991-2000年10年间广东省三大城市——广州、深圳和珠海的大气环境质量进行综合评价，发现20世纪末，大气环境质量最好的是珠海市，最差的是广州市。与广东省平均水平相比、三大城市中。广州市污染程度高于全省平均水平。从空气环境质量的变化趋势来看，自1997年开始，三大城市大气环境质量趋于好转，为新世纪城市空气环境质量改善提供了良好的基础。     

Tectonic relation between northeastern China and the Korean peninsula revealed by interpretation of GRACE satellite gravity data

The major continental blocks in northeastern Asia are the North China block and the South China block, which have collided starting from the Korean peninsula. Geologic and geophysical interpretations reveal a well defined suture zone in northeastern China from Qinling through Dabie to Jiaodong. The discovery of high-pressure metamorphic rocks in the Hongseong area of the Korean peninsula, prominent evidence for the collision zone, indicates extension of the collision zone in northeastern China into the Korean peninsula. Interpretation of the GRACE satellite gravity dataset shows two prominent structural boundaries in the Yellow Sea. One extends from the Jiaodong Belt in eastern China to the Imjingang Belt in the Korean peninsula. The other extends from near Nanjing, eastern China, to Hongseong. Tectonic movement in or near the suture zone may be responsible for seismic activity in the western Korean peninsula and the development of the Yellow Sea sedimentary basin.     

Semi-analytical solutions of a contaminant transport equation with nonlinear sorption in 1D

A new method to determine semi-analytical solutions of one-dimensional contaminant transport problem with nonlinear sorption is described. This method is based on operator splitting approach where the convective transport is solved exactly and the diffusive transport by finite volume method. The exact solutions for all sorption isotherms of Freundlich and Langmuir type are presented for the case of piecewise constant initial profile and zero diffusion. Very precise numerical results for transport with small diffusion can be obtained even for larger time steps (e.g., when the Courant-Friedrichs-Lewy (CFL) condition failed).     

Measurement of the total nitrogen in lake sediments： Comparison and its environmental significance

Nitrogen cycle is an important bio-geochemical process in the environment. Measurement of the total nitrogen （TN） is a routine experiment in agriculture, biology and environmental sciences. The Kjeldahl method （KM） and elemental analyzer method （EA） are both commonly used to determine TN. Total nitrogen by EA is the sum of nitrate （NO3）, nitrite （NO2）, organic nitrogen and ammonia. Total nitrogen by KM （TKN） is made up of both organic nitrogen and ammonia. A comparative study focused on the two methods is conducted by analysis of TN in 97 samples from the sediment sequence of Gouchi, a salt lake in North China. KM presents a higher degree of accuracy than EA with a standard deviation of 0.007 vs. 0.024. With the presence of nitrate and/or nitrite nitrogen, however, measurement by KM is considerably lower than that by EA. Therefore, for samples from lake sediment sequences or soils in North China, KM is inapplicable to determining TN because of usually high contents of nitrous salt. Despite the inconsistency, use of both methods to the same samples makes sense in reconstructions of climatic and environmental changes from lake sediments. In Lake Gouchi, the contents of nitrate and nitrite nitrogen vary from 1.40% in the lower part of the sequence to 14.77% in the uppermost part, suggesting a gradual evolution process from a fresh water lake to the present-day salt lake.     

Engineering geological studies and tunnel support design at Sulakyurt dam site, Turkey

This paper presents the results of engineering geological investigations and tunnel support design studies, carried out at the Sulakyurt dam site, northeast of Ankara, Turkey. The Sulakyurt dam will be used for flow control and water storage for irrigation projects. Studies were carried out both in the field and the laboratory. Field studies include engineering geological mapping, intensive discontinuity surveying, core drilling and sampling for laboratory testing. The diversion tunnel will be driven in rock mass, consisting of granite and diorite. Empirical, analytical and numerical methods were combined for safe tunnel design. Rock mass rating (RMR), Rock mass quality (Q) and Geological strength index (GSI) systems were used for empirical rock mass quality determination, site characterization and support design. The convergence–confinement method was used as analytical method and software called Phase², a 2D finite element program, was utilized as numerical method. According to the results acquired from the empirical, analytical and numerical methods, tunnel stability problems were expected in both granite and diorite rock masses. The support system, suggested by empirical methods, was applied and the performance of suggested support system was evaluated by means of numerical modelling. It was concluded that the suggested support systems were adequate, since after applying the suggested support system to granite and diorite, tunnel deformation and the yielded elements around the tunnel decreased significantly. Thus, it is suggested that for more reliable support design empirical, numerical and analytical methods should be combined.     

Modeling the seismic response of basin-edge structure by a sequential waveform method

A sequential waveform method is developed to simulate the seismic response of basin-edge structure excited by a plane incident P-wave. The full procedure involves: (a) a previous parameterization of the investigated model using the seismic wave velocities and depths of the sedimentary stratifications; (b) an input motion determined from the records at stations installed on hard rock; (c) forward computation of the P-SV elastic wave field by means of a two-dimensional finite difference (FD) method; (d) the optimization of the model vector using simulated annealing technique and comparing the simulated seismic response of the tested structure with the observed wave field; (e) the correction of the initial model by trial-and-error by testing the differences between synthetics and observed data, and (f) the final solution obtained by iteration using the conjugate gradient algorithm. The search of an optimal basin-edge model has been parallel processed by varying the shapes and velocities of strata on the basis of the fitting of relative timing, amplitude and phase between the output and the observed data. The input motion and sensitivity have been checked and the validity of the method has been demonstrated by numeric analysis. Using the teleseismic records generated by 7 earthquakes recorded at 26 broadband seismic stations, we have studied the seismic velocity structure of the southern edge of the Jiyang depression located in the Bohai Bay basin, northern China. Two cross sections show an agreement between the velocity results and the geological sections available in the region. In addition, we obtain evidence of three hidden faults under the sections and features that suggest major extensions at the Paleogene.     

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.