首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 281 毫秒
1.
A possible effective stress variable for wet granular materials is numerically investigated based on an adapted discrete element method (DEM) model for an ideal three‐phase system. The DEM simulations consider granular materials made of nearly monodisperse spherical particles, in the pendular regime with the pore fluid mixture consisting of distinct water menisci bridging particle pairs. The contact force‐related stress contribution to the total stresses is isolated and tested as the effective stress candidate for dense or loose systems. It is first recalled that this contact stress tensor is indeed an adequate effective stress that describes stress limit states of wet samples with the same Mohr‐Coulomb criterion associated with their dry counterparts. As for constitutive relationships, it is demonstrated that the contact stress tensor used in conjunction with dry constitutive relations does describe the strains of wet samples during an initial strain regime but not beyond. Outside this so‐called quasi‐static strain regime, whose extent is much greater for dense than loose materials, dramatic changes in the contact network prevent macroscale contact stress‐strain relationships to apply in the same manner to dry and unsaturated conditions. The presented numerical results also reveal unexpected constitutive bifurcations for the loose material, related to stick‐slip macrobehavior.  相似文献   

2.
Computational fluid dynamics and discrete element method (CFD–DEM) is extended with the volume of fluid (VOF) method to model free‐surface flows. The fluid is described on coarse CFD grids by solving locally averaged Navier–Stokes equations, and particles are modelled individually in DEM. Fluid–particle interactions are achieved by exchanging information between DEM and CFD. An advection equation is applied to solve the phase fraction of liquid, in the spirit of VOF, to capture the dynamics of free fluid surface. It also allows inter‐phase volume replacements between the fluid and solid particles. Further, as the size ratio (SR) of fluid cell to particle diameter is limited (i.e. no less than 4) in coarse‐grid CFD–DEM, a porous sphere method is adopted to permit a wider range of particle size without sacrificing the resolution of fluid grids. It makes use of more fluid cells to calculate local porosities. The developed solver (cfdemSolverVOF) is validated in different cases. A dam break case validates the CFD‐component and VOF‐component. Particle sedimentation tests validate the CFD–DEM interaction at various Reynolds numbers. Water‐level rising tests validate the volume exchange among phases. The porous sphere model is validated in both static and dynamic situations. Sensitivity analyses show that the SR can be reduced to 1 using the porous sphere approach, with the accuracy of analyses maintained. This allows more details of the fluid phase to be revealed in the analyses and enhances the applicability of the proposed model to geotechnical problems, where a highly dynamic fluid velocity and a wide range of particle sizes are encountered. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

3.
Analytical solutions for the steady‐state response of an infinite beam resting on a visco‐elastic foundation and subjected to a concentrated load moving with a constant velocity are developed in this paper. The beam responses investigated are deflection, bending moment, shear force and contact pressure. The mechanical resistance of the foundation is modeled using two parameters ks and ts — ks accounts for soil resistance due to compressive strains in the soil and ts accounts for the resistance due to shear strains. Since this model represents the ground behavior more accurately than the Winkler spring model, the developed solutions produce beam responses that are closer to reality than those obtained using the existing solutions for Winkler model. The dynamic beam responses depend on the damping present in the system and on the velocity of the moving load. Based on the study, dynamic amplification curves are developed for beam deflection. Such amplification curves for deflection, bending moment, shear force and contact pressure can be developed for any beam‐foundation system and can be used in design. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

4.
By incorporating the nonlinear variation of a soil's compressibility and permeability during the process of consolidation, an analytical solution for the radial consolidation of vertical drains has been developed for a general time‐variable loading. The general solution was verified for the cases of instantaneous loading and ramp loading. Detailed solutions were further derived for two special loading schemes: multistage loading and preloading–unloading–reloading. The nonlinear consolidation behavior of a vertical drain subjected to these two types of loading schemes was then investigated by a parametric study. The results show that the loading rate, the ratio of the compressibility index to the permeability index (Cc/Ck), and the initial stress state have a significant influence on the consolidation rate. A smaller value of Cc/Ck, a larger initial stress, or a fast loading rate always leads to a rapid consolidation rate. During the unloading period, a negative excess pore water pressure may occur, and a slower unloading rate may reduce this negative value. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

