Micromechanical parameters in bonded particle method for modelling of brittle material failure   总被引：1，自引：0，他引：1  

T. Kazerani  J. Zhao 《国际地质力学数值与分析法杂志》2010,34(18):1877-1895

Bonded particle modelling (BPM) is nowadays being extensively used for simulating brittle material failure. In BPM, material is modelled as a dense assemblage of particles (grains) connected together by contacts (cement). This sort of modelling seriously depends on the mechanical properties of particle and contact, which are named here as micro‐parameters. However, a definite calibration methodology to obtain micro‐parameters has not been so far established; and many have reported some serious problems. In this research, a calibration procedure to find a unique set of micro‐parameters is established. To attain this purpose, discrete element code of UDEC is used to perform BPM. This code can be conveniently developed by the user. The proposed BPM is composed of rigid polygonal particles interacting at their contact points. These contacts can undergo a certain amount of tension, and their shear resistance is provided by cohesion and friction angle. The results demonstrate that each material macro‐property (i.e. Young's modulus, Poisson's ratio, internal friction angel, internal cohesion, and tensile strength) is directly originated from and distinctly related to the contact properties (i.e. normal and shear stiffness, friction angel, cohesion, and tensile strength). Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

The integration of thermal infrared imaging,discharge measurements and numerical simulation to quantify the relative contributions of freshwater inflows to small estuaries in Atlantic Canada     

Serban Danielescu  Kerry T. B. MacQuarrie  Russell N. Faux 《水文研究》2009,23(20):2847-2859

Nutrient fluxes from developed catchments are often a significant factor in the declining water quality and ecological functioning in estuaries. Determining the relative contributions of surface water and groundwater discharge to nutrient‐sensitive estuaries is required because these two pathways may be characterized by different nutrient concentrations and temporal variability, and may thus require different remedial actions. Quantifying the volumetric discharge of groundwater, which may occur via diffuse seepage or springs, remains a significant challenge. In this contribution, the total discharge of freshwater, including groundwater, to two small nutrient‐sensitive estuaries in Prince Edward Island (Canada) is assessed using a unique combination of airborne thermal infrared imaging, direct discharge measurements in streams and shoreline springs, and numerical simulation of groundwater flow. The results of the thermal infrared surveys indicate that groundwater discharge occurs at discrete locations (springs) along the shoreline of both estuaries, which can be attributed to the fractured sandstone bedrock aquifer. The discharge measured at a sub‐set of the springs correlates well with the area of the thermal signal attributed to each discharge location and this information was used to determine the total spring discharge to each estuary. Stream discharge is shown to be the largest volumetric contribution of freshwater to both estuaries (83% for Trout River estuary and 78% for McIntyre Creek estuary); however, groundwater discharge is significant at between 13% and 18% of the total discharge. Comparison of the results from catchment‐scale groundwater flow models and the analysis of spring discharge suggest that diffuse seepage to both estuaries comprises only about 25% of the total groundwater discharge. The methods employed in this research provide a useful framework for determining the relative volumetric contributions of surface water and groundwater to small estuaries and the findings are expected to be relevant to other fractured sandstone coastal catchments in Atlantic Canada. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

Application of the Takagi–Sugeno fuzzy system for combination forecasting of river flow in semiarid mountain regions     

Lijie Zhang  Wenzhi Zhao  Zhibin He  Hu Liu 《水文研究》2009,23(10):1430-1436

Because of high spatial heterogeneity and the degree of uncertainty about hydrological processes in large‐scale catchments of semiarid mountain areas, satisfactory forecasting of daily discharge is seldom available using a single model in many practical cases. In this paper the Takagi–Sugeno fuzzy system (TS) and the simple average method (SAM) are applied to combine forecasts of three individual models, namely, the simple linear model (SLM), the seasonally based linear perturbation model (LPM) and the nearest neighbour linear perturbation model (NNLPM) for modelling daily discharge, and the performance of modelling results is compared in five catchments of semiarid areas. It is found that the TS combination model gives good predictions. The results confirm that better prediction accuracy can be obtained by combining the forecasts of different models with the Takagi–Sugeno fuzzy system in semi‐arid mountain areas. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

