基于三维裂隙网络的裂隙岩体表征单元体研究

岩体表征单元体（REV）的研究一直是岩土工程中的重要问题。裂隙岩体REV是其内部复杂结构在尺寸效应上的体现,归因于岩体内裂隙系统的随机性。在三维网络模型的基础上通过计算体积节理数建立了确定岩体REV的方法。这种方法考虑了岩体内裂隙的三维空间密度的特征,采用体积节理数这一参数的收敛状况确定岩体REV的大小。建立了实际岩体的模拟裂隙系统,用一定大小的立方体来分割裂隙系统,统计不同立方体内的体积节理数,并将均值作为最终裂隙岩体的参数。将此参数作为裂隙岩体的属性,通过判别其收敛状况获得了岩体REV的尺寸。考虑了岩体不连续面的特征,建立了REV的分析模型,得出工程岩体的REV的尺寸是裂隙平均迹长的4倍的结论。     

石灰改良膨胀土的应力-应变-强度特征与本构描述

为探讨石灰改良膨胀土的变形特征与破坏机制,以压实度为95 %的荆门石灰改良膨胀土为研究对象,开展了单轴、侧限、三轴压缩应力状态下的力学性质试验。试验结果表明：即使湿化饱和后石灰土也具有较高的刚度与强度;单轴压缩状态下,无论饱和还是非饱和状态,石灰土的破坏都为典型的脆性破坏;三轴压缩状态下石灰土破坏前剪缩,同时伴随应变强化,即将破坏时开始剪胀,随后表现为应变软化;围压对石灰土的脆性破坏与破坏后的剪胀有一定的抑制作用,但即使在200 kPa围压下,试样仍发生脆性破坏。在辨明石灰土应力、变形机制的基础上,选用Duncan模型描述其脆性破坏前表现出的压硬性、剪缩性与应变强化特性,标定了相应模型参数,通过数值模拟与平行试验的对比验证了模型的适用性与参数的可靠性。     

利用旁压试验的砂土初始状态反演分析

讨论了砂土旁压试验的反分析问题。建议利用旁压试验实测数据反演确定土体的初始状态而非材料参数,材料参数可以通过实验室常规试验确定。主要利用了新型的亚塑性本构模型参数不受土体状态影响的特点。主要内容包括：建立有限元数值模型,利用文献中的试验数据确定了特定砂土的材料参数;通过模拟砂土的实验室旁压试验,验证了数值模型;通过改变土体初始应力状态和相对密实度进行数值计算,并利用数值计算成果,建立了可适用于砂土初始状态反演分析的实测旁压荷载和土体中应力与相对密度的数学关系。     

Modeling the effects of historical vegetation change on near-surface atmosphere in the northern Chihuahuan Desert

Our goal was to evaluate effects of broad-scale changes in vegetation from grasslands to shrublands over the past 150 years on near-surface atmosphere over the Jornada Experimental Range in the northern Chihuahuan Desert, using a regional climate model. Simulations were conducted using 1858 and 1998 vegetation maps, and data collected in the field. Overall, the vegetation shift led to small changes in sensible heat (SH) and an increase in latent heat (LH). The impacts of shrub encroachment depended on shrubland type: conversion from grass to mesquite cools the near-surface atmosphere and from grass to creosotebush warms it. Higher albedo of mesquite relative to grasses reduced available energy, which was dissipated mainly as LH due to the deeper root system in mesquite. In creosotebush-dominated areas, a decrease in albedo, an increase in roughness length and displacement height contributed to the SH increase and warmer temperatures. Sensitivity simulations showed that an increase in soil moisture content enhanced shrub LH and a reduction in mesquite cover enhanced the temperature differences. The observed shift in vegetation led to complex interactions between land and surface fluxes, demonstrating that vegetation itself is a weather and climate variable as it significantly influences temperature and humidity.     

Simulating the 1999 Capbreton canyon turbidity current with a Cellular Automata model

A numerical model has been developed for the simulation of turbidity currents driven by nonuniform, non cohesive sediment and flowing over a complex three dimensional submarine topography. The model is based on an alternative approach known as Cellular Automata paradigm. The model is validated by comparing a simulation with a reported field-scale event. The chosen case is a turbidity current which occurred in Capbreton Canyon and was initiated by a storm in December 1999. Using data from recent oceanographic cruises, the deposit of the event has been precisely described, which constrain values of model parameters. The model simulates the 1999 turbidity current over the actual canyon topography and related turbidite using three different types of particle. The model successfully simulates areas of erosion and deposition in the canyon. It predicts the vertical and longitudinal grain size evolution, and shows that the fining-up sequence can be deposited by several phases of deposition and erosion related to the current energetic variation during its evolution. This result could explain the presence of intrabed contacts or the frequent lack of facies in Bouma sequences.     

基于紧凑存储的任意四边形网格有限元法地震波场数值模拟

