MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 343842 10.1155/2013/343842 343842 Research Article XFEM for Thermal Crack of Massive Concrete http://orcid.org/0000-0003-2308-7578 Liu Guowei Hu Yu Li Qingbin Zuo Zheng Cen Song State Key Laboratory of Hydroscience and Engineering Tsinghua University Beijing 100084 China tsinghua.edu.cn 2013 9 10 2013 2013 28 06 2013 02 09 2013 02 09 2013 2013 Copyright © 2013 Guowei Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Thermal cracking of massive concrete structures occurs as a result of stresses caused by hydration in real environment conditions. The extended finite element method that combines thermal fields and creep is used in this study to analyze the thermal cracking of massive concrete structures. The temperature field is accurately simulated through an equivalent equation of heat conduction that considers the effect of a cooling pipe system. The time-dependent creep behavior of massive concrete is determined by the viscoelastic constitutive model with Prony series. Based on the degree of hydration, we consider the main properties related to cracking evolving with time. Numerical simulations of a real massive concrete structure are conducted. Results show that the developed method is efficient for numerical calculations of thermal cracks on massive concrete. Further analyses indicate that a cooling system and appropriate heat preservation measures can efficiently prevent the occurrence of thermal cracks.

1. Introduction

In massive hardening concrete, the thermal gradients under the heat of hydration during hardening can cause thermal cracking. The thermal cracking of concrete structures is a serious concern in construction. Some cracks with wide openings can particularly result in durability problems.

To address this concern, various numerical models have been proposed to investigate the nature of thermal cracks for massive concrete by means of finite element calculations. de Borst and van den Boogaard  studied deformation and cracking in early-age concrete modeled by finite element method (FEM). Xiang et al.  conducted thermoelastic analysis by examining time-dependent evolutions of heat release caused by hydration and material properties. Mazars and Bournazel , Waller et al. , and Lackner and Mang  performed thermal-elastic (or plastic) analyses of massive concrete without considering creep strain. de Schutter  simulated crack initiation and propagation by employing a simple stress-based cracking criterion and a smeared crack model based on the heat of hydration in hardening massive concrete elements. Benboudjema and Torrenti  developed a numerical model to predict early-age cracking of massive concrete through Kelvin-Voigt elements and an elastic damage model that considered autogenous and thermal shrinkage. Briffaut et al. [8, 9] examined the early-age cracking of massive concrete structures caused by thermal restrained shrinkage by using the thermal active restrained shrinkage ring test and FE methods. Buffo-Lacarrière et al.  proposed FE modeling to predict thermal cracking through the evolution of the damage variable law, which considered creep and damage behavior based on nonlinear viscoelasticity.

Massive concrete has been widely used in civil engineering, such as in the construction of dams, buildings, and bridges. Traditional finite element methods have been widely applied in the thermal analysis of massive concrete. However, conventional FEM is formulated with continuous media, so additional remeshing is necessary to accurately predict irregular crack propagation . The extended finite element method (XFEM)  enhances conventional FEM capabilities by excluding the mesh requirement to conform to discontinuities. XFEM has been widely used in numerous fields with discontinuous problems, particularly in fracture mechanics, because XFEM is an excellent method of addressing discrete crack propagation in various types of materials. XFEM was applied to thermal problems in  and to shear band problems with thermal effects in . It was applied by Duflot for the first time in steady-state thermoelastic fracture mechanics . Zamani and Eslami [19, 20] investigated thermoelastic fractures by implementing XFEM for dynamic fractures and by improving accuracy through high-order crack tip enrichments. In these papers, the authors employed XFEM to characterize both the displacement discontinuity and the temperature discontinuity. In addition to XFEM, the embedded finite element method (EFEM) , meshfree methods (MMs) , and boundary elements (BEMs)  are also efficient computational methods based on the continuum-mechanics for discrete cracks.

An alternative pathway to simulate crack initiation and propagation offers methods based on the discrete-mechanics such as the particle simulation method  and the particle element method [38, 39]. The methods have been recently developed to discretize the continuum into particles so that crack generation at a contact between two particles is a natural process. Zhao et al.  simulated magma intrusion problems using solid particles for the surrounding rock and fluid particles for the intruded magma, respectively. The particle simulation method can efficiently simulate hydraulic fracturing because the fluid particles are allowed to move through the cracks between the solid particles. In addition, the particle simulation method has also been used to solve a broad range of scientific and engineering problems .