5.
This paper presents results of analysis of full-scale pile load test data of 14 piles embedded in either loose or medium dense sands. The analysis was performed using two methods, p–y curve approach and a more recently developed khmax approach. Comparison of the results obtained using both the methods is also presented. A step-by-step analysis procedure is presented for predicting lateral load deflection response of single piles in sand using the khmax approach. The results presented show that the khmax approach has promise over the p–y curve approach because of its simplicity and the fact that it provides upper- and lower-bound curves, which are valuable guides to making engineering decisions. For loose sands, a new range of khmax values is recommended to better predict the lateral load–deflection response of single piles.  相似文献   

6.
The Jiangzhuang ultrahigh‐pressure (UHP) metamorphic peridotite from south Sulu, eastern China occurs as a layer within gneiss with eclogite blocks, and consists of coarse‐grained garnet porphyroblasts and a fine‐grained matrix assemblage of garnet + forsterite + enstatite + diopside ± phlogopite ± Ti‐clinohumite ± magnesite. Both types of garnet are characterized by high MgO content and depletion of light rare earth element (LREE) and enrichment of heavy rare earth element, but the matrix garnet has lower MgO, TiO2 and higher Cr2O3 and REE contents. Diopside displays LREE enrichment, and has low but variable large‐ion lithophile element (LILE) contents. Phlogopite is a major carrier of LILE. Ti‐clinohumite contains high Nb, Ta, Cr, Ni, V and Co contents. The P–T conditions of 4.5–6.0 GPa and 850–950 °C were estimated for matrix mineral assemblages. Most peridotites are depleted in Al2O3, CaO and TiO2, and enriched in SiO2, K2O, REE and LILE. In contrast to phlogopite‐free peridotites, the phlogopite‐bearing peridotites have higher K2O, Zr, REE and LILE contents. Zircon occurs only in the phlogopite‐bearing peridotites, shows no zoning, with low REE contents and Th/U ratios, and yields tight U–Pb ages of 225–220 Ma, indicating the peridotites experienced consistent Triassic UHP metamorphism with subducted supercrustal rocks. These data demonstrate that the Jiangzhuang peridotites were derived from the depleted mantle wedge of the North China Craton, and experienced various degrees of metasomatism. The phlogopite‐free peridotites may have been subjected to an early cryptic metasomatism at UHP conditions of the mantle wedge, whereas the phlogopite‐bearing peridotites were subjected to a subsequent strong metasomatism, characterized by distinctly enrichment in LILE, LREE, Zr and K as well as the growth of zircon and volatile‐bearing minerals at UHP subduction conditions. The related metasomatism may have resulted from the filtration of fluids sourced mainly from deeply subducted supracrustal rocks.  相似文献   

7.
The vibrations of railway tracks on a poroelastic half‐space generated by moving trains are investigated through a vehicle–track–ground coupling model. The theoretical model incorporates a vehicle, a track, and a fully saturated poroelastic half‐space soil medium. The source of vibration excitation is divided into two components: the quasi‐static loads and the dynamic loads. The quasi‐static loads are related to the static component of the axle loads, whereas the dynamic loads are due to the dynamic wheel–rail interaction. A linear Hertizian contact spring is introduced between each wheelset and the rail to consider the dynamic loads. Biot's dynamic theory is used to characterize the poroelastic half‐space soil medium. Using the Fourier transform, the governing equations for the track–ground system are solved and the numerical results are presented for a single axle vehicle model. The different dynamic characteristics of the elastic soil medium and the saturated poroelastic medium are investigated. In addition, the different roles of the moving axle loads and the roughness‐induced dynamic loads are identified. It is concluded that the vibration level of the free field off the track predicted by the poroelastic soil medium is smaller than that predicted by the elastic soil medium for vehicle speed below the Rayleigh wave speed of the poroelastic half‐space, whereas it is larger for vehicle speed above the Rayleigh wave speed. The dynamic loads play an important role in the dynamic responses of the track–ground system. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

