On the use of discrete element modelling to simulate the viscoelastic deformation behaviour of an idealized asphalt mixture

The use of discrete-element modelling (DEM) to simulate the behaviour of a highly idealized bituminous mixture under uniaxial and triaxial compressive creep tests is investigated in this paper. The idealized mixture comprises single-sized spherical particles (sand) mixed with bitumen and was chosen so that the packing characteristics are known and the behaviour of the mixture is dominated by the bitumen. The bitumen is represented as shear and normal (tensile and compressive) contact stiffnesses. Numerical sample preparation procedures for specimens containing spherical particles or clumps have been developed to ensure that the final specimen is isotropic and has the correct volumetric proportions. An elastic contact was used for the compressive normal contact stiffness and a viscoelastic contact was used for shear and tensile normal contact stiffness. Simulation results show that the idealized mixture is found to dilate when the ratio of compressive to tensile contact stiffness increases as a function of loading time. Uniaxial and triaxial viscoelastic simulations have been performed to investigate the effect of stress ratio on dilation and the numerical results have been verified with experimental data. The effects of introducing a proportion of frictional contacts and a more complex particle shape (clump) on dilation have been examined.     

基于结构台阵记录的地震响应分析

为了研究土-结构相互作用对高层建筑地震响应数值模拟的影响．选择一栋具有地震观测记录的高层建筑（结构台阵）作为研究对象．利用Abaqus有限元软件和2009年姚安地震加速度记录,并假定了水平成层和有倾角成层的2种土介质模型．进行土-结构动力反应分析。结果表明：有倾角成层土-桩-结构相互作用模型的计算结果与观测资料相比有明显的放大现象,在NS方向体现的更加明显。这说明,就高层建筑的抗震分析而言,土介质模型对计算结果影响很大,在进行地震作用下土-结构地震响应分析时．土层模型应尽可能与真实土层相符。否则其计算结果会产生较大误差。     

剪切带断层泥特征及其在胶东金矿床研究中的应用

剪切带断层泥是中低温热液蚀变的产物。研究表明,断层泥的形成演化与剪切带水-力行为间存在强烈的作用与反馈。断层泥的粒度分布、矿物组成以及孔渗结构受断裂带构造变形和流体作用的控制。同时,由于具有组成和结构上的特性,断层泥可以通过增强岩石各向异性和储存或释放流体而改变断裂带剪切强度和流体压力。以胶东金矿集中区为例,讨论了剪切带构造变形、流体输运-反应与断层泥形成演化的耦合作用,及其对剪切带流变学行为及金矿化形成和分布的影响。     

南宁膨胀土非线性流变模型研究

对南宁非饱和膨胀土进行一系列的加载-卸载固结蠕变试验,得到了不同含水率的膨胀土在不同应力水平下应变-时间曲线和不同时刻的应力-应变等时曲线簇。通过对膨胀土各阶段蠕变变形的分析,得到了膨胀土的固结蠕变是包含有黏弹性变形和黏塑性变形成分的非线性蠕变。采用模型流变与经验流变理论相结合的方法,将膨胀土的蠕变变形分成线性和非线性两部分进行分析,蠕变柔量和线性黏弹性模量随时间减小,而线性黏塑性模量随时间增大。结合膨胀土的线性黏塑性模型和非线性黏塑性经验模型,建立能够描述膨胀土非线性流变的黏塑性本构方程。根据试验所得数据,通过非线性拟合方法得到各模型参数,用该模型进行数值分析得到的理论与蠕变试验结果对比非常吻合。研究成果为实际工程提供了流变变形计算分析的可靠依据。     

盐岩储气库运营期失效概率分析研究

随着能源地下储备的日益迫切需求,越来越多的盐岩地下储库开始投入使用。盐岩地下储库因其良好的使用功能,愈发被世界各国所采用。因此,对盐岩储气库运营期的安全稳定性研究具有重要的理论意义和工程应用价值。运营期间,储库腔体将受到随时间变化的储气内压的作用,同时,客观存在的诸多不确定性因素不可避免地使盐岩地下储库具有一定的风险性,而由于这种不确定性很难采用确定性方法对其进行研究。以金坛盐岩地下储气库为例,基于有限元软件ANSYS,建立储气库运营期腔体内压-时间历程曲线,对储气库长期流变进行了瞬态力学分析,并基于随机力学分析方法,以体积收缩作为风险控制指标,利用ANSYS-PDS得出储气库在整个运营期失效概率的变化规律。计算表明,在周期性变化的储气内压作用下,腔体位移呈波浪线形不断增长,但增幅逐渐降低;同时,储气库失效概率呈现先慢速增长再快速增长最后逐渐趋于稳定的变化规律。探讨了注采气速率对储气库稳定性的影响,适当降低注采速率可抑制腔体收缩,有利于提高储气库的稳定性。对类似盐岩储库运营期风险分析与评估具有借鉴意义。     

