Equivalent continuum strain calculations based on 3D particle kinematic measurements of sand

Constitutive modeling of granular materials has been a subject of extensive research for many years. While the calculation of the Cauchy stress tensor using the discrete element method has been well established in the literature, the formulation and interpretation of the strain tensor are not as well documented. According to Bagi, 1 researchers mostly adopt well‐known continuum or discrete microstructural approaches to calculate strains within granular materials. However, neither of the 2 approaches can fully capture the behavior of granular materials. They are considered complementary to each other where each has its own strengths and limitations in solving granular‐mechanics problems. Zhang and Regueiro 2 proposed an equivalent continuum approach to calculating finite strain measures at the local level in granular materials subjected to large deformations. They used three‐dimensional discrete element method results to compare the proposed strains measures. This paper presents an experimental application of the Zhang and Regueiro 2 approach using three‐dimensional synchrotron microcomputed tomography images of a sheared Ottawa sand specimen. Invariant Eulerian finite strain measures were calculated for representative element volumes within the specimen. The spatial maps of Eulerian octahedral shear and volumetric strain were used to identify zones of intense shearing within the specimen and compared well with maps of incremental particle translation and rotation for the same specimen. The local Eulerian volumetric strain was compared to the global volumetric strains, which also can be considered as an averaging of all local Eulerian volumetric strains.     

Numerical modeling of three‐phase dissolution of underground cavities using a diffuse interface model

Natural evaporite dissolution in the subsurface can lead to cavities having critical dimensions in the sense of mechanical stability. Geomechanical effects may be significant for people and infrastructures because the underground dissolution may lead to subsidence or collapse (sinkholes). The knowledge of the cavity evolution in space and time is thus crucial in many cases. In this paper, we describe the use of a local nonequilibrium diffuse interface model for solving dissolution problems involving multimoving interfaces within three phases, that is, solid–liquid–gas as found in superficial aquifers and karsts. This paper generalizes developments achieved in the fluid–solid case, that is, the saturated case [1]. On one hand, a local nonequilibrium dissolution porous medium theory allows to describe the solid–liquid interface as a diffuse layer characterized by the evolution of a phase indicator (e.g., porosity). On the other hand, the liquid–gas interface evolution is computed using a classical porous medium two‐phase flow model involving a phase saturation, that is, generalized Darcy's laws. Such a diffuse interface model formulation is suitable for the implementation of a finite element or finite volume numerical model on a fixed grid without an explicit treatment of the interface movement. A numerical model has been implemented using a finite volume formulation with adaptive meshing (e.g., adaptive mesh refinement), which improves significantly the computational efficiency and accuracy because fine gridding may be attached to the dissolution front. Finally, some examples of three‐phase dissolution problems including density effects are also provided to illustrate the interest of the proposed theoretical and numerical framework. Copyright © 2014 John Wiley & Sons, Ltd.     

Modelling the time‐dependent behaviour of granular material with hypoplasticity

This paper presents a constitutive model for time‐dependent behaviour of granular material. The model consists of 2 parts representing the inviscid and viscous behaviour of granular materials. The inviscid part is a rate‐independent hypoplastic constitutive model. The viscous part is represented by a rheological model, which contains a high‐order term denoting the strain acceleration. The proposed model is validated by simulating some element tests on granular soils. Our model is able to model not only the non‐isotach behaviour but also the 3 creep stages, namely, primary, secondary, and tertiary creep, in a unified way.     

对欧拉矢量附加不等式约束的RELSM模型及其应用

借鉴选权拟合思想,对区域地壳运动的整体旋转和线性应变模型中欧拉矢量部分附加不等式约束,将区域地壳运动的欧拉参数约束到由多年观测资料确定的平均欧拉参数值附近,建立附有不等式约束的整体旋转和线性应变模型。对中国大陆构造环境监测网络在环渤海区域的近几期GPS观测速度场数据进行拟合分析,结果表明,本文方法通过不等式约束可适当增强块体整体运动趋势的平滑性,获得对速度场更高的拟合精度,可进一步增强模型的物理解释能力。     

