A steady state model for flow in saturated–unsaturated soils

1984 ◽  
Vol 21 (3) ◽  
pp. 419-430 ◽  
Author(s):  
A. T. Papagianakis ◽  
D. G. Fredlund
Keyword(s):  
Finite Element ◽  
Steady State ◽  
Pore Water ◽  
Unsaturated Soils ◽  
Pore Water Pressure ◽  
Water Pressure ◽  
Pressure Head ◽  
Flow Equation ◽  
Coefficient Of Permeability ◽  
Steady State Model

A model is proposed describing continuous flow between saturated and unsaturated soil. The flow is assumed to be two dimensional and under steady state conditions. In the unsaturated zone, the coefficient of permeability is treated as a function of pore-water pressure head. The nonlinear differential equation governing the flow is solved using an iterative finite element scheme. The flow equation for an element is derived using the Galerkin weighed residuals method. Several example problems are solved and compared with flow net solutions. The proposed flow model is superior to traditional models, which consider flow only in the saturated zone. The results show that the zero pressure isobar is not an upper flow boundary. The finite element solution is shown to be relatively insensitive to the function used to express the relationship between the coefficient of permeability and the pore-water pressure head. Key words: saturated–unsaturated, pore-water pressure, head, phreatic line.

scholarly journals Finite element analysis for identifying locations of cracking and hydraulic fracturing in homogeneous earthen dams

2021 ◽  
Vol 12 (1) ◽  
Author(s):  
Priyanka Talukdar ◽  
Arindam Dey
Keyword(s):  
Finite Element ◽  
Steady State ◽  
Hydraulic Fracturing ◽  
Pore Water ◽  
Pore Water Pressure ◽  
Water Pressure ◽  
The Body ◽  
Principal Stresses ◽  
Crack Generation ◽  
Earthen Dams

AbstractThis paper reports a finite element study to identify the locations of crack initiation in a homogeneous earthen dam at its post-construction and reservoir operation stages. The steady state and transient analyses, including reservoir rise-up and drawdown conditions, are simulated to identify the favorable conditions and locations of crack generation. The behavioral response of the dam is represented in terms of the developed total stress, horizontal deformations of the faces, strain accumulation in the dam, and the differential settlement of the dam base. The locations of post-construction cracks are identified based on the negative minor principal stresses developed on the dam faces. For both single-lift and multiple-lift modeling techniques, the upstream face of the dam is found to be the most favorable location for the crack generation. For transient reservoir operations (rise-up and drawdown scenarios), it is identified that hydraulic fracturing may occur on either faces of the dam at specific heights, governed by the minor principal stresses becoming lesser than the developed pore-water pressure. Depending on whether the reservoir drawdown occurs before or after the attainment of steady-state phreatic level within the body of the dam, the pore-water pressure distribution within the dam are found notably different. This results in hydraulic fracturing occurring at different faces of the dam and at different heights. It is important to have a thorough understanding of the tentative location of the cracks developed in homogeneous earthen dams so that proper mitigation measure can be adopted as per requirement.

Numerical study of soil conditions under which matric suction can be maintained

2004 ◽  
Vol 41 (4) ◽  
pp. 569-582 ◽  
Author(s):  
L L Zhang ◽  
D G Fredlund ◽  
L M Zhang ◽  
W H Tang
Keyword(s):  
Steady State ◽  
Pore Water ◽  
Pore Water Pressure ◽  
Water Pressure ◽  
Water Storage ◽  
Matric Suction ◽  
Rainfall Infiltration ◽  
Soil Conditions ◽  
Soil System ◽  
Coefficient Of Permeability