In accordance with the literature review, this paper aims to develop a numerical method based on XFEM to analyze thermal cracks in massive concrete structures. A practical thermal model that considers an artificial water-cooling system is adopted to simulate the temperature field. Time-dependent viscoelastic behavior is implemented for creep. As input parameters for XFEM that combines the temperature field and the viscoelastic constitutive model, the main properties develop over time when thermal cracking is calculated. The reliability of the developed method is verified by a case study on a lift of an arch dam, which underwent a cold wave during its construction. Subsequently, further analysis on thermal cracking is conducted.

2. XFEM for Thermal Cracking of Massive Concrete 2.1. Briefing on XFEM

Compared with the classical FEM, XFEM provides significant benefits in modeling crack propagation. The crack geometry in XFEM does not need to be aligned with the element edges. XFEM is based on the partition of unity concept and introduces additional degrees of freedom which are tied to the nodes of the elements intersected by the crack [12, 13].

XFEM displacement approximation around the crack uh(x) is enriched with the step function to model the crack surface and with asymptotic displacement fields to model crack tips. The displacement approximation takes the following form : (1)uh(x)=nNnNn(x)ui+jNcrNj(x)h(x)aj+kNtipNk(l=14ψl(x)bkl), where the first term is the classical FE approximation, Nn is the FE shape function, and ui is the vector of the regular degrees of nodal freedom. The other terms are the enriched terms, where Ncr and Ntip stand for the set of nodes that have crack faces and crack tips in their support domain, respectively. aj and bkl denote the vectors of additional degrees of nodal freedom for modeling crack faces and crack tips, respectively.

Step function h(x) simulates the displacement at the crack location with the help of the signed distance function ϕ(x) to the crack, which is also called the normal level set [13, 15]: (2)h(x)=sign(ϕ(x)).

The four following branch functions, which define the discontinuity near the crack tip, enrich the nodes with a distance to the front inferior to a prescribed enrichment radius: (3)ψ(x)=[rsinθ2,rcosθ2,rsinθ2cosθ2,rcosθ2cosθ], where r and θ are polar coordinates in the local crack-tip coordinate system.

Duflot  applied the same procedure to temperature enrichment. The step function h(x) can account for the temperature jump. The leading term of the asymptotic field for the temperature of an adiabatic crack can be written as (4)T=-KTλ2rπsin(θ2), where coefficient KT provides the strength of flux singularity at the crack tip. The temperature field is discretized similar to the displacement field but with only the second branch function (4), which is the only discontinuous branch function. The approximation of the temperature field Th(x) is then written as (5)Th(x)=nNNn(x)Ti+j=NcrNj(x)h(x)cj+kNtipNkψ2(x)dk.

As in standard FEM, performing numerical integration over the element domain is necessary to compute the element stiffness matrix. The subdomain method , in which the crack is one of the subdomain boundaries to carry out the numerical integration, is adopted for the elements that are cut by the crack.

2.2. Equivalent Equation of Heat Conduction

An isotropic homogenous domain Ω bounded by boundary Γ is considered. The thermal properties of this domain are independent of temperature, so the equation of thermal conduction is (6)q=-λT-q+ρcT˙=Q,˙ where q is the heat flux, λ is the thermal conductivity, T is the temperature, ρ is the density, c is the specific heat, and Q˙ is the heat source.

For massive concrete structures, the hydration heat of concrete, which varies with age, is the main heat source. Temperature control is one of the critical problems in the construction of massive concrete, and embedding water-cooling pipes in massive concrete is one of the key methods to address this problem. Zhu  suggested that heat removal through a cooling system can be regarded as a negative heat source. Therefore, the heat source consists of the following two parts in consideration of the cooling pipe system: (7)Q=ρcξ(t+Δt2)+ρc[T(t)-Tw]φ(t+Δt2). In the equation, the first part is the heat source caused by hydration, and the second part is the “cool source” caused by the water cooling system. ξ(t) is the adiabatic temperature rise of concrete at time t, Δt is the increment step, Tw is the temperature of the cooling water inlet, and φ(t) is the cooling function illustrated by Yang et al. . Yang et al.  emphasized that although the temperature gradient near the cooling pipes cannot be obtained, uniformly dispersing heat throughout the solution domain is feasible.

For the adiabatic temperature rise of concrete, Zhu  proposed a compound exponential formula that is convenient for mathematical operation and is in accordance with experimental results. The hydration model can be written as follows: (8)ξ(t)=i=1nξi(1-e-mit), where ξi and mi are the parameters that can be obtained through the optimization method or the trial method whose details can be found in the literature .

2.3. Thermal Boundary Conditions

When the concrete is in contact with air, we assume that the heat flux on the concrete surface is proportional to the difference between concrete surface temperature Tsurf and air temperature Tair. Therefore, the thermal boundary conditions at the concrete surface can be written as (9)-λTn=hT(Tsurf-Tair), where hT is the surface heat convection coefficient and n represents the normal direction of outer surfaces.

