• Rocscience Settle 3d

Download Rocscience Settle 3D Full Patch Free Gratis - Rocscience Settle 3D adalah perngkat lunak 3 dimensi yang digunakan untuk analisis penyelesaian dan konsolidasi di bawah yayasan, tanggul dan penggalian permukaan. Rocscience Settle3D Overview Rocscience Settle3D is an excellent 3D program for the analysis of settlement and consolidation under the foundations and surface [.] The post Rocscience Settle3D Free Download appeared first on Get Into PC.

• Aubrie Malone
• 2 years ago
• Views:

Transcription

1 Settle3D Settlement and consolidation analysis Theory Manual Rocscience Inc.

2

3 Table of Contents Notes.. Stress.. Effective Stress.. Initial stress and pore pressure..3 Stress change due to external load Boussinesq : Multi-layer solution Westergaard solution Stress change due to change in water table elevation Stress change due to excavation Stress change due to rigid load Settlement Immediate settlement Loading Unloading Immediate settlement with mean stress Settlement due to consolidation Linear Non-linear Janbu Koppejan Overconsolidation and underconsolidation Secondary settlement Standard method Mesri method Hydroconsolidation settlement Pore pressures Initial pore pressure Excess pore pressure when load is applied Excess pore pressure above water table Pore pressure dissipation Vertical flow Horizontal flow due to drains Horizontal and vertical flow Permeability and Coefficient of Consolidation Variable permeability Degree of consolidation Buoyancy Buoyancy effect Buoyancy in Settle3D Assumptions Buoyancy examples Infinite embankment, water table at the ground surface i

4 5.4. Finite load, water table at depth Stress correction due to compaction Compaction effect Compaction Correction Examples Infinite embankment, water table at the surface Finite load, water table at depth Empirical Methods Schmertmann Approximation Settlement calculation Influence factors from Schmertmann (970) Modified influence factors from Schmertmann et al. (978) Elastic modulus Peck, Hanson and Thornburn SPT testing Settlement calculation Schultze and Sherif D'Appolonia Method References Table of symbols.. 5 ii

5 Notes Directional quantities such as displacement and stress are vertical unless otherwise noted. No subscripts are used when denoting vertical stress and vertical displacement. In keeping with geotechnical practice, all compressive stresses and strains are positive. The top of the soil horizon always has a vertical (z) coordinate of 0 and vertical coordinates below the surface are positive. Distances and displacements are positive downwards. Stress. Effective Stress Settlement depends on effective stress. The effective stress is the total stress due to gravity and external loads minus the pore water pressure. The vertical effective stress, σ at any point is simply: σ = σ u Where σ is the total vertical stress and u is the pore water pressure. Compressive stresses are positive. Initial stress and pore pressure The initial total stress is just the stress due to gravity loading. At any point this stress is found by summing the weights of the above material layers. The weight of any given layer is the unit weight, γ, times the layer thickness. Therefore the initial total stress at any point is: σ i = γ H Where H is the layer thickness. The unit weight γ is set to either the saturated unit weight or the moist unit weight depending on whether the layer is below or above the piezometric line (water table). The initial pore pressure at any point is the pressure due to the weight of overlying water. If the elevation of the point of interest is z and the water table is at elevation z wt, the initial pressure is: i ( z zwt ) water z zwt u = γ > u = 0 z <= i z wt 3 Where γ water is the unit weight of water.

6 The initial effective stress can then be calculated from equations, and 3: σ = σ u i i i 4 An example initial stress calculation is shown below. 4.3 m γ = 7.8 kn/m 3 γ = 8. kn/m 3. m γ water = 9.8 kn/m 3 Sand 6.0 m Clay γ = 6.0 kn/m 3 A At point A σ i = (. m)(7.8 kn/m 3 ) + (4.3 m. m) (8. kn/m 3 ) + (6 m)(6 kn/m 3 ) = 73.5 kpa u i = (4.3 m. m)(9.8 kn/m 3 ) + (6 m)(9.8 kn/m 3 ) = 89. kpa σ i = 73.5 kpa 89. kpa = 84.3 kpa

7 .3 Stress change due to external load External loads such as embankments or fills can be applied in any stage. These external loads cause a change in total stress depending on the geometry and magnitude of the external load. The stress change can be calculated by four different methods.3. Boussinesq The Boussinesq method uses the theory of elasticity to calculate the vertical stress under a point load in a homogeneous, semi-infinite half space: σ L = 3Q π z cos 5 θ 5 Where σ L is the loading stress at any point and the meaning of the other symbols are as shown. Q Y θ z X σ L Z Useful solutions for stresses under different footing shapes can be obtained by integrating over the area of the footing. 3

8 .3. : The : method assumes that the zone of influence for an applied rectangular load has a slope of : as shown. The vertical loading stress at some depth is then calculated by: Q σ L = 6 ( B + z)( L + z) For a non-rectangular footing, the stress is calculated by computing the area of the load at the surface. With increasing depth, the area over which the load is applied increases at a : ratio and the magnitude of the loading stress decreases correspondingly. L B Q z σ L.3.3 Multi-layer solution The complete elastic solution for an arbitrary shaped foundation resting on a layered elastic medium is computed by the integration of point load solution (Green s functions). Yue (995, 996) has recently provided a computational schema to efficiently compute the point load solutions through Hankel transforms. The elastic state of each homogeneous material body is governed by a set of partial differential equations and a set of boundary conditions. In this case of homogeneous elastic layers the boundary conditions are the set of matching stress and displacement conditions. Hankel transformation is used to convert the set of partial differential equations governing each material layer to a set of ordinary differential equations and the respective boundary conditions to algebraic expressions. Yue (995, 996) solved the systems of ordinary differential equations with the matching boundary conditions to a set of algebraic system of equations, unknowns being specially constructed functions in the Hankel transform domain. The stresses and displacements can be obtained through combinations of the 4

9 inverse transforms of these functions. Due to extreme complexity of the functions in the Hankel transform domain, we need to resort to numerical methods which are computationally intensive and therefore very slow. For our purpose, we need to perform two levels of numerical integration: the first level of computation evaluates the point load solutions, and the subsequent level of computation performs numerical quadrature to evaluate stresses and displacements due to the distributed loading over the given area. We have developed an extremely efficient computational method by combining the two computational levels in different order. We have also enhanced the numerical accuracy by converting the area integrals into boundary integrals by a method developed by Vijayakumar, Yacoub and Curran (000). This boundary conversion method also circumvents the numerical difficulty associated with the singularity of the point load solutions. For any rapidly varying stress locations such as in the vicinity of boundaries, a form of adaptive special subdivision scheme is used.3.4 Westergaard solution The Boussinesq equation is one the most commonly used methods for obtaining the stress distribution induced by an applied vertical load. In developing his point load formula, Boussinesq assumed that a soil medium is elastic, isotropic and homogeneous. However, most soils are neither isotropic nor homogeneous. In fact, many soils encountered in practice exist in layers with distinctive characteristics. In particular, when the soil medium is a layered stratum of fine and coarse material, such as material underneath a road pavement or alternating layers of clay and sand, many experts believe that the Westergaard equations give a better estimate for the stress. Westergaard assumed that an elastic soil medium is laterally reinforced by numerous, closely spaced, horizontal sheets of negligible thickness but of infinite rigidity. These sheets prevent the mass as a whole from experiencing any lateral strains. Hence, only downward vertical deformations are allowed. These assumptions can better represent stratified soils, in which soft layers are reinforced by strong layers. The effect of which is a greater lateral stress distribution than what is usually obtained from methods that assume isotropic elastic solids. Unlike the Boussinesq solution, the Westergaard equations include the Poisson s ratio of the soil medium. For large lateral restraint, the Poisson s ratio may be taken to be zero. In addition, the main effect of Westergaard s assumptions is to reduce the stresses obtained directly below the center of the load as compared to Boussinesq. However, at a certain distance away from the center, stresses obtained from Westergaard begin to exceed those calculated from Boussinesq. For example, for radius r and depth z with a Poisson s ratio ν = 0. 0, Westergaard stresses are about 3 that of Boussinesq when r z = 0. At r z =. 5, both methods give identical values, while for values where r z >. 5, the Westergaard stresses are larger. 5

10 Similar to Boussinesq, closed form expressions for the stress profiles of: ) the center of a uniformly loaded circular area and ) the corner of a uniformly loaded rectangular area are both available. Further, the contributions of adjacent loads can be superimposed to give the total stress at a given point. King of fighters wing 1.8 move list. In general, sedimentary soils like natural clay strata accentuate the non-isotropic condition of the soil medium. Hence, for these cases, the Westergaard equations serve as better models for reality. However, practicing geotechnical engineers often prefer Boussinesq primarily because this method gives more conservative results. In any case, the choice of analysis depends on how closely field conditions match a model s basic assumptions. As described in Venkatramaiah (006), for a soil medium with Poisson s ratio ν, the vertical stress σ due to a point load Q as obtained by Westergaard is given by: z σ = z Q z ν π ν 3 ν r + ν z 7 For large lateral restraint, ν may be taken as zero. The vertical stress below the center of a circular footing can be obtained analytically by integrating equation 7. The solution of which is given by: σ z = q 8 + a ηz where: ν η = and ν q is a uniform load. The vertical stress below the corner of a rectangular footing can be obtained analytically by integrating equation 7. The solution of which is given by: q ν ν σ z = cot π ν m n ν m n where: m = L / z n = W / z 6

11 L and W are the respective lengths and widths of the rectangle z is the depth and q is a uniform load. For the case of a square, L = B. Hence, m = n.4 Stress change due to change in water table elevation If the water table is lowered, then some soil will change from being saturated to unsaturated. This change in weight will cause a decrease in the total stress at all points below the original water table elevation. However, the water table drop will also cause a decrease in pore pressures, and therefore the effective stress will increase by equation. This effect is generally larger than the change in total stress; therefore there is generally a net increase in effective stress. The opposite will occur for a water table rise.5 Stress change due to excavation When material is excavated, the weight of excavated material is calculated from the unit weight and this is applied as an upward stress at the bottom of the excavation. The stress at any point below this is then calculated as if a stress was applied due to a negative load (section.3). There is no attempt to calculate true three-dimensional stress changes due to excavations. Also, no changes in water pressure are considered.6 Stress change due to rigid load We use a point load based (Green s function) error minimization hybrid method for the computation of the stresses resulting from a loaded rigid foundation. The rigid elastic foundation problem is one of the important modeling problems in Civil Engineering practice. The entire displacement of a rigid foundation will follow a linear function with only 3 hitherto unknown variables. On the other hand, the traction on the elastic body will follow an unknown, but much more complicated functionality. The coupling between the rigid body and the elastic body lying below is done by minimizing the average of the square (RMS) of the difference between the displacements. Here, a method using linear piece-wise functions with several parameters via triangular discretizations is adopted (Vijayakumar et al, 009). The elastic response of the piece-wise functions (with the unknown parameters) is calculated by the integration of point load solutions over the triangular discretizations. The point load solutions are available in the analytical form for a homogeneous elastic body and via numerically using the method outlined in Yue (995, 996) for a layered elastic body. The unknown parameters governing the linear functional form in each triangulation are calculated by the overall minimization of the RMS error as mentioned above. Once the parameters are known, the stress and displacement at any point in the elastic body can be calculated using the boundary quadrature method outlined in Vijayakumar et al, (000). By performing several numerical tests, it has been found that the error minimization method is a very efficient and accurate method for the rigid foundation problem. 7