Metamorphic evolution,partial melting and rapid exhumation above an ancient flat slab: insights from the San Emigdio Schist,southern California     

A. D. CHAPMAN  P. I. LUFFI  J. B. SALEEBY  S. PETERSEN 《Journal of Metamorphic Geology》2011,29(6):601-626

The San Emigdio and related Pelona, Orocopia, Rand and Sierra de Salinas schists of southern California were underplated beneath the southern Sierra Nevada batholith and adjacent southern California batholith along a shallow segment of the subducting Farallon plate in Late Cretaceous to early Tertiary time. These subduction accretion assemblages represent a regional, deeply exhumed, shallowly dipping domain from an ancient slab segmentation system and record the complete life cycle of the segmentation process from initial flattening and compression to final extensional collapse. An important unresolved question regarding shallow subduction zones concerns how the thermal structure evolves during the slab flattening process. New field relationships, thermobarometry, thermodynamic modelling and garnet diffusion modelling are presented that speak to this issue and elucidate the tectonics of underplating and exhumation of the San Emigdio Schist. We document an upsection increase in peak temperature (i.e. inverted metamorphism), from 590 to 700 °C, peak pressures ranging from 8.5 to 11.1 kbar, limited partial melting, microstructural evidence for large seismic events, rapid cooling (825–380 °C Myr^?1) from peak conditions and an ‘out and back’P–T path. While inverted metamorphism is a characteristic feature of southern California schists, the presence of partial melt and high temperatures (>650 °C) are restricted to exposures with maximum depositional ages between 80 and 90 Ma. Progressive cooling and tectonic underplating beneath an initially hot upper plate following the onset of shallow subduction provide a working hypothesis explaining high temperatures and partial melting in San Emigdio and Sierra de Salinas schists, inverted metamorphism in the schist as a whole, and the observed P–T trajectory calculated from the San Emigdio body. Lower temperatures in Pelona, Orocopia and Rand schists are likewise explained in the context of this overarching model. These results are consistent with an inferred tectonic evolution from shallow subduction beneath the then recently active Late Cretaceous arc to exhumation by rapid trench‐directed channelized extrusion in the subducted schist.  相似文献   

Tectono‐metamorphic history recorded in garnet porphyroblasts: insights from thermodynamic modelling and electron backscatter diffraction analysis of inclusion trails     

E. SKRZYPEK  K. SCHULMANN  P. ŠTÍPSKÁ  F. CHOPIN  J. LEHMANN  O. LEXA  J. HALODA 《Journal of Metamorphic Geology》2011,29(4):473-496

In a Barrovian metamorphic sequence, garnetiferous mica schists document a heterogeneously developed superposition of sub‐orthogonal fabrics and multiple garnet growth episodes. In the variably deformed domains, four types of garnet porphyroblasts have been defined based on inclusion trail patterns. Modelled garnet zoning in the MnNCKFMASHTO system indicates a prograde evolution from 4–4.5 kbar and 490–510 °C to 5–6 kbar and 520–550 °C in the earliest subhorizontal fabric progressing towards 6.5–7.5 kbar and 560–590 °C in the subsequent subvertical foliation. This fabric is heterogeneously deformed into a shallow‐dipping retrograde foliation associated with garnet resorption. In situ electron backscatter diffraction measurements of ilmenite inclusions in individual garnet grains yield precise data on included planar and linear elements. Consistent orientations of internal foliations, lineations and foliation intersection axis sets indicate a superposition of three sub‐orthogonal foliation systems. Weak variations of internal records with increasing intensity of deformation suggest that a moderate buckling stage occurred, but apparent lack of porphyroblast rotation is interpreted as a result of dominant passive flow. Coupling the orientation of internal fabric sets with P–T estimates is used to complement the tectono‐metamorphic evolution of the thickened crust. We demonstrate that garnet porphyroblasts preserve features which reflect large‐scale tectonic processes in orogens.  相似文献   

Modelling coastal ground‐ and surface‐water interactions using an integrated approach     

Dekui Yuan  Binliang Lin 《水文研究》2009,23(19):2804-2817