2.4. Viscoelastic Creep Model

Early-age concrete exhibits high creep, which reduces the stresses significantly more than any other influencing parameter, as well as assisting in mitigating thermal cracking by reducing thermal stress. Bažant  briefly summarized the method of calculating concrete creep.

The mechanism of concrete creep is not entirely clear, but the viscoelastic constitutive model composed of spring and dashpots can describe its performance. In this study, creep is included in time-dependent viscoelastic behavior. The constitutive relationship of viscoelastic materials can be represented in an integral form as follows : (10)σij(t)=20tG(t-t)εijdevtdt+30tK(t-t)εkk(t)tdt, where G and K are the shear and bulk moduli, respectively, which are functions of reduced time t. Superscript dots denote differentiation with respect to time t, and εijdev and εkk are the mechanical deviatory and volumetric strains, respectively. The bulk (K) and shear (G) moduli can be defined individually through Prony series representation, which is expressed as (11)G(t)=G0[1-i=1nGgi(1-e-t/τiG)],K(t)=K0[1-i=1nKki(1-e-t/τiK)], where gi, τiG, ki, and τiK are material constants. Variable numbers of Prony series parameters end in nK and nG. G0=G(0) and K0=K(0) are the instantaneous relaxation moduli obtained by the instantaneous Young’s moduli as follows: (12)G0=E02·(1+υ),K0=E03·(1-2υ), where υ is the Poisson’s ratio.

To obtain the relaxation modulus, the normalized shear and bulk moduli can be defined as functions of the creep coefficient, which is expressed as follows : (13)G(t)G0=1(1+φ(t,t0)),K(t)K0=1(1+φ(t,t0)), where φ(t,t0) is the creep coefficient, which can be calculated as follows : (14)φ(t,t0)=φ(t,t0)stan·β(t0)β(h0)β(fck)·βTβfβα with (15)φ(t,t0)stan=3.7419(1-(11520-t11520)21.528)0.263, where β(t0), β(h0), β(fck), βT, βf, and βα represent a series of correction coefficients that were described in detail in the literature . Through nonlinear least-squares fit, the Prony series parameters in (11) are calculated.

The developed XFEM method that combines the equivalent equation of heat conduction and the viscoelastic creep model, as previously described, can be used in the numerical calculation of the thermal cracking of massive concrete. Therefore, considering hydration and creep, we can simulate a few discrete thermal cracks on a massive concrete structure.

3. Fracture Criteria and Evolution of Material Properties

The temperature gradient in massive concrete leads to tensile stress, which in turn causes crack initiation and evolution. In modeling thermal cracking evolution, an appropriate law on crack initiation and damage evolution should be defined. The mechanical properties of concrete are essential to predict the development of thermal stress in massive concrete elements.

3.1. Criterion for Initiation

As a possible first approximation, the maximum principal stress criterion can be applied to simulate the onset of thermal cracks. This criterion can be defined as follows: (16)f={σmaxσmax0}, where f is the maximum principle stress ratio and σmax and σmax0 denote the maximum principle stress and maximum allowable principle stress (tensile strength), respectively. The symbol represents the so-called Macaulay brackets, which enable σmax to have an alternative value of 0 (when σmax<0) and σmax (when σmax0). The crack initiates when f=1.

3.2. Criterion for Evolution

In this study, the analysis of crack growth is based on the approach of linear elastic fracture mechanics and is conducted with the framework of the extended-FEM method. The power law criterion proposed by Wu and Reuter  can be extended to 3D problems as follows : (17)(GIGIc)α+(GIIGIIc)β+(GIIIGIIIc)χ=1, where α, β, and χ are parameters with values of 1 to 2, respectively. GI, GII, and GIII denote the strain energy release rate parameters of mode I, mode II, and mode III components, respectively; GIc, GIIc, and GIIIc are the critical values of the strain energy release rate (the fracture toughness) for the three fracture modes. The law forms an energy envelope surface, and if the energy release rate exceeds this surface, the crack propagates.

3.3. Evolution of Material Properties

During the construction period in which concrete changes from an almost liquid state to a solid state, most of the mechanical properties rapidly vary in relation to the age of concrete. Among these properties, modulus of elasticity, tensile strength, and fracture toughness are the key parameters used in analyzing thermal cracking.

