Statistical and stochastic description of microvoid evolution observed in dynamic failure of ductile materials

A brief review is presented of statistical approaches on microdamage evolution. An experimental study of statistical microdamage evolution in two ductile materials under dynamic loading is carried out. The observation indicates that there are large differences in size and distribution of microvoids between these two materials. With this phenomenon in mind, kinetic equations governing the nucleation and growth of microvoids in nonlinear rate-dependent materials are combined with the balance law of void number to establish statistical differential equations that describe the evolution of microvoids' number density. The theoretical solution provides a reasonable explanation of the experimentally observed phenomenon. The effects of stochastic fluctuation which is influenced by the inhomogeneous microscopic structure of materials are subsequently examined (i.e. stochastic growth model). Based on the stochastic differential equation, a Fokker-Planck equation which governs the evolution of the transition probability is derived. The analytical solution for the transition probability is then obtained and the effects of stochastic fluctuation is discussed. The statistical and stochastic analyses may provide effective approaches to reveal the physics of damage evolution and dynamic failure process in ductile materials.

Study Of The Damage-Induced Anisotropy Of Quasi-Brittle Materials Using The Component Assembling Model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

Internal Friction Of Bend-Deformed Nanocrystalline Nickel By Mechanical Spectroscopy

Internal friction of nanocrystalline nickel is investigated by mechanical spectroscopy from 360 K to 120 K. Two relaxation peaks are found when nanocrystalline nickel is bent up to 10% strain at room temperature and fast cooling. However, these two peaks disappear when the sample is annealed at room temperature in vacuum for ten days. The occurrence and disappearance of the two relaxation peaks can be explained by the interactions of partial dislocations and point defects in nanocrystalline materials.

Nonlinear Constitutive Behavior Of Ferroelectric Materials

The ferroelectric specimen is considered as an aggregation of many randomly oriented domains. According to this mechanism, a multi-domain mechanical model is developed in this paper. Each domain is represented by one element. The applied stress and electric field are taken to be the stress and electric field in the formula of the driving force of domain switching for each element in the specimen. It means that the macroscopic switching criterion is used for calculating the volume fraction of domain switching for each element. By using the hardening relation between the driving force of domain switching and the volume fraction of domain switching calibrated, the volume fraction of domain switching for each element is calculated. Substituting the stress and electric field and the volume fraction of domain switching into the constitutive equation of ferroelectric material, one can easily get the strain and electric displacement for each element. The macroscopic behavior of the ferroelectric specimen is then directly calculated by volume averaging. Meanwhile, the nonlinear finite element analysis for the ferroelectric specimen is carried out. In the finite element simulation, the volume fraction of domain switching for each element is calculated by using the same method mentioned above. The interaction between different elements is taken into account in the finite element simulation and the local stress and electric field for each element is obtained. The macroscopic behavior of the specimen is then calculated by volume averaging. The computation results involve the electric butterfly shaped curves of axial strain versus the axial electric field and the hysteresis loops of electric displacement versus the electric field for ferroelectric specimens under the uniaxial coupled stress and electric field loading. The present theoretical prediction agrees reasonably with the experimental results.

Rapidly Infrared-Assisted Cooperatively Self-Assembled Highly Ordered Multiscale Porous Materials

In this paper, cooperative self-assembly (CSA) of colloidal spheres with different sizes was studied. It was found that a complicated jamming effect makes it difficult to achieve an optimal self-assembling condition for construction of a well-ordered stacking of colloidal spheres in a relatively short growth time by CSA. Through the use of a characteristic infrared (IR) technique to significantly accelerate local evaporation on the growing interface without changing the bulk growing environment, a concise three-parameter (temperature, pressure, and IR intensity) CSA method to effectively overcome the jamming effect has been developed. Mono- and multiscale inverse opals in a large range of lattice scales can be prepared within a growth time (15-30 min) that is remarkably shorter than the growth times of several hours for previous methods. Scanning electron microscopy images and transmittance spectra demonstrated the superior crystalline and optical qualities of the resulting materials. More importantly, the new method enables optimal conditions for CSA without limitations on sizes and materials of multiple colloids. This strategy not only makes a meaningful advance in the applicability and universality of colloidal crystals and ordered porous materials but also can be an inspiration to the self-assembly systems widely used in many other fields, such as nanotechnology and molecular bioengineering.

Effects Of The Length Of A Cylindrical Solid Shield On The Entrainment Of Ambient Air Into Turbulent And Laminar Impinging Argon Plasma Jets

