A micromechanical analysis of elastoplastic behavior of porous materials

In this paper, we present a micromechanical analysis of elastoplastic behavior of porous materials. The non-uniform transformation field analysis (NTFA) is used and the non-uniform distribution of local plastic strain in the solid matrix is taken into account. Comparisons with the classical Gurson's model and standard FEM solution are presented.     

A micromechanical model for porous materials with a reinforced matrix

In this study a micromechanical model is proposed for ductile porous material whose matrix is reinforced by small inclusions. The solid phase is described by a pressure sensitive plastic model. Based on works of Maghous et al. [6], a macroscopic plastic criterion is firstly obtained by using a two-step homogenization procedure. The effect of porosity at the mesoscale and the influence of inclusions at the microscale are taken into account simultaneously by this criterion. With a non-associated plastic flow rule, the micro-macro model is applied to modeling of mechanical behavior of a cement paste. In particular, we have considered at the microscopic scale the formation of calcite grains by carbonation process in the solid matrix. The studied cement paste is then seen as a reinforced matrix–pore system. Comparisons between numerical results and experimental data show that the proposed model is able to capture the main features of the mechanical behavior of the studied material.     

Three-dimensional micromechanical modeling of voided polymeric materials

A three-dimensional micromechanical unit cell model for particle-filled materials is presented. The cell model is based on a Voronoi tessellation of particles arranged on a body-centered cubic (BCC) array. The three-dimensionality of the present cell model enables the study of several deformation modes, including uniaxial, plane strain and simple shear deformations, as well as arbitrary principal stress states.The unit cell model is applied to studies on the micromechanical and macromechanical behavior of rubber-toughened polycarbonate. Different load cases are examined, including plane strain deformation, simple shear deformation and principal stress states. For a constant macroscopic strain rate, the different load cases show that the macroscopic flow strength of the blend decreases with an increase in void volume fraction, as expected. The main mechanism for plastic deformation is broad shear banding across inter-particle ligaments. The distributed nature of plastic straining acts to reduce the amount of macroscopic strain softening in the blend as the initial void volume fraction is increased. In the case of plane strain deformation, the plastic flow is observed to initiate across inter-particle ligaments in the direction of constraint. This particular mode of deformation could not have been captured using a two-dimensional, plane strain idealization of cylindrical voids in a matrix.The potential for localized crazing and/or cavitation in the matrix is addressed. It is observed that the introduction of voids acts to relieve hydrostatic stress in the matrix material, compared to the homopolymer. It is also seen that the predicted peak hydrostatic stress in the matrix is higher under plane strain deformation than under triaxial tension (with equal lateral stresses), for the same macroscopic stress triaxiality.The effect of void volume fraction on the macroscopic uniaxial tension behavior of the different blends is examined using a Considère construction for dilatant materials. The natural draw ratio was predicted to decrease with an increase in void volume fraction.     

A micromechanical basis for partitioning the evolution of grain bridging in brittle materials

A micromechanical model is developed for grain bridging in monolithic ceramics. Specifically, bridge formation of a single, non-equiaxed grain spanning adjacent grains is addressed. A cohesive zone framework enables crack initiation and propagation along grain boundaries. The evolution of the bridge is investigated through a variance in both grain angle and aspect ratio. We propose that the bridging process can be partitioned into five distinct regimes of resistance: propagate, kink, arrest, stall, and bridge. Although crack propagation and kinking are well understood, crack arrest and subsequent “stall” have been largely overlooked. Resistance during the stall regime exposes large volumes of microstructure to stresses well in excess of the grain boundary strength. Bridging can occur through continued propagation or reinitiation ahead of the stalled crack tip. The driving force required to reinitiate is substantially greater than the driving force required to kink. In addition, the critical driving force to reinitiate is sensitive to grain aspect ratio but relatively insensitive to grain angle. The marked increase in crack resistance occurs prior to bridge formation and provides an interpretation for the rapidly rising resistance curves which govern the strength of many brittle materials at realistically small flaw sizes.     