Beach water table fluctuations have an impact on the transport of beach sediments and the exchange of solute and mass between coastal aquifer and nearby water bodies. Details are given of the refinement of a dynamically integrated ground‐ and surface‐water model, and its application to study ground‐ and surface‐water interactions in coastal regions. The depth‐integrated shallow‐water equations are used to represent the surface‐water flow, and the extended Darcy's equation is used to represent the groundwater flow, with a hydrostatic pressure distribution being assumed to apply for both these two types of flows. At the intertidal region, the model has two layers, with the surface‐water layer being located on the top of the groundwater layer. The governing equations for these two types of flows are discretized in a similar manner and they are combined to give one set of linear algebraic equations that can be solved efficiently. The model is used to predict water level distributions across sloping beaches, where the water table in the aquifer may or may not decouple from the free water surface. Five cases are used to test the model for simulating beach water table fluctuations induced by tides, with the model predictions being compared with existing analytical solutions and laboratory and field data published in the literature. The numerical model results show that the integrated model is capable of simulating the combined ground‐ and surface‐water flows in coastal areas. Detailed analysis is undertaken to investigate the capability of the model. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

A Linking Test that investigates the feasibility of inverse modelling: application to a simple rainfall interception model for Mt Gambier,southeast South Australia     

J. A. P. Pollacco  R. Angulo‐Jaramilo 《水文研究》2009,23(14):2023-2032

Interception loss has an important influence on the water yield of forested areas. Nevertheless, in most studies stemflow is not measured, therefore the question of how to determine the feasibility of optimizing interception and stemflow parameters simultaneously by matching daily simulated throughfall to fortnightly measurements of cumulative throughfall is an important one. By applying a daily empirical interception model, a goodness fit of 2·2 mm/day is obtained between observed and simulated cumulative throughfall. However, by applying the simple but robust Linking Test, it was shown that the parameters are non‐unique and falsely linked, i.e. inter‐relationships between different vegetation parameter sets give similar throughfall but non‐unique net precipitation. The Linking Test investigates the causes of obtaining falsely linked parameters and shows that objective equifinality is not the source of the problem. Objective equifinality occurs when an inappropriate objective function is used. The Linking Test also shows that falsely linked parameters are not caused by measuring throughfall on a non‐daily basis (termed frequency sampling equifinality). By expanding the interception model to the second degree, it was found that the non‐uniqueness is due to the inherent nature of interception and stemflow functions that behave similarly and therefore can easily compensate each other (termed similarity equifinality). It is also shown that a simple daily empirical exponential interception model developed for conifers in the uplands of the United Kingdom is suitable to model interception in Pinus radiata plantations in the Mediterranean climate of southern Australia by using only daily gross precipitation data and two parameters. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

Experimental and numerical behaviour of reinforced mortar plates subjected to accelerated corrosion     

Quang Thanh Nguyen  Sabine Caré  Yves Berthaud  Alain Millard  Frédéric Ragueneau 《国际地质力学数值与分析法杂志》2011,35(10):1141-1159

Most reinforced concrete structures are damaged due to corrosion of reinforcements in concrete. In normal conditions the pH near the reinforcements is around 12–13 which means that steel is in a passive state. But aggressive species, such as chloride ions or carbon dioxide, may penetrate into concrete and promote active corrosion. As a consequence (hydro)oxides are produced leading to degradation of concrete structures. For instance cracking of the concrete is generated due to the pressure induced by rust. In this paper, we study the inception and the propagation of cracking on reinforced mortar plates with rebars located either in the middle or at the corner. Additional experiments have been performed on cylindrical specimens to determine the local effect of rust pressure at the interface rust/mortar. The specimens have been subjected to imposed current density in order to enhance the corrosion and digital image intercorrelation has been used to determine displacement fields. The experiments have been compared to numerical modelling. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

Polymetamorphism in the mainland Lewisian complex,NW Scotland – phase equilibria and geochronological constraints from the Cnoc an t’Sidhean suite     

T. E. JOHNSON  R. W. WHITE  T. ZACK 《Journal of Metamorphic Geology》2012,30(8):865-885