8.
Both oceanic and continental HP rocks are juxtaposed in the Huwan shear zone in the western Dabie orogen, and thus provide a window for testing the buoyancy‐driven exhumation of dense oceanic HP rocks. The HP metamorphic age of the continental rocks in this zone has not been well constrained, and hence it is not known if they are of the same age as the exhumation of the HP oceanic rocks. In situ laser ablation (multiple collector) inductively coupled plasma mass spectrometry (LA‐(MC‐)ICP‐MS), U–Pb, trace element and Hf isotope analyses were made on zircon in a granitic gneiss and two eclogites from the Huwan shear zone. U–Pb age and trace element analysis of residual magmatic zircon in an eclogite constrain its protolith formation at 411 ± 4 Ma. The zircon in this sample displays εHf (t) values of +6.1 to +14.4. The positive εHf (t) values up to +14.4 suggest that the protolith was derived from a relatively depleted mantle source, most likely Palaeotethyan oceanic crust. A granitic gneiss and the other eclogite yield protolith U–Pb ages of 738 ± 6 and 700 ± 14 Ma, respectively, which are both the Neoproterozoic basement rocks of the Yangtze Block. The zircon in the granitic gneiss has low εHf (t) values of ?14.2 to ?10.5 and old TDM2 ages of 2528–2298 Ma, suggesting reworking of Palaeoproterozoic crust during the Neoproterozoic. The zircon in the eclogite has εHf (t) values of ?1.0 to +7.4 and TDM1 ages of 1294–966 Ma, implying prompt reworking of juvenile crust during its protolith formation. Metamorphic zircon in both eclogite samples displays low Th/U ratios, trace element concentrations, relatively flat heavy rare earth element patterns, weak negative Eu anomalies and low 176Lu/177Hf ratios. All these features suggest that the metamorphic zircon formed in the presence of garnet but in the absence of feldspar, and thus under eclogite facies conditions. The metamorphic zircon yields U–Pb ages of 310 ± 3 and 306 ± 7 Ma. Therefore, both the oceanic‐ and continental‐type eclogites share the same episode of Carboniferous eclogite facies metamorphism. This suggests that high‐pressure continental‐type metamorphic rocks might have played a key role in the exhumation and preservation of oceanic‐type eclogites through buoyancy‐driven uplift.  相似文献   

9.
  相似文献   

10.
We present major and trace element data for eighteen 1.71–1.66 Ga granitoid samples, and Sm–Nd whole‐rock isotope data for eleven of these samples, in a transect across the border between the Transscandinavian Igneous Belt (TIB) and the Eastern Segment in central southern Sweden. The geochemistry of the granitoids varies from alkalic to alkali‐calcic and peraluminous in the east to predominantly calc‐alkaline and metaluminous in the west. Rocks in the west also have lower SiO2 contents. Trace element signatures favour formation in an active continental margin setting. Nd isotope data are completely overlapping along the transect and initial εNd values are mildly depleted in the range +0.3 to +2.6. The combined data suggest that the magmas were derived mainly from juvenile, pre‐existing crust, increasingly mafic and less alkaline towards the west. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

11.
Within the framework of compressible multi‐phase flow through deformable porous media, wave–soil interactions in the near‐shore region during wave runup and drawdown are modeled. Critical non‐dimensional parameters governing the interaction processes are identified. Within the context of wave basin and centrifuge wave tank facilities, we propose scaling relations for the experimental investigations of the transient and steady‐state responses of wave–soil systems. Numerical simulations are conducted to illustrate and confirm the theoretical and scaling analyses. Based on the simulations results, the implications on the design of experiments and interpretation of results are discussed. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

12.
  相似文献   

