Adaptive Mesh Refinement FEM for Damage Evolution of Heterogeneous Brittle Media

Damage evolution of heterogeneous brittle media involves a wide range of length scales. The coupling between these length scales underlies the mechanism of damage evolution and rupture. However, few of previous numerical algorithms consider the effects of the trans-scale coupling effectively. In this paper, an adaptive mesh refinement FEM algorithm is developed to simulate this trans-scale coupling. The adaptive serendipity element is implemented in this algorithm, and several special discontinuous base functions are created to avoid the incompatible displacement between the elements. Both the benchmark and a typical numerical example under quasi-static loading are given to justify the effectiveness of this model. The numerical results reproduce a series of characteristics of damage and rupture in heterogeneous brittle media.

Phase Evolution and its Effect on Magnetic Properties of Nd60Al10Fe20Co10 Bulk Metallic Glass

The thermal stability of nanocrystalline clusters, the phase evolution, and their effects on magnetic Propertieswere studied for as-cast Nd60Al10Fe20Co10 alloy using differential scanning calorimetry curves, x-ray diffraction patterns, scanning electron microscopy, and high-resolution transition electron microscopy. Thermomagnetic curves and hysteresis loops of the bulk metallic glass were measured during the annealing process. The high thermostability of the hardmagnetic properties of the samples observed is attributed to the stability of the nanocrystalline clusters upon annealing, while the slight enhancement in the magnetization is due to the precipitation of some Nd-rich metastable phases. The mechanism of thermostability of the nanocrystalline clusters and the formation of the metastable phases are discussed.

Numerical simulation of temperature field evolution of metallic alloy surface due to pulsed laser treatment

This paper performed a numerical simulation on temperature field evolution for the surface layer of a metallic alloy subjected to pulsed Nd:YAG laser treatment. The enthalpy method was adopted to solve the moving boundary problem, I.e. Stefan problem. Computational results were obtained to show the temperature field evolution. Effects of latent heat and mushy zone width on the temperature field were investigated. The results also show very high values of temperature gradient and cooling rate, which are typical characteristics during the solidification process.

Electron Beam Surface Treatment. Part II: Microstructure Evolution of Stainless Steel and Aluminum Alloy During Electron Beam Rapid Solidification

Surface rapid solidification microstructures of AISI 321 austenitic stainless steel and 2024 aluminum alloy have been investigated by electron beam remelting process and optical microscopy observation. It is indicated that the morphologies of the melted layer of both stainless steel and aluminum alloy change dramatically compared to the original materials. Also, the microstructures were greatly refined after the electron beam irradiation.

Microdamage Evolution, Energy Dissipation and Their Trans-Scale Effects on Macroscopic Failure

The process of damage evolution concerns various scales, from micro- to macroscopic. How to characterize the trans-scale nature of the process is on the challenging frontiers of solid mechanics. In this paper, a closed trans-scale formulation of damage evolution based on statistical microdamage mechanics is presented. As a case study, the damage evolution in spallation is analyzed with the formulation. Scaling of the formulation reveals that the following dimensionless numbers: reduced Mach number M, damage number S, stress wave Fourier number P, intrinsic Deborah number D*, and the imposed Deborah number De*, govern the whole process of deformation and damage evolution. The evaluation of P and the estimation of temperature increase show that the energy equation can be ignored as the first approximation in the case of spallation. Hence, apart from the two conventional macroscopic parameters: the reduced Mach number M and damage number S, the damage evolution in spallation is mainly governed by two microdamage-relevant parameters: the Deborah numbers D* and De*. Higher nucleation and growth rates of microdamage accelerate damage evolution, and result in higher damage in the target plate. In addition, the mere variation in nucleation rate does not change the spatial distribution of damage or form localized rupture, while the increase of microdamage growth rate localizes the damage distribution in the target plate, which can be characterized by the imposed Deborah number De*.

Dislocation evolution in epitaxial multilayers and graded composition buffers

This paper presents models to describe the dislocation dynamics of strain relaxation in an epitaxial uniform layer, epitaxial multilayers and graded composition buffers. A set of new evolution equations for nucleation rate and annihilation rate of threading dislocations is developed. The dislocation interactions are incorporated into the kinetics process by introducing a resistance term, which depends only on plastic strain. Both threading dislocation nucleation and threading dislocation annihilation are characterized. The new evolution equations combined with other evolution equations for the plastic strain rate, the mean velocity and the dislocation density rate of the threading dislocations are tested on GexSi1-x/Si(100) heterostructures, including epitaxial multilayers and graded composition buffers. It is shown that the evolution equations successfully predict a wide range of experimental results of strain relaxation and threading dislocation evolution in the materials system. Meanwhile, the simulation results clearly signify that the threading dislocation annihilation plays a vital role in the reduction of threading dislocation density.

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.

Numerical Investigation On Detonation Cell Evolution In A Channel With Area-Changing Cross Section

The two-dimensional cellular detonation propagating in a channel with area-changing cross section was numerically simulated with the dispersion-controlled dissipative scheme and a detailed chemical reaction model. Effects of the flow expansion and compression on the cellular detonation cell were investigated to illustrate the mechanism of the transverse wave development and the cellular detonation cell evolution. By examining gas composition variations behind the leading shock, the chemical reaction rate, the reaction zone length, and thermodynamic parameters, two kinds of the abnormal detonation waves were identified. To explore their development mechanism, chemical reactions, reflected shocks and rarefaction waves were discussed, which interact with each other and affect the cellular detonation in different ways.

