Unified elastoplastic finite difference and its application
|
|
- Reynard Wiggins
- 6 years ago
- Views:
Transcription
1 Appl. Math. Mech. -Engl. Ed., 344, DOI /s c Shanghai University and Springer-Verlag Berlin Heidelberg 013 Applied Mathematics and Mechanics English Edition Unified elastoplastic finite difference and its application Zong-yuan MA 1, Hong-jian LIAO, Fa-ning DANG 1 1. School of Civil Engineering and Architecture, Xi an University of Technology, 5 South Jinhua Road, Xi an , P. R. China;. Department of Civil Engineering, Xi an Jiaotong University, 8 Xianning West Road, Xi an , P. R. China Abstract Two elastoplastic constitutive models based on the unified strength theory UST are established and implemented in an explicit finite difference code, fast Lagrangian analysis of continua FLAC/FLAC3D, which includes an associated/nonassociated flow rule, strain-hardening/softening, and solutions of singularities. Those two constitutive models are appropriate for metallic and strength-different SD materials, respectively. Two verification examples are used to compare the computation results and test data using the two-dimensional finite difference code FLAC and the finite element code ANSYS, and the two constitutive models proposed in this paper are verified. Two application examples, the large deformation of a prismatic bar and the strain-softening behavior of soft rock under a complex stress state, are analyzed using the three-dimensional code FLAC3D. The two new elastoplastic constitutive models proposed in this paper can be used in bearing capacity evaluation or stability analysis of structures built of metallic or SD materials. The effect of the intermediate principal stress on metallic or SD material structures under complex stress states, including large deformation, three-dimensional and non-association problems, can be analyzed easily using the two constitutive models proposed in this paper. Key words elastoplastic constitutive model, unified strength theory, explicit finite difference, effect of intermediate principal stress Chinese Library Classification O Mathematics Subject Classification 74C05 1 Introduction As one of the most important constitutive relations, strength theory yield or strength criterion was implemented in various nonlinear computer codes based on the finite element method FEM and the finite difference method FDM. Strength theory also plays an important role in material nonlinearity analysis, e.g., bearing capacity evaluation for structures. Metallic materials e.g., steel, copper, and aluminum have equal strength both in tension and compression, and hydrostatic stress has little influence on the strength of metallic materials. Nonmetallic Received Sept. 14, 01 / Revised Nov. 0, 01 Project supported by the National Natural Science Foundation of China No and the Central Financial Funds for the Development of Characteristic Key Disciplines in Local Universities Nos X101 and 106-5X105 Corresponding author Zong-yuan MA, Ph.D., mzy gogo@hotmail.com
2 458 Zong-yuan MA, Hong-jian LIAO, and Fa-ning DANG materials e.g., concrete, rock, and soil have different strengths in tension and compression i.e., the strength-different SD effect, and hydrostatic stress has a marked influence on the strength of SD materials. Plasticity theory can be divided into two components for strength theory, classical plasticity for metallic materials and generalized plasticity for SD materials. Three forms of work-hardening/softening or strain-hardening/softening law are used in plasticity theory, i.e., the isotropic hardening/softening case, the kinematic hardening/softening case, and the anisotropic hardening/softening case. The anisotropic hardening/softening law is basically discussed in generalized plasticity for SD materials. Many complex stress tests, such as the combined tension-torsion test and the true triaxial test, have verified that the intermediate principal stress has big influence on the mechanical behavior of metallic materials and SD materials. Lode [1] performed combined-stress experiments on thin-walled tubes made of iron, copper, and nickel. Taylor and Quinney [] did experiments on mild steel, copper, and aluminum thin-walled tubes that were subjected to tension and torsion. The effect of the intermediate principal stress was verified for all of the tested materials. Ivey [3] and Mair and Pugh [4] performed experiments with combined tension and torsion stress experiments of thin-walled tubes made of aluminum and pure copper and demonstrated the effect of the intermediate principal stress. Michelis [5] and Mogi [6] obtained experimental results from true triaxial testing of many rock samples and showed that the effect of the intermediate principal stress must be considered for rock. The results of complex stress tests of concrete indicate that the intermediate principal stress also has a marked effect on the strength of concrete [7 8]. The results from true triaxial tests show that the intermediate principal stress has a marked effect on the characteristics of soil under three different principal stresses [9 13]. Tresca and von Mises are two widely used yield criteria for metallic materials with the same yield stress in both tension and compression. The Mohr-Coulomb criterion is a widely used strength criterion for nonmetallic materials concrete or geomaterials, and it considers the SD and hydrostatic stress effects. However, the effect of the intermediate principal stress is not taken into account by the Tresca and the Mohr-Coulomb criterion, and some metals strengths exceed the yield surface of the von Mises criterion under a complex stress state [3 4]. Several strength criteria that take the intermediate principal stress and nonlinear yield surface into account, such as the William and Warnke [7], Lade and Duncan [14], and Matsuoka and Sun [15] criteria, which are proposed for geomaterials. Moreover, a smooth model for the Mohr-Coulomb criterion was established by Zienkiewicz and Pande [16], using the shape function suggested by Gudehus and Argyris, which used in many computer programs. However, those nonlinear yield surface criteria are appropriated for particular materials only. For example, the William-Warnke criterion is appropriate for concrete, and the Lade-Duncan and Matsuoka-Nakai criteria are both appropriate for sand. In addition, nonlinear yield surface criteria are difficult to implement in computer programming or to apply in theoretical analysis. The advances in strength theories of materials under complex stress states in the 0th century have been summarized by Yu [17]. Theories and methods.1 Unified strength theory UST The twin-shear stress element and the multi-shear stresses element, proposed by Yu [18], are shown in Fig. 1g and Fig.1h. The stress elements shown in Fig.1a, Fig.1b, and Fig. 1c are principal stress element, maximum shear stress element, and octahedral shear stress element, respectively. The multi-shear element see Fig. 1d is a rhombic dodecahedral multiple slip element differing from that of the three dimensional principal stress used in common continuum mechanics. There are three sets of principal shear stresses τ 13, τ 1, and τ 3 and normal stresses σ 13, σ 1, and σ 3 acting on the same sections. The maximum principal shear stress τ 13 equals the sum of the other two, i.e., τ 13 = τ 1 +τ 3. Thus, there are only two independent components among the three principal shear stresses. If the states of the three principal stresses change,
3 Unified elastoplastic finite difference and its application 459 e.g., if the intermediate principal stress σ increases from the value of the minimum principal stress σ 3 to that of the maximum principal stress σ 1, the intermediate principal shear stress τ 1 will be decreased, and τ 3 will be the intermediate principal shear stress. The Tresca and Mohr-Coulomb criteria consider only the maximum principal shear stress τ 13. Thus, the Tresca and Mohr-Coulomb criteria are both single-shear strength criteria and do not take into account the effect of the intermediate principal stress. The UST which takes the influence of the intermediate principal shear stress and the intermediate principal stress into account with a bilinear function, was proposed by Yu [19]. The expression of the UST for metallic materials equal strength in tension and compression can be formulated as follows: { f = τ13 + bτ 1 C, τ 1 bτ 3, f = τ 13 + bτ 3 C, τ 1 bτ 3. 1 Fig. 1 Various polyhedral element action by different stresses The expression of the UST for SD materials different strength in tension and compression can be formulated as follows: { f = τ13 + bτ 1 + βσ 13 + bσ 1 C, τ 1 + βσ 1 bτ 3 + βσ 3, f = τ 13 + bτ 3 + βσ 13 + bσ 3 C, τ 1 + βσ 1 bτ 3 + βσ 3, where β is a coefficient reflecting the effect of the normal stress on the strength of materials, b is a coefficient reflecting the effect of the intermediate principal shear stress τ 1 or τ 3 on the strength of materials, and σ 13, σ 1, and σ 3 are the normal stresses corresponding to τ 13, τ 1, and τ 3 that act on the sides of a multi-shear stress element. The stresses τ 13, τ 1, τ 3, σ 13, σ 1, and σ 3 can be obtained from the principal stresses σ 1, σ, and σ 3 as follows: τ 13 = σ 1 σ 3 σ 13 = σ 1+σ 3, τ 1 = σ 1 σ, σ 1 = σ 1 + σ, τ 3 = σ σ 3,, σ 3 = σ + σ 3. 3
4 460 Zong-yuan MA, Hong-jian LIAO, and Fa-ning DANG Therefore, the expressions of the UST for metallic materials and SD materials can be written in the forms of principal stresses as follows: f = σ b bσ + σ 3 σ s when σ 1 σ 1 + σ 3, f = b σ 1 + bσ σ 3 σ s when σ 1 σ 1 + σ 3, f = bσ + σ bn ϕ σ 1 + c f = σ 3 σ 1 + bσ + c N ϕ 1 + b when σ 1 + sin ϕ Nϕ when σ 1 + sinϕ Nϕ σ sin ϕ σ 3, σ sinϕ σ 3, where σ s is the tensile yield stress for metal, c and ϕ are the cohesion and friction angle for SD material, and N ϕ = 1 + sinϕ/1 sinϕ. The order of the three principal stresses follows σ 1 σ σ 3. For the non-associated flow rule [0], the plastic potential function g can be written as follows: g = bσ + σ 3 σ 1 when σ 1 + sin ϕ σ sin ϕ σ 3, 1 + bn ψ g = σ 3 σ 1 + bσ N ψ 1 + b when σ 1 + sin ϕ σ sin ϕ σ 3, where N ψ = 1+sin ψ 1 sin ψ, and ψ is the dilation angle for SD material. If ψ < ϕ, then the plastic flow rule is non-associated, while for the associated flow rule, ψ = ϕ. Direct comparisons with the analytical method can be made when the associated flow rule is used in numerical methods. We can see that the effect of the intermediate principal stress is an inherent property for materials. The different effects of the intermediate principal stress for different materials can be reflected by the coefficient b in the UST. The effect of the intermediate principal stress on the strength of materials is increased with increasing values of the coefficient b. A series of convex strength criteria can be obtained when the value of the coefficient b varies in the range of 0 b 1, and more materials can be described. The UST offers the possibility to choose a reasonable strength criterion for studies or applications. In addition, the UST can also introduce a family of non-convex strength criteria when b < 0 or b > 1. The limit loci of the UST on the deviatoric plane or principal stress space are shown in Fig., where θ b is the stress angle for the junction of two yield surfaces and depends only on the friction angle of the material. The yield surfaces of the UST cover the entire region of convex theory from the lower bound UST b = 0.0 to the upper bound UST b = 1.0. The Mohr-Coulomb strength criterion is the special case of the UST when the coefficient b is equal to zero. Other nonlinear yield surface criteria can be encompassed or linearly approximated by the UST when 0.0 < b < 1.0. Comparisons between the yield surfaces of the UST and test data for metallic materials [1 4], granite [1], and sand [10] are plotted in Fig. 3. The results shown in Fig. 3 indicate that the intermediate principal stress has certain influences on the strengths of different metallic materials or geomaterials under complex stress states, and the potential strengths of different materials can be predicted by the UST with various values of the parameter b. Yu et al. [1] proposed an elastoplastic associated and non-associated constitutive model based on the UST for SD materials and implemented it in a two-dimensional finite element program. The objective of this study is to establish two elastoplastic constitutive models for metal and SD materials based on the UST and implement those two models in the finite difference code fast Lagrangian analysis of continua FLAC/FLAC3D. Those two new elastoplastic constitutive models for 4 5 6
5 Unified elastoplastic finite difference and its application 461 FLAC/FLAC3D are programed easily in computer code and can be widely used in material nonlinearity analysis under complex stress states for both metallic and SD materials. Fig. Limit loci of UST on deviatoric plane Fig. 3 Effects of intermediate principal stress for different materials. Finite difference and its elastoplastic constitutive model FLAC/FLAC3D is a two/three-dimensional finite difference code with an explicit computation scheme. The behavior of structures built of metallic materials or SD materials that may undergo plastic flow when their yield limits are reached can be simulated easily by
6 46 Zong-yuan MA, Hong-jian LIAO, and Fa-ning DANG FLAC/FLAC3D [ 3]. An approach of mixed discretization is applied in FLAC/FLAC3D, as described by Marti and Cundall [4]. The components of the strain-rate tensor for a tetrahedral element are expressed as follows: ε ij = 1 6V 4 l=1 v l i nl j + v l j nl i S l. 7 The critical time step corresponding to the explicit finite difference scheme is given by the following expression [5] : t = T π = m k, 8 where T is the period of the system, k is the stiffness between two nodes, and m is the point mass. An important factor for the convergence of FLAC/FLAC3D is the out-of-balance force of the whole system. The out-of-balance force approaches zero when the medium reaches equilibrium or steady flow, and the calculation process of FLAC/FLAC3D is converged. Two elastoplastic constitutive models based on the UST, which are appropriate for metallic and SD materials, respectively, are written in C++ as user-written constitutive models and compiled as dynamic link library DLL files that can be loaded into the FLAC or FLAC3D code. The direction of the plastic strain increment is not uniquely on the additional corners of bilinear yield surfaces in the UST. Thus, corner singularity treatment must be processed. The method of dividing the stress space is used to address the problem of corner singularity in FLAC/FLAC3D. Figure 4 shows the stress space of the bilinear yield surfaces divided by a function σ = σ1+σ3 or σ = 1+sin ϕ σ sin ϕ σ 3 on the deviatoric stress plane. In this paper, the vector summation method is used to eliminate the corner singularity on the corner of two yield surfaces in the UST [6 8]. Fig. 4 Stress space dividing and corner singularity treatment for UST on deviatoric plane The elastic relations between elastic principal strain increments and principal stress increments in the constitutive model for FLAC/FLAC3D are as follows: σ i = S i ε e n, i = 1, n, 9 where S i is a linear function of the elastic principal strain increments ε e n. The plastic principal strain increments can be written as follows: ε p i = λ g σ i, 10
7 Unified elastoplastic finite difference and its application 463 where λ is a non-negative multiplier if plastic loading occurs. The expression of the elastoplastic constitutive model for FLAC/FLAC3D can be formulated as follows: g σi N = σi I λ S i, 11 σ n where λ S i g σ n is the plastic principal strain component in the matrix of the constitutive model, g is the plastic potential function, σi I are the stress components obtained from the incremental elastic law, and σi N are the new stress components obtained from the plastic flow rule. The plastic factor λ can be written as follows: λ = fσ I n f S n g σ n f0. 1 After substituting the expressions 4, 5, and 6 into 10 and 11, the elastoplastic constitutive model based on the UST can be obtained. For the metallic materials, the term λs i f σ n is expressed as follows: When σ 1 σ 1 + σ 3, S 1 λ f, λ f σ, λ f = λα α 1, S λ f, λ f σ, λ f = λb 1 + b α 1 α, S 3 λ f, λ f σ, λ f = λ 1 + b α 1 α. 13a When σ 1 σ 1 + σ 3, S 1 λ f, λ f σ, λ f = λ 1 + b α α 1, S λ f, λ f σ, λ f = λb 1 + b α α 1, S 3 λ f, λ f σ, λ f = λα 1 α, 13b where α 1 = K+ 4G 3, α = K G 3, K is the elastic bulk modulus, G is the elastic shear modulus, and the plastic factor λ is expressed as follows: λ = fσ1 I, σi, σi 3 f S g 1 σ n, g S σ n, g when σ 1 S3 σ n f0 σ 1 + σ 3, f σ1, I σ, I σ3 I f S g 1 σ n, g S σ n, g S3 σ n f 0 when σ 1 σ 1 + σ 3. For SD materials, the term λs i f σ n is expressed as follows: 14
8 464 Zong-yuan MA, Hong-jian LIAO, and Fa-ning DANG When σ < 1+sin ϕ σ sin ϕ σ 3, S 1 λ g, λ g, λ g α = λ α 1, σ N ψ S λ g, λ g, λ g bα 1 1 = λ + α 1, σ 1 + bn ψ 1 + bn ψ S 3 λ g, λ g, λ g α 1 b = λ + α 1. σ 1 + bn ψ 1 + bn ψ 15a When σ > 1+sin ϕ σ sin ϕ σ 3, S 1 λ g, λ g, λ g = λ σ α 1 N ψ b α 1, 1 + b 1 + b S λ g, λ g, λ g 1 = λ α 1 bα 1, σ N ψ 1 + b 1 + b S 3 λ g, λ g, λ g 1 = λ α 1 α, σ N ψ where the plastic factor λ is expressed as follows: fσ1 I, σi, σi 3 λ = f S 1 g, g σ, g, S g, g σ, g, S3 g, g σ, g f0 when σ < 1 + sin ϕ σ sin ϕ σ 3, f σ I 1, σi, σi 3 f S 1 g, g σ, g, S g, g σ, g, S3 g, g σ, g f 0 when σ > 1 + sin ϕ σ sin ϕ σ 3. 15b In the stage of strain-hardening/softening, the shear hardening parameter e ps is used to characterize the plastic shear strain, whose incremental form is defined as follows: e ps = 1 εps 1 εps m + ε ps εps m + ε ps 3 εps m 1 16, 17 where ε ps 1, εps, and εps 3 are the three plastic principal strains of the shear strength envelope, and ε ps m = ε ps 1 + εps + εps 3 /3. It is observed from 17 that the hardening parameter eps equals the square root of the second invariant of the strain. For the plastic hardening/softening stage, the yield function f for metallic materials can be written as follows: fσ n, σ s e ps = 0, belongs to the isotropic hardening/softening case of classical plasticity theory. The yield function f for SD materials can be written as follows: fσ n, ce ps, ϕe ps = 0, belongs to the anisotropic hardening/softening case, if the value of the strength parameter ϕ is varied with the hardening parameter e ps. 19 changes to the isotropic hardening/softening case, if only the value of the strength parameter c is varied with the hardening parameter e ps.
9 3 Verification problems Unified elastoplastic finite difference and its application Perforated plate in tension An experiment on a thin perforated plate in tension was produced by Theocaris and Marketos [9]. The material they tested was aluminum alloy with linear strain-hardening characteristics. The yield function f for linear strain-hardening can be written as follows: fσ n, σ s + He ps = 0, 0 where H is the plastic modulus for linear strain-hardening material, and the isotropic hardening law is used to analyze this problem. The material properties for the perforated plate were reported by Theocaris and Marketos as follows: Young s modulus E = MPa, Poisson s ratio ν = 0., the yield stress σ s = 43 MPa, and the linear hardening modulus H= 40 MPa. This problem is analyzed by Zienkiewicz [30] using FEM with the variable stiffness iterative method, assuming plane stress conditions and applying the von Mises criterion. In this paper, the same problem is analyzed using FLAC with the UST and ANSYS with the von Mises criterion. The FD or FE meshes, with dimensions and boundary conditions are given in Fig. 5. A quadrilateral element with linear interpolation Plane4 and the variable stiffness iterative method are used in ANSYS. A tensile load σ is applied to the top of the perforated plate, and the tensile load σ is divided into nine load increments a single load increment is 0.1 σ s. The relationship of average tensile stress σ mean along a section AB versus maximum tensile strain ε yy at the first yield point is shown in Fig. 6. The results from FLAC based on the UST b = 0.5 and b = and from FEM using the von Mises criterion show excellent agreement. Comparison of the tensile stress σ yy along AB at the load σ/σ s = obtained by numerical simulation and experiment is shown in Fig. 7. A satisfactory agreement is achieved among the tensile stress distributions obtained from FDM and FEM with the von Mises criterion. A slight discrepancy exists between the experimental data and the model predictions in the plastic region. Fig. 5 FD and FE meshes with boundary conditions for perforated plate plane stress
10 466 Zong-yuan MA, Hong-jian LIAO, and Fa-ning DANG Fig. 6 Stress-strain relationship for perforated plate Fig. 7 Stress distribution at root of plate with load σ/σ s = Bearing capacity of strip footing An example was given by Humpheson and Naylor [31] to analyze the influence of different forms of strength criteria i.e., Mohr-Coulomb, Drucker-Prager, William-Warnke, and Gudehus- Argyris on the bearing capacity of a strip footing using FEM with 3 parabolic elements. This example was examined again by Zienkiewicz and Pande [16]. Some differences in the values of the ultimate bearing capacity obtained using different strength criteria were detected for a plane strain and flexible footing on a weightless SD material based on the associated flow rule ψ = ϕ. The same problem is analyzed using the two-dimensional explicit finite difference code FLAC with the UST. The foundation soil is considered to be a linearly elasticperfectly plastic and homogeneous material with the following properties: Young s modulus E= 40.0 MPa, Poisson s ratio ν = 0., the cohesion c = 10.0 kpa, the friction angle ϕ = 0, and the dilation angle ψ = ϕ. The uniform vertical pressure is gradually increased the incremental pressure is 0. MPa until the algorithm could not converge. The ultimate bearing capacity is reached when the numerical nonconvergence occurres. Figure 8 shows the limit loci of the UST, Drucker-Prager, and nonlinear yield surface criteria William-Warnke and Gudehus- Argyris on the deviatoric plane. Figure 9 shows the FD meshes with the dimensions and boundary conditions used in analysis. The values of the ultimate bearing capacity P u /c calculated by Humpheson and Naylor and obtained using FLAC are presented in Table 1. Figure 10 shows a comparison of the values of ultimate bearing capacity P u /c obtained from FLAC using the UST with those calculated by Humpheson and Naylor. The values of P u /c calculated by Humpheson and Naylor with the William-Warnke and Gudehus-Argyris criteria lie between the values obtained using the UST for b = 0.5 and b = 0.5, and the values of P u /c calculated by Humpheson and Naylor with the Gudehus-Argyris criterion are close to the UST solution for b = Application problems 4.1 Large deformation and non-association problem Large deformation is a common phenomenon for metallic or nonmetallic structures under static or dynamic loading. There are two ways to describe the movement in continuum mechanics, the Lagrangian description and the Eulerian description [30]. In the Lagrangian formulation, which is often applied in solid mechanics, the nodes of the Lagrangian mesh move together with the material. Large deformations may lead to problematic mesh and large element distortions
11 Unified elastoplastic finite difference and its application 467 in Lagrangian simulations. In the Eulerian formulation, an Eulerian reference mesh is needed to trace the motion of the material in the Eulerian domain. This is a field description that is often applied in fluid mechanics. Fig. 8 Limit loci of UST and other criteria on deviatoric plane Fig. 9 FD meshes with boundary conditions for bearing capacity of strip footing plane strain Table 1 Ultimate bearing capacity P u/c values from FLAC and Humpheson & Naylor P u/c Yield criteria Humpheson & Naylor FEM FLAC FDM Drucker-Prager Mohr-Coulomb/UST b = Drucker-Prager extension cone UST b = Williams-Warnke 1.1 Gudehus-Argyris.6 UST b = UST b = Drucker-Prager compression cone
12 468 Zong-yuan MA, Hong-jian LIAO, and Fa-ning DANG Fig. 10 Comparison bearing capacity P u/c for flexible strip footing values for different criteria The updated Lagrangian method is used for large deformation analysis in FLAC/FLAC- 3D [ 3]. Incremental displacements are added to the coordinates so that the grid moves and deforms with the material it represents. Considering the large displacements, displacement gradients, and rotations, the node coordinates and element stresses are updated after each step calculation, and the rotation components must be taken into account in calculation of the element stresses. The updates of node coordinates and element stresses are expressed as follows: x l i t + t = x l i t + tv l i t + t, 1 σ ij = σ ij + σ C ij, σ C ij = ω ik σ kj σ ik ω kj t, where x l i and v l i are the components of coordinate and velocity for nodes in large deformation analysis and σij C is the stress correction for displacement gradients and rotations and not taken into consideration in small deformation mode. The fundamental hypothesis of associated plasticity, which is applicable in plasticity theory for metallic materials, does not hold for some SD materials, e.g., concrete, rock, or soil [0]. The plastic volume dilation of SD materials is much higher than that predicted by the associative flow rule, and the plastic volume dilation can be reduced by the non-associative flow rule. Considering the influence of large deformations and the non-associative flow rule, an example of a prismatic bar with a square cross section subject to axial compression load is simulated in FLAC3D. The updated Lagrangian formulation is used to analyze the large deformation of the prismatic bar in compression. The material of the prismatic bar has the following properties, Young s modulus E = 50.0 MPa, Poisson s ratio ν = 0., the cohesion c = 10.0 kpa, the friction angle ϕ = 30, the dilation angle ψ = ϕ, ψ = ϕ, and ψ = 0. The boundary conditions and dimensions used in the analysis are shown in Fig. 11. The relationship between the axial compression load and the axial displacement of the prismatic bar is shown in Fig. 1. For small deformations, the uniaxial stress state is closely approximated by the prismatic bar being subject to pure axial compression, and the intermediate principal stress has little influence on the strength and deformation of the prismatic bar. However, the uniaxial stress state is changed to a complex stress state when the deformation is large. As Fig.1 shows, a difference induced by the effect of the intermediate principal stress is detected from the results of the Mohr-Coulomb criterion and UST b = 1.0. The influence of the intermediate principal stress is decreased with decreasing values of the dilation angle for the non-association problem. However, little effect of the intermediate principal stress is indicated for small deformations.
13 Unified elastoplastic finite difference and its application 469 Fig. 11 Dimensions and boundary conditions for large deformation analysis three-dimensional Fig. 1 Relationship between axial compression load and axial displacement 4. Strength of soft rock under complex stress states The hump curve function for the friction angle ϕ is used to simulate the strain-softening behavior of soft rock in this study. The hump curve function is shown as follows: ϕe ps = eps H + Re ps H + Pe ps, 3 where the slope of the hump curve is controlled by the parameter H, and the peak and residual values of the hump curve are controlled by the parameters P and R, and the anisotropic hardening/softening law is used to analyze this problem. A conventional triaxial shear test on diatomaceous soft rock was conducted in Refs. [3 33]. The strength parameter c for diatomaceous soft rock varies slightly with increasing value of ε 1, but the strength parameter ϕ has a peak and residual value. Figure 13 shows the variation in friction angle for diatomaceous soft rock with the axial strain increased. Based on the variation in friction angle determined with the conventional triaxial test, the values of the parameters H, P, and R in 3 can be determined. A hexahedral element with unite length is used to simulate a cubical specimen of diatomaceous soft rock, and the conventional triaxial and true triaxial tests are simulated using the
14 470 Zong-yuan MA, Hong-jian LIAO, and Fa-ning DANG three-dimensional finite difference code FLAC3D. The confining pressure σ 3 is applied by the stress boundary condition on all sides of the cubical specimen element to generate the consolidated stress. For the conventional triaxial test simulation simple stress state, the three principal stresses follow σ 1 σ = σ 3, a vertical velocity load m per step is applied to the top of the cubical specimen element to simulate the maximum principal stress σ 1, and the confining pressure σ 3 is also applied on both sides of the cubical specimen element. The bottom of the cubical specimen element is fixed in the vertical direction, so that both ends of the cubical specimen element are rigid and smooth. For the true triaxial test simulation complex stress state, the three principal stresses follow σ 1 σ σ 3 and the maximum principal stress σ 1 and confining pressure σ 3 also depend on the velocity and the stress boundary condition. The intermediate principal stress σ depends on the stress boundary condition according to the principal stress ratio of σ σ 3 /σ 1 σ 3. The magnitude of σ is controlled by the program to observe the principal stress ratio of σ σ 3 /σ 1 σ 3 during the true triaxial test simulation. The technology of the true triaxial test for a cubical specimen based on the principal stress ratio of σ σ 3 /σ 1 σ 3 is widely used in complex stress tests for geomaterials [9 10]. The values of the material property parameters used in the numerical analysis are presented in Table. Because the cubical specimen of diatomaceous soft rock has little volume dilation after the conventional triaxial test, the dilation angle ψ is set to be zero non-associated flow rule throughout the whole calculation process for the triaxial and true triaxial test simulations. Figure 14 shows the relationships of stress σ 1 σ 3 versus axial strain ε 1 for diatomaceous soft rock as measured in the triaxial test and as predicted by FLAC3D. The comparison between the model prediction and the measured results shows that the stress-strain relationship of diatomaceous soft rock can be simulated by the strain-hardening/softening model with the hump curve function. Fig. 13 Variation of friction angle ϕ with axial strain ε 1 for diatomaceous soft rock Fig. 14 Relationship of stress σ 1 σ 3 to axial strain ε 1 for diatomaceous soft rock measured by conventional triaxial test or computed by FLAC3D Table Material parameters used in true triaxial tests simulation of diatomaceous soft rock Name Diatomaceous soft rock Linearly elastic-perfectly plastic Strain softening Elastic modulus Poisson s ratio Cohesion E/MPa ν c/kpa P H R
15 Unified elastoplastic finite difference and its application 471 The strain-hardening/softening constitutive model based on the UST is then used to simulate a true triaxial test for diatomaceous soft rock. The hardening/softening rule for the friction angle ϕ follows. The stress-strain relationships for diatomaceous soft rock predicted by the true triaxial test simulation for various values of the principal stress ratio σ σ 3 /σ 1 σ 3 are shown in Fig.15. The results shown in Fig.15 indicate that the peak strength q p and the residual strength q s for diatomaceous soft rock under a complex stress state are higher than under a simple stress state. The stress-strain curve obtained using the Mohr-Coulomb criterion is insensitive to the variations in the intermediate principal stress. The peak strength q p and residual strength q s values obtained from the true triaxial test simulation for various values of the principal stress ratio σ σ 3 /σ 1 σ 3 are shown in Fig.16. The peak strength q p and residual strength q s first increase and then decrease with increasing values of the principal stress ratio σ σ 3 /σ 1 σ 3. The values of q p and q s under the stress state of σ = σ 1 σ σ 3 /σ 1 σ 3 = 1.0 are lower than those under the stress state of σ = σ 3 σ σ 3 /σ 1 σ 3 = 0.0. Fig. 15 Stress-strain curves yield by true triaxial test simulation under different principal stress ratios
16 47 Zong-yuan MA, Hong-jian LIAO, and Fa-ning DANG Fig. 16 Variation of peak and residual strength from true triaxial test simulation with principal stress ratio of σ σ 3/σ 1 σ 3 5 Conclusions Two new elastoplastic constitutive models, appropriate for metallic materials and SD materials, are presented and implemented in the explicit finite difference code FLAC/FLAC3D. Two verification examples are used to compare test data and computation results. Good agreement is achieved between the test results and model predictions, and the two constitutive models proposed in this paper are verified. Two application examples, a prismatic bar subject to axial compression and soft rock under a complex stress state, are analyzed. Based on the results of this study, the following conclusions can be drawn: i The two elastoplastic constitutive models proposed in this paper, based on the unified strength theory with bilinear expression and yield surfaces, can be easily implemented and used in the explicit finite difference code FLAC/FLAC3D. ii These two elastoplastic constitutive models are verified using two classical examples and may be widely used in bearing capacity evaluation or stability analysis of structures built of metallic or SD materials in situations, where the effect of the intermediate principal stress must be taken into account. iii The results from the analysis of the large deformation and non-association problem indicate that the effect of the intermediate principal stress on large deformation problems is related to the plastic volume dilation of SD materials. The true triaxial test simulation results for diatomaceous soft rock suggest that the intermediate principal stress has a significant influence on the strength and strain-softening behavior of diatomaceous soft rock under a complex stress state. iv Considering the entire region between the lower bound using either the Tresca or the Mohr-Coulomb criterion and the upper bound UST b = 1.0 of the convex theory, the effect of the intermediate principal stress on metallic or SD materials or structures under complex stress states, including large deformation, three-dimensional or non-association problems, can be easily analyzed using the two new elastoplastic constitutive models proposed in this paper. Acknowledgements The first author expresses gratitude to Prof. Mao-hong YU of Xi an Jiaotong University, People s Republic of China. Discussion with Prof. YU is of great help and value.
17 Unified elastoplastic finite difference and its application 473 References [1] Lode, W. Versuche ueber den Einlfuss der mittleren Hauptspannung auf dsa fliessen der metals eisen kupfer und nickel. Z. Physik., 36, [] Taylor, G. I. and Quinney, H. The plastic distortion of metals. Phil. Trans. Roy. Soc., 30, [3] Ivey, H. J. Plastic stress-strain relations and yield surfaces for aluminium alloys. J. Mech. Eng. Sci., 3, [4] Mair, W. M. and Pugh, H. L. D. Effect of pre-strain on yield surfaces in copper. J. Mech. Eng. Sci., 6, [5] Michelis, P. Polyaxial yielding of granular rock. J. Eng. Mech., 1118, [6] Mogi, K. Effect of the intermediate principal stress on rock failure. J. Geophys. Res., 70, [7] William, K. J. and Warnke, E. P. Constitutive model for the triaxial behavior of concrete. Proceedings of Seminar on Concrete Structures Subjected to Triaxial Stresses, [8] Du, X. L., Lu, D. C., Gong, Q. M., and Zhao, M. Nonlinear unified strength criterion for concrete under three-dimensional stress states. J. Eng. Mech., 1361, [9] Green, G. E. and Bishop, A. W. A note on the drained strength of sand under generalized strain conditions. Geotechnique, 191, [10] Lade, P. V. and Duncan, J. M. Cubical triaxial tests on cohesionless soil. J. Soil Mech. Found. Div., 9910, [11] Ochiai, H. and Lade, P. V. Three-dimensional behavior of sand with anisotropic fabric. Journal of Geotechnical Engineering, 10910, [1] Matsuoka, H. and Nakai, T. Stress-deformation and strength characteristics of soil under three difference principal stress. Proceeding of Japan Society of Civil Engineers, 3, [13] Xing, Y. C., Liu, Z. D., and Zheng, Y. R. A failure criterion of loess in Chinese. Journal of Hydraulic Engineering, 1, [14] Lade, P. V. and Duncan, J. M. Elastoplastic stress-strain theory for cohesionless soil. J. Geo. Eng. Div., 10110, [15] Matsuoka, H. and Sun, D. A. Extension of spatially mobilized plane SMP to frictional and cohesive materials and its application to cemented sands. Soils and Foundations, 354, [16] Zienkiewicz, O. C. and Pande, G. N. Some Useful Forms of Isotropic Yield Surfaces for Soil and Rock Mechanics, Wiley, New York, [17] Yu, M. H. Advances in strength theories for materials under complex stress state in the 0th Century. Appl. Mech. Rev., 553, [18] Yu, M. H. Twin shear stress yield criterion. Int. J. Mech. Sci., 51, [19] Yu, M. H. Unified Strength Theory and Its Applications, Springer, Berlin 004 [0] Vermeer, P. A. and de Borst, R. Non-associated plasticity for soils concrete and rock. Heron, 93, [1] Yu, M. H., Yang, S. Y., Fan, S. C., and Ma, G. W. Unified elasto-plastic associated and nonassociated constitutive model and its engineering applications. Comput. Struct., 716, [] Itasca Consulting Group. FLAC-Fast Lagrangian Analysis of Continua, User s manual, Minneapolis 005 [3] Itasca Consulting Group. FLAC 3D-Fast Lagrangian Analysis of Continua in 3 Dimensions, User s manual, Minneapolis 005 [4] Marti, J. and Cundall, P. A. Mixed discretization procedure for accurate modelling of plastic collapse. Int. J. Numer. Anal. Meth. Geomech., 6, [5] Bathe, K. J. and Wilson, E. L. Numerical Methods in Finite Element Analysis, Prentice-Hall, New Jersey 1976
18 474 Zong-yuan MA, Hong-jian LIAO, and Fa-ning DANG [6] Koiter, W. T. Stress-strain relations, uniqueness and variational theorems for elastic-plastic materials with singular yield surface. Quart. Appl. Math., 11, [7] Sloan, S. W. and Booker, J. R. Removal of singularities in Tresca and Mohr-Coulomb yield functions. Commun. Appl. Numer. Mech.,, [8] Yu, H. S. A close-form solution of stiffness matrix for Tresca and Mohr-Coulomb plasticity models. Comput. Struct., 533, [9] Theocaris, P. S. and Marketos, E. Elastic-plastic analysis of perforated thin strips of a strainhardening material. J. Mech. Phys. Solids, 1, [30] Zienkiewicz, O. C. The Finite Element Method in Engineering Science, McGraw-Hill, New York 1971 [31] Humpheson, C. and Naylor, D. J. The Importance of the Form of the Failure Criterion, University of Wales, Swansea 1975 [3] Liao, H. J., Su, L. J., Pu, W. C., and Yin, J. H. Test and numerical analysis of the constitutive relation of a diatomaceous soft rock. Marine Georesources & Geotechnology, 1, [33] Liao, H. J., Pu, W. C., and Yin, J. H. Numerical modeling of the effect of strain rate on the stress-strain relation for soft rock using a 3-D elastic visco-plastic model. Int. J. Rock. Mech. Min., 413,
A return mapping algorithm for unified strength theory model
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2015; 104:749 766 Published online 10 June 2015 in Wiley Online Library (wileyonlinelibrary.com)..4937 A return mapping
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 informationPlasticity R. Chandramouli Associate Dean-Research SASTRA University, Thanjavur
Plasticity R. Chandramouli Associate Dean-Research SASTRA University, Thanjavur-613 401 Joint Initiative of IITs and IISc Funded by MHRD Page 1 of 9 Table of Contents 1. Plasticity:... 3 1.1 Plastic Deformation,
More informationINFLUENCE OF NONASSOCIATIVITY ON THE BEARING CAPACITY
INFLUENCE OF NONASSOCIATIVITY ON THE BEARING CAPACITY OF A STRIP FOOTING By Jian-Hua Yin, 1 Yu-Jie Wang, and A. P. S. Selvadurai 3 ABSTRACT: This paper examines the ultimate bearing capacity of a strip
More informationTrue Triaxial Tests and Strength Characteristics Study on Silty Sand Liang MA and Ping HU
217 2 nd International Conference on Test, Measurement and Computational Method (TMCM 217) ISBN: 978-1-6595-465- True Triaxial Tests and Strength Characteristics Study on Silty Sand Liang MA and Ping HU
More informationNUMERICAL ANALYSIS OF PASSIVE EARTH PRESSURES WITH INTERFACES
III European Conference on Computational Mechanics Solids, Structures and Coupled Problems in Engineering C.A. Mota Soares et.al. (eds.) Lisbon, Portugal, 5-8 June 2006 NUMERICAL ANALYSIS OF PASSIVE EARTH
More informationCONSIDERATIONS CONCERNING YIELD CRITERIA INSENSITIVE TO HYDROSTATIC PRESSURE
CONSIDERATIONS CONCERNING YIELD CRITERIA INSENSITIVE TO HYDROSTATIC PRESSURE ADRIAN SANDOVICI, PAUL-DORU BARSANESCU Abstract. For distinguishing between pressure insensitive and pressure sensitive criteria,
More informationClassical fracture and failure hypotheses
: Chapter 2 Classical fracture and failure hypotheses In this chapter, a brief outline on classical fracture and failure hypotheses for materials under static loading will be given. The word classical
More informationUniversity of Sheffield The development of finite elements for 3D structural analysis in fire
The development of finite elements for 3D structural analysis in fire Chaoming Yu, I. W. Burgess, Z. Huang, R. J. Plank Department of Civil and Structural Engineering StiFF 05/09/2006 3D composite structures
More informationGEO E1050 Finite Element Method Mohr-Coulomb and other constitutive models. Wojciech Sołowski
GEO E050 Finite Element Method Mohr-Coulomb and other constitutive models Wojciech Sołowski To learn today. Reminder elasticity 2. Elastic perfectly plastic theory: concept 3. Specific elastic-perfectly
More informationNonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield Surface
Transaction A: Civil Engineering Vol. 16, No. 6, pp. 512{519 c Sharif University of Technology, December 2009 Research Note Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield
More informationEnhancing Prediction Accuracy In Sift Theory
18 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS Enhancing Prediction Accuracy In Sift Theory J. Wang 1 *, W. K. Chiu 1 Defence Science and Technology Organisation, Fishermans Bend, Australia, Department
More informationFINITE ELEMNT ANALYSIS FOR EVALUATION OF SLOPE STABILITY INDUCED BY CUTTING
FINITE ELEMNT ANALYSIS FOR EVALUATION OF SLOPE STABILITY INDUCED BY CUTTING Toshinori SAKAI Department of Environmental Science and Technology, Mie University, Tsu, Japan Tadatsugu TANAKA Graduate School
More informationDiscrete Element Modelling of a Reinforced Concrete Structure
Discrete Element Modelling of a Reinforced Concrete Structure S. Hentz, L. Daudeville, F.-V. Donzé Laboratoire Sols, Solides, Structures, Domaine Universitaire, BP 38041 Grenoble Cedex 9 France sebastian.hentz@inpg.fr
More informationFLAC3D analysis on soil moving through piles
University of Wollongong Research Online Faculty of Engineering - Papers (Archive) Faculty of Engineering and Information Sciences 211 FLAC3D analysis on soil moving through piles E H. Ghee Griffith University
More informationMODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Uniaxial Model: Strain-Driven Format of Elastoplasticity
MODELING OF CONCRETE MATERIALS AND STRUCTURES Kaspar Willam University of Colorado at Boulder Class Meeting #3: Elastoplastic Concrete Models Uniaxial Model: Strain-Driven Format of Elastoplasticity Triaxial
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 informationGeneral Strength Criterion for Geomaterials
TECHNICAL NOTES General Strength Criterion for Geomaterials M. D. Liu 1 and J. P. Carter 2 Abstract: This note presents the results of a study on the variation of the strength of geomaterials. A general
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 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 informationEstimation of the Passive Earth Pressure with Inclined Cohesive Backfills: the Effect of Intermediate Principal Stress is Considered
108 The Open Mechanical Engineering Journal, 011, 5, 108-116 Open Access Estimation of the Passive Earth Pressure with Inclined Cohesive Backfills: the Effect of Intermediate Principal Stress is Considered
More informationEffect of the intermediate principal stress on fault strike and dip - theoretical analysis and experimental verification
Effect of the intermediate principal stress on fault strike and dip - theoretical analysis and experimental verification B. Haimson University of Wisconsin, USA J. Rudnicki Northwestern University, USA
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 informationStress and Strains in Soil and Rock. Hsin-yu Shan Dept. of Civil Engineering National Chiao Tung University
Stress and Strains in Soil and Rock Hsin-yu Shan Dept. of Civil Engineering National Chiao Tung University Stress and Strain ε 1 1 2 ε 2 ε Dimension 1 2 0 ε ε ε 0 1 2 ε 1 1 2 ε 2 ε Plane Strain = 0 1 2
More informationPrediction of torsion shear tests based on results from triaxial compression tests
Prediction of torsion shear tests based on results from triaxial compression tests P.L. Smith 1 and N. Jones *2 1 Catholic University of America, Washington, USA 2 Geo, Lyngby, Denmark * Corresponding
More informationNumerical Assessment of the Influence of End Conditions on. Constitutive Behavior of Geomaterials
Numerical Assessment of the Influence of End Conditions on Constitutive Behavior of Geomaterials Boris Jeremić 1 and Zhaohui Yang 2 and Stein Sture 3 ABSTRACT In this paper we investigate the behavior
More informationUnified Constitutive Model for Engineering- Pavement Materials and Computer Applications. University of Illinois 12 February 2009
Unified Constitutive Model for Engineering- Pavement Materials and Computer Applications Chandrakant S. Desai Kent Distinguished i Lecture University of Illinois 12 February 2009 Participation in Pavements.
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 informationFEM for elastic-plastic problems
FEM for elastic-plastic problems Jerzy Pamin e-mail: JPamin@L5.pk.edu.pl With thanks to: P. Mika, A. Winnicki, A. Wosatko TNO DIANA http://www.tnodiana.com FEAP http://www.ce.berkeley.edu/feap Lecture
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 information2 CONSTITUTIVE MODELS: THEORY AND IMPLEMENTATION
CONSTITUTIVE MODELS: THEORY AND IMPLEMENTATION 2-1 2 CONSTITUTIVE MODELS: THEORY AND IMPLEMENTATION 2.1 Introduction There are twelve basic constitutive models provided in, arranged into null, elastic
More informationA circular tunnel in a Mohr-Coulomb medium with an overlying fault
MAP3D VERIFICATION EXAMPLE 9 A circular tunnel in a Mohr-Coulomb medium with an overlying fault 1 Description This example involves calculating the stresses and displacements on a fault overlying a 5 m
More informationLecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity
Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity by Borja Erice and Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling
More informationMODELING GEOMATERIALS ACROSS SCALES
MODELING GEOMATERIALS ACROSS SCALES JOSÉ E. ANDRADE DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING AFOSR WORKSHOP ON PARTICULATE MECHANICS JANUARY 2008 COLLABORATORS: DR XUXIN TU AND MR KIRK ELLISON
More informationMODELING GEOMATERIALS ACROSS SCALES JOSÉ E. ANDRADE DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING EPS SEMINAR SERIES MARCH 2008
MODELING GEOMATERIALS ACROSS SCALES JOSÉ E. ANDRADE DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING EPS SEMINAR SERIES MARCH 2008 COLLABORATORS: DR XUXIN TU AND MR KIRK ELLISON THE ROADMAP MOTIVATION
More informationConstitutive models: Incremental plasticity Drücker s postulate
Constitutive models: Incremental plasticity Drücker s postulate if consistency condition associated plastic law, associated plasticity - plastic flow law associated with the limit (loading) surface Prager
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 informationSOIL MODELS: SAFETY FACTORS AND SETTLEMENTS
PERIODICA POLYTECHNICA SER. CIV. ENG. VOL. 48, NO. 1 2, PP. 53 63 (2004) SOIL MODELS: SAFETY FACTORS AND SETTLEMENTS Gabriella VARGA and Zoltán CZAP Geotechnical Department Budapest University of Technology
More informationOn the mechanical behaviour of dry cohesionless soils by DEM simulations
On the mechanical behaviour of dry cohesionless soils by DEM simulations Dorival M. Pedroso, Sergio Galindo-Torres & David J. Williams The University of Queensland, Brisbane, Australia This paper presents
More informationStress Integration for the Drucker-Prager Material Model Without Hardening Using the Incremental Plasticity Theory
Journal of the Serbian Society for Computational Mechanics / Vol. / No., 008 / pp. 80-89 UDC: 59.74:004.0 Stress Integration for the Drucker-rager Material Model Without Hardening Using the Incremental
More informationABSTRACT INTRODUCTION
Optimization of soil anchorages M. Teschner, C. Mattheck Kernforschungszentrum Karlsruhe GmbH, Institut fur Materialforschung II, W-7500 Karlsruhe 1, Postfach 3640, Germany ABSTRACT A new optimization
More informationA Simple and Accurate Elastoplastic Model Dependent on the Third Invariant and Applied to a Wide Range of Stress Triaxiality
A Simple and Accurate Elastoplastic Model Dependent on the Third Invariant and Applied to a Wide Range of Stress Triaxiality Lucival Malcher Department of Mechanical Engineering Faculty of Tecnology, University
More informationMODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam
MODELING OF CONCRETE MATERIALS AND STRUCTURES Class Meeting #1: Fundamentals Kaspar Willam University of Colorado at Boulder Notation: Direct and indicial tensor formulations Fundamentals: Stress and Strain
More informationA simple elastoplastic model for soils and soft rocks
A simple elastoplastic model for soils and soft rocks A SIMPLE ELASTO-PLASTIC MODEL FOR SOILS AND SOFT ROCKS by Roberto Nova Milan University of Technology 1. MODEL HISTORY The model is the result of the
More informationStresses Analysis of Petroleum Pipe Finite Element under Internal Pressure
ISSN : 48-96, Vol. 6, Issue 8, ( Part -4 August 06, pp.3-38 RESEARCH ARTICLE Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure Dr.Ragbe.M.Abdusslam Eng. Khaled.S.Bagar ABSTRACT
More informationAuthor(s) Okajima, Kenji; Tanaka, Tadatsugu; Symposium on Backwards Problem in G.
Title Backwards Analysis for Retaining Wa based upon ateral Wall Displacemen Author(s) Okajima, Kenji; Tanaka, Tadatsugu; Proceeding of TC302 Symposium Osaka Citation Symposium on Backwards Problem in
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 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 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 informationNumerical Modeling of Interface Between Soil and Pile to Account for Loss of Contact during Seismic Excitation
Numerical Modeling of Interface Between Soil and Pile to Account for Loss of Contact during Seismic Excitation P. Sushma Ph D Scholar, Earthquake Engineering Research Center, IIIT Hyderabad, Gachbowli,
More informationMathematical Relations Related to the Lode. Parameter for Studies of Ductility
Advanced Studies in Theoretical Physics Vol. 0, 06, no., - 4 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/astp.06.50 Mathematical Relations Related to the Lode Parameter for Studies of Ductility
More informationTheory of Shear Strength
MAJ 1013 ADVANCED SOIL MECHANICS Theory of Shear Strength Prepared by, Dr. Hetty 1 Strength of different materials Steel Concrete Soil Tensile strength Compressive strength Shear strength Complex behavior
More informationEffect Of The In-Situ Stress Field On Casing Failure *
Effect Of The In-Situ Stress Field On Casing Failure * Tang Bo Southwest Petroleum Institute, People's Republic of China Lian Zhanghua Southwest Petroleum Institute, People's Republic of China Abstract
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationApplication of Three Dimensional Failure Criteria on High-Porosity Chalk
, 5-6 May 00, Trondheim, Norway Nordic Energy Research Programme Norwegian U. of Science and Technology Application of Three Dimensional Failure Criteria on High-Porosity Chalk Roar Egil Flatebø and Rasmus
More informationTheory of Shear Strength
SKAA 1713 SOIL MECHANICS Theory of Shear Strength Prepared by, Dr. Hetty 1 SOIL STRENGTH DEFINITION Shear strength of a soil is the maximum internal resistance to applied shearing forces The maximum or
More informationUSER S MANUAL. 1D Seismic Site Response Analysis Example. University of California: San Diego.
USER S MANUAL 1D Seismic Site Response Analysis Example http://www.soilquake.net/ucsdsoilmodels/ University of California: San Diego August 2, 2017 Table of Contents USER'S MANUAL TABLE OF CONTENTS Page
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 informationCONTENTS. Lecture 1 Introduction. Lecture 2 Physical Testing. Lecture 3 Constitutive Models
CONTENTS Lecture 1 Introduction Introduction.......................................... L1.2 Classical and Modern Design Approaches................... L1.3 Some Cases for Numerical (Finite Element) Analysis..........
More informationLaboratory Testing Total & Effective Stress Analysis
SKAA 1713 SOIL MECHANICS Laboratory Testing Total & Effective Stress Analysis Prepared by: Dr. Hetty Mohr Coulomb failure criterion with Mohr circle of stress 2 ' 2 ' ' ' 3 ' 1 ' 3 ' 1 Cot Sin c ' ' 2
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 informationModule 5: Failure Criteria of Rock and Rock masses. Contents Hydrostatic compression Deviatoric compression
FAILURE CRITERIA OF ROCK AND ROCK MASSES Contents 5.1 Failure in rocks 5.1.1 Hydrostatic compression 5.1.2 Deviatoric compression 5.1.3 Effect of confining pressure 5.2 Failure modes in rocks 5.3 Complete
More informationSeismic stability safety evaluation of gravity dam with shear strength reduction method
Water Science and Engineering, 2009, 2(2): 52-60 doi:10.3882/j.issn.1674-2370.2009.02.006 http://kkb.hhu.edu.cn e-mail: wse@hhu.edu.cn Seismic stability safety evaluation of gravity dam with shear strength
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 informationELASTOPLASTICITY THEORY by V. A. Lubarda
ELASTOPLASTICITY THEORY by V. A. Lubarda Contents Preface xiii Part 1. ELEMENTS OF CONTINUUM MECHANICS 1 Chapter 1. TENSOR PRELIMINARIES 3 1.1. Vectors 3 1.2. Second-Order Tensors 4 1.3. Eigenvalues and
More informationPLAXIS 3D TUNNEL. Material Models Manual version 2
PLAXIS 3D TUNNEL Material Models Manual version 2 TABLE OF CONTENTS TABLE OF CONTENTS 1 Introduction...1-1 1.1 On the use of different models...1-1 1.2 Limitations...1-2 2 Preliminaries on material modelling...2-1
More information1 Introduction. Abstract
Abstract This paper presents a three-dimensional numerical model for analysing via finite element method (FEM) the mechanized tunneling in urban areas. The numerical model is meant to represent the typical
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 information13 Plastic Flow in a Punch Problem
Plastic Flow in a Punch Problem 13-1 13 Plastic Flow in a Punch Problem 13.1 Problem Statement Difficulties are sometimes reported in the modeling of plastic flow where large velocity gradients exist.
More informationModule 4 Lecture 20 Pore water pressure and shear strength - 4 Topics
Module 4 Lecture 20 Pore water pressure and shear strength - 4 Topics 1.2.6 Curvature of the Failure Envelope Effect of angularity of soil particles Effect of rate of loading during the test 1.2.7 Shear
More informationParticle flow simulation of sand under biaxial test
5th International Conference on Civil Engineering and Transportation (ICCET 2015) Particle flow simulation of sand under biaxial test Xiao-li Dong1,2, a *,Wei-hua Zhang1,a 1 Beijing City University, China
More informationTHE BEHAVIOUR OF REINFORCED CONCRETE AS DEPICTED IN FINITE ELEMENT ANALYSIS.
THE BEHAVIOUR OF REINFORCED CONCRETE AS DEPICTED IN FINITE ELEMENT ANALYSIS. THE CASE OF A TERRACE UNIT. John N Karadelis 1. INTRODUCTION. Aim to replicate the behaviour of reinforced concrete in a multi-scale
More informationProceedings of the ASME th International Conference on Ocean, Offshore and Arctic Engineering OMAE2016 June 19-24, 2016, Busan, South Korea
Proceedings of the ASME 26 35th International Conference on Ocean, Offshore and Arctic Engineering OMAE26 June 9-24, 26, Busan, South Korea OMAE26-54554 LOCAL STRAIN AND STRESS CALCULATION METHODS OF IRREGULAR
More informationFinite Element analysis of Laterally Loaded Piles on Sloping Ground
Indian Geotechnical Journal, 41(3), 2011, 155-161 Technical Note Finite Element analysis of Laterally Loaded Piles on Sloping Ground K. Muthukkumaran 1 and N. Almas Begum 2 Key words Lateral load, finite
More informationReference material Reference books: Y.C. Fung, "Foundations of Solid Mechanics", Prentice Hall R. Hill, "The mathematical theory of plasticity",
Reference material Reference books: Y.C. Fung, "Foundations of Solid Mechanics", Prentice Hall R. Hill, "The mathematical theory of plasticity", Oxford University Press, Oxford. J. Lubliner, "Plasticity
More informationNUMERICAL MODELING OF INSTABILITIES IN SAND
NUMERICAL MODELING OF INSTABILITIES IN SAND KIRK ELLISON March 14, 2008 Advisor: Jose Andrade Masters Defense Outline of Presentation Randomized porosity in FEM simulations Liquefaction in FEM simulations
More informationOn generalization of constitutive models from two dimensions to three dimensions
University of Wollongong Research Online Faculty of Engineering and Information Sciences - Papers: Part A Faculty of Engineering and Information Sciences 008 On generalization of constitutive models from
More informationA rate-dependent Hosford-Coulomb model for predicting ductile fracture at high strain rates
EPJ Web of Conferences 94, 01080 (2015) DOI: 10.1051/epjconf/20159401080 c Owned by the authors, published by EDP Sciences, 2015 A rate-dependent Hosford-Coulomb model for predicting ductile fracture at
More informationSimulation of footings under inclined loads using different constitutive models
Simulation of footings under inclined loads using different constitutive models J. Hintner, P.A. Vermeer Institute of Geotechnical Engineering, University of Stuttgart, Germany P.-A. von Wolffersdorff
More informationAnalysis of Controlling Parameters for Shear behavior of Rock Joints with FLAC3D
International Journal of Structural and Civil Engineering Research Vol. 5, No. 1, February 2016 Analysis of Controlling Parameters for Shear behavior of Rock Joints with FLAC3D Prasoon Tiwari and Hakan
More informationFE Prediction of Bearing Capacity of Reinforced Soil under Plane Strain
The First Pan American Geosynthetics Conference & Exhibition 2-5 March 28, Cancun, Mexico FE Prediction of Bearing Capacity of Reinforced Soil under Plane Strain Nogueira, C., Oliveira, R., Zornberg, J.G.,
More informationNumerical and experimental analyses for bearing capacity of rigid strip footing subjected to eccentric load
DOI: 10.1007/s11771-014-2386-5 Numerical and experimental analyses for bearing capacity of rigid strip footing subjected to eccentric load LU Liang( 卢谅 ) 1, 2, WANG Zong-jian( 王宗建 ) 3, K. ARAI 4 1. School
More informationProf. B V S Viswanadham, Department of Civil Engineering, IIT Bombay
56 Module 4: Lecture 7 on Stress-strain relationship and Shear strength of soils Contents Stress state, Mohr s circle analysis and Pole, Principal stressspace, Stress pathsin p-q space; Mohr-Coulomb failure
More informationStress and fabric in granular material
THEORETICAL & APPLIED MECHANICS LETTERS 3, 22 (23) Stress and fabric in granular material Ching S. Chang,, a) and Yang Liu 2 ) Department of Civil Engineering, University of Massachusetts Amherst, Massachusetts
More informationModule-4. Mechanical Properties of Metals
Module-4 Mechanical Properties of Metals Contents ) Elastic deformation and Plastic deformation ) Interpretation of tensile stress-strain curves 3) Yielding under multi-axial stress, Yield criteria, Macroscopic
More informationExpansion of circular tubes by rigid tubes as impact energy absorbers: experimental and theoretical investigation
Expansion of circular tubes by rigid tubes as impact energy absorbers: experimental and theoretical investigation M Shakeri, S Salehghaffari and R. Mirzaeifar Department of Mechanical Engineering, Amirkabir
More informationNumerical evaluation of the bearing capacity factor N of ring footings for associated soils
Numerical evaluation of the bearing capacity factor N of ring footings for associated soils BENMEBAREK Sadok Numerical Modeling and Instrumentation Laboratory, Biskra University, Algeria, benmebareks@yahoo.fr
More informationTHE INTERPRETATION OF LABORATORY SOIL TESTS
57ième CONGRÈS CANADIEN DE GÉOTECHNIQUE 5ième CONGRÈS CONJOINT SCG/AIH-CNN 57TH CANADIAN GEOTECHNICAL CONFERENCE 5TH JOINT CGS/IAH-CNC CONFERENCE THE INTERPRETATION OF LABORATORY SOIL TESTS P. Guo Department
More informationON THE FACE STABILITY OF TUNNELS IN WEAK ROCKS
33 rd 33 Annual rd Annual General General Conference conference of the Canadian of the Canadian Society for Society Civil Engineering for Civil Engineering 33 e Congrès général annuel de la Société canadienne
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 informationElastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount
Elastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount Min Ding* 1, and Jian Lin 2 1 MIT/WHOI Joint Program, 2 Woods Hole Oceanographic Institution *Woods Hole Oceanographic
More informationEDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d
Institute of Industrial Science, University of Tokyo Bulletin of ERS, No. 48 (5) A TWO-PHASE SIMPLIFIED COLLAPSE ANALYSIS OF RC BUILDINGS PHASE : SPRING NETWORK PHASE Shanthanu RAJASEKHARAN, Muneyoshi
More informationME 2570 MECHANICS OF MATERIALS
ME 2570 MECHANICS OF MATERIALS Chapter III. Mechanical Properties of Materials 1 Tension and Compression Test The strength of a material depends on its ability to sustain a load without undue deformation
More informationExample-3. Title. Description. Cylindrical Hole in an Infinite Mohr-Coulomb Medium
Example-3 Title Cylindrical Hole in an Infinite Mohr-Coulomb Medium Description The problem concerns the determination of stresses and displacements for the case of a cylindrical hole in an infinite elasto-plastic
More informationEffect of embedment depth and stress anisotropy on expansion and contraction of cylindrical cavities
Effect of embedment depth and stress anisotropy on expansion and contraction of cylindrical cavities Hany El Naggar, Ph.D., P. Eng. and M. Hesham El Naggar, Ph.D., P. Eng. Department of Civil Engineering
More informationPLASTICITY FOR CRUSHABLE GRANULAR MATERIALS VIA DEM
Plasticity for crushable granular materials via DEM XIII International Conference on Computational Plasticity. Fundamentals and Applications COMPLAS XIII E. Oñate, D.R.J. Owen, D. Peric and M. Chiumenti
More informationMicromechanics-based model for cement-treated clays
THEORETICAL & APPLIED MECHANICS LETTERS 3, 216 (213) Micromechanics-based model for cement-treated clays Zhenyu Yin, 1, a) 2, b) and Pierre Yves Hicher 1) Department of Civil Engineering, Shanghai Jiao
More informationGeology 229 Engineering Geology. Lecture 5. Engineering Properties of Rocks (West, Ch. 6)
Geology 229 Engineering Geology Lecture 5 Engineering Properties of Rocks (West, Ch. 6) Common mechanic properties: Density; Elastic properties: - elastic modulii Outline of this Lecture 1. Uniaxial rock
More informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationConstitutive Equations
Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module
More information