Modelling heterogeneity of concrete using 2D lattice network for concrete fracture and comparison with AE study

In this paper, numerical modelling of fracture in concrete using two-dimensional lattice model is presented and also a few issues related to lattice modelling technique applicable to concrete fracture are reviewed. A comparison is made with acoustic emission (AE) events with the number of fractured elements. To implement the heterogeneity of the plain concrete, two methods namely, by generating grain structure of the concrete using Fuller's distribution and the concrete material properties are randomly distributed following Gaussian distribution are used. In the first method, the modelling of the concrete at meso level is carried out following the existing methods available in literature. The shape of the aggregates present in the concrete are assumed as perfect spheres and shape of the same in two-dimensional lattice network is circular. A three-point bend (TPB) specimen is tested in the experiment under crack mouth opening displacement (CMOD) control at a rate of 0.0004 mm/sec and the fracture process in the same TPB specimen is modelled using regular triangular 2D lattice network. Load versus crack mouth opening isplacement (CMOD) plots thus obtained by using both the methods are compared with experimental results. It was observed that the number of fractured elements increases near the peak load and beyond the peak load. That is once the crack starts to propagate. AE hits also increase rapidly beyond the peak load. It is compulsory here to mention that although the lattice modelling of concrete fracture used in this present study is very similar to those already available in literature, the present work brings out certain finer details which are not available explicitly in the earlier works.

Verification of the applicability of lattice model to concrete fracture by AE study

Notched three-point bend specimens (TPB) were tested under crack mouth opening displacement (CMOD) control at a rate of 0.0004 mm/s and the entire fracture process was simulated using a regular triangular two-dimensional lattice network only over the expected fracture proces zone width. The rest of the beam specimen was discretised by a coarse triangular finite element mesh. The discrete grain structure of the concrete was generated assuming the grains to be spherical. The load versus CMOD plots thus simulated agreed reasonably well with the experimental results. Moreover, acoustic emission (AE) hits were recorded during the test and compared with the number of fractured lattice elements. It was found that the cumulative AE hits correlated well with the cumulative fractured lattice elements at all load levels thus providing a useful means for predicting when the micro-cracks form during the fracturing process, both in the pre-peak and in the post-peak regimes.

Fracture process zone size and true fracture energy of concrete using acoustic emission

Acoustic emission (AE) energy, instead of amplitude, associated with each of the event is used to estimate the fracture process zone (FPZ) size. A steep increase in the cumulative AE energy of the events with respect to time is correlated with the formation of FPZ. Based on the AE energy released during these events and the locations of the events, FPZ size is obtained. The size-independent fracture energy is computed using the expressions given in the boundary effect model by least squares method since over-determined system of equations are obtained when data from several specimens are used. Instead of least squares method a different method is suggested in which the transition ligament length, measured from the plot of histograms of AE events plotted over the un-cracked ligament, is used directly to obtain size-independent fracture energy. The fracture energy thus calculated seems to be size-independent.

Fracture analysis of stiffened panels under combined tensile, bending, and shear loads

The objective is to present the formulation of numerically integrated modified virtual crack closure integral technique for concentrically and eccentrically stiffened panels for computation of strain-energy release rate and stress intensity factor based on linear elastic fracture mechanics principles. Fracture analysis of cracked stiffened panels under combined tensile, bending, and shear loads has been conducted by employing the stiffened plate/shell finite element model, MQL9S2. This model can be used to analyze plates with arbitrarily located concentric/eccentric stiffeners, without increasing the total number of degrees of freedom, of the plate element. Parametric studies on fracture analysis of stiffened plates under combined tensile and moment loads have been conducted. Based on the results of parametric,studies, polynomial curve fitting has been carried out to get best-fit equations corresponding to each of the stiffener positions. These equations can be used for computation of stress intensity factor for cracked stiffened plates subjected to tensile and moment loads for a given plate size, stiffener configuration, and stiffener position without conducting finite element analysis.

Fracture behavior of concrete–concrete interface using acoustic emission technique

The fracture behavior of concrete–concrete interface is characterized using acoustic emission (AE). Beams of different sizes having jointed interface between two different strengths of concrete are tested. The results of load, displacement, CMOD, AE-events and AE-energy are analyzed. The width of fracture process zone and damage zone are computed using AE-data and are found to be independent of size. It is observed that, as the difference in compressive strength of concrete on either side of interface increases, the load carrying capacity, number of AE-events, AE-energy, width of fracture process zone and damage zone decreases.

The fracture toughness of bulk metallic glasses