Coupling between experimental measurements and polycrystal finite element calculations for micromechanical study of metallic materials

This paper presents a methodology for multiscale coupling between the morphology and texture of a microstructure as has been characterised experimentally, and the results of mechanical strain field analysis. This methodology is based on a coupling between experimental characterisation of the microstructure, in situ and/or ex situ mechanical tests, local strain field measurements performed at the grain scale, and finite element simulations. First, with orientation imaging microscopy, a map of the microstructure is generated that can be meshed. Then, finite element calculations are carried out on this mesh, using a constitutive law which takes into account the crystallographic orientation of each grain, as has been determined by the orientation imaging itself. These numerical results are then compared to the experimental strain field as obtained by digital image correlation at the scale of the grains.     

A micromechanical model for predicting thermal properties and thermo-viscoelastic responses of functionally graded materials

This study introduces a micromechanical model for predicting effective thermo-viscoelastic behaviors of a functionally graded material (FGM). The studied FGM consists of two constituents with varying compositions through the thickness. The microstructure of the FGM is idealized as solid spherical particles spatially distributed in a homogeneous matrix. The mechanical properties of each constituent can vary with temperature and time, while the thermal properties are allowed to change with temperature. The FGM model includes a transition zone where the inclusion and matrix constituents are not well defined. At the transition zone, an interchange between the two constituents as inclusion and matrix takes place such that the maximum inclusion volume contents before and after the transition zone are less than 50%. A micromechanical model is used to determine through-thickness effective thermal conductivity, coefficient of thermal expansion, and time-dependent compliance/stiffness of the FGM. The material properties at the transition zone are assumed to vary linearly between the two properties at the bounds of the transition zone. The micromechanical model is designed to be compatible with finite element (FE) scheme and used to analyze heat conduction and thermo-viscoelastic responses of FGMs. Available experimental data and analytical solutions in the literature are used to verify the thermo-mechanical properties of FGMs. The effects of time and temperature dependent constituent properties on the overall temperature, stress, and displacement fields in the FGM are also examined.     

Evolution of the micromechanical properties of impacted granular materials

The impact of a projectile on a granular material induces important changes in the micromechanical properties of the impacted material. These changes are studied using a Discrete Element Method model of the impact. The numerical results show that the impact first entails an energy propagation from the impact point to the limit of the sample through the existing force chains. A significant kinetic energy and a total breakage of the existing contact force chains are the main consequences of the energy propagation. During the long recovery balance phase observed after the energy propagation phase, frictional processes cause the kinetic energy dissipation. The motions of the particles and the numerous contact openings first prevent the formation of stable force chains. However, for long durations after the beginning of the impact, contact openings stop. The balance recovery phase therefore finally results in the creation of new stable contact force networks.     

A one dimensional analysis of convective drying of porous materials

A one-dimensional finite element analysis has been carried out to study temperature and moisture variation of porous materials during convective drying. The predicted results for brick and mortar with constant transport properties agree with available experimental data. Henderson’s modified equation has been used to obtain the desorption data for the materials.     

A study of fracture in fiber-reinforced cement-based composites using laser holographic interferometry

Micromechanisms of matrix fracture in Portland cement-based fiber composites were studied by means of reflection holographic interferometry and quantitative image analysis. An experimental investigation was conducted with different volume concentrations of polypropylene fibers. Uniaxial tensile specimens were loaded within a strain range of up to three percent. The deformation history was recorded as interference holograms. The holograms were acquired in an image-analysis system. After enhancement, they were analyzed for crack density, length, opening profile, and spacing. The evolution of microcracks, their propagation, distributed microcracking and the material response beyond the characteristic damage state are also discussed.     

Some micromechanical aspects of failure in granular materials based on second-order work

