虚拟样机技术在地质工程机械领域的应用研究

基于对国内地质工程机械的研发现状的考虑，根据虚拟样机技术的特点，讨论了在地质工程机械领域实施虚拟样机技术的必要件，并提出了建立地质工程机械虚拟样机的策略。联合PRO／E、ANSYS及ADAMS建立了气动潜孔锤的虚拟样机模型，对气动潜孔锤虚拟样机进行了仿真，阐明了虚拟样机技术在地质工程机械领域的应用前景。     

Atomistic calculations of structural and elastic properties of serpentine minerals: the case of lizardite

The physical properties of the hydrous phyllosilicate lizardite have been investigated by atomistic simulation using the GULP code based on transferable semi-empirical interatomic potentials. Lizardite behavior was first investigated during structure relaxation at room temperature. The Helmholtz free energy is minimum for an equilibrium structure that is in agreement with experiment. The bulk, shear, and Young modulii for lizardite were calculated along with the Poisson ratio. From the shear and bulk modulii, we also calculated translational and longitudinal acoustic wave velocities that are important quantities for tectonophysics models. As expected, lizardite is stiffer in the a direction parallel to the layers than in the c perpendicular direction; the variation of the unit cell parameters with pressure is in good agreement with experiment. The cohesive energy between two successive layers along c direction was calculated at 0.33 eV (i.e., 0.11 eV per OH bond) in good agreement with recent ab initio calculations. Upon pressure and temperature variations, we evidenced that structural changes are mainly pressure induced; pressure being accommodated by a decrease of the c parameter up to 10 GPa. We also found that the change of slope in the derivative of the c cell parameter with respect to pressure occurring around 2 GPa originates from the bending of the interlayer hydroxyl groups with respect to the layer normal direction.     

含预制软弱带的岩石破裂过程的数值模拟

考虑到岩石脆性破坏过程中介质的不均匀性之特点，对岩石样品中含预制弱介质条带的岩石样品破坏过程进行二维有限元数值模拟，并对弱介质带的破坏贯通过程、新的断层的产生和有关的地震活动进行了研究．数值模型展示了岩石从变形、微观破坏到整体破坏的全过程以及微震活动的时空分布特征．应力、应变和微震活动的时空分布形象地描绘了岩石变形的局部化和时空迁移等现象，这与实际地壳中所观测到的现象是一致的．此外，模拟结果与实验观测结果也是一致或相似的．      

Atomistic simulation of the structure and elastic properties of pyrite (FeS2) as a function of pressure

Interatomic potential parameters have been derived at simulated temperatures of 0 K and 300 K to model pyrite FeS₂. The predicted pyrite structures are within 1% of those determined experimentally, while the calculated bulk modulus is within 7%. The model is also able to simulate the properties of marcasite, even though no data for this phase were included in the fitting procedure. There is almost no difference in results obtained for pyrite using the two potential sets; however, when used to model FeS₂ marcasite, the potential fitted at 0 K performs better. The potentials have also been used to study the high-pressure behaviour of pyrite up to 44 GPa. The calculated equation of state gives good agreement with experiment and shows that the Fe–S bonds shorten more rapidly that the S–S dimer bonds. The behaviour of marcasite at high pressure is found to be similar to that of pyrite.     

从118.75 GHz附近六通道亮温反演大气温度廓线的数值模拟研究

为了能在静止气象卫星上实现微波被动遥感探测大气温度廓线,并保持一定的地面空间分辨率(如视场小于60 km),就需要使用高频微波及大天线.欧洲和美国下一代静止气象卫星上都已考虑采用118.75 GHz附近通道.为了充分了解118.75 GHz附近通道遥感反演温度廓线的能力,为仪器研制及今后资料的解释反演提供必要的基础数据,作者开展了采用118.75 GHz附近六个通道遥感反演大气温度廓线的数值模拟研究.统计反演的数值试验表明,118.75 GHz附近六通道对温度垂直分布有一定的遥感反演能力;温度反演较好的层次对应于权重函数峰值所在的位置.     

