Finite element modeling of slab–foundation interaction on rigid pavement applications

Bending stresses produced in pavement concrete slabs are influenced by the contact condition along the slab–foundation interface. A finite element model for the analysis of frictional contact between unbonded pavement layers is presented. A 2D plate element is proposed for the idealization of composite slabs. Interface elements are used to connect the unbonded layers within the composite slab and to model the interface between the bottom slab layer and the foundation layer. Interface elements have the ability to capture the separation and sliding between pavement layers, due to thermal loads, and to calculate the frictional traction at their interface.     

Transversely isotropic frost heave modeling with heat–moisture–deformation coupling

The transversely isotropic (TISO) constitutive and frost heave models for the freezing of fine-grained soils are more accurate than the isotropic model and simpler than the orthotropic models. First, in combination with the mesoscopic composition of freezing soils, a mechanical model for the interaction between the equivalent ice lens and the soil in frozen soils is established based on the series and parallel models in the theory of composite mechanics. Second, the TISO constitutive model together with the analytic expression of five elastic constants is provided for analysis of the freezing soils. Third, a preliminary elastoplastic model for TISO freezing soil is established based on the Hill plastic model. Fourth, the heat–moisture–deformation coupling TISO model and the hydrodynamic frost heave model are derived according to a thermodynamics equation, a soil water motion equation, and generalized Hooke’s law. Synchronization and uniformity of the TISO constitutive model and the TISO frost heave model are realized for analyzing the interaction between permafrost soils and buildings. Finally, an indoor standard frost heave test and the frost heave of a prototype canal are simulated based on the above models. The numerical results indicated that the models presented in this paper accurately described the frost heave and revealed the interaction between permafrost and buildings.

     

Application of genetic algorithm technique to inverse modeling of tide–aquifer interaction

Modeling of tide–aquifer interaction plays a vital role in the management of coastal aquifer systems. A novel and robust methodology is presented in this paper for estimating aquifer parameters of coastal aquifers from tide–aquifer interaction data using tide–aquifer interaction model and genetic algorithm (GA). Two stand-alone computer programs were developed to optimize hydraulic diffusivities of unconfined and confined coastal aquifers at multiple sites using GA technique and tide–aquifer interaction model and considering two approaches (‘lumped tidal component approach’ and ‘multi-tidal component approach’). Five sets of real-world tide–aquifer interaction data at two sites of an unconfined aquifer and one set of tide–aquifer interaction data at three sites of a confined aquifer were used to demonstrate the efficacy of the methodology. The analysis of the GA-based inverse modeling results indicated that the ‘multi-tidal component approach’ yields more accurate and reliable hydraulic diffusivities for the unconfined aquifer (RMSE = 0.0129–0.0521 m, NSE = 0.70–0.97, and d1 = 0.91–0.99) as well as for the confined aquifer (RMSE = 0.0204–0.0545 m, NSE = 0.95–0.97, and d1 = 0.99) compared with the ‘lumped tidal component approach’. A comparative evaluation of data-size revealed that the short-duration datasets of the unconfined aquifer provide more reliable estimate of hydraulic diffusivity than the long-duration datasets. Further, it was found that the spring and neap tidal data yield unreasonable values of hydraulic diffusivity with considerably high values of RMSE and very low values of r ², NSE, and d1, thereby suggesting that spring and neap tidal data are not suitable for aquifer parameter estimation. Overall, it is concluded that the GA-based tide–aquifer interaction model following ‘multi-tidal component approach’ is the most efficient tool for estimating aquifer parameters of unconfined and confined aquifers from tide–aquifer interaction data. The developed methodology is also applicable to other coastal basins of the world irrespective of hydrogeological settings.     

Three dimensional Finite Element modeling of seismic soil–structure interaction in soft soil

Earthquakes in regions underlain by soft clay have amply demonstrated the detrimental effects of soil–structure interaction (SSI) in such settings. This paper describes a new three dimensional Finite Element model utilizing linear elastic single degree of freedom (SDOF) structure and a nonlinear elasto-plastic constitutive model for soil behavior in order to capture the nonlinear foundation–soil coupled response under seismic loadings. Results from an experimental SSI centrifuge test were used to verify the reliability of the numerical model followed by parametric studies to evaluate performance of linear elastic structures underlain by soft saturated clay. The results of parametric study demonstrate that rigid slender (tall) structures are highly susceptible to the SSI effects including alteration of natural frequency, foundation rocking and excessive base shear demand. Structure–foundation stiffness and aspect ratios were found to be crucial parameters controlling coupled foundation–structure performance in flexible-base structures. Furthermore, frequency content of input motion, site response and structure must be taken into account to avoid occurrence of resonance problem.     

Swell–consolidation characteristics of artificial sand clay mixes

