An adaptive singular ES-FEM for mechanics problems with singular field of arbitrary order

This paper presents a singular edge-based smoothed finite element method (sES-FEM) for mechanics problems with singular stress fields of arbitrary order. The sES-FEM uses a basic mesh of three-noded linear triangular (T3) elements and a special layer of five-noded singular triangular elements (sT5) connected to the singular-point of the stress field. The sT5 element has an additional node on each of the two edges connected to the singular-point. It allows us to represent simple and efficient enrichment with desired terms for the displacement field near the singular-point with the satisfaction of partition-of-unity property. The stiffness matrix of the discretized system is then obtained using the assumed displacement values (not the derivatives) over smoothing domains associated with the edges of elements. An adaptive procedure for the sES-FEM is proposed to enhance the quality of the solution with minimized number of nodes. Several numerical examples are provided to validate the reliability of the present sES-FEM method. (C) 2012 Elsevier B.V. All rights reserved.

Ductile fracture by multiple void growth and interaction ahead of a notch tip in polycrystalline plastic solids

The objectives of this paper are to study the effects of plastic anisotropy and evolution in crystallographic texture with deformation on the ductile fracture behaviour of polycrystalline solids. To this end, numerical simulations of multiple void growth and interaction ahead of a notch tip are performed under mode I, plane strain, small scale yielding conditions using two approaches. The first approach is based on the Hill yield theory, while the second employs crystal plasticity constitutive equations and a Taylor-type homogenization in order to represent the ductile polycrystalline solid. The initial textures pertaining to continuous cast Al-Mg AA5754 sheets in recrystallized and cold rolled conditions are considered. The former is nearly-isotropic, while the latter displays pronounced anisotropy. The results indicate distinct changes in texture in the ligaments bridging the voids ahead of the notch tip with increase in load level which gives rise to retardation in porosity evolution and increase in tearing resistance for both materials.

On the mechanism and the length scales involved in the ductile fracture of a bulk metallic glass

Some bulk metallic glasses (BMGs) exhibit high crack initiation toughness due to shear band mediated plastic flow at the crack tip and yet do not display additional resistance to crack growth due to the lack of a microstructure. Thus, at crack initiation, the fracture behavior of BMGs transits from that of ductile alloys to that of brittle ceramics. In this paper, we attempt to understand the physics behind the characteristic length from the notch root at which this transition occurs, through testing of four-point bend specimens made of a nominally ductile Zr-based BMG in three different structural states. In the as-cast state, both symmetric (mode I) and asymmetric (mixed mode) bend specimens are tested. The process of shear band mediated plastic flow followed by crack initiation at the notch root was monitored through in situ imaging. Results show that stable crack growth occurs inside a dominant shear band through a distance of, similar to 60 mu m, irrespective of the structural state and mode mixity, before attaining criticality. Detailed finite element simulations show that this length corresponds to the distance from the notch root over which a positive hydrostatic stress gradient prevails. The mean ridge heights on fractured surfaces are found to correlate with the toughness of the BMG. The Argon and Salama model, which is based on the meniscus instability phenomenon at the notch root, is modified to explain the experimentally observed physics of fracture in ductile BMGs. (C) 2012 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Fracture processes of self compacting Concrete using image analysis

The mode I fracture toughness of concrete can be experimentally determined using three point bend beam in conjunction with digital image correlation (DIC). Three different geometrically similar sizes of beams are cast for this study. To study the influence of fly ash and silica fume on fracture toughness of SCC, three SCC mixes are prepared with and without mineral additions. The scanning electron microscope (SEM) images are taken on the fractured surface to add information on fracture process in SCC. From this study, it is concluded that the fracture toughness of SCC with mineral addition is higher when compared to those without mineral addition.

The Fracture Characteristics of a Near Eutectic Al-Si Based Alloy Under Compression

The fracture of eutectic Si particles dictates the fracture characteristics of Al-Si based cast alloys. The morphology of these particles is found to play an important role in fracture initiation. In the current study, the effects of strain rate, temperature, strain, and heat treatment on Si particle fracture under compression were investigated. Strain rates ranging from 3 x 10(-4)/s to 10(2)/s and three temperatures RT, 373 K, and 473 K (100 A degrees C and 200 A degrees C) are considered in this study. It is found that the Si particle fracture shows a small increase with increase in strain rate and decreases with increase in temperature at 10 pct strain. The flow stress at 10 pct strain exhibits the trend similar to particle fracture with strain rate and temperature. Particle fracture also increases with increase in strain. Large and elongated particles show a greater tendency for cracking. Most fracture occurs on particles oriented nearly perpendicular to the loading axis, and the cracks are found to occur almost parallel to the loading axis. At any strain rate, temperature, and strain, the Si particle fracture is greater for the heat-treated condition than for the non-heat-treated condition because of higher flow stress in the heat-treated condition. In addition to Si particle fracture, elongated Fe-rich intermetallic particles are also seen to fracture. These particles have specific crystallographic orientations and fracture along their major axis with the cleavage planes for their fracture being (100). Fracture of these particles might also play a role in the overall fracture behavior of this alloy since these particles cleave along their major axis leading to cracks longer than 200 mu m.