12 3 Settlement The total settlement is the sum of three components: Immediate or initial settlement Settlement due to consolidation Secondary settlement (creep) Each of these is calculated as follows. 3. Immediate settlement 3. Loading The immediate settlement occurs instantly after load is applied and is assumed to be linear elastic. The strain for each element in a string can then be easily calculated from the D modulus and the effective stress. The D modulus (or constrained modulus) (E s ) is input directly into Settle3D. The relationship between the D modulus and 3D Young s modulus (E) is: E = E s ( + ν )( ν ) ν 0 Where ν is the Poisson s Ratio. The vertical strain in each sublayer is calculated by: ε = Δσ E s Where Δ σ is the change in vertical total stress. Initial settlement is then calculated from these strains. For each string, the bottom point is assumed to be fixed (non-moving). The vertical displacement of the point second from bottom is then: δ = Δz = εh Where h is the original thickness of the bottom sublayer. The settlement of the i th point is then the settlement of the point below (i+) plus the settlement in sublayer i: δ + ε h i = δ i + i i 3 8

13 h δ = δ + ε h h δ = δ 3 + ε h h 3 δ 3 = ε 3 h Unloading In Settle3D, the user may supply an unload/reload modulus. If unloading occurs (due to decreasing the magnitude of an existing load, or adding an excavation) then the unload/reload modulus (E sur ) is used in place of E s in equation. If the soil is then reloaded, E sur continues to be used until the soil sublayer reaches its previous stress state, at which time E s is again used in equation to compute strain as shown. σ Loading Slope = E s Unload/reload Slope = E sur γ 3.3 Immediate settlement with mean stress By default, settlement is calculated using only the vertical stress. A more accurate analysis can be performed by using the three-dimensional mean stress in the calculations. The mean stress at any point is the average of the volumetric stress components: ( σ + σ + σ ) or σ = ( σ + σ σ ) σ M = xx yy zz M

14 The Immediate Settlement is then calculated by computing the strain in each sublayer: ε = ( + ν )( Δσ) ν ( Δσ ) 3 M E 5 Where ν is the Poisson s Ratio, Δ σ is the change in total vertical stress, Δ σ M is the change in mean stress and E is the Young s modulus. When mean stress is being used, Young's modulus is input directly for each material (not the D modulus). The unload/reload modulus E ur may be used in place of E as described above. The mean stress is not used in the calculation of Consolidation Settlement and Secondary Settlement because the relationship between strain and mean stress is not clearly defined for non-linear materials. 3. Settlement due to consolidation Settlement due to consolidation progresses gradually as pore pressures dissipate and effective stresses increase. The settlements are calculated from the strains in each sublayer as in equation 3, however the calculation of the strains depends on the type of material as follows. 3. Linear Linear material is assumed to be linear elastic. Therefore the change in vertical strain for any given element for a change in vertical effective stress Δσ is simply: Δε = m Δσ v 6 Where m v is the one-dimensional compressibility. During unload/reload cycles, m v is replaced with the unload/reload compressibility m vur. 3. Non-linear In non-linear material, the modulus is not constant but is a function of stress. This relationship is most commonly shown on a graph of void ratio versus the logarithm of effective stress as shown. 0

15 . Void Ratio Recompression curve Slope = C r Unload-reload curve Slope = C r Virgin curve Slope = C c.6 0 P c Vertical effective stress (kpa) The void ratio in a soil is the ratio of the volume of voids to the volume of solids V e = V v s 7 The stress level P c is the preconsolidation stress and represents the maximum effective stress experienced by the soil in the past. If a soil is experiencing an effective stress less than P c it is overconsolidated. The relationship between the void ratio and the logarithm of the effective stress is given by the recompression index, C r. A soil that has a stress greater than or equal to P c is normally consolidated and its deformation is dictated by the compression index C c. A soil that is unloading is always considered overconsolidated as shown. For a stress change in an overconsolidated soil layer, the change in void ratio, Δe, can be calculated from the initial effective stress, σ i and the final effective stress, σ f by: Δe = C r σ f log σ i 8 Vertical strain is related to void ratio by:

16 Δe ε = + e 0 9 Where e 0 is the initial void ratio. Therefore, Δε = Cr + e 0 σ f log σ i 0 For a normally consolidated soil, the equation is the same with C r replaced with C c. It is also possible for a soil layer to start out overconsolidated and end up normally consolidated if the final stress is greater than P c. In this case, the vertical strain can be calculated by: Δε = Cr + e Pc Cc log + + e σ i 0 0 σ f log P c Note that the initial stress does not necessarily have to refer to the initial in-situ stress due to gravitational loading. For a multi-stage analysis, the change in strain is computed for each stage by using the effective stress at the start of the stage and the effective stress at the end of the stage in equation 0 or. Strains may also be calculated using strain based versions of the compression indices, C cε and C rε. In this case, stress is related directly to strain, rather than void ratio so for a completely overconsolidated material equation 0 becomes: Δε = C r σ f ε log σ i Equations involving normally consolidated material are formed similarly Janbu The Janbu approach (Janbu, 963, 965) can be linear or non-linear depending on the stress exponent, a. The D modulus, M is given by: M σ m σ = = σ r ε σ r a

17 Where m is a modulus number and σ r is a reference stress, usually equal to 00 kpa. The change in strain for any given change in stress is then found by integrating equation. This yields: For a > 0 Δε = ma σ f σ r a a σ i σ r 3 For a = 0 Δε = σ f m ln σ i 4 Where σ i and σ f are the initial and final effective stress for the calculation stage. For a =, the Janbu method is the same as the linear method with m = m σ v r a = When a = 0, the Janbu method is the same as the non-linear method with m + e ln0 C 0 = a = 0 c These equations are written for normally consolidated soil. For overconsolidated soils, m can be replaced by m r, the recompression modulus number. Alternatively, the user can specify a constant modulus M oc for overconsolidated soil. If this option is chosen, then: M = M oc 5 and Δε = M oc σ f σ i 6 3

18 3.4 Koppejan The Koppejan method (Koppejan, 948) computes strain with the following equation: U ε = C p + log C s () t σ f ln σ i 7 Where U is the degree of consolidation, C p and C s are compression constants, t is time, σ i is the initial effective stress and σ f is the final effective vertical stress after loading and consolidation. The degree of consolidation U is computed for any location at any time by: U = u e ue + σ σ i 8 Where u e is the excess pore water pressure and σ is the effective stress. For overconsolidated material, C p and C s are replaced by C p and C s. In Settle3D, an incremental version is used to get the change in strain for any time step, Δt: ΔU Δε = Cp t + Δt σ f + log ln C s t σ i 9 For the first time step, when t equals 0, the logarithm is replaced by log(δt). This equation allows the settlement to be calculated at intermediate times during a settlement analysis. The problem with the Koppejan method is how to deal with changing loads. Traditionally, some type of superposition is used (see van Baars, 003). Since Settle3D computes the strain in an incremental fashion, this is not really necessary. Instead, Settle3D performs the following steps when a change in loading occurs: The current effective vertical stress is recorded and used as σ i in subsequent time steps. The degree of consolidation is calculated using the new σ i Overconsolidation and underconsolidation The Non-linear and Janbu methods use the preconsolidation stress, P c, to determine if the soil is over-consolidated or normally consolidated. Instead of specifying P c you may specify either the overconsolidation ratio (OCR) or overconsolidation margin (OCM). The OCR is related to P c by 4

19 P OCR = c σ 30 Where σ is the current effective stress. The OCM is related to P c by OCM = P c σ (Coduto, 999, p39). 3 Therefore if OCR = or OCM = 0 and there is no load applied, then the material is all normally consolidated and each point in the soil has a P c equal to the in-situ stress due to gravity. When a load is applied, the material remains normally consolidated and P c increases with stress. If OCR > or OCM > 0 then the material is all overconsolidated and at each point in the soil, P c equals OCR times the in-situ stress, or OCM plus the in-situ stress. If a load is then applied, the material remains overconsolidated and P c remains constant until the stress in the soil surpasses P c at which point the material becomes normally consolidated and P c increases as stress increases. If the material is then unloaded, it will then become overconsolidated again and P c will remain constant and equal to the highest stress attained by the soil. If OCR is less than or if OCM is negative, then the soil is underconsolidated. This means that there is still excess pore pressure in the soil due to a previous load that has not yet dissipated. To account for this, Settle3D adds an excess pore pressure to all points in the soil such that u or e ( OCR ) σ = i 3 u e = OCM 33 Where σ i is the initial (in-situ) effective stress due to gravity. P c is then set according to equation 30 or equation 3 and the material behaves as if it is normally consolidated (i.e. C c is used to compute displacements in the non-linear model). As the excess pore pressure dissipates, P c increases as the effective stress increases. Note that a material can only be underconsolidated if the Time-dependent Consolidation Analysis option is turned on in the Project Settings. 3.3 Secondary settlement Secondary settlement, or creep, can occur at a constant effective stress for some types of soils. There are several different methods for computing the magnitude of secondary 5

20 settlement especially when a surcharge load is applied. Settle3D implements two different methods as described in the subsequent sections Standard method The magnitude of the secondary settlement is assumed to vary linearly with the logarithm of time as shown below Void ratio.45 Secondary compression Slope = C α.4 e p.35.3 t p Time Using the void ratio at the end of primary consolidation, e p, the change in secondary vertical strain for a change in time from t to t is Δε S Cα t = log + e t p 34 The secondary displacement can then be calculated by multiplying by the layer thickness as in equations and 3. For linear and Janbu materials, there is no void ratio associated with them. Therefore instead of C α, the user must input the modified secondary compression index or strain based secondary compression index C αε Cα = + e p 35 6

