Fracture criteria for combined cleavage and dislocation emission

A general theory of fracture criteria for mixed dislocation emission and cleavage processes is developed based on Ohr's model. Complicated cases involving mixed-mode loading are considered. Explicit formulae are proposed for the critical condition of crack cleavage propagation after a number of dislocation emissions. The effects of crystal orientation, crack geometry and load phase angle on the apparent critical energy release rates and the total number of the emitted dislocations at the initiation of cleavage are analysed in detail. In order to evaluate the effects of nonlinear interaction between the slip displacement and the normal separation, an analysis of fracture criteria for combined dislocation emission and cleavage is presented on the basis of the Peierls framework. The calculation clearly shows that the nonlinear theory gives slightly high values of the critical apparent energy release rate G(c) for the same load phase angle. The total number N of the emitted dislocations at the onset of cleavage given by nonlinear theory is larger than that of linear theory.

Molecular dynamics simulation of interaction of a dislocation array from a crack tip with grain boundaries

The interaction of a dislocation array emitted from a crack tip under mode II loading with asymmetric tilt grain boundaries (GBs) is analysed by the molecular dynamics method. The GBs can generally be described by planar and linear matching zones and unmatching zones. All GBs are observed to emit dislocations. The GBs migrated easily due to their planar and linear matching structure and asymmetrical type. The diffusion induced by stress concentration is found to promote the GB migration. The transmissions of dislocations are either along the matched plane or along another plane depending on tilt angle theta. Alternate processes of stress concentration and stress relaxation take place ahead of the pileup. The stress concentration can be released either by transmission of dislocations, by atom diffusion along GBs, or by migration of GBs by formation of twinning bands. The simulated results also unequivocally demonstrate two processes, i.e. asymmetrical GBs evolving into symmetrical ones and unmatching zones evolving into matching ones during the loading process.

Formation, microstructure and development of the localized shear deformation in low-carbon steels

An investigation has been made into the effect of microstructural parameters on the propensity for forming shear localization produced during high speed torsional testing by split Hopkinson bar with different average rates of 610, 650 and 1500 s(-1) in low carbon steels. These steels received the quenched, quenched and tempered as well as normalized treatments that provide wide microstructural parameters and mechanical properties. The results indicate that the occurrence of the shear localization is susceptible to the strength of the steels. In other words, the tendency of the quenched steel to form a shear band is higher than that of the other two steels. It is also found that there is a critical strain at which the shear localization occurs in the steels. The critical strain value is strongly dependent on the strength of the steels. Before arriving at this point, the material undergoes a slow work-hardening. After this point, the material suffers work-softening, corresponding to a process during which the deformation is gradually localized and eventually becomes spatially correlated to form a macroscopic shear band. Examinations by SEM reveal that the shear localization within the band involves a series of sequential crystallographic and non-crystallographic events including the change in crystal orientation, misorientation, generation and even perhaps damage in microstructures such as the initiation, growth and coalescence of the microcracks. It is expected that the sharp drop in the load-carrying capacity is associated with the growth and coalescence of the microcracks rather than the occurrence of the shear localization, but the shear localization is seen to accelerate the growth and coalescence of the microcracks. The thin foil observations by TEM reveal that the density of dislocations in the band is extremely high and the tangled arrangement and cell structure of dislocations tends to align along the shear direction. The multiplication and interaction of dislocations seems to be responsible for work-hardening of the steels. The avalanche of the dislocation cells corresponds to the sharp drop in shear stress at which the deformed specimen is broken. Double shear bands and kink bands are also observed in the present study. The principal band develops first and its width is narrower than that of the secondary band.

Molecular dynamics studies on the dislocation gliding near a tilt boundary

The gliding behavior of edge dislocation near a grain boundary(QB) in copper under pure shear stresses is simulated by using molecular dynamics(MD) method. Many-body potential incorporating the embedded atom method (EAM) is used. The critical shear stresses for a single disocation to pass across GB surface are obtained at values of sigma(c)=23MPa similar to 68 MPa and 137 MPa similar to 274 MPa for Sigma=165 small angle tilt GB at 300 K and 20 K, respectively. The first result agrees with the experimental yield stress sigma(y)(=42 MPa) quite well. It suggests that there might be one of the reasons of initial plastic yielding caused by single dislocation gliding across GB. In addition, there might be possibility to obtain yield strength from microscopic analysis. Moreover, the experimental value of sigma(y) at low temperature is generally higher than that at room temperature. So, these results are in conformity qualitatively with experimental fact. On the other hand, the Sigma=25 GB is too strong an obstacle to the dislocation. In this case, a dislocation is able to pass across GB under relatively low stress only when it is driven by other dislocations. This is taken to mean that dislocation pile-up must be built up in front of this kind of GB, if this GB may take effect on the process of plastic deformation.