The effect of negative pore-water pressure is often ignored in slope stability studies. There is a perception among geotechnical engineers that negative pore-water pressures will dissipate with rainfall infiltration and cannot be relied upon in design considerations. The objective of this paper is to illustrate that under certain conditions soil suction can be maintained. Based on the theory of infiltration and seepage through a saturated–unsaturated soil system, steady state and transient finite element seepage analyses were conducted using Seep/W on a 20 m high slope inclined at 30°. The results of the analysis showed that under steady state conditions, the most important factor influencing the permanency of matric suction in the soil is the magnitude of rainfall flux expressed as a percentage of the saturated coefficient of permeability of soil. For the analysis under transient seepage conditions, the results showed that the pore-water pressure profile depends on the magnitude of the rainfall flux, the saturated coefficient of permeability, the soil-water characteristic curve, and the water storage function. For a soil with a low coefficient of permeability and a large water storage capacity, the matric suction needs a substantial amount of time to dissipate and thus may be maintained over a longer time period than the rain is likely to fall, even if the ground surface flux is equal to or greater than the saturated coefficient of permeability. Engineers should address more appropriate engineering design assumptions that can be related to the permanence of matric suction in soil slopes based on the numerical analysis. Measures such as slope cover or surface recompaction can be taken into consideration to minimize the rainfall infiltration and thus maintain active matric suction in slopes.Key words: unsaturated soils, slope, rainfall infiltration, matric suction, permeability.

scholarly journals Effect of Reduced Zone on Time-Dependent Analysis of Tunnels

2011 ◽  
Vol 2011 ◽  
pp. 1-12
Author(s):  
Mohammed Y. Fattah ◽  
Kais T. Shlash ◽  
Nahla M. Salim
Keyword(s):  
Finite Element ◽  
Pore Water ◽  
Pore Water Pressure ◽  
Water Pressure ◽  
Outer Diameter ◽  
Transient Effects ◽  
Excess Pore Water Pressure ◽  
Modified Cam Clay ◽  
Tunnel Diameter ◽  
Excess Pore

The problem of the proposed “Baghdad metro line” which consists of two routes of 32 km long and 36 stations is analyzed. The tunnel is circular in cross-section with a 5.9 m outer diameter. The finite element analyses were carried out using elastic-plastic and modified Cam clay models for the soil. The excavation has been used together with transient effects through a fully coupled Biot formulation. All these models and the excavation technique together with Biot consolidation are implemented into finite-element computer program named “Modf-CRISP” developed for the purpose of these analyses. The results indicate that there is an inward movement at the crown and this movement is restricted to four and half tunnel diameters. A limited movement can be noticed at spring line which reaches 0.05% of tunnel diameter, while there is a heave at the region below the invert, which reaches its maximum value of about 0.14% of the diameter and is also restricted to a region extending to 1.5 diameters. The effect of using reduced zone on excess pore water pressure and surface settlement (vertical and horizontal) was also considered and it was found that the excess pore water pressure increases while the settlement trough becomes deeper and narrower using reduced .

Aspects of the behaviour of compacted clayey soils on drying and wetting paths

2002 ◽  
Vol 39 (6) ◽  
pp. 1341-1357 ◽  
Author(s):  
Jean-Marie Fleureau ◽  
Jean-Claude Verbrugge ◽  
Pedro J Huergo ◽  
António Gomes Correia ◽  
Siba Kheirbek-Saoud
Keyword(s):  
Water Content ◽  
Pore Water ◽  
Unsaturated Soils ◽  
Pore Water Pressure ◽  
Water Pressure ◽  
Dry Density ◽  
Clayey Soils ◽  
Maximum Dry Density ◽  
Optimum Water Content ◽  
Optimum Water

A relatively large number of drying and wetting tests have been performed on clayey soils compacted at the standard or modified Proctor optimum water content and maximum density and compared with tests on normally consolidated or overconsolidated soils. The results show that drying and wetting paths on compacted soils are fairly linear and reversible in the void ratio or water content versus negative pore-water pressure planes. On the wet side of the optimum, the wetting paths are independent of the compaction water content and can be approached by compaction tests with measurement of the negative pore-water pressure. Correlations have been established between the liquid limit of the soils and such properties as the optimum water content and negative pore-water pressure, the maximum dry density, and the swelling or drying index. Although based on a limited number of tests, these correlations provide a fairly good basis to model the drying–wetting paths when all the necessary data are not available.Key words: compaction, unsaturated soils, clays, drying, wetting, Proctor conditions.

Mini suction probe for matric suction measurements

2002 ◽  
Vol 39 (6) ◽  
pp. 1427-1432 ◽  
Author(s):  
Inge Meilani ◽  
Harianto Rahardjo ◽  
Eng-Choon Leong ◽  
Delwyn G Fredlund
Keyword(s):  
Pore Water ◽  
Unsaturated Soils ◽  
Pore Water Pressure ◽  
Water Pressure ◽  
Matric Suction ◽  
Lateral Direction ◽  
Triaxial Apparatus ◽  
Specimen Height ◽  
Pressure Changes ◽  
Accurate Reading

