LINE SEARCH TECHNIQUES FOR ELASTO-PLASTIC FINITE ELEMENT COMPUTATIONS IN GEOMECHANICS
|
|
- Joel Beverly Jefferson
- 6 years ago
- Views:
Transcription
1 REPORT NO. UCD/CGM-01/04 CENTER FOR GEOTECHNICAL MODELING LINE SEARCH TECHNIQUES FOR ELASTO-PLASTIC FINITE ELEMENT COMPUTATIONS IN GEOMECHANICS BY B. JEREMIC DEPARTMENT OF CIVIL & ENVIRONMENTAL ENGINEERING COLLEGE OF ENGINEERING UNIVERSITY OF CALIFORNIA AT DAVIS SEPTEMBER 2001
2 University of Califorina at Davis, Center for Geotechnical Modeling Report UCD-CGM-00 Published in the Communications in Numerical Methods in Engineering, Vol. 17, issue 2, pages , Line Search Techniques for Elasto Plastic Finite Element Computations in Geomechanics Boris Jeremić Department of Civil and Environmental Engineering, University of California, Davis, CA 95616, phone , fax , Abstract In this paper we present a globally convergent modification of Newton's method for integrating constitutive equations in elasto plasticity of geomaterials. Newton's method is known to be q-quadratically convergent when the current solution approximation is adequate. Unfortunately, it is not unusual to expend significant computational time in order to achieve satisfactory results. We will present a technique which can be used when the Newton step is unsatisfactory. This scheme can be considered as a modified version of the traditional concept of backtracking along the Newton direction if a full Newton step provides unsatisfactory results. The method is also known as line search technique. The technique is applied to the fully implicit Newton algorithm for a hardening or softening general isotropic geomaterials at the constitutive level. Various solution details and visualizations are presented, which emerge from the realistic modeling of highly nonlinear constitutive behavior observed in the analysis of cohesionless granular materials. Key Words: Elasto plasticity, Geomechanics, Constitutive integrations, Newton iterative method, Line search technique 1
3 1 Introduction Problem of accurately following the equilibrium path in numerical modeling of geomaterials elasto plasticity has been researched for some time now. It has been shown (Runesson and Samuelson [12], Simo and Taylor [14]) that consistent use of Newton iterative algorithm on the global, finite element and local constitutive levels provides for fast convergence. However, complexities of elasto plastic constitutive models for geomaterials put high demand on the constitutive driver. In particular, low confinement region, usually found near the soil surface or during liquefaction behavior of sands, with highly nonlinear hardening and softening behavior, high dilatancy angles (non associative behavior) and highly curved yield surface, can lead to the numerical failure (divergence) of the constitutive driver. This leads to the interruption of finite element computations and consequently requires rerun of the problem with smaller loading steps. While the use of full Newton scheme improves rate of convergence, in some cases the iterative algorithm does not converge at all. The necessary condition for good convergence behavior of Newton methods is that current solution approximation is good enough, i.e. it is within the convergence region of Newton method. If the current solution approximation is not good enough, Newton method will not converge but rather bounce through solution space until erroneous iterations are interrupted. One of the strategies used to prevent failure of Newton methods is the line search technique. It should be noted that line search will not cure all of the problems associated with local, constitutive level Newton iterations. However, as long as the function (and their derivatives) used in iterations (yield function, potential function, hardening/softening functions) are analytic, line search will provide for global convergence in the sense defined by Dennis and Schnabel [3]. Equilibrium iterations for material nonlinear finite element computations can be in general separated in two levels. First iteration level is tied to the constitutive, elasto plastic computations. On this level, constitutive driver, for a given strain increment 1 iterates in stress and internal variable space until convergence criteria is met. On the second level, nonlinear finite element system of equations solver is iterating until balance of internal and external forces is achieved. Global level iteration to achieve balance of internal and external forces have seen use of line search techniques (eg. Crisfield [2], Larsson et al. [10], Simo and Meschke [13]). In their recent paper, Dutko et al. [4], 1 Assuming displacement based implementation of FEM. 2
4 used a variant of line search algorithm for constitutive level iterations. They applied this technique to the biaxial anisotropic yield criterion (Barlat et al. [1]) and presented interesting and useful results. Of particular interest are statistics on number of line search iterations for a sharply curved region of the yield line. Unfortunately, the actual line search technique used was just briefly described. The paper is organized as follows. Section 2 briefly describes elastic plastic algorithmic formulation based on the fully implicit, Newton procedure and the B Material Model used in computations. Section 3 describes theory behind the line search technique and section 4 describes application of the line search techniques to the constitutive integration problems. Section 5 presents numerical examples that illustrate described developments. It is worthwhile noting that while the presented techniques is applied to small deformation elasto plastic problems it has been used in the large deformation elasto plastic algorithms as well. However, inherent complexities of LDEP algorithms developments might hide the basic ideas of using line search techniques so we restrict our presentation to small deformation elasto plasticity. Although developments presented here are simplified to small deformation format, the generality oftheapproach is not lost. 2 Elastic Plastic Geomechanics In this section we briefly present Backward Euler algorithms for the solutions of elastic plastic problems in geomechanics. To this end, we use the additive decomposition of the strain increment into elastic and plastic parts together with flow theory of plasticity and Karush Kuhn Tucker conditions to formulate the elastic plastic problem. Detailed description of the formulations can be found elsewhere (eg. Jeremić and Sture [8]). The fully implicit Backward Euler algorithm is developed in general stress tensor internal variable tensor setting and is based on the following equations ff ij = E ijkl ffl kl ffl p kl ; ffl p ij = n ffl p ij + m ij q Λ = n q Λ + h Λ ; F =0 (1) where, ffl ij ; ffl e ij and ffl p ij are the total, elastic and plastic strain tensors respectively, ff ij is the Cauchy stress tensor, q Λ represents suitable set of internal variables, h Λ is the plastic moduli and F is the yield surface function at the final position. The asterisk 3
5 in the place of indices in q Λ replaces n indices, so that for example in the case of isotropic hardening q Λ is a scalar while for kinematic hardening q Λ is second order tensor. Equations (1), are the nonlinear algebraic equations to be solved for the unknowns ff ij, ffl p ij, q Λ and. Newton iterativescheme is used to solve for the single vector return in stress tensor internal variable space. We can circumvent the problem of finding the solution in the softening region, which, in general stress tensor internal variable space belongs to non convex space, by defining the tensor of stress residuals. Then, by working on the problem of minimizing the stress residual we convert the non convex problem to the convex one. Of course, as usual, by using Newton iterative method, we have to provide for continuation of iterative procedure even if the Newton step is not satisfactory. That is precisely the point where the line search technique shows it usefulness. The Backward Euler algorithm is based on the elastic predictor plastic corrector strategy: ff ij = pred ff ij E ijkl m kl (2) where pred ff ij = E ijkl ffl kl is the elastic trial stress state and m kl = (@Q=@ff kl )j is the gradient to the plastic potential function in stress space at the final position. We define a tensor of residuals r ij = ff ij pred ff ij E ijkl m kl that represents the difference between the current stress state ff ij (3) and the Backward Euler stress state pred ff ij E ijkl m kl. By using first order Taylor series expansion of the tensor of residuals r ij and the first order Taylor expansion of the yield function F, and after some algebraic manipulations we obtain the change in consistency parameter d = F old n mn old r ij T 1 ijmn n mn E ijkl H kl T 1 ijmn ο Λ h Λ (4) We have also introduced the fourth order tensors T ijmn and H ijmn : T ijmn = ffi im ffi nj + E mn fi ; H kl = m kl Λ fi h Λ (5) where n mn mn, ο Λ Λ and dq Λ = d h Λ (ff ij ;q Λ ). With the solutions for d we can write the iterative solution for dff mn and dq Λ, in the extended stress internal variable space as: ψ! dff mn = old F old n old mn r ij T 1 ijmn r ij + E n mn E ijkl H kl T 1 ijkl H kl T 1 ijmn ijmn (6) ο Λ h Λ 4
6 dq Λ = ψ! F old n old mn r ij T 1 ijmn h n mn E ijkl H kl T 1 Λ (7) ijmn ο Λ h Λ Iterative procedure is continued until the objective function kr ij k =0is satisfied given a certain tolerance. In order to start the Newton iterative procedure, we use forward Euler solution as a starting point start ff mn = E mnpq dffl pq E mnpq cross n rs E rstu dffl tu cross n ab E abcd cross m cd ο Λ h Λ cross m pq (8) start q Λ = previous q Λ + ψ cross n mn E mnpq dffl pq cros n mn E mnpq cros m pq ο Λ h Λ! h Λ (9) where cross () denotes the point where trial state crosses yield surface. Consistent, Newton iterations on the constitutive level are reflected on the global, finite element level through use of algorithmic tangent stiffness (ATS) tensor AT S E ep pqmn defined as dff pq = AT S E ep pqmn dffl mn with AT S E ep R pqkl H kl n ij R ijmn pqmn = R pqmn (10) n ot R otpq H pq + ο Λ h Λ where we have used reduced stiffness tensor R mnkl =( T ijmn ) 1 E ijkl. It is important to note that the line search technique described later does not alter the ATS tensor AT S E ep pqmn. A detailed derivation is given in Jeremić and Sture [6]. We use small deformation version of the B material model (Jeremić et al. [5]) for our computations. The model relies on the development behind the so called MRS Lade model (Sture et al. [16]). The B Model is a single surface model, with uncoupled cone portion and cap portion hardening. Very low confinement region was carefully modeled and the yield surface was shaped in such a way to mimic recent findings obtained during Micro Gravity Mechanics tests aboard Space Shuttle (Sture et al. [15]). Figure (1) depicts meridian and deviatoric traces and a full yield surface. It also depicts cone and cap hardening/softening functions. Detailed description of the model is given by Jeremić et al. [5]. 3 Line Search Technique In this section we develop theory behind the line search techniques. We base our developments on concept of minimization theory for a scalar function f(x Λ ). 5 We show
7 q η σ 1 low confinment trace a) b 1 η init η cone p t 0.71(p c - p t ) p c p σ 2 σ 3 b) high confinment trace c) η peak d) κ peak η start η residual κ cone e) p c Figure 1: (a) Meridian trace of yield, ultimate or potential surface. (b) Change of the deviatoric trace of yield surface changes along the mean stress axes. (c) Yield and/or potential surface in principal stress space. (d) Cone hardening function. (e) Cap hardening function. later that in our case, function to be minimized will be the energy norm of the vector of residuals r ij (Eq. 3). The minimization problem is defined as (eg. Dennis and Schnabel [3].): min x Λ2R n f(x Λ):R n!r (11) The basic idea of a global method for minimization is to take the step that lead downhill for the function f(x Λ ). One chooses a direction p Λ from the current point x c Λ in which f(x Λ ) decreases initially and a new point x + Λ in this direction from x c Λ is such that f(x + Λ ) < f(x c Λ). Such a direction p Λ is called a descent direction. From the mathematical point of view, p Λ is a descent direction from x c Λ if the directional derivative of f(x Λ ) at x c Λ in the direction p Λ is (x c Λ p Λ < 0 (12) 6 k cap
8 If (12) holds, then it is guaranteed that for a small positive, f(x c Λ + p Λ ) <f(x c Λ). The idea of line search algorithm can be described as follows: ffl At iteration k do: calculate a descent direction p k Λ, set x k+1 Λ ψ (x k Λ + k p k Λ) for some k that makes x k+1 Λ iterate. Figure (2) shows the basic concept: select x k+1 Λ an acceptable next by considering the half of a one dimensional cross section of f(x Λ ) in which f(x Λ ) decreases initially from x k Λ. f( x k + * ζ p k ) * Quadratic approximation slope d f( x k ) * d x * p * < 0 ζ=0 ζ Figure 2: A cross section of f(x Λ ) from x k Λ in the direction p k Λ. The term line search" refers to the procedure of choosing the acceptable k. The so called exact line search" accounts for finding the exact solution of the one dimensional minimization problem, i.e. finding the exact k so that f(x k Λ + k p k Λ) attains minimum. This was the preferred approach to the problem until mid 1960s. More careful computational testing has led to the use of slack line search" which has a weak acceptance criteria for k as a more computationally efficient procedure. The common procedure now is to try the full Newton step first (with k =1)and then, if k =1fails to satisfy criterion, to reduce k in a systematic way, along the direction defined by that step. Systematic reduction of k along the descent direction can be achieved by applying the line search techniques through backtracking algorithm: Given ff 2 (0; 1): 2 k =1; 7
9 while f(x k Λ + k p k Λ) >f(x k Λ)+ff Λ) p Λ Λ do: k ψ ρ k where, usually ρ = 1 2 ; x k+1 Λ ψ x k Λ + k p k Λ; Dennis and Schnabel [3] provide rather powerful convergence results for properly chosen steps. 4 Application to Elasto Plastic Computations in Geomechanics In this section we describe the application of line search technique to the elasto plastic constitutive level equilibrium iterations. The objective function that is followed is the Euclidean norm of tensor of residuals r ij. To this end we rewrite the elastic predictor plastic corrector Eq. (2) as: ff ij = pred ff ij E ijkl m kl (13) It is important to note that the addition of scalar line search parameter does not change the initial equations. Moreover, as is only used in line search improvement and is not a function of any state variable (stress, internal variables or displacements) first order Taylor expansions used to come up with the iterative steps (Eq. (4) (7) are not changed. The only difference is that now solution for the change in consistency parameter d reads: d = with changed tensors T ijmn and H kl T ijmn = ffi im ffi nj + E mn fi ; F old n mn old r ij T 1 ijmn n mn E ijkl H kl T 1 ijmn ο Λ h Λ (14) H kl = m kl Λ fi h Λ (15) With this changes the line search algorithm can now be specialized to elasto plastic implicit computations. In each iterative step k we do: ffl Given fi 2 (0; 1 2 ): ffl k =1;(Full Newton step) 8
10 ffl while ψ kr ij ( k ff ij + k d k ff ij ; k q Λ + k d k q Λ )k > kr ij ( k ff ij ; k q Λ )k + fi k fl k ff mn Λ do: k ψ ρ k where, usually ρ = 1 2 ; k+1 ff ij ψ k ff ij + k d k ff ij and k+1 q Λ ψ k q Λ + k d k q Λ k q Λ!fl We always start with k =1(that is, a full Newton step) and only if that iteration fail, we apply the backtrack algorithm. As we are following value of the Euclidean norm of tensor of residuals r ij, failure of Newton step is defined as divergence of that scalar variable. In other words, increase the value of Euclidean norm of tensor of residuals r ij means that our Newton step in the extended stress internal variable space is not valid. It also means that convergence of the Newton iterative procedure is questionable and that it will most probably fail in this iterative step. 5 Numerical Example Previous theoretical developments has been implemented in FEI (Finite Element Interpreter) finite element program. FEI is our experimental software platform, written in C++ with some Fortran modules. It is build on top of ndarray and FEMtools class libraries (Jeremić and Sture [9]). The line search algorithm is implemented on both global, finite element level as well as on the local, constitutive level, described in this paper. We follow nonlinear finite element iterations in analyzing directional shear cell (DSC) experiments done at the University of Colorado at Boulder. For example McFadden ([11]) used DSC apparatus in order to investigate shear behavior of various sands. Here we present modeling of particular plane strain experiment K 0 14:5 30 in figure (3). Loading is divided in two stages. The first stage comprises isotropic compression to target confinement state (in this case p = 180:0kPa). The second stage is a shear loading until instability occurs (instability in the experimental test, numerically we can follow the specimen beyond limit point). It is not clear if the instabilities in DSC experiments (loss of control over loading process) were due to the bifurcation phenomena inside the specimen or to the global rotation of the specimen. Since DSC is a load controlled device, only shear deformation of fl xy ' 3:5 %was reached in the laboratory experiment. 9
11 Shear stress τ [kn/m^2] volumetric strain [%] a) Shear deformation γ [%] b) Shear deformation γ [%] Figure 3: Numerical modeling of a shear test K 0 14:5 30 a) Shear stress axial strain curves. b) Volumetric strain axial strain curves. As the numerical experiment is proceeding, the length of the loading steps is controlled by using the variable hyperspherical arc length constraint ([7]). In this particular case, the Newton iterative algorithm failed at constitutive level after stress state advanced into elastic plastic regime. This part of the elastic plastic computations is sometimes problematic, since there is a sudden activation of the hardening mechanism. In this particular case, we have stopped the computations at the beginning of the problematic incremental step. Then, we resumed iterations by manually decreasing the step size from 0.5% to 0.1%. Figure 4(a) shows values of the norm of the tensor of residual kr ij k in such a successful iteration. It can be seen from Fig. 4(a) that convergence is quite fast, leading to the conclusion that the initial estimate of the stress and internal variable state was within the Newton convergence region. Similar to the previous numerical experiment, we resumed computations at the beginning of the problematic step, but this time with larger incremental deformation of 0.5% and with the line search algorithm turned off. Figure 4(b) shows values of the norm of the tensor of residual kr ij k for the first incremental step. We can observe that the initial, uncorrected iteration produces non converging set of values for the norm of the tensor of residual kr ij k. It is interesting to note that the very first iterative step takes the value of kr ij k to higher than initial values. Moreover, Figure 5(a) shows that although value of evaluated yield surface F initially point downward from the initial predictor, it is eventually diverging too. Finally, once again we resume computations at the beginning of the problematic incremental step, now with the line search algorithm turned on. It can be seen from 10
12 norm of residual r 1e-03 1e-04 1e-05 norm of residual r 1e-01 1e-02 1e-03 1e-04 1e-05 erratic iterations (line search off) first iterative step corrected iterations (line search on) a) 1e number of iterations b) 1e number of iterations Figure 4: (a) Residual norm kr ij k values for a successful iteration with small incremental steps (one step of 0.1% shear deformation). (b) Residual norm values kr ij k values for erratic and corrected iterations (one step of 0.5% shear deformation). Fig. 4(b) that not only the algorithm helps advance the iterations without any signs of divergence, it also produces a very good convergence rate. This, of course, is the result of our algorithm always first trying the full Newton step (for which we set k = 1). Figure 5(b) depicts the problematic, first iterative step from Figure 4(b) in some more details. We follow the value of the objective function kr ij k through 10 subincrements within the problematic iterative step. At this point it is important to remember that Newton iterative algorithm makes a quadratic approximation (in multidimensional space) of the objective function, here the Euclidean norm kr ij k. In this particular step, the approximation for kr ij k does a poor job, and although the objective function initially decreases, eventually its final value is higher than that it had at the initial state. It is interesting to note that at 5th subincrement, objective function kr ij k actually does attain a minimum value, but full Newton step (at 10th subincrement) leads to a higher than initial values for kr ij k. Figure (6) depicts failed and corrected iterations in the space of stress invariants p, q and p = 1 3 ff kk q = s s ijs ij cos 3 = 3p s 3 ijs jk s q ki ( 1 2 s ijs ij ) 3 s ij = ff ij 1 3 ff kkffi ij (16) In this particular view, we are looking at the B model yield surface from the tensile region (negative p). Divergent iteration is oscillating through the stress space and there is no sign of recovery. On the other hand, corrected oscillation, after modifying one 11
13 1e+04 1e-01 erratic iterations 1e+02 1e-02 yield function F 1e+00 1e-02 1e-04 correct iterations norm of residual r 1e-03 function quadratic approximation a) 1e number of iterations b) 1e subincremented iterative step Figure 5: (a) Yield surface values F for erratic and corrected iterations (one step of 0.5% shear deformation). (b) Dissection of failed iteration. iterative step by the line search algorithm, converges successfully. 6 Summary In this paper we havepresented a globally convergent modification of Newton's method for integrating constitutive equations in elasto plasticity of geomaterials. The method has been developed in rigorous mathematical framework and implemented in experimental finite element program FEI. The practical use of method was illustrated in details. The application of line search techniques in numerical modeling of elasto plasticity of geomaterials should provide for robust following of equilibrium path on both, global, finite element and local, constitutive levels. Acknowledgment The author gratefully acknowledge partial support by National Science Foundation's Pacific Earthquake Engineering Research Center Grant NSF-PEER
14 q predictor point φ correct iterations erratic iterations yield surface p Figure 6: Illustration of erratic and correct iterations in stress space. 13
15 References [1] Barlat, F., Lege, D. J., and Brem, J. C. A six component yield function for anisotropic materials. International Journal for Plasticity 7 (1991), [2] Crisfield, M. A. Non - Linear Finite Element Analysis of Solids and Structures Volume 1: Essentials. John Wiley and Sons, Inc. New York, [3] Dennis, Jr., J. E., and Schnabel, R. B. Numerical Methods for Unconstrained Optimization and Nonlinear Equations. Prentice Hall, Engelwood Cliffs, New Jersey , [4] Dutko, M., Perić, D., and Owen, D. R. J. Universal anisotropic yield criterion based on superquadratic functional representation: Part1. algorithmic issues and accuray analysis. Computer methods in applied mechanics and engineering 109 (1993), [5] Jeremić, B., Runesson, K., and Sture, S. A model for elastic plastic pressure sensitive materials subjected to large deformations. International Journal of Solids and Structures 36, 31/32 (1999), [6] Jeremić, B., and Sture, S. Implicit integration rules in plasticity: Theory and implementation. Report to NASA Marshall Space Flight Center, Contract: NAS , University of Colorado at Boulder, June [7] Jeremić, B., and Sture, S. Refined finite element analysis of geomaterials. In Proceedings of 11th Engineering Mechanics Conference (Fort Lauderdale, Florida, May 1996), Y. K. Lin and T. C. Su, Eds., Engineering Mechanics Division of the American Society of Civil Engineers, pp [8] Jeremić, B., and Sture, S. Implicit integrations in elasto plastic geotechnics. International Journal of Mechanics of Cohesive Frictional Materials 2 (1997), [9] Jeremić, B., and Sture, S. Tensor data objects in finite element programming. International Journal for Numerical Methods in Engineering 41 (1998), [10] Larsson, R., Runesson, K., and Sture, S. Embedded localization band in undrained soil based on regularized strong discontinuity theory and FE analysis. International Journal of Solids and Structures 33, (1996), [11] McFadden, J. J. Experimental response of sand during pricipal stress rotations. Master of Science thesis, University od Colorado at Boulder, December
16 [12] Runesson, K., and Samuelsson, A. Aspects on numerical techniques in small deformation plasticity. In NUMETA 85 Numerical Methods in Engineering, Theory and Applications (1985), G. N. P. J. Middleton, Ed., AA.Balkema., pp [13] Simo, J. C., and Meschke, G. A new class of algorithms for classical plasticity extended to finite strain. application to geomaterials. Computational Mechanics 11 (1993), [14] Simo, J. C., and Taylor, R. L. Consistent tangent operators for rateindependent elastoplasticity. Computer Methods in Applied Mechanics and Engineering 48 (1985), [15] Sture, S., Costes, N., Batiste, S., Lankton, M., AL-Shibli, K., Jeremić, B., Swanson, R., and Frank, M. Mechanics of granular materials at low effective stresses. ASCE Journal of Aerospace Engineering 11, 3 (July 1998), [16] Sture, S., Runesson, K., and Macari-Pasqualino, E. J. Analysis and calibration of a three invariant plasticity model for granular materials. Ingenieur Archiv 59 (1989),
Abstract The paper presents derivation of a fully implicit Newton algorithm for direct integration of constitutive equations, in extended stress inter
Published in: International Journal for Mechanics of Cohesive Frictional Materials and Structures Volume 2, Pages 165-183, 1997 Implicit Integrations in Elasto Plastic Geotechnics Boris Jeremić z Stein
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 informationPublished in the International Journal of Solids and Structures, vol. 36 No. 31/32 pages , A Model for Elastic Plastic Pressure Sensiti
Published in the International Journal of Solids and Structures, vol. 36 No. 31/32 pages 491-4918, 1999. A Model for Elastic Plastic Pressure Sensitive Materials Subjected to Large Deformations Boris Jeremić
More informationISSUES IN MATHEMATICAL MODELING OF STATIC AND DYNAMIC LIQUEFACTION AS A NON-LOCAL INSTABILITY PROBLEM. Ronaldo I. Borja Stanford University ABSTRACT
ISSUES IN MATHEMATICAL MODELING OF STATIC AND DYNAMIC LIQUEFACTION AS A NON-LOCAL INSTABILITY PROBLEM Ronaldo I. Borja Stanford University ABSTRACT The stress-strain behavior of a saturated loose sand
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 information3D and 2D Formulations of Incremental Stress-Strain Relations for Granular Soils
Archives of Hydro-Engineering and Environmental Mechanics Vol. 55 (2008), No. 1 2, pp. 45 5 IBW PAN, ISSN 121 726 Technical Note D and 2D Formulations of Incremental Stress-Strain Relations for Granular
More informationDEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS
DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS Mohsen Safaei, Wim De Waele Ghent University, Laboratory Soete, Belgium Abstract The present work relates to the
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 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 informationStochastic Elastic-Plastic Finite Element Method for Performance Risk Simulations
Stochastic Elastic-Plastic Finite Element Method for Performance Risk Simulations Boris Jeremić 1 Kallol Sett 2 1 University of California, Davis 2 University of Akron, Ohio ICASP Zürich, Switzerland August
More informationTransactions on Engineering Sciences vol 6, 1994 WIT Press, ISSN
Large strain FE-analyses of localized failure in snow C.H. Liu, G. Meschke, H.A. Mang Institute for Strength of Materials, Technical University of Vienna, A-1040 Karlsplatz 13/202, Vienna, Austria ABSTRACT
More informationMODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD
MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD Andrzej Skrzat, Rzeszow University of Technology, Powst. Warszawy 8, Rzeszow, Poland Abstract: User-defined material models which can be used
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 information1. Background. is usually significantly lower than it is in uniaxial tension
NOTES ON QUANTIFYING MODES OF A SECOND- ORDER TENSOR. The mechanical behavior of rocks and rock-like materials (concrete, ceramics, etc.) strongly depends on the loading mode, defined by the values and
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 Probabilistic Yielding of Materials
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2; :1 6 [Version: 22/9/18 v1.2] On Probabilistic Yielding of Materials Boris Jeremić and Kallol Sett Department of Civil and
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 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 informationCourse in. Geometric nonlinearity. Nonlinear FEM. Computational Mechanics, AAU, Esbjerg
Course in Nonlinear FEM Geometric nonlinearity Nonlinear FEM Outline Lecture 1 Introduction Lecture 2 Geometric nonlinearity Lecture 3 Material nonlinearity Lecture 4 Material nonlinearity it continued
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 informationSTUDY OF THE BARCELONA BASIC MODEL. INFLUENCE OF SUCTION ON SHEAR STRENGTH
STUDY OF TH BARCLONA BASIC MODL. INFLUNC OF SUCTION ON SHAR STRNGTH Carlos Pereira ABSTRACT The Barcelona Basic Model, BBM, is one of the most used elasto-plastic models for unsaturated soils. This summary
More informationDiscrete Element Modeling of Soils as Granular Materials
Discrete Element Modeling of Soils as Granular Materials Matthew R. Kuhn Donald P. Shiley School of Engineering University of Portland National Science Foundation Grant No. NEESR-936408 Outline Discrete
More informationA 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 informationReturn Mapping Algorithms and Stress Predictors for Failure Analysis in Geomechanics
Return Mapping Algorithms and Stress Predictors for Failure Analysis in Geomechanics Jinsong Huang 1 and D. V. Griffiths 2 Abstract: Two well-known return mapping algorithms, the closest point projection
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 informationCOUPLED FINITE-INFINITE ELEMENTS MODELING OF BUILDING FRAME-SOIL INTERACTION SYSTEM
VOL. 4, NO. 10, DECEMBER 009 SSN 1819-6608 006-009 Asian Research Publishing Network (ARPN. All rights reserved. COUPLED FNTE-NFNTE ELEMENTS MODELNG OF BULDNG FRAME-SOL NTERACTON SYSTEM Ramakant Agrawal
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 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 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 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 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 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 informationBAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS
Journal of Computational and Applied Mechanics, Vol.., No. 1., (2005), pp. 83 94 BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Vladimír Kutiš and Justín Murín Department
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 informationNon-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises
Non-linear and time-dependent material models in Mentat & MARC Tutorial with Background and Exercises Eindhoven University of Technology Department of Mechanical Engineering Piet Schreurs July 7, 2009
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 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 informationA Study on Modelling the Plane Strain Behaviour of Sand and its Stability
Archives of Hydro-Engineering and Environmental Mechanics Vol. 59 (2012), No. 3 4, pp. 85 100 DOI: 10.2478/heem-2013-0001 IBW PAN, ISSN 1231 3726 A Study on Modelling the Plane Strain Behaviour of Sand
More informationPropagation of Seismic Waves through Liquefied Soils
Propagation of Seismic Waves through Liquefied Soils Mahdi Taiebat a,b,, Boris Jeremić b, Yannis F. Dafalias b,c, Amir M. Kaynia a, Zhao Cheng d a Norwegian Geotechnical Institute, P.O. Box 393 Ullevaal
More informationOn the Numerical Modelling of Orthotropic Large Strain Elastoplasticity
63 Advances in 63 On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity I. Karsaj, C. Sansour and J. Soric Summary A constitutive model for orthotropic yield function at large strain
More informationFROM THEORY TO PRACTICE IN GEOTECHNICAL ENGINEERING: THE MISSING LINKS. Yannis F. Dafalias
FROM THEORY TO PRACTICE IN GEOTECHNICAL ENGINEERING: THE MISSING LINKS Yannis F. Dafalias Department of Civil and Environmental Engineering, University of California at Davis, & Department of Mechanics,
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 information3D MATERIAL MODEL FOR EPS RESPONSE SIMULATION
3D MATERIAL MODEL FOR EPS RESPONSE SIMULATION A.E. Swart 1, W.T. van Bijsterveld 2, M. Duškov 3 and A. Scarpas 4 ABSTRACT In a country like the Netherlands, construction on weak and quite often wet soils
More informationLETTERS TO THE EDITOR
NTERNATONAL JOURNAL FOR NUMERCAL AND ANALYTCAL METHODS N GEOMECHANCS, VOL. 13, 101-107 (1989) LETTERS TO THE EDTOR A LARGE STRAN THEORY AND TS APPLCATON N THE ANALYSS OF THE CONE PENETRATON MECHANSM (P.
More informationA MODEL FOR COUPLED DYNAMIC ELASTO-PLASTIC ANALYSIS OF SOILS
Journal of GeoEngineering, Vol. 7, No. 3, pp. 089-096, December 2012 Fattah et al.: A Model for Coupled Dynamic Elastic Plastic Analysis of Soils 89 A MODEL FOR COUPLED DYNAMIC ELASTO-PLASTIC ANALYSIS
More informationImprovement of a hypoplastic model to predict clay behaviour under undrained conditions
Improvement of a hypoplastic model to predict clay behaviour under undrained conditions Mašín, D. and Herle, I. September 3, 27 Accepted for publication in Acta Geotechnica ABSTRACT A shortcoming of the
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 informationA Krylov Subspace Accelerated Newton Algorithm
A Krylov Subspace Accelerated Newton Algorithm Michael H. Scott 1 and Gregory L. Fenves 2 ABSTRACT In the nonlinear analysis of structural systems, the Newton-Raphson (NR) method is the most common method
More informationSTRAIN LOCALIZATION AS BIFURCATION ELASTO-PLASTIC SOFTENING MATERIALS
Fracture Mechanics of Concrete Structures, Proceedings FRAMCOS-2, edited by Folker H. Wittmann, AEDIFICA TIO Publishers, D-79104 Frei burg (1995) STRAIN LOCALIZATION AS BIFURCATION ELASTO-PLASTIC SOFTENING
More informationTable of Contents. Foreword... xiii Introduction... xv
Foreword.... xiii Introduction.... xv Chapter 1. Controllability of Geotechnical Tests and their Relationship to the Instability of Soils... 1 Roberto NOVA 1.1. Introduction... 1 1.2. Load control... 2
More informationINVERSE ANALYSIS METHODS OF IDENTIFYING CRUSTAL CHARACTERISTICS USING GPS ARRYA DATA
Problems in Solid Mechanics A Symposium in Honor of H.D. Bui Symi, Greece, July 3-8, 6 INVERSE ANALYSIS METHODS OF IDENTIFYING CRUSTAL CHARACTERISTICS USING GPS ARRYA DATA M. HORI (Earthquake Research
More informationNumerical Analysis of Pile Behavior under Lateral Loads in. Layered Elastic Plastic Soils
INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS Int. J. Numer. Anal. Meth. Geomech. ; : 3 [Version: /3/ v.] Numerical Analysis of Pile Behavior under Lateral Loads in Layered
More informationObservations on Return Mapping Algorithms for Piecewise Linear Yield Criteria
Observations on Return Mapping Algorithms for Piecewise Linear Yield Criteria Jinsong Huang 1 and D. V. Griffiths 2 Abstract: This paper shows that for perfect plasticity, the closest point projection
More informationTransactions on Engineering Sciences vol 6, 1994 WIT Press, ISSN
Spatial discretization of strain localization E. Rizzi,* K. Willam University of Colorado Boulder, CEAE Department, Boulder, Colorado 80309-0428, USA * On leave from Politecnico di Milano, DIS, Milano,
More information16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations
6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =
More informationMGM-III. STS-107 IWG Presentation. Stein Sture, University of Colorado at Boulder Nicholas Costes, University of Colorado at Boulder
MGM-III STS-107 IWG Presentation Principal Investigator Co- Investigator Project Scientist Presented by Stein Sture, University of Colorado at Boulder Nicholas Costes, University of Colorado at Boulder
More informationHÅLLFASTHETSLÄRA, LTH Examination in computational materials modeling
HÅLLFASTHETSLÄRA, LTH Examination in computational materials modeling TID: 2016-28-10, kl 14.00-19.00 Maximalt 60 poäng kan erhållas på denna tenta. För godkänt krävs 30 poäng. Tillåtet hjälpmedel: räknare
More informationUCD CompGeoMech Contributions to OpenSees: Deliverables
UCD CompGeoMech Contributions to OpenSees: Deliverables Prof. Yang, Dr. Cheng, GSRs Jie, Sett, Taiebat, Preisig, Wu, Liu Boris Jeremić, University of California, Davis OpenSees Annual Meeting, San Francisco,
More informationLocalization in Undrained Deformation
Localization in Undrained Deformation J. W. Rudnicki Dept. of Civil and Env. Engn. and Dept. of Mech. Engn. Northwestern University Evanston, IL 621-319 John.Rudnicki@gmail.com Fourth Biot Conference on
More informationTangent Modulus in Numerical Integration of Constitutive Relations and its Influence on Convergence of N-R Method
Applied and Computational Mechanics 3 (2009) 27 38 Tangent Modulus in Numerical Integration of Constitutive Relations and its Influence on Convergence of N-R Method R. Halama a,, Z. Poruba a a Faculty
More informationMicro-mechanics in Geotechnical Engineering
Micro-mechanics in Geotechnical Engineering Chung R. Song Department of Civil Engineering The University of Mississippi University, MS 38677 Fundamental Concepts Macro-behavior of a material is the average
More 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 informationNumerical Simulation of Fully Saturated Porous Materials
INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS Numerical Simulation of Fully Saturated Porous Materials Boris Jeremić ( ) 1, Zhao Cheng ( ) 2 3 (مهدی طيبات ( Taiebat Mahdi and
More informationMODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Class Meeting #5: Integration of Constitutive Equations
MODELING OF CONCRETE MATERIALS AND STRUCTURES Kaspar Willam University of Colorado at Boulder Class Meeting #5: Integration of Constitutive Equations Structural Equilibrium: Incremental Tangent Stiffness
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 informationCONSOLIDATION BEHAVIOR OF PILES UNDER PURE LATERAL LOADINGS
VOL., NO., DECEMBER 8 ISSN 89-8 -8 Asian Research Publishing Network (ARPN). All rights reserved. CONSOLIDATION BEAVIOR OF PILES UNDER PURE LATERAL LOADINGS Qassun S. Mohammed Shafiqu Department of Civil
More informationContinuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms
Continuum mechanics office Math 0.107 ales.janka@unifr.ch http://perso.unifr.ch/ales.janka/mechanics Mars 16, 2011, Université de Fribourg 1. Constitutive equation: definition and basic axioms Constitutive
More informationDevelopment of Geotechnical Capabilities in OpenSees
Development of Geotechnical Capabilities in OpenSees PEER 2132-3 Boris Jeremić University of California Pacific Earthquake Engineering Research Center University of California, Berkeley September 21 ii
More informationThe Role of Slope Geometry on Flowslide Occurrence
American Journal of Environmental Sciences 3 (3): 93-97, 27 ISSN 1553-345X 27 Science Publications Corresponding Author: The Role of Slope Geometry on Flowslide Occurrence Chiara Deangeli DITAG, Politecnico
More informationfor material models with uncertain parameters. Developments in this paper are based on already derived general formulation presented in the companion
Acta Geotechnica manuscript No. (will be inserted by the editor) Probabilistic Elasto-Plasticity: Solution and Verification in 1D Kallol Sett 1, Boris Jeremić 2, M. Levent Kavvas 3 1 Graduate Research
More informationFE FORMULATIONS FOR PLASTICITY
G These slides are designed based on the book: Finite Elements in Plasticity Theory and Practice, D.R.J. Owen and E. Hinton, 1970, Pineridge Press Ltd., Swansea, UK. 1 Course Content: A INTRODUCTION AND
More informationNumerical analysis of pile behaviour under lateral loads in layered elastic plastic soils
INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS Int. J. Numer. Anal. Meth. Geomech., ; 6:1385 146 (DOI: 1.1/nag.5) Numerical analysis of pile behaviour under lateral loads in
More informationNumerical Modelling of Blockwork Prisms Tested in Compression Using Finite Element Method with Interface Behaviour
13 th International Brick and Block Masonry Conference Amsterdam, July 4-7, 2004 Numerical Modelling of Blockwork Prisms Tested in Compression Using Finite Element Method with Interface Behaviour H. R.
More informationDrained Against Undrained Behaviour of Sand
Archives of Hydro-Engineering and Environmental Mechanics Vol. 54 (2007), No. 3, pp. 207 222 IBW PAN, ISSN 1231 3726 Drained Against Undrained Behaviour of Sand Andrzej Sawicki, Waldemar Świdziński Institute
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 informationUnderstand basic stress-strain response of engineering materials.
Module 3 Constitutive quations Learning Objectives Understand basic stress-strain response of engineering materials. Quantify the linear elastic stress-strain response in terms of tensorial quantities
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 informationBifurcation Analysis in Geomechanics
Bifurcation Analysis in Geomechanics I. VARDOULAKIS Department of Engineering Science National Technical University of Athens Greece and J. SULEM Centre d'enseignement et de Recherche en Mecanique des
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 informationINTERPRETATION OF UNDRAINED SHEAR STRENGTH OF UNSATURATED SOILS IN TERMS OF STRESS STATE VARIABLES
INTERPRETATION OF UNDRAINED SHEAR STRENGTH OF UNSATURATED SOILS IN TERMS OF STRESS STATE VARIABLES S. K. Vanapalli and D.G. Fredlund Department of Civil Engineering University of Saskatchewan, Saskatoon
More informationAnalysis of Blocky Rock Slopes with Finite Element Shear Strength Reduction Analysis
Analysis of Blocky Rock Slopes with Finite Element Shear Strength Reduction Analysis R.E. Hammah, T. Yacoub, B. Corkum & F. Wibowo Rocscience Inc., Toronto, Canada J.H. Curran Department of Civil Engineering
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 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 informationGEOMETRIC NONLINEAR ANALYSIS
GEOMETRIC NONLINEAR ANALYSIS The approach for solving problems with geometric nonlinearity is presented. The ESAComp solution relies on Elmer open-source computational tool [1] for multiphysics problems.
More informationCh 4a Stress, Strain and Shearing
Ch. 4a - Stress, Strain, Shearing Page 1 Ch 4a Stress, Strain and Shearing Reading Assignment Ch. 4a Lecture Notes Sections 4.1-4.3 (Salgado) Other Materials Handout 4 Homework Assignment 3 Problems 4-13,
More informationENGN 2290: Plasticity Computational plasticity in Abaqus
ENGN 229: Plasticity Computational plasticity in Abaqus The purpose of these exercises is to build a familiarity with using user-material subroutines (UMATs) in Abaqus/Standard. Abaqus/Standard is a finite-element
More information1.8 Unconfined Compression Test
1-49 1.8 Unconfined Compression Test - It gives a quick and simple measurement of the undrained strength of cohesive, undisturbed soil specimens. 1) Testing method i) Trimming a sample. Length-diameter
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationSoil Uncertainty and Seismic Ground Motion
Boris Jeremić and Kallol Sett Department of Civil and Environmental Engineering University of California, Davis GEESD IV May 2008 Support by NSF CMMI 0600766 Outline Motivation Soils and Uncertainty Stochastic
More informationContinuous Surface Cap Model for Geomaterial Modeling: A New LS-DYNA Material Type
7 th International LS-DYNA Users Conference Material Technology () Continuous Surface Cap Model for Geomaterial Modeling: A New LS-DYNA Material Type Leonard E Schwer Schwer Engineering & Consulting Services
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 information(Refer Slide Time: 02:18)
Geology and Soil Mechanics Prof. P. Ghosh Department of Civil Engineering Indian Institute of Technology Kanpur Lecture 40 Shear Strength of Soil - C Keywords: Shear strength of soil, direct shear test,
More informationGeneral method for simulating laboratory tests with constitutive models for geomechanics
General method for simulating laboratory tests with constitutive models for geomechanics Tomáš Janda 1 and David Mašín 2 1 Czech Technical University in Prague, Faculty of Civil Engineering, Czech Republic
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 informationClassical Elasticity and Plasticity
Chapter 2 Classical Elasticity and Plasticity 2.1 Elasticity Fung (1965 provides elegant definitions for the different forms of elasticity theory, and we follow his terminology here. A material is said
More informationUseful Formulae ( )
Appendix A Useful Formulae (985-989-993-) 34 Jeremić et al. A.. CHAPTER SUMMARY AND HIGHLIGHTS page: 35 of 536 A. Chapter Summary and Highlights A. Stress and Strain This section reviews small deformation
More informationComputational models of diamond anvil cell compression
UDC 519.6 Computational models of diamond anvil cell compression A. I. Kondrat yev Independent Researcher, 5944 St. Alban Road, Pensacola, Florida 32503, USA Abstract. Diamond anvil cells (DAC) are extensively
More information