13.
At Bangriposi, variable stages in replacement of staurolite by chloritoid – Na–K–Ca mica shimmer aggregates in muscovite schists provides insight into the complex interplay between fluid flow, mass transfer, and dissolution–precipitation during pseudomorph growth. Idioblastic chloritoid growing into mica caps without causing visible deformation, and monomineralic chloritoid veins (up to 300 μm wide) within shimmer aggregates replacing staurolite attest to chloritoid nucleation in fluid‐filled conduits along staurolite grain boundaries and crystallographic planes. The growth of shimmer aggregates initiated along staurolite margins, and advanced inwards into decomposing staurolite along networks of crystallographically controlled fluid‐filled conduits. Coalescence among alteration zones adjacent to channel fills led to dismemberment and the eventual demise of staurolite. Mass balance calculation within a volume‐fixed, silica‐conserved reference frame indicate the shimmer aggregates grew via precipitation from fluids in response to mass transport that led to the addition of H2O, K2O, Na2O and CaO in the reaction zone, and Al2O3 was transported outward from the inward‐retreating margin of decomposing staurolite. This aided precipitation of chloritoid in veins and in the outer collars, and as disseminated grains in the shimmer aggregates at mid‐crustal condition (~520 ± 20 °C, 5.5 ± 2.0 kbar). Computation using one‐dimensional transport equation suggests that staurolite decomposition involved advection dominating over diffusive transport; the permeation of externally derived H2O caused flattening of chemical potential gradients in H2O and aqueous species, for example, and , computed using the Gibbs method. This suggests that staurolite decomposition was promoted by the infiltration of a large volume of H2O that flattened existing chemical potential gradients. In the initial stages of replacement, chloritoid super‐saturation in fluid caused preferential nucleation and growth of chloritoid at staurolite grain boundaries and in crystallographic planes. As reaction progressed, further chloritoid nucleation was halted, but chloritoid continued to grow as the 3‐mica aggregates continued to replace the remaining staurolite in situ, while the chloritoid‐compatible elements were transported in the water‐rich phase facilitating continued growth of the existing chloritoid grains.  相似文献   

14.
Biot theory predicts wave velocities in a saturated granular medium using the pore geometry, viscosity, densities, and elastic moduli of the solid skeleton and pore fluid, neglecting the interaction between constituent particles and local flow, which becomes essential as the wavelength decreases. Here, a hydro-micromechanical model, for direct numerical simulations of wave propagation in saturated granular media, is implemented by two-way coupling the lattice Boltzmann method (LBM) and the discrete element method (DEM), which resolve the pore-scale hydrodynamics and intergranular behavior, respectively. The coupling scheme is benchmarked with the terminal velocity of a single sphere settling in a fluid. In order to mimic a small amplitude pressure wave entering a saturated granular medium, an oscillating pressure boundary on the fluid is implemented and benchmarked with the one-dimensional wave equation. The effects of input waveforms and frequencies on the dispersion relations in 3D saturated poroelastic media are investigated with granular face-centered-cubic crystals. Finally, the pressure and shear wave velocities predicted by the numerical model at various effective confining pressures are found to be in excellent agreement with Biot analytical solutions, including his prediction for slow compressional waves.  相似文献   

15.
Analytical solutions for wave velocities and wave vectors are yielded for a continuously inhomogeneous cross‐anisotropic medium, in which Young's moduli (E, E′) and shear modulus (G′) varied exponentially as depth increased. However, for the rest moduli in cross‐anisotropic materials, ν and ν′ remained constant regardless of depth. We assume that cross‐anisotropy planes are parallel to the horizontal surface. The generalized Hooke's law, strain–displacement relationships, and equilibrium equations are integrated to constitute governing equations. In these equations, displacement components are fundamental variables and, hence, the solutions of three quasi‐wave velocities, VP, VSV, and VSH, and the wave vectors, $\mathop{\mathop{l}\limits^{\rightharpoonup}}\nolimits_{P}$ $\mathop{\mathop{l}\limits^{\rightharpoonup}}\nolimits_{\mathit{SV}}$, and $\mathop{\mathop{l}\limits^{\rightharpoonup}}\nolimits_{{\mathit{SH}}}$, can be generated for the inhomogeneous cross‐anisotropic media. The proposed solutions and those obtained by Daley and Hron, and Levin correlate well with each other when the inhomogeneity parameter, k, is 0. Additionally, parametric study results indicate that the magnitudes and directions of wave velocity are markedly affected by (1) the inhomogeneous parameter, k; (2) the type and degree of geomaterial anisotropy (E/E′, G′/E′, and ν/ν′); and (3) the phase angle, θ. Consequently, one must consider the influence of inhomogeneous characteristic when investigating the behaviors of wave propagation in a cross‐anisotropic medium. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

