Embedded cell model of three-dimensional two-phase composites

An embedded cell model is presented to obtain the effective elastic moduli and the elastic-plastic stress-strain relations of three-dimensional two-phase particulate composites. Each cell consists of an ellipsoidal inclusion surrounded by a finite ellipsoidal matrix that embedded in an infinite matrix. When both matrix and particle are elastic, the effective elastic moduli are derived which is an exact analytic formula without any simplified approximation that can be expressed in an explicit form. Further, the elastic-plastic stress-strain relations are obtained for spherical cells and oblate spheroid cells, in which the matrix is elastic and the particle is elastic-plastic. In addition, the macroscopic elastic-plastic constitutive relation of particle reinforced composites (PRC) is investigated by a systematic approach [1] in which the matrix is elastic-plastic and the particle is elastic.

Revelation of a Functional Dependence of the Sum of Two Uniaxial Strengths/Hardness on Elastic Work/Total Work of Indentation

Dimensional and finite element analyses were used to analyze the relationship between the mechanical properties and instrumented indentation response of materials. Results revealed the existence of a functional dependence of (engineering yield strength sigma(E,y) + engineering tensile strength sigma(E,b))/Oliver & Pharr hardness on the ratio of reversible elastic work to total work obtained from an indentation test. The relationship links up the Oliver & Pharr hardness with the material strengths, although the Oliver & Pharr hardness may deviate from the true hardness when sinking in or piling up occurs. The functional relationship can further be used to estimate the SUM sigma(E,y) + sigma(E,b) according to the data of an instrumented indentation test. The sigma(E,y) + sigma(E,b) value better reflects the strength of a material compared to the hardness value alone. The method was shown to be effective when applied to aluminum alloys. The relationship can further be used to estimate the fatigue limits, which are usually obtained from macroscopic fatigue tests in different modes.

Mesofracture mechanics: a necessary link

Classical fracture mechanics is based on the premise that small scale features could be averaged to give a larger scale property such that the assumption of material homogeneity would hold. Involvement of the material microstructure, however, necessitates different characteristic lengths for describing different geometric features. Macroscopic parameters could not be freely exchanged with those at the microscopic scale level. Such a practice could cause misinterpretation of test data. Ambiguities arising from the lack of a more precise range of limitations for the definitions of physical parameters are discussed in connection with material length scales. Physical events overlooked between the macroscopic and microscopic scale could be the link that is needed to bridge the gap. The classical models for the creation of free surface for a liquid and solid are oversimplified. They consider only the translational motion of individual atoms. Movements of groups or clusters of molecules deserve attention. Multiscale cracking behavior also requires the distinction of material damage involving at least two different scales in a single simulation. In this connection, special attention should be given to the use of asymptotic solution in contrast to the full field solution when applying fracture criteria. The former may leave out detail features that would have otherwise been included by the latter. Illustrations are provided for predicting the crack initiation sites of piezoceramics. No definite conclusions can be drawn from the atomistic simulation models such as those used in molecular dynamics until the non-equilibrium boundary conditions can be better understood. The specification of strain rates and temperatures should be synchronized as the specimen size is reduced to microns. Many of the results obtained at the atomic scale should be first identified with those at the mesoscale before they are assumed to be connected with macroscopic observations. Hopefully, "mesofracture mechanics" could serve as the link to bring macrofracture mechanics closer to microfracture mechanics.

Critical Sensitivity and Trans-Scale Fluctuations in Catastrophic Rupture

Rupture in the heterogeneous crust appears to be a catastrophe transition. Catastrophic rupture sensitively depends on the details of heterogeneity and stress transfer on multiple scales. These are difficult to identify and deal with. As a result, the threshold of earthquake-like rupture presents uncertainty. This may be the root of the difficulty of earthquake prediction. Based on a coupled pattern mapping model, we represent critical sensitivity and trans-scale fluctuations associated with catastrophic rupture. Critical sensitivity means that a system may become significantly sensitive near catastrophe transition. Trans-scale fluctuations mean that the level of stress fluctuations increases strongly and the spatial scale of stress and damage fluctuations evolves from the mesoscopic heterogeneity scale to the macroscopic scale as the catastrophe regime is approached. The underlying mechanism behind critical sensitivity and trans-scale fluctuations is the coupling effect between heterogeneity and dynamical nonlinearity. Such features may provide clues for prediction of catastrophic rupture, like material failure and great earthquakes. Critical sensitivity may be the physical mechanism underlying a promising earthquake forecasting method, the load-unload response ratio (LURR).

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.

Toughness Of Ni/Al2O3 Interfaces As Dependent On Micron-Scale Plasticity And Atomistic-Scale Separation