不同采样率下四分量钻孔应变观测同震变化特征

利用湟源台四分量钻孔应变观测的分钟值、1 sps、10 sps和100 sps四种不同采样率的观测数据,通过自检分析、同震应变阶分析和频谱分析等方法研究青海玛多M7.4地震同震变化特征。研究结果表明,分钟采样记录的地震波信息缺失严重,用分钟采样数据进行地震波初动、同震变化幅度等研究将会得到信度较低的结论;采样率越高,记录应变地震波信息的能力越强,但100 sps采样和10 sps采样结果相差不大,10 sps采样已能记录到比较全面的应变地震波信息;同震应变阶的变化性质和变化幅度与采样率无关;未来布设四分量钻孔应变仪时,建议应将采样率至少提高至1 sps。     

1998年1月10日张北-尚义6.2级地震前昌平地震台应变资料的异常变化分析

分析了昌平地震台体应变、线应变及差应变等观测资料在１９９８年１月１０日张北－尚义６．２级地震前的异常变化。结果表明：该次地震前昌平台的应变观测资料出现了明显的趋势异常;体应变的异常时间与震级有较发的相关性;差应变资料显示了短临异常;由线应变资料计算的主应变方向在地震前逐渐向震中偏转;Ｍ２波潮汐因子也出现了明显的异常变化。     

A Study on the Distribution of Wave Phases

Probability distributions of wave phases in association with distributions of surface elevations arestudied with wave records.Wave records of different nature are used for comparison.These are surface fluc-tuations acquired during wind wave flume experiments,representing wave generation under strong wind:andwave records measured in the northern part of Taiwan for waves in natural environments.Three probabilitymodels,the unifrom distribution,the beta distribution,and a model from Tayfun and Lo(1989)are adoptedto study the possible distributions of wave phases.It is found that when surface elevations become skewed,wave phases deviate from the usually assumed uniform distribution and a better model would be the beta dis-tribution.     

易县地震台四分量钻孔应变观测震前异常特征

通过对河北易县地震台钻孔应变资料进行处理,发现在2014年9月6日涿鹿4.3级和2015年9月14日昌黎4.2级地震前,该台四分量钻孔应变观测到显著前兆异常变化,与2012年5月28日唐山4.8级地震前出现的异常相似度较高,异常形态均表现为"来快去慢",均在极值附近维持数日后,逐渐恢复原值过程中发生地震。     

仪器标定格值变化对钻孔应变观测数据的影响

以巴南地震台钻孔应变某一时段观测资料异常变化为例,指出观测仪器元件标定格值是影响钻孔应变观测数据的重要参数,与固体潮汐幅度成反比例关系,且格值变化也可能是改变钻孔应变趋势形态的一种因素。文中分析格值与应变观测数据之间的物理机制,校正恢复异常时段观测曲线,为钻孔应变异常干扰识别工作提供参考。     

青藏高原东南缘活动断裂的地震应变能释放系统研究

本文以青藏高原东南缘为研究区域,利用G-R震级能量经验公式和Benioff地震应变能释放曲线,对该区域内1500年以来的历史地震应变能释放进行了系统性的研究。文中给出了各断裂带和断块区的地震应变能释放周期表,及相应的地震危险性。分析发现研究区域地震应变能的释放具有东强西弱,南强北弱的特征,整体上各断层断块区的历史地震应变能释放符合准周期模式,某些断层和断块区上的地震周期具有某种程度上的同步现象。青藏高原东南缘现今处于大释放期中,地震的危险性不能忽视。局部结果显示,安宁河-则木河断裂带、小江断裂带的危险性很高,对于这些危险区要重点跟踪研究。今后仍需结合不同研究方法来提高地震危险性评估的可靠性。