21 So that equation 34 becomes: Δ ε = C S αε log t t To get the time at the end of primary consolidation (t p ) some degree of primary consolidation must be assumed. A primary consolidation of 00% only occurs after infinite time. Therefore a default value of 95% primary consolidation is used as the starting point for secondary consolidation where degree of consolidation is u U = u e e0 36 Where u e is the current excess pore water pressure and u e0 is the initial excess pore water pressure at time of load application (if different loads have been applied at different times, then u e0 is the maximum excess pore water pressure attained). Note that this definition of degree of consolidation is only strictly valid for linear material however it is only being used here as an estimate for the start of secondary consolidation, so perfect accuracy is not crucial. The percent of primary consolidation at which to start secondary consolidation can be set by going to the Advanced options in the Project Settings dialog. Also in this dialog, the user can set the minimum level of stress required to start secondary consolidation. It is not realistic that soil layers subjected to very small stress changes will creep, therefore a minimum stress must be obtained for creep to be activated. By default, this is % of the initial stress due to gravity Mesri method With this method, a value of C / C c is specified. In this way, the secondary consolidation will reflect changes in C c. In Settle3D, C c remains constant, unless the material goes from overconsolidated to normally consolidated or vice versa. It is understood in the Mesri method, that C c in the ratio of C / C c can refer to C r if the material is overconsolidated. So, for simple loading cases, the Mesri method will give exactly the same amount of secondary consolidation as the standard method if C / C c is chosen carefully. The differences arise in the case of surcharge loading. Terzaghi, Peck and Mesri (996) describe how secondary consolidation follows a different pattern when a surcharge (preload) is applied and then removed. When load is removed at time t s, there is primary rebound up to a time of t pr. This is followed by a small amount of secondary rebound to a time of t l. After t l, secondary compression restarts. 7

22 Settlement versus time for a surcharge load, from Terzaghi et al., 996 The postsurcharge secondary compression index, C is not constant with time, as shown in the above figure. It generally starts small and then increases in magnitude. To make secondary settlement calculations easier, a secant index, C is proposed as shown. The secondary compression at any time t after t l can then be calculated from C. The method of calculating secondary compression for a surcharge load proceeds as follows: A surcharge is applied. Consolidation occurs and the maximum effective vertical stress before the surcharge removal is σ vs. The surcharge is removed, leaving behind the desired final load. Rebound occurs and the effective vertical stress at the end of primary rebound is σ vf. The effective surcharge ratio is then calculated: ( ) R = σ / σ s vs vf The value for R s is then used to compute the time to the start of secondary compression t l by: tl t pr = 00 Rs t = 0t R l pr s.7 Soft clays Peat Both t l and t pr are measured relative to the time of load removal, t s. These equations are obtained from laboratory studies as described in Terzaghi et al. (996), Mesri et al. (997) and Mesri et al. (00). It is assumed that no secondary settlement occurs until time t l (the secondary rebound is ignored). To calculate the secondary settlement at any time t after t l, the ratio of C /C is determined from one of the following charts, depending on the material type. Charts are obtained from Terzaghi et al. (996), Mesri et al. (997) and Mesri et al. (00). 8

23 Middleton Peat With the Mesri method, a value of C / C c is generally specified. Therefore the secondary strain at time t is: Δε S = C C C c α α + ep Cc Cα log t t l 39 9

24 3.4 Hydroconsolidation settlement Some soils are very sensitive to moisture content and may experience swelling or collapse due to wetting. In arid and semi-arid climates, natural soils and man-made fills have been known to experience considerable volume reduction as moisture content increases due to irrigation, urbanization etc. Hydroconsolidation in Settle3D accounts for these settlements. Soil collapse, or hydroconsolidation, occurs when a loose clayey sand is exposed to water and the clay bonds break causing significant volume reduction. The amount of collapse depends on the stress. At very low stresses, these soils may actually experience swelling. To quantify the hydroconsolidation behaviour, it is common to perform a series of oedometer tests in which the sample is loaded to a certain stress state and then saturated. The change in void ratio (or strain) due to wetting is then measured (see ε w in the figure below). From a series of such tests, a curve of compaction versus stress can be created. Alternatively, one could simply run two oedometer tests: one with dry material and one with wet. The difference in the two curves can then be used to get a curve of compaction versus stress (see ε w in the figure below). In Settle3D, a table of values for Δε w (change in strain) or Δe (change in void ratio) at different stress states is entered. A hydroconsolidation region is specified. Hydroconsolidation can then occur through two possible mechanisms: 0

25 . The water table is raised. When this happens, the soil will collapse if it is within the hydroconsolidation region that was previously dry, and then becomes wet. Note that soil that is below the water table in Stage is assumed to be collapsed already and will not collapse further. A wetting stage can be specified. This might simulate a broken water pipe or some other wetting mechanism that does not necessarily raise the water table. At this stage and within the specified region, the soil will collapse. The collapse in the soil is calculated by first calculating the stress at the centre of each sublayer in the collapsing region. From this stress, the change in strain is determined and applied according to the user-defined strain versus stress relationship. Therefore an instant collapse (or swelling) is observed. By default, Settle3D assumes unsaturated unit weight for a wetted material above the water table. However you can force the program to consider the saturated unit weight by checking the appropriate box in the hydroconsolidation material parameters dialog. If this option is chosen, the saturated unit weight is only used in calculating the hydroconsolidation settlement. It will not change the stresses plotted in graphs and contour plots and it will not affect calculation of other settlements (immediate, consolidation or secondary).

26 4 Pore pressures 4. Initial pore pressure The initial pore pressure at any point is the pressure due to the weight of overlying water as given in equation Excess pore pressure when load is applied By default, when a load is applied, the pore pressure at each node increases by an amount equal to the change in vertical stress at that node: Δu = Δσ 40 The excess pore pressure is then simply the pore pressure minus the initial pore pressure due to gravity: u = u e u i 4 A more accurate approach is to set the change in pore pressure equal to the change in undrained mean stress. This can be accomplished by choosing the Use mean 3D stress option in the Project Settings (Advanced). The undrained mean stress is the average of the volumetric stress components for an undrained (incompressible) material: σ undrained M undrained ( σ xx + σ yy + σ zz ) or σ M = ( σ + σ + σ 3) = 0. 5 = ν = 0.5 ν An undrained material has a Poisson s ratio ν = 0.5. Setting the initial pore pressure equal to the initial stress is only strictly true for a fluid of infinite stiffness in a fully saturated soil. For partially saturated material, the pore spaces are partly filled with air. In this case the pore fluid (air + water) cannot be considered infinitely stiff. Therefore some load is initially carried by the soil matrix. To account for this effect, the user can specify the Skempton pore pressure coefficient B such that Δu = B undrained ( Δσ ) M (Skempton, 954). The coefficient B ranges from 0 for dry soil to for fully saturated soil. If only vertical stresses are being used, then the parameter B is set such that Δu = B ( ) Δσ By default in Settle3D, B or B =

27 4. Excess pore pressure above water table By default, the soil is assumed to be dry above the water table, so an applied load does not generate excess pore pressure above the water table. If you want the soil to behave as if it were saturated above the water table, then you can select the Generate excess pore pressures above water table checkbox under the Groundwater tab of the Project Settings dialog. In this case, applied loads will generate excess pore pressure above the water table, as if the material were saturated with B =. 4.3 Pore pressure dissipation 4.3. Vertical flow Vertical consolidation is dictated by Terzaghi s D consolidation equation: u t = c v u z e 45 Where u e is the excess pore pressure, c v is the coefficient of consolidation and z is the vertical distance below the ground surface. The excess pressure at any time is calculated from this equation and is then used to calculate the effective stress (equation ). Strains can then be calculated depending on the material type and the strains are used to calculate settlement (equations and 3). This equation can be solved analytically for a single layer with linear stiffness. For multiple layers with differing coefficients of consolidation and different thicknesses, a finite difference approach can be used to solve for the pore pressure. This is the approach adopted by Settle3D. If the problem domain is discretized in time and space as shown below, equation 45 can be replaced by: u u = t c v Δt u i, t+ Δt i, t u c v = ( φ ) [ u + ] + ( ) i, t ui, t ui+, t φ x Δh c v u u + u ( ) [ ] i, t+δt t+δt i+ t+δt Δh i,, 46 Where φ is a time weighting parameter. In Settle3D, φ = 0.5. Note that we have dropped the subscript, e, that indicates excess pore pressure. These equations still refer to excess pore pressure but the subscript has been dropped to simplify the notation. 3

28 time Δt z u i-,t h i- u i,t u i,t+δt h i u i+,t Explicit solution This equation can then be solved implicitly or explicitly. The explicit solution calculates the pressure at time t+δt for each node sequentially. Initial excess pore pressures at time t are calculated for all nodes according to equation 43 or 44. The pore pressure at time t+δt for node i is then calculated by rearranging equation 46 to give: u c Δt ( u + u + u + u ) v i, t+δt = ui, t + i, t i+, t i, t+δt i+, t+ Δt + ( Δh) + Δtcv cvδt + ( Δh) ( Δh) cvδt u + Δtcv i, t Note that the solution for pressure at node i depends on the pressure at node i+, which has not yet been calculated, so the solution requires iterations. In addition, the timestep, Δt, must be small to ensure stability: ( Δh ) Δt < β explicit c v Where β is a dimensionless factor that must be less than 0.5. For accurate results, β is set to 0. in Settle3D. 4

29 Implicit solution For models with long time scales, the explicit solution can be slow, therefore an implicit solution is also implemented in Settle3D. To perform the implicit analysis, equation 46 is rearranged such that all the pressures at time t+δt are on the left hand side: u i, t+δt ( Δh) ( Δh) + ui, t+δt + ui+, t+δt = ui, t + ui, t ui+, t cv t Δ cvδt One equation of this form is written for each node. This leads to a system of linear algebraic equations that can be solved simultaneously using matrix inversion. The system is unconditionally stable, regardless of Δt, however results are poor if Δt is too low: ( h) Δ Δ t > implicit 3 c v (Abid and Pyrah, 988) Choice of solution method In Settle3D, explicit timesteps are executed until the time is greater than the minimum timestep required for the implicit approach. The solution then switches to the implicit approach for the remainder of the solution. As the implicit solution progresses, the timestep is gradually increased. This speeds up the solution without loss of accuracy because of decreasing pore pressure gradients. Non-uniform material Most models will not have a constant value for Δh or c v since different material layers will be present. To account for this, the equations are adjusted slightly. All u i+ terms are multiplied by a factor α i : ki hi α i = ki hi Where k i-, k i are the permeabilities of the sublayers above and below the node, h i-, h i are the thicknesses of the sublayers above and below the node. All u i terms are split into two pieces so for example, in the implicit solution, the term for u i,t+δt becomes: ( Δh) ( Δh) u i t Δt i c t, + + Δ ( c ) α ( Δh ) i+ u t t ( c ) t i +Δ Δ v v t, Δ i v i+ 5