Degree of hydration is a fundamental parameter that describes concrete behavior during hardening. de Schutter and Taerwe [54, 55] proposed a degree-of-hydration-based description of these material properties as follows: (18)E(γ)E(γ=1)=(γ-γ01-γ0)a0,σmax0(γ)σmax0(γ=1)=(γ-γ01-γ0)b0,GcGc(γ=1)=(γ-γ01-γ0)c0, where E(γ), σmax0(γ), and Gc(γ) are the Young’s modulus, tensile strength, and critical values of strain energy release rate at degree of reaction γ, respectively; and E(γ=1), σmax0(γ=1), and Gc(γ=1) correspond to these values when γ=1. γ0, a0, b0, and c0 are parameters that depend on the concrete composition.

A practical estimation of the degree of hydration at each moment t can be expressed as the ratio between the accumulated heat of hydration at time t and that at the end of reaction, which can be calculated through the adiabatic temperature rise: (19)γ(t)=0tξ(τ)dτ0tendξ(τ)d(τ), where γ(t) is the degree of hydration at time t and tend is the time at end of hydration reaction.

The coefficient of thermal expansion of concrete is another key parameter for thermal cracking analysis. Its measured values tend to be stable within two days , so a constant for it is adopted in this paper.

4. Numerical Examples 4.1. Introduction

A 286 m high parabolic double-curvature arch dam under construction is located in Southwest China at an altitude of 610 m above sea level (level of crest). The dam consists of 31 monoliths, each containing a series of lifts. The 10th lift of the 13th monolith was poured on September 5, 2009 and underwent a cold wave before the next lift was poured. This lift was selected for the thermal cracking analysis. The 9th lift was also used in our calculation as the lower boundary of the 10th lift.

Figure 1 illustrates the finite element mesh of the two lifts with cast corners (they were used to facilitate the construction and increase the rigidity of the bottom of the dam) at downstream; each lift (1.5 m high) contains four-layer elements. Each element in Figure 1 has an edge length of about 0.984 m and a thickness of 0.375 m. The numerical model includes a total of 15488 elements, in which 5544 elements (the 10th lift) are used in the XFEM mesh.

Element mesh of the two lifts (9th and 10th).

4.2. Boundary Conditions

During construction, the 10th lift of the monolith had been standing for 86 days before the next lift was poured. The increasing upstream and downstream surface of the poured concrete slabs and the top surface of the latest slab were assumed to be directly exposed to the environment, where heat transfer by convention occurred. Because the other lifts under the 9th lift and in the other monoliths were not considered, the boundary condition of the lower surface of the 9th lift was assumed to be adiabatic with fixed vertical displacement.

The calculations correspond to a period of 128 days (the duration between starting the 9th lift and starting the 11th lift), and the 10th lift started to be poured 42 days after the 9th lift was poured. The daily mean air temperature curve at the dam site, as shown in Figure 2, was used in this analysis.

Development of air temperature.

4.3. Properties of Materials

According to , the hydration model in the engineering case analyzed in this paper may be written as follows: (20)ξ(t)=16.4(1-e-0.516t)+10.9(1-e-0.052t). The equation is in accordance with the experimental data. Each concrete slab contains one-layer staggered heterogeneous iron cooling pipes in the middle for artificial cooling; related property values are illustrated in [46, 57].

Equations (18) to (19) show that the Young’s modulus, the maximum allowable principle stress, and the critical values of strain energy release rate can be written as time-dependent functions (E(t), σmax0(t), and Gc(t)) when (20) is used. In the simplest case, the fracture process can be isotropic, and parameters GIc(t)=GIIc(t)=GIIIc(t) and α=β=χ=1. The main properties of concrete are presented in Table 1. Figure 3 shows the inclusion of the creep coefficient in viscoelastic material behavior during the cracking analysis for the 10th lift.

Main properties of concrete.

Property Values Note
Density ρ 2663.0 kg/m3 Constant
Thermal conductivity λ 184.9 kJ/m·day °C Constant
Specific heat c 0.86 kJ/kg °C Constant
Heat convection coefficient hT 1200.0 kJ/m2·day °C Constant
a 0 1.40 Constant
b 0 0.62 Constant
c 0 0.88 Constant
γ 0 0.25 Constant
Elastic modulus E (28 d) 40 GPa Time-dependent
Poisson’s ratio υ 0.167 Constant
Coefficient of thermal expansion α 7 × 10−6/°C Constant
Allowable maximum principle σmax0 (28 d) 1.98 MPa Time-dependent
Critical energy release rate for model I GIc (28 d) 47.92 N/m Time-dependent

Development of the creep coefficient.

4.4. Results and Discussions

Accurate temperature fields are required to calculate thermal cracking. Figure 4 compares the temperature histories of the concrete in the calculated results and the actual measured temperature from thermometers buried at middle height of the lift. Overall, the calculated results obtained from the methods used in this study agree with the measured data. A sudden temperature drop occurred in the position of the buried thermometers five days after the cold wave, and an inflection point was observed in the computed curve two days in advance. These conditions can be related to the positions of the thermometers and the thermal conductivity.