16.
In this study, plagiogranites in the Diyanmiao ophiolite of the southeastern Central Asian Orogenic Belt(Altaids) were investigated for the first time. The plagiogranites are composed predominantly of albite and quartz, and occur as irregular intrusive veins in pillow basalts. The plagiogranites have high SiO_2(74.37–76.68 wt%) and low Al_2O_3(11.99–13.30 wt%), and intensively high Na_2O(4.52–5.49 wt%) and low K_2O(0.03–0.40 wt%) resulting in high Na_2O/K_2O ratios(11.3–183). These rocks are classified as part of the low-K tholeiitic series. The plagiogranites have low total rare earth element contents(∑REE)(23.62–39.77 ppm), small negative Eu anomalies(δEu=0.44–0.62), and flat to slightly LREE-depleted chondrite-normalized REE patterns((La/Yb)N=0.68–0.76), similar to N-MORB. The plagiogranites are also characterized by Th, U, Zr, and Hf enrichment, and Nb, P, and Ti depletion, have overall flat primitivemantle-normalized trace element patterns. Field and petrological observations and geochemical data suggest that the plagiogranites in the Diyanmiao ophiolite are similar to fractionation-type plagiogranites. Furthermore, the REE patterns of the plagiogranites are similar to those of the gabbros and pillow basalts in the ophiolite. In plots of SREE–SiO_2, La–SiO_2, and Yb–SiO_2, the plagiogranites, pillow basalts, and gabbros show trends typical of crystal fractionation. As such, the plagiogranites are oceanic in origin, formed by crystal fractionation from basaltic magmas derived from depleted mantle, and are part of the Diyanmiao ophiolite. LA–ICP–MS U–Pb dating of zircons from the plagiogranites yielded ages of 328.6±2.1 and 327.1±2.1 Ma, indicating an early Carboniferous age for the Diyanmiao ophiolite. These results provide petrological and geochronological evidence for the identification of the Erenhot–Hegenshan oceanic basin and Hegenshan suture of the Paleo-Asian Ocean.  相似文献   

17.
Dynamic two‐phase interaction of soil can be modelled by a displacement‐based, two‐phase formulation. The finite element method together with a semi‐implicit Euler–Cromer time‐stepping scheme renders a discrete equation that can be solved by recursion. By experience, it is found that the CFL stability condition for undrained wave propagation is not sufficient for the considered two‐phase formulation to be numerically stable at low values of permeability. Because the stability analysis of the two‐phase formulation is onerous, an analysis is performed on a simplified two‐phase formulation that is derived by assuming an incompressible pore fluid. The deformation of saturated porous media is now captured in a single, second‐order partial differential equation, where the energy dissipation associated with the flow of the fluid relative to the soil skeleton is represented by a damping term. The paper focuses on the different options to discretize the damping term and its effect on the stability criterion. Based on the eigenvalue analyses of a single element, it is observed that in addition to the CFL stability condition, the influence of the permeability must be included. This paper introduces a permeability‐dependent stability criterion. The findings are illustrated and validated with an example for the dynamic response of a sand deposit. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