The metamorphic evolution of rocks cropping out near Stoer, within the Assynt terrane of the central region of the mainland Lewisian complex of NW Scotland, is investigated using phase equilibria modelling in the NCKFMASHTO and MnNCKFMASHTO model systems. The focus is on the Cnoc an t’Sidhean suite, garnet‐bearing biotite‐rich rocks (brown gneiss) with rare layers of white mica gneiss, which have been interpreted as sedimentary in origin. The results show that these rocks are polymetamorphic and experienced granulite facies peak metamorphism (Badcallian) followed by retrograde fluid‐driven metamorphism (Inverian) under amphibolite facies conditions. The brown gneisses are inferred to have contained an essentially anhydrous granulite facies peak metamorphic assemblage of garnet, quartz, plagioclase and ilmenite (±rutile, K‐feldspar and pyroxene) with biotite, hornblende, muscovite, chlorite and/or epidote as hydrous retrograde minerals. P–T constraints imposed by phase equilibria modelling imply conditions of 13–16 kbar at >900 °C for the Badcallian granulite facies metamorphic peak, consistent with the field evidence for partial melting in most lithologies. The white mica gneiss comprises a muscovite‐dominated matrix containing porphyroblasts of staurolite, corundum, kyanite and rare garnet. Previous studies have suggested that staurolite, corundum, kyanite and muscovite all grew at the granulite facies peak, with partial melting and melt loss producing a highly aluminous residue. However, at the inferred peak P–T conditions, staurolite and muscovite are not predicted to be stable, suggesting they are retrograde phases that grew during amphibolite facies retrograde metamorphism. The large proportion of mica suggests extensive H₂O‐rich fluid‐influx, consistent with the retrograde growth of hornblende, biotite, epidote and chlorite in the brown gneisses. P–T conditions of 5.0–6.5 kbar at 520–550 °C are derived for the Inverian event. In situ dating of zircon from samples of the white mica gneiss yield apparent ages that are difficult to interpret. However, the data are permissive of granulite facies (Badcallian) metamorphism having occurred at c. 2.7–2.8 Ga with subsequent fluid driven (Inverian) retrogression at c. 2.5–2.6 Ga, consistent with previous interpretations.  相似文献   

Polymetamorphic evolution of pelitic schists and evidence for Permian low-pressure metamorphism in the Vepor Unit, West Carpathians     

P. JE&#;ÁBEK  M. JANÁK  S. W. FARYAD  F. FINGER  P. KONE&#;NÝ 《Journal of Metamorphic Geology》2008,26(4):465-485

Phase equilibrium modelling and monazite microprobe dating were used to characterize the polymetamorphic evolution of metapelites from the northern part of the Vepor Unit, West Carpathians. Three generations of garnet and associated metamorphic assemblages found in these rocks correspond to three distinct metamorphic events related to the Variscan orogeny, a Permian phase of crustal extension and the Alpine orogeny. Variscan staurolite‐bearing and Alpine chloritoid‐bearing assemblages record medium‐temperature and medium‐pressure regional metamorphisms reaching 540–570 °C/5–7.5 kbar and 530–550 °C/5–6.5 kbar respectively. The Permian metamorphic assemblage involves garnet, andalusite, sillimanite, biotite, muscovite, plagioclase and corundum and locally forms silica‐undersaturated andalusite‐biotite‐spinel coronas around older staurolite. The transition from andalusite to sillimanite indicates a prograde low‐pressure and medium‐temperature metamorphism characterized by temperature increase from 500 to 650 °C at ～3 kbar. As accessory monazite is abundant in the rocks, an attempt was made to derive its age of formation by means of electron microprobe‐based Th‐U‐Pb chemical dating. Despite the polymetamorphic nature of the metapelites, the monazite yielded uniform Permian ages. Microstructures confirm that monazite was formed in relation to the low‐pressure and medium‐temperature paragenesis, and the weighted average ages obtained for two different samples are 278 ± 5 and 275 ± 12 Ma respectively. The virtual lack of Variscan and Alpine monazite populations points to interesting aspects concerning the growth systematics of monazite in metamorphic rocks.  相似文献