采用有限元法精确模拟复杂介质条件下的地震波场,一般采用三角单元,但在节点数相同的情况下,三角单元的计算精度不如矩形单元高,采用三角单元模拟复杂界面时,编制确定结构刚度矩阵非零元素位置的程序也较为麻烦。采用矩形单元离散含有倾斜或起伏界面的地质模型时,无法避免绕射噪声,若加密网格又会增加计算量。为此,本文基于任意四边形单元模拟声波的传播,在倾斜或起伏界面条件下,可以有效避免因离散引起的“阶梯状”界面,在不增加计算量以及内存占用的前提下,有效地消除离散绕射噪声。采用对角的集中质量矩阵代替一致质量矩阵,避免矩阵的求逆运算,从而提高显式有限元法的计算效率;对结构刚度矩阵采用紧凑存储格式,每一行需要存储的元素最多为5个,同时零元素不参与运算,既减少内存的占用,又极大地提高计算效率。      

砂岩侵入体系模拟及形成机理分析

砂岩侵入体系是一种普遍的地质现象,近些年来越来越受到地学界的重视,但对其形成机理探讨较少。作者及研究团队通过设计一套室内模拟实验装置,来简化模拟3层地层结构中砂岩侵入体系的形成和演化过程。在之前模拟实验(实验变量为顶层沉积物的厚度、进水管的结构)结果的基础上,进一步考察了地形坡度对实验结果的影响,并对砂岩侵入体系的形成机理进行了分析。实验过程中,共观察到7种不同形态不同样式的管道,并对实验过程中压力变化进行了测试。对实验结果分析认为:盖层厚度越大,压力下降的速度越慢,越不利于压力的释放;随着实验的继续进行,砂岩侵入体最终会转化为砂岩喷出体并喷出地表;不稳定的地形有利于触发疏松沉积物发生变形。实验表明,超压是砂岩侵入体发生的最重要条件之一, 沉积盆地中的压实不均衡和生烃作用是可独立产生大规模超压的2种主要机制。     

Comparison of 1D linear,equivalent-linear,and nonlinear site response models at six KiK-net validation sites

Vertical seismometer arrays represent a unique interaction between observed and predicted ground motions, and they are especially helpful for validating and comparing site response models. In this study, we perform comprehensive linear, equivalent-linear, and nonlinear site response analyses of 191 ground motions recorded at six validation sites in the Kiban–Kyoshin network (KiK-net) of vertical seismometer arrays in Japan. These sites, which span a range of geologic conditions, are selected because they meet the basic assumptions of one-dimensional (1D) wave propagation, and are therefore ideal for validating and calibrating 1D nonlinear soil models. We employ the equivalent-linear site response program SHAKE, the nonlinear site response program DEEPSOIL, and a nonlinear site response overlay model within the general finite element program Abaqus/Explicit. Using the results from this broad range of ground motions, we quantify the uncertainties of the alternative site response models, measure the strain levels at which the models break down, and provide general recommendations for performing site response analyses. Specifically, we find that at peak shear strains from 0.01% to 0.1%, linear site response models fail to accurately predict short-period ground motions; equivalent-linear and nonlinear models offer a significant improvement at strains beyond this level, with nonlinear models exhibiting a slight improvement over equivalent-linear models at strains greater than approximately 0.05%.     

Newmark sliding block model for pile-reinforced slopes under earthquake loading

Recent studies have demonstrated that the use of a discretely-spaced row of piles can be effective in reducing the deformations of slopes in earthquakes. In this paper, an approximate strain-dependant Newmark sliding-block procedure for pile-reinforced slopes has been developed, for use in analysis and design of the piling scheme, and the model is validated against centrifuge test data. The interaction of the pile within the slipping soil was idealised using a non-linear elasto-plastic (P–y) model, while the interaction within the underlying stable soil was modelled using an elastic response model in which (degraded) soil stiffness is selected for an appropriate amount of shear strain. This combined soil–pile interaction model was incorporated into the improved Newmark methodology for unreinforced slopes presented by Al-defae et al. [1], so that the final method additionally incorporates strain-dependent geometric hardening (slope re-grading). When combined with the strain-dependent pile resistance, the method is therefore applicable to analysis of both the mainshock and subsequent aftershocks acting on the deformed slope. It was observed that the single pile resistance is mobilised rapidly at the start of a strong earthquake and that this and the permanent slope deformation are therefore strongly influenced by pile stiffness properties, pile spacing and the depth of the slip surface. The model shows good agreement with the centrifuge test data in terms of the prediction of permanent deformation at the crest of the slope (important in design for selecting an appropriate pile layout/spacing i.e. S/B) and in terms of the maximum permanent bending moments induced in the piles (important for appropriate structural detailing of the piles), so long as the slip surface depth can be accurately predicted. A method for doing this, based on limit analysis, is also presented and validated.