Effects of geology and human activity on the dynamics of salt-water intrusion in three coastal aquifers in southern Spain

 Mathematical modelling of salt-water intrusion processes in three aquifers on the southern coast of Spain (Río Verde, Río Vélez and Castell de Ferro) reveals that, although all three systems are subject to the same climate and seasonal over-exploitation, geological and human factors have very different effects on the dynamics of contamination. In the Río Verde aquifer, the most important influence is the high volume of extractions occurring during the dry season; in Río Vélez, the intrusion is strongly controlled by infiltration of water from the river to the aquifer, and, in the Castell de Ferro system, an intensely karstified carbonate massif lying in contact with both the sea water and the detrital aquifer represents the main entrance point for influx of sea water and subsequent washing of the aquifer. We have undertaken a mathematical simulation of various possible measures to counteract intrusion, according to the specific characteristics of the process in each aquifer. These measures include artificial recharge, use of natural recharge from the river as a hydraulic barrier, and the construction of a low-permeability barrier. Received: 5 December 1995 · Accepted: 12 April 1996     

Large-Eddy Simulation of Turbulent Flow Over a Rough Surface

A family of wall models is proposed that exhibits moresatisfactory performance than previousmodels for the large-eddy simulation (LES) of the turbulentboundary layer over a rough surface.The time and horizontally averaged statistics such asmean vertical profiles of windvelocity, Reynolds stress, turbulent intensities, turbulentkinetic energy and alsospectra are compared with wind-tunnel experimental data.The purpose of the present study is to obtain simulatedturbulent flows that are comparable with wind-tunnelmeasurements for use as the wind environment for thenumerical prediction by LES of source dispersion in theneutral atmospheric boundary layer.     

Turbulent Pressure Statistics in the Atmospheric Boundary Layer from Large-Eddy Simulation

We use large-eddy simulation (LES) to study the turbulent pressure field in atmospheric boundary layers with free convection, forced convection, and stable stratification. We use the Poisson equation for pressure to represent the pressure field as the sum of mean-shear, turbulence–turbulence, subfilter-scale, Coriolis, and buoyancy contributions. We isolate these contributions and study them separately. We find that in the energy-containing range in the free-convection case the turbulence–turbulence pressure dominates over the entire boundary layer. That part dominates also up to midlayer in the forced-convection case; above that the mean-shear pressure dominates. In the stable case the mean-shear pressure dominates over the entire boundary layer.We find evidence of an inertial subrange in the pressure spectrum in the free and forced-convection cases; it is dominated by the turbulence–turbulence pressure and has a three-dimensional spectral constant of about 4.0. This agrees well with quasi-Gaussian predictions but is a factor of 2 less than recent results from direct numerical simulations at moderate Reynolds numbers. Measurements of the inertial subrange pressure spectral constant at high Reynolds numbers, which might now be possible, would be most useful.     

Simulation of oasis breeze circulation in the arid region of the Northwestern China

In this paper, the process of oasis-desert circulation (ODC) is simulated by MM5V3.5 model through designing an ideal oasis-desert scheme and assuming that initial atmosphere is at rest (V = 0). The findings showed that the key of forming special oasis boundary structure is the difference of energy and water between oasis and desert. The evaporation of oasis surface consumes heat energy, and the low temperature of oasis causes an oasis breeze circulation (OBC), which drives an ODC with a downdraft over the oasis and an updraft over the desert. Later, the cold, dry and stable boundary over oasis is gradually formed, on the contrary,the atmospheric boundary over desert on the edge of oasis is hot, humid and unstable and its height is about 600 hPa. The updraft over the desert forms a wet ring that acts as a vertical wall weakening the low-level moisture exchange between the oasis and desert. The downdraft of OBC increases the atmospheric stability that reduces the oasis evaporation. The low-level outflow from the oasis (into the desert) prevents the dry, hot air flowing from the desert into the oasis.Thus an oasis self-preservation mechanism may be formed due to OBC. The horizontal area influenced by oasis is twice as oasis area and the vertical range is four times as oasis. The ODC is strong in the daytime and reaches the strongest at 17:00, and the influenced area is the largest at 20:00.