30 Boundary conditions Nodes on the boundaries can be drained or undrained. For a drained boundary, the pressure is set to 0 and the calculations are not performed at this node. An undrained node on the boundary is assumed to be impermeable. To account for this, a dummy node is generated to permit the finite difference calculation to proceed. For example, if the bottom node is impermeable, then a dummy node is created below the bottom node. The distance between the bottom node and the dummy node is the same as the distance between the bottom node and the node second from bottom. The pressure on the dummy node is assigned a pressure equal to the pressure of the node second from bottom. In this way, there will be no flow across the bottom node. Soil u n- Δh Δh u n u n- Impermeable boundary Dummy node 4.3. Horizontal flow due to drains Drains may be added to speed consolidation by permitting horizontal flow. If an array of drains is constructed, then the distance the water must flow to the nearest drain is short and the pore pressure dissipation will be greatly accelerated. In addition, the horizontal permeability is often greater than the vertical permeability so horizontal flow is faster than vertical flow. In general, an array of drains is constructed in a square or triangular pattern. 6

31 drain D D D e D e Square pattern Triangular pattern These drain patterns ensure that the entire area within the grid is being drained. Each drain therefore has associated with it a zone of influence. It is assumed that horizontal flow occurs only toward the drain within this zone of influence, i.e. the edge of each zone of influence is an impermeable boundary. For a triangular array, the zone of influence is hexagonal and for a square array, the zone is square. To simplify the mathematics, each zone of influence is assumed to be circular and the diameter of the equivalent circle is calculated as follows: D e =.3D D e =.05D Square pattern Triangular pattern Where D is the actual spacing between the drains as shown. For each drain, the pore pressure at any given point in the zone of influence at any given time can be calculated analytically using the radial flow equation. The average pore pressure within the zone of influence can then be calculated at any time. This average value will be the same for all drains, therefore everywhere within the drain array is assumed to have the same pore pressure at any given time. To make the equations more realistic, a smear zone can be included. When a drain is constructed, the soil close to the drain walls is disturbed forming a smear zone. The permeability in the smear zone may be lower than the permeability of the undisturbed soil and this will affect the pore pressure dissipation. The geometry of the problem for a single drain is shown below: 7

32 Drain Clay Smear zone d d s D e z l The excess pore pressure within the drain array is then calculated by (Barron, 948): u e = u e0 8T exp μ r 47 ch Where T r =, c h is the horizontal consolidation coefficient and μ is a function of the De drain geometry: μ = n = n D d e n S n S kh n S ln S 4n ks n ln S 48 8

33 S = d s d kh And is the ratio of horizontal permeability in the undisturbed zone to permeability in ks kh the smear zone. A common value of S is about.5 3 (Hansbo, 979). If S = and k = s, then there is no smear zone. Equivalent diameter of band shaped drains Wick drains are generally not cylindrical as assumed in the above equations. Kjellman (948) showed that the draining effect of a drain depends on the circumference of its cross section and not on the cross-sectional area. Therefore the equivalent diameter of a band shaped drain with width b and thickness t is d = ( b + t) π 49 Well resistance The above equations assume that the drain itself has an infinite permeability. In fact this is not the case and for long drains, the well resistance may slow down fluid flow. To account for this, μ in equation 47 should be replaced by μ r (Hansbo, 98) ( ) μ = μ+ πz l z k / q r h w Where l is the length of the well, z is the distance from the top of the well, k h is the horizontal permeability of the clay and q w is the discharge capacity of the well (volume / time). If the well drains at both the top and the bottom (i.e. the bottom intersects a highly permeable layer), the l is set to half the length of the drain. A Note on permeability Although equation 47 uses c h to compute the pore pressure decrease due to drains, most practitioners will likely not have this information. In Settle3D, the user may choose to enter the permeability, k h instead. The value for c h will then be calculated from the vertical stiffness (see section 4.4). This is not strictly correct since the horizontal stiffness should be used in the calculation, however it is thought that the error introduced with this method will be small compared to the uncertainty in the input permeability. 9

34 4.3.3 Horizontal and vertical flow The pore pressure at any time is a function of pressure loss due to vertical drainage and pressure loss due to horizontal flow into drains. To compute the total pore pressure at each time step, the pressure due to vertical flow is calculated first using the finite difference approach in equation 46. This pressure is then used as the starting pressure for the calculation of equation 47 that accounts for pressure loss due to vertical drains. The result of equation 47 is then the ultimate pressure for the timestep. 4.4 Permeability and Coefficient of Consolidation To solve for pore pressures in a multi-layered material, both the material permeability (k) and material coefficient of consolidation (c v ) are required. However in Settle3D, the user only needs to specify one or the other of these quantities. This is because these quantities are not independent and are related to each other through the material stiffness. The material stiffness depends on the material type. Therefore if the user enters c v, permeabilities are calculated as follows: k = c v m v γ w linear material k k cvccγ w =.3 ( + e ) σ c γ 0 zi v w = a zi m σ σ r σ r non-linear material Janbu material 50 cvγ w k = C σ p zi Koppejan material Where σ zi is the initial effective stress and γ w is the unit weight of pore water. For overconsolidated material, c v is replaced by c vr, C c is replaced by C r, m is replaced by m r and C p is replace by C p. If the user enters k instead, then equation 50 can be rearranged to obtain c v Variable permeability It is often observed that increasing stress in a soil causes a permeability decrease as porosity (and void ratio) decreases. Several relationships have been proposed to mathematically account for changing permeabilities with changing stress. Settle3D currently implements two different relationships as described below. 30

35 The changing permeability is accounted for in Settle3D by computing stress (and/or void ratio) at the start of each stage and setting the permeability for the duration of the stage. Permeability does not change during a stage. If the material permeability is very sensitive to stress changes and the stress is changing rapidly, then the user should specify many stages of short duration to ensure accurate results. Changing permeability with changing void ratio Terzaghi et al. (996) give the following relationship: Δe C k = 5 Δ logk Where e is the void ratio, k is the permeability and C k is a unitless parameter to be specified by the user. This equation cannot be used for Linear or Janbu material since no void ratio is specified for these material types. It is up to the user to specify C k, however Tavenas et al. (983) propose the following empirical relationship between C k and e for soft clay and silt deposits: C k e 0 = Permeability is a function of effective stress Vaughan (989) proposed that the permeability, k, varies exponentially with the mean effective stress, σ : ( σ ) k = k exp B 0 Where k 0 is the permeability at zero mean effective stress and B is a user-specified material property with units of AREA/FORCE. This B should not be confused with the Skempton pore pressure coefficient a completely different parameter. In Settle3D, vertical stress is used instead of mean stress in equation Degree of consolidation The degree of consolidation is a number between 0 00% that indicates the stage of consolidation (0 for completely unconsolidated to 00% for totally consolidated). Lambe and Whitman (969) give two definitions for degree of consolidation. The first is: compression at time t U = compression at end of consolidation 54 This is the equation used by Settle3D, where consolidation settlements are used in the calculation (not total settlements). 3

36 However, for some situations, namely when B <, then this may be unsatisfactory. When B <, or when the mean stress is used to calculate settlements, the applied load is not completely balanced by excess pore pressures so there is some instant settlement. In these cases you will see the degree of consolidation start from some non-zero value and progress to 00% as excess pore pressures dissipate. This may not be what you want, so Settle3D also provides another indicator for U called Average Degree of Consolidation. This follows the second method of calculation in Lambe and Whitman (969) and the calculation is performed as follows: First, at each point below the surface, the consolidation ratio is calculated: U z = u u e peak Where u e is the excess pore pressure and u peak is the maximum excess pore pressure in the past. These values are then summed over the thickness of a material layer as shown in the figure below. In this way, an average degree of consolidation can be calculated for each material layer. z given t shaded area U = total area U z For a linear material with B =, the values for U and U avg are equal. In other cases they are not. It is up to the user to decide which indicator is most appropriate. Keep in mind however, that U reflects the consolidation of all soil below the point of interest, whereas U avg reflects the consolidation of each material layer independently. 3

37 5 Buoyancy 5. Buoyancy effect In practice, as settlement occurs, the stresses at a given point in the soil will change due to a buoyancy effect. The buoyancy effect can be described as follows: As a point moves downwards, it descends further below the water table. This causes a pore pressure increase at the point and therefore an effective stress decrease. In addition, if soil that was originally above the water table settles below the water table, then the soil becomes saturated and its unit weight may increase. This causes an effective stress increase at all points below the submerged soil. 5. Buoyancy in Settle3D In Settle3D, the buoyancy effect can be accounted for by selecting the Include buoyancy effect checkbox found in the Advanced option under the General tab of the Project Settings dialog. When this option is turned on, Settle3D will account for buoyancy by following these steps for each stage:. Settlement is calculated due to changes in applied load, changes in the water table, excavations or changes in excess pore pressure (for time-dependent analysis). The change in stress due to buoyancy effects is computed. 3. This stress change is applied and the settlement is recomputed. 4. Steps and 3 are repeated until the difference in settlement between two iterations falls below some tolerance. In Settle3D, the total settlement is the sum of immediate settlement, consolidation settlement and secondary settlement (creep). The way that each of these is affected by buoyancy is described below: Immediate settlement: Immediate settlement only depends on total stress, not effective stress. Therefore pore pressure changes due to settlement have no effect on the immediate settlement. However, changing unit weights of material descending below the water table will affect the immediate settlement. Also, the immediate settlement affects the consolidation settlement since changing pore pressures do influence the consolidation settlement. Consolidation settlement: Changes in pore pressure due to immediate, consolidation and secondary settlement will change the effective stress and will therefore influence the consolidation settlement. Secondary settlement: Buoyancy does not affect secondary settlement because the amount of secondary settlement does not depend on stress only time. Note however, that the secondary settlement will affect the calculation of consolidation settlement since secondary settlement produces pore pressure (and therefore effective stress) changes. 33

38 5.3 Assumptions. When material settles below the water table and its unit weight changes, this stress change does not affect all underlying points equally since the settlement does not have an infinite horizontal extent. The stress change at each point below the water table is therefore calculated by multiplying the stress change by an influence factor. The influence factor is the influence factor of the most recently applied load. This assumption may lead to inaccuracies if the most recently applied load is not the load that is having the most effect on the settlement. In a time-dependent analysis, the buoyancy effect should increase gradually as settlement occurs with time. In Settle3D, the buoyancy effect is only calculated at the end of each stage, therefore small inaccuracies will result. To improve the accuracy of the time-dependent results, you could include more stages in the analysis. 3. The water table can never be higher than the current ground surface elevation (except for the case of embankment loads, see number 4). This means that a water table at the surface will actually move downwards as settlement occurs to maintain the water table at the top of the soil. The stress effect of this moving water table is multiplied by influence factors as in assumption. 4. For embankment loads, if the embankment settles partially below the water table, the water table will extend into the embankment material and therefore there will be no change in the water table height. If the embankment settles completely below the water table, the water table cannot be higher than the top surface of the embankment, therefore the water table is lowered as in assumption As settlement occurs, it is assumed that the density of material increases such that the compaction of a material layer does not change its weight. To account for decreasing soil weight due to compaction, see section 6. 34

