岩土工程数值模拟新方法

论述了岩土工程数值模拟的新进展和近年来新出现的数值模拟方法,重点介绍了无网络伽辽金法（ＥＦＧＭ）和数值流形法（ＭＭ）的原理及其在岩土工程中的应用,指出了今后研究的重点方向。     

水震波与地震面波的对比研究及其应用

利用同一地质单元内的水震波与地震面波资料进行对比研究,估算水位对含水层应力庆变的响应“格值”,选取唐山矿井水震波资料,唐山陆河地震台记录的地震面波资料,估算“格值”是９．１８２×１０＾－１０／ｃｍ,即唐山矿井观测含水层垂直应变量为９．１８２×１０＾－１０／ｃｍ时,井水位变化１ｃｍ。     

Explicit solution of a free boundary problem for a nonlinear absorption model of mixed saturated-unsaturated flow

In wet soils, zones of saturation naturally develop in the vicinity of impermeable strata, surface ponds and subterranean cavities. Hydrology must be then concerned with transient flow through coexisting unsaturated and saturated zones. The models of advancing saturated zones necessarily involve a nonlinear free boundary problem.A closed-form analytic solution is presented for a nonlinear diffusion model under conditions of ponding at the surface. The soil water diffusivity is restricted to the special functional form D(θ) = a/(b − θ)², where θ is the water content field to be determined and a, b are positive constants. The explicit solution depends on a parameter C (determined by the data of the problem), according to two cases: 1 < C < C₁ or C ≥ C₁, where C₁ is a constant which is obtained as the unique solution of an equation. This result complements the study given in P. Broadbridge, Water Resources Research, 1990, 26, 2435–2443, in order to established when the explicit solution is available. The behavior of the bifurcation parameter C₁ as a function of the driving potential is studied with the corresponding limits for small and large values. Moreover, the sorptivity is proven to be continuously differentiable function of the variable C.     

Heat flow in the Basin and Range province and thermal effects of tectonic extension

In regions of tectonic extension, vertical convective transport of heat in the lithosphere is inevitable. The resulting departure of lithosphere temperature and thickness from conduction-model estimates depends upon the mechanical mode of extension and upon how rapidly extension is (and has been) taking place. Present knowledge of these processes is insufficient to provide adequate constraints on thermal models. The high and variable regional heat flow and the intense local heat discharge at volcanic centers in the Basin and Range province of the United States could be accounted for by regional and local variations in extensional strain rate without invoking anomalous conductive heat flow from the asthenosphere. Anomalous surface heat flow typical of the province could be generated by distributed extension at average rates of about 1/2 to 1%/m.y., similar to rates estimated from structural evidence. To account for higher heat flow in subregions like the Battle mountain High, these rates would be increased by a factor of about 3, and locally at active bimodal volcanic centers, by an order of magnitude more.     

NEAR-FAULT DISPLACEMENT AND DEFORMATION OBTAINED FROM ONE-KILOMETER-LONG FAULT-CROSSING BASELINE MEASUREMENTS-A PRELIMINARY EXPERIMENT AT 2 SITES ON THE EASTERN BOUNDARY OF THE SICHUAN-YUNNAN BLOCK

