CDPM2: A damage-plasticity approach to modelling the failure of concrete
|
|
- Archibald Smith
- 6 years ago
- Views:
Transcription
1 CDPM2: A damage-plasticity approach to modelling the failure of concrete Peter Grassl 1, Dimitrios Xenos 1, Ulrika Nyström 2, Rasmus Rempling 2, Kent Gylltoft 2 arxiv: v1 [cond-mat.mtrl-sci] 26 Jul School of Engineering, University of Glasgow, Glasgow, UK 2 Department of Civil and Environmental Engineering, Chalmers University of Technology, Göteborg, Sweden Corresponding author: peter.grassl@glasgow.ac.uk, Phone: , Fax: Abstract A constitutive model based on the combination of damage mechanics and plasticity is developed to analyse the failure of concrete structures. The aim is to obtain a model, which describes the important characteristics of the failure process of concrete subjected to multiaxial loading. This is achieved by combining an effective stress based plasticity model with a damage model based on plastic and elastic strain measures. The model response in tension, uni-, bi- and triaxial compression is compared to experimental results. The model describes well the increase in strength and displacement capacity for increasing confinement levels. Furthermore, the model is applied to the structural analyses of tensile and compressive failure. Keywords: concrete; constitutive model; plasticity; damage mechanics; fracture; mesh dependence 1 Introduction Concrete is a strongly heterogeneous material, which exhibits a complex nonlinear mechanical behaviour. Failure in tension and low confined compression is characterised by softening which is defined as decreasing stress with increasing deformations. This softening response is accompanied by a reduction of the unloading stiffness of concrete, and irreversible (permanent) deformations, which are localised in narrow zones often called cracks or shear bands. On the other hand, the behaviour of concrete subjected to high confined compression is characterised by a ductile hardening response; that is, increasing stress with increasing deformations. These phenomena should be considered in a constitutive model for analysing the multiaxial behaviour of concrete structures. 1
2 There are many constitutive models for the nonlinear response of concrete proposed in the literature. Commonly used frameworks are plasticity, damage mechanics and combinations of plasticity and damage mechanics. Stress-based plasticity models are useful for the modelling of concrete subjected to triaxial stress states, since the yield surface corresponds at a certain stage of hardening to the strength envelope of concrete (Leon, 1935; Willam and Warnke, 1974; Pramono and Willam, 1989; Etse and Willam, 1994; Menétrey and Willam, 1995; Pivonka, 2001; Grassl et al., 2002; Papanikolaou and Kappos, 2007; Červenka and Papanikolaou, 2008; Folino and Etse, 2012). Furthermore, the strain split into elastic and plastic parts represents realistically the observed deformations in confined compression, so that unloading and path-dependency can be described well. However, plasticity models are not able to describe the reduction of the unloading stiffness that is observed in experiments. Conversely, damage mechanics models are based on the concept of a gradual reduction of the elastic stiffness (Kachanov, 1980; Mazars, 1984; Ortiz, 1985; Resende, 1987; Mazars and Pijaudier-Cabot, 1989; Carol et al., 2001; Tao and Phillips, 2005; Voyiadjis and Kattan, 2009). For strain-based isotropic damage mechanics models, the stress evaluation procedure is explicit, which allows for a direct determination of the stress state, without an iterative calculation procedure. Furthermore, the stiffness degradation in tensile and low confined compressive loading observed in experiments can be described. However, isotropic damage mechanics models are often unable to describe irreversible deformations observed in experiments and are mainly limited to tensile and low confined compression stress states. On the other hand, combinations of isotropic damage and plasticity are widely used for modelling both tensile and compressive failure and many different models have been proposed in the literature (Ju, 1989; Lee and Fenves, 1998; Jason et al., 2006; Grassl and Jirásek, 2006; Nguyen and Houlsby, 2008; Nguyen and Korsunsky, 2008; Voyiadjis et al., 2008; Grassl, 2009; Sánchez et al., 2011; Valentini and Hofstetter, 2012). One popular class of damage-plastic models relies on a combination of stress-based plasticity formulated in the effective (undamaged) stress space combined with a strain based damage model. The combined damage-plasticity model recently developed by Grassl and Jirásek (2006); Grassl and Jirásek (2006a) belongs to this group. This model, called here Concrete Damage Plasticity Model 1 (CDPM1), is characterised by a very good agreement with a wide range of experimental results of concrete subjected to multiaxial stress states. Furthermore, it has been used in structural analysis in combination with techniques to obtain mesh-independent results and has shown to be robust (Grassl and Jirásek, 2006a; Valentini and Hofstetter, 2012). However, CDPM1 is based on a single damage parameter for both tension and compression. This is sufficient for monotonic loading with unloading, but is not suitable for modelling the transition from tensile to compressive failure realistically. When the model was proposed, this limitation was already noted and a generalisation to isotropic formulations with several damage parameters was recommended. In the present work, CDPM1 is revisted to address this issue by proposing separate damage variables for tension and compression. The introduction of two isotropic damage variables for tension and compression was motivated by the work of Mazars (1984); Ortiz (1985); Fichant et al. (1999). Secondly, in CDPM1, a perfect plastic response in the nom- 2
3 inal post-peak regime is assumed for the plasticity part and damage is determined by a function of the plastic strain. For the nonlocal version of CDPM1 presented in Grassl and Jirásek (2006a), this perfect-plastic response resulted in mesh-dependent plastic strain profiles, although the overall load-displacement response was mesh-independent. Already in Grassl and Jirásek (2006a), it was suggested that the plastic strain profile could be made mesh-independent by introducing hardening in the plasticity model for the nominal post-peak regime. In the present model, the damage functions for tension and compression depend on both plastic and elastic strain components. Furthermore, hardening is introduced in the nominal post-peak regime. With these extensions, the damage laws can be analytically related to chosen stress-inelastic strain relations, which simplifies the calibration procedure. The extension to hardening is based on recent 1D damage-plasticity model developments in Grassl (2009), which are here for the first time applied to a 3D model. The present damage-plasticity model for concrete failure is an augmentation of CDPM1. Therefore, the model is called here CDPM2. The aim of this article is to present in detail the new phenomenological model and to demonstrate that this model is capable of describing the influence of confinement on strength and displacement capacity, the presence of irreversible displacements and the reduction of unloading stiffness, and the transition from tensile to compressive failure realistically. Furthermore, it will be shown, by analysing structural tests, that CDPM2 is able to describe concrete failure mesh independently. 2 Damage-plasticity constitutive model 2.1 General framework The damage plasticity constitutive model is based on the following stress-strain relationship: σ = (1 ω t ) σ t + (1 ω c ) σ c (1) where σ t and σ c are the positive and negative parts of the effective stress tensor σ, respectively, and ω t and ω c are two scalar damage variables, ranging from 0 (undamaged) to 1 (fully damaged). The effective stress σ is defined as σ = D e : (ε ε p ) (2) where D e is the elastic stiffness tensor based on the elastic Young s modulus E and Poisson s ratio ν, ε is the strain tensor and ε p is the plastic strain tensor. The positive and negative parts of the effective stress σ in (1) are determined from the principal effective stress σ p as σ pt = σ p + and σ pc = σ p, where + and are positive and negative part operators, respectively, defined as x + = max (0, x) and x = min (0, x). For instance, for a combined tensile and compressive stress state with principal effective stress components σ p = ( σ, 0.2 σ, 0.1 σ) T, the positive and negative principal stresses are 3
4 σ pt = (0, 0.2 σ, 0.1 σ) T and σ pc = ( σ, 0, 0) T, respectively. The plasticity model is based on the effective stress, which is independent of damage. The model is described by the yield function, the flow rule, the evolution law for the hardening variable and the loading-unloading conditions. The form of the yield function is f p ( σ, κ p ) = F ( σ, q h1, q h2 ) (3) where q h1 (κ p ) and q h2 (κ p ) are dimensionless functions controlling the evolution of the size and shape of the yield surface. The flow rule is ε p = λ g p σ ( σ, κ p) (4) where ε p is the rate of the plastic strain, λ is the rate of the plastic multiplier and g p is the plastic potential. The rate of the hardening variable κ p is related to the rate of the plastic strain by an evolution law. The loading-unloading conditions are f p 0, λ 0, λfp = 0 (5) A detailed description of the individual components of the plasticity part of the model are discussed in Section 2.2. The damage part of the model is described by the damage loading functions, loading unloading conditions and the evolution laws for damage variables for tension and compression. For tensile damage, the main equations are For compression, they are f dt = ε t ( σ) κ dt (6) f dt 0, κ dt 0, κ dt f dt = 0 (7) ω t = g dt (κ dt, κ dt1, κ dt2 ) (8) f dc = α c ε c ( σ) κ dc (9) f dc 0, κ dc 0, κ dc f dc = 0 (10) ω c = g dc (κ dc, κ dc1, κ dc2 ) (11) Here, f dt and f dc are the loading functions, ε t ( σ) and ε c ( σ) are the equivalent strains and κ dt, κ dt1, κ dt2, κ dc, κ dc1 and κ dc2 are damage history variables. Furthermore, α c is a variable that distinguishes between tensile and compressive loading. A detailed description of the variables is given in Section
5 2.2 Plasticity part The plasticity part of the model is formulated in a three-dimensional framework with a pressure-sensitive yield surface, hardening and non-associated flow. The main components are the yield function, the flow rule, the hardening law and the evolution law for the hardening variable Yield function The yield surface is described in terms of the cylindrical coordinates in the principal effective stress space (Haigh-Westergaard coordinates), which are the volumetric effective stress σ V = I 1 (12) 3 the norm of the deviatoric effective stress ρ = 2J 2 (13) and the Lode angle ( θ = 1 arccos 3 ) 3 J J 3/2 2 The foregoing definitions use the first invariant (14) I 1 = σ : δ = σ ij δ ij (15) of the effective stress tensor σ, and the second and third invariants J 2 = 1 2 s : s = 1 2 s2 : δ = 1 s 2 ij s ij (16) J 3 = 1 3 s3 : δ = 1 s 3 ij s jk s ki (17) of the deviatoric effective stress tensor s = σ δi 1 /3. The yield function { ( f p ( σ V, ρ, θ; ρ κ p ) = [1 q h1 (κ p )] + σ ) 2 V 3 ρ + 6fc f c 2 f c [ ρ + m 0 qh1(κ 2 p )q h2 (κ p ) r(cos θ) + σ ] V q 6fc f h1(κ 2 p )qh2(κ 2 p ) c } 2 (18) depends on the effective stress (which enters in the form of cylindrical coordinates) and on the hardening variable κ p (which enters through the dimensionless variables q h1 and q h2 ). Parameter f c is the uniaxial compressive strength. For q h2 = 1, the yield function is identical to the one of CDPM1. 5
6 Figure 1: The evolution of the meridional section of the yield surface during hardening. The meridians of the yield surface f p = 0 are parabolic, and the deviatoric sections change from triangular shapes at low confinement to almost circular shapes at high confinement. The shape of the deviatoric section is controlled by the function r(cos θ) = 4(1 e 2 ) cos 2 θ + (2e 1) 2 2(1 e 2 ) cos θ + (2e 1) 4(1 e 2 ) cos 2 θ + 5e 2 4e (19) proposed by Willam and Warnke (1974). The calibration of the eccentricity parameter e is described in Jirásek and Bažant (2002) and in section 5. The friction parameter m 0 is given by m 0 = 3 (f c 2 ft 2 ) e (20) f c f t e + 1 where f t is the tensile strength. The shape and evolution of the yield surface is controlled by the variables q h1 and q h2 (Figs. 1 and 2). If the two variables q h1 and q h2 in (18) are set equal to one and the resulting yield function is set equal to zero, the failure surface [ 3 ρ 2 ρ + m 0 r(cos 2 θ) + σ ] V 1 = 0 (21) 6fc f 2 c is obtained, which was originally proposed by Menétrey and Willam (1995). f c Flow rule In the present model, the flow rule in (4) is non-associated, which means that the yield function f p and the plastic potential g p do not coincide and, therefore, the direction of 6
7 Figure 2: The evolution of the deviatoric section of the yield surface during hardening for a constant volumetric stress of σ V = f c /3. the plastic flow g p / σ is not normal to the yield surface. The plastic potential is given as { ( ρ g p ( σ V, ρ; κ p ) = [1 q h1 (κ p )] + σ ) 2 } 2 V 3 ρ + 6fc f c 2 f c ( + qh1(κ 2 m0 ρ p ) + m ) (22) g( σ V, κ p ) 6fc f c where m g ( σ V, κ p ) = A g (κ p ) B g (κ p ) f c exp σ V q h2 (κ p )f t /3 B g (κ p ) f c (23) is a variable controlling the ratio of volumetric and deviatoric plastic flow. Here, A g (κ p ) and B g (κ p ), which depend on q h2 (κ p ), are derived from assumptions on the plastic flow in uniaxial tension and compression in the post-peak regime. The derivation of these two variables is illustrated in the following paragraphs. Here, the notation m g p is introduced. In the principal stress space, the plastic flow σ tensor m has three components, m 1, m 2 and m 3 associated with the three principal stress components. The flow rule (4) is split into a volumetric and a deviatoric part, i.e., the gradient of the plastic potential is decomposed as m = g σ = g σ V σ V σ + g ρ ρ σ (24) Taking into account that σ V / σ = δ/3 and ρ/ σ = s/ ρ, restricting attention to the post-peak regime (in which q h1 = 1) and differentiating the plastic potential (22), we 7
8 rewrite equation (24) as m = g σ = m ( g δ m ) 0 s (25) σ V 3f c f c 6 ρ f c Experimental results for concrete loaded in uniaxial tension indicate that the strains perpendicular to the loading direction are elastic in the softening regime. Thus, the plastic strain rate in these directions should be equal to zero (m 2 = m 3 = 0). Under uniaxial tension, the effective stress state in the post-peak regime is characterised by σ 1 = f t q h2, σ 2 = σ 3 = 0, σ V = f t q h2 /3, s 1 = 2f t q h2 /3, s 2 = s 3 = f t q h2 /3 and ρ = 2/3f t q h2. Substituting this into (25) and enforcing the condition m 2 = m 3 = 0, we obtain an equation from which m g σv =f σ tq h2 /3 = 3f tq h2 + m 0 V f c 2 (26) In uniaxial compressive experiments, a volumetric expansion is observed in the softening regime. Thus, the inelastic lateral strains are positive while the inelastic axial strain is negative. In the present approach, a constant ratio D f = m 2 /m 1 = m 3 /m 1 between lateral and axial plastic strain rates in the softening regime is assumed. The effective stress state at the end of hardening under uniaxial compression is characterised by σ 1 = f c q h2, σ 2 = σ 3 = 0, σ V = f c q h2 /3, s 1 = 2f c q h2 /3, s 2 = s 3 = f c q h2 /3 and ρ = 2/3f c q h2. Substituting this into (25) and enforcing the condition m 2 = m 3 = D f m 1, we get an equation from which m g σv = f σ cq h2 /3 = 2D f 1 ( 3q h2 + m ) 0 V D f Substituting the specific expression for m g / σ V constructed by differentiation of (23) into (26) and (27), we obtain two equations from which parameters (27) A g = 3f tq h2 + m 0 f c 2 (q h2 /3) (1 + f t /f c ) B g = ln A g ln (2D f 1) ln (3q h2 + m 0 /2) + ln (D f + 1) (28) (29) can be computed. The gradient of the dilation variable m g in (23) decreases with increasing confinement. The limit σ V corresponds to purely deviatoric flow. As in CDPM1, the plastic potential does not depend on the third Haigh-Westergaard coordinate (Lode angle θ), which increases the efficiency of the implementation and the robustness of the model. 8
9 Figure 3: The two hardening laws q h1 (solid line) and q h2 (dashed line) Hardening law The dimensionless variables q h1 and q h2 that appear in (18), (22) and (23) are functions of the hardening variable κ p. They control the evolution of the size and shape of the yield surface and plastic potential. The first hardening law q h1 is q h1 (κ p ) = { qh0 + (1 q h0 ) ( κ 3 p 3κ 2 p + 3κ p ) Hp ( κ 3 p 3κ 2 p + 2κ p ) if κp < 1 1 if κ p 1 (30) The second hardening law q h2 is given by q h2 (κ p ) = { 1 if κp < H p (κ p 1) if κ p 1 (31) The initial inclination of the hardening curve q h1 at κ p = 0 is positive and finite, and the inclinations of q h1 on the left of κ p = 1 and q h2 on the right of κ p = 1 are equal to H p, as depicted in Fig. 3. For H p = 0, the hardening law reduces to the one proposed in Grassl and Jirásek (2006) Hardening variable The evolution law for the hardening variable, κ p = ε p ( 2 λ m ( 2 2 cos θ) = 2 cos θ) (32) x h ( σ V ) x h ( σ V ) 9
10 sets the rate of the hardening variable equal to the norm of the plastic strain rate scaled by a hardening ductility measure A h (A h B h ) exp ( R h ( σ V )/C h ) if R h ( σ V ) 0 x h ( σ V ) = (33) E h exp(r h ( σ V )/F h ) + D h if R h ( σ V ) < 0 For pure volumetric stress states, θ in (32) is set to zero. The dependence of the scaling factor x h on the volumetric stress σ V is constructed such that the model response is more ductile under compression. The variable R h ( σ V ) = σ V f c 1 3 (34) is a linear function of the volumetric effective stress. Model parameters A h, B h, C h and D h are calibrated from the values of strain at peak stress under uniaxial tension, uniaxial compression and triaxial compression, whereas the parameters E h and F h are determined from the conditions of a smooth transition between the two parts of equation (33) at R h = 0: E h = B h D h (35) F h = (B h D h ) C h A h B h (36) This definition of the hardening variable is identical to the one in CDPM1 described in Grassl and Jirásek (2006), where the calibration procedure for this part of the model is described. 2.3 Damage part Damage is initiated when the maximum equivalent strain in the history of the material reaches the threshold ε 0 = f t /E. For uniaxial tension only, the equivalent strain could be chosen as ε = σ t /E, where σ t is the effective uniaxial tensile stress. Thus, damage initiation would be linked to the axial elastic strain. However, for general triaxial stress states a more advanced equivalent strain expression is required, which predicts damage initiation when the strength envelope is reached. This expression is determined from the yield surface (f p = 0) by setting q h1 = 1 and q h2 = ε/ε 0. From this quadratic equation for ε, the equivalent strain is determined as ε = ε 0m 0 2 ( ρ r (cos θ) + σ ) V + 6fc f c ( ε 2 0m 2 0 ρ 4 6fc r (cos θ) + σ V f c ) 2 + 3ε2 0 ρ 2 For uniaxial tension, the effective stress state is defined as σ 1 = σ t, σ 2 = σ 3 = 0, σ V = σ t /3, s 1 = 2 σ t /3, s 2 = s 3 = σ t /3, ρ = 2/3 σ t and r(cos θ) = 1/e. Setting this into (37) 2f 2 c (37) 10
11 and using the definition of m 0 in (20) gives ε = ε 0 σ t f t = σ t /E (38) which is suitable equivalent strain for modelling tensile failure. For uniaxial compression, the effective stress state is defined as σ 1 = σ c, σ 2 = σ 3 = 0, σ V = σ c /3, s 1 = 2/3 σ c, s 2 = s 3 = 1/3 σ c, ρ = 2/3 σ c, and r(cos θ) = 1. Here, σ c is the magnitude of the effective compressive stress. Setting this into (37), the equivalent strain is ε = σ cε 0 f c = σ cf t Ef c (39) If σ c = (f c /f t ) σ t, the equivalent strain is again equal to the axial elastic strain component in uniaxial tension. Consequently, the equivalent strain definition in (37) is suitable for both tension and compression, which is very convenient for relating the damage variables in tension and compression to stress-inelastic strain curves. The damage variables ω t and ω c in (1) are determined so that a prescribed stress-inelastic strain relation in uniaxial tension is obtained. Since, the damage variables are evaluated for general triaxial stress states, the inelastic strain in uniaxial tension has to be expressed by suitable scalar history variables, which are obtained from total and plastic strain components. To illustrate the choice of these components, a 1D damage-plastic stressstrain law of the form σ = (1 ω) σ = (1 ω) E (ε ε p ) (40) is considered. Here, ω is the damage variable. This law can also be written as σ = E {ε [ε p + ω (ε ε p )]} = E (ε ε i ) (41) where ε i is the inelastic strain which is subtracted from the total strain. The geometrical interpretation of the inelastic strain and its split for monotonic uniaxial tension, linear hardening plasticity and linear damage evolution are shown in Fig. 4. Furthermore, the way how the hardening influences damage and plasticity dissipation has been discussed in Grassl (2009). The part ω (ε ε p ) is reversible and ε p is irreversible. The damage variable is chosen, so that a softening law is obtained, which relates the stress to the inelastic strain, which is written here in generic form as σ = f s (ε i ) (42) Setting (41) equal with (42) allows for determining the damage variable ω. However, the inelastic strain ε i in (41) and (42) needs to be expressed by history variables, so that the expression for the damage variable can be used for non-monotonic loading. Furthermore, to be able to describe also the influence of multiaxial stress states on the damage evolution, the inelastic strain in (42) is replaced by different history variables than the inelastic strain in (41). The choice of the history variables for tension and compression 11
12 Figure 4: Geometrical meaning of the inelastic strain ε i for the combined damageplasticity model. The inelastic strain is composed of reversible ω (ε ε p ) and irreversible ε p parts. The dashed lines represent elastic unloading with the same stiffness as the initial elastic loading. 12
13 is explained in sections and History variables for tension The tensile damage variable ω t in (1) is defined by three history variables κ dt, κ dt1 and κ dt2. The variable κ dt is used in the definition of the inelastic strain in (41), while κ dt1 and κ dt2 enter the definition of the inelastic strain in (42). The history variable κ dt is determined from ε t using (6) and (7). Here, ε t is given implicitly in incremental form by ε t = ε (43) with ε given in (37). For κ dt1, the inelastic strain component related the plastic strain ε p is replaced by 1 ε κ dt1 = p if κ dt > 0 and κ dt > ε 0 x s (44) 0 if κ dt = 0 or κ dt < ε 0 Here, the pre-peak plastic strains do not contribute to this history variable, since κ dt1 is only nonzero, if κ dt > ε 0. Finally, the third history variable is related to κ dt as κ dt2 = κ dt x s (45) In (44) and (45), x s is a ductility measure, which describes the influence of multiaxial stress states on the softening response, see Sec History variables for compression The compression damage variable ω c is also defined by three history variables κ dc, κ dc1 and κ dc2. Analogous to the tensile case, the variable κ dc is used in the definition of the inelastic strain in (41), while κ dc1 and κ dc2 enter the definition of the equivalent strain in (42). In addition, a variable α c is introduced which distinguishes tensile and compressive stresses. It has the form 3 σ pci ( σ pti + σ pci ) α c = (46) σ p 2 i=1 where σ pti and σ pci are the components of the compressive and tensile part of the principal effective stresses, respectively, which were previously used for the general stress strain law in (1). The variable α c varies from 0 for pure tension to 1 for pure compression. For instance, for the mixed tensile compressive effective stress state σ p = { σ, 0.2 σ, 0.1 σ}, considered in Sec. 2.1, the variable is α c = The history variable κ dc is determined from ε c using (9) and (10), where, analogous to 13
14 the tensile case, the ε c is specified implicitly by ε c = α c ε (47) The other two history variables are α c β c ε κ dc1 = p if κ dt > 0 κ dt > ε 0 x s (48) 0 if κ dt = 0 κ dt < ε 0 and In (48), the factor β c is κ dc2 = κ dc x s (49) β c = f tq h2 2/3 ρ (50) 1 + 2Df 2 This factor provides a smooth transition from pure damage to damage-plasticity softening processes, which can occur during cyclic loading, as described in section Damage variables for bilinear softening With the history variables defined in the previous two sections, the damage variables for tension and compression are determined. The form of these damage variables depends on the type of softening law considered. For bilinear softening used in the present study, the stress versus inelastic strain in the softening regime is σ = f t f t σ 1 ε f1 ε i if 0 < ε i ε f1 σ 1 σ 1 ε f ε f1 (ε i ε f1 ) if ε f1 < ε i ε f 0 if ε f ε i where ε f is the inelastic strain threshold at which the uniaxial stress is equal to zero and ε f1 is the threshold where the uniaxial stress is equal to σ 1 as shown in Fig. 5. Furthermore, ε i is the inelastic strain in the post-peak regime only. Since damage is irreversible, the inelastic strain ε i in (51) is expressed by irreversible damage history variables as (51) ε i = κ dt1 + ω t κ dt2 (52) Furthermore, the term ε ε p in (40) is replaced by κ dt, which gives σ = (1 ω t ) Eκ dt (53) 14
15 Figure 5: Bilinear softening. Setting (51) with (52) equal to (53), and solving for ω t gives ω t = (Eκ dt f t )ε f1 (σ 1 f t )κ dt1 Eκ dt ε f1 + (σ 1 f t )κ dt2 Eκ dt (ε f ε f1 ) + σ 1 (κ dt1 ε f ) Eκ dt (ε f ε f1 ) σ 1 κ dt2 if 0 < ε i ε f1 if ε f1 < ε i ε f 0 if ε f < ε i (54) For the compressive damage variable, an evolution based on an exponential stress-inelastic strain law is used. The stress versus inelastic strain in the softening regime in compression is ( σ = f t exp ε ) i if 0 < ε i (55) ε fc where ε fc is an inelastic strain threshold which controls the initial inclination of the softening curve. The use of different damage evolution for tension and compression is one important improvement over CDPM1 as it will shown later on when the structural applications are discussed Ductility measure The history variables κ dt1, κ dt2, κ dc1 and κ dc2 in (44), (45), (48) and (49), respectively, depend on a ductility measure x s, which takes into account the influence of multiaxial stress states on the damage evolution. This ductility measure is given by x s = 1 + (A s 1) R s (56) 15
16 where R s is 6 σv if σ R s = V 0 ρ 0 if σ V > 0 and A s is a model parameter. For uniaxial compression σ V / ρ = 1/ 6, so that R s = 1 and x s = A s, which simplifies the calibration of the softening response in this case. (57) Constitutive response to cyclic loading The response of the constitutive model is illustrated by a quasi-static strain cycle (Fig. 6, solid line), before it is compared to a wide range of experimental results in the next section. The strain is increased from point 0 to point 1, where the tensile strength of the material is reached. Up to point 1, the material response is elastic-plastic with small plastic strains. With a further increase of the strain from point 1 to point 2, the effective stress part continues to increase, since H p > 0, whereas the nominal stress decreases, since the tensile damage variable ω t increases. A reverse of the strain at point 2 results in an reduction of the stress with an unloading stiffness, which is less than the elastic stiffness of an elasto-plastic model, but greater than the stiffness of an elasto-damage mechanics model, i.e. greater than the secant stiffness. At point 3, when the stress is equal to zero, a further reduction of the strain leads to a compressive response following a linear stress-strain relationship between the points 3 and 4 with the original Young s modulus E of the undamaged material. This change of stiffness is obtained by using two damage variables, ω t and ω c. At point 3, ω t > 0, but ω c = 0. Up to point 5, no further plastic strains are generated, since the hardening from point 0 to 2 has increased the elastic domain of the plasticity part, so that the yield surface is not reached. Thus, the softening from point 4 to 5 is only described by damage. Only at point 5, the plasticity surface is reached and a subsequent increase of strain results in hardening of the plasticity part, which corresponds to an increase of the effective stress. However, the nominal stress, shown in Fig. 6, decreases, since ω c increases. The continuous slopes of parts 4-5 and 5-6 are obtained, since the additional factor β c in (48) is introduced. A second reversal of the strain direction (point 6) changes the stress from compression to tension at point 7, which is again associated with a change of the stiffness. The above response is very different from the one obtained with CDPM1 with only one damage parameter, which is also shown in Fig. 6 by a dashed line. With CDPM1, the compressive response after point 3 is characterised by both a reduced stiffness and strength which would depend on the amount of damage accumulated in tension. For the case of damage equal to 1 in tension, both the strength and stiffness in compression would be zero, which is not realistic for the tension-compression transition in concrete. 16
17 1 0 stress [MPa] -1 CDPM1 CDPM strain [mm/m] Figure 6: Model response for cyclic loading with f t = 1 and f c = 3 for CDPM2 (solid line) and CDPM1 (dashed line). 3 Mesh adjusted softening modulus If the constitutive model described in the previous sections is straightaway used within the finite element method, the amount of dissipated energy might be strongly meshdependent. This mesh-dependence is caused by localisation of deformations in mesh-size dependent zones. The finer the mesh, the less energy would be dissipated. This is a well known limitation of constitutive laws with strain softening. One way to overcome this mesh-dependence is to adjust the softening modulus with respect to the element size. For the present model, this approach is applied for the tensile damage variable by replacing in the tensile damage law in (54) the strain thresholds ε f1 and ε f with w f1 /h and w f /h, respectively. Here, w f1 and w f are displacement thresholds and h is the finite element size. Thus, with this approach the damage variables for bilinear softening are ω t = (Eκ dt f t )w f1 (σ 1 f t )κ dt1 h Eκ dt w f1 + (σ 1 f t )κ dt2 h Eκ dt (w f w f1 ) + σ 1 (κ dt1 h w f ) Eκ dt (w f w f1 ) σ 1 κ dt2 h if 0 < hε i w f1 h if w f1 < hε i w f 0 if w f < hε i These expressions are used when the constitutive model is compared to experimental results in the next section. However, the evolution law for compressive damage is kept to be independent of the element size, as compressive failure is often accompanied by mesh-independent zones of localised displacements. (58) 17
18 4 Implementation The present constitutive model has been implemented within the framework of the nonlinear finite element method, where the continuous loading process is replaced by incremental time steps. In each step the boundary value problem (global level) and the integration of the constitutive laws (local level) are solved. For the boundary value problem on the global level, the usual incremental-iterative solution strategy is used, in the form of a modified Newton-Raphson iteration method. For the local problem, the updated values ( ) (n+1) of the stress and the internal variables at the end of the step are obtained by a fully implicit (backward Euler) integration of the rate form of the constitutive equations, starting from their known values ( ) (n) at the beginning of the step and applying the given strain increment ε = ε (n+1) ε (n). The integration scheme is divided into two sequential steps, corresponding to the plastic and damage parts of the model. In the plastic part, the plastic strain ε p and the effective stress σ at the end of the step are determined. In the damage part, the damage variables ω t and ω c, and the nominal stress σ at the end of the step are obtained. The implementation strategy for the local problem, described in detail in Grassl and Jirásek (2006) for CDPM1, applies to the present model as well. To improve the robustness of the model, a subincrementation scheme is employed for the integration of the plasticity part. 5 Comparison with experimental results In this section, the model response is compared to five groups of experiments reported in the literature. For each group of experiments, the physical constants Young s modulus E, Poisson s ratio ν, tensile strength f t, compressive strength f c and tensile fracture energy G Ft are adjusted to obtain a fit for the different types of concrete used in the experiments. The first four constants are model parameters. The last physical constant, G Ft, is directly related to model parameters. For the bilinear softening law in section 2.3.3, the tensile fracture energy is G Ft = f t w f1 /2 + σ 1 w f /2 (59) For σ 1 /f t = 0.3 and w f1 /w f = 0.15 (shown by Jirásek and Zimmermann (1998) to result in a good fit for concrete failure), the expression for the fracture energy reduces to G Ft = f t w f / The compressive energy is G Fc = f c ε fc l c A s, where l c is the length in which the compressive displacement are assumed to localise and A s is the ductility measure in Sec If no experimental results are available, the five constants can be determined using, for instance, the CEB-FIP Model Code (CEB, 1991). The other model parameters are set to their default values for all groups. The eccentricity constant e that controls the shape of the deviatoric section is evaluated using the formula 18
19 in Jirásek and Bažant (2002), p. 365: e = 1 + ɛ 2 ɛ, where ɛ = f t f bc f 2 bc f 2 c f 2 c f 2 t (60) where f bc is the strength in equibiaxial compression, which is estimated as f bc = 1.16f c according to the experimental results reported in Kupfer et al. (1969). Parameter q h0 is the dimensionless ratio q h0 = f c0 /f c, where f c0 is the compressive stress at which the initial yield limit is reached in the plasticity model for uniaxial compression. Its default value is q h0 = 0.3. For the hardening modulus the default value is H p = Furthermore, the default value of the parameter of the flow rule is chosen as D f = 0.85, which yields a good agreement with experimental results in uniaxial compression. The determination of parameters A h, B h, C h and D h that influence the hardening ductility measure is more difficult. The effective stress varies within the hardening regime, even for monotonic loading, so that the ratio of axial and lateral plastic strain rate is not constant. Thus, an exact relation of all four model parameters to measurable material properties cannot be constructed. In Grassl and Jirásek (2006), it has been shown that a reasonable response is obtained with parameters A h = 0.08, B h = 0.003, C h = 2 and D h = These values were also used in the present study. Furthermore, the element size h in the damage laws in Section 3 was chosen as h = 0.1 m. The first analysis is a uniaxial tensile setup with unloading. The model response is compared to the experimental results reported in Gopalaratnam and Shah (1985) (Fig. 7). The relevant model parameters for this experiment are E = 28 GPa, ν = 0.2, f c = 40 MPa, f t = 3.5 MPa, G Ft = 55 J/m 2. 4 constitutive model experiments 3 stress [MPa] strain [mm/m] Figure 7: Uniaxial tension: Model response compared to experimental results in Gopalaratnam and Shah (1985). 19
20 The next example is an uniaxial compression test with unloading, for which the model response is compared to experimental results reported in Karsan and Jirsa (1969) (Fig. 8). The model parameters are E = 30 GPa, ν = 0.2, f c = 28 MPa, f t = 2.8 MPa. Furthermore, the model constants for compression are A s = 5 and ε fc = The value of the tensile fracture energy G Ft does not influence the model response in compression, which also applies to all other compression tests considered in the following paragraphs. Therefore, only the compressive fracture energy is stated stress [MPa] constitutive model experiments strain [mm/m] Figure 8: Uniaxial compression: Model response compared to experimental results reported in Karsan and Jirsa (1969). Next, the model is compared to uniaxial and biaxial compression tests reported in Kupfer et al. (1969). For these experiments, the model parameters are set to E = 32 GPa, ν = 0.2, f c = 32.8 MPa, f t = 3.3 MPa. Furthermore, the model constants for compression are A s = 1.5 and ε fc = The comparison with experimental results is shown in Fig. 9 for uniaxial, equibiaxial and biaxial compression. For the biaxial compression case, the stress ratio of the two compressive stress components is σ 1 /σ 2 = 1/ 0.5. Furthermore, the performance of the model is evaluated for triaxial tests reported in Caner and Bažant (2000). The material parameters for this test are E = 25 GPa, ν = 0.2, f c = 45.7 MPa, f t = 4.57 MPa. Furthermore, the model constants for compression are A s = 15 and ε fc = The model response is compared to experimental results presented in Figs. 10. Next, the model response in triaxial compression is compared to the experimental results reported in Imran and Pantazopoulou (1996). The material parameters for this test are E = 30 GPa, ν = 0.2, f c = 47.4 MPa, f t = 4.74 MPa. Furthermore, the model constants for compression are A s = 15 and ε fc =
21 0-10 const. -1/0 const. -1/-0.5 const. -1/-1 exp. -1/0 exp. -1/-0.5 exp. -1/-1 stress [MPa] strain [mm/m] Figure 9: Uniaxial and biaxial compression: Model response compared to experimental results reported in Kupfer et al. (1969) stress [MPa] constitutive model experiments strain [mm/m] Figure 10: Confined compression: Model response compared to experiments used in Caner and Bažant (2000). Finally, the model response in hydrostatic compression is compared to the experimental results reported in Caner and Bažant (2000). The material parameters are the same as for the triaxial test shown in Fig. 10. Overall, the agreement of the model response with the experimental results is very good. The model is able to represent the strength of concrete in tension and multiaxial compression. In addition, the strains at maximum stress in tension and compression agree well 21
22 0 const. response experiments -50 stress [MPa] strain [mm/m] Figure 11: Confined compression: Model response compared to experiments reported in Imran and Pantazopoulou (1996) constitutive model experiments stress [MPa] strain [mm/m] Figure 12: Hydrostatic compression: Model response compared to experiments reported in Caner and Bažant (2000). with the experimental results. The bilinear stress-crack opening curve that was used results in a good approximation of the softening curve in uniaxial tension and compression. With the above comparisons, it is demonstrated that CDPM2, provides, very similar to CDPM1, a very good agreement with experimental results. 22
23 6 Structural analysis The performance of the proposed constitutive model is further evaluated by structural analysis of three fracture tests. The main objective of this part of the study is to demonstrate that the structural response obtained with the model is mesh-independent. This is achieved by adjusting the softening modulus with respect to the element size (section 3). 6.1 Three point bending test The first structural example is a three-point bending test of a single-edge notched beam reported by Kormeling and Reinhardt (1982). The experiment is modelled by triangular plane strain finite elements with three mesh sizes. The geometry and loading set up is shown in Fig. 13. The input parameters are chosen as E = 20 GPa, ν = 0.2, f t = 2.4 MPa, G ft = 100 N/m, f c = 24 MPa (Grassl and Jirásek, 2006a). All other parameters are set to their default values described in section 5. For this type of analysis, local stress-strain relations with strain softening are known to result in mesh-dependent load-displacement curves. The capability of the adjustment of the softening modulus approach presented in section 3 to overcome this mesh-dependence is assessed with this test. The global response Figure 13: Three point bending test: Geometry and loading setup. thickness is 0.1 m. The notch thickness is 5 mm. The out-of-plane in the form of load-crack Mouth Opening Displacement (CMOD) is shown in Fig. 14. The local response in the form of tensile damage patterns at loading stages marked in Fig. 14 for the three meshes is shown in Fig. 15. Overall, the load-cmod curves in Fig. 14 are in good agreement with the experimental results and almost mesh independent. On the other hand, the damage zones in Fig. 15 depend on the mesh size. 23
24 coarse mesh medium mesh fine mesh experimental bounds load [kn] displacement [mm] Figure 14: Load-CMOD curves of analyses with three mesh sizes compared to the experimental bounds reported in Kormeling and Reinhardt (1982). Figure 15: Tensile damage patterns for the coarse, medium and fine mesh for the three point bending test. Black indicates a tensile damage variable of Four point shear test The second structural example is a four point shear test of a single-edge notched beam reported in Arrea and Ingraffea (1982). Again, the experiment is modelled by triangular plane strain finite elements with three different mesh sizes. The geometry and loading setup are shown in Fig. 16. The input parameters are chosen as E = 30 GPa, ν = 0.18, f t = 3.5 MPa, G ft = 140 N/m, f c = 35 MPa (Jirásek and Grassl, 2008). All other parameters are set to their default values described in section 5. The global responses of analyses and experimental results are compared in the form of load-crack Mouth Sliding Displacement (CMSD) curves in Fig. 17. Furthermore, the damage patterns for the three meshes at loading stages marked in Fig. 17 are compared to the experimental crack patterns in Fig. 18. The load-cmsd curves obtained with the three meshes are in good agreement with the 24
25 Figure 16: Four point shear test: Geometry and loading setup. The out-of-plane thickness is 0.15 m. A zero notch thickness is assumed coarse medium fine experimental bounds Load [kn/m] CMSD [mm] Figure 17: Load-CMSD curves of analyses with three mesh sizes compared to the experimental bounds reported in Arrea and Ingraffea (1982). Figure 18: Four point shear test: Tensile damage patterns for the coarse, medium and fine mesh compared to the experimental crack patterns reported in Arrea and Ingraffea (1982). Black indicates a tensile damage variable of 1. 25
26 36.8 P ε (a) [mm] (b) Figure 19: a) Geometry and loading setup of the eccentric compression test. coarse finite element mesh. b) The experimental results. The coarse mesh overestimates the load levels obtained with the medium and fine mesh. However, the two finer meshes are in good agreement. Again, the width of the damaged zone depends on the element size. Furthermore, the damage zones are influenced by the mesh orientation. In particular, for the fine mesh the damage zone follows the regular element arrangement, so that the crack is less curved than reported in the experiments. This is a well known behaviour of models using the adjustement of the softening modulus approach, which has been studied in more detail in Jirásek and Grassl (2008); Grassl and Rempling (2007). 6.3 Eccentric compression test The third structural example studies the failure of a concrete prism subjected to eccentric compression, tested by Debernardi and Taliano (2001). The geometry and loading setup are shown in Figure 19a. The specimen with a relatively great eccentricity of 36.8 mm is modeled by a thin layer of linear 3D elements to reduce the computational time compared to a full 3D analysis. Three different mesh sizes with element lengths of 7.5, 5 and 2.5 mm were chosen (see Figure 19b for the coarse mesh). The model parameters were set to E = 30 GPa, ν = 0.2, f t = 4 MPa, f c = 46 MPa, G Ft = 100 N/m, A s = 10 and ε fc = The model response in terms of the overall load versus the mean compressive strain of the compressed side obtained on the fine mesh is compared to the experimental result in Figure 20. The load capacity and the strain at peak are underestimated by the model. The overall behaviour, however, is captured 26
27 coarse mesh medium mesh fine mesh experiment load [kn] average strain [mm/m] Figure 20: Comparison of the analysis of the eccentric compression test with the experiment. well. The comparison of the load-compressive strain relations for the analyses on meshes of different sizes indicates that the description of this type of compressive failure is nearly mesh-independent. The evolution of the damage zone for the analysis on the coarse mesh is depicted in Figure 21 for the final stage of the analyses in Figure 20. On the tensile side several zones of localized damage form, whereas the failure on the compressive side is described by a diffuse damage zone. 7 Conclusions The present damage plasticity model CDPM2, which combines a stress-based plasticity part with a strain based damage mechanics model, is based on an enhancement of an already exisiting damage-plasticity model called CDPM1 (Grassl and Jirásek (2006)). Based on the work presented in this manuscript, the following conclusions can be drawn on the improvements that this constitutive model provides: 1. The model is able to describe realistically the transition from tensile to compressive failure. This is achieved by the introduction of two separate isotropic damage variables for tension and compression. 2. The model is able to reproduce stress inelastic strain relations with varying ratios of reversible and irreversible strain components. The ratio can be controlled by the hardening modulus of the plasticity part. 27
28 (a) (b) (c) Figure 21: Contour plots of the damage variable for the a) coarse, b) medium and c) fine mesh of the eccentric compression test. 3. The model gives meshindependent load-displacement curves for both tensile and compressive failure. In addition, the model response is in good agreement with experimental results for a wide range of loading from uniaxial tension to confined compression. Acknowledgements This work have been performed partially within the project Dynamic behaviour of reinforced concrete structures subjected to blast and fragment impacts, which in turn is financially sponsored by MSB - the Swedish Civil Contigencies Agency. The first author would also like to thank Prof. Bořek Patzák of the Czech Technical University for kind assistance with his finite element package OOFEM (Patzák, 1999; Patzák and Bittnar, 2001). In addition, the first and second authors acknowledge funding received from the UKs Engineering and Physical Sciences Research Council (EPSRC) under grant EP/I036427/1. References M. Arrea, A.R. Ingraffea, Mixed-mode crack propagation in mortar and concrete, Department of Structural Engineering 81-83, Cornell University, Ithaca, NY,
29 F.C. Caner, Z.P. Bažant, Microplane model M4 for concrete: II. Algorithm and calibration, Journal of Engineering Mechanics, ASCE 126 (2000) I. Carol, E. Rizzi, K.J. Willam, On the formulation of anisotropic elastic degradation. I. Theory based on a pseudo-logarithmic damage tensor rate, International Journal of Solids and Structures 38 (2001) CEB, CEB-FIP Model Code 1990, Design Code, Thomas Telford, London, P.G. Debernardi, M. Taliano, Softening behaviour of concrete prisms under eccentric compressive forces, Magazine of Concrete Research 53 (2001) G. Etse, K.J. Willam, A fracture-energy based constitutive formulation for inelastic behavior of plain concrete, Journal of Engineering Mechanics, ASCE 120 (1994) S. Fichant, C.L. Borderie, G. Pijaudier-Cabot, Isotropic and anisotropic descriptions of damage in concrete structures, Mechanics of Cohesive-Frictional Materials 4 (1999) P. Folino, G. Etse, Performance dependent model for normal and high strength concretes, International Journal of Solids and Structures 49 (2012) V.S. Gopalaratnam, S.P. Shah, Softening response of plain concrete in direct tension, ACI Journal Proceedings 82 (1985). P. Grassl, On a damage-plasticity approach to model concrete failure, Proceedings of the ICE - Engineering and Computational Mechanics 162 (2009) P. Grassl, M. Jirásek, Damage-plastic model for concrete failure, International Journal of Solids and Structures 43 (2006) P. Grassl, M. Jirásek, A plastic model with nonlocal damage applied to concrete, International Journal for Numerical and Analytical Methods in Geomechanics 30 (2006a) P. Grassl, K. Lundgren, K. Gylltoft, Concrete in compression: a plasticity theory with a novel hardening law, International Journal of Solids and Structures 39 (2002) P. Grassl, R. Rempling, Influence of volumetric-deviatoric coupling on crack prediction in concrete fracture tests, Engineering Fracture Mechanics 74 (2007) I. Imran, S.J. Pantazopoulou, Experimental study of plain concrete under triaxial stress, ACI Materials Journal 93 (1996) L. Jason, A. Huerta, G. Pijaudier-Cabot, S. Ghavamian, An elastic plastic damage formulation for concrete: Application to elementary tests and comparison with an isotropic damage model, Computer Methods in Applied Mechanics and Engineering 195 (2006)
Card Variable MID RO E PR ECC QH0 FT FC. Type A8 F F F F F F F. Default none none none 0.2 AUTO 0.3 none none
Note: This is an extended description of MAT_273 input provided by Peter Grassl It contains additional guidance on the choice of input parameters beyond the description in the official LS-DYNA manual Last
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 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 informationEVALUATION OF NONLOCAL APPROACHES FOR MODELLING FRACTURE IN NOTCHED CONCRETE SPECIMENS
VIII International Conference on Fracture Mechanics of Concrete and Concrete Structures FraMCoS-8 J.G.M. Van Mier, G. Ruiz, C. Andrade, R.C. Yu and X.X. Zhang (Eds) EVALUATION OF NONLOCAL APPROACHES FOR
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 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 information6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS
6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS Blondet et al. [25] carried out a cyclic test on an adobe wall to reproduce its seismic response and damage pattern under in-plane loads. The displacement
More informationLecture #8: Ductile Fracture (Theory & Experiments)
Lecture #8: Ductile Fracture (Theory & Experiments) by Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling of Materials in Manufacturing 2015 1 1 1 Ductile
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 informationarxiv: v1 [cond-mat.mtrl-sci] 24 Apr 2014
A simple and robust elastoplastic constitutive model for concrete arxiv:144.645v1 [cond-mat.mtrl-sci] 24 Apr 214 F. Poltronieri 1, A. Piccolroaz 1, D. Bigoni 1, S. Romero Baivier 2 1 Department of Civil,
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 informationDynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models
Dynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models I. Rhee, K.J. Willam, B.P. Shing, University of Colorado at Boulder ABSTRACT: This paper examines the global
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 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 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 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 informationAn anisotropic continuum damage model for concrete
An anisotropic continuum damage model for concrete Saba Tahaei Yaghoubi 1, Juha Hartikainen 1, Kari Kolari 2, Reijo Kouhia 3 1 Aalto University, Department of Civil and Structural Engineering 2 VTT 3 Tampere
More information5 ADVANCED FRACTURE MODELS
Essentially, all models are wrong, but some are useful George E.P. Box, (Box and Draper, 1987) 5 ADVANCED FRACTURE MODELS In the previous chapter it was shown that the MOR parameter cannot be relied upon
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 informationDevelopment and numerical implementation of an anisotropic continuum damage model for concrete
Development and numerical implementation of an anisotropic continuum damage model for concrete Juha Hartikainen 1, Kari Kolari 2, Reijo Kouhia 3 1 Tampere University of Technology, Department of Civil
More informationNumerical Simulation of Reinforced Concrete Beam with Utilization of Elasto-plastic Material Model of Concrete
Numerical Simulation of Reinforced Concrete Beam with Utilization of Elasto-plastic Material Model of Concrete FILIP HOKEŠ, MARTIN HUŠEK, PETR KRÁL, JIŘÍ KALA Faculty of Civil Engineering, Institute of
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 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 informationModelling 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 informationAnálisis Computacional del Comportamiento de Falla de Hormigón Reforzado con Fibras Metálicas
San Miguel de Tucuman, Argentina September 14 th, 2011 Seminary on Análisis Computacional del Comportamiento de Falla de Hormigón Reforzado con Fibras Metálicas Antonio Caggiano 1, Guillermo Etse 2, Enzo
More informationChapter 7. Highlights:
Chapter 7 Highlights: 1. Understand the basic concepts of engineering stress and strain, yield strength, tensile strength, Young's(elastic) modulus, ductility, toughness, resilience, true stress and true
More informationDiscrete Element Modelling of a Reinforced Concrete Structure
Discrete Element Modelling of a Reinforced Concrete Structure S. Hentz, L. Daudeville, F.-V. Donzé Laboratoire Sols, Solides, Structures, Domaine Universitaire, BP 38041 Grenoble Cedex 9 France sebastian.hentz@inpg.fr
More informationModeling of Interfacial Debonding Induced by IC Crack for Concrete Beam-bonded with CFRP
Proceedings of the World Congress on Engineering 21 Vol II WCE 21, June 2 - July 1, 21, London, U.K. Modeling of Interfacial Debonding Induced by IC Crack for Concrete Beam-bonded with CFRP Lihua Huang,
More informationCONCRETE FATIGUE MODEL BASED ON CUMULATIVE INELASTIC SHEAR STRAINS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 11-1 June 218, Glasgow, UK COCRETE FATIGUE MODEL BASED O CUMULATIVE IELASTIC
More informationLecture #7: Basic Notions of Fracture Mechanics Ductile Fracture
Lecture #7: Basic Notions of Fracture Mechanics Ductile Fracture by Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling of Materials in Manufacturing
More informationNew methodology for calculating damage variables evolution in Plastic Damage Model for RC structures
New methodology for calculating damage variables evolution in Plastic Damage Model for RC structures B. Alfarah 1, F. López-Almansa 2, S. Oller 3 1 PhD Candidate, Technical University of Catalonia, Civil
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 informationINCREASING RUPTURE PREDICTABILITY FOR ALUMINUM
1 INCREASING RUPTURE PREDICTABILITY FOR ALUMINUM Influence of anisotropy Daniel Riemensperger, Adam Opel AG Paul Du Bois, PDB 2 www.opel.com CONTENT Introduction/motivation Isotropic & anisotropic material
More informationAN UNLOADING AND RELOADING STRESS-STRAIN MODEL FOR CONCRETE CONFINED BY TIE REINFORCEMENTS
AN UNLOADING AND RELOADING STRESS-STRAIN MODEL FOR CONCRETE CONFINED BY TIE REINFORCEMENTS Junichi SAKAI 1 And Kazuhiko KAWASHIMA SUMMARY This paper presents a series of uniaxial compressive loading tests
More informationA rate-dependent Hosford-Coulomb model for predicting ductile fracture at high strain rates
EPJ Web of Conferences 94, 01080 (2015) DOI: 10.1051/epjconf/20159401080 c Owned by the authors, published by EDP Sciences, 2015 A rate-dependent Hosford-Coulomb model for predicting ductile fracture at
More 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 informationTHE BEHAVIOUR OF REINFORCED CONCRETE AS DEPICTED IN FINITE ELEMENT ANALYSIS.
THE BEHAVIOUR OF REINFORCED CONCRETE AS DEPICTED IN FINITE ELEMENT ANALYSIS. THE CASE OF A TERRACE UNIT. John N Karadelis 1. INTRODUCTION. Aim to replicate the behaviour of reinforced concrete in a multi-scale
More informationPOST-PEAK BEHAVIOR OF FRP-JACKETED REINFORCED CONCRETE COLUMNS
POST-PEAK BEHAVIOR OF FRP-JACKETED REINFORCED CONCRETE COLUMNS - Technical Paper - Tidarut JIRAWATTANASOMKUL *1, Dawei ZHANG *2 and Tamon UEDA *3 ABSTRACT The objective of this study is to propose a new
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 informationMODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam
MODELING OF CONCRETE MATERIALS AND STRUCTURES Class Meeting #1: Fundamentals Kaspar Willam University of Colorado at Boulder Notation: Direct and indicial tensor formulations Fundamentals: Stress and Strain
More informationModelling the nonlinear shear stress-strain response of glass fibrereinforced composites. Part II: Model development and finite element simulations
Modelling the nonlinear shear stress-strain response of glass fibrereinforced composites. Part II: Model development and finite element simulations W. Van Paepegem *, I. De Baere and J. Degrieck Ghent
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 informationMICROPLANE MODEL M4 FOR CONCRETE. II: ALGORITHM AND CALIBRATION
MICROPLANE MODEL M4 FOR CONCRETE. II: ALGORITHM AND CALIBRATION By Ferhun C. Caner 1 and Zdeněk P. Bažant, 2 Fellow, ASCE ABSTRACT: This paper represents Part II of a two-part study in which a new improved
More informationNUMERICAL SIMULATIONS OF CORNERS IN RC FRAMES USING STRUT-AND-TIE METHOD AND CDP MODEL
Numerical simulations of corners in RC frames using Strut-and-Tie Method and CDP model XIII International Conference on Computational Plasticity. Fundamentals and Applications COMPLAS XIII E. Oñate, D.R.J.
More informationMesoscopic Simulation of Failure of Mortar and Concrete by 3D RBSM
Journal of Advanced Concrete Technology Vol., No., 85-4, October 5 / Copyright 5 Japan Concrete Institute 85 Scientific paper Mesoscopic Simulation of Failure of Mortar and Concrete by D RBSM Kohei Nagai,
More informationUniversity of Sheffield The development of finite elements for 3D structural analysis in fire
The development of finite elements for 3D structural analysis in fire Chaoming Yu, I. W. Burgess, Z. Huang, R. J. Plank Department of Civil and Structural Engineering StiFF 05/09/2006 3D composite structures
More informationAvailable online at ScienceDirect. Procedia Engineering 161 (2016 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 161 (2016 ) 892 897 World Multidisciplinary Civil Engineering-Architecture-Urban Planning Symposium 2016, WMCAUS 2016 Parameter
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 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 informationEquivalent localization element for crack band approach to mesh-sensitivity in microplane model
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2005; 62:700 726 Published online 3 December 2004 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/nme.1216
More informationModule 5: Failure Criteria of Rock and Rock masses. Contents Hydrostatic compression Deviatoric compression
FAILURE CRITERIA OF ROCK AND ROCK MASSES Contents 5.1 Failure in rocks 5.1.1 Hydrostatic compression 5.1.2 Deviatoric compression 5.1.3 Effect of confining pressure 5.2 Failure modes in rocks 5.3 Complete
More 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 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 informationEFFECTS OF CONFINED CONCRETE MODELS ON SIMULATING RC COLUMNS UNDER LOW-CYCLIC LOADING
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 1498 EFFECTS OF CONFINED CONCRETE MODELS ON SIMULATING RC COLUMNS UNDER LOW-CYCLIC LOADING Zongming HUANG
More informationSIZE EFFECTS IN THE BIAXIAL TENSILE-COMPRESSIVE BEHAVIOUR OF CONCRETE: PHYSICAL MECHANISMS AND MODELLING
Fracture Mechanics of Concrete Structures, Proceedings FRAMCOS-2, edited by Folker H. Wittmann, AEDIFICA TIO Publishers, D-79104 Frei burg (1995) SIZE EFFECTS IN THE BIAXIAL TENSILE-COMPRESSIVE BEHAVIOUR
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 informationLattice Discrete Particle Model (LDPM) for Failure Behavior of Concrete. II: Calibration and Validation.
Lattice Discrete Particle Model (LDPM) for Failure Behavior of Concrete. II: Calibration and Validation. By Gianluca Cusatis 1, Andrea Mencarelli 2, Daniele Pelessone 3, James Baylot 4 A Paper Submitted
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 informationStrain-Based Design Model for FRP-Confined Concrete Columns
SP-230 57 Strain-Based Design Model for FRP-Confined Concrete Columns by N. Saenz and C.P. Pantelides Synopsis: A constitutive strain-based confinement model is developed herein for circular concrete columns
More informationMODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Isotropic Elastic Models: Invariant vs Principal Formulations
MODELING OF CONCRETE MATERIALS AND STRUCTURES Kaspar Willam University of Colorado at Boulder Class Meeting #2: Nonlinear Elastic Models Isotropic Elastic Models: Invariant vs Principal Formulations Elastic
More informationNumerical investigation of EDZ development around a deep polymetallic ore mine
Paper No. 198 ISMS 2016 Numerical investigation of EDZ development around a deep polymetallic ore mine Mountaka Souley a *, Marwan Al Heib a, Vincent Renaud a a INERIS, c/o Ecole des Mines de Nancy, Campus
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 informationEnhancing Prediction Accuracy In Sift Theory
18 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS Enhancing Prediction Accuracy In Sift Theory J. Wang 1 *, W. K. Chiu 1 Defence Science and Technology Organisation, Fishermans Bend, Australia, Department
More informationSIMPLIFIED CONCRETE MODELING WITH *MAT_CONCRET_DAMAGE_REL3
SIMPLIFIED CONCRETE MODELING WITH *MAT_CONCRET_DAMAGE_REL3 Leonard E Schwer Schwer Engineering & Consulting Services, Windsor CA, USA and L. Javier Malvar Karagozian & Case Structural Engineers, Burbank
More informationNumerical Characterization of Concrete Heterogeneity
Vol. Materials 5, No. Research, 3, 2002Vol. 5, No. 3, Statistical 309-314, 2002. Characterization of the Concrete Numerical Modeling of Size Effect In Heterogeneity 2002 309 Numerical Characterization
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 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 informationConstitutive Modeling of Reinforced Concrete Panel Behavior under Cyclic Loading
Constitutive Modeling of Reinforced Concrete Panel Behavior under Cyclic Loading K. Orakcal & D. Ulugtekin Bogazici University, Istanbul, Turkey L. M. Massone University of Chile, Santiago, Chile SUMMARY:
More informationAnisotropic Damage Mechanics Modeling of Concrete under Biaxial Fatigue Loading
Open Journal of Civil Engineering, 2015, 5, 8-16 Published Online March 2015 in SciRes. http://www.scirp.org/journal/ojce http://dx.doi.org/10.4236/ojce.2015.51002 Anisotropic Damage Mechanics Modeling
More informationHeterogeneous structures studied by interphase elasto-damaging model.
Heterogeneous structures studied by interphase elasto-damaging model. Giuseppe Fileccia Scimemi 1, Giuseppe Giambanco 1, Antonino Spada 1 1 Department of Civil, Environmental and Aerospace Engineering,
More informationImpact and Crash Modeling of Composite Structures: A Challenge for Damage Mechanics
Impact and Crash Modeling of Composite Structures: A Challenge for Damage Mechanics Dr. A. Johnson DLR Dr. A. K. Pickett ESI GmbH EURO-PAM 99 Impact and Crash Modelling of Composite Structures: A Challenge
More informationMesoscopic Simulation of Failure of Mortar and Concrete by 2D RBSM
Journal of Advanced Concrete Technology Vol., No., 59-7, October / Copyright Japan Concrete Institute 59 Mesoscopic Simulation of Failure of Mortar and Concrete by D RBSM Kohei Nagai, Yasuhiko Sato and
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 informationDAMAGE MODEL FOR CONCRETE INCLUDING RESIDUAL HYSTERETIC LOOPS: APPLICATION TO SEISMIC AND DYNAMIC LOADING
Fracture Mechanics of Concrete Structures Proceedings FRAMCOS-3 AEDIFICA TIO Publishers, D-79104 Frei burg, Germany DAMAGE MODEL FOR CONCRETE INCLUDING RESIDUAL HYSTERETIC LOOPS: APPLICATION TO SEISMIC
More informationarxiv: v1 [cond-mat.mtrl-sci] 31 Oct 2008
Influence of aggregate size and volume fraction on shrinkage induced micro-cracking of concrete and mortar arxiv:0811.0019v1 [cond-mat.mtrl-sci] 31 Oct 2008 Abstract Peter Grassl* 1, Hong S. Wong 2 and
More informationMesh-sensitivity analysis of seismic damage index for reinforced concrete columns
Mesh-sensitivity analysis of seismic damage index for reinforced concrete columns Jun Won Kang 1a and Jeeho Lee 2 1 Department of Civil Engineering, Hongik University, 94 Wausan-ro, Mapo-gu, Seoul 04066,
More informationFailure behavior modeling of slender reinforced concrete columns subjected to eccentric load
520 Failure behavior modeling of slender reinforced concrete columns subjected to eccentric load Abstract This work presents a numerical model to simulate the failure behavior of slender reinforced concrete
More informationAnalysis of RC concrete structures subject to elevated temperatures in ZSoil v2018
Analysis of RC concrete structures subject to elevated temperatures in ZSoil v2018 Andrzej Truty and Krzysztof Podleś Cracow University of Technology & ZACE Services Ltd August, 2018 1 / 39 General goal
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 informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationMicroplane Model formulation ANSYS, Inc. All rights reserved. 1 ANSYS, Inc. Proprietary
Microplane Model formulation 2010 ANSYS, Inc. All rights reserved. 1 ANSYS, Inc. Proprietary Table of Content Engineering relevance Theory Material model input in ANSYS Difference with current concrete
More informationEQUIVALENT FRACTURE ENERGY CONCEPT FOR DYNAMIC RESPONSE ANALYSIS OF PROTOTYPE RC GIRDERS
EQUIVALENT FRACTURE ENERGY CONCEPT FOR DYNAMIC RESPONSE ANALYSIS OF PROTOTYPE RC GIRDERS Abdul Qadir Bhatti 1, Norimitsu Kishi 2 and Khaliq U Rehman Shad 3 1 Assistant Professor, Dept. of Structural Engineering,
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 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 informationMATHEMATICAL AND NUMERICAL MODEL OF ROCK/CONCRETE MECHANICAL BEHAVIOR IN A MULTI-PLANE FRAMEWORK
ROMAI J., 4, 2(2008), 146 168 MATHEMATICAL AND NUMERICAL MODEL OF ROCK/CONCRETE MECHANICAL BEHAVIOR IN A MULTI-PLANE FRAMEWORK Seyed Amirodin Sadrnejad Faculty of Civil Engineering, K.N.Toosi University
More informationPlastic Damage Analysis of Koyna Dam in Different Damping Mechanisms with Dam Water Interaction
Plastic Damage Analysis of Koyna Dam in Different Damping Mechanisms with Dam Water Interaction O. Omidi & V. Lotfi Amirkabir University of Technology, Iran S. Valliappan University of New South Wales,
More informationPrediction of torsion shear tests based on results from triaxial compression tests
Prediction of torsion shear tests based on results from triaxial compression tests P.L. Smith 1 and N. Jones *2 1 Catholic University of America, Washington, USA 2 Geo, Lyngby, Denmark * Corresponding
More informationNUMERICAL 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 informationFRACTURE IN HIGH PERFORMANCE FIBRE REINFORCED CONCRETE PAVEMENT MATERIALS
FRACTURE IN HIGH PERFORMANCE FIBRE REINFORCED CONCRETE PAVEMENT MATERIALS ERIK DENNEMAN A thesis submitted in partial fulfilment of the requirements for the degree of PHILOSOPHIAE DOCTOR (ENGINEERING)
More informationFinite Element Modelling with Plastic Hinges
01/02/2016 Marco Donà Finite Element Modelling with Plastic Hinges 1 Plastic hinge approach A plastic hinge represents a concentrated post-yield behaviour in one or more degrees of freedom. Hinges only
More informationOptimization as a tool for the inverse identification of parameters of nonlinear material models
Optimization as a tool for the inverse identification of parameters of nonlinear material models FILIP HOKES, JIRI KALA, ONDREJ KRNAVEK Institute of Structural Mechanics Brno University of Technology,
More informationCode_Aster. SSNP161 Biaxial tests of Kupfer
Titre : SSNP161 Essais biaxiaux de Kupfer Date : 10/10/2012 Page : 1/8 SSNP161 Biaxial tests of Kupfer Summary: Kupfer [1] is interested to characterize the performances of the concrete under biaxial loadings.
More informationModule-4. Mechanical Properties of Metals
Module-4 Mechanical Properties of Metals Contents ) Elastic deformation and Plastic deformation ) Interpretation of tensile stress-strain curves 3) Yielding under multi-axial stress, Yield criteria, Macroscopic
More informationLimit analysis of brick masonry shear walls with openings under later loads by rigid block modeling
Limit analysis of brick masonry shear walls with openings under later loads by rigid block modeling F. Portioli, L. Cascini, R. Landolfo University of Naples Federico II, Italy P. Foraboschi IUAV University,
More informationsottotitolo Department of Civil and Environmental Engineering
A thirty-year Titolo survey presentazione of Microplane Models sottotitolo Milano, Dr. Giovanni XX mese Di Luzio 20XX Department of Civil and Environmental Engineering Origin of microplane models Slip
More informationANSYS Mechanical Basic Structural Nonlinearities
Lecture 4 Rate Independent Plasticity ANSYS Mechanical Basic Structural Nonlinearities 1 Chapter Overview The following will be covered in this Chapter: A. Background Elasticity/Plasticity B. Yield Criteria
More informationReciprocal of the initial shear stiffness of the interface K si under initial loading; reciprocal of the initial tangent modulus E i of the soil
Appendix F Notation a b B C c C k C N C s C u C wt C θ D r D 1 D 2 D 10 D 30 Reciprocal of the initial shear stiffness of the interface K si under initial loading; reciprocal of the initial tangent modulus
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 informationA generalisation of the Hillerborg s model for the analytical evaluation of ductility of RC beams in bending
Magazine of Concrete Research, 21, 62, No. 8, August, 557 567 doi: 1.168/macr.21.62.8.557 A generalisation of the Hillerborg s model for the analytical evaluation of ductility of RC beams in bending E.
More informationFLEXURAL MODELLING OF STRAIN SOFTENING AND STRAIN HARDENING FIBER REINFORCED CONCRETE
Proceedings, Pro. 53, S.A.R.L., Cachan, France, pp.55-6, 7. FLEXURAL MODELLING OF STRAIN SOFTENING AND STRAIN HARDENING FIBER REINFORCED CONCRETE Chote Soranakom and Barzin Mobasher Department of Civil
More information