This paper discusses the notion of failure in a granular assembly by examining the key microstructural mechanisms which are most likely to trigger the nucleation and propagation of instabilities within a granular material. For this purpose, the key variable to predict the occurrence of failure, known as second-order work, is expressed from variables on the grain scale. The local behaviour incidents (where contacts may open or slide), compared to the global response of the assembly, are analysed by two approaches. First, numerical computations made by a discrete element model confirm the microscopic definition of the second-order work. Secondly, a micromechanical model, based on a homogenization procedure, relating the macroscopic behaviour to microscopic ingredients, namely contact planes, points to a close link between the occurrence of failure on the macroscopic scale as well as on the contact planes.     

An embedding method for modeling micromechanical behavior and macroscopic properties of composite materials

This paper presents a numerical method for modeling the micromechanical behavior and macroscopic properties of fiber-reinforced composites and perforated materials. The material is modeled by a finite rectangular domain containing multiple circular holes and elastic inclusions. The rectangular domain is assumed to be embedded within a larger circular domain with fictitious boundary loading represented by truncated Fourier series. The analytical solution for the complementary problem of a circular domain containing holes and inclusions is obtained by using a combination of the series expansion technique with a direct boundary integral method. The boundary conditions on the physical external boundary are satisfied by adopting an overspecification technique based on a least squares approximation. All of the integrals arising in the method can be evaluated analytically. As a result, the elastic fields and effective properties can be expressed explicitly in terms of the coefficients in the series expansions. Several numerical experiments are conducted to verify the accuracy and efficiency of the numerical method and to demonstrate its application in determination of the macroscopic properties of composite materials.     

A micromechanical model of surface steps

The surface of an epitaxial thin film typically consists of terraces separated by steps of atomic height and it evolves largely by the motion of steps. Steps are sources of stress that interact with other residual stress fields, and these interactions have a profound effect on surface evolution. A model of the elastic field arising from a two-dimensional step is presented as a departure from the commonly used half-plane point-multipole model. The field is calculated asymptotically for small step height up to second order in terms of ‘structural’ parameters that can be determined from empirical data or atomistic calculations. Effects of a lattice mismatch and surface stress are included. The model is shown to be in remarkable agreement with atomistic predictions. It is demonstrated that second-order terms are necessary for understanding non-trivial step-step interactions, and that these second-order fields cannot be described by point sources on a half-plane.     

Modeling the viscoplastic micromechanical response of two-phase materials using Fast Fourier Transforms

A viscoplastic approach using the Fast Fourier Transform (FFT) method for obtaining local mechanical response is utilized to study microstructure-property relationships in composite materials. Specifically, three-dimensional, two-phase digital materials containing isotropically coarsened particles surrounded by a matrix phase, generated through a Kinetic Monte Carlo Potts model for Ostwald ripening, are used as instantiations in order to calculate the stress and strain-rate fields under uniaxial tension. The effects of the morphology of the matrix phase, the volume fraction and the contiguity of particles, and the polycrystallinity of matrix phase, on the stress and strain-rate fields under uniaxial tension are examined. It is found that the first moments of the stress and strain-rate fields have a different dependence on the particle volume fraction and the particle contiguity from their second moments. The average stresses and average strain-rates of both phases and of the overall composite have rather simple relationships with the particle volume fraction whereas their standard deviations vary strongly, especially when the particle volume fraction is high, and the contiguity of particles has a noticeable effect on the mechanical response. It is also found that the shape of stress distribution in the BCC hard particle phase evolves as the volume fraction of particles in the composite varies, such that it agrees with the stress field in the BCC polycrystal as the volume of particles approaches unity. Finally, it is observed that the stress and strain-rate fields in the microstructures with a polycrystalline matrix are less sensitive to changes in volume fraction and contiguity of particles.     

Rate dependent properties of perovskite type tetragonal piezoelectric materials using micromechanical model

