A Stabilization Procedure for Soil-water Coupled Problems Using the Mesh-free Method
|
|
- Andrea Fitzgerald
- 5 years ago
- Views:
Transcription
1 The 12 th International Conference of International Association for Computer Methods and Advances in Geomechanics (IACMAG) 1-6 October, 28 Goa, India A Stabilization Procedure for Soil-ater Coupled Problems Using the Mesh-free Method T. Shibata Dept. of Civil and Environmental Engineering, Matsue National College of Technology, Shimane, Japan A. Murakami Graduate School of Environmental Science, Okayama University, Okayama, Japan Keyords: Mesh-free method, soil-ater coupled problem, stabilization procedure ABSTRACT: The development of stability problems related to classical mixed methods has recently been observed. In this study, a soil-ater coupled boundary-value problem, one type of stability problem, is presented using the Element-free Galerkin Method (EFG Method). In this soil-ater coupled problem, anomalous behavior appears in the pressure field unless stabilization techniques are used. The remedy to such numerical instability has generally been to adopt a higher interpolation order for the displacements than for the pore pressure. As an alternative, hoever, an added stabilization term is incorporated into the equilibrium equation. The advantages of this stabilization procedure are as follos: (1) The interpolation order for the pore pressure is the same as that for the displacements. Therefore, the interpolation functions in the pore pressure field do not reduce the accuracy of the numerical results. (2) The stabilization term consists of first derivatives. The first derivatives of the interpolation functions for the EFG Method are smooth, and therefore, the solutions for pore pressure are accurate.in order to validate the above stabilization technique, some numerical results are given. It can be seen from the results that a good convergence is obtained ith this stabilization term. 1 Introduction The development of numerical computation technologies has enabled a variety of engineering problems to be solved and has brought about remarkable progress in recent decades. Among the related findings, meshless and/or mesh-free methods in particular have been applied to some problems for hich the usual finite element method is ineffective in dealing ith significant mesh distortion brought about by large deformations, crack groth, and moving discontinuities. Various meshless and/or mesh-free methods have been used for geotechnical problems, instead of the finite element method, to overcome the above-mentioned difficulties. Consolidation phenomena have been analyzed by means of EFGM (Modaressi et al., 1998)(Nogami et al., 24)(Murakami et al., 25)(Wang et al., 26), the point/radial point interpolation method (PIM/RPIM) (Wang et al., 21)(Wang et al., 22), the local RPIM (Wang et al., 25), RKPM (Chen et al., 21)(Zhang et al., 25), and the natural neighbor method (Cai et al., 25), the transient response of saturated soil has been dealt ith under cyclic loading by means of EFGM (Karim et al., 22)(Sato et al., 26), ave-induced seabed response and instability have been examined by EFGM (Wang et al., 27) and RPIM (Wang et al., 24), slip lines have been modeled by geological materials using EFGM (Rabczuk et al., 26), and a Bayesian inverse analysis has been carried out in conjunction ith the meshless local Petrov-Galerkin method (Sheu, 26). Hoever, unless certain requirements are met in dealing ith soil-ater coupled problems for the finite element computation, based on the coupled formulation becoming ill-conditioned, numerical instabilities ill occur (Chapelle et al., 1993). In order to overcome these eaknesses, several strategies have been proposed (Pastor et al., 1999). For example, as a necessary condition for stability, the interpolation degree of the displacement field must be higher than that of the pore pressure field. An alternative means of stabilization as also proposed based on the Simo-Rifai enhanced strain method hich even allos an equal order of interpolation degree for both variables. Hoever, these strategies are not directly applicable to meshless/mesh-free methods, because all the nodal points simultaneously have the same degree of freedom for both the displacement field and the pore pressure field, and no information beteen the element and the nodes can be utilized. The purpose of this paper is to present a stabilization methodology for the mesh-free analysis of soil-ater coupled problems by incorporating the stabilizing term into the eak form. The folloing sections deal ith descriptions of the formulation, including the stabilization term. In Section 3, to applications of the strategy to 64
2 soil-ater coupled problems are analyzed, one being the saturated soil column test appearing in Mira et al.(23), to demonstrate the effectiveness of the strategy, and the other being the foundation behavior under displacement-controlled condition, for hich the feasibility of theanalysis ill be thoroughly discussed hile Table 1. Order of the interpolation function in FEM Number Pressure Order Displacement :Pressure :Displacement Figure 1. Rearrangement and allocation of background cells to cover the changing domain comparing the finite element solutions. Conclusions ill follo in Section 4. 2 The Stabilization term The governing equations for soil-ater coupled problems ith boundary conditions and initial conditions are given as follos: div S & δ v dv + ρ v v δ v d S = (1) V t Γ u ( ) V V Γ Γ ( tr ) δhdv + v gradδhdv qδhds β ( p p ) D δhds = (2) here S & t is the nominal stress rate, v is the velocity, v is the boundary value of the velocity, D is the stretching, h is the total head, v is the velocity of the pore ater, q is the discharge per unit area ith units of length per time, β is Lagrange multiplier, p is the pore ater pressure, and p is the boundary value of the pore ater pressure. As previously mentioned, mixed displacement-pressure formulations (e.g., finite element methods) produce locking phenomena in the pressure field unless stabilization techniques are used. The remedy for such numerical instability has generally been to adopt a higher interpolation order for the displacements than for the pore pressure. As an alternative, hoever, an added stabilization term is incorporated. It is shon in this study that the instability can be eliminated by adding the stabilization term to the continuity of pore ater (Equation (2)) hich consists of the square of the pore ater pressure of the first derivatives. q h + δ ( tr D) δhdv + v gradδh dv qδh ds β ( p p ) V q h 2 α ( p ) dv = V Γ Γ V δh ds (3) here α is the stability parameter, p is the differentiate p ith respect to z, and z is the vertical axis. In this study, the soil column test created by Zienkieicz et al. is performed in order to examine the numerical stability of this procedure, and the values for the pore ater pressure are illustrated along the vertical axis of the soil column. In order to smooth the values for the pore ater pressure along the vertical axis ith the continuity condition of the pore ater, the stabilization term is differentiated ith respect to z. The to advantages of this stabilization procedure are as follos: (1) The interpolation order for the pore pressure is the same as that for the displacements, namely, a loer interpolation order is not adopted for the pore 65
3 pressure. Therefore, the interpolation functions in the pore pressure field do not reduce the accuracy of the numerical results. (2) Table 1 summarizes the order of the interpolation functions for the pressure field and for the displacement field in FEM. The first derivatives of the interpolation functions in the pore pressure field are the zero order or the first order, as shon Table 1. Therefore, accuracy in the numerical results cannot be obtained. With EFGM, hoever, the interpolation functions are derived by the MLS approximant and a linearbased polynomial is used, in other ords, the resultant interpolation functions are smooth. The eak forms are integrated using the MLS interpolation functions in space and employing the explicit time scheme in time. In the manipulation of the stiffness matrix, a numerical integration is performed on the Table 2. Material parameters Compression index λ.11 Selling index κ.4 Critical state parameter M 1.42 Poisson s ratio ν.333 Initial void ratio e.83 Initial volume ratio v = 1+ e 1.83 Initial consolidation stress (kpa) p 294 z = 3 p = u = 3 m z = u x = u = z 1 m q p z p x z = z u = x = x Figure 2. Saturated column test 3 Depth ( m ) 2 1 k= ( cm/sec ) Pore ater pressure ( Pa ) [ 1 5 ] ithout stabilization term α =.1 α =.1 α =.1 Figure 3. Effect of the dimensionless parameter background cells. The background cells, hich intersect or make contact ith the boundary, are divided into four finer cells. The cells are then rearranged during the computation, according to changes in the domain, as shon in Figure 1. 3 Numerical Examples 3.1 The saturated column soil test For soil-ater coupled problems, numerical instabilities are often encountered at the initial stage under undrained conditions unless stabilization techniques are performed. In this chapter, the numerical stability of an EFG computation is examined using the proposed stabilization term described in the last chapter. Let s consider the 1D problem analyzed by Zienkieicz et al.(1986) in hich the saturated soil column test in Figure 2 uses the material parameters listed in Table 2. Here, a sand permeability of cm/s is employed. A model discretized by 62 nodal points is adopted. Linear-based polynomials are applied for the interpolation functions of the EFG Method. The functions have the same order for both the displacements and the pore pressure. The eight function is a quartic spline type of eight function and its radius of support is 1.. Here, 5 5 Gaussian points are used. The background cells are 1.m ithin the domain and.5m outside of the domain. The scale factor, hich is defined as the magnification of the support diameter to the side length of the square background cell, is 1.5. This study employs penalty methods to apply the boundary conditions, and the value of the penalty factor is In order to solve the stiffness equation, the forard difference is 66
4 adopted. Figure 3 shos the effect of the proposed stabilization term and dimensionless parameter α on the numerical profile of the pore pressure just after loading under undrained conditions, for hich a time difference of Δ t =.1 and a value of permeability of k = cm/s are adopted. There are 62 nodal points and the eight function is a quartic spline. From Figure 3, it can be concluded that the numerical solution has been improved in the case here the value of the dimensionless parameter is.1. In subsequent examinations, a permeability of cm/s and a time difference of.1 are adopted. 2. Here, the numerical instability is caused from the large difference of the value in the stiffness matrix in beteen displacements and pore ater pressures. Generally, the values of stiffness matrix in pore ater pressure are smaller than in displacements. Hoever, if the loer order interpolation function in pore ater pressures is adopted, the values of stiffness matrix in pore ater pressures become large, so that the difference of the values in stiffness matrix becomes small. Similarly, if the stabilization term is considered, the values of stiffness matrix Loading plate p = λ =.11 κ =.4 M = 1.42 ν =.333 e =.83 v = 1+ e =1.83 p = 294 (kpa) 5. m 4.m 2.m 4.m 1.m (a) Geometry and boundary conditions (b) Initial collocation of the nodal points Figure 4. Sketch of the problem [kpa] (a) Contours of the pore pressure (b) Contours of the normalized strain (c) Displacement distributions (d) Collocation of the nodal points Figure 5. Numerical results ithout stabilization term 67
5 in pore ater pressure also become large. Therefore, the stabilization term quells the anomalous pore pressure behavior because of the small difference. 3.2 Foundation Problem Subjected to Strip Loading In order to examine the numerical availability of the stabilization procedure to EFGM, a 2D soil-ater coupled problem is solved in relation to a soft soil foundation. In this problem, the p& in Equation (1) is the differentiate p ith respect to x and z, ith x being the horizontal axis and z being the vertical axis. The geometry and the boundary conditions are given in Figure 4(a), hile the initial collocation of the nodal points is shon in Figure 4(b). Vertical displacements are applied at the top face to simulate loading on the foundation, hile the [kpa] (a) Contours of the pore pressure (b) Contours of the normalized strain (c) Displacement distributions (d) Collocation of the nodal points Figure 6. Numerical results ith stabilization term Applied footing stress [MPa].2.1 Prandtl's solution EFG solution Settlement [m] Figure 7. Comparison EFG solution ith Prandtl s solution transverse direction is restrained. Thus, the boundary conditions under the loading surface are given by the settlements, and the other locations on the upper boundary are stress free. The bottom boundary is fixed in all directions, hile the side boundaries are fixed in only the horizontal direction, such that vertical displacement is alloed. Hydrostatic pressure is used here for the initial conditions in the pore pressure field. The model is generated ith 1,326 nodal points, in other ords, 26 vertical nodal points and 51 horizontal nodal points. The size of background cells is.2m, and the value of the dimensionless parameter α is same as the previous analysis. Figures 5 and 6 present the numerical results ithout and ith the stabilization term, respectively, in hich (a), (b), (c) and (d) sho the contours of the pore ater pressure, the contours of the normalized strain, the displacement distributions and collocation of the nodal points, respectively. In these results, the settlements 68
6 under the loading surface are.4m in Figures 5 and.4m in Figures 6, respectively. The normalized strain measure is given as t here ε = Ddt (Yatomi et al., 1989). T ( ) ε = tr εε (2) Spurious oscillations arise in the pore pressure field, as can be seen in Figure 5(a). In contrast to the results using the stabilization procedure, high values are obtained for the pore pressure and the normalized strain in Figures 5(a) and (b), respectively. In particular, anomalous behavior appears on left side in these figures. Moreover, the directions of the displacement vectors are disorderly in Figure 5(c) because of the interaction beteen the pore pressure field and the displacement field. If stabilization techniques are used, hoever, no oscillations appear in the solution, as observed in Figure 6. Figure 6(a) accurately displays the rise in the pore ater pressure belo the loading surface. In Figure 6(b), the prominently localized zones of the normalized strain occur just beneath the edge of the loading surface. The deformed pattern in Figure 6(c) is very similar to the classical slip line solution obtained by Prandtl. The shear bands are recognized as the localized deformation. Figure 7 compares the EFG solution ith Prandtl s solution. The Prandtl s solution q f is expressed as q f =5.14 c u (3) here c u is the undrained shear strength. From this figure, it is revealed that the EFG solution approaches Prandtl s solution. This result shos that the EFG method ith stabilizing procedure is very efficient for solving problems of computational geomechanics. 4 Conclusion In this paper, e have proposed a stabilization method for soil-ater coupled problems using the Element-free Galerkin Method. A stabilization term has been presented by the addition of a continuity condition. The proposed stabilization procedure has the folloing to characteristics: (1) The interpolation order for the pore pressure field is the same as that for the displacements, in other ords, a loer interpolation order for the pore pressure is not adopted. (2) The stabilization term consists of first derivatives in hich the interpolation functions are smooth because of the MLS approximation. The saturated column test and the foundation loading problem have been solved using the stabilization procedure. Numerical examples have shon that the stabilization method can indeed quell the anomalous pore pressure behavior. 5 References Modaressi H, Aubert P., 1998, Element-free Galerkin method for deforming multiphase porous media, International Journal for Numerical Methods Engineering, 42, Nogami T, Wang W, Wang JG., 24, Numerical method for consolidation analysis of lumpy clay fillings ith meshless method, Soils and Foundations, 44(1), Murakami A, Setsuyasu T, Arimoto S., 25, Mesh-free method for soil-ater coupled problem ithin finite strain and its numerical validity, Soils and Foundations, 45(2). Wang ZL, Li YC., 26, Analysis of factors influencing the solution of the consolidation problem by using an element-free Galerkin method, Computers & Geosciences, 32, Wang JG, Liu GR, Wu YG., 21, A point interpolation method for simulating dissipation process of consolidation, Computer Methods in Applied Mechanics and Engineering, 19, Wang JG, Liu GR, Lin P., 22, Numerical analysis of Biot s consolidation process by radial point interpolation method, International Journal of Solids and Structures, 39, Wang JG, Yan L, Liu GR., 25, A local radial point interpolation method for dissipation process of excess pore ater pressure, International Journal of Numerical Methods for Heat & Fluid Flo, 15(6), Chen JS, Wu CT, Chi L, Huck F., 21, A meshfree method for geotechnical materials, Journal of Engineering Mechanics, 127, Zhang JF, Zhang WP, Zheng Y., 25, A meshfree method and its applications to elasto-plastic problems, Journal of Zhejiang University SCIENCE, 6A(2), Cai YC, Zhu HH, Xia CC., 25, Meshless natural neighbor method and its application to complex geotechnical engineering, Yanshilixue Yu Gongcheng Xuebao/Chinese Journal of Rock Mechanics and Engineering, 24(11), Karim MR, Nogami T, Wang JG., 22, Analysis of transient response of saturated porous elastic soil under cyclic loading using element-free Galerkin method, International Journal of Solids and Structures, 39, Sato T, Matsumaru T., 26, Numerical simulation of liquefaction and flo process using mesh free method, Journal of Geotechnical Engineering, 813, (in Japanese). 69
7 Wang JG, Karim MR, Lin PZ., 27, Analysis of seabed instability using element free Galerkin method, Ocean Engineering, 34(2), Wang JG, Nogami T, Dasari GR, Lin PZ., 24, A eak coupling algorithm for seabed-ave interaction analysis, Computer Methods in Applied Mechanics and Engineering, 193, Rabczuk T, Areias PMA., 26, A ne approach for modeling slip lines in geological materials ith cohesive models, International Journal for Numerical and Analytical Methods in Geomechanics, 3, Sheu GY., 26, Direct back analysis by the meshless local Petrov-Galerkin method and Bayesian statistics, International Journal for Numerical and Analytical Methods in Geomechanics, 3(8), Chapelle D, Bathe KJ., 1993, The inf-sup test, Computers & Structures, 47, Pastor M, Li T, Liu X, Zienkieicz OC., 1999, Stabilized lo-order finite elements for failure and localization problems in undrained soils and foundations, Computer Methods in Applied Mechanics and Engineering, 174, Mira P, Pastor M, Li T, Liu X., 23, A ne stabilized enhanced strain element ith equal order of interpolation for soil consolidation problems, Computer Methods in Applied Mechanics and Engineering, 192, Zienkieicz OC, Qu S, Taylor RL, Nakazaa S., 1986, The patch test for mixed formulations, International Journal for Numerical Methods Engineering, 23, Yatomi C, Yamashita A, Iizuka A, Sano I., 1989, Shear bands formation numerically simulated by a non-coaxial Cam-Clay model, Soils and Foundations, 29,
AN ELASTO-PLASTIC LIQUEFACTION ANALYSIS METHOD BASED ON FINITE DEFORMATION THEORY
The Fifteenth KKCNN Symposium on Civil Engineering December 19-20, 2002, Singapore AN ELASTO-PLASTIC LIQUEFACTION ANALYSIS METHOD BASED ON FINITE DEFORMATION THEORY F. Oka, T. Kodaka, R. Morimoto and N.
More informationTIME-DEPENDENT BEHAVIOR OF PILE UNDER LATERAL LOAD USING THE BOUNDING SURFACE MODEL
TIME-DEPENDENT BEHAVIOR OF PILE UNDER LATERAL LOAD USING THE BOUNDING SURFACE MODEL Qassun S. Mohammed Shafiqu and Maarib M. Ahmed Al-Sammaraey Department of Civil Engineering, Nahrain University, Iraq
More informationMonitoring of underground construction
Monitoring of underground construction Geotechnical Aspects of Underground Construction in Soft Ground Yoo, Park, Kim & Ban (Eds) 2014 Korean Geotechnical Society, Seoul, Korea, ISBN 978-1-138-02700-8
More informationAnalysis of Steady State Heat Conduction Problem Using EFGM
International Journal of Engineering and Management Research, Vol.-2, Issue-6, December 2012 ISSN No.: 2250-0758 Pages: 40-47 www.ijemr.net Analysis of Steady State Heat Conduction Problem Using EFGM Manpreet
More informationINFLUENCE OF SOIL NONLINEARITY AND LIQUEFACTION ON DYNAMIC RESPONSE OF PILE GROUPS
INFLUENCE OF SOIL NONLINEARITY AND LIQUEFACTION ON DYNAMIC RESPONSE OF PILE GROUPS Rajib Sarkar 1 and B.K. Maheshwari 2 1 Research Scholar, Dept. of Earthquake Engineering, IIT Roorkee, India, e-mail:
More informationFinite Deformation Analysis of Dynamic Behavior of Embankment on Liquefiable Sand Deposit Considering Pore Water Flow and Migration
6 th International Conference on Earthquake Geotechnical Engineering 1-4 November 215 Christchurch, New Zealand Finite Deformation Analysis of Dynamic Behavior of Embankment on Liquefiable Sand Deposit
More informationEFFECTS OF GROUND WATER ON SEISMIC RESPONSES OF BASIN
EFFECTS OF GROUND WATER ON SEISMIC RESPONSES OF BASIN Huei-Tsyr CHEN And Jern-Chern HO 2 SUMMARY It has long been recognized that the local soil and geology conditions may affect significantly the nature
More informationNUMERICAL ANALYSIS OF DAMAGE OF RIVER EMBANKMENT ON SOFT SOIL DEPOSIT DUE TO EARTHQUAKES WITH LONG DURATION TIME
Proceedings of the International Symposium on Engineering Lessons Learned from the 2011 Great East Japan Earthquake, March 1-4, 2012, Tokyo, Japan NUMERICAL ANALYSIS OF DAMAGE OF RIVER EMBANKMENT ON SOFT
More informationApplication of cyclic accumulation models for undrained and partially drained general boundary value problems
Application of cyclic accumulation models for undrained and partially drained general boundary value problems A. M. Page Risueño Yngres Dag 2014, May 15 th 2014 Introduction Cyclic loads in geotechnical
More informationInfluences of material dilatancy and pore water pressure on stability factor of shallow tunnels
Influences of material dilatancy and pore water pressure on stability factor of shallow tunnels YANG Xiao-li( ), HUANG Fu( ) School of Civil and Architectural Engineering, Central South University, Changsha
More informationNon-linear consolidation of soil with vertical and horizontal drainage under time-dependent loading
Uniersity of Wollongong Research Online Faculty of Engineering - Papers (Archie) Faculty of Engineering and Information Sciences Non-linear consolidation of soil ith ertical and horizontal drainage under
More informationNumerical model comparison on deformation behavior of a TSF embankment subjected to earthquake loading
Numerical model comparison on deformation behavior of a TSF embankment subjected to earthquake loading Jorge Castillo, Yong-Beom Lee Ausenco, USA Aurelian C. Trandafir Fugro GeoConsulting Inc., USA ABSTRACT
More informationFEM MODEL OF BIOT S EQUATION FREE FROM VOLUME LOCKING AND HOURGLASS INSTABILITY
he 14 th World Conference on Earthquake Engineering FEM MODEL OF BIO S EQUAION FREE FROM OLUME LOCKING AND HOURGLASS INSABILIY Y. Ohya 1 and N. Yoshida 2 1 PhD Student, Dept. of Civil and Environmental
More informationEFFECTS OF PLASTIC POTENTIAL ON THE HORIZONTAL STRESS IN ONE-DIMENSIONAL CONSOLIDATION
Journal of GeoEngineering, Vol. 11, No. 1, pp. 27-31, April Iinuma 216 et al.: Effects of Plastic Potential on the Horizontal Stress in One-Dimensional Consolidation 27 http://dx.doi.org/1.631/jog.216.11(1).3
More informationA MODEL FOR COUPLED DYNAMIC ELASTO-PLASTIC ANALYSIS OF SOILS
Journal of GeoEngineering, Vol. 7, No. 3, pp. 089-096, December 2012 Fattah et al.: A Model for Coupled Dynamic Elastic Plastic Analysis of Soils 89 A MODEL FOR COUPLED DYNAMIC ELASTO-PLASTIC ANALYSIS
More informationCONSOLIDATION BEHAVIOR OF PILES UNDER PURE LATERAL LOADINGS
VOL., NO., DECEMBER 8 ISSN 89-8 -8 Asian Research Publishing Network (ARPN). All rights reserved. CONSOLIDATION BEAVIOR OF PILES UNDER PURE LATERAL LOADINGS Qassun S. Mohammed Shafiqu Department of Civil
More informationPSEUDO ELASTIC ANALYSIS OF MATERIAL NON-LINEAR PROBLEMS USING ELEMENT FREE GALERKIN METHOD
Journal of the Chinese Institute of Engineers, Vol. 27, No. 4, pp. 505-516 (2004) 505 PSEUDO ELASTIC ANALYSIS OF MATERIAL NON-LINEAR PROBLEMS USING ELEMENT FREE GALERKIN METHOD Raju Sethuraman* and Cherku
More informationA Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials
Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical
More information1.8 Unconfined Compression Test
1-49 1.8 Unconfined Compression Test - It gives a quick and simple measurement of the undrained strength of cohesive, undisturbed soil specimens. 1) Testing method i) Trimming a sample. Length-diameter
More informationLateral Load Capacity of Rock Socketed Piers Using Finite Difference Approach
Journal of Civil Engineering esearch. 11; 1(1): 1-8 DOI: 1. 593/.ce.1111.1 ateral oad Capacity of ock Socketed Piers Using Finite Difference Approach Mayank Suhirid 1, Kiran B adhane 1, Mahendra Singh
More informationNonlinear Time-Dependent Soil Behavior due to Construction of Buried Structures
Journal of Earth Sciences and Geotechnical Engineering, vol. 4, no. 1, 214, 71-88 ISSN: 172-4 (print), 172- (online) Scienpress Ltd, 214 Nonlinear Time-Dependent Soil Behavior due to Construction of Buried
More informationAdvanced model for soft soils. Modified Cam-Clay (MCC)
Advanced model for soft soils. Modified Cam-Clay (MCC) c ZACE Services Ltd August 2011 1 / 62 2 / 62 MCC: Yield surface F (σ,p c ) = q 2 + M 2 c r 2 (θ) p (p p c ) = 0 Compression meridian Θ = +π/6 -σ
More informationLandslide FE Stability Analysis
Landslide FE Stability Analysis L. Kellezi Dept. of Geotechnical Engineering, GEO-Danish Geotechnical Institute, Denmark S. Allkja Altea & Geostudio 2000, Albania P. B. Hansen Dept. of Geotechnical Engineering,
More informationAnalytical and Numerical Investigations on the Vertical Seismic Site Response
Analytical and Numerical Investigations on the Vertical Seismic Site Response Bo Han, Lidija Zdravković, Stavroula Kontoe Department of Civil and Environmental Engineering, Imperial College, London SW7
More informationNumerical analysis of dissipation of pore water pressure using Natural Element Method
Bulletin of Environment, Pharmacology and Life Sciences Bull. Env.Pharmacol. Life Sci., Vol 4 [Spl issue 1] 2015: 371-375 2014 Academy for Environment and Life Sciences, India Online ISSN 2277-1808 Journal
More informationModelling Progressive Failure with MPM
Modelling Progressive Failure with MPM A. Yerro, E. Alonso & N. Pinyol Department of Geotechnical Engineering and Geosciences, UPC, Barcelona, Spain ABSTRACT: In this work, the progressive failure phenomenon
More informationInstability analysis and numerical simulation of the dissociation process of methane hydrate bearing soil
Computer Methods and Recent Advances in Geomechanics Oka, Murakami, Uzuoka & Kimoto (Eds.) 2015 Taylor & Francis Group, London, ISBN 978-1-138-00148-0 Instability analysis and numerical simulation of the
More informationAnalysis of fluid induced vibration of cryogenic pipes in consideration of the cooling effect
Journal of Mechanical Science and Technology (8) 375~385 Journal of Mechanical Science and Technology.springerlink.com/content/1738-494x DOI 1.17/s16-8-55-7 Analysis of fluid induced vibration of cryogenic
More informationPoint interpolation method based on local residual formulation using radial basis functions
Structural Engineering and Mechanics, Vol. 14, No. 6 (2002) 713-732 713 Point interpolation method based on local residual formulation using radial basis functions G.R. Liu, L. Yan, J.G. Wang and Y.T.
More informationGEO-SLOPE International Ltd, Calgary, Alberta, Canada Wick Drain
1 Introduction Wick Drain This example is about modeling the behavior of a wick drain. The primary purpose here is to illustrate how interface elements can conveniently be used to include the effects of
More informationPLAXIS. Scientific Manual
PLAXIS Scientific Manual 2016 Build 8122 TABLE OF CONTENTS TABLE OF CONTENTS 1 Introduction 5 2 Deformation theory 7 2.1 Basic equations of continuum deformation 7 2.2 Finite element discretisation 8 2.3
More informationAnalyzing Plane-plate Bending with EFGM
Journal of Mathematics Research March, 2009 Analyzing Plane-plate Bending with EFGM Yajing Zhang, Maohui Xia & Yun Zhai Department of Science, YanShan University Qinhuangdao 066004, China E-mail: zhang
More informationDYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION
October 1-17,, Beijing, China DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION Mohammad M. Ahmadi 1 and Mahdi Ehsani 1 Assistant Professor, Dept. of Civil Engineering, Geotechnical Group,
More informationTikrit University. College of Engineering Civil engineering Department CONSOILDATION. Soil Mechanics. 3 rd Class Lecture notes Up Copyrights 2016
Tikrit University CONSOILDATION College of Engineering Civil engineering Department Soil Mechanics 3 rd Class Lecture notes Up Copyrights 2016 Stresses at a point in a soil mass are divided into two main
More informationOPTIMAL SENSOR LOCATION FOR PARAMETER IDENTIFICATION IN SOFT CLAY
2 th International Conference on the Application of Computer Science and Mathematics in Architecture and Civil Engineering K. Gürlebeck and T. Lahmer (eds.) Weimar, Germany, 2-22 July 215 OPTIMAL SENSOR
More informationHypoplastic Cam-clay model
Hypoplastic model David Mašín Charles University in Prague corresponence address: David Mašíın Charles University in Prague Faculty of Science Albertov 6 12843 Prague 2, Czech Republic E-mail: masin@natur.cuni.cz
More information5. TWO-DIMENSIONAL FLOW OF WATER THROUGH SOILS 5.1 INTRODUCTION
5. TWO-DIMENSIONAL FLOW OF WATER TROUG SOILS 5.1 INTRODUCTION In many instances the flo of ater through soils is neither one-dimensional nor uniform over the area perpendicular to flo. It is often necessary
More informationThe Preliminary Study of the Impact of Liquefaction on Water Pipes
The Preliminary Study of the Impact of Liquefaction on Water Pipes Jerry J. Chen and Y.C. Chou ABSTRACT Damages to the existing tap-water pipes have been found after earthquake. Some of these damages are
More informationCompressibility & Consolidation
CHAPTER Compressibility & Consolidation Settlement If a structure is placed on soil surface, then the soil will undergo an elastic and plastic deformation. In engineering practice, the deformation or reduction
More informationBenefits of Collaboration between Centrifuge Modeling and Numerical Modeling. Xiangwu Zeng Case Western Reserve University, Cleveland, Ohio
Benefits of Collaboration between Centrifuge Modeling and Numerical Modeling Xiangwu Zeng Case Western Reserve University, Cleveland, Ohio ABSTRACT There is little doubt that collaboration between centrifuge
More informationVerification Manual GT
Verification Manual GT Written by: The SoilVision Systems Ltd. Team Last Updated: Tuesday, February 20, 2018 SoilVision Systems Ltd. Saskatoon, Saskatchewan, Canada Software License The software described
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationDeformation And Stability Analysis Of A Cut Slope
Deformation And Stability Analysis Of A Cut Slope Masyitah Binti Md Nujid 1 1 Faculty of Civil Engineering, University of Technology MARA (Perlis), 02600 Arau PERLIS e-mail:masyitahmn@perlis.uitm.edu.my
More information4 Undrained Cylindrical Cavity Expansion in a Cam-Clay Medium
Undrained Cylindrical Cavity Expansion in a Cam-Clay Medium 4-1 4 Undrained Cylindrical Cavity Expansion in a Cam-Clay Medium 4.1 Problem Statement The stress and pore pressure changes due to the expansion
More informationCavity Expansion Methods in Geomechanics
Cavity Expansion Methods in Geomechanics by Hai-Sui Yu School of Civil Engineering, University of Nottingham, U. K. KLUWER ACADEMIC PUBLISHERS DORDRECHT / BOSTON / LONDON TABLE OF CONTENTS Foreword Preface
More informationTheoretical Model and Numerical Analysis on Unsaturated Expansive Soil Slope during Digging and Climate Change Courses (II)-Numerical Analysis
The 12 th International Conference of International Association for Computer Methods and Advances in Geomechanics (IACMAG) 1-6 October, 2008 Goa, India Theoretical Model and Numerical Analysis on Unsaturated
More information4 One-dimensional consolidation
4 One-dimensional consolidation In the previous chapter e studied briefly the conditions of the steady flo of ater through a stationary soil structure; in this chapter e shall be concerned ith the more
More informationTable of Contents Chapter 1 Introduction to Geotechnical Engineering 1.1 Geotechnical Engineering 1.2 The Unique Nature of Soil and Rock Materials
Table of Contents Chapter 1 Introduction to Geotechnical Engineering 1.1 Geotechnical Engineering 1.2 The Unique Nature of Soil and Rock Materials 1.3 Scope of This Book 1.4 Historical Development of Geotechnical
More informationEvaluation of dynamic behavior of culverts and embankments through centrifuge model tests and a numerical analysis
Computer Methods and Recent Advances in Geomechanics Oka, Murakami, Uzuoka & Kimoto (Eds.) 2015 Taylor & Francis Group, London, ISBN 978-1-138-00148-0 Evaluation of dynamic behavior of culverts and embankments
More informationDynamic consolidation problems in saturated soils solved through u w formulation in a LME meshfree framework
Dynamic consolidation problems in saturated soils solved through u w formulation in a LME meshfree framework Pedro Navas a, Rena C. Yu a, Susana López-Querol b,, and Bo Li c a E. T. S. de Ingenieros de
More information10th Asian Regional Conference of IAEG (2015)
0th Asian Regional Conference of IAEG (05) Normalized Pore Water Pressure Ratio and Post-Cyclic Settlement of Saturated Clay Subjected to Undrained Uni-Directional and Multi-Directional Cyclic Shears TRAN
More information3-D FINITE ELEMENT NONLINEAR DYNAMIC ANALYSIS FOR SOIL-PILE-STRUCTURE INTERACTION
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-, 4 Paper No. 157 3-D FINITE ELEMENT NONLINEAR DYNAMIC ANALYSIS FOR SOIL-PILE-STRUCTURE INTERACTION B.K. MAHESHWARI 1,
More informationIntro to Soil Mechanics: the what, why & how. José E. Andrade, Caltech
Intro to Soil Mechanics: the what, why & how José E. Andrade, Caltech The What? What is Soil Mechanics? erdbaumechanik The application of the laws of mechanics (physics) to soils as engineering materials
More informationSettlement and Bearing Capacity of a Strip Footing. Nonlinear Analyses
Settlement and Bearing Capacity of a Strip Footing Nonlinear Analyses Outline 1 Description 2 Nonlinear Drained Analysis 2.1 Overview 2.2 Properties 2.3 Loads 2.4 Analysis Commands 2.5 Results 3 Nonlinear
More informationEXTENDED ABSTRACT. Dynamic analysis of elastic solids by the finite element method. Vítor Hugo Amaral Carreiro
EXTENDED ABSTRACT Dynamic analysis of elastic solids by the finite element method Vítor Hugo Amaral Carreiro Supervisor: Professor Fernando Manuel Fernandes Simões June 2009 Summary The finite element
More informationTowards Efficient Finite Element Model Review Dr. Richard Witasse, Plaxis bv (based on the original presentation of Dr.
Towards Efficient Finite Element Model Review Dr. Richard Witasse, Plaxis bv (based on the original presentation of Dr. Brinkgreve) Journée Technique du CFMS, 16 Mars 2011, Paris 1/32 Topics FEA in geotechnical
More informationExamining the Soil Responses during the Initiation of a Flow Landslide by Coupled Numerical Simulations
The 2012 World Congress on Advances in Civil, Environmental, and Materials Research (ACEM 12) Seoul, Korea, August 26-30, 2012 Examining the Soil Responses during the Initiation of a Flow Landslide by
More informationSoil strength. the strength depends on the applied stress. water pressures are required
Soil Strength Soil strength u Soils are essentially frictional materials the strength depends on the applied stress u Strength is controlled by effective stresses water pressures are required u Soil strength
More informationMeshfree Inelastic Frame Analysis
Theory & Results Louie L. Yaw, Sashi Kunnath and N. Sukumar University of California, Davis Department of Civil and Environmental Engineering Minisymposium 47 Recent Advances in Modeling of Engineering
More information2D Liquefaction Analysis for Bridge Abutment
D Liquefaction Analysis for Bridge Abutment Tutorial by Angel Francisco Martinez Integrated Solver Optimized for the next generation 64-bit platform Finite Element Solutions for Geotechnical Engineering
More informationCYCLIC LIQUEFACTION POTENTIAL OF LACUS- TRINE CARBONATE SILT FROM JULIAN ALPS
CYCLIC LIQUEFACTION POTENTIAL OF LACUS- TRINE CARBONATE SILT FROM JULIAN ALPS BOJAN ŽLENDER and STANISLAV LENART About the authors Bojan Žlender University of Maribor, Faculty of Civil Engineering Smetanova
More informationCalculation of 1-D Consolidation Settlement
Calculation of 1-D Consolidation Settlement A general theory for consolidation, incorporating threedimensional flow is complicated and only applicable to a very limited range of problems in geotechnical
More informationFINITE ELEMENT SIMULATION OF RETROGRESSIVE FAILURE OF SUBMARINE SLOPES
FINITE ELEMENT SIMULATION OF RETROGRESSIVE FAILURE OF SUBMARINE SLOPES A. AZIZIAN & R. POPESCU Faculty of Engineering & Applied Science, Memorial University, St. John s, Newfoundland, Canada A1B 3X5 Abstract
More informationA truly meshless Galerkin method based on a moving least squares quadrature
A truly meshless Galerkin method based on a moving least squares quadrature Marc Duflot, Hung Nguyen-Dang Abstract A new body integration technique is presented and applied to the evaluation of the stiffness
More informationApplicability of Multi-spring Model Based on Finite Strain Theory to Seismic Behavior of Embankment on Liquefiable Sand Deposit
Applicability of Multi-spring Model Based on Finite Strain Theory to Seismic Behavior of Embankment on Liquefiable Sand Deposit Kyohei Ueda Railway Technical Research Institute, Kokubunji, Tokyo, Japan
More informationCalculation types: drained, undrained and fully coupled material behavior. Dr Francesca Ceccato
Calculation types: drained, undrained and fully coupled material behavior Dr Francesca Ceccato Summary Introduction Applications: Piezocone penetration (CPTU) Submerged slope Conclusions Introduction Porous
More informationConstitutive modelling of fabric anisotropy in sand
Geomechanics from Micro to Macro Soga et al. (Eds) 2015 Taylor & Francis Group, London, ISBN 978-1-138-02707-7 Constitutive modelling of fabric anisotropy in sand Z.W. Gao School of Engineering, University
More informationTime Rate of Consolidation Settlement
Time Rate of Consolidation Settlement We know how to evaluate total settlement of primary consolidation S c which will take place in a certain clay layer. However this settlement usually takes place over
More information3-D FINITE ELEMENT MODELING OF THE REMOTE FIELD EDDY CURRENT EFFECT
3-D FINITE ELEMENT MODELING OF THE REMOTE FIELD EDDY CURRENT EFFECT Y. Sun and H. Lin Department of Automatic Control Nanjing Aeronautical Institute Nanjing 2116 People's Republic of China Y. K. Shin,
More informationA thermo-hydro-mechanically coupled analysis of clay using a thermo-elasto-viscoplastic model
JHUWS05 A thermo-hydro-mechanically coupled analysis of clay using a thermo-elasto-viscoplastic model by F. Oka, S. Kimoto, Y.-S. Kim, N. Takada Department of Civil & Earth Resources Engineering, Kyoto
More informationEARTH PRESSURES AGAINST RIGID RETAINING WALLS IN THE MODE OF ROTATION ABOUT BASE BY APPLYING ZERO-EXTENSION LINE THEORY
EARTH PRESSURES AGAINST RIGID RETAINING WALLS IN THE MODE OF ROTATION ABOUT BASE BY APPLYING ZERO-EXTENSION LINE THEORY Morshedi S.M., Ghahramani A. 2, Anvar S.A. 2 and Jahanandish M. 2 Department of Geotechnical
More information3-D Numerical simulation of shake-table tests on piles subjected to lateral spreading
3-D Numerical simulation of shake-table tests on piles subjected to lateral spreading M. Cubrinovski 1, H. Sugita 2, K. Tokimatsu 3, M. Sato 4, K. Ishihara 5, Y. Tsukamoto 5, T. Kamata 5 1 Department of
More informationCommon pitfalls while using FEM
Common pitfalls while using FEM J. Pamin Instytut Technologii Informatycznych w Inżynierii Lądowej Wydział Inżynierii Lądowej, Politechnika Krakowska e-mail: JPamin@L5.pk.edu.pl With thanks to: R. de Borst
More information3D ANALYSIS OF H-M COUPLED PROBLEM WITH ZERO-THICKNESS INTERFACE ELEMENTS APPLIED TO GEOMECHANICS
Environmental 3D analysis of H-M and Geosciences coupled problem with zero-thickness interface elements applied to Geomechanics XIII International Conference on Computational Plasticity. Fundamentals and
More informationMicro-mechanics in Geotechnical Engineering
Micro-mechanics in Geotechnical Engineering Chung R. Song Department of Civil Engineering The University of Mississippi University, MS 38677 Fundamental Concepts Macro-behavior of a material is the average
More informationCubzac-les-Ponts Experimental Embankments on Soft Clay
Cubzac-les-Ponts Experimental Embankments on Soft Clay 1 Introduction In the 197 s, a series of test embankments were constructed on soft clay at Cubzac-les-Ponts in France. These full-scale field tests
More informationPOST-CYCLIC RECOMPRESSION CHARACTERISTICS OF A CLAY SUBJECTED TO UNDRAINED UNI-DIRECTIONAL AND MULTI-DIRECTIONAL CYCLIC SHEARS
NCEE Tenth U.S. National Conference on Earthquake Engineering Frontiers of Earthquake Engineering July -5, 4 Anchorage, Alaska POST-CYCLIC RECOMPRESSION CHARACTERISTICS OF A CLAY SUBJECTED TO UNDRAINED
More informationThe Effects of Different Surcharge Pressures on 3-D Consolidation of Soil
The Effects of Different Surcharge Pressures on 3-D Consolidation of Soil Arpan Laskar *1 and Sujit Kumar Pal 2 *1 Department of Civil Engineering, National Institute of Technology Agartala, Tripura, India.
More informationEstimation of Multi-Directional Cyclic Shear-Induced Pore Water Pressure on Clays with a Wide Range of Plasticity Indices
Proceedings of the 2 nd International Conference on Civil, Structural and Transportation Engineering (ICCSTE 16) Ottawa, Canada May 5 6, 216 Paper No. 116 Estimation of Multi-Directional Cyclic Shear-Induced
More informationDynamic Analysis Contents - 1
Dynamic Analysis Contents - 1 TABLE OF CONTENTS 1 DYNAMIC ANALYSIS 1.1 Overview... 1-1 1.2 Relation to Equivalent-Linear Methods... 1-2 1.2.1 Characteristics of the Equivalent-Linear Method... 1-2 1.2.2
More informationThe Study of Large-Scale Triaxial Texts on Metamorphic Soft Rock Embankment Filler
1 st International Conference on Transportation Infrastructure and Materials (ICTIM 2016) ISBN: 978-1-60595-367-0 The Study of Large-Scale Triaxial Texts on Metamorphic Soft Rock Embankment Filler Yanyan
More informationD1. A normally consolidated clay has the following void ratio e versus effective stress σ relationship obtained in an oedometer test.
(d) COMPRESSIBILITY AND CONSOLIDATION D1. A normally consolidated clay has the following void ratio e versus effective stress σ relationship obtained in an oedometer test. (a) Plot the e - σ curve. (b)
More informationVapor Pressure Prediction for Stacked-Chip Packages in Reflow by Convection-Diffusion Model
Vapor Pressure Prediction for Stacked-Chip Packages in Reflo by Convection-Diffusion Model Jeremy Adams, Liangbiao Chen, and Xuejun Fan Lamar University, PO Box 10028, Beaumont, TX 77710, USA Tel: 409-880-7792;
More informationReinforced Soil Structures Reinforced Soil Walls. Prof K. Rajagopal Department of Civil Engineering IIT Madras, Chennai
Geosynthetics and Reinforced Soil Structures Reinforced Soil Walls continued Prof K. Rajagopal Department of Civil Engineering IIT Madras, Chennai e-mail: gopalkr@iitm.ac.inac in Outline of the Lecture
More informationThe method by which this error has to be combined with other sources of error is I
The error in the discharge coefficient (including C,) of a triangular broad-crested (truncated) eir, hich has been constructed ith reasonable care and skill, may be deduced form the equation X, = f (3
More informationNumerical Modeling of Direct Shear Tests on Sandy Clay
Numerical Modeling of Direct Shear Tests on Sandy Clay R. Ziaie Moayed, S. Tamassoki, and E. Izadi Abstract Investigation of sandy clay behavior is important since urban development demands mean that sandy
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationThe Preliminary Study of the Impact of Liquefaction on Water Pipes
The Preliminary Study of the Impact of Liquefaction on Water Pipes Jerry J. Chen and Y.C. Chou Geotechnical Engineer, Dept. of Geotechnical Engineering, CECI Engineering Consultants, Inc. CONTENT 1. Introduction
More informationSHEAR STRENGTH OF SOIL
Soil Failure Criteria SHEAR STRENGTH OF SOIL Knowledge about the shear strength of soil important for the analysis of: Bearing capacity of foundations, Slope stability, Lateral pressure on retaining structures,
More informationNUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS
IGC 009, Guntur, INDIA NUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS Mohammed Younus Ahmed Graduate Student, Earthquake Engineering Research Center, IIIT Hyderabad, Gachibowli, Hyderabad 3, India.
More informationPANCRAS MUGISHAGWE BUJULU AND GUSTAV GRIMSTAD
DOI: http://dx.doi.org/1.4314/gjer.v11i1.7 GLOBAL JOURNAL OF ENGINEERING RESEARCH VOL 11, NO. 1, 212: 57-66 COPYRIGHT BACHUDO SCIENCE CO. LTD PRINTED IN NIGERIA. ISSN 1596-292x.globaljournalseries.com;
More informationSoft Ground Coupled Consolidation
1 Introduction Soft Ground Coupled Consolidation This example is about constructing an embankment in delayed stages on a soft foundation so that some of the excess pore-pressure is allowed to dissipate
More informationFormation rock failure mode and its recognition method under high pressure squeezing condition in horizontal wells
01 36 Vol. 36 No. 5 5 Journal of China University of Petroleum Oct. 01 1673-5005 01 05-0105-05 1 3 4 1 1 1. 66580. 30045 3. 710018 4. 73600 TE 358. 1 A doi 10. 3969 /j. issn. 1673-5005. 01. 05. 019 Formation
More informationShear Strength of End Web Panels
Paper 10 Shear Strength of End Web Panels Civil-Comp Press, 2012 Proceedings of the Eleventh International Conference on Computational Structures Technology, B.H.V. Topping, (Editor), Civil-Comp Press,
More informationStatic Pile Head Impedance using 3D Nonlinear FEM Analysis
Static Pile Head Impedance using 3D Nonlinear FEM Analysis Ichiro NAGASHIMA Technology Center, Taisei Corporation, 344-1 Nasecho, Totsuka-ku, Yokohama 245-51, Japan, ichiro.nagashima@sakura.taisei.co.jp
More informationOn equal settlement plane height in piled reinforced embankments
Available online www.jocpr.com Journal of Chemical and Pharmaceutical Research, 2014, 6(4):23-29 Research Article ISSN : 0975-7384 CODEN(USA) : JCPRC5 On equal settlement plane height in piled reinforced
More informationComputation Time Assessment of a Galerkin Finite Volume Method (GFVM) for Solving Time Solid Mechanics Problems under Dynamic Loads
Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 215 Paper o. 31 Computation Time Assessment of a Galerkin Finite Volume
More informationUSER S MANUAL 1D Seismic Site Response Analysis Example University of California: San Diego August 30, 2017
USER S MANUAL 1D Seismic Site Response Analysis Example http://www.soilquake.net/ucsdsoilmodels/ University of California: San Diego August 30, 2017 Table of Contents USER'S MANUAL TABLE OF CONTENTS Page
More informationFinite Element Instability Analysis of the Steel Joist of Continuous Composite Beams with Flexible Shear Connectors
Journal of Solid Mechanics Vol. 9, No. (07) pp. 5-69 Finite Element Instability Analysis of the Steel Joist of Continuous Composite Beams ith Fleible Shear Connectors H. Bakhshi,*, H.R. Ronagh Engineering
More informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
More information