Development of concrete temperature.

Both elastic and viscoelastic models were applied in the stress analysis. Four elements of different heights at the center of the slab (see Figure 1) were selected, with the time-dependent maximum principal stress at the centroid of these elements shown in Figure 5. Overall creep can reduce the stress levels by approximately 20%. At the beginning, the bottom element bears low tensile stress levels because of the effects from the lower lift. The three other elements bear compressive stresses because of the restraint of thermal dilatation and the rise to their peaks seven days after the lift casting. After the temperature decreases in the core as a result of artificial water cooling, compressive stresses decrease gradually, and then, tensile stress occurs. Affected by the ambient temperature, the stresses fluctuate over time and increase suddenly as a result of the cold wave. The stress of the surface element that contacts with air has the largest varied amplitude and exceeds the tensile strength eventually.

Development of stresses of different heights.

The cold wave results in a thermal gradient (between the core and the skin of the lift), which causes tensile stress that is particularly serious on the surface. The cracks occur at discrete locations where the maximum principle stress exceeds the allowable values and propagates based on the criterion. The cracks simulated by XFEM and observed on the slab are presented in Figure 6. This figure illustrates that the numerical results are generally consistent with real crack occurrence. Most of the main cracks are initiated at the edges of the lift where high-temperature gradients occur, and this condition favors crack initiation. The evolution paths of the simulated cracks are almost close to the lines because of the isotropic assumption for the material and the crack evolution law. The discrepancy between the simulated and real crack directions in the middle of the lift surface is a result of excluding the other lifts (except the 9th and 10th) and the foundation in the calculation. In addition, crack crossing cannot be simulated because of the limitations of XFEM. The numerical simulation indicates that ambient temperature is among the main reasons for the thermal cracking of massive concrete structures.

Comparison between simulated and actual thermal cracks: (a) simulated via XFEM and (b) actual crack occurrence.

The casting temperature obviously affects the cracking risk if no cooling system is used . The concrete pouring temperature in the preceding analysis is 6.7°C, and then, we gradually increase the casting temperature to the ambient one (30.3°C). The maximum principle stresses with different pouring temperatures are presented in Figure 7. The result indicates that an increase in casting temperature does not exert significant effects on the stresses during the early stage because the embedded water-cooling system decrease the concrete temperature effectively. A negligible discrepancy is observed between these stresses at a later stage.

Development of stresses with different pouring temperatures.

Heat preservation measures can be adopted during the intermission between the pouring of two lifts to reduce the thermal gradient of concrete. In this study, we used different heat convection coefficients to simulate different levels of thermal insulation effects. The stress evolution with a series of heat convection coefficients is illustrated in Figure 8. The stress levels decrease significantly with a decrease in the coefficients, particularly from 400 kJ/m·day °C to 200 kJ/m·day °C. In contrast to the case during the early stage, a small convection coefficient assumes a critical function when the cold wave arrives. The results demonstrate that appropriate heat insulation measures can effectively reduce the risk of thermal cracking.

Development of stresses with different heat convection coefficients.

5. Conclusions

This study presents an effective XFEM method to simulate the thermal cracking of massive concrete. The temperature field in the slab can be obtained by the equivalent equation of heat conduction, which includes the hydration heat of concrete and the effect of cooling pipe systems. The time-dependent viscoelastic behavior presented by the Prony series for concrete is included in the prediction of stresses to describe the stress relaxation that results from creep. The Young’s modulus, tensile strength, and critical values of strain energy release rate which evolve over time are the key parameters that affect the likelihood of thermal cracking in massive concrete. Thermal cracking of massive concrete can be described through a combination of XFEM with the temperature field and the viscoelastic creep model.

Compared with existing numerical methods to simulate thermal cracks on massive concrete, the developed XFEM method in this study has advantages in efficiently simulating multiple discrete thermal crack propagations without presupposing crack path and remeshing. Therefore, this method can be conveniently applied to analyze thermal cracking for massive concrete structures during the design and construction phases.

The precision of the developed method was verified by an engineering application. The simulated and actual results were in agreement, as shown in two comparisons: the comparison between the calculated time-dependent temperature curves and the measured ones and the comparison between the numerical simulated cracks and the actual observed crack occurrence. If the structure is embedded with a cooling system and the ambient temperature is appropriate, the casting temperature for concrete exerts a limited effect on thermal cracking risk that is different for the case without a cooling system. However, heat preservation measures can exert a significant effect on thermal crack control for massive concrete, particularly when a cold wave occurs. These measures deserve attention during the construction of massive concrete structures.

Acknowledgments

