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 Engineering 2 VTT 3 Tampere University of Technology, Department of Mechanical Engineering and Industrial Systems 15th International Conference on Fracture and Damage Mechanics, 14-16 September 2016, Alicante, Spain
Outline 2 Ottosen s 4 parameter model formulation 5 Some results and future work Anisotropic damage R. Kouhia 15.9.2016 2/18
2 Ottosen s 4 parameter model formulation 5 Some results and future work Anisotropic damage R. Kouhia 15.9.2016 3/18
Introduction The non-linear behaviour of quasi-brittle materials under loading is mainly due to damage and micro-cracking rather than plastic deformation. Damage of such materials can be modelled using scalar, vector or higher order damage tensors. Failure of rock-like materials in tension is mainly due to the growth of the most critical micro-crack Failure of rock-like materials in compression can be seen as a cooperative action of a distributed microcrack array http://mps-il.com Anisotropic damage R. Kouhia 15.9.2016 4/18
2 Ottosen s 4 parameter model formulation 5 Some results and future work Anisotropic damage R. Kouhia 15.9.2016 5/18
Ottosen s 4 parameter model A J 2 + Λ J 2 + BI 1 σ c = 0, σ c { k1 cos[ 1 Λ = 3 arccos(k 2 cos 3θ)] if cos 3θ 0 k 1 cos[ 1 3 π 1 3 arccos( k 2 cos 3θ)] if cos 3θ 0. cos 3θ = 3 3 2 J 3 J 3/2 2, : Lode angle σ c : the uniaxial compressive strength I 1 = trσ: the first invariant of the stress tensor J 2 = 1 2 s : s, J 3 = det s = 1 3 trs3 : deviatoric invariants A, B, k 1, k 2 : material constants Anisotropic damage R. Kouhia 15.9.2016 6/18
(a) (b) Meridian plane & plane stress Figure 1: Comparison of the failure surfaces on the deviatoric plane. Green line indicates the Mohr-Coulomb criteria with tension cut-off, blue line the Ottosen model and red line the Barcelona model. (a) π-plane, (b) σm = fc. θ = 0 σe/fc θ = 60 7 5 4 3 2 1 0 1 (a) 6 5 4 3 2 1 σm/fc 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.2 Green line = Mohr-Coulomb with tension cut-off Blue line = Ottosen s model Red line = Barcelona model6 (b) 0.4 0.6 0.8 1.0 1.2 1.4 Figure 2: Comparison of the failure surfaces in the (a) meridian plane and in (b)plane stress. Line colours as in Fig.1. Black dots indicate the experimental results (Kupfer et al., 1969). σ2/fc σ1/fc Anisotropic damage R. Kouhia 15.9.2016 7/18
Deviatoric plane σ 1 σ 1 σ 2 σ 3 σ 2 σ 3 π plane σ m = f c (a) (b) Green line = Mohr-Coulomb with tension cut-off Blue line = Ottosen s model Red line = Barcelona model Kuva 5.11: Murtopintojen vertailu deviatorisella tasolla. Vihreä viiva on vedossa katkaistu Mohr-Coulomb, sininen Ottosenin malli ja punainen viiva kuvaa Lublinerin mallia. (a) π- taso, (b) σ m = f c. σ e /f c θ = 60 7 Anisotropic damage R. Kouhia σ 2 /f c 15.9.2016 8/18 σ 1 /f c
2 Ottosen s 4 parameter model formulation 5 Some results and future work Anisotropic damage R. Kouhia 15.9.2016 9/18
Thermodynamic formulation Two potential functions ψ c = ψ c (S), S = (σ, D, κ) Specific Gibbs free energy γ = ρ 0 ψ c σ : ɛ, γ 0 Clausius-Duhem inequality ϕ(w ; S), W = (Y, K) Dissipation potential γ B Y : Y + B K K ψ c Define Y = ρ 0 D K = ρ ψ c 0 κ ψ (ρ c ) ) 0 σ ɛ : σ + (Ḋ BY : Y + ( κ B K ) K = 0 ɛ = ρ 0 ψ c σ, Ḋ = B Y, κ = B K Anisotropic damage R. Kouhia 15.9.2016 10/18
2 Ottosen s 4 parameter model formulation 5 Some results and future work Anisotropic damage R. Kouhia 15.9.2016 11/18
Specific model Specific Gibbs free energy ρ 0 ψ c (σ, D, κ) = 1 + ν 2E Elastic domain [ trσ 2 + tr(σ 2 D) ] ν 2E (1 + 1 3 trd)(trσ)2 + ψ c,κ (κ) Σ = {(Y, K) f(y, K; σ) 0} where the damage surface is defined as f(y, K; σ) = A J 2 σ c0 + Λ J 2 + BI 1 (σ c0 + K) = 0 Anisotropic damage R. Kouhia 15.9.2016 12/18
Invariants in terms of Y J 2 = 1 1 + ν J 3 = 2 3(1 + ν) [ EtrY 1 (1 2ν)(trσ)2] 6 { E [tr(σy ) trσtry ] + 1 (1 2ν)(trσ)3} 9 ϕ(y, K; σ) = I Σ(Y, K; σ) where I Σ is the indicator function { 0 if (Y, K) Σ I Σ(Y, K; σ) = + if (Y, K) / Σ (0, ( 0), if f(y, K; σ) < 0 (B Y, B K) = λ f ) f, λ, λ 0, if f(y, K; σ) = 0 Y K Ḋ = λ f Y, f κ = λ K Anisotropic damage R. Kouhia 15.9.2016 13/18
2 Ottosen s 4 parameter model formulation 5 Some results and future work Anisotropic damage R. Kouhia 15.9.2016 14/18
Some results Uniaxial compression - ultimate compressive stength σ c = 32.8 MPa σ c0 = 18 MPa, σ t0 = 1 MPa, (I 1, J 2 ) = ( 5 3σ c0, 4σ c0 / 2) A = 2.694, B = 5.597, k 1 = 19.083, k 2 = 0.998 K = [a 1 (κ/κ max) + a 2 (κ/κ max) 2 ]/[1 + b(κ/κ max) 2 ] a 1 = 85.3 MPa, a 2 = 12.65 MPa, b = 0.7032 (a) σ11/σc 1.25 1 0.75 0.5 0.25 model exp. Damage 0.4 0.3 0.2 0.1 (b) D 11 D 22 = D 33 (a) Damage 0.04 0.02 D D 22 = D 0 0 0.5 1 1.5 2 0 0 0.5 1 0 0 0.02 ε11/εc ε 11/ε c Experimental results from Kupfer et al. 1969. Figure 6: (a) Predicted damage-strain curves of the constitutive model for the concrete men under uniaxial compression, (b) Predicted damage-strain curves of the constitutive Figure 2. (a) Stress-strain diagram in uniaxial compression. for the concreteexperimental specimen under data uniaxial fromtension. Ref. [9], (b) Damage evolution in uniaxial compression. Notice Anisotropic that damage damage is larger R. Kouhia in the planes parallel 15.9.2016 to the 15/18 loading direction indicating splitting failure mode.
E (GPa) 30 20 10 Young s modulus (a) and apparent Poisson s ratio (b) 0 0 0.5 1 1.5 2 2.5 σ11/σc 1.2 1 0.8 0.6 0.4 0.2 Const. Exp. 0 0 0.1 0.2 0.3 0.4 0.5 ε 11/ε c ν app Biaxial compression Figure 5: (a) The predicted reduce of the Young s modulus E,(b) Apparent Poisson s ratio νapp = ε22/ε11 of the concrete (a) (a) (b) specimen under uniaxial compression plotted using the (b) constitutive model compared to experimental results (Chen, 1982). 0.012 0.012 1.2 D D 11 = 11 D = D 22 22 D 22 = D 33. This is a characteristic property of elastic-brittle materials (Basista, D D 33 33 2003; Murakami and Kamiya, 1997; Nemat-Nasser and Hori, 0.008 0.008 1993). In0.8 Fig.7 the volumetric strain-stress curves obtained by the solution of the constitutive equations and implementation of the hardening variables of equations (43) and (44) are illustrated and compared to the experimental observations (Kupfer 0.4 et al., 1969). Eq. model Although (43) some deviation from0.004 experimental 0.004 data is 0.4 observed, the model with Eq. hardening exper. (44) variable (43) qualitatively predicts the volumetric change of the material. Exp. Applying the hardening variable (44), however, causes the 0 model to fail to predict the volume increase which will eventually result in positive 0 0 values for V/V, after the material reaches its minimum volume. 0 0.4 0.4 0.8 0.8 1.2 1.2 1.6 1.6 0 0 0.2 0.2 0.4 0.4 0.6 0.6 The predicted behaviour ε 11/ε of cc the concrete specimen subjected to equibiaxial ε 11/ε ε 11/ε c c compression is compared to the experimental data (Kupfer et al., 1969), as can Figure be9: seen 8: (a) (a) incomparison Fig.8a. between Although between the predicted the predicted model stress-strain with stress-strain both curves hardening curves for concrete forvariables concrete specimen specimen gives under under equibiaxial an acceptable equibiaxial compression compression approximation with withof experimental hardening the ultimate variables data equibiaxial (Kupfer of equations etcompressive al., 1969), (43) and (b) stress, Predicted (44) and theex- perimental curves of the constitutive model for the concrete specimen under equibiaxial damage-strain magnitude data of (Kupfer the predicted et al., 1969), strain (b) ε Predicted 11 corresponding damage-strain Anisotropic damage to the curves ultimate of R. thekouhia stress constitutive 15.9.2016 16/18 is compression. model for the concrete specimen under equibiaxial compression. approximately 60% smaller than expected. Moreover, as shown in Fig.8b, pre- σ11/σc σ11/σc Damage Damage
2 Ottosen s 4 parameter model formulation 5 Some results and future work Anisotropic damage R. Kouhia 15.9.2016 17/18
Conclusions and future work Continuum damage formulation of the Ottosen s 4 parameter model Can model axial splitting Implementation into FE software (own codes, ABAQUS) Development of directional hardening model Regularization by higher order gradients Thank you for your attention! Juana Francés (1926-1990) 1960 Anisotropic damage R. Kouhia 15.9.2016 18/18