Influence of stochastic mesoscopic structure on macroscopic mechanical behavior of brittle material

A newly developed numerical code, MFPA(2D) (Material Failure Process Analysis), is applied to study the influence of stochastic mesoscopic structure on macroscopic mechanical behavior of rock-like materials. A set of uniaxial compression tests has been numerically studied with numerical specimens containing pre-existing crack-like flaw. The numerical results reveal the influence of random mesoscopic structure on failure process of brittle material, which indicates that the variation of failure mode is strongly sensitive to the local disorder feature of the specimen. And the patterns of the crack evolution in the specimens are very different from each other due to the random mesoscopic structure in material. The results give a good explanation for various kinds of fracture modes and peak strength variation observed in laboratory studies with specimens made from the same rock block being statistically homogenous in macro scale. In addition, the evolution of crack is more complicated in heterogeneous cases than in homogeneous cases.

Evolution induced catastrophe in a nonlinear dynamical model of material failure

In order to study the failure of disordered materials, the ensemble evolution of a nonlinear chain model was examined by using a stochastic slice sampling method. The following results were obtained. (1) Sample-specific behavior, i.e. evolutions are different from sample to sample in some cases under the same macroscopic conditions, is observed for various load-sharing rules except in the globally mean field theory. The evolution according to the cluster load-sharing rule, which reflects the interaction between broken clusters, cannot be predicted by a simple criterion from the initial damage pattern and even then is most complicated. (2) A binary failure probability, its transitional region, where globally stable (GS) modes and evolution-induced catastrophic (EIC) modes coexist, and the corresponding scaling laws are fundamental to the failure. There is a sensitive zone in the vicinity of the boundary between the GS and EIC regions in phase space, where a slight stochastic increment in damage can trigger a radical transition from GS to EIC. (3) The distribution of strength is obtained from the binary failure probability. This, like sample-specificity, originates from a trans-scale sensitivity linking meso-scopic and macroscopic phenomena. (4) Strong fluctuations in stress distribution different from that of GS modes may be assumed as a precursor of evolution-induced catastrophe (EIC).

Micro- and macroscopic characterizations of the viscoelastic fracture of resin-based fibre composites

Micro- and macroscopic characterizations of the viscoelastic fracture of a unidirectional carbon-fibre-reinforced epoxy composite are presented. First, the micro-cracking behavior of the material is studied by the use of scanning electron microscopy; the in situ creep cracking process is observed and the crack propagation is measured. In order to obtain insight into the mechanisms of the observed creep cracking, macroscopic investigations were also carried out. Finite-element method simulations were carried out to calculate the stress distribution and the variation of stresses with time. A theoretical analysis of the orthotropy of viscoelastic fracture behavior of the material is also conducted.

Influence of secondary void damage in the matrix material around voids

The mechanism of ductile damage caused by secondary void damage in the matrix around primary voids is studied by large strain, finite element analysis. A cylinder embedding an initially spherical void, a plane stress cell with a circular void and plane strain cell with a cylindrical or a flat void are analysed under different loading conditions. Secondary voids of smaller scale size nucleate in the strain hardening matrix, according to the requirements of some stress/strain criteria. Their growth and coalescence, handled by the empty element technique, demonstrate distinct mechanisms of damage as circumstances change. The macroscopic stress-strain curves are decomposed and illustrated in the form of the deviatoric and the volumetric parts. Concerning the stress response and the void growth prediction, comparisons are made between the present numerical results and those of previous authors. It is shown that loading condition, void growth history and void shape effect incorporated with the interaction between two generations of voids should be accounted for besides the void volume fraction.

Material response and failure mechanism of unidirectional metal matrix composites under impulsive shear loading

The material response and failure mechanism of unidirectional metal matrix composite under impulsive shear loading are investigated in this paper. Both experimental and analytical studies were performed. The shear strength of unidirectional C-f/A356.0 composite and A356.0 aluminum alloy at high strain rate were measured with a modified split Hopkinson torsional bar technique. The results indicated that the carbon fibers did not improve the shear strength of aluminum matrix if the fiber orientation aligned with the shear loading axis. The microscopic inspection of the fractured surface showed a multi-scale zigzag feature which implied a complicated shear failure mechanism in the composite. In addition to testing, the micromechanical stress field in the composite was analyzed by the generalized Eshelby equivalent method (GEEM). The influence of cracking in matrix on the micromechanical stress field was investigated as well. The results showed that the stress distribution in the composite is quite nonhomogeneous and very high shear stress concentrations are found in some regions in the matrix. The high shear stress concentration in the matrix induces tensile cracking at 45 degrees to the shear direction. This in turn aggravates the stress concentration at the fiber/matrix interface and finally leads to a catastrophic failure in the composite. From the correlation between the analysis and experimental results, the shear failure mechanism of unidirectional C-f/A356.0 composite can be elucidated qualitatively.

