A novel coupled FDM-DEM modelling method for flexible membrane boundary in laboratory tests     

Hua-Xiang Zhu  Zhen-Yu Yin  Qiang Zhang 《国际地质力学数值与分析法杂志》2020,44(3):389-404

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.  相似文献   

Geostatistical features of streambed vertical hydraulic conductivities in Frenchman Creek Watershed in Western Nebraska     

Olufemi P. Abimbola  Aaron R. Mittelstet  Troy E. Gilmore 《水文研究》2020,34(16):3481-3491

This study evaluated the spatial variability of streambed vertical hydraulic conductivity (K_v) in different stream morphologies in the Frenchman Creek Watershed, Western Nebraska, using different variogram models. Streambed K_v values were determined in situ using permeameter tests at 10 sites in Frenchman, Stinking Water and Spring Creeks during the dry season at baseflow conditions. Measurements were taken both in straight and meandering stream channels during a 5 day period at similar flow conditions. Each test site comprised of at least three transects and each transect comprised of at least three K_v measurements. Linear, Gaussian, exponential and spherical variogram models were used with Kriging gridding method for the 10 sites. As a goodness-of-fit statistic for the variogram models, cross-validation results showed differences in the median absolute deviation and the standard deviation of the cross-validation residuals. Results show that using the geometric means of the 10 sites for gridding performs better than using either all the K_v values from the 93 permeameter tests or 10 K_v values from the middle transects and centre permeameters. Incorporating both the spatial variability and the uncertainty involved in the measurement at a reach segment can yield more accurate grid results that can be useful in calibrating K_v at watershed or sub-watershed scales in distributed hydrological models.  相似文献   

Evolutionary Processes in Low-Mass Binary Systems     

Peter P. Eggleton  Ludmila Kisseleva-Eggleton 《Astrophysics and Space Science》2005,296(1-4):327-336

We consider the evolution of certain low-mass binaries, incorporating models of (a) internal evolution, (b) tidal friction, (c) dynamo activity driven by an elementary α,Ω dynamo, (d) stellar wind driven by the activity, and (e) magnetic braking as a consequence of wind and poloidal dynamo-generated magnetic field. In some circumstances the stellar wind is found to remove mass on a nuclear timescale, as is necessary to explain some observed systems. We can hope that various uncertainties in the model may be clarified by a careful comparison of the models with such observed quantities as rotation periods. These are modified by processes (a), (b) and (e). Assuming that stellar evolution is slow, rotation rate should in some circumstances represent a balance between magnetic braking trying to slow the star down and tidal friction trying to spin it up. Preliminary attempts are promising, but indicate that some fine tuning is necessary. When there is a third body present, in an orbit which is inclined but not necessarily of short period, the eccentricity of a close binary can be strongly modified by ‘Kozai cycles’. We show that this may complicate attempts to account for spin rates of stars in close binaries.  相似文献   

Smeared crack approach: back to the original track     

M. Cervera  M. Chiumenti 《国际地质力学数值与分析法杂志》2006,30(12):1173-1199

This paper briefly reviews the formulations used over the last 40 years for the solution of problems involving tensile cracking, with both the discrete and the smeared crack approaches. The paper focuses on the smeared approach, identifying as its main drawbacks the observed mesh‐size and mesh‐bias spurious dependence when the method is applied ‘straightly’. A simple isotropic local damage constitutive model is considered, and the (exponential) softening modulus is regularized according to the material fracture energy and the element size. The continuum and discrete mechanical problems corresponding to both the weak discontinuity (smeared cracks) and the strong discontinuity (discrete cracks) approaches are analysed and the question of propagation of the strain localization band (crack) is identified as the main difficulty to be overcome in the numerical procedure. A tracking technique is used to ensure stability of the solution, attaining the necessary convergence properties of the corresponding discrete finite element formulation. Numerical examples show that the formulation derived is stable and remarkably robust. As a consequence, the results obtained do not suffer from spurious mesh‐size or mesh‐bias dependence, comparing very favourably with those obtained with other fracture and continuum mechanics approaches. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

液核地球自转速率变化的带谐潮效应Ⅱ--数值计算与比较     

张捍卫  郑勇  杜兰  盘关松 《天文研究与技术》2003,(2):67-74

地球自转速率的潮汐变化可由无量纲参数k/cm)(k和cm分别为壳幔的有效Love数和有效极转动惯量)来表征。对于一个具有弹性地幔、平衡海潮和核幔不耦合的地球k/cm=0.944,且与潮波频率无关。海潮的非平衡扰动使k/cm为复数,且与频率有关。大气对自转速率有效勒夫数的贡献约为Δkat=0.0075。同时地幔滞弹性对勒夫数也产生扰动。利用本文得到的理论公式和最新的潮汐数据计算了地球自转速率的潮汐变化,及其有关地球物理机制的影响。  相似文献   

Anthropogenic trace metals in an ice core at Vestfonna,Svalbard, Norway     

Sumito Matoba  Hideaki Motoyama  Hideki Narita  Okitsugu Watanabe 《极地研究(英文版)》2003,14(1)

1　IntroductionAnicecoreobtainedfrom polarglaciersoricesheetsisoneofthemostimportantarchivestoreconstructpaleoclimaticandpaleoatmosphericcondition .Informationonpale o environmentcanbeextractedfromicecoresaschemicaland/orphysicalsignals.Amongthechemicalsignals,heavymetalsarenotedassignalsofterrestrialenvironmentalchangeandanthropogenicpollution (e.g .Murozumietal.1 969;NgandPatterson 1 981 ;Hongetal.1 994) .SinceconcentrationsofmostofthemetalsinpolarsnowincentralGreen landareatorbelowthepptl…  相似文献   