Thermal vibration analysis of monolayer graphene embedded in elastic medium based on nonlocal continuum mechanics

This paper presents the thermal vibration analysis of single-layer graphene sheet embedded in polymer elastic medium, using the plate theory and nonlocal continuum mechanics for small scale effects. The graphene is modeled based on continuum plate theory and axial stress caused by the thermal effects is also considered. Nonlocal governing equation of motion for this graphene sheet system is derived from the principle of virtual displacements. The closed form solution for thermal-vibration frequencies of a simply supported rectangular nanoplate has been obtained by using the Navier's method of solution. Numerical results obtained by the present theory are compared with available solutions in the literature and the molecular dynamics results. The influences of the small scale coefficient, the room or low temperature, the high temperature, the half wave number and the aspect ratio of nanoplate on the natural frequencies are considered and discussed in detail. The thermal vibration analysis of single- and double-layer graphene sheets are considered for the analysis. The mode shapes of the respective graphene system are also captured in this work. The present analysis results can be used for the design of the next generation of nanodevices that make use of the thermal vibration properties of the graphene.

Determination of size-independent specific fracture energy of concrete mixes by two methods

The RILEM work-of-fracture method for measuring the specific fracture energy of concrete from notched three-point bend specimens is still the most common method used throughout the world, despite the fact that the specific fracture energy so measured is known to vary with the size and shape of the test specimen. The reasons for this variation have also been known for nearly two decades, and two methods have been proposed in the literature to correct the measured size-dependent specific fracture energy (G(f)) in order to obtain a size-independent value (G(F)). It has also been proved recently, on the basis of a limited set of results on a single concrete mix with a compressive strength of 37 MPa, that when the size-dependent G(f) measured by the RILEM method is corrected following either of these two methods, the resulting specific fracture energy G(F) is very nearly the same and independent of the size of the specimen. In this paper, we will provide further evidence in support of this important conclusion using extensive independent test results of three different concrete mixes ranging in compressive strength from 57 to 122 MPa. (c) 2013 Elsevier Ltd. All rights reserved.

Experimental Studies on the Mechanics of Cohesive Frictional Granular Media

The mechanical behaviour of cohesive-frictional granular materials is a combination of the strength pervading as intergranular friction (represented as an angle of internal friction - Phi), and the cohesion (C) between these particles. Most behavioral or constitutive models of this class of granular materials comprise of a cohesion and frictional component with no regard to the length scale i.e. from the micro structural models through the continuum models. An experimental study has been made on a model granular material, viz. angular sand with different weights of binding agents (varying degrees of cohesion) at multiple length scales to physically map this phenomenon. Cylindrical specimen of various diameters - 10, 20, 38, 100, 150 mm (and with an aspect ratio of 2) are reconstituted with 2, 4 and 8% by weight of a binding agent. The magnitude of this cohesion is analyzed using uniaxial compression tests and it is assumed to correspond to the peak in the normalized stress-strain plot. Increase in the cohesive strength of the material is seen with increasing size of the specimen. A possibility of ``entanglement'' occurring in larger specimens is proposed as a possible reason for deviation from a continuum framework.

Dynamics and fracture of ligaments from a droplet on a vibrating surface

A droplet residing on a vibrating surface and in the pressure antinode of an asymmetric standing wave can spread radially outward and atomize. In this work, proper orthogonal decomposition through high speed imaging is shown to predict the likelihood of atomization for various viscous fluids based on prior information in the droplet spreading phase. Capillary instabilities are seen to affect ligament rupture. Viscous dissipation plays an important role in determining the wavelength of the most unstable mode during the inception phase of the ligaments. However, the highest ligament capillary number achieved was less than 1, and the influence of viscosity in the ligament growth and breakup phases is quite minimal. It is inferred from the data that the growth of a typical ligament is governed by a balance between the inertial force obtained from the inception phase and capillary forces. By including the effect of acoustic pressure field around the droplet, the dynamics of the ligament growth phase is revealed and the ligament growth profiles for different fluids are shown to collapse on a straight line using a new characteristic time scale.

Bilinear tension softening diagrams of concrete mixes corresponding to their size-independent specific fracture energy