Experimental Studies of the Material Properties of the Forewing of Cicada (Homoptera, Cicadidae)

Detailed investigations on the structural and mechanical properties of the forewing of the cicada were carried out. Measurement of the structures of the wings showed that the thickness of the membrane of each cell and the diameter of each vein were non-uniform in both the longitudinal and transverse directions, and their means were approximately 12.2 and 133.3 mum, respectively. However, the aspect ratios of the wings and the bodies were quite uniform and were approximately equal to 2.98 and 2.13, respectively. Based on the measured thickness, mass and area of the membranes of the cells, the mean density and the mean area density of the wing were approximately 2.3 g cm(-3) and 2.8 x 10(-3) g cm(-2), respectively. In addition, the diameters of the veins of the wings, including the diameters of the holes in the vein of the leading edge, were examined. The mechanical properties of the wing were investigated separately by nanoindentation and tensile testing. The results indicated that the mean Young's modulus, hardness and yield stress of the membranes of the wings were approximately 3.7 Gpa, 0.2 Gpa and 29 Mpa, respectively, and the mean Young's modulus and strength of the veins along the direction of the venation of wings were approximately 1.9 Gpa and 52 Mpa, respectively. Finally, the relevant results were briefly analyzed and discussed, providing a guideline to the biomimetic design of the aerofoil materials of micro air vehicles.

Elastic-plastic constitutive relation of particle reinforced composites

In this paper, a systematic approach is proposed to obtain the macroscopic elastic-plastic constitutive relation of particle reinforced composites (PRC). The strain energy density of PRC is analyzed based on the cell model, and Che analytical formula for the macro-constitutive relation of PRC is obtained. The strength effects of volume fraction of the particle and the strain hardening exponent of matrix material on the macro-constitutive relation are investigated, the relation curve of strain versus stress of PRC is calculated in detail. The present results are consistent; with the results given in the existing references.

Determination of Proper Processing Conditions of Material Surface Heat Treatment Through Finite Element Method

A two-dimensional model has been developed based on the experimental results of stainless steel remelting with the laminar plasma technology to investigate the transient thermo-physical characteristics of the melt pool liquids. The influence of the temperature field, temperature gradient, solidification rate and cooling rate on the processing conditions has been investigated numerically. Not only have the appropriate processing conditions been determined according to the calculations, but also they have been predicted with a criterion established based on the concept of equivalent temperature area density (ETAD) that is actually a function of the processing parameters and material properties. The comparison between the resulting conditions shows that the ETAD method can better predict the optimum condition.

Influences of strain rate and temperature variation on material intrinsic length of plasticity strain gradient theory

Cowper-Symonds and Johnson-Cook dynamic constitutive relations are used to study the influence of both strain rate effect and temperature variation on the material intrinsic length scale in strain gradient plasticity. The material intrinsic length scale decreases with increasing strain rates, and this length scale increases with temperature.

Three-Dimensional Modeling of Heat Transfer and Fluid Flow in Laminar-Plasma Material Re-Melting Processing

Modeling study is performed concerning the heat transfer and fluid flow for a laminar argon plasma jet impinging normally upon a flat workpiece exposed to the ambient air. The diffusion of the air into the plasma jet is handled by using the combined-diffusion-coefficient approach. The heat flux density and jet shear stress distributions at the workpiece surface obtained from the plasma jet modeling are then used to study the re-melting process of a carbon steel workpiece. Besides the heat conduction within the workpiece, the effects of the plasma-jet inlet parameters (temperature and velocity), workpiece moving speed, Marangoni convection, natural convection etc. on the re-melting process are considered. The modeling results demonstrate that the shapes and sizes of the molten pool in the workpiece are influenced appreciably by the plasma-jet inlet parameters, workpiece moving speed and Marangoni convection. The jet shear stress manifests its effect at higher plasma-jet inlet velocities, while the natural convection effect can be ignored. The modeling results of the molten pool sizes agree reasonably with available experimental data.

Implication of scaling hierarchy associated with nonequilibrium: field and particulate

The advent of nanotechnology has necessitated a better understanding of how material microstructure changes at the atomic level would affect the macroscopic properties that control the performance. Such a challenge has uncovered many phenomena that were not previously understood and taken for granted. Among them are the basic foundation of dislocation theories which are now known to be inadequate. Simplifying assumptions invoked at the macroscale may not be applicable at the micro- and/or nanoscale. There are implications of scaling hierrachy associated with in-homegeneity and nonequilibrium. of physical systems. What is taken to be homogeneous and equilibrium at the macroscale may not be so when the physical size of the material is reduced to microns. These fundamental issues cannot be dispensed at will for the sake of convenience because they could alter the outcome of predictions. Even more unsatisfying is the lack of consistency in modeling physical systems. This could translate to the inability for identifying the relevant manufacturing parameters and rendering the end product unpractical because of high cost. Advanced composite and ceramic materials are cases in point. Discussed are potential pitfalls for applying models at both the atomic and continuum levels. No encouragement is made to unravel the truth of nature. Let it be partiuclates, a smooth continuum or a combination of both. The present trend of development in scaling tends to seek for different characteristic lengths of material microstructures with or without the influence of time effects. Much will be learned from atomistic simulation models to show how results could differ as boundary conditions and scales are changed. Quantum mechanics, continuum and cosmological models provide evidence that no general approach is in sight. Of immediate interest is perhaps the establishment of greater precision in terminology so as to better communicate results involving multiscale physical events.