39 5.4 Buoyancy examples 5.4. Infinite embankment, water table at the ground surface Assume an embankment of infinite extent, thickness = t, unit weight = γ fill The soil has a saturated and unsaturated unit weight both equal to γ soil The water table is at the surface and the water unit weight is γ w γ fill t δ γ soil A h δ At point A prior to settlement, u = hγ w At point A after settlement, ( h ) γ w u = + δ We assume no change in weight of overlying material, therefore the change in effective stress due to settlement is simply the change in pore pressure: Δ ( u u ) = δ γ w σ = 35

40 5.4. Finite load, water table at depth Assume a circular load of radius r and magnitude F. The soil has a moist unit weight of γ moist above the water table, and a saturated unit weight of γ sat. The water table is at a depth of h and the water unit weight is γ w Load F, radius = r γ moist h γ sat A δ h δ B At point B the change in pore pressure is, Δ u = δ γ w And the change in stress due to changing weight of overlying material is: σ = δ ( γ sat γ moist )I Δ where I represents the influence factor that accounts for the finite load. (Recall that we are assuming the density of material increases to offset the decreased thickness) Therefore the change in effective stress due to settlement is: Δ = Δσ Δu = δ ( γ sat γ moist ) I δ γ w σ At a distance, z beneath the centre of a circular load of radius r, the influence factor is I () z = + ( r / z) 3 36

41 6 Stress correction due to compaction 6. Compaction effect As settlement occurs, the stress at a given point in the soil may decrease as the above layers compact and their weight decreases. In Settle3D, you can account for this effect by selecting the Include vertical stress reduction due to settlement above a point checkbox found in the Advanced option under the General tab of the Project Settings dialog. The compaction effect also includes the effect of soil above the water table settling below the water table and therefore experiencing a change in unit weight. Settle3D corrects for compaction by following these steps for each stage:. Settlement is calculated due to changes in applied load, changes in the water table, excavations or changes in excess pore pressure (for time-dependent analysis). The change in stress due to compaction effects is computed. 3. This change in stress is applied and the settlement is recomputed. 4. Steps and 3 are repeated until the difference in settlement between two iterations falls below some tolerance. The same analysis is performed for immediate settlement and consolidation settlement. Compaction does not affect secondary settlement because the amount of secondary settlement does not depend on stress only time. Note however, that the secondary settlement WILL affect the consolidation settlement. In Settle3D, the following assumptions are made:. If the compaction correction is turned on, it is assumed that the density of soil does not change as it compacts. In reality there will probably be some increase in density therefore Settle3D will tend to overcorrect for compaction. For a load of finite extent, the stress change does not affect all underlying points equally since the settlement does not have an infinite horizontal extent. The stress change at each point below the load is therefore calculated by multiplying the stress change by an influence factor. The influence factor is the influence factor of the most recently applied load. This assumption may lead to inaccuracies if the most recently applied load is not the load that is having the most effect on the settlement. 3. In a time-dependent analysis, the compaction effect should increase gradually as settlement occurs with time. In Settle3D, the compaction effect is only calculated at the end of each stage, therefore small inaccuracies will result. To improve the accuracy of the time-dependent results, you could include more stages in the analysis. 37

42 6. Compaction Correction Examples 6. Infinite embankment, water table at the surface Assume an embankment of infinite extent, thickness = t, unit weight = γ fill The soil has a saturated and unsaturated unit weight both equal to γ soil The water table is at the surface and the water unit weight is γ w γ fill t δ γ soil A h δ At point A prior to settlement, σ = tγ + hγ fill soil At point A after settlement, ( h + δ δ ) γ soil σ = tγ fill + We assume there is no change in pore pressure, therefore the change in effective stress due to settlement is: ( δ δ ) γ soil Δ σ = σ σ = 38

43 6. Finite load, water table at depth Assume a circular load of radius r and magnitude F. The soil has a moist unit weight of γ moist above the water table, and a saturated unit weight of γ sat. The water table is at a depth of h Load F, radius = r γmoist h δ γ sat h δ A At point A prior to settlement, σ = FI + hγ moist + h γ sat Where I is the influence factor that accounts for the finite load At point A after settlement, we can calculate the stress state as if the settlement was infinite in horizontal extent, ( h δ) γ moist + ( h δ ) γ sat inf inite σ = FI + + The difference between these two stresses is inf inite σ σ = δγ moist + δ γ sat We assume there is no change in pore pressure. Therefore to get the change in effective stress due to settlement, multiply the above difference by the influence factor: Δ = ( δ γ δ )I σ sat γ moist At a distance, z beneath the centre of a circular load of radius r, the influence factor is I () z = + ( r / z) 3 39

44 7 Empirical Methods A suite of empirical methods exist for calculating immediate settlement in cohesionless soil. These methods are generally based on data from field tests (Standard Penetration Test, Cone Penetration Test, etc.) and require little user input. In general, settlement can only be calculated for regular shaped loads (rectangles and circles). Several of these methods are implemented in Settle3D as described below. 7. Schmertmann Approximation The Schmertmann method (970) calculates settlement from layer stiffness data or cone tip bearing resistances, q c obtained from a Cone Penetration Test (CPT). The method proposed a simplified triangular strain distribution and calculates the settlement accordingly. A time factor can also be included to account for time dependent (creep) effects. 7. Settlement calculation The equation for settlement is: δ = C C ΔP t n i = Δzi I E si zi Where od C = the correction to account for strain relief from excavated soil, σ ΔP σ od = effective overburden pressure at bottom of the footing ΔP = the net applied footing pressure (σ L σ od ) C t = correction for time-dependent creep, + 0. log( 0t) t = time (years) E si = one-dimensional elastic modulus of soil layer i Δz i = thickness of soil layer i I zi = the influence factor at the centre of soil layer i as described below. 55 The influence factor, I z is based on an approximation of strain distributions below the footing. Two different ways have been proposed to calculate these factors. 7. Influence factors from Schmertmann (970) The simplest approach is to assume that the strain influence factor I z increases linearly from zero at the bottom of the footing to a maximum of 0.6 at a depth of ½B below the footing where B is the footing width. The strain influence factor then decreases linearly to zero at a depth of B below the footing bottom. This distribution is shown below. 40

45 B 0 D B/ B Strain Influence factor, I z Strain influence factors from Schmertmann (970) 7.3 Modified influence factors from Schmertmann et al. (978) A modified version of the Shertmann approach considers the shape of the load when calculating the strain influence factors. With this approach, the peak value for I z is calculated by I zp ΔP = σ op 56 Where σ op is the overburden pressure at the depth of I zp. The depth of I zp depends on the shape of the load. For an axisymmetric load (a circle or a square), I zp occurs at a depth of B/. For the plane strain case (length of the load is > 0 the width), the I zp occurs at a depth of B. The values for I z (shown in the figure below) are calculated as follows: Axisymmetric: I z varies linearly from 0. at the bottom of the footing to I zp at a depth of B/. The strain influence factor then decreases to zero at a depth of B. Plane strain: I z varies linearly from 0. at the bottom of the footing to I zp at a depth of B. The strain influence factor then decreases to zero at a depth of 4B. 4

46 In Settle3D, I z is calculated using the axisymmetric equations for circle and square loads. For rectangular loads in which the length is greater than ten times the width, the plane strain approach is used. For rectangular loads in which the length is less than ten times the width, a linear interpolation between the axisymmetric and plane strain case is performed, dependent on the length to width ratio. B 0 D B/ B Axisymmetric B Plane strain 4B Strain Influence factor, I z Strain influence factors from Schmertmann et al. (978) Elastic modulus The elastic modulus E s can be estimated from the results of a Cone Penetration test: E s =.0q c Schmertmann (970) E s =.5q c E s = 3.5q c Modified Schmertmann (978), Axisymmetric Footing Modified Schmertmann (978), Plane strain footing Where q c is the cone tip bearing resistance. As with the calculation of the strain influence factor, the value for E s is calculated to be between the axisymmetric case and plane strain 4

47 case if the length of the load is less than ten times the width (for the modified Schmertmann calculations). 7. Peck, Hanson and Thornburn The method of Peck, Hanson and Thornburn (974) uses the results of a Standard Penetration Test (SPT) to obtain settlement. The water table and overburden pressure at the location of the SPT are taken into account. Settlement results are obtained from matching the problem geometry and corrected SPT results to empirical curves. 7. SPT testing The Standard Penetration Test produces a value, N, equal to the number of hammer blows needed to drive an SPT sample foot through the soil. Standards dictate the details of the procedure. To use the Peck, Hanson and Thornburn method, the N value should be corrected to account for the efficiency of your testing system. Peck et al., assume that a standard test is being performed in which 60% of the energy is transmitted to the soil. If your system is more or less efficient, a correction should be made: N 60 E m = N Where, N 60 = the SPT N value corrected for field testing procedures N is the recorded number of blows per foot, E m is the hammer efficiency in percent (generally between 45-95) Skempton (986) proposes further corrections to account for the rod length and borehole diameter: Em N60 = CbCrN Where, C b is the borehole diameter correction (C b =.0 for borehole diameter less than 5 mm,.05 for borehole diameters between 5 and 50 mm and.5 if the diameter is > 50 mm) C r is the rod length correction (C r = 0.75 for less than 4 m of drill rods, 0.85 for 4-6 m of drill rods, 0.95 for 6-0 m, and.0 for drill rods more than 0 m) Other corrections are often performed to account for overburden pressure and the water table but these are generally method dependent and will be described separately for each method. 43

48 7. Settlement calculation The settlement in inches is calculated by: δ = σ L Cw ΔP Where σ L is the applied load, ΔP is the load required to produce a settlement of inch and C w is the correction factor for water table depth. The water table correction C w is given by: 59 C zw * C D + B w = w 60 Where z w is the depth to the water table D is the depth to the bottom of the load B is the load width. The value for ΔP is obtained from the empirical charts shown below. To use the charts, a corrected blow count is required: ( N60 ) = N 60Cn 6 Where C n is the correction for overburden pressure. The overburden correction is given by: 0 C 0.77log n = C n σ od 6 Where σ od is the initial effective stress in tons/ft due to overburden soil. The value of (N 60 ) represents the average value between the bottom of the footing and a depth of B below the bottom of the footing. So, if there are multiple layers, the overburden stress is calculated for the midpoint of each layer and the correction is calculated by equation 6. (N 60 ) is then calculated for each layer. The average (N 60 ) between depths of D and D+B is then calculated and used to get ΔP from the charts below. Once you have obtained (N 60 ), use this value along with your footing width B and ratio of footing depth to width D/B to obtain ΔP from the charts below. This is then used in equation 59 to get the settlement. Settle3D automatically interpolates between the curves to obtain ΔP. 44