Stiffness, strength, and toughness are the three primary attributes of a material, in terms of its mechanical properties. Bulk metallic glasses (BMGs) are known to exhibit elastic moduli at a fraction lower than crystalline alloys and have extraordinary strength. However, the reported values of fracture toughness of BMGs are highly variable; some BMGs such as the Zr-based ones have toughness values that are comparable to some high strength steels and titanium alloys, whereas there are also BMGs that are almost as brittle as silicate glasses. Invariably, monolithic BMGs exhibit no or low crack growth resistance and tend to become brittle upon structural relaxation. Despite its critical importance for the use of BMGs as structural materials, the fracture toughness of BMGs is relatively poorly understood. In this paper, we review the available literature to summarize the current understanding of the mechanics and micromechanisms of BMG toughness and highlight the needs for future research in this important area.

A comparative study of dynamic,ductile fracture initiation in two specimen configurations

In this work a single edge notched plate (SEN(T)) subjected to a tensile stress pulse is analysed, using a 2D plane strain dynamic finite element procedure. The interaction of the notch with a pre-nucleated hole ahead of it is examined. The background material is modelled by the Gurson constitutive law and ductile failure by microvoid coalescence in the ligament connecting the notch and the hole is simulated. Both rate independent and rate dependent material behaviour is considered. The notch tip region is subjected to a range of loading rates j by varying the peak value and the rise time of the applied stress pulse. The results obtained from these simulations are compared with a three point bend (TPB) specimen subjected to impact loading analysed in an earlier work [3] The variation of J at fracture initiation, J(c), with average loading rate j is obtained from the finite element simulations. It is found that the functional relationship between J(c) and j is fairly independent of the specimen geometry and is only dependent on material behaviour.

Effect of copper addition on the fracture and fatigue crack growth behavior of solution heat-treated SUS 304H austenitic steel

Small additions of Cu to the SUS 304H, a high temperature austenitic stainless steel, enhance its high temperature strength and creep resistance. As Cu is known to cause embrittlement, the effect of Cu on room temperature mechanical properties that include fracture toughness and fatigue crack threshold of as-solutionized SUS 304H steel were investigated in this work. Experimental results show a linear reduction in yield and ultimate strengths with Cu addition of up to 5 wt.% while ductility drops markedly for 5 wt.% Cu alloy. However, the fracture toughness and the threshold stress intensity factor range for fatigue crack initiation were found to be nearly invariant with Cu addition. This is because the fracture in this alloy is controlled by the debonding from the matrix of chromium carbide precipitates, as evident from fractography. Cu, on the other hand, remains either in solution or as nano-precipitates and hence does not influence the fracture characteristics. It is concluded that small additions of Cu to 304H will not have adverse effects on its fracture and fatigue behavior. (C) 2010 Elsevier B.V. All rights reserved.

Effect of Prestrain on Fracture Toughness of Rolled Mild Steel with Significant Inclusion Content

The basic method of JIc calculation using a single specimen is discussed. Dokouipil's approach for evaluating the JIc value is extended further and the effect of prestrain on rolled mild steel with significant inclusions is studied using this modified approach. Although this method does not give an accurate value of JIc, it is quite effective for a comparative study. While the fracture toughness of annealed and 7% prestrained materials are about the same, the fracture toughness of 3% prestrained material is significantly lower.

Nonlinear fracture properties of concrete-concrete interfaces

The fracture properties of different concrete-concrete interfaces are determined using the Bazant's size effect model. The size effect on fracture properties are analyzed using the boundary effect model proposed by Wittmann and his co-workers. The interface properties at micro-level are analyzed through depth sensing micro-indentation and scanning electron microscopy. Geometrically similar beam specimens of different sizes having a transverse interface between two different strengths of concrete are tested under three-point bending in a closed loop servo-controlled machine with crack mouth opening displacement control. The fracture properties such as, fracture energy (G(f)), length of process zone (c(f)), brittleness number (beta), critical mode I stress intensity factor (K-ic), critical crack tip opening displacement CTODc (delta(c)), transitional ligament length to free boundary (a(j)), crack growth resistance curve and micro-hardness are determined. It is seen that the above fracture properties decrease as the difference between the compressive strength of concrete on either side of the interface increases. (C) 2010 Elsevier Ltd. All rights reserved.

Alloy design for fracture resistance

Several methods for improving the strength of metallic materials are available and correlations between strength and various microstructural features have been established. The purpose of this paper is to review parallel developments favouring improved fracture resistance. Resistance to fracture in monotonie loading, cyclic loading and when fracture is environment-aided have been considered in steels, aluminium alloys and anisotropic materials. Finally, the question of optimising alloy behaviour is discussed.

A finite element analysis of dynamic fracture initiation by ductile failure mechanisms in a 4340 steel

