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.     

地震前后垂直形变场动态演化的量化指标

提出一种表达垂直形变场动态演化过程的量化指标——区域应变率、应变集中度。在此基础上，对南北地震带各水准监测区近30年的垂直形变资料进行了实际计算．并结合具体震例进行了对比研究。结果表明：量化指标在一定程度上能够反映地震前后形变场的动态演化过程和地壳运动状态，对地震的中长期预报有一定的积极意义。     

有限元法与伪谱法混合求解弹性波动方程

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

差分观测与伪距观测的定位精度

严格论证了在对空间或地面目标的定位中,利用伪距测量观测模型的几何因子大于用台站间一次差分观测模型的几何因子.此外,台站间完全差分测量观测的定位精度优于不完全测量情形的定位精度.这些结果给人们在实际工作中采取何种定位方式提供了可靠的理论依据.     

拉格朗日插值法在渔业资源评估上的应用

应用拉格朗日插值法拟合鱼类体长与体重之间的函数关系 ,使用计算机VB语言进行编写程序 ,与常用的线性回归法、Ricker法进行比较 ,得出体长与体重的关系式为Pn(x) =∑nk=0yklk(x)=∑nk=0(∏nj=0j≠kx-xjxk-xj)yk 。     

An approach to quantitative analysis of regional economic links

Aimed at promoting regional coalitions and expanding the approach to economic links, this paper puts forward some new concepts such as link intensity and receiving coefficient, expounds the indexes of quantitative analysis of economic links and establishes the quantitative-analysis model of economic links. With help of the model, this paper calculates the values of the link intensities between Su-Xi-Chang (Suzhou-Wuxi-Changzhou) region and Shanghai, and analyzes the regional difference of the economic links.     

原子钟噪声的蒙特卡络模拟方法

讨论了原子钟噪声的蒙特卡络模拟方法．具体过程是：先产生模拟所需的正态分布随机数，从实际钟的阿仑方差估算出噪声电平，然后用１组递点函数产生高质量的模拟钟时差的随机序列．在对１台商品铯钟进行具体模拟后，给出了真实数据与模拟值的阿仑方差之问的相符程度，最后对实际问题进行了分析．     

Travel time based approach for the assessment of vulnerability of karst groundwater: the Transit Time Method

Karst aquifers represent important water resources in many parts of the world. Unfortunately, karst aquifers are characterised by high contamination risks. This paper presents a travel time based method for the estimation of karst groundwater vulnerability. It considers (1) physics-based lateral flow within the uppermost weathered zone (epikarst) in a limestone-dominated region and (2) high velocities of vertical infiltration at discrete infiltration points (e.g. sinkholes) or lines (e.g. dry valleys, faults). Consequently, the Transit Time Method honours the actual flow path within the unsaturated zone of a karst aquifer system. A test site in Northern Jordan was chosen for the demonstration of the assessment technique, i.e. the catchment area of the Qunayyah Spring north of the capital Amman. The results demonstrate that zones of highest vulnerability lie within valleys and nearby main fault zones. It also reveals that regions, categorised as protected areas by other methods due to thick unsaturated zones, contribute to a major degree to the total risk.     

The double solid reactant method: II. An application to the shallow groundwaters of the Porto Plain,Vulcano Island (Italy)

This paper presents an example of application of the double solid reactant method (DSRM) of Accornero and Marini (Environmental Geology, 2007a), an effective way for modeling the fate of several dissolved trace elements during water–rock interaction. The EQ3/6 software package was used for simulating the irreversible water–rock mass transfer accompanying the generation of the groundwaters of the Porto Plain shallow aquifer, starting from a degassed diluted crateric steam condensate. Reaction path modeling was performed in reaction progress mode and under closed-system conditions. The simulations assumed: (1) bulk dissolution (i.e., without any constraint on the kinetics of dissolution/precipitation reactions) of a single solid phase, a leucite-latitic glass, and (2) precipitation of amorphous silica, barite, alunite, jarosite, anhydrite, kaolinite, a solid mixture of smectites, fluorite, a solid mixture of hydroxides, illite-K, a solid mixture of saponites, a solid mixture of trigonal carbonates and a solid mixture of orthorhombic carbonates. Analytical concentrations of major chemical elements and several trace elements (Cr, Mn, Fe, Ni, Cu, Zn, As, Sr and Ba) in groundwaters were satisfactorily reproduced. In addition to these simulations, similar runs for a rhyolite, a latite and a trachyte permitted to calculate major oxide contents for the authigenic paragenesis which are comparable, to a first approximation, with the corresponding data measured for local altered rocks belonging to the silicic, advanced argillic and intermediate argillic alteration facies. The important role played by both the solid mixture of trigonal carbonates as sequestrator of Mn, Zn, Cu and Ni and the solid mixture of orthorhombic carbonates as scavenger of Sr and Ba is emphasized.

A method to estimate the forces exerted by internal solitons on cylindrical piles

Internal soliton is the large amplitude wave existing in the pycnocline, induced by internal tide in the condition of special bottom topography. During its propagation process, the induced disturbance can bring about strong convergence of sea water and sudden strong current (wave-induced-current), which can cause severe threat to the ocean engineering structures, such as oil drilling platform and pipeline. In this paper, Morison’s empirical method, modal separation and regression analyses are introduced to estimate the forces and torques exerted by internal soliton on cylindrical piles. As an example, a limited set of observational data recording a passage of the internal soliton near Dongsha Islands is used to estimate the horizontal velocity and its acceleration in a vertical section for computing the force and torque on a supposed pile, and the estimation results are reasonable. It is shown that, the higher number of the modes retained in the calculation, the better the estimation of velocity profile fits the observational one. A better overall approximation to the real solution can be reached if there are more observational current data acquired in a whole vertical profile.