49 Charts used for obtaining ΔP, the load required to induce a settlement of inch. The N value refers to the blow counts corrected for efficiency, overburden and water table. Reproduced from Peck, Hanson and Thornburn, Schultze and Sherif This method, proposed by Schultze and Sherif (973) is based on the equation for settlement in an elastic isotropic half-space: σ LB δ = E s f n 63 Where σ L is the loading stress, B is the load width, E s is the one-dimensional modulus and f n is an influence factor that is a function of: ν Poisson s ratio L/B load length divided by load width H/B thickness of compressible layer divided by load width The geometry of the problem is shown below. The contribution of Schultze and Sherif was to write this equation in terms of the blow counts, N 60 instead of E s. Using empirical data, the following equation was derived: δ = σ L f n B 60 + / 0.87 ( N ) ( 0.4D B) 64 The value for f n is taken from the charts shown below. This assumes that the Poisson s ratio is 0. It is assumed that below a depth of B the influence of the load is minimal so for values of H/B >, the value of f n at H/B = is used. 45

When to Use Immediate Settlement in Settle 3D

When to Use Immediate Settlement in Settle 3D Most engineers agree that settlement is made up of three components: immediate, primary consolidation and secondary consolidation (or creep). Most engineers

More information

9. COMPRESSIBILITY AND SETTLEMENT

9-1 9. COMPRESSIBILITY AND SETTLEMENT 9.1 COMPRESSIBILITY OF AN ELASTIC SOLID Structures such as buildings, bridges, etc. are frequently erected on soil as the founding material. Since soil is a compressible

More information

INTRODUCTION TO SOIL MODULI. Jean-Louis BRIAUD 1

INTRODUCTION TO SOIL MODULI By Jean-Louis BRIAUD 1 The modulus of a soil is one of the most difficult soil parameters to estimate because it depends on so many factors. Therefore when one says for example:

More information

4.3 Results.. 27 4.3.1 Drained Conditions.. 27 4.3.2 Undrained Conditions.. 28 4.4 References.. 30 4.5 Data Files.. 30 5 Undrained Analysis of

Table of Contents 1 One Dimensional Compression of a Finite Layer.. 3 1.1 Problem Description.. 3 1.1.1 Uniform Mesh.. 3 1.1.2 Graded Mesh.. 5 1.2 Analytical Solution.. 6 1.3 Results.. 6 1.3.1 Uniform

More information

Stress Analysis Verification Manual

Settle3D 3D settlement for foundations Stress Analysis Verification Manual 007-01 Rocscience Inc. Table of Contents Settle3D Stress Verification Problems 1 Vertical Stresses underneath Rectangular Footings

More information

Secondary Consolidation and the effect of Surcharge Load

Secondary Consolidation and the effect of Surcharge Load Thuvaragasingam Bagavasingam University of Moratuwa Colombo, Sri Lanka International Journal of Engineering Research & Technology (IJERT) Abstract

More information

Applying a circular load. Immediate and consolidation settlement. Deformed contours. Query points and query lines. Graph query.

Quick Start Tutorial 1-1 Quick Start Tutorial This quick start tutorial will cover some of the basic features of Settle3D. A circular load is applied to a single soil layer and settlements are examined.

More information

FUNDAMENTALS OF CONSOLIDATION

FUNDAMENTALS OF CONSOLIDATION SAND (Vertical Stress Increase) CLAY CONSOLIDATION: Volume change in saturated soils caused by the expulsion of pore water from loading. Saturated Soils: causes u to increase

More information

White Paper: Establishing and Investigating Foundation Zones of Influence

White Paper: Establishing and Investigating Foundation Zones of Influence Prepared by: Bracken Engineering, Inc. 2701 W. Busch Blvd, Ste 200 Tampa, Florida 33618 December 1, 2012 (Revised: December 2014)

More information

Estimation of Compression Properties of Clayey Soils Salt Lake Valley, Utah

Estimation of Compression Properties of Clayey Soils Salt Lake Valley, Utah Report Prepared for the Utah Department of Transportation Research Division by Steven F. Bartlett, PhD. P.E. Assistant Professor

More information

Embankment Consolidation

Embankment Consolidation 36-1 Embankment Consolidation In this tutorial, RS2 is used for a coupled analysis of a road embankment subject to loading from typical daily traffic. Model Start the RS2 9.0 Model

More information

Consolidation and Settlement Analysis

19 Consolidation and Settlement Analysis Patrick J. Fox Purdue University 19.1 Components of Total Settlement 19.2 Immediate Settlement 19.3 Consolidation Settlement Total Consolidation Settlement Rate

More information

UU Triaxial Test. Please refer the same materials you got in Unconfined Compression test.

1 Civil & Environmental Engineering Department UU Triaxial Test Concept of Shear Strength Please refer the same materials you got in Unconfined Compression test. a. Triaxial Shear Test Figure 1 triaxial

More information

CE 366 SETTLEMENT (Problems & Solutions)

CE 366 SETTLEMENT (Problems & Solutions) P. 1) LOAD UNDER A RECTANGULAR AREA (1) Question: The footing shown in the figure below exerts a uniform pressure of 300 kn/m 2 to the soil. Determine vertical

More information

Module 1 (Lecture 3) GEOTECHNICAL PROPERTIES OF SOIL AND OF REINFORCED SOIL

Module 1 (Lecture 3) GEOTECHNICAL PROPERTIES OF SOIL AND OF REINFORCED SOIL Topics 1.1 CAPILLARY RISE IN SOIL 1.2 CONSOLIDATIONS-GENERAL 1.3 CONSOLIDATION SETTLEMENT CALCULATION 1.4 TIME RATE OF CONSOLIDATION

More information

Period #16: Soil Compressibility and Consolidation (II)

Period #16: Soil Compressibility and Consolidation (II) A. Review and Motivation (1) Review: In most soils, changes in total volume are associated with reductions in void volume. The volume change of the

More information

Drained and Undrained Conditions. Undrained and Drained Shear Strength

Drained and Undrained Conditions Undrained and Drained Shear Strength Lecture No. October, 00 Drained condition occurs when there is no change in pore water pressure due to external loading. In a drained

More information

Groundwater Analysis Overview

Finite Element Groundwater Seepage (Overview) 7-1 Groundwater Analysis Overview Introduction Within the Slide program, Slide has the capability to carry out a finite element groundwater seepage analysis

More information

VERTICAL STRESS INCREASES IN SOILS TYPES OF LOADING

14.53 THEORETICAL SOIL MECHANICS VERTICAL STRESS INCREASES IN SOILS TYPES OF LOADING Point Loads (P) Line Loads (q/unit length) Figure 6.11. Das FGE (005). Examples: -Posts Figure 6.1. Das FGE (005). Examples:

More information

Flat Plate Dilatometer Test & TEXAM Pressuremeter Tests in Soil Investigations with Special Emphasis on Evaluating Settlements

Flat Plate Dilatometer Test & TEXAM Pressuremeter Tests in Soil Investigations with Special Emphasis on Evaluating Settlements David L. Tucker, P.E., Regional Technical Expert Deep Foundations U.S. Army

More information

VERTICAL STRESS INCREASES IN SOIL TYPES OF LOADING. Point Loads (P) Line Loads (q/unit length) Examples: - Posts. Examples: - Railroad track

VERTICAL STRESS INCREASES IN SOIL Point Loads (P) TYPES OF LOADING Line Loads (q/unit length) Revised 0/015 Figure 6.11. Das FGE (005). Examples: - Posts Figure 6.1. Das FGE (005). Examples: - Railroad

More information

Total water potential is sum of gravitational and soil water pressure potential, or. Column of water

SSC107 Fall 000 Chapter 3 Page 3-1 Chapter 3. Saturated Water Flow All pores are filled with water, i.e., volumetric water content is equal to porosity ( θ = θ s with θ s = φ ) Nonequilibrium. Water flows

More information

ENCE 4610 Foundation Analysis and Design

This image cannot currently be displayed. ENCE 4610 Foundation Analysis and Design Shallow Foundations Total and Differential Settlement Schmertmann s Method This image cannot currently be displayed. Strength

More information

INDIAN INSTITUTE OF TECHNOLOGY GANDHINAGAR Department of Civil Engineering Soil Mechanics Laboratory

THEORY: INDIAN INSTITUTE OF TECHNOLOGY GANDHINAGAR CONSOLIDATION TEST- (OEDOMETER TEST) (IS 2720-PART-15-1986) Reaffirmed-2002 When a compressive load is applied to soil mass, a decrease in its volume

More information

A Theoretical Solution for Consolidation Rates of Stone Column-Reinforced Foundations Accounting for Smear and Well Resistance Effects

The International Journal of Geomechanics Volume, Number, 135 151 (00) A Theoretical Solution for Consolidation Rates of Stone Column-Reinforced Foundations Accounting for Smear and Well Resistance Effects

More information

Prof. B V S Viswanadham, Department of Civil Engineering, IIT Bombay

37 Module 6: Lecture -2: Buried Structures Content in this module: Load on Pipes, Marston s load theory for rigid and flexible pipes, Trench and Projection conditions, minimum cover, Pipe floatation and

More information

LIST OF SYMBOLS. - diameter of drilled shaft or width of spread foundation; dilatometer B reading

LIST OF SYMBOLS ENGLISH LETTERS - UPPER CASE A A tip A v B CIUC CU CPT COV COV F COV G COV K COV Ko COV Q COV R COV e COV su COV V COV w COV ε COV ξa - dilatometer A reading - tip or base area - area exposed

More information

Axisymmetry Tutorial. A few representations of simple axisymmetric models are shown below. Axisymmetry Tutorial 6-1 0, 18 28, MPa 0, 12.

Axisymmetry Tutorial 6-1 Axisymmetry Tutorial 0, 18 28, 18 0, 12 10 MPa 10 MPa 0, 6 0, 0 4, 0 x = 0 (axis of symmetry) user-defined external boundary 4, 24 12, 24 20, 24 28, 24 This tutorial will illustrate

More information

11. THE STABILITY OF SLOPES

11-1 11. THE STABILITY OF SLOPES 11.1 INTRODUCTION The quantitative determination of the stability of slopes is necessary in a number of engineering activities, such as: (a) (b) (c) (d) the design of earth

More information

Materials & Loading Tutorial