In this paper, a three-dimensional micromechanical model is presented for simulation of the rate dependent properties of certain perovskite type tetragonal piezoelectric materials. The model is based on linear constitutive, nonlinear domain switching, and linear kinetics theories. The simulation starts with a virgin bulk material of randomly oriented grains. Then the material is electrically loaded with an alternating voltage of various frequencies, which are in the order of 0.01 Hz to 1 Hz. An energy equation in combination with a probability function is used to determine the onset of the domain switching inside the grains. Such a probability function leads to a better phenomenological model for the domain switching even for electrical loadings, which are in a range far below the coercive fields. The propagation of the domain wall during the domain switching process in grains is modeled by means of linear kinetics relations after domain nucleation. The response of the bulk ceramic is predicted by averaging the response of individual grains using Euler angles for the transformation from local coordinates of the grains to global coordinate. Electric displacement hysteresis loops for different frequencies and amplitudes of the alternating electric fields are simulated. A simple micromechanical model without the probabilistic approach is compared with the one that takes it into account. Both models give important insights into the rate dependency of piezoelectric materials, which was observed in some experiments reported in the literature.     

A micromechanical model based on hypersingular integro-differential equations for analyzing micro-crazed interfaces between dissimilar elastic materials

The current work models a weak(soft) interface between two elastic materials as containing a periodic array of micro-crazes. The boundary conditions on the interfacial micro-crazes are formulated in terms of a system of hypersingular integro-differential equations with unknown functions given by the displacement jumps across opposite faces of the micro-crazes. Once the displacement jumps are obtained by approximately solving the integro-differential equations, the effective stiffness of the micr...     

A micromechanical study of residual stress and its effect on transverse failure in polymer–matrix composites

Cure residual stress and its effect on damage in unidirectional fibre-reinforced polymer–matrix composites under transverse loading were studied using a micromechanical unit cell model and the finite element method. The overall residual stress introduced from curing was determined by considering two contributions: volume shrinkage of matrix resin from the crosslink polymerization during isothermal curing and thermal contraction of both resin and fibre as a result of cooling from the curing temperature to room temperature. To examine the effect of residual stress on failure, a model based on the Maximum Principal Stress criterion and stiffness degradation technique was used for damage analysis of the unit cell subjected to mechanical loading after curing. Predicted damage initiation and evolution are clearly influenced by the inclusion of residual stress. Residual stress is always detrimental for transverse compressive loading and pure shear loading. For transverse tensile loading, residual stress is detrimental for relatively low resin strength and beneficial for relatively high resin strength. Failure envelopes were obtained for both biaxial normal loading and combined shear and normal loading and the results show that residual stress results in a shifting and contraction of the failure envelopes.     

A micromechanical theory of grain-size dependence in metal plasticity

The effect of grain-size on the elastoplastic behavior of metals is investigated from the micromechanics standpoint. First, based on the observations that dislocation pile-ups, formation of cell structures, and other inelastic activities influenced by the presence of grain boundary actually take place transcrystallinely, a grain-size dependent constitutive equation is proposed for the slip deformation of slip systems. By means of a modified Hill's self-consistent relation the local stress of a grain is calculated, and used in conjunction with this constitutive equation to evaluate the plastic strain of each constituent grain. The grain-size effect on the plastic flow of polycrystals then can be determined by an averaging process. To check the validity of the proposed theory it was finally applied to predict the stress-strain curves and flow stresses of a copper at various grain-sizes. The obtained results were found to be in good agreement with experimental data.     

A micromechanical prediction of localization in a granular material

We predict localization in a strained, random aggregate of identical elasto-frictional particles. As triaxial compression proceeds, pairs of particles deform subject to local equilibrium and anisotropy develops because of contact deletion and the dependence of contact stiffness on the average strain. The combination of the deviations of the particle displacements from the average strain and the induced anisotropy results in a relation between increments in average strain and average stress that does not possess the major symmetry. This leads to the possibility of a discontinuity in a component of the incremental strain at a predicted value of the shear strain and a predicted orientation relative to the axis of greatest compression.