A modified triaxial apparatus with mini suction probes was fabricated to study the matric suction along the specimen height during unsaturated triaxial testing. Three mini suction probes were placed at 3/4, 1/2, and 1/4 height of the specimen, each at 120° apart in the lateral direction. This paper presents the development of the mini probe for matric suction measurements. Evaluation of the performance shows that the fabricated mini probe provides a rapid response and accurate reading under negative and positive pore-water pressure changes. Matric suctions as high as 400 kPa were successfully measured on soil specimens over a time span of 15 h. On the other hand, the mini suction probes were also found to be able to measure a matric suction of 200 kPa for a longer period of 155 h.Key words: matric suction, mini suction probe, triaxial, unsaturated soils, mid-height pore-water pressure measurement.

Experimental verification of the theory of consolidation for unsaturated soils

1995 ◽  
Vol 32 (5) ◽  
pp. 749-766 ◽  
Author(s):  
Harianto Rahardjo ◽  
Delwyn G. Fredlund
Keyword(s):  
Water Content ◽  
Pore Pressure ◽  
Pore Water ◽  
Unsaturated Soils ◽  
Pore Water Pressure ◽  
Water Pressure ◽  
Matric Suction ◽  
Water Volume ◽  
Loading Tests ◽  
Undrained Loading

An experimental program was designed to study the behavior of unsaturated soils during undrained loading and consolidation. A Ko cylinder was designed and built for the testing program. Simultaneous measurements of pore-air and pore-water pressures could be made throughout a soil specimen using this Ko cylinder. Four types of tests were performed on a silty sand. These are (1) undrained loading tests where both the air and water are not allowed to drain, (2) constant water content tests where only the water phase is not allowed to drain, (3) consolidation tests where both the air and water phases are allowed to drain, and (4) increasing matric suction tests. Undrained loading tests or constant water content loading tests were conducted for measuring the pore pressure parameters for the unsaturated soil. Drained tests consisting of either consolidation tests or increasing matric suction tests were conducted to study the pore pressure distribution and volume change behavior throughout an unsaturated soil during a transient process. The experimental pore pressure parameters obtained from the undrained loadings and constant water content leadings agreed reasonably well with theory. The pore-air pressure was found to dissipate instantaneously when the air phase is continuous. The pore-water pressure dissipation during the consolidation test was found to be faster than the pore-water pressure decrease during the increasing matric suction test. The differing rates of dissipation were attributed to the different coefficients of water volume change for each of the tests. The water volume changes during the consolidation test were considerably smaller than the water volume changes during the increasing matric suction tests for the same increment of pressure change. Key words : consolidation, Ko loading, matric suction, pore-air pressures, pore-water pressures, unsaturated soils

scholarly journals Numerical Simulation Of The Treatment Of Soil Swelling Using Grid Geocell Columns

2015 ◽  
Vol 23 (2) ◽  
pp. 9-18 ◽  
Author(s):  
Mohammed Y. Fattah ◽  
Raid R. Al-Omari ◽  
Haifaa A. Ali
Keyword(s):  
Finite Element ◽  
Pore Water ◽  
Pore Water Pressure ◽  
Water Pressure ◽  
Expansive Soil ◽  
Matric Suction ◽  
Degree Of Saturation ◽  
Axial Forces ◽  
Column Material ◽  
Cap Model

Abstract In this paper, a method for the treatment of the swelling of expansive soil is numerically simulated. The method is simply based on the embedment of a geogrid (or a geomesh) in the soil. The geogrid is extended continuously inside the volume of the soil where the swell is needed to be controlled and orientated towards the direction of the swell. Soils with different swelling potentials are employed: bentonite base-Na and bentonite base-Ca samples in addition to kaolinite mixed with bentonite. A numerical analysis was carried out by the finite element method to study the swelling soil's behavior and investigate the distribution of the stresses and pore water pressures around the geocells beneath the shallow footings. The ABAQUS computer program was used as a finite element tool, and the soil is represented by the modified Drucker-Prager/cap model. The geogrid surrounding the geocell is assumed to be a linear elastic material throughout the analysis. The soil properties used in the modeling were experimentally obtained. It is concluded that the degree of saturation and the matric suction (the negative pore water pressure) decrease as the angle of friction of the geocell column material increases due to the activity of the sand fill in the dissipation of the pore water pressure and the acceleration of the drainage through its function as a drain. When the plasticity index and the active depth (the active zone is considered to be equal to the overall depth of the clay model) increase, the axial movement (swelling movement) and matric suction, as a result of the increase in the axial forces, vary between this maximum value at the top of the layer and the minimum value in the last third of the active depth and then return to a consolidation at the end of the depth layer.