Materials & Loading Tutorial 2-1 Materials & Loading Tutorial This tutorial will demonstrate how to model a more complex multimaterial slope, with both pore water pressure and an external load. MODEL FEATURES:

More information

EML 5526 FEA Project 1 Alexander, Dylan. Project 1 Finite Element Analysis and Design of a Plane Truss

Problem Statement: Project 1 Finite Element Analysis and Design of a Plane Truss The plane truss in Figure 1 is analyzed using finite element analysis (FEA) for three load cases: A) Axial load: 10,000

More information

Seismic Design of Shallow Foundations

Shallow Foundations Page 1 Seismic Design of Shallow Foundations Reading Assignment Lecture Notes Other Materials Ch. 9 FHWA manual Foundations_vibrations.pdf Homework Assignment 10 1. The factored forces

More information

Constant Rate of Strain Consolidation with Radial Drainage

Geotechnical Testing Journal, Vol. 26, No. 4 Paper ID GTJ10173_264 Available online at: www.astm.org Tian Ho Seah 1 and Teerawut Juirnarongrit 2 Constant Rate of Strain Consolidation with Radial Drainage

More information

3 Numerical Methods for Convection

3 Numerical Methods for Convection The topic of computational fluid dynamics (CFD) could easily occupy an entire semester indeed, we have such courses in our catalog. The objective here is to eamine some

More information

Laboratory 3 Hydraulic Conductivity of a Porous Media

Laboratory Hydraulic Conductivity of a Porous Media INTRODUCTION Permeability refers to the propensity of a material to allow a gas or fluid to move through its pores or interstices. A material is permeable

More information

Finite Element Analysis of Elastic Settlement of Spreadfootings Founded in Soil

Finite Element Analysis of Elastic Settlement of Spreadfootings Founded in Soil Jae H. Chung, Ph.D. Bid Bridge Software Institute t University of Florida, Gainesville, FL, USA Content 1. Background 2.

More information

DETERMINING POTENTIAL VERTICAL RISE

Test Procedure for DETERMINING POTENTIAL VERTICAL RISE TxDOT Designation: Tex-124-E Effective Dates: August 1999 February 2014. 1. SCOPE 1.1 This procedure determines the potential vertical rise (PVR)

More information

Finite Element Method (FEM) Introduction

Engineering manual No. 20 Updated: 06/2016 Finite Element Method (FEM) Introduction The objective of this engineering manual is to explain the basic terms of a particular field of problems and the practical

More information

c. Borehole Shear Test (BST): BST is performed according to the instructions published by Handy Geotechnical Instruments, Inc.

Design Manual Chapter 6 - Geotechnical 6B - Subsurface Exploration Program 6B-2 Testing A. General Information Several testing methods can be used to measure soil engineering properties. The advantages,

More information

HWR 431 / 531 HYDROGEOLOGY LAB SECTION LABORATORY 2 DARCY S LAW

HWR 431 / 531 HYDROGEOLOGY LAB SECTION LABORATORY 2 DARCY S LAW Introduction In 1856, Henry Darcy, a French hydraulic engineer, published a report in which he described a series of experiments he had performed

More information

CHAPTER 9 FEM MODELING OF SOIL-SHEET PILE WALL INTERACTION

391 CHAPTER 9 FEM MODELING OF SOIL-SHEET PILE WALL INTERACTION 9.1 OVERVIEW OF FE SOIL-STRUCTURE INTERACTION Clough and Denby (1969) introduced Finite Element analysis into the soil-structure interaction

More information

Copyright 2011 Casa Software Ltd. www.casaxps.com. Centre of Mass

Centre of Mass A central theme in mathematical modelling is that of reducing complex problems to simpler, and hopefully, equivalent problems for which mathematical analysis is possible. The concept of

More information

The First Law of Thermodynamics: Closed Systems. Heat Transfer

The First Law of Thermodynamics: Closed Systems The first law of thermodynamics can be simply stated as follows: during an interaction between a system and its surroundings, the amount of energy gained

More information

Fluid Mechanics Prof. T. I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay. Lecture No. # 36 Pipe Flow Systems

Fluid Mechanics Prof. T. I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay Lecture No. # 36 Pipe Flow Systems Welcome back to the video course on Fluid Mechanics. In today

More information

RESILIENT MODULUS: WHAT, WHY, AND HOW?

RESILIENT MODULUS: WHAT, WHY, AND HOW? Shane Buchanan, Vulcan Materials Company, 08/1/07 Resilient modulus (M r ) continues to be a key design parameter for pavement systems. Accurate knowledge of the

More information

Chapter 3 Geotechnical Conditions

PART II ACTIONS AND MATERIAL STRENGTH REQUIREMENTS, CHAPTER 3 GEOTECHNICAL CONDITIONS Chapter 3 Geotechnical Conditions Public Notice Geotechnical Conditions Article 13 Geotechnical conditions shall be

More information

GEO-SLOPE International Ltd, Calgary, Alberta, Canada Seepage Exit Gradients

1 Introduction Seepage Exit Gradients Most Soil Mechanics text books present and discuss the concept of seepage exit gradients and state that the exit gradients should not be greater than 1.0 Applying

More information

Soil Mechanics Prof. B.V.S. Viswanadham Department of Civil Engineering Indian Institute of Technology, Bombay Lecture - 2

Soil Mechanics Prof. B.V.S. Viswanadham Department of Civil Engineering Indian Institute of Technology, Bombay Lecture - 2 In the previous lecture, we have seen the origin of soils and different types

More information

Neural network for constitutive modelling in finite element analysis

Neural network for constitutive modelling in finite element analysis A. A. Javadi, T. P. Tan Department of Engineering, School of Engineering and Computer Science University of Exeter, Exeter, Devon, EX4

More information

Chapter 8 Steady Incompressible Flow in Pressure Conduits

Chapter 8 Steady Incompressible Flow in Pressure Conduits Outline 8.1 Laminar Flow and turbulent flow Reynolds Experiment 8.2 Reynolds number 8.3 Hydraulic Radius 8.4 Friction Head Loss in Conduits of

More information

Shallow Foundation Analysis on Heterogeneous Soil Formation in the Niger Delta of Nigeria

Shallow Foundation Analysis on Heterogeneous Soil Formation in the Niger Delta of Nigeria Akpila, S.B Department of Civil Engineering Rivers State University of Science & Technology, PortHarcourt. Eluozo,

More information

Foundation Engineering Prof. Mahendra Singh Department of Civil Engineering Indian Institute of Technology, Roorkee

Foundation Engineering Prof. Mahendra Singh Department of Civil Engineering Indian Institute of Technology, Roorkee Module - 03 Lecture - 09 Stability of Slopes Welcome back to the classes of on this Stability

More information

Module 2 Lecture 7 Permeability and Seepage -3 Topics

Module 2 Lecture 7 Permeability and Seepage -3 Topics 1.1.6 Determination of Coefficient of Permeability in the Field Pumping from wells Packer test 1.1.7 Theoretical Solution for Coefficient of Permeability

More information

Stress and Deformation Analysis. Representing Stresses on a Stress Element. Representing Stresses on a Stress Element con t

Stress and Deformation Analysis Material in this lecture was taken from chapter 3 of Representing Stresses on a Stress Element One main goals of stress analysis is to determine the point within a load-carrying

More information

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 01

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 01 Welcome to the series of lectures, on finite element analysis. Before I start,

More information

FLUID FORCES ON CURVED SURFACES; BUOYANCY

Rocscience Settle 3d

FLUID FORCES ON CURVED SURFCES; BUOYNCY The principles applicable to analysis of pressure-induced forces on planar surfaces are directly applicable to curved surfaces. s before, the total force on the

More information

Laterally Loaded Piles

Laterally Loaded Piles 1 Soil Response Modelled by p-y Curves In order to properly analyze a laterally loaded pile foundation in soil/rock, a nonlinear relationship needs to be applied that provides soil

More information

Worked Examples of mathematics used in Civil Engineering

Worked Examples of mathematics used in Civil Engineering Worked Example 1: Stage 1 Engineering Surveying (CIV_1010) Tutorial - Transition curves and vertical curves. Worked Example 1 draws from CCEA Advanced

More information

SOIL IMPROVEMENT BY PRECOMPRESSION AT A TANK FARM SITE IN CENTRAL JAVA, INDONESIA

SOIL IMPROVEMENT BY PRECOMPRESSION AT A TANK FARM SITE IN CENTRAL JAVA, INDONESIA Kul Bhushan, President, Group Delta Consultants, Inc., Aliso Viejo, CA, USA Carlos V. Amante, Project Engineer, Group Delta

More information

CBE 6333, R. Levicky 1. Potential Flow

CBE 6333, R. Levicky Part I. Theoretical Background. Potential Flow Potential Flow. Potential flow is irrotational flow. Irrotational flows are often characterized by negligible viscosity effects. Viscous

More information

CE-632 Foundation Analysis and Design

CE-63 Foundation Analysis and Design Shallow Foundations 1 SUMMARY of Terminology Gross Loading Intensity Total pressure at the level of foundation including the weight of superstructure, foundation, and

More information

Soil Mechanics. Outline. Shear Strength of Soils. Shear Failure Soil Strength. Laboratory Shear Strength Test. Stress Path Pore Pressure Parameters

Soil Mechanics Shear Strength of Soils Chih-Ping Lin National Chiao Tung Univ. cplin@mail.nctu.edu.tw 1 Outline Shear Failure Soil Strength Mohr-Coulomb Failure Criterion Laboratory Shear Strength Test

More information

Differential Relations for Fluid Flow. Acceleration field of a fluid. The differential equation of mass conservation

Differential Relations for Fluid Flow In this approach, we apply our four basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of

More information

GEOTECHNICAL ENGINEERING FORMULAS. A handy reference for use in geotechnical analysis and design

GEOTECHNICAL ENGINEERING FORMULAS A handy reference for use in geotechnical analysis and design TABLE OF CONTENTS Page 1. SOIL CLASSIFICATION..3 1.1 USCS: Unified Soil Classification System..3 1.1.1

More information

Types of Elements