The current and conventional fault-crossing short baseline measurement has a relatively high precision, but its measurement arrays usually fail to or cannot completely span major active fault zones due to the short length of the baselines, which are only tens to 100 meters. GNSS measurement has relatively low resolution on near-fault deformation and hence is not suitable for monitoring those faults with low motion and deformation rates, due to sparse stations and relatively low accuracy of the GNSS observation. We recently built up two experimental sites on the eastern boundary of the active Sichuan-Yunnan block, one crossing the Daqing section of the Zemuhe Fault and the other crossing the Longshu section of the Zhaotong Fault, aiming to test the measurement of near-fault motion and deformation by using fault-crossing arrays of one-kilometer-long baselines. In this paper, from a three-year-long data set we firstly introduce the selection of the sites and the methods of the measurement. We then calculate and analyze the near-field displacement and strain of the two sites by using three hypothetical models, the rigid body, elastic and composed models, proposed by previous researchers. In the rigid body model, we assume that an observed fault is located between two rigid blocks and the observed variances in baseline lengths result from the relative motion of the blocks. In the elastic model, we assume that a fault deforms uniformly within the fault zone over which a baseline array spans, and in the array baselines in different directions may play roles as strainmeters whose observations allow us to calculate three components of near-fault horizontal strain. In the composed model, we assume that both displacement and strain are accumulated within the fault zone that a baseline array spans, and both contribute to the observed variances in baseline lengths. Our results show that, from the rigid body model, variations in horizontal fault-parallel displacement component of the Zemuhe Fault at the Daqing site fluctuate within 3mm without obvious tendencies. The displacement variation in the fault-normal component keeps dropping in 2015 and 2016 with a cumulative decrease of 6mm, reflecting transverse horizontal compression, and it turns to rise slightly(suggesting extension)in 2017. From the elastic model, the variation in horizontal fault-normal strain component of the fault at Daqing shows mainly compression, with an annual variation close to 10^-5, and variations in the other two strain components are at the order of 10^-6. For the Longshu Fault, the rigid-body displacement of the fault varies totally within a few millimeters, but shows a dextral strike-slip tendency that is consistent with the fault motion known from geological investigation, and the observed dextral-slip rate is about 0.7mm/a on average. The fault-parallel strain component of the Longshu Fault is compressional within 2×10^-6, and the fault-normal strain component is mainly extensional. Restricted by the assumption of rigid-body model, we have to ignore homolateral deformation on either side of an observed fault and attribute such deformation to the fault displacement, resulting in an upper limit estimate of the fault displacement. The elastic model emphasizes more the deformation on an observed fault zone and may give us information about localizations of near-fault strain. The results of the two sites from the composed model suggest that it needs caution when using this model due to that big uncertainty would be introduced in solving relevant equations. Level surveying has also been carried out at the meantime at the two sites. The leveling series of the Daqing site fluctuates within 4mm and shows no tendency, meaning little vertical component of fault motion has been observed at this site; while, from the rigid-body model, the fault-normal motion shows transverse-horizontal compression of up to 6mm, indicating that the motion of the Zemuhe Fault at Daqing is dominantly horizontal. The leveling series of the Longshu site shows a variation with amplitude comparable with that observed from the baseline series here, suggesting a minor component of thrust faulting; while the baseline series of the same site do not present tendencies of fault-normal displacement. Since the steep-dip faults at the two sites are dominantly strike-slip in geological time scale, we ignore probable vertical movement temporarily. In addition, lengths of homolateral baselines on either side of the faults change somewhat over time, and this makes us consider the existence of minor faults on either side of the main faults. These probable minor faults may not reach to the surface and have not been identified through geological mapping; they might result in the observed variances in lengths of homolateral baselines, fortunately such variations are small relative to those in fault-crossing baselines. In summary, the fault-crossing measurement using arrays with one-kilometer-long baselines provides us information about near-fault movement and strain, and has a slightly higher resolution relative to current GNSS observation at similar time and space scales, and therefore this geodetic technology will be used until GNSS networks with dense near-fault stations are available in the future.     

包含局地自由对流的Monin—Obukhov相似性函数

本文从湍流产生和发展的本质出发,以湍流动能和湍流混合长为湍流输送的决定因子,用因次分析的方法得到了一套可包括局地自由对流在内的不稳定大气近地层Monin－Obukhov相似性函数。这套相似性函数不仅有较充分的理论依据,而且还有较明确的物理含意。它能在极端情况下自动蜕化出局地自由对流相似性预测形式。     

Performance of mixed and enhanced finite elements for strain localization in hypoplasticity