In some recent dropweight impact experiments [5] with pre-notched bend specimens of 4340 steel, it was observed that considerable crack tunneling occurred in the interior of the specimen prior to gross fracture initiation on the free surfaces. The final failure of the side ligaments happened because of shear lip formation. The tunneled region is characterized by a flat, fibrous fracture surface. In this paper, the experiments of [5] (corresponding to 5 m/s impact speed) are analyzed using a plane strain, dynamic finite element procedure. The Gurson constitutive model that accounts for the ductile failure mechanisms of micro-void nucleation, growth and coalescence is employed. The time at which incipient failure was observed near the notch tip in this computation, and the value of the dynamic J-integral, J d, at this time, compare reasonably well with experiments. This investigation shows that J-controlled stress and deformation fields are established near the notch tip whenever J d , increases with time. Also, it is found that the evolution of micro-mechanical quantities near the notch root can be correlated with the time variation of J d .The strain rate and the adiabatic temperature rise experienced at the notch root are examined. Finally, spatial variations of stresses and deformations are analyzed in detail.

Influence of strain softening on the fracture of plain concrete beams

Fracture behaviour of notched and un-notched plain concrete slender beams subjected to three-point or four-point bending is analyzed through a one-dimensional model, also called Softening Beam Model. Fundamental equations of equilibrium are used to develop the model. The influence of structural size in altering the fracture mode from brittle fracture to plastic collapse is explained through the stress distribution across the uncracked ligament obtained by varying the strain softening modulus. It is found that at the onset of fracture instability, stress at the crack tip is equal to zero. The maximum load and fracture load are found to be different and a unique value for the fracture load is obtained. It is shown that the length of the fracture process zone depends on the value of the strain softening modulus. Theoretical limits for fracture process zone length are also calculated. Several nonlinear fracture parameters, such as, crack tip opening displacement, crack mouth opening displacement and fracture energy are computed for a wide variety of beam specimens reported in the literature and are found to compare very well with experimental and theoretical results. It is demonstrated that by following a simple procedure, both pre-peak and post-peak portions of load versus crack mouth opening displacement curve can be obtained quite accurately. Further, a simple procedure to calculate the maximum load is also developed. The predicted values of maximum load are found to agree well with the experimental values. The Softening Beam Model (SBM), proposed in this investigation is very simple and is based on rational considerations. It can completely describe the fracture process from the beginning of formation of the fracture process zone till the onset of fracture instability.A l'aide d'un modèle unidimensionnel dit ldquoSoftening Beam Modelrdquo (SBM), on analyse le comportement à rupture de poutres élancées pleines entaillées ou non, soumises en flexion en trois ou quatre points. Des équations fondamentales d'équilibre sont utilisées pour développer le modèle. On explique l'influence de la taille du composant sur l'altération du mode de rupture en rupture fragile et en effondrement plastique par la distribution par la distribution des contraintes sur le ligament non fissuré lorsque varie le module d'adoucissement. On trouve que la contrainte à l'extrémité de la fissure est nulle est nulle au début de l'instabilité de la rupture. La charge maximum et la charge à la rupture sont trouvées différentes, et on obtient une valeur unique de la charge à la rupture. On montre que la longueur de la zone concernée par le processus de rupture d'pend de la valeur du module d'adoucissement. On calcule également les limites théoriques de longueur de cette zone. Divers paramètres de rupture non linéaire sont calculés pour une large gamme d'éprouvettes en poutres reprises dans la littérature; on trouve qu'il existe une bonne concordance avec les résultats expérimentaux et théoriques. On démontre qu'en suivant une procédure simple on peut obtenir avec une bonne précision la courbe reliant les portions avant et après le pic de sollicitation en fonction du COD de la fissure. En outre, on développe une procédure simple pour calculer la charge maximum. Les valeurs prédites sont en bon accord avec les valeurs expérimentales. Le modèle SBM proposé est très simple et est basé sur des considérations rationnelles. Il est susceptible de décrire complètement le processus de rupture depuis le début de la formation de la zone intéressée jusqu'à l'amorçage de la rupture instable.

Dynamic growth of tensile cracks by ductile and brittle fracture mechanisms in a viscoplastic material

In this paper, a finite element analysis of steady-state dynamic crack growth under Mode I, plane strain, small-scale yielding conditions is performed in a rate dependent plastic material characterized by the over-stress model. The main objective of the paper is to obtain theoretically the dependence of dynamic fracture toughness on crack speed. Crack propagation due to a ductile (micro-void) mechanism or a brittle (cleavage) mechanism, as well as transition from one mode to another are considered. The conversion from ductile to brittle has been observed experimentally but has received very little attention using analytical methods. Local fracture criteria based on strains and stresses are used to describe ductile and brittle fracture mechanisms. The results obtained in this paper are in general agreement with micro-structural observations of mode conversion during fracture initiation. Finally, the particular roles played by material rate sensitivity and inertia are examined in some detail.