chapter : Modeling and Simulation 439 142 20 600 Then from the first equation, P 1 = 140(0.0714) = 9.996 kn. 280 = MPa =, psi The structure pushes on the wall with a force of 9.996 kn. (Note: we could

More information

Program COLANY Stone Columns Settlement Analysis. User Manual

User Manual 1 CONTENTS SYNOPSIS 3 1. INTRODUCTION 4 2. PROBLEM DEFINITION 4 2.1 Material Properties 2.2 Dimensions 2.3 Units 6 7 7 3. EXAMPLE PROBLEM 8 3.1 Description 3.2 Hand Calculation 8 8 4. COLANY

More information

Benchmarking Multi-Dimensional Large Strain Consolidation Analyses D. Priestley 1, M.D. Fredlund 2 and D. van Zyl 3

Benchmarking Multi-Dimensional Large Strain Consolidation Analyses D. Priestley 1, M.D. Fredlund 2 and D. van Zyl 3 1,3 University of British Columbia 6350 Stores Road Vancouver, BC, V6T 1Z4 2 SoilVision

More information

Soil Mechanics SOIL STRENGTH page 1

Soil Mechanics SOIL STRENGTH page 1 Contents of this chapter : CHAPITRE 6. SOIL STRENGTH..1 6.1 PRINCIPAL PLANES AND PRINCIPAL STRESSES..1 6.2 MOHR CIRCLE..1 6.2.1 POLE METHOD OF FINDING STRESSES ON

More information

Math 018 Review Sheet v.3

Math 018 Review Sheet v.3 Tyrone Crisp Spring 007 1.1 - Slopes and Equations of Lines Slopes: Find slopes of lines using the slope formula m y y 1 x x 1. Positive slope the line slopes up to the right.

More information

Integration of a fin experiment into the undergraduate heat transfer laboratory

Integration of a fin experiment into the undergraduate heat transfer laboratory H. I. Abu-Mulaweh Mechanical Engineering Department, Purdue University at Fort Wayne, Fort Wayne, IN 46805, USA E-mail: mulaweh@engr.ipfw.edu

More information

New building plans on the surface above the existing underground Willemspoortunnel in Rotterdam

New building plans on the surface above the existing underground Willemspoortunnel in Rotterdam Ragavan Appiah, Aronsohn raadgevende ingenieurs bv, Rotterdam, The Netherlands The development of new projects

More information

METHODS FOR ACHIEVEMENT UNIFORM STRESSES DISTRIBUTION UNDER THE FOUNDATION

International Journal of Civil Engineering and Technology (IJCIET) Volume 7, Issue 2, March-April 2016, pp. 45-66, Article ID: IJCIET_07_02_004 Available online at http://www.iaeme.com/ijciet/issues.asp?jtype=ijciet&vtype=7&itype=2

More information

10-1 10. CONSOLIDATION

10-1 10. CONSOLIDATION 10.1 INFLUENCE OF DRAINAGE ON RATE OF SETTLEMENT When a saturated stratum of sandy soil is subjected to a stress increase, such as that caused by the erection of a building on the

More information

ESTIMATION OF UNDRAINED SETTLEMENT OF SHALLOW FOUNDATIONS ON LONDON CLAY

International Conference on Structural and Foundation Failures August 2-4, 2004, Singapore ESTIMATION OF UNDRAINED SETTLEMENT OF SHALLOW FOUNDATIONS ON LONDON CLAY A. S. Osman, H.C. Yeow and M.D. Bolton

More information

The interpretation of the flat dilatometer tests in Ariake clay, southwest Japan

Mem. Fac. Sci. Eng. Shimane Univ. Series A 37, pp. 1-6 2003 The interpretation of the flat dilatometer tests in Ariake clay, southwest Japan Takeshi KAMEI 1) and Masanori TANAKA 2) 1) Department of Geoscience,

More information

Eurocode 7 - Geotechnical design - Part 2 Ground investigation and testing

Brussels, 18-20 February 2008 Dissemination of information workshop 1 Eurocode 7 - Geotechnical design - Part 2 Ground investigation and testing Dr.-Ing. Bernd Schuppener, Federal Waterways Engineering

More information

Chapter (1) Fluids and their Properties

Chapter (1) Fluids and their Properties Fluids (Liquids or gases) which a substance deforms continuously, or flows, when subjected to shearing forces. If a fluid is at rest, there are no shearing forces

More information

SOME OBSERVATIONS ON THE OEDOMETRIC CONSOLIDATION STRAIN RATE BEHAVIORS OF SATURATED CLAY

Journal of GeoEngineering, Vol. 5, No. 1, pp. 1-7, Feng: April Some 2010 Observations on the Oedometric Consolidation Strain Rate Behaviors of Saturated Clay 1 SOME OBSERVATIONS ON THE OEDOMETRIC CONSOLIDATION

More information

DESIGN OF PILES AND PILE GROUPS CONSIDERING CAPACITY, SETTLEMENT, AND NEGATIVE SKIN FRICTION

DESIGN OF PILES AND PILE GROUPS CONSIDERING CAPACITY, SETTLEMENT, AND NEGATIVE SKIN FRICTION Introduction Bengt H. Fellenius, Dr.Tech., P.Eng. Background Notes for Demo Example for UniPile at www.unisoftltd.com

More information

4 The strength of the lithosphere

4 The strength of the lithosphere Constitutive equations specify the relations between stress and strain or the time derivatives of strain and thus can be used in conjunction with the equations of motion

More information

4. Subgrade : Foundation Considerations

4. Subgrade : Foundation Considerations 1. Introduction The performance of any structure is strongly dependent on the foundations. Road pavements are no exceptions and their foundations need to be fundamentally

More information

Module 2. The Science of Surface and Ground Water. Version 2 CE IIT, Kharagpur

Module 2 The Science of Surface and Ground Water Lesson 2 Runoff and Infiltration Instructional Objectives At the end of this lesson, the student shall be able to learn: 1. The importance of runoff and

More information

Cement Stabilization of Organic Subgrades for Pavement Preservation

Tenth LACCEI Latin American and Caribbean Conference (LACCEI 2012), Megaprojects: Building Infrastructure by fostering engineering collaboration, efficient and effective integration and innovative planning,

More information

CHAPTER: 6 FLOW OF WATER THROUGH SOILS

CHAPTER: 6 FLOW OF WATER THROUGH SOILS CONTENTS: Introduction, hydraulic head and water flow, Darcy s equation, laboratory determination of coefficient of permeability, field determination of coefficient

More information

A Theoretical Model for Mutual Interaction between Coaxial Cylindrical Coils Lukas Heinzle

A Theoretical Model for Mutual Interaction between Coaxial Cylindrical Coils Lukas Heinzle Page 1 of 15 Abstract: The wireless power transfer link between two coils is determined by the properties of the

More information

A Strategy for Teaching Finite Element Analysis to Undergraduate Students

A Strategy for Teaching Finite Element Analysis to Undergraduate Students Gordon Smyrell, School of Computing and Mathematics, University of Teesside The analytical power and design flexibility offered

More information

Gravity waves on water

Phys374, Spring 2006, Prof. Ted Jacobson Department of Physics, University of Maryland Gravity waves on water Waves on the surface of water can arise from the restoring force of gravity or of surface tension,

More information

Senior Secondary Australian Curriculum

Senior Secondary Australian Curriculum Mathematical Methods Glossary Unit 1 Functions and graphs Asymptote A line is an asymptote to a curve if the distance between the line and the curve approaches zero

More information

CONSTANT HEAD AND FALLING HEAD PERMEABILITY TEST

CONSTANT HEAD AND FALLING HEAD PERMEABILITY TEST 1 Permeability is a measure of the ease in which water can flow through a soil volume. It is one of the most important geotechnical parameters. However,

More information

DETERMINATION OF SOIL STRENGTH CHARACTERISTICS PERFORMING THE PLATE BEARING TEST

III Międzynarodowa Konferencja Naukowo-Techniczna Nowoczesne technologie w budownictwie drogowym Poznań, 8 9 września 005 3rd International Conference Modern Technologies in Highway Engineering Poznań,

More information

Fluid Mechanics. Fluid Statics [3-1] Dr. Mohammad N. Almasri. [3] Fall 2010 Fluid Mechanics Dr. Mohammad N. Almasri [3-1] Fluid Statics

1 Fluid Mechanics Fluid Statics [3-1] Dr. Mohammad N. Almasri Fluid Pressure Fluid pressure is the normal force exerted by the fluid per unit area at some location within the fluid Fluid pressure has the

More information

Soil Mechanics Prof. B.V.S. Viswanathan Department of Civil Engineering Indian Institute of Technology, Bombay Lecture 19 Effective Stress-I I I

Soil Mechanics Prof. B.V.S. Viswanathan Department of Civil Engineering Indian Institute of Technology, Bombay Lecture 19 Effective Stress-I I I Welcome to the third lecture on effective stress. In the

More information

Module 7 (Lecture 24 to 28) RETAINING WALLS

Module 7 (Lecture 24 to 28) RETAINING WALLS Topics 24.1 INTRODUCTION 24.2 GRAVITY AND CANTILEVER WALLS 24.3 PROPORTIONING RETAINING WALLS 24.4 APPLICATION OF LATERAL EARTH PRESSURE THEORIES TO DESIGN 24.5

More information

Centre, Bandar Tasik Selatan, Kuala Lumpur, Malaysia.

Methodology for Design of Piled Raft for Five-Storeys Buildings on Very Soft Clay Tan, Y.C. 1, Cheah, S.W. 1 and Taha, M.R. 2 Abstract Conventional piled foundation is usually designed for buildings to

More information

VARIOUS ASPECTS OF EXPANSIVE SOILS RELEVENT TO GEOENGINEERING PRACTICE. Simple Correlations Between Soil Plasticity and Expansion Potential

GE 441 Advanced Engineering Geology & Geotechnics Spring 2004 VARIOUS ASPECTS OF EXPANSIVE SOILS RELEVENT TO GEOENGINEERING PRACTICE Simple Correlations Between Soil Plasticity and Expansion Potential

More information

International Society for Helical Foundations (ISHF)

QUICK DESIGN GUIDE For Screw-Piles and Helical Anchors in Soils Ver. 1.0 Prepared by Dr. Alan J. Lutenegger, P.E., F. ASCE for International Society for Helical Foundations (ISHF) Copyright 2015 (ISHF)

More information

Probabilistic Analysis

Probabilistic Analysis Tutorial 8-1 Probabilistic Analysis This tutorial will familiarize the user with the basic probabilistic analysis capabilities of Slide. It will demonstrate how quickly and easily

More information

MECHANICS OF SOLIDS - BEAMS TUTORIAL 3 THE DEFLECTION OF BEAMS

MECHANICS OF SOLIDS - BEAMS TUTORIAL THE DEECTION OF BEAMS This is the third tutorial on the bending of beams. You should judge your progress by completing the self assessment exercises. On completion

More information

Calculating Areas Section 6.1

A B I L E N E C H R I S T I A N U N I V E R S I T Y Department of Mathematics Calculating Areas Section 6.1 Dr. John Ehrke Department of Mathematics Fall 2012 Measuring Area By Slicing We first defined

More information

1 Drained and Undrained Triaxial Tests on Sand

1 Drained and Undrained Triaxial Tests on Sand The behavior of sands in triaxial state of stress is the focus of attention in this verification. Sands, while shearing in triaxial test, show a hardening

More information