Many remedial measures have been devised to lessen the damage caused by expansive soils. Physical alteration, chemical stabilization, innovative foundation techniques like belled piers, drilled piers, under-reamed piles and granular pile anchors are some of these remedial measures. Mixing a non swelling material such as gravel or sand to expansive soil is one of the methods of physical alteration. This paper presents experimental data on artificially prepared sand-clay mixes. Swell and consolidation characteristics of these artificially prepared sand-clay mixes were studied in one dimensional consolidometer. Fine sand content and fines content in the expansive soil were arbitrarily varied in the investigation. The fines content was varied as 425–300 μm and 150–75 μm, separated from the same expansive soil based on the grain size. Swell potential and swelling pressure decreased with increasing fine sand content but increased with increasing fines content. Coefficient of compressibility, coefficient of volume compressibility and compression index of the samples decreased initially up to a sand content of 15% and thereafter increased at higher sand contents.     

An analytical model of soil–structure interaction with swelling soils during droughts

Lightly loaded structures constructed on expansive soils may develop structural damage as a result of changes in the soil’s moisture content. This study investigated an analytical model of soil–structure interaction to assess the settlement of dwellings built on swelling soils when droughts occur. The building behavior was investigated with the Euler–Bernoulli beam theory, and the ground behavior was investigated with a Winkler-derived model based on the state surface approach. The analytical model results were compared to those of a finite element analysis using the Barcelona Expansive Model (BExM) performed with Code_Bright.The analytical model was then used to assess the settlement transmission ratio for a typology of clayey soils and different parameters of building. The results indicated that the final deflection of the building increased with the building length and soil suction. The building deflection due to the suction variations was inversely proportional to the load, the rigidity of the building and the embedding depth of the foundation. Increasing these parameters made the building less vulnerable to shrinkage and swelling action.     

Development of an adaptive CTM–RPIM method for modeling large deformation problems in geotechnical engineering

Acta Geotechnica - In this paper, a meshfree method called adaptive CTM–RPIM is developed to model geotechnical problems with large deformation. The developed adaptive CTM–RPIM is a...     

Sediment–well interaction during depressurization

Depressurization gives rise to complex sediment–well interactions that may cause the failure of wells. The situation is aggravated when high depressurization is imposed on sediments subjected to an initially low effective stress, such as in gas production from hydrate accumulations in marine sediments. Sediment–well interaction is examined using a nonlinear finite element simulator. The hydro-mechanically coupled model represents the sediment as a Cam-Clay material, uses a continuous function to capture compressibility from low to high effective stress, and recognizes the dependency of hydraulic conductivity on void ratio. Results highlight the critical effect of hydro-mechanical coupling as compared to constant permeability models: A compact sediment shell develops against the screen, the depressurization zone is significantly smaller than the volume anticipated assuming constant permeability, settlement decreases, and the axial load on the well decreases; in the case of hydrates, gas production will be a small fraction of the mass estimated using a constant permeability model. High compressive axial forces develop in the casing within the production horizon, and the peak force can exceed the yield capacity of the casing and cause its collapse. Also tensile axial forces may develop in the casing above the production horizon as the sediment compacts in the depressurized zone and pulls down from the well. Well engineering should consider: slip joints to accommodate extensional displacement above the production zone, soft telescopic/sliding screen design to minimize the buildup of compressive axial force within the production horizon, and enlarged gravel pack to extend the size of the depressurized zone.     

Mantle melting in equilibrium with an Iron–Wüstite–Graphite buffered COH-fluid

Partial melting experiments on a San Carlos peridotite were done in a Walker type multi-anvil press at pressures from 5 to 12.5 GPa. Experiments were done in the presence of a COH-fluid and at oxygen fugacity controlled by the Fe–FeO buffer. Olivine, clinopyroxene, garnet and orthopyroxene are stable in all but the highest temperature 10 GPa experiments where olivine and garnet coexist, and the highest temperature 5 GPa experiments where olivine is the single crystalline phase. The solidus at 5 GPa was found to be at approximately 1,200°C and the liquidus was estimated to be at 1,325°C, which is ∼500°C lower than has been reported for dry melting of peridotite. The aluminum concentration of the melts decreases with increasing melt fraction and decreases also with increasing pressure. At 5 GPa the melts have a CaO/Al₂O₃-ratio of 0.85–1.0, which is similar to that of undepleted komatiites; major element concentrations are also identical to those of undepleted komatiites such as the Munro komatiites. At 10 and 12.5 GPa the partial melts have CaO/Al₂O₃-ratios above 1.5 and major element composition almost identical to aluminum depleted komatiites such as the Barberton komatiites. We therefore conclude that in the presence of a reducing COH-fluid both aluminum-depleted and -undepleted komatiites could have formed at temperatures much lower than generally accepted.     

Three-dimensional elastic–viscoplastic consolidation behavior of transversely isotropic saturated soils

Acta Geotechnica - This paper investigates three-dimensional elastic–viscoplastic consolidation behaviors of transversely isotropic saturated soils. The Drucker–Prager yield criterion...     

Coupled thermal–hydraulic modeling of artificial ground freezing with uncertainties in pipe inclination and thermal conductivity

Acta Geotechnica - Artificial ground freezing (AGF) has been widely used as a temporary soil stabilization and waterproofing technique in geotechnical practices (e.g., tunnel construction). Many...     

Physicochemical modeling of formation of Ag–Au–Hg solid solutions: Kyuchyus deposit (Yakutia,Russia) as an example