Earthquake-induced deformation analyses of the Upper San Fernando Dam under the 1971 San Fernando Earthquake

2001 ◽  
Vol 38 (1) ◽  
pp. 1-15 ◽  
Author(s):  
Guoxi Wu
Keyword(s):  
Finite Element ◽  
Pore Water ◽  
Effective Stress ◽  
Pore Water Pressure ◽  
Ground Deformation ◽  
Water Pressure ◽  
Soil Liquefaction ◽  
Shear Stresses ◽  
Cyclic Shear ◽  
San Fernando

A nonlinear effective stress finite element approach for dynamic analysis of soil structure is described in the paper. Major features of this approach include the use of a third parameter in the two-parameter hyperbolic stress-strain model, a modified expression for unloading–reloading modulus in the Martin–Finn–Seed pore-water pressure model, and an additional pore-water pressure model based on cyclic shear stress. The additional pore-water pressure model uses the equivalent number of uniform cyclic shear stresses for the assessment of pore-water pressure. Dynamic analyses were then conducted to simulate the seismically induced soil liquefaction and ground deformation of the Upper San Fernando Dam under the 1971 San Fernando Earthquake. The analyses were conducted using the finite element computer program VERSAT. The computed zones of liquefaction and deformation are compared with the measured response and with results obtained by others.Key words: effective stress method, finite element analysis, Upper San Fernando Dam, earthquake deformation, VERSAT.

scholarly journals Finite Element Analysis of Vertical and Horizontal Drainage Structures under Vacuum Combined Surcharge Preloading

2021 ◽  
Vol 2021 ◽  
pp. 1-11
Author(s):  
Daqing Wang ◽  
Dong Wei ◽  
Guoyi Lin ◽  
Jiannan Zheng ◽  
Zhiting Tang ◽  
...  
Keyword(s):  
Finite Element ◽  
Pore Water ◽  
Pore Water Pressure ◽  
Water Pressure ◽  
Horizontal Displacement ◽  
Reinforcement Effect ◽  
Vacuum Suction ◽  
Prefabricated Vertical Drains ◽  
Surcharge Preloading ◽  
Prefabricated Vertical Drain

Combined vacuum and surcharge preloading has gradually been widely used because of its advantages of low cost, green environmental protection, and good treatment effect. The conventional prefabricated vertical drain presents obvious defects in vacuum preloading treatment, such as obvious silting, serious bending of the drainage board, large attenuation of vacuum degree of drainage board along the depth, long construction period, and so on, which affect the final reinforcement effect. In this paper, the MIDAS finite element simulation of combined vacuum and surcharge preloading of prefabricated vertical drains (PVDs) and prefabricated horizontal drains (PHDs) is established through the comparative experiment of the engineering field. The comparative experimental study is carried out from the aspects of the vertical settlement, horizontal displacement, and pore water pressure. The results show that under combined vacuum and surcharge preloading, the consolidation effect of soft soil with PHDs is better than that with PVDs. When PHDs are used, the vertical settlement increases by 7.2% compared with PVDs; the horizontal displacement is larger; and the pore water pressure dissipates faster. This is because when the PHDs are adopted, the consolidation direction of the soil is consistent with the direction of the vacuum suction, which is mainly caused by vertical settlement. With the consolidation, the spacing between PHDs is gradually shortened, and the drainage distance is reduced, which can effectively reduce the consolidation time and improve the reinforcement effect of the soil. In addition, the PHDs can move downward uniformly with the soil during the consolidation process and have almost no bending deformation, which makes the vacuum transfer more uniform and effective.

Sign in / Sign up

Export Citation Format

Share Document