Shear Localization In Dynamic Deformation: Microstructural Evolution

Investigations made by the authors and collaborators into the microstructural aspects of adiabatic shear localization are critically reviewed. The materials analyzed are low-carbon steels, 304 stainless steel, monocrystalline Fe-Ni-Cr, Ti and its alloys, Al-Li alloys, Zircaloy, copper, and Al/SiCp composites. The principal findings are the following: (a) there is a strain-rate-dependent critical strain for the development of shear bands; (b) deformed bands and white-etching bands correspond to different stages of deformation; (c) different slip activities occur in different stages of band development; (d) grain refinement and amorphization occur in shear bands; (e) loss of stress-carrying capability is more closely associated with microdefects rather than with localization of strain; (f) both crystalline rotation and slip play important roles; and (g) band development and band structures are material dependent. Additionally, avenues for new research directions are suggested.

Numerical Investigation On Evolution Of Cylindrical Cellular Detonation

Cylindrical cellular detonation is numerically investigated by solving two-dimensional reactive Euler equations with a finite volume method on a two-dimensional self-adaptive unstructured mesh. The one-step reversible chemical reaction model is applied to simplify the control parameters of chemical reaction. Numerical results demonstrate the evolution of cellular cell splitting of cylindrical cellular detonation explored in experimentas. Split of cellular structures shows different features in the near-field and far-field from the initiation zone. Variation of the local curvature is a key factor in the behavior of cell split of cylindrical cellular detonation in propagation. Numerical results show that split of cellular structures comes from the self-organization of transverse waves corresponding to the development of small disturbances along the detonation front related to detonation instability.

Damage Evolution In Cement Mortar Due To Erosion Of Sulphate

Ultrasonic technique is used to detect the velocity change of stress wave propagated in the cement mortar immersed in the solution of sodium sulfate for 425 days. Also the density change of specimens at different erosion time is measured. By curve fitting, the effect of solutions' concentration and water/cement ratio on the damage evolution is analyzed. The SEM observation on the growth of delayed ettringite is also performed. It shows that the damage evolution of specimens attacked by sulphate solution is dominantly induced by the nucleation and growth of delayed ettringite, and the average size of microvoids in cement mortar affects the damage evolution significantly. (c) 2008 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.

Evolution of thermoplastic shear localization and related microstructures in Awl/SiCp composites under dynamic compression

The localized shear deformation in the 2024 and 2124 Al matrix composites reinforced with SiC particles was investigated with a split Hopkinson pressure bar (SHPB) at a strain rate of about 2.0x10(3) s(-1). The results showed that the occurrence of localized shear deformation is sensitive to the size of SiC particles. It was found that the critical strain, at which the shear localization occurs, strongly depends on the size and volume fraction of SiC particles. The smaller the particle size, the lower the critical strain required for the shear localization. TEM examinations revealed that Al/SiCp interfaces are the main sources of dislocations. The dislocation density near the interface was found to be high and it decreases with the distance from the particles. The Al matrix in shear bands was highly deformed and severely elongated at low angle boundaries. The Al/SiCp interfaces, particularly the sharp corners of SiC particles, provide the sites for microcrack initiation. Eventual fracture is caused by the growth and coalescence of microcracks along the shear bands. It is proposed that the distortion free equiaxed grains with low dislocation density observed in the center of shear band result from recrystallization during dynamic deformation.

Long-range automaton models of earthquakes: Power-law accelerations, correlation evolution, and mode-switching

We introduce a conceptual model for the in-plane physics of an earthquake fault. The model employs cellular automaton techniques to simulate tectonic loading, earthquake rupture, and strain redistribution. The impact of a hypothetical crustal elastodynamic Green's function is approximated by a long-range strain redistribution law with a r(-p) dependance. We investigate the influence of the effective elastodynamic interaction range upon the dynamical behaviour of the model by conducting experiments with different values of the exponent (p). The results indicate that this model has two distinct, stable modes of behaviour. The first mode produces a characteristic earthquake distribution with moderate to large events preceeded by an interval of time in which the rate of energy release accelerates. A correlation function analysis reveals that accelerating sequences are associated with a systematic, global evolution of strain energy correlations within the system. The second stable mode produces Gutenberg-Richter statistics, with near-linear energy release and no significant global correlation evolution. A model with effectively short-range interactions preferentially displays Gutenberg-Richter behaviour. However, models with long-range interactions appear to switch between the characteristic and GR modes. As the range of elastodynamic interactions is increased, characteristic behaviour begins to dominate GR behaviour. These models demonstrate that evolution of strain energy correlations may occur within systems with a fixed elastodynamic interaction range. Supposing that similar mode-switching dynamical behaviour occurs within earthquake faults then intermediate-term forecasting of large earthquakes may be feasible for some earthquakes but not for others, in alignment with certain empirical seismological observations. Further numerical investigation of dynamical models of this type may lead to advances in earthquake forecasting research and theoretical seismology.

Micromechanics analysis on evolution of crack in viscoelastic materials

A preliminary analysis on crack evolution in viscoelastic materials was presented. Based on the equivalent inclusion concept of micro-mechanics theory, the explicit expressions of crack opening displacement delta and energy release rate G were derived, indicating that both delta and G are increasing with time. The equivalent modulus of the viscoelastic solid comprising cracks was evaluated. It is proved that the decrease of the modulus comes from two mechanisms: one is the viscoelasticity of the material; the other is the crack opening which is getting larger with time.