This study was financially supported by National Basic Research Program of China (2013CB035902), Research Project of State Key Laboratory of Hydroscience and Engineering of Tsinghua University (2011-KY-4), and the National Natural Science Foundation of China (51279087).

de Borst R. van den Boogaard A. H. Finite-element modeling of deformation and cracking in early-age concrete Journal of Engineering Mechanics 1994 120 12 2519 2534 2-s2.0-0028669682 Xiang Y. Zhang Z. He S. Dong G. Thermal-mechanical analysis of a newly cast concrete wall of a subway structure Tunnelling and Underground Space Technology 2005 20 5 442 451 2-s2.0-20844456448 10.1016/j.tust.2005.02.005 Mazars J. Bournazel J. P. Modelling of damage processes due to volumic variations for maturing and matured concrete Proceedings of the International RILEM Conference on Concrete: From Material to Structure 1998 43 54 Waller V. d'Aloïa L. Cussigh F. Lecrux S. Using the maturity method in concrete cracking control at early ages Cement and Concrete Composites 2004 26 5 589 599 2-s2.0-2542505835 10.1016/S0958-9465(03)00080-5 Lackner R. Mang H. A. Chemoplastic material model for the simulation of early-age cracking: from the constitutive law to numerical analyses of massive concrete structures Cement and Concrete Composites 2004 26 5 551 562 2-s2.0-2542449469 10.1016/S0958-9465(03)00071-4 de Schutter G. Finite element simulation of thermal cracking in massive hardening concrete elements using degree of hydration based material laws Computers and Structures 2002 80 27–30 2035 2042 2-s2.0-0036857433 10.1016/S0045-7949(02)00270-5 Benboudjema F. Torrenti J. M. Early-age behaviour of concrete nuclear containments Nuclear Engineering and Design 2008 238 10 2495 2506 2-s2.0-49849092377 10.1016/j.nucengdes.2008.04.009 Briffaut M. Benboudjema F. Torrenti J. M. Nahas G. A thermal active restrained shrinkage ring test to study the early age concrete behaviour of massive structures Cement and Concrete Research 2011 41 1 56 63 2-s2.0-78649329010 10.1016/j.cemconres.2010.09.006 Briffaut M. Benboudjema F. Torrenti J. M. Nahas G. Numerical analysis of the thermal active restrained shrinkage ring test to study the early age behavior of massive concrete structures Engineering Structures 2011 33 4 1390 1401 2-s2.0-79951911196 10.1016/j.engstruct.2010.12.044 Buffo-Lacarrière L. Sellier A. Turatsinze A. Escadeillas G. Finite element modelling of hardening concrete: application to the prediction of early age cracking for massive reinforced structures Materials and Structures 2011 44 10 1821 1835 2-s2.0-84856227056 10.1617/s11527-011-9740-y Abdelaziz Y. Hamouine A. A survey of the extended finite element Computers and Structures 2008 86 11-12 1141 1151 2-s2.0-43049165623 10.1016/j.compstruc.2007.11.001 Belytschko T. Black T. Elastic crack growth in finite elements with minimal remeshing International Journal for Numerical Methods in Engineering 1999 45 5 601 620 2-s2.0-0032649861 Moës N. Dolbow J. Belytschko T. A finite element method for crack growth without remeshing International Journal for Numerical Methods in Engineering 1999 46 1 131 150 2-s2.0-0033349534 Sukumar N. Chopp D. L. Moës N. Belytschko T. Modeling holes and inclusions by level sets in the extended finite-element method Computer Methods in Applied Mechanics and Engineering 2001 190 46-47 6183 6200 MR1857695 ZBL1029.74049 2-s2.0-0035860267 10.1016/S0045-7825(01)00215-8 Fries T.-P. Belytschko T. The extended/generalized finite element method: an overview of the method and its applications International Journal for Numerical Methods in Engineering 2010 84 3 253 304 MR2732652 ZBL1202.74169 2-s2.0-78650831150 10.1002/nme.2904 Merle R. Dolbow J. Solving thermal and phase change problems with the eXtended finite element method Computational Mechanics 2002 28 5 339 350 2-s2.0-0036577760 10.1007/s00466-002-0298-y Areias P. M. A. Belytschko T. Two-scale method for shear bands: thermal effects and variable bandwidth International Journal for Numerical Methods in Engineering 2007 72 6 658 696 2-s2.0-35448983166 10.1002/nme.2028 Duflot M. The extended finite element method in thermoelastic fracture mechanics International Journal for Numerical Methods in Engineering 2008 74 5 827 847 MR2406405 ZBL1195.74170 2-s2.0-42949170644 10.1002/nme.2197 Zamani A. Eslami M. R. Implementation of the extended finite element method for dynamic thermoelastic fracture initiation International Journal of Solids and Structures 2010 47 10 1392 1404 2-s2.0-77949492664 10.1016/j.ijsolstr.2010.01.024 Zamani A. Gracie R. Eslami M. R. Higher order tip enrichment of eXtended Finite Element Method in thermoelasticity Computational Mechanics 2010 46 6 851 866 MR2680104 ZBL05797418 2-s2.0-77956743302 10.1007/s00466-010-0520-2 Dvorkin E. N. Cuitino A. M. Gioia G. Finite elements with displacement interpolated embedded localization lines insensitive to mesh size and distortions International Journal for Numerical Methods in Engineering 1990 30 3 541 564 2-s2.0-0025473074 Linder C. Armero F. Finite elements with embedded strong discontinuities for the modeling of failure in solids International Journal for Numerical Methods in Engineering 2007 72 12 1391 1433 MR2372143 ZBL1194.74431 2-s2.0-37249060701 10.1002/nme.2042 Linder C. Armero F. Finite elements with embedded branching Finite Elements in Analysis and Design 2009 45 4 280 293 MR2508483 2-s2.0-60949083679 10.1016/j.finel.2008.10.012 Belytschko T. Lu Y. Y. Gu L. Element-free Galerkin methods International Journal for Numerical Methods in Engineering 1994 37 2 229 256 10.1002/nme.1620370205 MR1256818 ZBL0796.73077 2-s2.0-0028259955 Liu W. K. Jun S. Zhang Y. F. Reproducing kernel particle methods International Journal for Numerical Methods in Fluids 1995 20 8-9 1081 1106 10.1002/fld.1650200824 MR1333917 ZBL0881.76072 2-s2.0-0029102512 Atluri S. N. The Meshless Local Petrov-Galerkin (MLPG) Method 2002 Encino, Calif, USA Tech Science Press Bordas S. Rabczuk T. Zi G. Three-dimensional crack initiation, propagation, branching and junction in non-linear materials by an extended meshfree method without asymptotic enrichment Engineering Fracture Mechanics 2008 75 5 943 960 2-s2.0-36849039473 10.1016/j.engfracmech.2007.05.010 Aliabadi M. H. Boundary element formulations in fracture mechanics Applied Mechanics Reviews 1997 50 2 83 96 2-s2.0-0031077274 Pan E. Yuan F. G. Boundary element analysis of three-dimensional cracks in anisotropic solids International Journal for Numerical Methods in Engineering 2000 48 2 211 237 Sfantos G. K. Aliabadi M. H. Multi-scale boundary element modelling of material degradation and fracture Computer Methods in Applied Mechanics and Engineering 2007 196 7 1310 1329 2-s2.0-33750941585 10.1016/j.cma.2006.09.004 Simpson R. N. Bordas S. P. A. Trevelyan J. Rabczuk T. A two-dimensional Isogeometric boundary element method for elastostatic analysis Computer Methods in Applied Mechanics and Engineering 2012 209–212 87 100 MR2877957 ZBL1243.74193 2-s2.0-83455181910 10.1016/j.cma.2011.08.008 Zhao C. Hobbs B. E. Ord A. Fundamentals of Computational Geoscience: Numerical Methods and Algorithms 2009 Berlin, Germany Springer Zhao C. Hobbs B. E. Ord A. Hornby P. Peng S. Liu L. Particle simulation of spontaneous crack generation problems in large-scale quasi-static systems International Journal for Numerical Methods in Engineering 2007 69 11 2302 2329 2-s2.0-33947123334 10.1002/nme.1851 Zhao C. Hobbs B. E. Ord A. Peng S. Liu L. An upscale theory of particle simulation for two-dimensional quasi-static problems International Journal for Numerical Methods in Engineering 2007 72 4 397 421 2-s2.0-35348956367 10.1002/nme.2018 Zhao C. Nishiyama T. Murakami A. Numerical modelling of spontaneous crack generation in brittle materials using the particle simulation method Engineering Computations 2006 23 5 566 584 2-s2.0-33745869489 10.1108/02644400610671144 Zhao C. Hobbs B. E. Ord A. Robert P. A. Hornby P. Peng S. Phenomenological modelling of crack generation in brittle crustal rocks using the particle simulation method Journal of Structural Geology 2007 29 6 1034 1048 2-s2.0-34248507208 10.1016/j.jsg.2007.02.012 Zhao C. Hobbs B. E. Ord A. Peng S. Critical contact stiffness concept and simulation of crack generation in particle models of large length-scales Computers and Geotechnics 2009 36 1-2 81 92 2-s2.0-56949099458 10.1016/j.compgeo.2008.03.001 Oñate E. O. Idelsohn S. R. del Pin F. Aubry R. The particle finite element method: an overview International Journal of Computational Methods 2004 1 2 267 307 Idelsohn S. R. Oñate E. del Pin F. Calvo N. Fluid-structure interaction using the particle finite element method Computer Methods in Applied Mechanics and Engineering 2006 195 17-18 2100 2123 MR2202916 ZBL1178.76230 2-s2.0-32044443135 10.1016/j.cma.2005.02.026 Zhao C. Hobbs B. E. Ord A. Peng S. Particle simulation of spontaneous crack generation associated with the laccolithic type of magma intrusion processes International Journal for Numerical Methods in Engineering 2008 75 10 1172 1193 2-s2.0-52649173325 10.1002/nme.2287 Zhao C. Peng S. Liu L. Hobbs B. E. Ord A. Effective loading algorithm associated with explicit dynamic relaxation method for simulating static problems Journal of Central South University of Technology 2009 16 1 125 130 2-s2.0-60849094782 10.1007/s11771-009-0021-7 Zhao C. Computational simulation of frictional drill-bit movement in cemented granular materials Finite Elements in Analysis and Design 2011 47 8 877 885 2-s2.0-79957508804 10.1016/j.finel.2011.03.003 Zhao C. Peng S. Liu L. Hobbs B. E. Ord A. Computational modeling of free-surface slurry flow problems using particle simulation method Journal of Central South University 2013 20 6 1653 1660 Belytschko T. Moës N. Usui S. Parimi C. Arbitrary discontinuities in finite elements International Journal for Numerical Methods in Engineering 2001 50 4 993 1013 Zhu B.-F. Equivalent equation of heat conduction in mass concrete considering the effect of pipe cooling Journal of Hydraulic Engineering 1991 3 28 34 2-s2.0-0037898841 Yang J. Hu Y. Zuo Z. Jin F. Li Q. Thermal analysis of mass concrete embedded with double-layer staggered heterogeneous cooling water pipes Applied Thermal Engineering 2012 35 1 145 156 2-s2.0-81855207239 10.1016/j.applthermaleng.2011.10.016 Zhu B.-F. Compound exponential formula for variation of thermal and mechanical properties with age of concrete Journal of Hydraulic Engineering 2011 42 1 1 7 2-s2.0-79953083239 Bažant Z. P. Prediction of concrete creep and shrinkage: past, present and future Nuclear Engineering and Design 2001 203 1 27 38 2-s2.0-0343290215 10.1016/S0029-5493(00)00299-5 Belytschko T. Liu W. K. Moran B. Nonlinear Finite Elements for Continua and Structures 2000 Chichester, UK John Wiley & Sons xvi+650 MR1897447 Malm R. Sundquist H. Time-dependent analyses of segmentally constructed balanced cantilever bridges Engineering Structures 2010 32 4 1038 1045 2-s2.0-77649179922 10.1016/j.engstruct.2009.12.030 Guo L. Zhu Y.-M. Zhu M.-D. Zhang X.-Y. The formula design of creep coefficient for high performance concrete Journal of Sichuan University 2010 42 2 75 81 2-s2.0-77951996333 Wu E. M. Reuter R. C. Jr. Crack extension in fiberglass reinforced plastics 1965 275 University of Illinois Goyal V. K. Johnson E. R. Dávila C. G. Irreversible constitutive law for modeling the delamination process using interfacial surface discontinuities Composite Structures 2004 65 3-4 289 305 2-s2.0-3042682704 10.1016/j.compstruct.2003.11.005 de Schutter G. Taerwe L. Degree of hydration-based description of mechanical properties of early age concrete Materials and Structures 1996 29 6 335 344 2-s2.0-0000816515 de Schutter G. Taerwe L. Fracture energy of concrete at early ages Materials and Structures 1997 30 196 67 71 2-s2.0-0000149311 Cusson D. Hoogeveen T. An experimental approach for the analysis of early-age behaviour of high-performance concrete structures under restrained shrinkage Cement and Concrete Research 2007 37 2 200 209 2-s2.0-33846607320 10.1016/j.cemconres.2006.11.005 Zuo Z. Hu Y. Duan Y. Yang J. Simulation of the temperature field in mass concrete with double layers of cooling pipes during construction Journal of Tsinghua University 2012 52 2 186 228 2-s2.0-84860335938 Wu S. Huang D. Lin F.-B. Zhao H. Wang P. Estimation of cracking risk of concrete at early age based on thermal stress analysis Journal of Thermal Analysis and Calorimetry 2011 105 1 171 186 2-s2.0-79959375265 10.1007/s10973-011-1512-y