18.
New activity–composition (a–x) relations for minerals commonly occurring in metapelites are presented for use with the internally consistent thermodynamic dataset of Holland & Powell ( 2011 , Journal of Metamorphic Geology, 29 , 333–383). The a–x relations include a broader consideration of Fe2O3 in minerals, changes to the formalism of several phases and order–disorder in all ferromagnesian minerals where Fe–Mg mixing occurs on multiple sites. The a–x relations for chlorite, biotite, garnet, chloritoid, staurolite, cordierite, orthopyroxene, muscovite, paragonite and margarite have been substantially reparameterized using the approach outlined in the companion paper in this issue. For the first time, the entire set of a–x relations for the common ferromagnesian minerals in metapelitic rocks is parameterized simultaneously, with attention paid to ensuring that they can be used together to calculate phase diagrams of geologically appropriate topology. The a–x relations developed are for use in the Na2O–CaO–K2O–FeO–MgO–Al2O3–SiO2–H2O–TiO2–O2 (NCKFMASHTO) system for both subsolidus and suprasolidus conditions. Petrogenetic grids in KFMASH and KFMASHTO are similar in topology to those produced with earlier end‐member datasets and a–x relations, but with some notable differences. In particular, in subsolidus equilibria, the FeO/(FeO + MgO) of garnet is now greater than in coexisting staurolite, bringing a number of key staurolite‐bearing equilibria into better agreement with inferences from field and petrographic observations. Furthermore, the addition of Fe3+ and Ti to a number of silicate phases allows more plausible equilibria to be calculated in relevant systems. Pseudosections calculated with the new a–x relations are also topologically similar to equivalent diagrams using earlier a–x relations, although with many low variance fields shifting in P–T space to somewhat lower pressure conditions.  相似文献   

19.
Sun  Meng  Cao  Junnan  Cao  Jingjing  Zhang  Shuai  Chen  Yunmin  Bate  Bate 《Acta Geotechnica》2022,17(7):2633-2649

Bioremediation is widely used to improve ground soil by introducing calcium carbonate (CaCO3). Shear wave velocity (Vs) is usually adopted to evaluate effect but the microscopic mechanism is unclear. The discrete element method (DEM), a promising tool for simulating the behaviors of cohesive and noncohesive materials, was used in this study to simulate Vs evolution and wave propagation path of sand reinforced by calcite precipitates. Two basic calcite precipitate forms are proposed for representing individual calcite precipitation (CaCO3-P) and calcite aggregation (CaCO3-C). Contact cementation between adjacent sand grain pairs was the primary association pattern for calcite precipitates at a low calcite content. At a higher calcite content, the preferential shear wave propagation pathway is the clusters cemented by CaCO3-C. With calcite content increasing from 0 to 9%, the coordination number and average contact force increased. Vs increased from 169.73 to 2132.64 m/s but had high variability due to the spatial distribution. The results suggest that the calibrated DEM model can elucidate the microscopic mechanisms and evaluate the enhancement effect of microorganism-reinforced soil.

  相似文献   

20.
Understanding the Cenozoic vertical-axis rotation in the Tibetan Plateau is crucial for continental dynamic evolution. Paleomagnetic and rock magnetic investigations were carried out for the Oligocene and Miocene continental rocks of the Hoh Xil basin in order to better understand the tectonic rotations of central Tibet. The study area was located in the Tongtianhe area located in the southern part of the Hoh Xil basin and northern margin of the Tanggula thrust system in central-northern Tibet. A total of 160 independently oriented paleomagnetic samples were drilled from the Tongtianhe section for this study. The magnetic properties of magnetite and hematite have been recognized by measurements of magnetic susceptibility vs. temperature curves and unblocking temperatures. The mean directions of the Oligocene Yaxicuo Group in stratigraphic coordinates(Declination/Inclination = 354.9°/29.3°, k = 33.0, α_(95) = 13.5°, N =5 Sites) and of the Miocene Wudaoliang Group in stratigraphic coordinates(Declination/Inclination = 3.6°/36.4°, k = 161.0, α_(95) = 9.7°, N =3 Sites) pass reversal tests, indicating the primary nature of the characteristic magnetizations. Our results suggested that the sampled areas in the Tuotuohe depression of the Hoh Xil basin have undergone no paleomagnetically detectable rotations under single thrusting from the Tanggula thrust system. Our findings, together with constraints from other tectonic characteristics reported by previous paleomagnetic studies, suggest tectonic rotations in the Cuoredejia and Wudaoliang depressions of the Hoh Xil basin were affected by strike-slip faulting of the Fenghuo Shan-Nangqian thrust systems. A closer examination of geological data and different vertical-axis rotation magnitudes suggest the tectonic history of the Hoh Xil basin may be controlled by thrust and strike-slip faulting since the Eocene.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号