When materials processing is conducted in air surroundings by use of an impinging plasma jet, the ambient air will be entrained into the materials processing region, resulting in unfavorable oxidation of the feedstock metal particles injected into the plasma jet and of metallic substrate material. Using a cylindrical solid shield may avoid the air entrainment if the shield length is suitably selected and this approach has the merit that expensive vacuum chamber and its pumping system are not needed. Modeling study is thus conducted to reveal how the length of the cylindrical solid shield affects the ambient air entrainment when materials processing (spraying, remelting, hardening, etc.) is conducted by use of a turbulent or laminar argon plasma jet impinging normally upon a flat substrate in atmospheric air. It is shown that the mass flow rate of the ambient air entrained into the impinging plasma jet cannot be appreciably reduced unless the cylindrical shield is long enough. In order to completely avoid the air entrainment, the gap between the downstream-end section of the cylindrical solid shield and the substrate surface must be carefully selected, and the suitable size of the gap for the turbulent plasma jet is appreciably larger than that for the laminar one. The overheating of the solid shield or the substrate could become a problem for the turbulent case, and thus additional cooling measure may be needed when the entrainment of ambient air into the turbulent impinging plasma jet is to be completely avoided.

Research Progress On Thermal Protection Materials And Structures Of Hypersonic Vehicles

Hypersonic vehicles represent future trends of military equipments and play an important role in future war. Thermal protection materials and structures, which relate to the safety of hypersonic vehicles, are one of the most key techniques in design and manufacture of hypersonic vehicles. Among these materials and structures, such as metallic temperature protection structure, the temperature ceramics and carbon/carbon composites are usually adopted in design. The recent progresses of research and application of ultra-high temperature materials in preparation, oxidation resistance, mechanical and physical characterization are summarized.

Non-Slipping Jkr Model For Transversely Isotropic Materials

A generalized plane strain JKR model is established for non-slipping adhesive contact between an elastic transversely isotropic cylinder and a dissimilar elastic transversely isotropic half plane, in which a pulling force acts on the cylinder with the pulling direction at an angle inclined to the contact interface. Full-coupled solutions are obtained through the Griffith energy balance between elastic and surface energies. The analysis shows that, for a special case, i.e., the direction of pulling normal to the contact interface, the full-coupled solution can be approximated by a non-oscillatory one, in which the critical pull-off force, pull-off contact half-width and adhesion strength can be expressed explicitly. For the other cases, i.e., the direction of pulling inclined to the contact interface, tangential tractions have significant effects on the pull-off process, it should be described by an exact full-coupled solution. The elastic anisotropy leads to an orientation-dependent pull-off force and adhesion strength. This study could not only supply an exact solution to the generalized JKR model of transversely isotropic materials, but also suggest a reversible adhesion sensor designed by transversely isotropic materials, such as PZT or fiber-reinforced materials with parallel fibers. (c) 2007 Elsevier Ltd. All rights reserved.

Threshold diversity and trans-scales sensitivity in a finite nonlinear evolution model of materials failure

We present a slice-sampling method and study the ensemble evolution of a large finite nonlinear system in order to model materials failure. There is a transitional region of failure probability. Its size effect is expressed by a slowly decaying scaling law. In a meso-macroscopic range (similar to 10(5)) in realistic failure, the diversity cannot be ignored. Sensitivity to mesoscopic details governs the phenomena. (C) 1997 Published by Elsevier Science B.V.

Numerical verification for instrumented spherical indentation techniques in determining the plastic properties of materials

 Instrumented indentation tests have been widely adopted for elastic modulus determination. Recently, a number of indentation-based methods for plastic properties characterization have been proposed, and rigorous verification is absolutely necessary for their wide application. In view of the advantages of spherical indentation compared with conical indentation in determining plastic proper-ties, this study mainly concerns verification of spherical indentation methods. Five convenient and simple models were selected for this purpose, and numerical experiments for a wide range of materials are carried out to identify their accuracy and sensitivity characteristics. The verification results show that four of these five methods can give relatively accurate and stable results within a certain material domain, which is defined as their validity range and has been summarized for each method.

Plane strain asymptotic fields for cracks terminating at the interface between elastic and pressure-sensitive dilatant materials

The plane strain asymptotic fields for cracks terminating at the interface between elastic and pressure-sensitive dilatant material are investigated in this paper. Applying the stress-strain relation for the pressure-sensitive dilatant material, we have obtained an exact asymptotic solution for the plane strain tip fields for two types of cracks, one of which lies in the pressure-sensitive dilatant material and the other in the elastic material and their tips touch both the bimaterial interface. In cases, numerical results show that the singularity and the angular variations of the fields obtained depend on the material hardening exponent n, the pressure sensitivity parameter mu and geometrical parameter lambda.

On the mechanism of void growth and the effect of straining mode in ductile materials