The effect of nucleation time on the growth of a microvoid in a viscoelastic material

In this paper, discussions are focused on the growth of a nucleated void in a viscoelastic material. The in situ tensile tests of specimens made of high-density polyethylene, filled with spherical glass beads (HDPE/GB) are carried out under SEM. The experimental result indicates that the microvoid nucleation is induced by the partially interfacial debonding of particles. By means of the Laplace transform and the Eshelby's equivalent inclusion method, a new analytical expression of the void strain at different nucleation times is derived. It can be seen that the strain of the nucleated void depends not only on the remote strain history, but also on the nucleation time. This expression is also illustrated by numerical examples, and is found to be of great usefulness in the study of damage evolution in viscoelastic materials.

Weibull modulus for diverse strength due to sample-specificity

By sample specificity it is meant that specimens with the same nominal material parameters and tested under the same environmental conditions may exhibit different behavior with diversified strength. Such an effect has been widely observed in the testing of material failure and is usually attributed to the heterogeneity of material at the mesoscopic level. The degree with which mesoscopic heterogeneity affects macroscopic failure is still not clear. Recently, the problem has been examined by making use of statistical ensemble evolution of dynamical system and the mesoscopic stress re-distribution model (SRD). Sample specificity was observed for non-global mean stress field models, such as the duster mean field model, stress concentration at tip of microdamage, etc. Certain heterogeneity of microdamage could be sensitive to particular SRD leading to domino type of coalescence. Such an effect could start from the microdamage heterogeneity and then be magnified to other scale levels. This trans-scale sensitivity is the origin of sample specificity. The sample specificity leads to a failure probability Phi (N) with a transitional region 0 < (N) < 1, so that globally stable and catastrophic modes could co-exist. It is found that the scatter in strength can fit the Weibull distribution very well. Hence, the Weibull modulus is indicative of sample specificity. Numerical results obtained from the SRD for different non-global mean stress fields show that Weibull modulus increases with increasing sample span and influence region of microdamage.

Dynamic stress intensity factor of a functionally graded material under antiplane shear loading

The dynamic response of a finite crack in an unbounded Functionally Graded Material (FGM) subjected to an antiplane shear loading is studied in this paper. The variation of the shear modulus of the functionally graded material is modeled by a quadratic increase along the direction perpendicular to the crack surface. The dynamic stress intensity factor is extracted from the asymptotic expansion of the stresses around the crack tip in the Laplace transform plane and obtained in the time domain by a numerical Laplace inversion technique. The influence of graded material property on the dynamic intensity factor is investigated. It is observed that the magnitude of dynamic stress intensity factor for a finite crack in such a functionally graded material is less than in the homogeneous material with a property identical to that of the FGM crack plane.

Weak Zone Related Seismic Cycles in Progressive Failure Leading to Collapse in Brittle Crust

Until quite recently our understanding of the basic mechanical process responsible for earthquakes and faulting was not well known. It can be argued that this was partly a consequence of the complex nature of fracture in crust and in part because evidence of brittle phenomena in the natural laboratory of the earth is often obliterated or obscured by other geological processes. While it is well understood that the spatial and temporal complexity of earthquakes and the fault structures emerge from geometrical and material built-in heterogeneities, one important open question is how the shearing becomes localized into a band of intense fractures. Here the authors address these questions through a numerical approach of a tectonic plate by considering rockmass heterogeneity both in microscopic scale and in mesoscopic scale. Numerical simulations of the progressive failure leading to collapse under long-range slow driving forces in the far-field show earthquake-like rupture behavior. $En Echelon$ crack-arrays are reproduced in the numerical simulation. It is demonstrated that the underlying fracturing induced acoustic emissions (or seismic events) display self-organized criticality------from disorder to order. The seismic cycles and the geometric structures of the fracture faces, which are found greatly depending on the material heterogeneity (especially on the macroscopic scale), agree with that observed experimentally in real brittle materials. It is concluded that in order to predict a main shock, one must have extremely detailed knowledge on very minor features of the earth's crust far from the place where the earthquake originated. If correct, the model proposed here seemingly provides an explanation as to why earthquakes to date are not predicted so successfully. The reason is not that the authors do not understand earthquake mechanisms very well but that they still know little about our earth's crust.