地下洞库勘察中的深斜孔钻探和压水试验技术          下载免费PDF全文

吉孟瑞  孙桂明  吴昌庆  原华山 《山东国土资源》2001,17(2):56-60

青岛市黄岛区拟建设海底地下油气储藏库,前期地质勘察工作要求施工倾角45°、深220m的取心钻孔4个,并在孔内分段进行压水试验.本文介绍了如何利用XY-4型金刚石小口径钻机施工基岩深斜孔及压水试验技术.  相似文献   

Metamorphic evolution of exhumed middle to lower crustal rocks in the Mojave Extensional Belt, southern California, USA     

D. J. HENRY  R. K. DOKKA 《Journal of Metamorphic Geology》1992,10(3):347-364

The Waterman Metamorphic Complex of the central Mojave Desert was exposed as a consequence of early Miocene detachment-dominated extension. However, it has evidence consistent with a more extensive geological history that involves collision of a crustal fragment(s), tectonic thickening by overthrusting and two periods of extension. The metamorphic complex contains granitoid intrusives and felsic mylonitic gneisses as well as polymetamorphic rocks that include marble, calc-silicate, quartzite. mafic granulite, pyribolite, amphibolite, migmatite and biotite schist. The latter group of rocks was affected by an initial series of high-grade metamorphic events (M1 and M2) and a localized lower grade overprint (M3). The initial metamorphism (M1) can be separated into two stages along its high-grade P–T path: M1a, a granulite facies metamorphism at 800–850° C and 7.5–9 kbar and Mlb, an upper amphibolite facies overprint at 750–800° C and 10–12 kbar. M1a developed mineral assemblages and textures consistent with granulite facies conditions at a reduced activity of H₂O and is associated with intense ductile deformation (D1) and minor local partial melting. M1b overprinted the granulite assemblages with a series of hydrous phases under conditions of increasing pressure and H₂O activity and is accompanied by little or no deformation. M2 developed at lower pressures and temperatures (650–750° C, 4.5–5.5 kbar) and is distinguished by a second local overprint of hydrous phases that reflects an input of aqueous fluids probably associated with the intrusion of a series of granitic dykes and veins. Effects of M3 are confined to the Mitchel detachment zone, an anastomosing early Miocene detachment fault, and are characterized by local ductile/brittle deformation (D2) of the pre-existing high-grade rocks and granitoid intrusives and by the production of mylonites and mylonitic gneisses under greenschist facies conditions (300–350° C, 3–5 kbar). The initial overprint (M1a) represents metamorphism, devolatilization and minor partial melting of supracrustal rocks under granulite facies conditions as a consequence of tectonic and, possibly, magmatic thickening. The increasing pressure transition of M1a to M1b reflects a period of continued compressional tectonism, thrusting and influx of H₂O, in part, locally related to crystallization of partial melts. The near isothermal decompression between M1b and M2 probably represents a pre-112-Ma extensional episode that may have been the result of a decompressional readjustment of a thickened crust. Following the initial extensional event, the metamorphic complex remained at depths of 10–17 km for at least 90 Ma until it was uplifted following Miocene extension. M3 develops locally in response to this second extensional period resulting from the early Miocene detachment faulting.  相似文献   

Simulation of double cold cores of the 35°N section in the Yellow Sea with a wave-tide-circulation coupled model   总被引：1，自引：0，他引：1  

Xia Changshui  Qiao Fangli  Zhang Mengning  Yang Yongzeng  Yuan Yeli 《中国海洋湖沼学报》2004,22(3):292-298

Based on the MASNUM wave-tide-circulation coupled numerical model, the temperature structure along 35°N in the Yellow Sea was simulated and compared with the observations. One of the notable features of the temperature structure along 35°N section is the double cold cores phenomena during spring and summer. The double cold cores refer to the two cold water centers located near 122°E and 125°E from the depth of 30m to bottom. The formation, maintenance and disappearance of the double cold cores are discussed. At least two reasons make the temperature in the center (near 123°E) of the section higher than that near the west and east shores in winter. One reason is that the water there is deeper than the west and east sides so its heat content is higher. The other is invasion of the warm water brought by the Yellow Sea Warm Current (YSWC) during winter. This temperature pattern of the lower layer (from 30m to bottom) is maintained through spring and summer when the upper layer (0 to 30m) is heated and strong thermocline is formed. Large zonal span of the 35°N section (about 600 km) makes the cold cores have more opportunity to survive. The double cold cores phenomena disappears in early autumn when the west cold core vanishes first with the dropping of the thermocline position. Supported by the National Basic Research Program of China (No. G1999043809) and the National Science Foundation of China (No. 49736190).  相似文献   

变形分析的Bayes估计与检验     

寇新建 《大地测量与地球动力学》1994,(3)

由于观测条件的变化，观测值方差将产生随机波动。本文根据Ｂａｙｅｓ估计理论，推导了线性模型待估参数的分布函数及观测值方差估计公式，并将其用于位移检验。计算结果表明，在同样的置信度下，本方法更符合实际。文中还对参数估值的统计性质进行了讨论。  相似文献