岩体流变情况下隧道合理支护时机的数值模拟

以广梧高速公路牛车顶隧道为工程背景,探讨岩体流变情况下隧道二衬支护时机确定的理论和数值计算方法。对比理论方法预测拱顶下沉时程曲线与现场实测结果,建立数值计算模型,分析隧道围岩应力释放系数对围岩稳定性的影响表明,理论方法预测拱顶下沉时程曲线与现场实测表现出相同的发展规律,并且二者得到的二衬支护时间基本相同,验证了理论方法的正确性。围岩应力随着掌子面推进而不断释放,III级围岩应力释放比较缓慢,空间效应比较明显;当开挖释放应力系数大于或等于60%时,整个围岩的松动区都大于锚喷加固范围,最终确定出二衬的合理支护时机,为工程实践提供参考。     

颗粒劣化效应的堆石料流变细观数值模拟

堆石料流变的主要机理是由于水位变化、降雨入渗、日晒雨淋等环境因素导致堆石料发生性质劣化,与此同时,颗粒发生高接触应力-破碎和重新排列-应力释放、调整和转移,这一过程由于颗粒的持续劣化而不断重复。在考虑颗粒破碎的堆石体不连续变形分析方法(SGDD)中,引入颗粒强度劣化模型反映颗粒强度随外界环境的持续劣化。应用该方法进行堆石料三轴流变数值试验,模拟结果与室内试验所观察到的规律一致,表明考虑流变效应的SGDD方法抓住堆石料流变的主要机理,适合模拟堆石料的流变变形这一复杂的、非线性演化问题。数值试验结果表明,堆石料随外界环境的劣化程度、劣化速率、母岩强度对宏观流变变形有较大影响。     

Analysis of the role of diffraction in topographic site effects using boundary element techniques

The role played by the diffraction field on the problem of seismic site effects is studied. For that purpose we solve and analyze simple scattering problems under P and SV in-plane wave assumptions, using two well known direct boundary-element-based numerical methods. After establishing the difference between scattered and diffracted motions, and introducing the concept of artificious and physically based incoming fields, we obtain the amplitude of the Fourier spectra for the diffracted part of the response: this is achieved after establishing the connection between the spatial distribution of the transfer function over the studied simple topographies and the diffracted field. From the numerical simulations it is observed that this diffracted part of the response is responsible for the amplification of the surface ground motions due to the geometric effect. Furthermore, it is also found that the diffraction field sets in a fingerprint of the topographic effect in the total ground motions. These conclusions are further supported by observations in the time-domain in terms of snapshots of the propagation patterns over the complete computational model. In this sense the geometric singularities are clearly identified as sources of diffraction and for the considered range of dimensionless frequencies it is evident that larger amplifications are obtained for the geometries containing a larger number of diffraction sources thus resulting in a stronger topographic effect. The need for closed-form solutions of canonical problems to construct a robust analysis method based on the diffraction field is identified.     

Seismic coda due to non-linear elasticity

Non-linear elastic response of rocks has been widely observed in laboratory, but very few seismic studies are reported in the literature, even though it is the most natural environment where this feature could be observed. Analytic solutions to the non-linear wave propagation phenomena are not readily available, and there is a need to use approximated techniques. It is clear that when a seismic wave propagates through a homogeneous non-linear elastic media, it will be perturbed by the non-linearity. This perturbation can be treated as a source of scattering, spreading the energy of the primary wave in space and time, contributing to the seismic coda. This is in some sense similar to the effect of heterogeneities. The properties of the coda due to the non-linearity depend on the amount of non-linearity and the seismic moment. Using a perturbation approach we calculate the amplitude of the scattered waves, and show that it can describe reasonably well the main features of real seismic codas.     

分形褶皱复杂性的流变因素研究

根据Ramberg的纵弯褶皱粘性力学实验,在褶皱形态的分形分析基础上,利用分形理论和褶皱的流变学理论导出了褶皱的分数维(D)与岩层厚度(h)和粘度(μ)间的关系式,并探讨了褶皱复杂性对褶皱分数维的影响,从中获得有关复杂褶皱的流变学信息。影响分形褶皱复杂程度的因素很多,主要因素包括岩层的厚度和粘度。因此,对褶皱的分形测量和岩层厚度及粘度的分析,可以定量分析分形褶皱形成的流变机理。这一研究是褶皱的非线性流变学理论研究的一个尝试。