The study of the fracture behaviour of concrete structures using the fictitious crack model requires two fracture properties of the concrete mix, namely, the size-independent specific fracture energy G(F). and the corresponding tension softening relation sigma(w) between the residual stress carrying capacity sigma and the crack opening w in the fracture process zone ahead of a real crack. In this paper, bi-linear tension softening diagrams of three different concrete mixes, ranging in compressive strength from 57 to 122 MPa whose size-independent specific fracture energy has been previously determined, have been constructed in an inverse manner based on the concept of a non-linear hinge from the load-crack mouth opening plots of notched three-point bend beams. (C) 2013 Elsevier Ltd. All rights reserved.

A systematic investigation of fracture mechanisms in Al-Si based eutectic alloy-Effect of Si modification

The fracture characteristics of Al-Si based eutectic alloy are investigated in the unmodified and modified conditions under compression. The investigations are carried out at different strain rates and temperatures. Fracture of the alloy starts with eutectic Si particle fracture and modification plays an important role in particle fracture. The fraction of fractured particles is found to be always lesser in the modified condition than in the unmodified condition. Particle fracture increases with increase in strain. It is found that the Si particle fracture shows an increase with increase in strain rate and decreases with increase in temperature at 10% strain. Large and elongated particles show a greater tendency for fracture in the unmodified and modified conditions. Particle orientation plays an important role on fracture and the cracks are found to occur almost in a direction normal to the tensile strain imposed upon the particles by the deforming matrix in the unmodified alloy. The modified alloy shows a random distribution of fractured particles and crack orientation. The criteria of fracture based on dislocation pile-up mechanism and fiber loading explain the observed difference in particle fracture characteristics due to modification. The particle fracture for the modified alloy is also discussed in terms of Weibull statistics and the existing models of dispersion hardening. Particle/matrix interface decohesion is observed at higher strain rates and temperatures in the modified alloy. Dendritic rotation of 10 degrees is also observed at higher strain rates, which can increase the amount of particle fracture. (C) 2013 Elsevier B.V. All rights reserved.

Incorporation of Interfacial Intermetallic Morphology in Fracture Mechanism Map for Sn-Ag-Cu Solder Joints

A fracture mechanism map (FMM) is a powerful tool which correlates the fracture behavior of a material to its microstructural characteristics in an explicit and convenient way. In the FMM for solder joints, an effective thickness of the interfacial intermetallic compound (IMC) layer (t (eff)) and the solder yield strength (sigma (ys,eff)) are used as abscissa and ordinate axes, respectively, as these two predominantly affect the fracture behavior of solder joints. Earlier, a definition of t (eff), based on the uniform thickness of IMC (t (u)) and the average height of the IMC scallops (t (s)), was proposed and shown to aptly explain the fracture behavior of solder joints on Cu. This paper presents a more general definition of t (eff) that is more widely applicable to a range of metallizations, including Cu and electroless nickel immersion gold (ENIG). Using this new definition of t (eff), mode I FMM for SAC387/Cu joints has been updated and its validity was confirmed. A preliminary FMM for SAC387/Cu joints with ENIG metallization is also presented.

Experimental study of fracture behavior of magnesium single crystals

In this work, the fracture behavior of magnesium single crystals is studied by conducting experiments with notched three point bend specimens of three crystallographic orientations. In the first and second orientations, the c-axis is along the normal to the flat surface of the notch, while in the third it is aligned with the notch front. For all the orientations, in situ electron back scattered diffraction observations made around the notch root show profuse tensile twinning of {10 (1) over bar2} type. Further, in the first two orientations basal and prismatic slip traces are identified from optical metallography. The width of the most prominent twin saturates at around 120-150 mu m, while twins continue to nucleate farther away to accommodate plastic deformation. In all the orientations, crack initiation occurs before the attainment of peak load and the crack grows stably along twin-matrix interface before deflecting at twin-twin intersections. Results show that profuse tensile twinning is an important energy dissipating mechanism that enhances the fracture toughness. (C) 2013 Elsevier B.V. All rights reserved.

Closed-form nonlinear sectional analysis of pretwisted anisotropic beam with multiple delaminations

Fiber reinforced laminated composite open-section beams are widely used as bearingless rotor flex beams because of their high specific strength and stiffness as well as fatigue life. These laminated composite structures exhibit a number of different failure modes, including fiber-matrix debonding within individual layers, delamination or separation of the layers, transverse cracks through one or more layers and fiber fracture. Delamination is a predominant failure mode in continuous fiber reinforced laminated composites and often initiate near the free edges of the structure. The appearance of delaminations in the composite rotorcraft flexbeams can lead to deterioration of the mechanical properties and, in turn, the helicopter performance as well as safety. Understanding and predicting the influence of free-edge delamination on the overall behavior of the laminates will provide quantitative measures of the extent of the damage and help ensure their damage tolerance.