MINISTRY OF TRAINING AND EDUCATION MINISTRY OF CONSTRUCTION HANOI ARCHITECTURAL UNIVERSITY NGO QUOC TRINH A STUDY ON PILE WORKING UNDER HORIZONTAL LOAD AND SEISMIC LOAD Major: Civil and Industrial Construction Engineering Code : 62 58 02 08 SUMMARY OF Eng.D THESIS HÀ NOI - 2014 The study was completed at: HANOI ARCHITECTURAL UNIVERSITY Scientific supervisors: 1. Supervisor: Asst. Prof. Dr Vuong Van Thanh Dr. Tran Huu Ha Reviewer 1: Prof. DSc Nguyen Dang Bich Reviewer 2: Prof.Ph.D Do Nhu Trang Reviewer 3: Asst. Prof. Dr Trinh Minh Thu The thesis will be defended in front of Eng.D Assessment Council at University level held in Hanoi Architectural University At: date month year 2014. Thesis can be found at: • National Library of Vietnam • Library of Hanoi Architectural University LIST OF WORKS BY AUTHOR 1. Ngo Quoc Trinh (2008), Study on the interaction problem between shallow foundation and ground deformation, Vietnam Road and Bridge Magazine 2. Ngo Quoc Trinh, Vuong Van Thanh, Tran Huu Ha (5/2012). Study on the interaction problem between ground mass and elastic foundation underhorizontal static load. Vietnam Road and Bridge Magazine. 3. Ngo Quoc Trinh, Vuong Van Thanh, Tran Huu Ha (6/2012). Study on the interaction problem between single pile and elastic foundation underhorizontal static load. Vietnam Road and Bridge Magazine, 4. Ngo Quoc Trinh, Vuong Van Thanh, Tran Huu Ha (11/2012), Using the solution of Mindlin to build the interaction problem between the pile and the elastic foundation under horizontal static load, Collection of Conference on materials science, structure and construction technology in 2012 (MSC2012), Hanoi Architectural University. 5. Ngo Quoc Trinh (12/2012), Using the comparison method to study the interaction problem between the pile and the elastic foundation under horizontal static load, 9th Nationwide Mechanics Conference. 6. Ngo Quoc Trinh (3/2013), Study on Love wave transmission problem in the foundation when earthquake occcurs. Tranposrt Journal. 1 INTRODUCTION 1. Background Though Vietnam is not located in the Ring of Fire of large earthquake areas in the world, it is still affected by strong earthquakes because the territory of Vietnam exists many faults acting complexly such as Lai Chau - Dien Bien fault, Ma River fault, Son La fault, Hong River fault zone, Ca River fault zone (in the history, there was a strong earthquake on 6.8 Richter scale). In order to design earthquake resistance for buildings, our country now is using some foreign translated standards: TCXDVN 375: 2006; 22 TCN 221-95 ; TCXD 205-1998 ; 22TCN 272- 05; however there is little detailed guidance on calculating the interaction between the building and the foundation. The biggest difficulty when designing pile foundations under horizontal load and seismic load is assessing the interaction between the pile and the foundation. Because the interaction between the pile and the foundation is too complex, the current calculating methods are often simplified by the model (Winkler model; continuous elastic model).Therefore, it is difficult to determine the interaction coefficients between the pile and the foundation (spring coefficient, viscosity coefficient), it is difficult to ensure the boundary conditions as well as endless radiation conditions; and the interaction between the pile and the foundation is insufficient and it is considered only in plane strain problem From above analysis, the study on the work of the pile, in which the study on the interaction between the pile and the foundation under horizontal load and seismic load is necessary, scientific and practical, contributing a fuller review of the calculating method of the pile foundation of constructions in Vietnam. 2 Study objectives Formulate theoretical methods to study interaction problems between pile – soil foundation and formulate calculating softwares to determine the status of the stress- strain of pile under horizontal load and seismic load. 3. Subject and scope of the study 2 The thesis studies vertical single pile located within an infinite elastic space under the effect of horizontal static load, horizontal dynamic load and seismic load. The thesis does not calculate in other foundation models (elastoplastic, visco-elastic), does not consider the liquefy phenomenon in the foundation when an earthquake occurs, does not consider the effect of pore water pressure in saturated foundation and does not study the limit state problem of pile. 4. Study content Study stress-strain state of the soil mass under horizontal static load. Study static interaction problem between pile and foundation under horizontal static load. Study dynamic interaction problem between the pile and the foundation under horizontal dynamic load and seismic load in the frequency domain and time domain. Formulate calculating softwares for above cases of study. 5 Study method Formulate theoretical problem by using comparison method of extreme principle Gauss (hereinafter refering to EP Gauss) when using the static solution of the semi-infinite elastic space (for all static interaction problem) and dynamic solution of infinite elastic space (for the dynamic interaction problem) as comparison system. Using the finite limit element method to solve and basing on the numerical results obtains the results proving the correctness and reliability of theoretical calculation. Chapter 1 OVERVIEW OF STUDY METHODS ON THE INTERACTION BETWEEN PILE AND FOUNDATION UNDER HORIZONTAL LOAD Basing on the annalysis of study methods on the interaction between pile and foundation under horizontal load, we can draw some comments as follows: + First: It’s dificult to determine stiffness coefficient “linear springs”, “nonlinear springs” (p-y curve), viscosity coefficient . 3 + Second: It’s dificult to determine the boundary conditions at infinity, especially for the wave transmission problem when an earthquake occurs. + Third: The interaction between the pile and the soil has not been adequately considered, only consider the effect of soil on the pile without considering the effect of pile on the soil. + Fourth: Mainly study on plane strain problem. From above issues, the author has based on the method of using comparison system of extreme principle Gauss method to formulate the interaction problem between pile and foundation under static load, dynamic horizontal load and seismic load with full consideration of the boundary conditions and endless radiation conditions as well as consideration of the full interaction between the pile and the soil and consider the 3-dimensional problem. Chapter 2 STUDY ON STRESS-STRAIN OF FOUNDATION UNDER STATIC HORIZONTAL LOAD 2.1 The basic equations and the wave propagation equation of elastic medium. 2.1.1 The basic contact of elastic medium 2.1.2 Formulate diferential balance equation and wave propagation equations according to EP Gauss 2.1.2.1 Extreme principle Gauss Extreme principle Gauss (EP Gauss) is an extreme principle of mechanics stated by Gauss K.F (1777 - 1855) in 1829 with the following content [5][6],[61]: “The motion of the system of points, optionally linked influenced by any force, in every moment occurred in accordance with the highest possible ability to move but the quality can be done if they are completely free, which means it occurs with the smallest amount of coercion if forced volume measurements in very little time period taken by the total volume of the mass of each point of the square deviation of the position they compared with when they are free”. Expression of forced amount in geometry form of EP Gauss is written as follows: 4 2 ii i i CBm ∑ ⇒ min! (2.4) Here 2 ii CB is the distance between 2 points B i and C i of material point i with the volume m i . B i is the position which material point i obtain when moving freely and C i is the position when that material point moves with the connection after extremely little time dt. Symbol Σ is the total number of points obtained from the system. EP Gauss is applied for material point system. Basing on this principle, in 1979, Prof.Ph.D Ha Huy Cuong recommended using Extreme principle Gauss to solve the problem of mechanical deformation of solid objects. 2.1.2.2 Formulate diferential balance equation Departs from Helmholtz theorem [60], for continuous environment, establishing three movements: forward motion, motion and rotation distortion. From EP Gauss for mechanics material points, applying EP Gauss method for moving elastic deformation elements distributed 3D, the author received three diferential balance equations of elastic systems (Navier equation) like diferential balance equations presented in many documents about elastic theory [26],[46],[53],[60]. 2.1.2.3 Formulate wave propagation equation Applying EP Gauss method for motion of volumetric strain and rotation like abrasive of distributed elements around the axis of x, y, z, the author got 4 of wave propagation equations (2.25), (2.33), (2.36), (2.37). Thus Navier equations or wave propagation equatations to study the motion of elastic medium. 2.2 The solution for an infinite elastic space and half-infinite elastic space 2.2.1 The solution for an infinite elastic space (solution of Kelvin) 2.2.2 The solution for half-infinite elastic space (solution of Mindlin ) 2.3 Formulate the interaction problem between elastic soil mass and half-infinite elastic space. 5 2.3.1 Comparison system is infinite elastic semi-space. In terms of rectangular soil mass V with elastic parameters E 1 , ν 1 in the elastic half-space with the elastic parameters E 0 , ν 0 . Horizontal force P affects in or outside the ground. Considering the comparision system as the infinite half space elastic with the elastic parameters E 0 , ν 0 , also under P horizontal force affected as the system needs calculating (Figure 2.5). Figure 2.4 Model of problem ofcalculating the elastic soil mass in the infinite elastic half-space Figure 2.5 Comparison system is half-infinite elastic space Note that on the boundary of the soil mass needed calculating, there is prestress σ ij affect (figure 2.4) and on the boundary of the soil mass of comparison system, there is prestress σ ij 0 affects (figure 2.5). Using prestress state σ ij 0 of the known comparison system to calculate the prestress state σ ij of the system needed calculating by wrtiting forced density as follows: Z V = ⌡ ⌠ V* (σ x -σ x 0 ) ε x dV * + ⌡ ⌠ V* (σ y -σ y 0 ) ε y dV * + ⌡ ⌠ V* (σ z -σ z 0 ) ε z dV * + ⌡ ⌠ V* (τ xy -τ xy 0 ) γ xy dV * + ⌡ ⌠ V* (τ xz -τ xz 0 ) γ xz dV * + ⌡ ⌠ V* (τ yz -τ yz 0 ) γ yz dV * →min (2.50) In (2.50), V * is volume of extended soil to consider the boundary conditions; V is volume of soil mass needed calculating (V * > V); ε x , ε y , ε z , γ xy , γ xz , γ yz are the strain of the soil mass; σ x 0 , σ y 0 , σ z 0 , τ xy 0 , τ xz 0 , τ yz 0 are the stress-strain of definite comparison system according to the solution of (figure 2.5); prestresses σ x , σ y , σ z , τ xy , τ xz , τ yz are the stress-strain of the soil mass of calculated system (figure 2.4). Replace the distortions by the contacts (2.1). EP Gauss considers actual displacement u, v, w in (2.51) as the virtual displacement, i.e whether the distortion is independent on the prestress, the extreme conditions of (2:51) is written as follows: Extended area to examine boundary condition P Soil mass E 0 , E 1 , ν c Compared soil mass P E 0 , ν E 0 , ν c 6 δZ V = ⌡ ⌠ V* (σ x -σ x 0 ) δ( ∂u ∂x )dV * + ⌡ ⌠ V* (σ y -σ y 0 ) δ( ∂v ∂y )dV * + ⌡ ⌠ V* (σ z -σ z 0 ) δ( ∂w ∂z )dV * + ⌡ ⌠ V* (τ xy -τ xy 0 ) δ ( ∂u ∂y + ∂v ∂x ) dV * + ⌡ ⌠ V* (τ xz -τ xz 0 ) δ ( ∂u ∂z + ∂w ∂x ) dV * + ⌡ ⌠ V* (τ yz -τ yz 0 ) δ ( ∂v ∂z + ∂w ∂y ) dV * = 0 (2.52) in which δ is marks obtaining variational. Note that soil mass here contains three implicits u, v, w, so from (2.52) we obtained 3 equation system: ⌡ ⌠ V* (σ x -σ x 0 ) δ( ∂u ∂x )dV * + ⌡ ⌠ V* (τ xy -τ xy 0 ) δ( ∂u ∂y ) dV * + ⌡ ⌠ V* (τ xz -τ xz 0 ) δ( ∂u ∂z )dV * = 0 ⌡ ⌠ V* (σ y -σ y 0 ) δ( ∂v ∂y )dV * + ⌡ ⌠ V* (τ xy -τ xy 0 ) δ( ∂v ∂x )dV * + ⌡ ⌠ V* (τ yz -τ yz 0 ) δ( ∂v ∂z )dV * =0 (2.53) ⌡ ⌠ V* (σ z - σ z 0 ) δ( ∂w ∂z )dV * + ⌡ ⌠ V* (τ xz -τ xz 0 ) δ( ∂w ∂x )dV * + ⌡ ⌠ V* (τ yz -τ yz 0 ) δ( ∂w ∂y )dV * = 0 Make variational calculations [34] for the (2.53) we received three following equations: ∂σ x ∂x + ∂τ xy ∂y + ∂τ xz ∂z = ∂σ x 0 ∂x + ∂τ xy 0 ∂y + ∂τ xz 0 ∂z ∂σ y ∂y + ∂τ xy ∂x + ∂τ yz ∂z = ∂σ y 0 ∂y + ∂τ xy 0 ∂x + ∂τ yz 0 ∂z (2.54) ∂σ z ∂z + ∂τ xz ∂z + ∂τ yz ∂y = ∂σ z 0 ∂z + ∂τ xz 0 ∂z + ∂τ yz 0 ∂y The right side (2.54) satisfies the balance equation when a horizontal force P in comparison system causes (Figure 2.4), so the left side of (2.54) is also the balance equation when a horizontal force P affecting in the calculated system (Figure 2.3) causes. Thus, by using the comparison system, we obtained three equilibrium equations of the system needed calculating. 2.3.2 Comparison system is infinite elastic space - Considering the case in which force P effects horizontally on soil V (Figure 2.8a). AB surface is the free surface. Let horizontal force P acting on the elastic space, using Kelvin solution to calculate the prestress state σ ij 0 in it.Because the system needed calculating lies in half-space (Figure 2.8a), we should only use the bottom half of infinite space (Figure 2.8b). 7 (a) (b) (c) Figure 2.8 The problem model of soil mass under horizontal force when using the comparison system as infinite elastic space The stress state σ ij 0 is equivalent to the force P / 2, so we have to put 2 forces P to calculate prestress σ ij 0 according to Kelvin solution. In case of the horizontal force P placed at the depth c compared with the free surface, we use 2 forces P symmetrically placed symmetrically through surface AB (Figure 2.8c). When calculating the above diagram, on the surface AB there is also the effect of prestress σ z 0 . Mindlin solution for the elastic half-space under horizontal force P derives from Kelvin solution with calculating diagram as in Figure 2.8c and finds the way to ensure σ z 0 = 0 on the surface AB. The solution obtained is the calculus solution. The author uses the diagram in figure 2.8c to calculate σ ij 0 . Due to the effect of prestress σ z 0 on the surface AB of the bottom half, it is necessary to consider the effect of this variable by writing forced amount as follows: Z AB = ⌡ ⌠ Ω AB [(σ z -σ z 0 ) w dΩ AB → min (2.55) with Ω AB the surface area of AB. Besides, it is necessary to ensure the condition of σ z = 0 on the surface AB. In a nutshell, the problem determining the prestress state of soil mass V when using Kelvin solution is written as follows: Z = Z V + Z AB → min (2.56) With constraint σ z = 0 on the surface AB. Z V = ⌡ ⌠ V* (σ x -σ x 0 ) ε x dV * + ⌡ ⌠ V* (σ y -σ y 0 ) ε y dV * + ⌡ ⌠ V* (σ z -σ z 0 ) ε z dV * + ⌡ ⌠ V* (τ xy -τ xy 0 ) γ xy dV * + ⌡ ⌠ V* (τ xz -τ xz 0 ) γ xz dV * + ⌡ ⌠ V* (τ yz -τ yz 0 ) γ yz dV * → min (2.57) P P P c c P A B A B Examined soil A B σ ij 0 σ z 0 8 In (2.57) prestresses σ x 0 , σ y 0 , σ z 0 , τ xy 0 , τ xz 0 , τ yz 0 are the prestress state of the comparison system determined according to Kelvin solution with two forces P (Figure 2.8c). By writing extended functional Lagrange, we turn constrained extreme problem into unconstrained extreme problem as follows: F = Z V + Z AB + λσ z → min (2.58) λ = λ(x,y) is a Lagrange factor as a new implicit function of the problem. Extreme conditions of F would be: δF = δZ V + δZ AB + δλσ z = 0 (2.59) 2.4 Solve the problem by using the finite element method The soil mass to be calculated as well as the soil mass of the comparison system is divided into rectangular elements ( 3D problem ) having any particle size. In order to consider the boundary conditions on the system to be calculated, the comparison system has 1 more number of elements than the system to be calculated according to the depth z and direction x, direction y . We can use rectangular elements with 8 nodes [38], but to get a better approximation, the author used rectangular elements with 20 nodes in the natural coordinate system with particle size ∆x = ∆y = ∆z = 2 and using displacements as unknowns . Each node has 3 parameters (unknown) to be determined as the displacement u according to direction x, v according to direction y, w according to direction z. Thus, the elements has 3 x 20 = 60 displacement parameters (60 unknowns) to be determined. Knowing the displacement of the nodes, the displacements at any point within the element is determined according to the interpolation function [39], [60] 2.5 Check the results and review 2.5.1 Problem using the comparison system as infinite elastic semi-space Considering the interaction problem between soil mass V having elastic parameters E 1 , ν 1 with infinite elastic half space having elastic parameters E 0 ,ν 0 (Figure 2.12). Based on Matlab software, the author formulated the calulating program Mstatic1 to survey following cases: 9 * Case 1: Put E 1 = E 0 , ν 1 = ν 0 (a) (b) Figure 2.13 Horizontal displacement chart of the soil mass when horizontal force P affects on surface (a) andbottom (b) of the land mass, in case E 1 = E 0 ; ν 1 = ν 0 . Recognizing that the results calculated by EP Gauss completely coincide with the results of the calculus solution of Mindlin (see Appendix 1) When changing the volume of block V, even in case that block V only has 1 element, we still get accurate results. * Case 2: Put ν 1 = ν 0 ; E 1 ≠ E 0 (E1 retained as in case 1, E 0 changed of the comparison system) (a) (b) Figure 2.14 Horizontal displacement chart of the soil mass when horizontal force P affects on surface (a) andbottom (b) of the soil mass, in case ν 1 = ν 0 ; E 1 ≠ E 0 Here we find a entirely concurrence between two results according to the EP Gauss solution in case 1 and case 2. When the volume V changes, we still get accurate results as above. Thus, through two survey cases we found that, though the comparison system has the same or different elastic module compared with the elastic module of comparison system to be 10 calculated, the displacement results of the system to be calculated is constant. This shows the correctness and reliability of theoretical calculations. 2.5.2 The problem of comparison system is infinite elastic space Surveying the land mass having E 1 = E 0 , μ 1 = μ 0 when let the horizontal force P affect alternately on 3 locations: c = 0 (the free surface of the land mass), c = 3m c = 5.4 m (bottom of the land mass) by two calculating programs Mstatic1 (comparison system is the infinite semi-space); Kstatic1 (comparison system is infinite space) (a) (b) (c) Hình 2.18 Horizontal displacement chart of the soil mass calculated according to 2 programs Mstatic1 and Kstatic1 when horizontal load P affects on the position c=0 (a); c=3m (b); c=5.4m (c) Calculating results show that displacement of the land mass when calculating according to Kstatic1 is approximately equal to the displacement of the land mass calculating according to Mstatic1 with the largest error of about 6% and the deeper the force is put compared with the free surface, the smaller the error between the two results is and almost overlap. Thus, through the numerical solution by finite element method, we can turn the solution of infinite elastic space (Kelvin solution) into the solution of infinite elastic senmi- space (Mindlin solution). 2.6 Conclusion of chapter 2 11 1- Formulating the interaction problem between the land mass with the remain infinite elastic semi-space under static horizontal load. With the conditions of displacement and prestress on the boundary surfaces of the land mass automatically satisfied exactly, there is no need to add extra links (i.e spring links) as the current methods and the condition at infinity is automatically satisfied. 2 - Formulating calculating program by using the finite element method in Matlab environment to calculate the land mass. Here use the rectangular element with 3-D, 20 nodes. Checking the numerical solution, we found a good fit between the calculated results with calculus solutions. 3 - Through numerical solution, we can turn the solution of infinite elastic space (Kelvin solution) into solution of half infinite elastic semi-space (Mindlin solution). Chapter 3 STUDY ON THE INTERACTION PROBLEM BETWEEN PILE AND FOUNDATION UNDER THE HORIZONTAL STATIC LOAD 3.1 Timoshenko beam theory Timoshenko beam theory is the bending beam theory considering the horizontal shear strain. Beam theory which considers current horizontal shear strain used two implicit functions u c (z); φ c (z) is independent implicit function often leads to Shear locking phenomenon (Shear locking). In the thesis, the author used Timoshenko beam theory but using two implicit functions which are the deflection u c (z) and shear force Q(z) in the pile. According to this method, there will be no longer Shear locking phenomenon. 3.2 Formulating bending beam problem considering horizontal shear strain according to EP Gauss By EP Gauss, the author formulated properly the deflection equation of bending beam considering the horizontal shear strain. 3.3 Finite element method for beam considering the horizontal shear strain Because there are two implicit functions, displacement function and shear force function of the beam, there are two types of elements: displacement element and shear force element. 12 Displacement element includes 2 nodes, each node has two unknown displacement and rotation, shear force element includes 3 nodes, each node has one unknown shear force. And ground element is rectangular element with 20 nodes, each node has three displacements u, v, w. 3.4 Formulating interaction problem between single pile with foundation under horizontal static load 3.4.1 In case of using comparison system as infinite elastic semi-space (a) System to be calculated (b) System to be compared According EP Gauss, forced density Z of the problem includes two components: Z = Z d + Z c → min Z d : forced account considering prestress state of the land mass of the comparison system which affects on the land mass containing piles of the system to be calculated (2.50). Z c considering the forced amount of bending pile considering horizontal shear strain γc in the pile. Z c = ⌡ ⌠ l Mχ c dz + ⌡ ⌠ l Qγ c dz (3.46) The condition ensuring the simultaneous work of the pile under the horizontal force compared with the foundation is that the horizontal displacement of the pile u c is equal to the horizontal displacement of the foundation u at the heart of the pile. We have: u c (z, x c , y c ) = u(z, x c , y c ) (3.49) We can lead the constrained extreme problem to unconstrained extreme problem by using Lagrange factor λ (z). The function λ (z) is implicit function to be calculated which changes according to the length of the pile. The Lagrange extended function F is now written as follows: F= Z d + Z c + ⌡ ⌠ l λ (z) (u c -u)dz → min (3.50) Extended area to examine boundary condition P Land mass contains pile Pile E 0 , ν 0 E 1 , ν 1 P E 0 , ν 0 E 0 , ν 0 Comparative soil mass 13 3.4.2 In case of using comparison system as infinite elastic space According EP Gauss, forced density Z of the problem includes two components: Z = Z d + Z c → min (3.53) Z d forced account considering prestress state of the land mass of the comparison system which affects on the land mass containing piles of the system to be calculated: Z d = Z V + Z AB ; Z V is the forced amount to calculate land mass V (formula 2.50); Z AB is the forced amount considering surface condition AB of the bottom half of the land mass: Z AB = ⌡ ⌠ Ω AB (σ z -σ z 0 ) w dΩ AB (3.57) Z c is forced amount (motion) of bending pile considering horizontal shear strain γ c in piles (formula 3.46). Constrained conditions u c (z, x c , y c ) = u(z, x c , y c ) and σ z = 0 on the free surface. We can lead the constrained extreme problem (3.53) to unconstrained extreme problem by using Lagrange factor λ as follows: F = Z d + Z c + ⌡ ⌠ l λ 1 (z) (u c -u)dz + ⌡ ⌠ Ω AB λ 2 (x,y) σ z d Ω AB → min (3.59) 3.5 Surveying some cases to test the reliability of the calculating program 3.5.1 Compare the results when let the elastic modules of comparison system be diferent a) b) Figure 3.9 Horizontal displacement diagram (a), bending moment (b) of the pile calculated according to two cases that the comparison system has E 0 = 10MPa; E 0 = 20MPa Realizing that the results of two cases are the same.Thus, the Such internal force displacement of pile in the system to be calculated does 14 not depend on the elastic module of the comparison systems, which proves that the algorithms provided is absolutely right. 3.5.3 Survey the problem compared with the method of Zavriev (1962) based on the local deformation foundation model [16] The author used input parameters of example V.5 in [16] calculated according to the method of Zavriev to calculate according to EP Gauss, then compared their results with each other. Table 3.5 Displacement value, maximum bending moment according to the method of Zavriev and EP Gauss Result Method Maximum displacement on the head of the pile (m) Maximum bending moment (kN.m) Pile under the load P, M Pile under the load P Pile under the load P, M Pile under the load P Zavriev 0,0116 0,0093 123,9 80,571 EP Gauss 0,0097 80,340 Comment: The displacement on the head of the pile, the maximum bending moment calculated according to the method of Zavriev is approximately equal to the results of The displacement on the head of the pile, the maximum bending moment calculated according to Gauss (about 4.1% error) 3.5.4 Survey the problem compared with the method of Poulos (1971) basing on continuous elastic foundation model [50] The author used the input parameters of the example 6.10 in [50] calculated according to the method of Poulos to calculate according to EP Gauss and then compared their results with each other. Table 3.6 Maximum displacement value on the head of the pile according to the method of Poulos and EP Gauss Result Method Maximum displacement on the head of the pile (cm) Pile under the load P, M Pile under the load P Poulos 5,8 4,2 EP Gauss 4,7 The displacement on the head of the pile calculated according to the method of Poulos is nearly equal to the displacement on the head of the pile calculated according to EP Gauss (12.7% error) 15 3.5.5 Survey the problem compared with Kim's research results based on the method of using p-y curve [45] The author used the input parameters in the study of soft pile of Kim [45] to formulate the KstaticPLs software calculated according to EP Gauss then compared their results with each other. (a) (b) Figure 3.14 Horizontal displacement chart, bending moment of piles calculated according to KstaticPLs (a); Kim, O’Neill, Matlock [45] (b) under the horizontal force effects in turn: 200kN, 400kN, 600kN, 800 kN. Comment: Results of displacement, bending moment of the problem based on the author's solution (KstaticPLs) are consistent with findingf of Kim, O'Neill, Matlock in the different cases of putting force both in shape, values and position that reaches the maximum value, minimum value, bending point. 3.6 Survey the parameters affecting the working of single pile under horizontal static load 3.6.1 Survey the change of pile length in homogeneous elastic foundation. Survey the short piles, long piles with reinforced concrete cross section (40x40) cm with the elastic module E c = 30.000MPa. Piles 16 have 2 different lengths: l =4m and l=16m and under the horizontal force P = 20kN on the head of the pile (Figure 3.8). (a) L = 4m (b) L = 16m Figure 3.18 Chart horizontal displacement, bending moment of pile length L= 4m (a); L = 16m(b) Survey results consistent with the results calculated by Matlock and Reese (1956); Zavriev (1962); Broms (1964) for short pile and long pile. However, according to the author's methodology, just a computer program can get results directly consider both short pile and long pile without sorting through piles step short, long pile; the single assumption simplification in calculations. 3.6.2 Single pile examination on hard rock layer The survey of single pile by reinforced concrete by section (30x30) cm, length l = 6m with elastic modulus Ec = of 30,000 MPa, Poisson's ratio ν c = 0.25. Pile effects of horizontal force P = 100kN at the poles. Calculation in two cases: uniform piles in the background; foot pile is plugged in tight limestone with 0.6 m thickness. (a) (b) Figure 3.19 Horizontal displacement diagram (a), bending moment (b) of the pile in uniform elastic foundation and is located in the elastic, hard truth. 17 Comment: when pile leans on the hard rock, horizontal displacement at pile foot is zero and turning point does not appear near the pile foot (Figure 3.19a), and bending moment in piles near the foot up from the pile in uniform background (Figure 3.19b). Thus, the method of comparison used Gauss system can also calculate similar pile against pile (pile tip is contraindicated in hard rock) 3.8 Conclusion of chapter 3 1 - By using method of comparative system has been built by the bending beam deflection equation under transverse shear strain. 2 - Develop a fully interactive problem between the pile and soil. Therefore, no need to add additional links contained soil piles at the edge, so this approach not only ensures the boundary conditions on the surface of the soil containing piles but also ensure the boundary conditions at infinity, boundary conditions between pile and soil. 3 - The problems were solved with the use of the Kelvin and Mindlin solution as comparative system to solve the problem on the horizontal force placed at the top of the pile, the pile foot or in different depths in the range including outside piles. Thereby the construction of pile load problem will be research content in the next chapter of the thesis. 4 - Results of the problem are compared with the results of a number of traditional methods helps improve the reliability of calculation theories. Chapter 4 STUDY INTERACTION PROBLEM BETWEEN PILE AND FOUNDATION UNDER HORIZONTAL LOAD AND SEISMIC LOAD 4.1 Solution of pulse units of infinite elastic space 4.2 Hysteretic damping coefficient of soil materials In the calculation of construction as well as foundation always consider the energy consumption in the fluctuations process and energy consumption which is described by viscous drag. Viscous drag is by viscous drag coefficient multiple with velocity. Under the loads, the foundation may appear deformed plastic, but plastic deformation does not depend on the frequency of the load, so this time instead of the usual viscous drag coefficient, people often [...]... the pile group, multiple spectral acceleration, etc, to receive generic and more accurate conclusions, as the basis for the calculation of seismic resistant design for pile foundation of construction 4.7 Conclusion of Chapter 4 1 - Develop the problem of interactive dynamics of the pile and the soil under horizontal dynamic loads and dynamic loads of soil that is allowed to mention the boundary conditions... problem, dynamic interaction between piles with soil under static load, horizontal dynamic load at any position Using the finite element method with soil tobe 3-D element, 20 buttons; pile with 2 button elements for displacement and 3 buttons for shear force to solve This method automatically satisfies the boundary conditions at infinity, the soil boundary conditions as well as conditions between pile and. .. add extra links as springs, viscous box on the contact surface between the pile ground, on the border of the land mass storing piles In addition, the study also investigates the parameters affecting the work of the piles such as pile length, stiffness of the pile, piles placed on a layer of hard rock and the effect of piles to the work of the soil 3 In the calculation of the construction dynamics and. .. waves in the horizontal plane and in vertical plane Based on the numerical solution of the finite element method, phenomenon of fluctuating surface amplification vertically to the wave propagation is studied, consistently with the theory of Love wave propagation 5 Formulating the interactive dynamic problem of the pile under the seismic load Using the convolution integral Duhamel provides solution in the... location of head and butt of pile (a) Horizontal displacement diagram (b), shear force (c), moment (d) following the length of pile at time 3,12s - Survey the tremor of time of earthquake T = 10,24s - The case mentions hysteretic damping coefficient (a) (b) Figure 4.16 The displacement diagram with frequency at the location of head of pile, middle and butt of pile (a) Horizontal displacement diagram following... of earthquake T = 5.12s (a) (b) (c) (d) Figure 4.19 Horizontal displacement diagram over time at the location of head and butt of pile (a) Horizontal displacement diagram (b), shear force (c), moment (d) following the length of pile at time 8,24s - Survey the tremor of time of earthquake T = 20,48s (a) (b) (a) (b) 24 (b) 25 (d) Figure 4.20 Horizontal displacement diagram over time at the location of... the constraint conditions as well as static interaction problem between the pile and foundation is presented in Chapter 3 When solving the problems in stable domain, meaning that there are no initial conditions Using the method of finite element is similar to the problem of static interactions between piles and foundation are presented in chapter 3 In this section, 20 nodes rectangular elements are... domain, then Fourier transform, reverse 26 (IFFT) provides results in the time domain Use the acceleration map of a real earthquake (El Centro 1940) as input parameters to survey, identify displacement parameters, torque, shear force of the pile at any time 6 Based on Matlab programming language to develop the calculating software program to serve the case studies and surveys: Mstatic1; Kstatic1; MstaticP1;... soil under horizontal dynamic loads and dynamic loads of soil, author received some main results as follows: 1 Through numerical solution using the finite element method, Kelvin solution about Mindlin solution can be provided, that means getting the solution of elastic infinite semi-space from the solution of elastic infinite space with load placed in any position 2 Formulating the static interactive... are used, each node has three displacements (u, v, w) in three axes x, y, z An element has a volume of 1 and is splitted into 8 nodes at the corner 4.5 Survey fluctuations of the pile of soil and piles under horizontal load 4.5.1 Survey Fluctuation of the pile of soil - The case does not mention hysteretic damping coefficient (a) (b) Figure 4.7 The chart of horizontal displacement of surface layers, the . formulate the interaction problem between pile and foundation under static load, dynamic horizontal load and seismic load with full consideration of the boundary conditions and endless radiation. BETWEEN PILE AND FOUNDATION UNDER HORIZONTAL LOAD Basing on the annalysis of study methods on the interaction between pile and foundation under horizontal load, we can draw some comments as follows:. of pile. 4. Study content Study stress-strain state of the soil mass under horizontal static load. Study static interaction problem between pile and foundation under horizontal static load.