A non-local model with tension and compression damage mechanisms
|
|
- Louise Rodgers
- 6 years ago
- Views:
Transcription
1 Eur. J. Mech. A/Solids ) Éditions scientifiques et médicales Elsevier SAS. All rights reserved S ) /FLA A non-local model with tension and compression damage mechanisms Claudia Comi 1 Department of Structural Engineering, Politecnico di Milano, P. Leonardo da inci, Milano, Italy Received 15 November 1999; revised and accepted 14 June 2000) Abstract A non-local isotropic damage model is proposed which can be used to predict the behaviour of rock-like materials up to failure. Two isotropic damage variables account for the progressive degradation of mechanical properties under stress states of prevailing tension and compression and two internal lengths, one for tension and the other for compression, are introduced as localisation limiters. A linear bifurcation analysis highlights the regularisation properties of the non-local model. An iterative scheme for the numerical solution of the finite-step problem consisting of a linear global predictor, an averaging phase and a non-linear local corrector is presented. Some illustrative examples of tension and splitting tests show the effectiveness of the proposed model Éditions scientifiques et médicales Elsevier SAS damage / non-local model / localisation 1. Introduction The correct modelling of mechanical behaviour of quasi-brittle materials such as concrete and rocks subject to multiaxial tension-compression stress states still represents a challenging theme, especially when the prediction of failure is in point. A key dissipative phenomenon for such materials is certainly microcracking which results in non-symmetric, progressive degradation of mechanical properties under tension and compression conditions. Strength and stiffness degradation can be effectively modelled in the framework of continuous damage mechanics: isotropic damage models with a single damage variable see, e.g., Mazars, 1984)) or with a tension damage variable and a compression damage variable see, e.g., Mazars, 1986; Mazars and Pijaudier-Cabot, 1989)), anisotropic damage models see, e.g., Ramtani et al., 1992; Halm and Dragon, 1996)), damage coupled to plasticity models see, e.g., Ortiz, 1985; Simo and Ju, 1987; Lee and Fenves, 1998)) have been proposed, to cite a few. However, as widely studied see, e.g., Maier and Hueckel, 1979; Benallal et al., 1988; Desoyer and Cormery, 1994; Comi et al., 1995)) the boundary value problem in the presence of damage-induced softening may become ill-posed and severe difficulties arise in numerical analyses. arious regularisation techniques have been proposed in the literature, leading to enhanced formulations; namely, micropolar Steinmann and Stein, 1994), non-local integral Pijaudier-Cabot and Bažant, 1987; Pijaudier-Cabot, 1995), rate dependent Dubé et al., 1994) and gradient-dependent damage formulations Frémond and Nedjar, 1996; Peerlings et al., 1996; Comi, 1999). In all these models the explicit or implicit introduction of a material characteristic length fixes the width of the zone in which damage localises, thus preventing strain localisation into a line with consequent zero energy dissipation. The experimental measure of this length is far from being trivial. Under uniaxial tension conditions concretelike materials exhibit a localised failure which is governed by the fracture energy. Comparison of experimental 1 comi@stru.polimi.it
2 2 C.Comi results concerning tension tests on concrete specimens in which localised failure occurs and tension tests in which a particular experimental apparatus forces the damage to diffuse in the whole specimen allows the definition of a material internal length related to the aggregates size Bažant and Pijaudier-Cabot, 1989). More recently, experimental studies of Jansen and Shah 1997) concerning simple compression of concrete showed that also the compression failure is localised and that a compression fracture energy and a compression internal length can be measured. This latter length turns out to be very different from the tension internal length. To simulate this behaviour which is non-symmetric both for the homogeneous and localised response, in this paper a two-surface non-local damage model is proposed. In this model the variation in load carrying capacity hardening followed by softening) and the elastic degradation are described by two scalar damage variables D t and D c, ranging from zero to one. The first one measures damage for stress states of prevailing tension and the second one for stress states of prevailing compression. The current elastic domain in stress space is limited by two surfaces F t and F c, evolving with damage and non-locality is introduced by means of the weighted averages of the strain invariants which affect the loading functions. The presence in this model of two loading surfaces, besides a correct definition of the elastic domain and of its evolution, gives the possibility of introducing two different internal material lengths, related to the width of the localised band for stress states of prevailing compression or tension. This model represents the non-local enhancement of the fracture-energy based, local model originally proposed in Comi and Perego, 2000). The paper is organised as follows. In section 2 the underlying local model is briefly presented and the non-local model is formulated. In section 3 a bifurcation analysis for the rate problem is performed and the regularising role of non-locality is highlighted. The numerical aspects related to the use of this model in finite element analyses are discussed in section 4. To solve the time-step problem, a predictor averaging corrector scheme is proposed in which the correction can be performed at the local Gauss-point level and only a procedure to evaluate the non-local value of the predicted strain invariants is required in addition to the customary procedure for non-linear local models. In section 5 numerical results obtained with the proposed non-local model are compared with experimental results concerning tension and splitting tests on concrete specimens. As the mesh is refined, convergence of numerical results is observed both in terms of the damage pattern and of the global response. 2. Two-surfaces isotropic damage model 2.1. Local model The local model considered and enhanced in what follows is the bi-dissipative isotropic damage model proposed in Comi and Perego, 2000). The model is characterised by two scalar damage variables D t and D c, ranging from zero to one, which measure damage under stress states of prevailing tension and of prevailing compression. These variables enter in the definition of the free energy ψ: ψ = 1 2 { 2µ0 1 D t )1 D c )e : e + K 0 1 D t )tr + ε) 2 + K 0 1 D c )tr ε) 2}, 1) where ε is the small-strain tensor, e is the strain deviator tensor, µ 0 and K 0 are the initial, undamaged, shear and bulk elastic moduli, respectively. The first addend represents the deviatoric strain energy which is affected by both damage variables. The volumetric free energy is split into two parts second and third addends), one depending on the positive part of the volumetric strain tr + ε tr + trε+ tr ε ε = ) and affected by the damage in 2
3 A non-local damage model 3 tension and the other depending on the negative part of the volumetric strain tr ε tr ε = trε trε ) and affected 2 by the damage in compression. The assumed expression of the free energy allows to partially model, while preserving isotropy, the unilateral effect due to fracture closure under compression experimentally observed see, e.g., Ramtani et al., 1992)). A discussion on the anisotropy induced by a different description of the unilateral effect can be found in Fichant et al., 1999). The stress strain relation is obtained from 1) as: σ = ψ ε = 2µ e + K + tr + ε I + K tr ε I 2) σ being the stress tensor, I the second-order identity tensor, µ µ 0 1 D t )1 D c ) the current damaged value of the shear modulus, K + K 0 1 D t ) the current damaged value of the bulk modulus when tr ε 0, and K K 0 1 D c ) the current damaged value of the bulk modulus when tr ε < 0. The activation and evolution of damage are governed by the following loading unloading conditions: with f t 0, Ḋ t 0, f t Ḋ t = 0, f c 0, Ḋ c 0, f c Ḋ c = 0 3) f t = J 2 a t I1 2 + b tr t D t )I 1 k t rt 2 D t)1 αd c ), 4) f c = J 2 + a c I1 2 + b cr c D c )I 1 k c rc 2 D c), 5) 1 1 σ e σ 0 )i D 2 D 0i D i ) 2 for D i <D 0i 0i r i D i ) = [ ) Di D ci ] 0.75 i = t or c. 6) 0i 1 for D i D 0i 1 D 0i The loading functions f t and f c are defined by 4) 5) in the stress space: J 2 = 1 s : s is the second invariant of 2 the stress deviator s, I 1 = tr σ is the first invariant of the stress tensor, a t,b t,k t,a c,b c,k c and α are non-negative material parameters, r t D t ) and r c D c ) are the hardening softening functions defined in equation 6). In these latter, D 0t and D 0c are the values of damage at the peak σ 0t and σ 0c of the stress strain curve in tension and compression, respectively, σ et and σ ec are the elastic limit values, c t > 1andc c > 1 are parameters defining the slope of the softening branch. The term 1 αd c ) in f t allows the modelling of the decrease of the load carrying capacity in tension due to compression damage. The tension loading function f t is a hyperboloid around the hydrostatic axis in the Haigh Westergaard space of principal stresses, with circular deviatoric sections and hyperbolic meridian sections. The compression loading function f c is an ellipsoid see figure 1). The elastic domain delimited by these functions initially expands due to damage hardening and then shrinks due to damage softening. The boundary of the failure domain, understood as the set of stress states corresponding to the peak of stress strain curves recorded in tests with radial stress paths, can be easily computed from equations 4), 5) by substituting, in f t and f c, r t D t ) = 1andr c D c ) = 1, i.e. the values of the hardening functions corresponding to the peak values of damage D 0t and D 0c. Under plane stress conditions the resulting rupture domain is represented in figure 2a) together with the experimental points concerning biaxial tests on concrete specimens taken from Kupfer et al., 1969). Figure 2b)
4 4 C.Comi Figure 1. Loading surfaces and resulting elastic domain. Figure 2. Biaxial tests. a) Failure domain: experimental points and numerical curves; b) Stress strain curves under biaxial compression; bold lines: numerical results, light lines: experimental results. shows a comparison between the experimental stress strain curves thin lines) and the curves obtained with the present model thick lines) for compression tests with different confinements. Using equation 2), the loading functions can be expressed equivalently in the strain space as: f t = f t ε,d t,d c ) = 4µ 2 J ε 9a t K + tr + ε + K tr ε) 2 + 3b t r t D t )K + tr + ε + K tr ε) k t r 2 t D t)1 αd c ), 7) f c = f c ε,d t,d c ) = 4µ 2 J ε + 9a c K + tr + ε + K tr ε) 2 + 3b c r c D c )K + tr + ε + K tr ε) k c rc 2 D c), 8)
5 A non-local damage model 5 where J ε = 1 e : e is the second invariant of the strain deviator e. These latter expressions of the loading functions 2 will be used in the rest of the paper. In this model the limit values of damage D t = 1andD c = 1 can only be asymptotically reached for infinite strains, however, the specific fracture energy, i.e. the area below the stress strain curve in one-dimensional conditions, is finite owing to the chosen exponents of the softening function. Remark 2.1: The values D 0t, σ 0t and D 0c, σ 0c can in principle be directly measured in homogeneous uniaxial tension and compression tests, respectively. In practice, unless special experimental apparatus are used Bažant and Pijaudier-Cabot, 1989), the specimen response is not homogeneous due to strain localization. However, the pre-peak experimental behaviour is not significantly affected by strain localisation and therefore the direct identification of D 0t,D 0c,σ 0t and σ 0c do make sense. Also the results of localisation analysis for the proposed model Comi and Rizzi, 2000) under uniaxial conditions confirm that the inception of non-homogeneous responses is around the peak or in the softening regime. Since non-locality affects the overall behaviour only after localisation the meaning of the above parameters also holds true for the non-local model of the next section Non-local model The proposed enhanced version of the above model is based on the concept of non-local damage Pijaudier- Cabot and Bažant, 1987). As pointed out in Bažant and Pijaudier-Cabot, 1988), non-locality in damage models can be introduced both through the definition of a weighted average of the damage variable itself or of a strain measure. This second possibility is here exploited in view of a simpler numerical treatment. The basic non-local variables at a point x are here assumed to be the average strain invariants, i.e. the weighted averages over a volume of the local strain invariants: J ε x) = Wx s)j ε s) ds, 9) tr + εx) = Wx s) tr + εs) ds, tr εx) = Wx s) tr εs) ds, 10) where Wx s) is the weighting function, adequately defined to normalise the averaging. In equations 9), 10) and in what follows the symbol denotes the weighted average value of the quantity. In this work Wx s) is assumed as the normalised Gauss function and the average is extended to the whole body so that coincides with the body volume: Wx s) = 1 ) W 0 x) exp x s 2 2l 2 ) x s 2 with W 0 x) = exp ds. 11) 2l 2 The length l is a material parameter which can be related to the width of the zone in which damage phenomena localise. No particular provisions need to be introduced for points near the boundary of the body since W 0 x) in equation 11) already normalises the averaging. For concrete-like materials the localised behaviour is very different in tension and compression see, e.g., Bažant and Pijaudier-Cabot, 1989)) for experimental results in tension and Jansen and Shah, 1997) for experimental results in compression) and it seems reasonable to introduce two different internal lengths l t and l c, one for prevailing tension states and the other for prevailing compression states, to take into account such non-symmetry. Therefore, two different weighting functions are used in 9) 10) to obtain the mean values of
6 6 C.Comi strain invariants to be introduced in the loading function in tension and compression, respectively: W t x s) = exp ) x s 2 2lt 2 exp ), W c x s) = exp ) x s 2 2lc 2 x s 2 2l ds t 2 exp ). 12) x s 2 2l ds c 2 The non-local model is then obtained by replacing in the loading functions the local strain invariants with their averages; in the loading function f t the averages should be performed with the weight W t, while in loading function f c the averages should be performed with the weight W c. Denoting by t and c the averages of the quantity performed with weights W t and W c, respectively, the resulting non-local loading functions F t and F c are: F t = 4µ 2 J ε t 9a t K+ tr + ε t + K tr ε t ) 2 + 3b t r t D t ) K + tr + ε t + K tr ε t ) kt r 2 t D t)1 αd c ), 13) F c = 4µ 2 J ε c + 9a c K+ tr + ε c + K tr ε c ) 2 + 3bc r c D c ) K + tr + ε c + K tr ε c ) kc r 2 c D c). 14) The non-local model is completed by the loading unloading conditions 3) in which f t and f c are replaced by F t and F c and by the local stress strain relation 2). Remarks: 2.2 In equations 13) 14) the damage variables are locally defined, but the damage at a point is influenced by the damage of the neighbouring points. In fact if damage tends to grow at a single point, say P, the strains increase at that point; as a result the average values of strain invariants which intervene in the loading functions of points in a proper surrounding of P increase, thus inducing damage also in the neighbouring points of P. Even though strains of all points of the body intervene in the average values defined in equations 9) 10) the shape of the weight functions allows to practically drop contributions outside a proper region, related to the characteristic lengths l t or l c. 2.3 Only the inelastic behaviour, involving damage activation, is treated as non-local, while the elastic relations remain local, thus satisfying the requirement already pointed out in Pijaudier-Cabot and Bažant, 1987). According to this basic assumption only the loading functions are made non-local: the activation of tensile damage is modified by the existence of a tensile internal length and the activation of compression damage is modified by the existence of a compression internal length. For this reason all averages are made using l t in F t and using l c in F c. The separation between states with positive volumetric strains and negative volumetric strains concerns the linear elastic behaviour during unloading; this behaviour is local and therefore the internal length does not need to change across the interface tr ε = With the assumed non-local definitions of F t and F c, equations 13) 14), the damage variables, locally defined, remain always limited between zero and one as happens for the local model. On the contrary, the introduction in the loading functions of a non-local definition of the damage variables would require ad hoc provisions to check that locally the damage does not exceed the limit value one. 2.5 A rigorous thermodynamic formulation of the proposed model can be obtained following the same lines proposed in Borino et al., 1999) for non-local plasticity models.
7 A non-local damage model 7 3. Bifurcation analysis 3.1. Rate problem Assuming that at each point of the solid full loading occurs with activation of both modes, i.e. considering the linear comparison solid Hill, 1959), the rate constitutive equations can be obtained from equation 2) and the consistency conditions. The stress strain relation in rate form reads: σ = Ẽ : ε AḊ t BḊ c 15) with [ Ẽ = 2µ I I + K + Htr ε) + K H tr ε) 2µ 3 ] I I, 16) A = 2µ 0 1 D c ) e + K 0 tr + ε I, 17) B = 2µ 0 1 D t ) e + K 0 tr ε I. 18) In the above relations Ẽ is the damaged isotropic elastic tensor, I I is the fourth order symmetric identity tensor such that its component ijhk is 1 2 δ ihδ jk + δ ik δ jh ), denotes the tensor product, H )is the Heavyside function Hx) = 1forx>0andHx) = 0forx<0). Relation 16) follows from 2) noting that: d d dt tr+ ε) = Htr ε) I : ε, dt tr ε) = H tr ε) I : ε, 19) ė = I I 13 ) I I : ε. 20) The damage increments are then evaluated from the consistency conditions. Consider first the local model. These conditions are: f t = f t D t Ḋ t + f t Ḋ c + f t : ε = 0, D c ε 21) f c = f c D t Ḋ t + f c Ḋ c + f c : ε = 0, D c ε 22) where the partial derivatives are computed from equations 7) and 8): f t = 8µµ 0 1 D c )J ε + 3K 0 tr + ε [ 6a t K+ tr + ε + K tr ε ) b t r t D t ) ] D t + 3b t dr t dd t K+ tr + ε + K tr ε ) 2k t r t D t ) dr t dd t 1 αd c ), 23) f t = 8µµ 0 1 D t )J ε + 3K 0 tr ε [ 6a t K+ tr + ε + K tr ε ) b t r t D t ) ] + αk t rt 2 D D t), c 24) f c = 8µµ 0 1 D c )J ε K 0 tr + ε [ 2a c K+ tr + ε + K tr ε ) + b c r c D c ), D t 25) f c D c = 8µµ 0 1 D t )J ε 3K 0 tr ε [ 6a c K+ tr + ε + K tr ε ) + b c r c D c ) ]
8 8 C.Comi + 3b c dr c dd c K+ tr + ε + K tr ε ) 2k c r c D c ) dr c dd c, 26) f t ε = 4µ2 e + 3 [ 6a t K+ tr + ε + K tr ε ) + b t r t D t ) ][ K + Htr ε) + K H tr ε) ] I, 27) f c ε = 4µ2 e + 3 [ 6a c K+ tr + ε + K tr ε ) + b c r c D c ) ][ K + Htr ε) + K H tr ε) ] I. 28) Solving the system 21) 22) and substituting the result into 15) one obtains the rate constitutive equations for the local model: σ = E t : ε; E t = Ẽ A f t f c D c ε f c f t ) D c ε B f c h f t D t ε f t f c D t ε with: h = f t f c f t f c. 30) D t D c D c D t Since two modes have been considered, the tangent tensor E t is a two rank one update of the damaged elastic tensor. If only one loading function is active the second or the third term in equation 29b) vanishes and the classical form of the tangent operator of single-surface models is recovered. Consider now the non-local model. The damage increments are again evaluated from the consistency conditions but using in this case the non-local loading functions 13) 14). These conditions are expressed by the following integral differential equations: h ) 29) where: Ḟ t = F t Ḋ t + F t Ḋ c + G t ε) = 0, 31) D t D c Ḟ c = F c Ḋ t + F c Ḋ c + G c ε) = 0, 32) D t D c F t D t = 8µµ 0 1 D c ) J ε t + 3K 0 tr + ε t [ 6at K+ tr + ε t + K tr ε t ) bt r t D t ) ] + 3b t dr t dd t K+ tr + ε t + K tr ε t ) 2kt r t D t ) dr t dd t 1 αd c ), 33) F t D c = 8µµ 0 1 D t ) J ε t + 3K 0 tr ε t [ 6at K+ tr + ε t + K tr ε t ) bt r t D t ) ] + αk t r 2 t D t), 34) F c D t = 8µµ 0 1 D c ) J ε c 3K 0 tr + ε c [ 6ac K+ tr + ε c + K tr ε c ) + bc r c D c ) ], 35) F c D c = 8µµ 0 1 D t ) J ε c 3K 0 tr ε c [ 6ac K+ tr + ε c + K tr ε c ) + bc r c D c ) ] + 3b c dr c dd c K+ tr + ε c + K tr ε c ) 2kc r c D c ) dr c dd c, 36) G t ε) = 4µ 2 W t e : ε d 3 [ 6a t K+ tr + ε t + K tr ε t ) bt r t D t ) ]
9 A non-local damage model 9 ] [K + W t Htr ε) I : ε d + K W t H tr ε) I : ε d, 37) G c ε) = 4µ 2 W c e : ε d + 3 [ 6a c K+ tr + ε c + K tr ) ε c + bc r c D c ) ] [K + W c Htr ε) I : ε d + K ] W c H tr ε) I : ε d. 38) Substituting the damage increments into equation 15) leads to the rate constitutive equation: with: σ = Ẽ : ε 1 H [ Ft A G c ε) F ) c Fc G t ε) + B G t ε) F )] t G c ε) D c D c D t D t H = F t F c F t F c. 40) D t D c D c D t Note that, due to the non-local nature of the model, a local tangent operator E t such that σ = E t : ε can no longer be defined. The rate problem describing the quasi-static evolution of an elasto-damageable solid is given by the rate constitutive equations equation 29) for the local model and equation 39) for the non-local one), the field equations of equilibrium and compatibility: div σ = 0, 41) 39) ε = 1 2 grad u + gradt u ), 42) where u is the displacement vector and small strains and zero rates of body forces have been assumed, and boundary conditions on displacements and traction rates Bifurcation conditions Starting from a homogeneous equilibrium state, possible bifurcated solutions of the rate problem with the non-local model of section 2.2 are sought in the form: where q is the wave number and n denotes the normal direction. The corresponding strain rate is: u = v 0 e iq n x, 43) ε = 1 2 iqv 0 n + n v 0 ) e iq n x. 44) This perturbation is admissible if it satisfies the rate equilibrium equations, with σ given by equation 39). As shown in the Appendix, this admissibility condition can be expressed in the form: where Lq) is a kind of non-local tangent operator defined as: Lq) = Ẽ 1 H n Lq) n) v 0 = 0, 45) { A W c q) f t f c D c ε W tq) f ) c f t + B W t q) f c f t D c ε D t ε W cq) f )} t f c. 46) D t ε
10 10 C. Comi The functions W t q) and W c q) for an infinite body are the Fourier transforms of the weight functions W t and W c, respectively: W t q) = W t ξ) e iq n ξ dξ = e l2 t q 2 2, 47) W c q) = W c ξ) e iq n ξ dξ = e l2 c q ) Thus, a perturbed solution of the form 43) exists if, and only if, the coefficient matrix of system 45) is singular, i.e.: detn Lq) n) = 0. 49) If one considers the local model, standard arguments lead to the following critical condition: detn E t n) = 0. 50) The tensor n Lq) n plays, for the non-local model here considered, a role similar to the one of the acoustic tensor n E t n for the local model, but the wave number q enters in the definition 46) through the Fourier transforms W t q) and W c q) of the weighting functions. Therefore unlike the local medium, in the non-local one the wavelength is not arbitrary at a bifurcation point. Localisation into a band of zero width, which occurs in the local medium as soon as condition 50) is fulfilled Rice, 1976; Borré and Maier, 1989), cannot occur in the non-local medium. In fact for zero wavelength q ) the Fourier transforms W t q) and W c q) of the weighting functions go to zero and Lq) tends to the elastic damaged tensor Ẽ which is positive definite for 0 D t < 1and0 D c < 1. Therefore the ill-posedness of the boundary value problem arising when using local damage models is avoided by the proposed non-local enhancement. The tensor n Lq) n reduces to the local acoustic tensor when W t q) = W c q) = 1, i.e. when the weights are Dirac distributions. To compare the bifurcation conditions of the local and non-local model assume that the internal length under compression l c be equal to the internal length under tension l t,thenw t q) = W c q) = Wq) for any q and one has: Lq) = Ẽ Wq) { ft f c A H D c ε f ) c f t fc f t + B D c ε D t ε f )} t f c. 51) D t ε Compare this tensor with the local tangent tensor of equation 29b): the two tensors are identical except for the presence of the term H/Wq) instead of h. Therefore, if perturbated solutions exist, for the local and the non-local models the direction n is the same. Denoting by H cr the maximum value of H, at varying n, which fulfils condition 49) and by h cr the maximum value of h, at varying n, which fulfils condition 50), being Wq) 1foranyq, one has: h cr H cr. 52) Note that, due to the presence of two loading surfaces, h and H are not directly related to a single scalar hardening softening parameter monotonically decreasing during loading. Therefore from equation 52) one cannot conclude that the criterion of bifurcation for the local continuum gives a lower bound for the corresponding non-local continuum in general. However, this is the case when bifurcation is sought starting from a state in which only one mode is active and the results obtained in Benallal and Pijaudier-Cabot, 1993) are recovered.
11 A non-local damage model Algorithmic aspects 4.1. Finite-step problem Consider the evolution of a damageable body of volume and boundary S subject to body forces Ft), surface forces ft) on the free part of the boundary S f, endowed with an outward normal m, and imposed displacements Ut) on the constrained part of the boundary S u and let the material behaviour be described by the non-local damage model of section 2.2. Numerical solution of this evolutive problem requires time and space discretizations. Time dicretization is performed by dividing the time interval of interest into finite steps. At the beginning, say t n, of a time step t all quantities are known and the solution must be computed at the end of the step, say t n+1, for given increments of external actions F, f and U. The finite step problem is then obtained by enforcing equilibrium and compatibility at the end of the step and integrating the constitutive law according to a backward difference scheme subscript n + 1 marks quantities computed at t n+1 and indicates the increment of a variable over the step): div σ n+1 + F n+1 = 0 in ; m σ n+1 = f n+1 on S f, 53) 1 grad un+1 + grad T ) u n+1 = εn+1 2 in ; u n+1 = U n+1 on S u, 54) σ n+1 = 2µ 0 1 D tn+1 )1 D cn+1 ) e n+1 + K 0 1 D tn+1 ) tr + ε n+1 I + K 0 1 D cn+1 ) tr ε n+1 I in,55) F t ) n+1 = F t ε n+1,d tn+1,d cn+1 ) 0; D t 0; F t ) n+1 D t = 0 in, 56) F c ) n+1 = F c ε n+1,d tn+1,d cn+1 ) 0; D c 0; F c ) n+1 D c = 0 in. 57) Space discretization is performed following a standard finite element displacement approach in which only displacements are modelled in terms of their nodal values ū ux) = N u x)ū) and strains are computed through compatibility. Note that even though the model is non-local, no additional variables should be modelled nor additional continuity requirements arise as it happens in gradient enhanced damage models see, e.g., Comi, 1999)) Iterative procedure The finite-step problem concerning the time interval t = t n+1 t n is non-linear and can be solved by an iterative scheme. The scheme here adopted is similar to the linear predictor non-linear corrector scheme used for local models except that an intermediate averaging phase should be added between predictor and corrector. In what follows the index n + 1 of the time step will be omitted and the superscripts i and i + 1 will refer to the iteration. Predictor phase At the iteration i + 1 for fixed damage variables values, compute the estimate of the nodal displacement vector ū i+1 from the linearized equations of the finite element model. Denoting by R i+1 the vector of residual forces and by S i+1 the prediction stiffness matrix, the equations to be solved read: S i+1 ū i+1 ū i) = R i+1. 58) The prediction matrix can be chosen in different ways. If in this phase the damage variables are set equal to the values reached at the end of the previous step, D t = D tn, D c = D cn, the procedure outlined above corresponds to a secant prediction matrix.
12 12 C. Comi Averaging phase At each Gauss point x g, compute the weighted mean values of the strain invariants by adding contributions of strains at all the NG Gauss points of the mesh: NG J ε x g ) = ω i Wx g x i )J ε x i ), 59) NG i=1 tr + εx g ) = ω i Wx g x i ) tr + εx i ), tr εx g ) = ω i Wx g x i ) tr εx i ). 60) i=1 In equations 59) 60) ω i is the Gauss weight of the ith Gauss point and W is the weighting function of the non-local model which is equal to W t equation 12a)) when computing average quantities to be introduced in the tension loading function and equal to W c equation 12b)) when computing average quantities to be introduced in the compression loading function. Note that the terms ω i W t x g x i ) and ω i W c x g x i ) can be computed and stored once for all at the beginning of the analysis. Corrector phase At the Gauss point x g, with the current of iteration i + 1), predicted values of the mean strain invariants and the values of damage of the previous time step, compute the trial loading functions F triali+1 t and F triali+1 c : F triali+1 t = F t ε i+1,d tn,d cn ), F trial i+1 c = F c ε i+1,d tn,d cn ) and calculate damage using the loading unloading conditions 56) 57). i) if Ft triali+1 0andFc triali+1 0, unloading occurs, the damage increments are zero, then Dt i+1 = D tn and Dc i+1 = D cn ; ii) if Ft triali+1 > 0 and/or Fc triali+1 > 0, compute the damage increments; the following three possibilities at most should be considered: tension activation: D t > 0and D c = 0, with D t computed from F t ε i+1,d tn + D t,d cn ) = 0; this solution is valid if F c ε i+1,d tn + D t,d cn ) 0; compression activation: D t = 0and D c > 0, with D c computed from F c ε i+1,d tn,d cn + D c ) = 0; this solution is valid if F t ε i+1,d tn,d cn + D c ) 0; vertex activation: D t > 0and D c > 0, with D t and D c computed from F t ε i+1,d tn + D t,d cn + D c ) = 0andF c ε i+1,d tn + D t,d cn + D c ) = 0. With the calculated values of damage variables, update the stresses. Convergence test Compute residual force vectors R i+1 with the current value of stress and check convergence comparing a norm of the residuals with a fixed tolerance. If the convergence test is not passed, perform a new iteration. The scheme of the iteration is summarized in box 1. It is worth noting that, despite the non-local nature of the material model, the corrector phase can be performed locally, at the Gauss point level. This is due to the particular choice to introduce non-local values of strain measures keeping damage variables locally defined. Thus in a strain driven process the non-local quantities are the input of the local corrector phase which furnishes updated estimates of the locally defined damage variables. NG i=1 61) 5. Numerical tests In order to investigate the predictive capability and the regularisation properties of the proposed model and the effectiveness of the numerical strategy, two simple preliminary numerical tests have been performed.
13 A non-local damage model 13 Box 1. Flow chart of the finite-step iterative scheme. In both examples, within a given time step, the secant prediction matrix computed at the end of the previous step has been used for all iterations. This choice results in a large number of iterations, but is quite efficient since it does not require the updating of the stiffness at each iteration. The proper definition of a consistent tangent matrix as proposed in Jirásek, 1999) is rather cumbersome for this model and further work is required in this direction Tension test on a four edge notched specimen The first example concerns the simulation of the experimental tension test performed by Hassanzadeh 1991) on four edge notched concrete specimen figure 3a)). This test has already been considered in the literature as a benchmark di Prisco et al., 2000). The material experimental compressive strength given by the author is 50 MPa, which corresponds to a concrete grade C40 in the CEB-FIP model code; therefore, the other data used to identify the model parameters have been taken from this code. Note that during the test only damage in tension is activated, therefore only the material parameters defining the loading function in tension are of interest. Both material and resulting model parameters are collected in figure 3. Even though the presence of the notch on the four edges of the specimen makes the response typically threedimensional, the finite element analyses have been performed under plane strain conditions; to account for the three dimensionality of the problem, the total reaction R has been obtained assuming axial-symmetric-like conditions as follows. From the computed nodal forces F i see figure 4a)), forces per unit area q i have been computed dividing F i by half of the support b i of node i see figure 4b)); these load densities have then been multiplied by suitable influence areas the square frames of width b i /2infigure 4c)) instead of the areas
14 14 C. Comi Figure 3. Tension test: a) geometrical data; b) 472-elements mesh; c) 908-elements mesh. Figure 4. Tension test: evaluation of the total reaction: a) nodal forces; b) forces per unit area; c) influence areas adopted; b) influence areas in plane strain conditions.
15 A non-local damage model 15 Figure 5. Tension test: experimental dotted line) and numerical global response. of 70 mm long strips see figure 4d)). The specimen has been discretized by two meshes of constant strain triangular elements figure 3b) c)). The computed total reactions R versus the imposed displacement of the top of the specimen u are plotted in figure 5 together with the experimental curve dotted line) taken from Hassanzadeh, 1991). One can observe that the two meshes give similar responses which are in quite good agreement with the experimental results. Note that the bump experimentally observed in the softening branch is due to rotational instability, occurring when the test apparatus has a low rotational stiffness Hillerborg, 1989). The two parts of the specimen, on both sides of the damaged zone rotate with respect to each other. This instability effect is not reproduced by the numerical analysis in which equal displacements are imposed to all points on the top face of the specimen. Also the final pattern of damage is almost mesh-independent, see figure 6, the non-local model allows for a spread of the damage over all the finite elements belonging to a zone whose dimensions are fixed by the a) b) Figure 6. Tension test: final damage distribution a) 472-elements mesh b) 908-elements mesh.
16 16 C. Comi a) b) c) d) Figure 7. Tension test: damage pattern evolution during the analysis 908-elements mesh). material length and not by the element size. The evolution of the damaged zone, simulating the formation of macro-cracks leading to the failure of the specimen, is shown in figure 7 for the refined mesh Splitting test on concrete prism The second example concerns the splitting test on prisms of plain micro-concrete, a test which is carried out under monotonically increasing imposed displacement. This example is chosen because both damage mechanisms in tension and compression are activated during this simple test; therefore its simulation requires the correct modelling of not only damage due to tension, but also damage due to compression. A similar test has been simulated with a two-surfaces plasticity model by Feenstra and de Borst 1996). The geometrical data, taken from Rocco et al., 1995), are shown in figure 8 together with the model parameters employed in the numerical analyses; these latter are identified on the basis of the material data given in that reference E = MPa, f c = 38 MPa, G f = N/mm) and of the data taken from the CEB-FIP
17 A non-local damage model 17 Figure 8. Splitting test: geometrical and material data. code. The analyses have been performed discretizing only a quarter of the specimen shown in figure 8)bytwo meshes of and triangular plane strain elements. The only experimental result given in Rocco et al., 1995) is the value of the reaction at the peak 26.7 kn) which is in good agreement with the computed value 27.2 kn). At the peak, both damage mechanisms are activated: figure 9 shows the corresponding contour plots of the damage variables D t and D c obtained with the coarse left hand side of the contour plots) and refined mesh right hand side of the contour plots). One can see that damage in tension develops at the centre of the specimen, where a crack forms, while damage in compression is activated below the bearing strips, giving rise to a punch. Again no spurious mesh dependence is visible. A map of the active modes during the analysis is displayed in figures 10a) e). Only the right, upper quarter of the specimen is represented. In black are represented regions where compression is active, in grey are represented regions where tension is active; light grey corresponds to elastic behaviour. Following the sequence, one can see that damage in tension begins at the middle of the specimen figure 10a)), then it spreads in both vertical and horizontal directions figures 10b) c)) and finally it tends to localize in a narrow vertical zone, thus simulating a macrocrack, and in the regions outside this zone elastic unloading occurs figures 10d) e)). Damage in compression develops below the bearing strips figures 10a) and b)) and it forms a punch figures 10c) e)). Figure 11 shows at each Gauss point, the active modes at the same step of the analysis of figure 10d) u = mm); the different post-treatment of the data allows the visualisation of the Gauss points where the vertex is active, marked by black circles. It can be observed, and this is true also for all the other steps, that the vertex is active in a very small number of points, differently to what would happen using multi-surfaces plasticity models. This is due to the fact that in the damage models there is no constraint of normality to be fulfilled, thus the line of vertex or the points of vertex in two dimensions) is not a privileged set of points where the constraint is relaxed, but it is only a single line on a surface of potentially active points.
18 18 C. Comi a) b) Figure 9. Splitting test: a) tension damage distribution at the peak, left elements mesh, right elements mesh; b) compression damage distribution at the peak, left elements mesh, right elements mesh. Figure 10. Splitting test: active modes at different steps of the analysis.
19 A non-local damage model 19 Figure 11. Splitting test: activation of vertex, compression or tension at the peak. 6. Conclusions A non-local isotropic damage model with two damage mechanisms has been developed. Two internal lengths have been introduced as localisation limiters of damage in tension and compression. For the proposed non-local model a bifurcation analysis has shown that the criterion for bifurcation is very similar to that of the underlying local continuum, but with the non-local model the boundary value problem remains well-posed at the onset of bifurcation. In fact while a perturbation with zero wavelength is an admissible solution to the rate problem in the local case, in the non-local case the wavelength is fixed and related to the material internal lengths. For the numerical solution of structural problems in the presence of this enhanced damage model, a simple iterative solution scheme of the finite time-step has been proposed. The effectiveness of the non-local enhancement, which allows the avoidance of spurious mesh dependence of numerical results, has been shown by some numerical tests. The correct identification of the internal, material lengths from experimental tests still remains an open question which, maybe, can be tackled by properly defined inverse problems. The simplicity of the underlying local model, in particular its isotropy, and the simplicity of the non-local enhancement, in which the damage variables rest locally defined, make the proposed model apt to be used in structural problems involving large computations as in dam engineering. The damage induced anisotropy linked to the formation of aligned micro-cracks is not introduced in the material model and is obtained only as a structural effect at the macro-scale. Even though a provision to describe the unilateral effect has been already included, the model is not able to correctly predict the behaviour in situations where loading reversal occurs, due to the lack of permanent strains in its present version. The coupling of the proposed damage model with plasticity and thus the inclusion of permanent strains is now in progress.
20 20 C. Comi Acknowledgements This work has been carried out in the framework of the project: Structural safety assessment of large dams supported by the Italian Ministry for Universities and Scientific Research. Appendix Consider the rate equilibrium equations for the non-local medium: { div Ẽ : ε 1 [ Ft A G c ε) F ) c Fc G t ε) + B G t ε) F )]} t G c ε) = 0. A.1) H D c D c D t D t Note first that, if bifurcation of the strain rates in the form 44) is sought starting from a homogeneous state, the averages of variables at this state do coincide with the corresponding local values, furthermore, only rate quantities give rise to non-zero terms in the divergence. Therefore equation A.1) can be rewritten as: div Ẽ : ε 1 [ Ft A div G c ε) F ) c Fc div G t ε) + B div G t ε) F )] t div G c ε) = 0. A.2) H D c D c D t D t To compute the term div G t ε), taking into account the definitions 37) and 44), one has to evaluate the following terms: term 1 = div [ 12 ] iqw t x s)e : v 0 n + n v 0 ) e iqn s ds, A.3) term 2 = div [ 12 ] iqw t x s)h±tr ε) I : v 0 n + n v 0 ) e iqn s ds. A.4) Analogous terms, with W c instead of W t, intervene in the evaluation of div G c ε). Perform the change of variables x s = ξ and note that, being W t symmetric, W t ξ) = W t ξ), from equations A.3) and A.4) one has: term 1 = iq e : v 0 n) W t ξ) e iqn ξ dξ div e iqn x = q 2 W t q)n e v 0 )n e iqn x, A.5) term 2 = iqh±tr ε) I : v 0 n) W t ξ) e iqn ξ dξ div e iqn x = q 2 W t q)h±tr ε)n I v 0 )n e iqn x, A.6) where W t q) is the Fourier transform of W t. Using equations A.5), A.6) and 37), remembering that the starting state is homogeneous and, hence, tr + ε t tr + ε c tr + ε, etc.), one can compute the divergence of G t ε) as: div G t ε) = q 2 W t q) { 4µ 2 n e v 0 ) [ 2a t K+ tr + ε + K tr ε ) b t r t D t ) ] [ K + Htr ε) + K H tr ε) ] n I v 0 ) } n e iqn x. A.7) Using the definition 27) of f t, one can rewrite A.7) as: ε div G t ε) = q 2 W t n f ) t ε v 0 n e iqn x. A.8)
21 A non-local damage model 21 Fully analogous passages allow to compute the divergence of G c ε) as: div G c ε) = q 2 W c n f ) c ε v 0 n e iqn x. A.9) Substituting A.8) and A.9) into A.2), condition 45) for the admissibility of the perturbation is arrived at. References Bažant Z.P., Pijaudier-Cabot G., Non-local continuum damage, localisation instability and convergence. J. Appl. Mech. 55, Bažant Z.P., Pijaudier-Cabot G., Measurement of characteristic length of non-local continuum. J. Eng. Mech. ASCE 115, Benallal A., Billardon R., Geymonat G., Some mathematical aspects of the damage softening problem. In: Mazars J., Bažant Z.P. Eds.), Cracking and Damage, Elsevier, Amsterdam, pp Benallal A., Pijaudier-Cabot G., Strain localisation and bifurcation in a non-local continuum. Int. J. Solids Structures 30 13), Borino G., Fuschi P., Polizzotto C., A thermodynamic approach to non-local plasticity and related variational principles. J. Appl. Mech. 66 4), Borré G., Maier G., On linear versus nonlinear flow rules in strain localisation analysis. Meccanica 24, CEB-FIP Model code 1990, Bullettin d information, CEB, Lausanne. Comi C., Berthaud Y., Billardon R., On localisation in ductile-brittle materials under compressive loadings. Eur. J. Mech., A/Solids 14 1), Comi C., Computational modelling of gradient-enhanced damage in quasi-brittle materials. Mech. Cohes.-Frict. Mater. 4 1), Comi C., Perego U., A bi-dissipative damage model for concrete with applications to dam engineering. In: Proc. ECCOMAS 2000, Barcelona, Spain. Comi C., Rizzi E., On bifurcation in local and non-local materials with tension and compression damage. In: Proc. ECCOMAS 2000, Barcelona, Spain. Desoyer T., Cormery F., On uniqueness and localisation in elastic-damage materials. Int. J. Solids Structures 31, di Prisco M., Ferrara L., Meftah F., Pamin J., de Borst R., Mazars J., Reynouard J.M., Mixed mode fracture in plain and reinforced concrete: some results on benchmark tests. Int. J. Fracture 103 2), Dubé J.-F., Pijaudier-Cabot G., La Borderie C., Glynn J., Rate dependent damage model for concrete- Wave propagation and localisation. In: Mang H., Bicanic N. and de Borst R. Eds.), Computer Modelling of Concrete Structures, pp Feenstra P.H., de Borst R., A composite plasticity model for concrete. Int. J. Solids Structures 33, Fichant S., La Borderie C., Pijaudier-Cabot G., Isotropic and anisotropic descriptions of damage in concrete structures. Mech. Cohes.-Frict. Mater. 4, Frémond M., Nedjar B., Damage, gradient of damage and principle of virtual power. Int. J. Solids Structures 33, Halm D., Dragon A., A model of anisotropic damage by mesocrack growth; unilateral effect. Int. J. Damage Mech. 5, Hassanzadeh M., Behaviour of fracture process zones in concrete influenced by simultaneously applied normal and shear displacements. PhD Thesis, Lund Institute of Technology, Sweden. Hill R., Some basic principles in the mechanics of solids without a natural time. J. Mech. Phys. Solids 7, Hillerborg A., Stability problems in fracture mechanics testing. In: Shah S.P., Swartz S.E., Barr B. Eds.), Fracture of concrete and rock, Elsevier Applied Science, pp Jansen D.C., Shah S.P., Effect of length on compressive strain softening of concrete. J. Eng. Mech. ASCE 123, Jirásek M., Computational aspects of non-local models. In: Proc. European Conf. Comp. Mech., München. Kupfer H., Hilsdorf H.K., Rush H., Behavior of concrete under biaxial stress. ACI J. 66, Lee J., Fenves G.L., Plastic-damage model for cyclic loading of concrete structures. J. Eng. Mech. ASCE 124, Maier G., Hueckel T., Nonassociated and coupled flow rules of elastoplasticity for rock-like materials. Int. J. Rock Mech. Min. Sci. & Geomech. Abstr. 16, Mazars J., Application de la mécanique de l endommagement au comportement non linéaire et à la rupture du béton de structure. Dr of Science thesis University of Paris 6. Mazars J., A model of a unilateral elastic damageable material and its application to concrete. In: Fracture Toughness and Fracture Energy of Concrete, Elsevier, New York. Mazars J., Pijaudier-Cabot G., Continuum damage theory application to concrete. J. Eng. Mech. ASCE 115, Ortiz M., A constitutive theory for the inelastic behavior of concrete. Mech. of Materials 4, Peerlings R.H.J., de Borst R., Brekelmans W.A.M., de ree J.H.P., Gradient-enhanced damage for quasi-brittle materials. Int. J. Num. Meth. Eng. 39, Pijaudier-Cabot G., Non local damage. In: Mülhaus H.-B. Ed.), Continuum Models for Materials with Microstructure, Wiley, pp Pijaudier-Cabot G., Bažant Z.P., Non-local damage theory. J. Eng. Mech. ASCE 113,
22 22 C. Comi Ramtani S., Berthaud Y., Mazars J., Orthotropic behavior of concrete with directional aspects: modelling and experiments. Nuclear Eng. Design 133, Rice J.R., The localisation of plastic deformation. In: Koiter W.T. Ed.), Theoretical and Applied Mechanics, North-Holland, Amsterdam, Rocco C., Guinea G.., Planas J., Elices M., The effect of the boundary conditions on the cylinder splitting strength. In: Wittman F.H. Ed.), Fracture Mechanics of Concrete Structures, AEDIFICATIO Publishers, Freiburg, pp Simo J.C., Ju J.W., Strain- and stress-based continuum damage models I. Formulation. Int. J. Solids Structures 23, Steinmann P., Stein E., Finite element localisation analysis of micropolar strength degrading materials. In: Mang H., Bicanic N., de Borst R. Eds.), Computer Modelling of Concrete Structures, pp
strain 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 informationCDPM2: A damage-plasticity approach to modelling the failure of concrete
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:1307.6998v1 [cond-mat.mtrl-sci] 26
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationDEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS
DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS Mohsen Safaei, Wim De Waele Ghent University, Laboratory Soete, Belgium Abstract The present work relates to the
More informationCode_Aster. Reaction of the steel subjected to corrosion
Titre : Comportement de l acier soumis à la corrosion Date : 16/0/01 Page : 1/1 Reaction of the steel subjected to corrosion Summarized: This documentation the model presents steel reaction of subjected
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 informationVersion default Titre : Modèle d endommagement de MAZARS Date : Page : 1/15. Responsable : HAMON François Clé : R Révision : 3ec eb
Titre : Modèle d endommagement de MAZARS Date : Page : 1/15 Model of damage of MAZARS Summary: This documentation presents the model of behavior of MAZARS who allows to describe the behavior rubber band-endommageable
More informationNumerical modeling of standard rock mechanics laboratory tests using a finite/discrete element approach
Numerical modeling of standard rock mechanics laboratory tests using a finite/discrete element approach S. Stefanizzi GEODATA SpA, Turin, Italy G. Barla Department of Structural and Geotechnical Engineering,
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 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 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 informationCard 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 informationTable of Contents. Foreword... xiii Introduction... xv
Foreword.... xiii Introduction.... xv Chapter 1. Controllability of Geotechnical Tests and their Relationship to the Instability of Soils... 1 Roberto NOVA 1.1. Introduction... 1 1.2. Load control... 2
More 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 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 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 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 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 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 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 informationExperiments and Numerical Simulations on Stress-State-Dependence of Ductile Damage Criteria
Experiments and Numerical Simulations on Stress-State-Dependence of Ductile Damage Criteria Michael Brünig, Steffen Gerke and Daniel Brenner Abstract The paper deals with a series of new experiments and
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 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 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 informationA FINITE ELEMENT MODEL FOR SIZE EFFECT AND HETEROGENEITY IN CONCRETE STRUCTURES
A FINITE ELEMENT MODEL FOR SIZE EFFECT AND HETEROGENEITY IN CONCRETE STRUCTURES Roque Luiz Pitangueira 1 and Raul Rosas e Silva 2 1 Department of Structural Engineering -Federal University of Minas Gerais
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationReference material Reference books: Y.C. Fung, "Foundations of Solid Mechanics", Prentice Hall R. Hill, "The mathematical theory of plasticity",
Reference material Reference books: Y.C. Fung, "Foundations of Solid Mechanics", Prentice Hall R. Hill, "The mathematical theory of plasticity", Oxford University Press, Oxford. J. Lubliner, "Plasticity
More informationSimulations of necking during plane strain tensile tests
J. Phys. IV France 134 (2006) 429 434 C EDP Sciences, Les Ulis DOI: 10.1051/jp4:2006134066 Simulations of necking during plane strain tensile tests F. Dalle 1 1 CEA/DAM Île de France, BP. 12, 91680 Bruyères-le-Châtel,
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 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 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 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 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 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 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 informationMasonry strength domain by homogenization in generalized plane state
Masonry strength domain by homogenization in generalized plane state Enrica Bernardini, Nicola Cavalagli, Federico Cluni, Vittorio Gusella Department of Civil and Environmental Engineering, University
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 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 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 informationMulti-scale representation of plastic deformation in fiber-reinforced materials: application to reinforced concrete
!!1 Multi-scale representation of plastic deformation in fiber-reinforced materials: application to reinforced concrete Abstract Here we present a multi-scale model to carry out the computation of brittle
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 informationFRACTURE MECHANICS OF COMPOSITES WITH RESIDUAL STRESSES, TRACTION-LOADED CRACKS, AND IMPERFECT INTERFACES
Proc. 2 nd ESIS TC4 Conference on Polymers and Composites, in press, 1999 Author prepared reprint FRACTURE MECHANICS OF COMPOSITES WITH RESIDUAL STRESSES, TRACTION-LOADED CRACKS, AND IMPERFECT INTERFACES
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 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 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 informationExcavation Damaged Zone Modelling in Claystone with Coupled Second Gradient Model
Excavation Damaged Zone Modelling in Claystone with Coupled Second Gradient Model Frédéric Collin * and Benoît Pardoen ** Argenco Department University of Liège 1, Chemin des Chevreuils B-4000 Liège, Belgium
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 informationTransactions on Modelling and Simulation vol 9, 1995 WIT Press, ISSN X
Elastic-plastic model of crack growth under fatigue using the boundary element method M. Scibetta, O. Pensis LTAS Fracture Mechanics, University ofliege, B-4000 Liege, Belgium Abstract Life of mechanic
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 informationSize effect in the strength of concrete structures
Sādhanā Vol. 27 Part 4 August 2002 pp. 449 459. Printed in India Size effect in the strength of concrete structures B L KARIHALOO and Q Z XIAO Division of Civil Engineering School of Engineering Cardiff
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 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 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 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 informationSTRAIN LOCALIZATION AS BIFURCATION ELASTO-PLASTIC SOFTENING MATERIALS
Fracture Mechanics of Concrete Structures, Proceedings FRAMCOS-2, edited by Folker H. Wittmann, AEDIFICA TIO Publishers, D-79104 Frei burg (1995) STRAIN LOCALIZATION AS BIFURCATION ELASTO-PLASTIC SOFTENING
More informationA C 0 FINITE ELEMENT IN GRADIENT PLASTICITY FOR LOCALIZED FAILURE MODES ANALYSIS
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
More informationDetermination of size-independent specific fracture energy of concrete from three-point bend and wedge splitting tests
Magazine of Concrete Research, 3, 55, No. 2, April, 133 141 Determination of size-independent specific fracture energy of concrete from three-point bend and wedge splitting tests H. M. Abdalla and B. L.
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 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 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 informationSiping Road 1239, , Shanghai, P.R. China
COMPARISON BETWEEN LINEAR AND NON-LINEAR KINEMATIC HARDENING MODELS TO PREDICT THE MULTIAXIAL BAUSCHINGER EFFECT M.A. Meggiolaro 1), J.T.P. Castro 1), H. Wu 2) 1) Department of Mechanical Engineering,
More informationComposite models 30 and 131: Ply types 0 and 8 calibration
Model calibration Composite Bi-Phase models 30 and 3 for elastic, damage and failure PAM-CRASH material model 30 is for solid and 3 for multi-layered shell elements. Within these models different ply types
More informationMESOSCOPIC MODELLING OF MASONRY USING GFEM: A COMPARISON OF STRONG AND WEAK DISCONTINUITY MODELS B. Vandoren 1,2, K. De Proft 2
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm MESOSCOPIC MODELLING OF MASONRY USING GFEM: A COMPARISON OF STRONG AND WEAK DISCONTINUITY
More informationELASTOPLASTICITY THEORY by V. A. Lubarda
ELASTOPLASTICITY THEORY by V. A. Lubarda Contents Preface xiii Part 1. ELEMENTS OF CONTINUUM MECHANICS 1 Chapter 1. TENSOR PRELIMINARIES 3 1.1. Vectors 3 1.2. Second-Order Tensors 4 1.3. Eigenvalues and
More information