Finite element analysis is employed to investigate void growth embedded in elastic-plastic matrix material. Axisymmetric and plane stress conditions are considered. The simulation of void growth in a unit cell model is carried out over a wide range of triaxial tensile stressing or large plastic straining for various strain hardening materials to study the mechanism of void growth in ductile materials. Triaxial tension and large plastic strain encircling around the void are found to be of most importance for driving void growth. The straining mode of incremental loading which favors the necessary strain concentration around void for its growth can be characterized by the vanishing condition of a parameter called "the third invariant of generalized strain rate". Under this condition, it accentuates the internal strain concentration and the strain energy stored/dissipated within the material layer surrounding the void. Experimental results are cited to justify the effect of this loading parameter. (C) 2000 Elsevier Science Ltd. All rights reserved.

Materials selection in mechanical design for microsensors and microactuators

Microsensors and microactuators are vital organs of microelectromechanical systems (MEMS), forming the interfaces between controller and environment. They are usually used for devices ranging in size at sub-millimeter or micrometer level, transforming energy between two or more domains. Presently, most of the materials used in MEMS devices belong to the silicon material system, which is the basis of the integrated circuit industry. However, new techniques are being explored and developed, and the opportunities for MEMS materials selection are getting broader. The present paper tries to apply 'performance index' to select the material best suited to a given application, in the early stage of MEMS design. The selection is based on matching performance characteristics to the requirements. A series of performance indices are given to allow a wide range comparison of materials for several typical sensing and actuating structures, and a rapid identification of candidates for a given task. (C) 2002 Elsevier Science Ltd. All rights reserved.

Electromechanical influence of crack velocity at bifurcation for poled ferroelectric materials

The piezoelastodynamic field equations are solved to determine the crack velocity at bifurcation for poled ferroelectric materials where the applied electrical field and mechanical stress can be varied. The underlying physical mechanism, however, may not correspond to that assumed in the analytical model. Bifurcation has been related to the occurrence of a pair of maximum circumferential stress oriented symmetrically about the moving crack path. The velocity at which this behavior prevails has been referred to as the limiting crack speed. Unlike the classical approach, bifurcation will be identified with finite distances ahead of a moving crack. Nucleation of microcracks can thus be modelled in a single formulation. This can be accomplished by using the energy density function where fracture initiation is identified with dominance of dilatation in relation to distortion. Poled ferroelectric materials are selected for this study because the microstructure effects for this class of materials can be readily reflected by the elastic, piezoelectic and dielectric permittivity constants at the macroscopic scale. Existing test data could also shed light on the trend of the analytical predictions. Numerical results are thus computed for PZT-4 and compared with those for PZT-6B in an effort to show whether the branching behavior would be affected by the difference in the material microstructures. A range of crack bifurcation speed upsilon(b) is found for different r/a and E/sigma ratios. Here, r and a stand for the radial distance and half crack length, respectively, while E and a for the electric field and mechanical stress. For PZT-6B with upsilon(b) in the range 100-1700 m/s, the bifurcation angles varied from +/-6degrees to +/-39degrees. This corresponds to E/sigma of -0.072 to 0.024 V m/N. At the same distance r/a = 0.1, PZT-4 gives upsilon(b) values of 1100-2100 m/s; bifurcation angles of +/-15degrees to +/-49degrees; and E/sigma of -0.056 to 0.059 V m/N. In general, the bifurcation angles +/-theta(0) are found to decrease with decreasing crack velocity as the distance r/a is increased. Relatively speaking, the speed upsilon(b) and angles +/-theta(0) for PZT-4 are much greater than those for PZT-6B. This may be attributed to the high electromechanical coupling effect of PZT-4. Using upsilon(b)(0) as a base reference, an equality relation upsilon(b)(-) < upsilon(b)(0) < upsilon(b)(+) can be established. The superscripts -, 0 and + refer, respectively, to negative, zero and positive electric field. This is reminiscent of the enhancement and retardation of crack growth behavior due to change in poling direction. Bifurcation characteristics are found to be somewhat erratic when r/a approaches the range 10(-2)-10(-1) where the kinetic energy densities would fluctuate and then rise as the distance from the moving crack is increased. This is an artifact introduced by the far away condition of non-vanishing particle velocity. A finite kinetic energy density prevails at infinity unless it is made to vanish in the boundary value problem. Future works are recommended to further clarify the physical mechanism(s) associated with bifurcation by means of analysis and experiment. Damage at the microscopic level needs to be addressed since it has been known to affect the macrocrack speeds and bifurcation characteristics. (C) 2002 Published by Elsevier Science Ltd.

Fracture statistics of brittle materials: Weibull or normal distribution

The fit of fracture strength data of brittle materials (Si3N4, SiC, and ZnO) to the Weibull and normal distributions is compared in terms of the Akaike information criterion. For Si3N4, the Weibull distribution fits the data better than the normal distribution, but for ZnO the result is just the opposite. In the case of SiC, the difference is not large enough to make a clear distinction between the two distributions. There is not sufficient evidence to show that the Weibull distribution is always preferred to other distributions, and the uncritical use of the Weibull distribution for strength data is questioned.