Ceramic/metal interfaces were studied that fail by atomistic separation accompanied by plastic dissipation in the metal. The macroscopic toughness of the specific Ni alloy/Al2O3 interface considered is typically on the order of ten times the atomistic work of separation in mode I and even higher if combinations of mode I and mode II act on the interface. Inputs to the computational model of interface toughness are: (i) strain gradient plasticity applied to the Ni alloy with a length parameter determined by an indentation test, and (ii) a potential characterizing mixed mode separation of the interface fit to atomistic results. The roles of the several length parameters in the strain gradient plasticity are determined for indentation and crack growth. One of the parameters is shown to be of dominant importance, thus establishing that indentation can be used to measure the relevant length parameter. Recent results for separation of Ni/Al2O3 interfaces computed by atomistic methods are reviewed, including a set of results computed for mixed mode separation. An approximate potential fit to these results is characterized by the work of separation, the peak separation stress for normal separation and the traction-displacement relation in pure shearing of the interface. With these inputs, the model for steady-state crack growth is used to compute the toughness of the interface under mode I and under the full range of mode mix. The effect of interface strength and the work of separation on macroscopic toughness is computed. Fundamental implications for plasticity-enhanced toughness emerge.

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.

Stability Analysis Of The Moving Interface In Piston- And Non-Piston-Like Displacements

From the macroscopic point of view, expressions involving reservoir and operational parameters are established for investigating the stability of moving interface in piston- and non-piston-like displacements. In the case of axi-symmetrical piston-like displacement, the stability is related to the moving interface position and water to oil mobility ratio. The capillary effect on the stability of moving interface depends on whether or not the moving interface is already stable and correlates with the wettability of the reservoir rock. In the case of non-piston-like displacement, the stability of the front is governed by both the relative permeability and the mobility ratio.

Role Of Convection Flow On The Pattern Formation In The Drying Process Of Colloidal Suspension

We study the macroscopic drying patterns of aqueous suspensions of colloidal silica spheres. It was found that convection strength can influence pattern formation. Uniformed films are obtained at weaker convection strength. In addition, we make clear that it is not reasonable to discuss individually the effect of temperature and humidity on the colloid self-assembly. The physical mechanism is that these factors have relationship with the evaporation rate, which can affect the convection strength.

Reconstructing The Vertical Distribution Of The Aeolian Saltation Mass Flux Based On The Probability Distribution Of Lift-Off Velocity

The probability distribution of lift-off velocity of the saltating grains is a bridge to linking microscopic and macroscopic research of aeolian sand transport. The lift-off parameters of saltating grains (i.e., the horizontal and vertical lift-off velocities, resultant lift-off velocity, and lift-off angle) in a wind tunnel are measured by using a Phase Doppler Particle Analyzer (PDPA). The experimental results show that the probability distribution of horizontal lift-off velocity of saltating particles on a bed surface is a normal function, and that of vertical lift-off velocity is an exponential function. The probability distribution of resultant lift-off velocity of saltating grains can be expressed as a log-normal function, and that of lift-off angle complies with an exponential function. A numerical model for the vertical distribution of aeolian mass flux based on the probability distribution of lift-off velocity is established. The simulation gives a sand mass flux distribution which is consistent with the field data of Namikas (Namikas, S.L., 2003. Field measurement and numerical modelling of acolian mass flux distributions on a sandy beach, Sedimentology 50, 303-326). Therefore, these findings are helpful to further understand the probability characteristics of lift-off grains in aeolian sand transport. (c) 2007 Elsevier B.V. All rights reserved.

Deformation Behavior Of Fe-Based Bulk Metallic Glass During Nanoindentation

Fe-based bulk metallic glasses (BMGs) normally exhibit super high strength but significant brittleness at ambient temperature. Therefore, it is difficult to investigate the plastic deformation behavior and mechanism in these alloys through conventional tensile and compressive tests due to lack of distinct macroscopic plastic strain. In this work, the deformation behavior of Fe52Cr15Mo9Er3C15B6 BMG was investigated through instrumented nanoindentation and uniaxial compressive tests. The results show that serrated flow, the typical plastic deformation feature of BMGs, could not be found in as-cast and partially crystallized samples during nanoindentation. In addition, the deformation behavior and mechanical properties of the alloy are insensitive to the applied loading rate. The mechanism for the appearance of the peculiar deformation behavior in the Fe-based BMG is discussed in terms of the temporal and spatial characteristics of shear banding during nanoindentation.

Theory Of Phase Transformation And Reorientation In Single Crystalline Shape Memory Alloys

A constitutive model, based on an (n + 1)-phase mixture of the Mori-Tanaka average theory, has been developed for stress-induced martensitic transformation and reorientation in single crystalline shape memory alloys. Volume fractions of different martensite lattice correspondence variants are chosen as internal variables to describe microstructural evolution. Macroscopic Gibbs free energy for the phase transformation is derived with thermodynamics principles and the ensemble average method of micro-mechanics. The critical condition and the evolution equation are proposed for both the phase transition and reorientation. This model can also simulate interior hysteresis loops during loading/unloading by switching the critical driving forces when an opposite transition takes place.

A Modified Free Volume Model For Characterizing Of Rate Effect In Bulk Metallic Glasses

We investigate the plastic deformation and constitutive behaviour of bulk metallic glasses (BMGs). A dimensionless Deborah number De(ID) = t(r)/t(i) is proposed to characterize the rate effect in BMGs, where t(r) is the structural relaxing characteristic time of BMGs under shear load, t(i) is the macroscopic imposed characteristic time of applied stress or the characteristic time of macroscopic deformation. The results demonstrate that the modified free volume model can characterize the strain rate effect in BMGs effectively.

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.

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.