New twelve node serendipity quadrilateral plate bending element based on Mindlin-Reissner theory using Integrated Force Method

The Integrated Force Method (IFM) is a novel matrix formulation developed for analyzing the civil, mechanical and aerospace engineering structures. In this method all independent/internal forces are treated as unknown variables which are calculated by simultaneously imposing equations of equilibrium and compatibility conditions. This paper presents a new 12-node serendipity quadrilateral plate bending element MQP12 for the analysis of thin and thick plate problems using IFM. The Mindlin-Reissner plate theory has been employed in the formulation which accounts the effect of shear deformation. The performance of this new element with respect to accuracy and convergence is studied by analyzing many standard benchmark plate bending problems. The results of the new element MQP12 are compared with those of displacement-based 12-node plate bending elements available in the literature. The results are also compared with exact solutions. The new element MQP12 is free from shear locking and performs excellent for both thin and moderately thick plate bending situations.

Fatigue crack propagation model for plain concrete - An analogy with population growth

In this work, an analytical model is proposed for fatigue crack propagation in plain concrete based on population growth exponential law and in conjunction with principles of dimensional analysis and self-similarity. This model takes into account parameters such as loading history, fracture toughness, crack length, loading ratio and structural size. The predicted results are compared with experimental crack growth data for constant and variable amplitude loading and are found to capture the size effect apart from showing a good agreement. Using this model, a sensitivity analysis is carried out to study the effect of various parameters that influence fatigue failure. (C) 2010 Elsevier Ltd. All rights reserved.

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.

Effects of temper level on the dependence of fatigue crack growth threshold and crack closure on the prior austenitic grain size

An attempt has been made to systematically investigate the effects of microstructural parameters, such as the prior austenite grain size (PAGS), in influencing the resistance to fatigue crack growth (FCG) in the near-threshold region under three different temper levels in a quenched and tempered high-strength steel. By austenitizing at various temperatures, the PAGS was varied from about 0.7 to 96 μm. The microstructures with these grain sizes were tempered at 200 °C, 400 °C, and 530 °C and tested for fatigue thresholds and crack closure. It has been found that, in general, three different trends in the dependence of both the total threshold stress intensity range, ΔK th , and the intrinsic threshold stress intensity range, ΔK eff, th , on the PAGS are observable. By considering in detail the factors such as cyclic stress-strain behavior, environmental effects on FCG, and embrittlement during tempering, the present observations could be rationalized. The strong dependence of ΔK th and ΔK eff, th on PAGS in microstructures tempered at 530 °C has been primarily attributed to cyclic softening and thereby the strong interaction of the crack tip deformation field with the grain boundary. On the other hand, a less strong dependence of ΔK th and ΔK eff, th on PAGS is suggested to be caused by the cyclic hardening behavior of lightly tempered microstructures occurring in 200 °C temper. In both microstructures, crack closure influenced near-threshold FCG (NTFCG) to a significant extent, and its magnitude was large at large grain sizes. Microstructures tempered at the intermediate temperatures failed to show a systematic variation of ΔKth and ΔKeff, th with PAGS. The mechanisms of intergranular fracture vary between grain sizes in this temper. A transition from “microstructure-sensitive” to “microstructure-insensitive” crack growth has been found to occur when the zone of cyclic deformation at the crack tip becomes more or less equal to PAGS. Detailed observations on fracture morphology and crack paths corroborate the grain size effects on fatigue thresholds and crack closure.

Deformational Behavior of Fiber-Reinforced Concrete Beams in Bending

An experimental study aimed at understanding the deformational behavior of conventionally reinforced steel fiber concrete beams in pure bending is reported in this paper. One group of beams has steel fibers dispersed in the entire volume of the beam and the second has fibers dispersed over half the depth of the beam on the tension side. A comparative study of the deformational characteristics of these beams has been made. Half-depth fiber inclusion, requiring only half the quantity of fibers of full-depth inclusion, is found to be equally effective in improving the deformational behavior of beams. Thus, by such modes of inclusion of fibers, an economical and efficient use of expensive steel fibers can be realized.