Displacement and mixed finite element formulations of shear localization in materials are presented. The formulations are based on hypoplastic constitutive laws for soils and the mixed enhanced treatment involving displacement, strain and stress rates as independently varied fields. Included in these formulations are the standard displacement method, the three‐field mixed formulation, the enhanced assumed strain method and the mixed enhanced strain method. Several numerical examples demonstrating the capability and performance of the different finite element formulations are presented. The numerical results are compared with available experimental data for Hostun RF sand and numerical results for Karlsruhe sand on biaxial tests. Copyright © 2011 John Wiley & Sons, Ltd.     

Microstructural evidence of low‐strain,time‐transgressive subglacial deformation

Microstructural analyses were used to investigate the formation of a macroscale‐massive till at Knud Strand in Denmark. More than 100 thin sections were examined and microstructures mapped and counted for quantitative comparison and interpretation. Microstructures indicative of both brittle (grain lineations, edge‐to‐edge crushed grains) and ductile (turbate structures) deformation are evenly distributed in vertical profiles through the till, suggesting that strain contributed to its formation. Discrete shears (grain lineations and plasmic fabric) probably accommodated most deformation, whereas rotational deformation was less prominent. The microshear geometry fits the predicted Coulomb–Mohr failure criterion, indicating that till behaves as a plastic material. Strain estimate of ca. 10¹ from micromorphological proxies is two–three orders of magnitude lower than expected if the till was subjected to pervasive deformation. A hybrid of lodgement and time‐transgressive deformation is envisaged as the till‐forming processes. Our data suggest that even abundant evidence of microscale deformation at continuing high levels of strain may only record the latest process of deposition and deformation and therefore not fully reflect the complexity of till genesis. Copyright © 2007 John Wiley & Sons, Ltd.     

Strain fields and mechanical response of a highly to medium fissured bentonite clay

This article presents the developments of an ongoing research aimed at modelling the influence of fissuring on the behaviour of clays. In particular, it recalls the main results of an extensive laboratory investigation on a fissured bentonite clay from the south of Italy and presents the data of a new investigation on the evolution with shearing of the strain fields developing within the clay, resulting from Digital Image Correlation (DIC). Element test results are analysed in the framework of continuum mechanics and linked to the clay fissuring features, once characterised using the Fissuring IDentity (F‐ID) chart. This article compares the bentonite behaviour with that of other fissured clays of different F‐IDs, highlighting the common behavioural features. Thereafter, the soil response at the macro level is related to the DIC‐derived strain fields evolving within the clay with loading. For this purpose, DIC was successfully used to investigate the deformation processes active in the fissured clay and the sources of the localisation phenomena. DIC is shown to provide indications of the extent to which highly to medium fissured clays element test results can be of use to model the clay behaviour according to continuum mechanics. Copyright © 2012 John Wiley & Sons, Ltd.     

大姚、民乐-山丹地震前的中短期前兆

利用川滇及青藏块体北缘地区上世纪80年代至2003年的跨断层流动形变观测资料，借助应变强度比时序曲线，通过丽江Ms7．0、景泰Ms5．9等震例研究，获取与2003年7月大姚Ms6．2、Ms6．1地震和10月民乐-山丹Ms6．1地震孕育有关的中短期前兆异常。分析表明：①大姚、民乐-山丹地震前1年内异常场地数远超过场地总数的1／3，相对近源区异常在发震地区占一半；②大姚Ms6．2地震前半年内，盐源-丽江地块及其邻近区域以及北东向延伸的川滇交界地区异常明显增多，相对集中。2003年7～9月大幅减弱；而鲜水河断裂、楚雄-通海断裂南段有所增强；③民乐-山丹震前十余天至数月内相对近源区及其邻近构造区域呈现异常相对分散-相对集中且局部增强-一临震前减弱的时空演变特征。震后一个月内青藏块体北缘地区仍有近1／3的场地呈现异常(海原一六盘山断裂带尤为突出)。