The «Selektor-C» software package and standard thermodynamic functions of ternary Ag–Au–Hg solid solutions were used for developing physicochemical models in natural processes with participation of gold, silver and mercury. On the example of the Kyuchyus Au–Sb–Hg deposit we have worked out hypogene and hypergene models of formation of native gold, including mercuric gold. We obtained thermodynamic evidence that the Kyuchyus deposit ores formed with the origin of electrum at the early main productive quartz–sulfide stage and ternary Au–Ag–Hg solid solutions at the late non-substantive Au-bearing stages.     

Galvanic interaction between galena and pyrite in an open system

1IntroductionMost sulfide minerals in nature have perfect con-ductivity and possess semiconducting properties.Gal-vanic corrosion may occur when two sulfide mineralswith different potentials are connected together in thesolution.The mineral with the highe…     

Noble-metal ore genesis and mantle–crust interaction

The mineral and geochemical compositions of noble-metal (first of all, gold) deposits of the Fennoscandian, Siberian, and Northeast Asian orogenic belts are considered. These deposits are of several types: Au (disseminated Au–sulfide and Au–quartz), Au–Bi, Au–Ag, Au–Sb, Ag–Sb, Au–Sb–Hg, and Ag–Hg. They formed in different geodynamic settings as a result of the active motion of crustal tectonic blocks of different nature. Subduction processes (both at the front and at the rear of continent-marginal and island-arc magmatic arcs) resulted in Au–Ag, Ag–Sb, Ag–Hg, Au–Sb–Hg, and Au–Bi deposits. Collision events gave rise to Au and Au–Bi deposits. Intraplate continental rifting and formation of orogenic belts along the boundaries of block (plate) sliding led to the origin of Au and Au–Bi ores in association with Au–Ag, Au–Sb–Hg, and complex ores. In all cases, the formation of noble-metal mineralization was accompanied by magmatism of different types and metamorphism. Because of this diversity of ores, there is no single concept of the genesis of noble-metal mineralization. Several competing models of genesis exist: hydrothermal-metamorphic, pluton-metamorphic, plutonic, activity of mantle fluid flows, and multistage concentration during the crust–mantle interaction with the leading role of sedimentary complexes.     

Leveraging HPC accelerator architectures with modern techniques — hydrologic modeling on GPUs with ParFlow

Computational Geosciences - Rapidly changing heterogeneous supercomputer architectures pose a great challenge to many scientific communities trying to leverage the latest technology in...     

The North Pole Region: First Data on the Snow–Sea Ice–Ice Water Sedimentation System

Doklady Earth Sciences - The dispersed sedimentary material in the snow–ice cover near the North Pole was analyzed directly for the first time. The composition of sediments in the...     

Geologically controlled agricultural contamination and water–rock interaction in an alluvial aquifer: results from a hydrochemical study

Environmental Earth Sciences - Hydrogeochemistry data collected from three multi-level monitoring wells in a sandy alluvial aquifer located in the Keum River watershed, South Korea, are used in...     

Trace element behaviour during interaction between basalt and crustal rocks at 0.5–0.8 GPa: an experimental approach

We experimentally investigate the major and trace elements behavior during the interaction between two partially molten crustal rocks (meta-anorthosite and metapelite) and a basaltic melt at 0.5–0.8 GPa. Results show that a hybrid melt is formed at the basalt-crust contact, where plagioclase crystallizes. This contact layer is enriched in trace elements which are incompatible with plagioclase crystals. Under these conditions, the trace element diffusion coefficients are one order of magnitude larger than those expected. Moreover, the HFSE diffusivity in the hybrid melt is surprisingly higher than the REE one. Such a feature is related to the plagioclase crystallization that changes the trace elements liquid-liquid partitioning (i.e. diffusivity) over a transient equilibrium that will persist as long as the crystal growth proceeds. These experiments suggests that the behaviour of the trace elements is strongly dependent on the crystallization at the magma-crust interface. Diffusive processes like those investigated can be invoked to explain some unusual chemical features of contaminated magmatic suites.     

Consistent modeling of a catastrophic flowslide at the Shenzhen landfill using a hydro-elasto-plastic model with solid–fluid transition

Flow-type landslides are an important hazard that can cause great destruction due to the rapid flow velocity and large disaster area. This paper presents a catastrophic flowslide that recently occurred at a landfill in Shenzhen, China. This disaster involved an area about 1100 m in length and 630 m in maximum width, and caused the death of 77 people and the destruction of 33 buildings. The precise reason for the landfill’s failure is still unknown, and therefore we try to contribute an increased understanding of the event for future prevention. In this study, the failure mechanism of the studied slope was analyzed and described under partially saturated condition. The solid–fluid transition during the flowslide occurrence was described using a unified constitutive model. The model was used to perform the hydro-elasto-plastic modeling in the pre-failure stage, the viscous modeling in the post-failure stage, and the second-order work criterion was introduced in between to model the solid–fluid transition. The consistent evolution of the flowslide, including initiation, propagation, and deposit stages, was simulated and analyzed using the finite element method with Lagrangian integration points after careful calibration of the viscous parameters. The numerical results were compared with the real case and used to explain the failure mechanism.