Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

A novel coupled FDM-DEM modelling method for flexible membrane boundary in laboratory tests

Difficulties are involved in discrete element method (DEM) modelling of the flexible boundary, that is, the membranes covering the soil sample, which can be commonly found in contemporary laboratory soil tests. In this paper, a novel method is proposed wherein the finite difference method (FDM) and DEM are coupled to simulate the rubber membrane and soil body, respectively. Numerical plane strain and triaxial tests, served by the flexible membrane, are implemented and analysed later. The effect of the membrane modulus on the measurement accuracy is considered, with analytical formulae derived to judge the significance of this effect. Based on an analysis of stress-strain responses and the grain rotation field, the mechanical performances produced by the flexible and rigid lateral boundaries are compared for the plane strain test. The results show that (1) the effect of the membrane on the test result becomes more significant at larger strain level because the membrane applies additional lateral confining pressure to the soil body; (2) the tested models reproduce typical stress and volumetric paths for specimens with shear bands; (3) for the plane strain test, the rigid lateral boundary derives a much higher peak strength and larger bulk dilatation, but a similar residual strength, compared with the flexible boundary. The latter produces a more uniform (or ‘diffuse') rotation field and more mobilised local kinematics than does the former. All simulations show that the proposed FDM-DEM coupling method is able to simulate laboratory tests with a flexible boundary membrane.     

Revisiting the existence of an effective stress for wet granular soils with micromechanics

A possible effective stress variable for wet granular materials is numerically investigated based on an adapted discrete element method (DEM) model for an ideal three‐phase system. The DEM simulations consider granular materials made of nearly monodisperse spherical particles, in the pendular regime with the pore fluid mixture consisting of distinct water menisci bridging particle pairs. The contact force‐related stress contribution to the total stresses is isolated and tested as the effective stress candidate for dense or loose systems. It is first recalled that this contact stress tensor is indeed an adequate effective stress that describes stress limit states of wet samples with the same Mohr‐Coulomb criterion associated with their dry counterparts. As for constitutive relationships, it is demonstrated that the contact stress tensor used in conjunction with dry constitutive relations does describe the strains of wet samples during an initial strain regime but not beyond. Outside this so‐called quasi‐static strain regime, whose extent is much greater for dense than loose materials, dramatic changes in the contact network prevent macroscale contact stress‐strain relationships to apply in the same manner to dry and unsaturated conditions. The presented numerical results also reveal unexpected constitutive bifurcations for the loose material, related to stick‐slip macrobehavior.     

基于噪声谱比法的内蒙古西部地区台站场地响应研究

利用内蒙古西部12个台站的地脉动噪声数据，采用噪声谱比法研究台站的场地响应情况。分析表明，台站场地响应按曲线形态可分为3类，且可能受地形地貌、局部构造和台基状况等条件影响。对比分析噪声谱比法与Moya方法的场地响应结果发现，场地响应曲线形态基本一致，只有极少数台站存在明显差异。结果表明，内蒙古西部区域大部分台站的台基状况总体较好，场地响应曲线较为平坦，无明显频率放大点。     

合并会计报表的电算编制方法探讨

利用用友软件对如何合并会计报表进行了初步研究，并以案例的形式阐明了采用用友软件实现合并会计报表的步骤，提出了合并会计报表的编制方法。     

Gamma-Hadron Separation using Čerenkov Photon Density Fluctuations

In the atmospheric Čerenkov technique γ-rays are detected against the abundant background produced by hadronic showers. In order to improve the signal to noise ratio of theexperiment, it is necessary to reject a significant fraction of hadronic showers. Traditional background rejection methods based on image shape parameters have been extensively used for the data from imaging telescopes. However, non-imaging Čerenkov telescopes have to develop very different means of statistically identifying and removing cosmic ray events. Some of the parameters, which could be potentially important for non-imaging arrays, are the temporal and spectral differences, the lateral distributions and density fluctuations of Čerenkov photons generated by γ-ray and hadron primaries. Here we study the differences in fluctuations of Čerenkov photon density in the light pool at the observation level from showers initiated by photons and those initiated by protons or heavier nuclei. The database of simulated events for the PACT array has been used to evaluate the efficiency of the new technique. Various types of density fluctuations like the short range and medium range fluctuations as well as flatness parameter are studied. The estimated quality factors reflect the efficiencies with which the hadrons can be rejected from the data. Since some of these parameters are independent, the cuts may be applied in tandem and we demonstrate that the proton rejection efficiency of ∼90% can be achieved. Use of density fluctuations is particularly suited for wavefront sampling observations and it seems to be a good technique to improve the signal to noise ratio. This revised version was published online in July 2006 with corrections to the Cover Date.     

RECENT TECTONIC DEFORMATION ANOMALY AND EARTHQUAKES IN GANSU-NINGXIA-QINGHAI AREA

1　IntroductionAccordingtogeologicalstructure ,theGansu Ningxia Qinghaiareabelongstothenortheastmar ginofQinghai Tibetblock .Thisareahasbeenpaidmuchattentionby geo specialistsinChinaandabroadbecauseofitssignificanttectonicmovement,itsintensiveseismicity ,anditsimportanceinearth quakehazardmitigation .IntheDevelopmentPro gramonNationalKeyBasicResearchesundertheProject“MechanismandPredictionofContinentalStrongEarthquakes”,themechanismsofcontinen talstrongearthquakesarestudied ,usinghypo…     

轴对称稳恒磁场的一种解法

本文依据稳恒磁场中的高斯定理和安培环路定理，用迭代法推导出对称轴上的磁感应强度与空间任意一点磁感应强度的关系式，从而可得出轴对称稳恒磁场的无电流分布空间的解式。     

Influence of the Choice of Core-Envelope Transition Point on the Binary Merger of Two Main-sequence Components

We have studied the influence of different choices of core-envelope transition point on the final merger of contact binaries with two main-sequence components. A binary of 1.00 + 0.90M⊙ with an initial orbital period of 0.35d is examined. The mass fraction of the primary mixed with the matter of the secondary, qmix, determined by the chosen core-envelope transition point, ranges from 0.04 to 1.00 in our analysis. If as qmix< 0.8, none of the     

观驾山—黑岛一带信息量法找金研究

本文以成矿地质模型为依据，对观驾山—黑山岛一带的化探、物探、遥感及地质等信息变量进行了详细的分析，选取了１０个信息变量２４个标志状态；在此基础上，运用信息量法对研究区内４５８个单元的信息量值进行计算，并由此确定了研究范围内的成矿远景区。事实表明，研究区信息量成矿预测是非常有效的