Experimental and theoretical study on numerical density evolution of short fatigue cracks

Fatigue testing was performed using a kind of triangular shaped specimen to obtain the characteristics of numerical density evolution for short cracks at the primary stage of fatigue damage. The material concerned is a structural alloy steel. The experimental results show that the numerical density of short cracks reaches the maximum value when crack length is slightly less than the average grain diameter, indicating grain boundary is the main barrier for short crack extension. Based on the experimental observations and related theory, the expressions for growth velocity and nucleation rate of short cracks have been proposed. With the solution to phase space conservation equation, the theoretical results of numerical density evolution for short cracks were obtained, which were in agreement with our experimental measurements.

The effect of thermal-activation on dislocation processes at an atomistic crack-tip

The effects of thermal activation on the dislocation emission from an atomistic crack tip are discussed, Molecular dynamics simulations at different constant temperatures are carried out to investigate the thermal effects. The simulated results show that the processes of the partial dislocation generation and emission are temperature dependent. As the temperature increases, the incipient duration of the partial dislocation nucleation becomes longer, the critical stress intensity factor for partial dislocation emission is reduced and, at the same loading level, more dislocations are emitted. The dislocation velocity moving away from the crack tip and the separations of partial dislocations are apparently not temperature dependent. The simulated results also show that, as the temperature increases, the stress distribution along the crack increases slightly. Therefore stress softening at the crack tip induced by thermal activation does not exist in the present simulation. A simple model is proposed to evaluate the relation of the critical stress intensity factor versus temperature. The obtained relation is in good agreement with our molecular dynamics results.

Behavior of ar plasma formed in a high-density plasma source - an ecr reactor

In order to develop the ultra-large scale integration(ULSI), low pressure and high density plasma apparatus are required for etching and deposit of thin films. To understand critical parameters such as the pressure, temperature, electrostatic potential and energy distribution of ions impacting on the wafer, it is necessary to understand how these parameters are influenced by the power input and neutral gas pressure. In the present work, a 2-D hybrid electron fluid-particle ion model has been developed to simulate one of the high density plasma sources-an Electron Cyclotron Resonance (ECR) plasma system with various pressures and power inputs in a non-uniform magnetic field. By means of numerical simulation, the energy distributions of argon ion impacting on the wafer are obtained and the plasma density, electron temperature and plasma electrostatic potential are plotted in 3-D. It is concluded that the plasma density depends mainly on both the power input and neutral gas pressure. However, the plasma potential and electron temperature can hardly be affected by the power input, they seem to be primarily dependent on the neutral gas pressure. The comparison shows that the simulation results are qualitatively in good agreement with the experiment measurements.

Potential energy model for turbulent entrainment across density interface in a linearly stratified fluid

A potential energy model is developed for turbulent entrainment in the absence of mean shear in a linearly stratified fluid. The relation between the entrainment distance D and the time t and the relation between dimensionless entrainment rate E and the local Richardson number are obtained. An experiment is made for examination. The experimental results are in good agreement with the model, in which the dimensionless entrainment distance D is given by DBAR = A(i)(SBAR)-1/4(fBAR)1/2(tBAR)1/8, where A(i) is the proportional coefficient, S is the dimensionless stroke, fBAR is the dimensionless frequency of the grid oscillation, tBAR the dimensionless time.

Influence of suspended particles and particle cover over the density interface

The effect of the particle cover over the density interface between two layers of fluids and of the suspended solid particles in the upper turbulcnt layer on the turbulent entrainment has been studied experimentally. The entrainment distance D is a function of the time of power: D=kt, where =0.200-0.130p. For suspended particles in the upper layer and pure 2-layer fluid is equal to 0.200, but the value of k for the suspended particles is smaller than that for the pure 2-layer fluid. The non-dimensional entrainment velocity is E=KRiln, where n=1.50+0.93 p. It is shown that the particle cover over the interface changes the power of Ril in the entrainment and hinders the turbulent entrainment. The variation rule of E for the suspended particles is the same as that for the pure 2-layer fluid, but the K value of the former is smaller than that of the latter. The turbulent mixing mechanism has been discussed.