Monitoring Fatigue Crack Propagation in Compact Tension Specimens via Remote Sensing of Back Face Strain

Back face strain (BFS) measurement is now well-established as an indirect technique to monitor crack length in compact tension (CT) fracture specimens [1,2]. Previous work [2] developed empirical relations between fatigue crack propagation (FCP) parameters. BFS, and number of cycles for CT specimens subjected to constant amplitude fatigue loading. These predictions are experimentally validated in terms of the variations of mean values of BFS and load as a function of crack length. Another issue raised by this study concerns the validity of assigning fixed values for the Paris parameters C and n to describe FCP in realistic materials.

Kinetic approach to dislocation bending in low-mobility films

A kinetic model has been developed for dislocation bending at the growth surface in compressively stressed low-mobility films such as III-V nitrides. It is based on a reduction in the number of atoms at the growth surface. Stress and nonstress sources of driving force for such a reduction are discussed. A comparison between the derived equations and experimentally measured stress evolution data yields good agreement between the predicted and observed angles through which dislocations bend.

Influence of crack tip constraint on void growth in ductile FCC single crystals

Effect of constraint (stress triaxiality) on void growth near a notch tip in a FCC single crystal is investigated. Finite element simulations within the modified boundary layer framework are conducted using crystal plasticity constitutive equations and neglecting elastic anisotropy. Displacement boundary conditions based on model, elastic, two term K-T field are applied on the outer boundary of a large circular domain. A pre-nucleated void is considered ahead of a stationary notch tip. The interaction between the notch tip and the void is studied under different constraints (T-stress levels) and crystal orientations. It is found that negative T-stress retards the mechanisms of ductile fracture. However, the extent of retardation depends on the crystal orientation. Further, it is found that there exists a particular orientation which delays the ductile fracture processes and hence can potentially improve ductility. This optimal orientation depends on the constraint level. (C) 2010 Published by Elsevier B.V.

Fracture Properties of Concrete-Concrete Interfaces Using Digital Image Correlation

The mode I and mode II fracture toughness and the critical strain energy release rate for different concrete-concrete jointed interfaces are experimentally determined using the Digital Image Correlation technique. Concrete beams having different compressive strength materials on either side of a centrally placed vertical interface are prepared and tested under three-point bending in a closed loop servo-controlled testing machine under crack mouth opening displacement control. Digital images are captured before loading (undeformed state) and at different instances of loading. These images are analyzed using correlation techniques to compute the surface displacements, strain components, crack opening and sliding displacements, load-point displacement, crack length and crack tip location. It is seen that the CMOD and vertical load-point displacement computed using DIC analysis matches well with those measured experimentally.

A general procedure for modified crack closure integral in 3D problems with cracks

In linear elastic fracture mechanics (LEFM), Irwin's crack closure integral (CCI) is one of the signficant concepts for the estimation of strain energy release rates (SERR) G, in individual as well as mixed-mode configurations. For effective utilization of this concept in conjunction with the finite element method (FEM), Rybicki and Kanninen [Engng Fracture Mech. 9, 931 938 (1977)] have proposed simple and direct estimations of the CCI in terms of nodal forces and displacements in the elements forming the crack tip from a single finite element analysis instead of the conventional two configuration analyses. These modified CCI (MCCI) expressions are basically element dependent. A systematic derivation of these expressions using element stress and displacement distributions is required. In the present work, a general procedure is given for the derivation of MCCI expressions in 3D problems with cracks. Further, a concept of sub-area integration is proposed which facilitates evaluation of SERR at a large number of points along the crack front without refining the finite element mesh. Numerical data are presented for two standard problems, a thick centre-cracked tension specimen and a semi-elliptical surface crack in a thick slab. Estimates for the stress intensity factor based on MCCI expressions corresponding to eight-noded brick elements are obtained and compared with available results in the literature.

Modified Crack Closure Integral (MCCI) For 3-D Problems Using 20-Noded Brick Elements

