A C 0 FINITE ELEMENT IN GRADIENT PLASTICITY FOR LOCALIZED FAILURE MODES ANALYSIS
|
|
- Doris Lawson
- 6 years ago
- Views:
Transcription
1 COMPUAIONAL MECHANICS New rends and Applications S. Idelsohn, E. Oñate and E. Dvorin (Eds.) CIMNE, Barcelona, Spain 1998 A C 0 FINIE ELEMEN IN GRADIEN PLASICIY FOR LOCALIZED FAILURE MODES ANALYSIS Féri Meftah, Gilles Pijaudier-Cabot, and Jean M. Reynouard L.M.-Cachan, E.N.S de Cachan / C.N.R.S / Université P. et M. Curie 61, Av. du Président Wilson, Cachan Cedex, France meftah@lmt.ens-cachan.fr L.M.-Cachan, E.N.S de Cachan / C.N.R.S / Université P. et M. Curie 61, Av. du Président Wilson, Cachan Cedex, France also at Institut Universitaire de France Gilles-Pijaudier.Cabot@lmt.ens-cachan.fr U.R.G.C.-Structures / Institut National des Sciences Appliquées de Lyon 20, Av. A. Einstein, Villeurbanne Cedex, France reynouar@gcu-beton.insa-lyon.fr Key words: Strain softening, Localization, Regularization, Gradient plasticity, Finite element, Failures analysis. Abstract. A new variational approach in gradient plasticity is derived for the finite element analysis of quasi-brittle failure modes. A consistent integration algorithm for updating the stress and the internal state variables at integration points is formulated in an original nonlocal sense. he model includes implicitly second gradient of the effective plastic strain in the yield condition, which reduces the continuity requirements on field interpolations to the use of C 0 shape functions. A one dimensional example is presented to illustrate the method. 1
2 1 INRODUCION Localization of deformation into narrow bands of intense straining caused by strain softening is a characteristic feature of plastic deformation. It has been experimentally observed in many engineering materials, such as concrete, rocs and soils. Localization phenomena are often associated with a significant reduction of the load-carrying capacity of the structures, hence, the onset of localization is naturally considered as the inception of failure of engineering structures. Considerable efforts have been devoted over the last decade to obtain a comprehensive understanding of the problem and to describe this behavior quantitatively. However, numerous attempts to simulate the behavior with softening plasticity or damage theories failed in the sense that the solution appeared to be fully determined by fineness and the direction of the finite element discretization. he underlying reason of this pathologically mesh dependency is the a local change of character of the governing equations, which results in a loss of well-posedness of the boundary value problem. More details and references on these aspects can be found in the papers by Pijaudier-Cabot et al. 14 and Lasry and Belytscho 7. o remedy the situation, generalized continuum theories should be adopted. hese models incorporate an internal length parameter which plays the role of localization limiter, that is, a parameter that allows to control the localization zone size by preventing loss of ellipticity of governing equations. his internal length parameter can be introduced through the incorporation of higher order derivatives leading to the so-called gradient continuum theories. he capabilities of gradient dependence have been investigated for both softening plasticity and damage models 4,10,12. Recently, formulations and algorithms for gradient dependent models have been presented in a finite element context 3,5,6,9,11,12,15. In a gradient plasticity model the yield strength depends not only on the effective plastic strain but also on its Laplacian. herefore, even if gradient dependent models bear the significant advantage of being local in a mathematical sense, the increment of the plastic strain can not be obtained at a local level since the consistency condition which governs the plastic flow becomes a second order partial differential equation. One possibility is to use a finite difference method. he algorithm is then a sequence of separate approximate solutions of the equilibrium problem using finite elements and the plastic yielding problem using finite differences 1. A more general approach is to use only finite elements and to solve the two coupled problems simultaneously 4. For this purpose, it was required to satisfy wealy the yield condition and to discretize the plastic strain field in addition to the standard discretization of the displacement field. herefore, C 1 continuous shape functions, for the plastic strain field interpolation, are needed in order to compute properly second gradients of the equivalent plastic strain. As a consequence, this mixed formulation leads to an oversized discretized problems maing the calculations unreasonable. With the use of the gradient plastic model including implicitly second gradient of the effective plastic strain in the yield condition, a finite element method for quasi-brittle fracture will be derived in our wor. he consistent integration algorithm to update the stress and the internal state variable at integration points for the gradient plasticity are formulated in an 2
3 original nonlocal sense. he new algorithm could require either a wea or point-wise satisfaction of the yield function together with the consistency condition, while only a C 0 continuity interpolation functions of the plastic problem related field is needed. he methodology to derive the finite element formulations of this model is presented. Finally, the pertinence of the finite element and the elaborated algorithm will be illustrated on a one dimensional numerical example. 2 CONSIUIVE RELAIONS he main feature of gradient plasticity is that the yield function F depends, beside the stress tensor σ, on both the equivalent plastic strain κ and its Laplacean 2 κ j+1. It can tae the following form 4,8, F( σ,κ κ) = f( σ) τ( κ κ), -, =0 (1) where f gives the stress related part and τ is the gradient dependent yield strength which reads (, ) ( )-g( ) 2 2 τκ κ = τκ κ κ. (2) In relation (2) τ is the classical, local, yield strength and g is a weight function which gives the effective nonlocal contribution. t n λ Σ = σ ij Elastic domain (e) Plastic domain (λ) λ Σ: stress tensor Figure 1: Evolution of the plastic domain he variational principle proposed previously by de Borst and Muhlhaus 4 for gradient plasticity allows to this theory to be used as a localization limiter in a large domain of finite 3
4 element calculations 8,11,15. his variational approach consists in a wea satisfaction of the yield function 2 ( σ +1, +1 +1) δλ F λ, λ d = 0 +1 κ besides the classical wea satisfaction of the equilibrium equation. In relation (3), δλ is a virtual variation of the plastic multiplier which is assumed to be related to κ by a linear relation. he above relation means that the yield function should only be satisfied at the end of the current loading step ( +1) and in a distributed manner. It should also be noticed that integration in Eq.3 is restricted to the plastic domain +1 κ since F < 0 in the elastic part of the body. Of course, due to Khun-ucer conditions λ 0, F 0, λ F = 0, (4) one can expect to extend the integration in Eq.3 to the whole domain, however the test function δλ should remain a general function satisfying only smoothness requirement and the boundary conditions of the approximated field. hereafter, since the plastic domain is a priori unnown, a predictor-corrector schema is proposed 4,8,11 for the solution of the algebraic equations resulting from the finite element discretization of the functional. Relation (3) becomes j + j+ 1 j+ 1 2 ( σ +1, +1 +1) δλ F j 1 λ, λ d = 0 +1 κ in which (j) is equilibrium iterations indicator in a Newton -Raphson algorithm. his means that we determine plastic flow at iteration (j+1) by integrating over the plastic domain given by all the integration points which yield (for which the stress state has violated the yield function) at the previous iteration. By this way of proceeding, one neglects the term giving the evolution process of plastic domain leading to a description of the motion of the elasto-plastic boundary λ which is not entirely correct. he consequence of the delay in the evolution of this boundary, is an overconstrain on the plastic flow which leads to important values of the Laplacean of the equivalent plastic strain and then causes numerical troubles during equilibrium iterations. Moreover, this variational principle gives the conditions on λ at the boundary of the plastically deforming part of the body λ, which must either be δλ = 0 or ( λ ) (3) (5) n = 0 (6) n λ is the unit outward normal at the elasto-plastic boundary (Fig.1). herefore, these conditions should be imposed at a moving boundary, a fortiori unnown, which renders the procedure impossible in a classical finite element code. Currently, the conditions are imposed λ 4
5 on only the external part of the elasto-plastic boundary. Finally, in order to avoid oscillations of the yield function values, the Laplacean of the plastic multiplier 2 λ should be properly computed by adopting C 1 continuous interpolation polynomials. his continuity requirement increases significantly the number of degrees of freedom and, thus, the calculation cost. In this contribution, we propose a new variational approach that allows to prevent the main drawbacs of the previous algorithm. For this purpose we start from an earlier definition 13 of a nonlocal quantity. We consider the nonlocal yield strength value τ in a material point x as weighted average of the local yield strength τ over the surrounding volume s in which φ( ξ) 1 τ = φ τ = ξ = exp 2 l 2 2 s ( ξ) ( x + ξ) d 0, with ( ) = φ ξ (7) s d is a gaussian weight function and ξ denotes the relative position vector pointing to the infinitesimal volume d. he parameter l is the internal length governing the contribution on the yield strength (additional carrying capacity) of the surrounding volume at the given material point x. A gradient formulation can be derived from this non local theory 10. he local yield strength is then expanded into aylor series according to τ( x + ) τ+ τ τ ( 2) 1 3 τ () ξ ξ ξ + ξ + τ ξ + 2! 3! 4! 4 (8) where, n is the n th order gradient operator, the dot denotes inner product between n th order tensors and ξ ( n) denotes the n factor dyadic product ξ n - times ξ. By substituting Eq.8 into Eq.7 and considering that the domain s remains small comparing to, an exact integration of relation (7) is feasible leading to the vanishing of the odd terms such that 4 () l c () l. (9) 2 τ τ+c 1 τ+ 2 τ+ he gradient parameters c i are of the dimension of a length to an even power, so that the internal length scale is present in this gradient formulation. hus, neglecting higher order terms in expression (9), we simply assume that the yield condition reads τ l τ τ = 0. (10) In the next section, we will show that this definition is suitable for numerical analysis since it enables a straightforward C 0 -continuous finite element interpolation. Furthermore, a very important feature of this relation is that it holds on the whole domain. Indeed, in this differential equation τ acts as a source term that drives the variation of the local quantity τ. hereafter, a change in the values of τ, notably at the elasto-plastic boundary, will provide 5
6 non zero values of τ in former elastic integration points close to the plastic domain. hus, the localization zone can extend and the regularization effect is ensured. In this contribution, a wea form of this equation will then be considered instead of Eq.1 (which is valid in the plastic domain only) to derive the set of equations governing both of equilibrium and plastic flow. It is important to note that the model does not differ fundamentally from the usual gradient plasticity one derived by de Borst and Muhlhaus 4. Nonlocality is expressed in a slightly different way only. 3 FINIE ELEMEN IMPLEMENAION he wea form of the partial differential equations which govern deformation processes are derived according to the weighted residuals approach 16. he equilibrium equation at the current load step ( +1) L σ +1 = 0 (11) in which L is a differential matrix operator 8 and where body forces have not been included for simplicity, is multiplied by the vectorial weight function δu (virtual displacement field) and subsequently integrated on the domain ( L σ ) δu +1 d = 0. (12) he wea form of the averaging equation (10) is obtained in a similar way. Multiplication with the weight function δτ and integration yields and 2 2 ( + 1 l ) δτ τ + τ τ d = 0. (13) With the aid of the divergence theorem, the resulting equation can be transformed to [( Lu δ ) σ] = [ δu t 1] ( Lu δ ) + [ σ ] d d d Σ d. (14) ( l ( ) ) = l ( ) ( ) 2 2 δτ d τ δτ d τ d δτ d τ d δτ τ d δτ τ δτ τ d(15) in which d denotes increment of the variable ( + 1 = λ + σd, τ + 1 = τ + dτ, τ + 1 = τ + dτ ), and t is the external stress vector acting on the configuration boundary, on which the non standard conditions δτ or ( ) 0 = 0 σdτ n = (16) 6
7 have been considered in equation (15). he last condition can be interpreted as zero history parameter flux through the boundary of the body. herefore, if these extra-boundary conditions are considered, Eq.15 shows no need for C 1 interpolation for the local yield strength. It remains however to treat the source term dτ such that its explicit dependency on the second order gradient can be avoided. During plastic flow, the stress tensor increment dσ can be related to both the strain and the local yield strength increments dε and dτ, respectively. For the sae of simplicity, an associated flow law is considered e [ ε n] d σ = D d d λ (17) in which D e is the elastic stiffness matrix and n = F σ is the gradient vector to the yield function. Further, a linear relation between and d λ = ηdκ (18) and the maximum plastic wor assumption in a local manner dτκ dκ = (19) h are adopted in which h is the softening modulus. his last equation is valid during the softening flow only. We do not consider here the case where the material strength has been completely exhausted. By considering Eqs together with Eq.14, one obtains the following equation e e η [( Lu δ ) D ε] ( Lu δ ) D n τ [ δu t 1] ( Lu δ ) = + ( ) [ σ ] d d d d d Σ d(20) h in which the strain and the local yield strength appear as the unnown fields. In the same way, the source term dτ in equation (15) should be expressed as a function of these two quantities. For this purpose, we mae call to the non local yield function that reads F ( σ ) f( σ),τ = τ = 0 (21) and in which the gradient dependency is made implicit through τ. herefore, the consistency condition gives F F d F( σ, τ) = d σ + d dσ d σ τ τ = n τ (22) and substitution of Eqs into Eq.22 provides the relation giving the increment of the nonlocal yield strength e η e dτ = df + n D d ε n D n dτ (23) h 7
8 in which df represents the difference of values of the yield function at the current (j+1) and previous (j) iteration within the load step (+1) ( σ j, j ) ( σ j, j) d F = F τ F τ. (24) he implementation of expression (23) in the second variational equation (15) governing the problem leads to e η e 2 δτ n D d ε d + δτ 1+ n D n d τ l ( δτ) d τ d h. (25) = δτ τ d δτ τ δτ τ d δτ d Fd 2 ( j l ( ) j) j where it can be clearly noticed that all the implied fields only require a C 0 continuity. Departing from the wea forms given by Eqs. (20) and (25), the finite element discretization is rather straightforward. Following the Galerin approach, the displacement and strain fields are discretized according to the normal finite element procedure u = Na, ε = Ba. (26) where N is the interpolation matrix containing the shape functions, B is the matrix of the shape function derivatives ( B = L N) and a is the nodal displacement vector. Further, a separate interpolation of the local yield strength is introduced τ=h Τ, τ=q Τ (27) where H denotes the interpolation functions, is the nodal values vector of the local yield strength and Q = [ H i ] containing the gradient of interpolation functions. It is emphasized that the interpolation polynomials of u and t do not need to be of the same order. Both discretization only need to satisfy C 0 -continuity requirements. o avoid stress oscillations, the use of an interpolation for the displacements which is one order higher than that of the local strength seems advisable (the Babusa-Brezzi conditions for mixed finite elements in compressible solids). Substituting the above identities in Eqs. (20) and (25) and requiring that these equations hold for any admissible variations δ a and δ Τ lead to the following set of algebraic equations which describe the incremental process in the discretized gradient enhanced elastoplastic continuum K K aa τa K K aτ da fe +1 fi d = (28) Τ f f ~ ττ τ τ where the elastic stiffness matrix K aa, the external force vector f e and the internal force vector f i are defined conventionally as 8
9 K = B D Bd, f N t aa the off-diagonal matrices K aτ and K τa are e e +1 = + 1 d i Σ, f = B σ d (29) K aτ and the gradient-dependent matrix K ττ is η e e = B D n Hd, Kτa H n D B d h =, (30) η e 2 Kττ = 1 + n D n H H l Q Q h d. (31) he vectors of non-standard residual forces which emerge from the inexact fulfillment of the equation (10), defining the non local yield function, read f τ ( τ ) = H d F d and f~ τ = [ H τ Q τ] l 2 d Note that the tangent stiffness operator in the system (28) is non-symmetric due to the expression of the off-diagonal matrices K aτ and K τa. 4 ALGORIHM his algorithm is base on a return mapping schema using the well nown elastic predictorplastic corrector procedure. We give here the different steps followed during the iterative process. 1. Solve for the nodal increments d ε = B da and d τ = H dτ according to Eq. (28) and update using a total incremental approach ε = ε + d ε, τ = τ +dτ (33) j+1 2. Compute the strain and the local yield strength at integration points j j+1 j (32) ε = ε + ε, τ = τ + τ (34) j+1 0 j+1 j+1 0 j+1 3. Compute the elastic predictor and estimate the loading function tr j+1 e tr tr = 0 +D j+1, Fj+1 F( j+1, j ) σ σ ε = σ τ (35) tr 4. If F j+1 0 then compute 4.1 he new stress state 9
10 e η σ j+1 = D ε j+ 1 τ h ( τ j+ 1) j+1 F σ σ=σ 4.2 the approximate increment of nonlocal yield strength else ( σ j+1, j+1) ( σ j, j) tr ; σ j+1 = σ0 + σ j+1 (36) e η e F = F τ F τ ; τ j+ 1 = Fj+ 1+ n D ε j+ 1 n D n τj+ 1 (37) h In a first approximation the value of the local yield strength is used for evaluating F j+1. τ = τ + τ, n 0 j+1 j j+1 =, η ( τ) h = 0, F j+ 1 = 0 (38) For elastic points the source term τ varies by leading to a variation of the local yield strength in the surrounding elastic domain near the plastic zone and thus provides an extension of the localization zone. Indeed, as a consequence of Eq.38, in the set of equations (28) the submatrices K aτ K τa vanish, but the term fτ f~ τ 0 such that non zero values of the local yield strength are obtained. 5. Update the internal and non standard residual forces. If the process has not converged, compute the new tangential stiffness operator and enter a new iteration. 5 NUMERICAL EXAMPLE y d σ x L Figure 2: Imperfect bar in tension he performance of the model and the algorithm is illustrated with a simple one dimensional problem of bar in tension (Fig.2). In the calculations, the length of the bar is L = 100 mm, Young s modulus E = N mm 2 2, the tensile strength f t = 2 N mm and the softening modulus h = 2000 N mm 2 (a linear softening is considered). he internal length scale is set to l = 5 mm. At the center of the bar, an imperfection zone ( d = 5 mm) with a 10% smaller value of f t is assumed. Calculations have been carried out for meshes consisting of 40, 80 and 160 one dimensional elements with a quadratic interpolation for the displacements and a linear interpolation for the local yield strength. hey are stopped before the local yield strength vanishes. 10
11 2.0 σ [ N mm 2 ] 9.0 ε [ 10 4 ] elements 80 elements 160 elements [ mm ] u τ [ N mm 2 ] x [mm] τ [ N mm 2 ] x [mm] x [mm] Figure 3: Load elongation diagrams for 40, 80 and 160 elements (a). Evolution of the strain (b), local (c) and non local (d) yield strengths along the bar for progressive loading. he load-elongation curves for the three discretizations are depicted in Fig.3 and show a convergence towards a solution with a mesh independent, finite energy dissipation. he same slope of the stress-displacement diagram is obtained upon mesh refinement. Furthermore, this figure shows the evolution of the strain, local and non local yield strengths in the bar for progressive loading. It can be observed that the algorithm allows to the source term τ to vary outside the imperfect zone such that the local yield strength follows. he localization zone can therefore extend to reach a finite size. 11
12 6. CONCLUSION A new variational approach in gradient plasticity has been derived. It consists in a wea satisfaction of the partial differential equation giving the non local history parameter, which is obtained from the averaging integral procedure. An advantage of this approach is that boundary conditions with regard to the second field can be taen into account in a natural way, and at the fixed external boundary of the body. Further, this algorithm requires only a C 0 continuous interpolation of the two fields. he incorporation of the gradient term in the mathematical description introduces an internal length scale in the model. he one dimensional example illustrates that this internal length scale preserves the well-posedness of the problem in the softening regime, same as in usual gradient plasticity models. Finally, some extension of the model should be considered, namely the incorporation of non linear softening flow and the treatment of the elastic zone and when the material strength has been exhausted. REFERENCES [1]. Belytscho and D. Lasry, A study of localization limiters for strain-softening in statics and dynamics. Comp. & Struct., 33, (1989). [2] Z.P. Bazant and G. Pijaudier-Cabot, Non-local continuum damage, localization instability and convergence. ASME J. Appl. Mech., 55, (1988). [3] A. Benallal, R. de Borst and O.M. Heeres, A gradient enhanced damage model: theory and computation. In : D.R.J. Owen, E. Onate and E. Hinlton (Eds.) : Proc. fifth. Int. Conf. on Computational Plasticity, heory and Applications. Barcelona : CIMNE, (1997). [4] R de Borst and H.B. Muhlhaus, Gradient-dependent plasticity: Formulation and algorithmic aspects. Int. J. Num. Meth. Eng., 35, (1992). [5] R. de Borst, L.J. Sluys, H.B. Muhlhaus, and J. Pamin, Fundamental issues in finite element analyses of localization of deformation. Eng. Comput., 10, (1993). [6] C. Comi and U. Perego, A generalized variable formulation for gradient dependent plasticity softening plasticity. Int. J. Num. Meth. Eng., 39, (1996). [7] D. Lasry and. Belytscho, Localization limiters in transient problems. Int. J. Solids Struct., 24, (1988). [8] F. Meftah, Contribution à l etude numerique des modes localisés de rupture dans les structures en bétons de type poutre. approche multicouches par la plasticité au gradient. Ph.D. hesis, INSA of Lyon, France, 220, (1997). [9] F. Meftah and J. M. Reynouard, A multilayered beam element in gradient plasticity for the analysis of localized failure modes. In Int. J. Mech. of Coh. and Frict. Mat., Accepted for publication September [10] H.B. Muhlhaus and E.C. Aifantis, A variational principle for gradient plasticity. Int. J. Solids Struct., 28, (1991). 12
13 [11] J. Pamin, Gradient-dependent plasticity in numerical simularion of localization phenomena. Dissertation, Delft University of echnology, Delft, he Netherlands, 134, (1994). [12] R.H.J Peerlings, R. de Borst, W.A.M. Breelmans and J.H.P. de Vree, Computational modelling of gradient enhanced damage for fracture and fatigue problems. In : D.R.J. Owen, E. Onate and E. Hinlton (Eds.) : Proc. Fourth Int. Conf. on Computational Plasticity, heory and Applications. Swansea : Pineridge Press, (1995). [13] G. Pijaudier-Cabot and Z.P. Bazant, Nonlocal damage theory. ASCE J. Eng. Mech., 113, (1987). [14] G. Pijaudier-Cabot and Z.P. Bazant and M. abbara, Comparison of various models for strain softening. Eng. Comput., 5, (1988). [15] L.J. Sluys, Wave propagation, localization and dispersion in softening solids. Dissertation, Delft University of echnology, Delft, 173 (1992). [16] O.C. Zieniewicz O.C. and R.L. aylor, he Finite Element Method, Fourth edition. London : McGraw-Hill, Vol. 1, (1991). 13
Modelling of ductile failure in metal forming
Modelling of ductile failure in metal forming H.H. Wisselink, J. Huetink Materials Innovation Institute (M2i) / University of Twente, Enschede, The Netherlands Summary: Damage and fracture are important
More informationTransactions on Engineering Sciences vol 6, 1994 WIT Press, ISSN
Significance of the characteristic length for micromechanical modelling of ductile fracture D.-Z. Sun, A. Honig Fraunhofer-Institut fur Werkstoffmechanik, Wohlerstr. 11, D-79108 Freiburg, Germany ABSTRACT
More informationA Method for Gradient Enhancement of Continuum Damage Models
TECHNISCHE MECHANIK, Band 28, Heft 1, (28), 43 52 Manuskripteingang: 31. August 27 A Method for Gradient Enhancement of Continuum Damage Models B. J. Dimitrijevic, K. Hackl A method for the regularization
More informationAbstract. 1 Introduction
Contact analysis for the modelling of anchors in concrete structures H. Walter*, L. Baillet** & M. Brunet* *Laboratoire de Mecanique des Solides **Laboratoire de Mecanique des Contacts-CNRS UMR 5514 Institut
More information1 Introduction Fracture of engineering components is often preceded by considerable changes in the microstructure of the material they are made of. Ac
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2000 Barcelona, 11-14 September 2000 cfleccomas ENHANCED DAMAGE MODELLING OF QUASI-BRITTLE AND FATIGUE FRACTURE COMPUTATIONAL
More informationNUMERICAL MODELLING AND DETERMINATION OF FRACTURE MECHANICS PARAMETERS FOR CONCRETE AND ROCK: PROBABILISTIC ASPECTS
NUMERICAL MODELLING AND DETERMINATION OF FRACTURE MECHANICS PARAMETERS FOR CONCRETE AND ROCK: PROBABILISTIC ASPECTS J. Carmeliet Catholic University of Leuven, Department of Civil Engineering, Belgium
More informationNonlocal computational methods applied to composites structures
Nonlocal computational methods applied to composites structures Norbert Germain, Frédéric Feyel, Jacques Besson To cite this version: Norbert Germain, Frédéric Feyel, Jacques Besson. Nonlocal computational
More information1 Nonlinear deformation
NONLINEAR TRUSS 1 Nonlinear deformation When deformation and/or rotation of the truss are large, various strains and stresses can be defined and related by material laws. The material behavior can be expected
More informationCracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics
Cracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics Tobias Gasch, PhD Student Co-author: Prof. Anders Ansell Comsol Conference 2016 Munich 2016-10-12 Contents Introduction Isotropic damage
More informationAn orthotropic damage model for crash simulation of composites
High Performance Structures and Materials III 511 An orthotropic damage model for crash simulation of composites W. Wang 1, F. H. M. Swartjes 1 & M. D. Gan 1 BU Automotive Centre of Lightweight Structures
More informationCracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics
Cracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics Tobias Gasch *1 and Anders Ansell 1 1 KTH Royal Institute of Technology, Department of Civil and Architectural Engineering *Corresponding
More informationA RATE-DEPENDENT MULTI-SCALE CRACK MODEL FOR CONCRETE
VIII International Conference on Fracture echanics of Concrete and Concrete Structures FraCoS-8 J.G.. Van ier, G. Ruiz, C. Andrade, R.C. Yu and X.X. Zhang (Eds) A RATE-DEPENDENT ULTI-SCALE CRACK ODEL FOR
More informationstrain appears only after the stress has reached a certain critical level, usually specied by a Rankine-type criterion in terms of the maximum princip
Nonlocal damage models: Practical aspects and open issues Milan Jirasek LSC-DGC, Swiss Federal Institute of Technology at Lausanne (EPFL), Switzerland Milan.Jirasek@ep.ch Abstract: The purpose of this
More informationHOW TO SIMPLIFY RETURN-MAPPING ALGORITHMS IN COMPUTATIONAL PLASTICITY: PART 2 - IMPLEMENTATION DETAILS AND EXPERIMENTS
How to simplify return-mapping algorithms in computational plasticity: Part 2 Implementation details and experiments XIII International Conference on Computational Plasticity. Fundamentals and Applications
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 information3D crack propagation with XFEM cohesive elements
Transactions, SMiRT-23, Paper ID 472 3D crack propagation with XFEM cohesive elements Patrick Massin 1, Guilhem Ferté 1,2, and Nicolas Moës 2 1 Institute of Mechanical Sciences and Industrial Applications,
More informationTENSILE CRACKING VIEWED AS BIFURCATION AND INSTABILITY IN A DISCRETE INTERFACE MODEL
Fracture Mechanics of Concrete Structures Proceeding FRAMCOS-3 AEDIFICATIO Publishers, D-79104 Frei burg, Germany TENSILE CRACKING VIEWED AS BIFURCATION AND INSTABILITY IN A DISCRETE INTERFACE MODEL A.
More informationCombined Isotropic-Kinematic Hardening Laws with Anisotropic Back-stress Evolution for Orthotropic Fiber-Reinforced Composites
Combined Isotropic-Kinematic Hardening Laws with Antropic Back-stress Evolution for Orthotropic Fiber- Reinforced Composites Combined Isotropic-Kinematic Hardening Laws with Antropic Back-stress Evolution
More informationA comparison between the Perzyna viscoplastic model and the Consistency viscoplastic model
European Journal of Mechanics A/Solids 21 (2002) 1 12 A comparison between the Perzyna viscoplastic model and the Consistency viscoplastic model Otto M. Heeres, Akke S.J. Suiker, René de Borst Koiter Institute
More informationFinite Element Solution of Nonlinear Transient Rock Damage with Application in Geomechanics of Oil and Gas Reservoirs
Finite Element Solution of Nonlinear Transient Rock Damage with Application in Geomechanics of Oil and Gas Reservoirs S. Enayatpour*1, T. Patzek2 1,2 The University of Texas at Austin *Corresponding author:
More informationTIME-DEPENDENT MESOSCOPIC MODELLING OF MASONRY USING EMBEDDED WEAK DISCONTINUITIES
XI International Conference on Computational Plasticity. Fundamentals and Applications COMPLAS 2011 E. Oñate and D.R.J. Owen (Eds) TIME-DEPENDENT MESOSCOPIC MODELLING OF MASONRY USING EMBEDDED WEAK DISCONTINUITIES
More informationA MULTIFIELD VARIATIONAL FE APPROACH FOR NONLOCAL DAMAGE
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau (eds.) R. Owen and M. Mikkola (assoc. eds.) Jyväskylä,
More informationMilan Jirasek 1 Introduction In the context of standard continuum theories, strain-softening constitutive models typically lead to ill-posed boundary
ECCM '99 European Conference on Computational Mechanics August 31 { September 3 Munchen, Germany Computational Aspects of Nonlocal Models Milan Jirasek Laboratory of Structural and Continuum Mechanics
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 informationCommon pitfalls while using FEM
Common pitfalls while using FEM J. Pamin Instytut Technologii Informatycznych w Inżynierii Lądowej Wydział Inżynierii Lądowej, Politechnika Krakowska e-mail: JPamin@L5.pk.edu.pl With thanks to: R. de Borst
More 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 informationMHA042 - Material mechanics: Duggafrågor
MHA042 - Material mechanics: Duggafrågor 1) For a static uniaxial bar problem at isothermal (Θ const.) conditions, state principle of energy conservation (first law of thermodynamics). On the basis of
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 informationCracked concrete structures under cyclic load
Cracked concrete structures under cyclic load Fabrizio Barpi & Silvio Valente Department of Structural and Geotechnical Engineering, Politecnico di Torino, Torino, Italy ABSTRACT: The safety of cracked
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 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 informationDamage and crack modeling in single-edge and doubleedge notched concrete beams
Engineering Fracture Mechanics 65 (2000) 247±261 www.elsevier.com/locate/engfracmech Damage and crack modeling in single-edge and doubleedge notched concrete beams M.G.D. Geers a,b, *, R. de Borst c, R.H.J.
More informationNONLOCAL PLASTICITY APPLIED TO BENDING OF BEAMS
IX International Conference on Computational Plasticity COMPLAS IX E. Oñate and D. R. J. Owen (Eds) CIMNE, Barcelona, 2007 NONLOCAL PLASTICITY APPLIED TO BENDING OF BEAMS L. Strömberg Div of Structural
More informationA Performance Modeling Strategy based on Multifiber Beams to Estimate Crack Openings ESTIMATE in Concrete Structures CRACK
A Performance Modeling Strategy based on Multifiber Beams to Estimate Crack Openings ESTIMATE in Concrete Structures CRACK A. Medjahed, M. Matallah, S. Ghezali, M. Djafour RiSAM, RisK Assessment and Management,
More informationAnalytical formulation of Modified Upper Bound theorem
CHAPTER 3 Analytical formulation of Modified Upper Bound theorem 3.1 Introduction In the mathematical theory of elasticity, the principles of minimum potential energy and minimum complimentary energy are
More informationFracture Mechanics of Concrete Structures Proceedings FRAMCOS-3 AEDIFICATIO Publishers, D Freiburg, Germany
Fracture Mechanics of Concrete Structures Proceedings FRAMCOS-3 AEDIFICATIO Publishers, D-79104 Freiburg, Germany MODELLING OF CRACK PROPAGATION WITH EMBEDDED DISCONTINUITY ELEMENTS L.J. Sluys, Koiter
More informationCorrection of local-linear elasticity for nonlocal residuals: Application to Euler-Bernoulli beams
Correction of local-linear elasticity for nonlocal residuals: Application to Euler-Bernoulli beams Mohamed Shaat* Engineering and Manufacturing Technologies Department, DACC, New Mexico State University,
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 informationPLAXIS. Scientific Manual
PLAXIS Scientific Manual 2016 Build 8122 TABLE OF CONTENTS TABLE OF CONTENTS 1 Introduction 5 2 Deformation theory 7 2.1 Basic equations of continuum deformation 7 2.2 Finite element discretisation 8 2.3
More informationThe concept of Representative Volume for elastic, hardening and softening materials
The concept of Representative Volume for elastic, hardening and softening materials Inna M. Gitman Harm Askes Lambertus J. Sluys Oriol Lloberas Valls i.gitman@citg.tudelft.nl Abstract The concept of the
More information3D Finite Element analysis of stud anchors with large head and embedment depth
3D Finite Element analysis of stud anchors with large head and embedment depth G. Periškić, J. Ožbolt & R. Eligehausen Institute for Construction Materials, University of Stuttgart, Stuttgart, Germany
More informationALGORITHM FOR NON-PROPORTIONAL LOADING IN SEQUENTIALLY LINEAR ANALYSIS
9th International Conference on Fracture Mechanics of Concrete and Concrete Structures FraMCoS-9 Chenjie Yu, P.C.J. Hoogenboom and J.G. Rots DOI 10.21012/FC9.288 ALGORITHM FOR NON-PROPORTIONAL LOADING
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 informationStructural behaviour of traditional mortise-and-tenon timber joints
Structural behaviour of traditional mortise-and-tenon timber joints Artur O. Feio 1, Paulo B. Lourenço 2 and José S. Machado 3 1 CCR Construtora S.A., Portugal University Lusíada, Portugal 2 University
More informationVirtual distortions applied to structural modelling and sensitivity analysis. Damage identification testing example
AMAS Workshop on Smart Materials and Structures SMART 03 (pp.313 324) Jadwisin, September 2-5, 2003 Virtual distortions applied to structural modelling and sensitivity analysis. Damage identification testing
More informationA non-local model with tension and compression damage mechanisms
Eur. J. Mech. A/Solids 20 2001) 1 22 2001 Éditions scientifiques et médicales Elsevier SAS. All rights reserved S0997-753800)01111-6/FLA A non-local model with tension and compression damage mechanisms
More informationModelling Localisation and Spatial Scaling of Constitutive Behaviour: a Kinematically Enriched Continuum Approach
Modelling Localisation and Spatial Scaling of Constitutive Behaviour: a Kinematically Enriched Continuum Approach Giang Dinh Nguyen, Chi Thanh Nguyen, Vinh Phu Nguyen School of Civil, Environmental and
More informationContemporary Research in Engineering Science
Bazant, Z.P., and Prat, P.C. (1995). 'Elastic material with systems of growing or closing cracks: Tangential Stiffness.' in Contemporary Research in Engineering Science (Proc., Eringen Medal Symp. honoring
More informationEnhanced coupled elasto-plastic-damage models to describe concrete behaviour in cyclic laboratory tests: comparison and improvement
Arch. Mech., 64, 3, pp. 227 259, Warszawa 2012 Enhanced coupled elasto-plastic-damage models to describe concrete behaviour in cyclic laboratory tests: comparison and improvement I. MARZEC, J. TEJCHMAN
More informationModeling issues of the FRP detachment phenomenon
Modeling issues of the FRP detachment phenomenon Elio Sacco in collaboration with: J. Toti, S. Marfia and E. Grande 1 Dipartimento di ngegneria Civile e Meccanica Università di Cassino e del Lazio Meridionale
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 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 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 informationNUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS
NUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS ABSTRACT : P Mata1, AH Barbat1, S Oller1, R Boroschek2 1 Technical University of Catalonia, Civil Engineering
More informationAVOIDING FRACTURE INSTABILITY IN WEDGE SPLITTING TESTS BY MEANS OF NUMERICAL SIMULATIONS
Damage, Avoiding fracture Fracture instability and Fatigue in wedge splitting tests by means of numerical simulations XIV International Conference on Computational Plasticity. Fundamentals and Applications
More informationArchetype-Blending Multiscale Continuum Method
Archetype-Blending Multiscale Continuum Method John A. Moore Professor Wing Kam Liu Northwestern University Mechanical Engineering 3/27/2014 1 1 Outline Background and Motivation Archetype-Blending Continuum
More informationEXTENDED ABSTRACT. Dynamic analysis of elastic solids by the finite element method. Vítor Hugo Amaral Carreiro
EXTENDED ABSTRACT Dynamic analysis of elastic solids by the finite element method Vítor Hugo Amaral Carreiro Supervisor: Professor Fernando Manuel Fernandes Simões June 2009 Summary The finite element
More informationAn Evaluation of Simplified Methods to Compute the Mechanical Steady State
An Evaluation of Simplified Methods to Compute the Mechanical Steady State T. Herbland a,b, G. Cailletaud a, S. Quilici a, H. Jaffal b, M. Afzali b a Mines Paris Paris Tech, CNRS UMR 7633, BP 87, 91003
More informationNUMERICAL ALGORITHMS AND MESH DEPENDENCE IN THE DISTURBED STATE CONCEPT
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, VOL. 40, 3059 3083 (1997) NUMERICAL ALGORITHMS AND MESH DEPENDENCE IN THE DISTURBED STATE CONCEPT CHANDRA S. DESAI, * CEMAL BASARAN AND WU ZHANG
More informationINVESTIGATION OF MESH DEPENDENCE OF STOCHASTIC SIMULATIONS OF QUASIBRITTLE FRACTURE
9th International Conference on Fracture Mechanics of Concrete and Concrete Structures FraMCoS-9 V. Saouma, J. Bolander, and E. Landis (Eds) DOI 10.21012/FC9.019 INVESTIGATION OF MESH DEPENDENCE OF STOCHASTIC
More informationCohesive Band Model: a triaxiality-dependent cohesive model inside an implicit non-local damage to crack transition framework
University of Liège Aerospace & Mechanical Engineering MS3: Abstract 131573 - CFRAC2017 Cohesive Band Model: a triaxiality-dependent cohesive model inside an implicit non-local damage to crack transition
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 informationElsevier Editorial System(tm) for Computer Methods in Applied Mechanics and Engineering Manuscript Draft
Elsevier Editorial System(tm) for Computer Methods in Applied Mechanics and Engineering Manuscript Draft Manuscript Number: Title: Modelling of Strain Softening Materials Based on Equivalent Damage Force
More informationVORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS
The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS K.
More informationFinite element analysis of diagonal tension failure in RC beams
Finite element analysis of diagonal tension failure in RC beams T. Hasegawa Institute of Technology, Shimizu Corporation, Tokyo, Japan ABSTRACT: Finite element analysis of diagonal tension failure in a
More informationPlane Strain Test for Metal Sheet Characterization
Plane Strain Test for Metal Sheet Characterization Paulo Flores 1, Felix Bonnet 2 and Anne-Marie Habraken 3 1 DIM, University of Concepción, Edmundo Larenas 270, Concepción, Chile 2 ENS - Cachan, Avenue
More informationMultiscale modeling of failure in ABS materials
Institute of Mechanics Multiscale modeling of failure in ABS materials Martin Helbig, Thomas Seelig 15. International Conference on Deformation, Yield and Fracture of Polymers Kerkrade, April 2012 Institute
More informationMASONRY MICRO-MODELLING ADOPTING A DISCONTINUOUS FRAMEWORK
MASONRY MICRO-MODELLING ADOPTING A DISCONTINUOUS FRAMEWORK J. Pina-Henriques and Paulo B. Lourenço School of Engineering, University of Minho, Guimarães, Portugal Abstract Several continuous and discontinuous
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 informationGeometric nonlinear sensitivity analysis for nonparametric shape optimization with non-zero prescribed displacements
0 th World Congress on Structural and Multidisciplinary Optimization May 9-24, 203, Orlando, Florida, USA Geometric nonlinear sensitivity analysis for nonparametric shape optimization with non-zero prescribed
More information3D ANALYSIS OF H-M COUPLED PROBLEM WITH ZERO-THICKNESS INTERFACE ELEMENTS APPLIED TO GEOMECHANICS
Environmental 3D analysis of H-M and Geosciences coupled problem with zero-thickness interface elements applied to Geomechanics XIII International Conference on Computational Plasticity. Fundamentals and
More informationA MIXED TRIANGULAR ELEMENT FOR THE ELASTO-PLASTIC ANALYSIS OF KIRCHHOFF PLATES
COMPUTATIONAL MECHANICS New Trends and Applications S. Idelsohn, E. Oñate and E. Dvorkin (Eds.) c CIMNE, Barcelona, Spain 1998 A MIXED TRIANGULAR ELEMENT FOR THE ELASTO-PLASTIC ANALYSIS OF KIRCHHOFF PLATES
More informationA BEM implicit formulation applied to damage mechanics models
A BEM implicit formulation applied to damage mechanics models A. S. ~otta', W. S. ~enturini' & A. ~enallal~ 'S& Carlos School oj'engineering, University of $60 Paulo, Brazil 2 Laboratoire de Me'canique
More informationNeuro -Finite Element Static Analysis of Structures by Assembling Elemental Neuro -Modelers
Neuro -Finite Element Static Analysis of Structures by Assembling Elemental Neuro -Modelers Abdolreza Joghataie Associate Prof., Civil Engineering Department, Sharif University of Technology, Tehran, Iran.
More informationFluid driven cohesive crack propagation in quasi-brittle materials
Fluid driven cohesive crack propagation in quasi-brittle materials F. Barpi 1, S. Valente 2 Department of Structural and Geotechnical Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129
More informationCheckerboard instabilities in topological shape optimization algorithms
Checkerboard instabilities in topological shape optimization algorithms Eric Bonnetier, François Jouve Abstract Checkerboards instabilities for an algorithm of topological design are studied on a simple
More informationNonlinear Static - 1D Plasticity - Isotropic and Kinematic Hardening Walla Walla University
Nonlinear Static - 1D Plasticity - Isotropic and Kinematic Hardening Walla Walla University Professor Louie L. Yaw c Draft date January 18, 2017 Copyright is reserved. Individual copies may be made for
More informationLocking in the incompressible limit: pseudo-divergence-free element free Galerkin
Locking in the incompressible limit: pseudodivergencefree element free Galerkin Yolanda Vidal Pierre Villon and Antonio Huerta Departament de Matemàtica Aplicada III, ES de Ingenieros de Caminos, Canales
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 informationComputational Inelasticity FHLN05. Assignment A non-linear elasto-plastic problem
Computational Inelasticity FHLN05 Assignment 2017 A non-linear elasto-plastic problem General instructions A written report should be submitted to the Division of Solid Mechanics no later than October
More informationPIANO SOUNDBOARD UNDER PRESTRESS: A NUMERICAL APPROACH
9th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, -7 SEPTEMBER 7 PIANO SOUNDBOARD UNDER PRESTRESS: A NUMERICAL APPROACH PACS: 43.75.Mn Mamou-Mani, Adrien ; Frelat, Joël ; Besnainou, Charles Institut Jean
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 informationChapter 3 Variational Formulation & the Galerkin Method
Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation
More informationOn the linearization of the microplane model
MECHANICS OF COHESIVE-FRICTIONAL MATERIALS Mech. Cohes.-Frict. Mater. 3, 343 364 (1998) On the linearization of the microplane model Ellen Kuhl and Ekkehard Ramm Institute of Structural Mechanics, University
More informationTIME-DEPENDENT ANALYSIS OF PARTIALLY PRESTRESSED CONTINUOUS COMPOSITE BEAMS
2nd Int. PhD Symposium in Civil Engineering 998 Budapest IME-DEPENDEN ANAYSIS OF PARIAY PRESRESSED CONINUOUS COMPOSIE BEAMS M. Sar and Assoc. Prof. J. apos 2 Slova echnical University, Faculty of Civil
More informationarxiv: v2 [cond-mat.mtrl-sci] 2 Jan 2009
On a damage-plasticity approach to model concrete failure arxiv:8.776v [cond-mat.mtrl-sci] Jan 9 Abstract Peter Grassl Department of Civil Engineering University of Glasgow, Glasgow, United Kingdom Email:
More informationMicrostructural effect on the cavity expansion in a soil cylinder
Microstructural effect on the cavity expansion in a soil cylinder J.D. Zhao, D.C. Sheng & S.W. Sloan Centre for Geotechnical and Materials Modelling, School of Engineering University of Newcastle, Callaghan,
More informationBRIDGES BETWEEN DAMAGE AND FRACTURE MECHANICS
BRIDGES BETWEEN DAMAGE AND FRACTURE MECHANICS Jacky Mazars and Gilles Pijaudier-Cabot,,. Laboratoire de Mecanique et Technologie - Ecole Normale Supeneure 94235 Cachan-France Abstract Fracture mechanics
More informationMIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING
144 MIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING J. N. Reddy* and Chen-Shyh-Tsay School of Aerospace, Mechanical and Nuclear Engineering, University of Oklahoma, Norman, Oklahoma The paper describes
More informationNonlinear Analysis of Reinforced Concrete Shells Subjected to Impact Loads
Transactions of the 7 th International Conference on Structural Mechanics in Reactor Technology (SMiRT 7) Prague, Czech Republic, August 7, 00 Paper # J0- Nonlinear Analysis of Reinforced Concrete Shells
More informationSome Aspects of a Discontinuous Galerkin Formulation for Gradient Plasticity at Finite Strains
Some Aspects of a Discontinuous Galerkin Formulation for Gradient Plasticity at Finite Strains Andrew McBride 1 and B. Daya Reddy 1,2 1 Centre for Research in Computational and Applied Mechanics, University
More informationComputational Inelasticity FHLN05. Assignment A non-linear elasto-plastic problem
Computational Inelasticity FHLN05 Assignment 2018 A non-linear elasto-plastic problem General instructions A written report should be submitted to the Division of Solid Mechanics no later than November
More informationA plane stress softening plasticity model for orthotropic materials Lourenco, P.B.; de Borst, R.; Rots, J.G.
A plane stress softening plasticity model for orthotropic materials Lourenco, P.B.; de Borst, R.; Rots, J.G. Published in: International Journal for Numerical Methods in Engineering DOI: 10.1002/(SICI)1097-0207(19971115)40:213.0.CO;2-0
More informationDomain optimisation using Trefftz functions application to free boundaries
Computer Assisted Mechanics and Engineering Sciences, 4: 317 326, 1997. Copyright c 1997 by Polska Akademia Nauk Domain optimisation using Trefftz functions application to free boundaries Mohamed Bouberbachene,
More informationA FINITE ELEMENT STUDY OF ELASTIC-PLASTIC HEMISPHERICAL CONTACT BEHAVIOR AGAINST A RIGID FLAT UNDER VARYING MODULUS OF ELASTICITY AND SPHERE RADIUS
Proceedings of the International Conference on Mechanical Engineering 2009 (ICME2009) 26-28 December 2009, Dhaka, Bangladesh ICME09- A FINITE ELEMENT STUDY OF ELASTIC-PLASTIC HEMISPHERICAL CONTACT BEHAVIOR
More informationSTRAIN ASSESSMENT USFOS
1 STRAIN ASSESSMENT IN USFOS 2 CONTENTS: 1 Introduction...3 2 Revised strain calculation model...3 3 Strain predictions for various characteristic cases...4 3.1 Beam with concentrated load at mid span...
More informationA novel approach to predict the growth rate of short cracks under multiaxial loadings
A novel approach to predict the growth rate of short cracks under multiaxial loadings F. Brugier 1&2, S. Pommier 1, R. de Moura Pinho 2, C. Mary 2 and D. Soria 2 1 LMT-Cachan, ENS Cachan / CNRS / UPMC
More informationNon-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 17, 2017, Lesson 5
Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 17, 2017, Lesson 5 1 Politecnico di Milano, February 17, 2017, Lesson 5 2 Outline
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 informationMechanics of Materials Primer
Mechanics of Materials rimer Notation: A = area (net = with holes, bearing = in contact, etc...) b = total width of material at a horizontal section d = diameter of a hole D = symbol for diameter E = modulus
More information