The Modified Crack Closure Integral (MCCI) technique based on Irwin's crack closure integral concept is very effective for estimation of strain energy release rates G in individual as well as mixed-mode configurations in linear elastic fracture mechanics problems. In a finite element approach, MCCI can be evaluated in the post-processing stage in terms of nodal forces and displacements near the crack tip. The MCCI expressions are however, element dependent and require a systematic derivation using stress and displacement distributions in the crack tip elements. Earlier a general procedure was proposed by the present authors for the derivation of MCCI expressions for 3-dimensional (3-d) crack problems modelled with 8-noded brick elements. A concept of sub-area integration was proposed to estimate strain energy release rates at a large number of points along the crack front. In the present paper a similar procedure is adopted for the derivation of MCCI expressions for 3-d cracks modelled with 20-noded brick elements. Numerical results are presented for centre crack tension and edge crack shear specimens in thick slabs, showing a comparison between present results and those available in the literature.

Modified crack closure integral (MCCI) for 3-d problems using 20 noded brick elememnts

The Modified Crack Closure Integral (MCCI) technique based on Irwin's crack closure integral concept is very effective for estimation of strain energy release rates G in individual as well as mixed-mode configurations in linear elastic fracture mechanics problems. In a finite element approach, MCCI can be evaluated in the post-processing stage in terms of nodal forces and displacements near the crack tip. The MCCI expressions are however, element dependent and require a systematic derivation using stress and displacement distributions in the crack tip elements. Earlier a general procedure was proposed by the present authors for the derivation of MCCI expressions for 3-dimensional (3-d) crack problems modelled with 8-noded brick elements. A concept of sub-area integration was proposed to estimate strain energy release rates at a large number of points along the crack front. In the present paper a similar procedure is adopted for the derivation of MCCI expressions for 3-d cracks modelled with 20-noded brick elements. Numerical results are presented for centre crack tension and edge crack shear specimens in thick slabs, showing a comparison between present results and those available in the literature.

Fatigue crack propagation model and size effect in concrete using dimensional analysis

It is well known that fatigue in concrete causes excessive deformations and cracking leading to structural failures. Due to quasi-brittle nature of concrete and formation of a fracture process zone, the rate of fatigue crack growth depends on a number of parameters, such as, the tensile strength, fracture toughness, loading ratio and most importantly the structural size. In this work, an analytical model is proposed for estimating the fatigue crack growth in concrete by using the concepts of dimensional analysis and including the above parameters. Knowing the governed and the governing parameters of the physical problem and by using the concepts of self-similarity, a relationship is obtained between different parameters involved. It is shown that the proposed fatigue law is able to capture the size effect in plain concrete and agrees well with different experimental results. Through a sensitivity analysis, it is shown that the structural size plays a dominant role followed by loading ratio and the initial crack length in fatigue crack propagation. (C) 2010 Elsevier Ltd. All rights reserved.

Fracture Analysis of Concrete Structural Components Accounting for Tension Softening Effect

This paper presents methodologies for fracture analysis of concrete structural components with and without considering tension softening effect. Stress intensity factor (SIF) is computed by using analytical approach and finite element analysis. In the analytical approach, SW accounting for tension softening effect has been obtained as the difference of SIP obtained using linear elastic fracture mechanics (LEFM) principles and SIP due to closing pressure. Superposition principle has been used by accounting for non-linearity in incremental form. SW due to crack closing force applied on the effective crack face inside the process zone has been computed using Green's function approach. In finite element analysis, the domain integral method has been used for computation of SIR The domain integral method is used to calculate the strain energy release rate and SIF when a crack grows. Numerical studies have been conducted on notched 3-point bending concrete specimen with and without considering the cohesive stresses. It is observed from the studies that SW obtained from the finite element analysis with and without considering the cohesive stresses is in good agreement with the corresponding analytical value. The effect of cohesive stress on SW decreases with increase of crack length. Further, studies have been conducted on geometrically similar structures and observed that (i) the effect of cohesive stress on SW is significant with increase of load for a particular crack length and (iii) SW values decreases with increase of tensile strength for a particular crack length and load.