Molecular dynamics simulation of crack growth under cyclic loading

The mechanical behaviors around a crack tip for a system including both a crack and two tilt grain boundaries under cyclic loading are examined using a molecular dynamics simulation. Not only a phase transition but also the emission of edge dislocations is observed in order to relax stress concentration around a crack tip during the first loading. Then, a dislocation pile-up is formed near the grain boundary after the edge dislocations reach the grain boundary, because they cannot move beyond the grain boundary. During the first unloading, the edge dislocations emitted from the crack tip return to the crack tip and disappear in the system. We observe several vacancies generated around the crack tip and crack growth corresponding to an atomic scale during cyclic loading. Conclusively, we propose the fatigue crack growth mechanism for the initial phase of the fatigue fracture. That is, a fatigue crack propagates due to coalescence of the crack and the vacancies caused by the emission and absorption of dislocations.     

A PROPOSED CRITERION FOR FATIGUE THRESHOLD: DISLOCATION SUBSTRUCTURE APPROACH

Abstract— A model for fatigue threshold has been proposed based on the dislocation subgrain cell structure that evolves at the crack tip in steels during the fatigue deformation process. The stabilized subgrain cells that develop in the material act as impenetrable barriers to dislocations in slip band pile-ups that emanate from the fatigue crack tip. The blocking of these dislocations tends to limit crack growth that occurs by crack tip emission of dislocations, thereby leading ultimately to the fatigue threshold condition. The grain size effect on threshold is deduced to be an indirect effect as it is proposed that the subgrain cell size is the controlling substructural parameter at the threshold stress intensity level. The subgrain cell size is shown to be proportional to the one-third power of the initial grain size.     

A FIRST STAGE IN THE DEVELOPMENT OF MICROMECHANICAL SIMULATIONS OF THE CRYSTALLOGRAPHIC PROPAGATION OF FATIGUE CRACKS UNDER MULTIAXIAL LOADING

Simulations of the nucleation of dislocations, glide and annihilation ahead of a fatigue crack growing along a localized slip band (a 'long' Stage I crack or a Stage II crack with a K value close to the threshold) are performed for the case of push–pull or reversed torsion loadings, ignoring, in a first approach, the effect of grain boundaries. The crack growth rates are deduced from the dislocation flux at the crack tip. An influence of the normal stress on the friction between the crack flanks as well as on the condition for dislocation emission is introduced. A slower Stage I growth rate is then predicted for reversed torsion, consistent with experimental data.     

Micromechanical simulations of microstructure-sensitive Stage I fatigue crack growth

Simulations of dislocation dynamics at the tip of a Stage I crack are performed, taking into account the influence of the normal stress on the friction of the crack flanks and on the condition for dislocation emission at the crack tip. The interactions of the emitted dislocations with microstructural obstacles are analysed. The repeated decelerations and sometimes arrests that characterize Stage I crack growth are properly described by the model, and the differences in Stage I kinetics observed in reversed torsion and push–pull are analysed in terms of crack tip–grain boundary interactions.     

Comparison of elastic interaction of a dislocation and a crack for four bonding conditions of the crack plane

A comparison of elastic interaction of a dislocation and a crack for four bonding conditions of the crack plane was made. Four cases of single crystalline material, sliding grain boundary, perfectly bonded interface, and sliding interface were considered. The stress intensity factors arising from edge and screw dislocations and their image forces for the above four cases were compared. The stress intensity factor at a crack tip along the perfectly bonded interface arising from screw dislocation can be obtained from that in a single crystalline material if the shear modulus in the single crystalline material is replaced by the harmonic mean of both shear moduli in the bimaterial. The stress intensity factor at a crack tip along the sliding interface arising from edge dislocation in the bimaterial can be obtained from that along the sliding grain boundary in the single material if the μ/(1−ν) in the single material is substituted by the harmonic mean of μ/(1− ν) in the bimaterial where μ and ν are the shear modulus and Poisson's ratio, respectively. The solutions of screw dislocation near a crack along the sliding grain boundary and sliding interface are the same as that of screw dislocation and its mirror image. Generally, the effect of edge dislocation for perfectly bonded interface on the crack propagation is more pronounced than that for the sliding interface. The effect of edge dislocation on the crack propagation is mixed mode for the cases of perfectly bonded interface and single crystalline material, but mode I fracture for the cases of sliding interface and sliding grain boundary. All curves of Fx versus distance r from the dislocation at interface to the right-hand crack tip are similar to one another regardless of dislocation source for both sliding interface and perfectly bonded interface. The level of Fx for m=0 is larger than that for m=−1. This revised version was published online in August 2006 with corrections to the Cover Date.     

Grain boundary influence on short fatigue crack growth rate

The influence from different grain boundary configurations on the crack growth rate of a microstructurally short edge crack, located within one grain and subjected to remote fatigue loading, is studied. The study is performed using a dislocation formulation, were the geometry is described by dislocation dipole elements in a boundary element approach and the plasticity by discrete dislocations, located along specific slip planes in the material. Plane strain and quasi-static conditions are assumed. The crack is assumed to grow in a single shear mechanism due to nucleation, glide and annihilation of discrete dislocations. Different grain boundary configurations in front of the growing crack are considered, including both high angle and low angle grain boundaries. It is shown that both grain boundary configuration and distance between the crack and a grain boundary has a pronounced influence on the crack growth rate.     

SiCw/Al复合材料在TEM下的拉伸行为

在备有拉伸台的透射电镜上原位观察了SiCw/Al复合材料的变形、裂纹萌生及其扩展过程.结果发现,变形初期,位错主要从晶界及SiCw-Al界面处开动,在晶界与晶须周围产生大量位错.裂纹萌生源主要为:(1)材料制备过程中受损伤或破断的晶须;(2)晶须富集区.晶须与晶界对裂纹的扩展有阻碍作用.当裂纹取向与晶须大致垂直时,裂尖在晶须附近钝化然后绕过晶须扩展,当裂纹取向与晶须接近平行时,裂尖在晶须端部附近钝化后沿其一侧扩展.     

The role of local strains from prior cold work on stress corrosion cracking of α-brass in Mattsson's solution

The dynamic strain rate ahead of a crack tip formed during stress corrosion cracking (SCC) under a static load is assumed to arise from the crack propagation. The strain surrounding the crack tip would be redistributed as the crack grows, thereby having the effect of dynamic strain. Recently, several studies have shown cold work to cause accelerated crack growth rates during SCC, and the slip-dissolution mechanism has been widely applied to account for this via a supposedly increased crack-tip strain rate in cold worked material. While these interpretations consider cold work as a homogeneous effect, dislocations are generated inhomogeneously within the microstructure during cold work. The presence of grain boundaries results in dislocation pile-ups that cause local strain concentrations. The local strains generated from cold working α-brass by tensile elongation were characterized using electron backscatter diffraction (EBSD). The role of these local strains in SCC was studied by measuring the strain distributions from the same regions of the sample before cold work, after cold work, and after SCC. Though, the cracks did not always initiate or propagate along boundaries with pre-existing local strains from the applied cold work, the local strains surrounding the cracked boundaries had contributions from both the crack propagation and the prior cold work.     

Influence from geometrical features on the growth rate of a micro-structurally short fatigue crack

The influence on the crack growth rate on a micro-structurally short edge crack subjected to fatigue loading from changes in crack length, distance to grain boundaries and applied load has been investigated. The crack is assumed to grow in a single shear mechanism due to nucleation, glide and annihilation of dislocations along preferred slip planes in the material. The external geometry is modelled by distributed dislocation dipole elements in a boundary element approach under quasi-static and plane strain conditions. The evolving plasticity is described by individual discrete dislocations along a slip plane emanating from the crack in the crack direction. The crack growth rate is shown to be controlled by the plasticity, which in turn is controlled by geometrical parameters in combination with the external load.     

Shielding effects and residual stresses at cleavage due to pre-existing dislocations

The motion of pre-existing edge dislocations in an infinite linear elastic body is studied at initiation of crack growth and at quasi-static steady-state crack growth. Dislocation nucleation is assumed not to occur. Thus, the study concerns only dislocations that are present in the virgin material. A dislocation is assumed to glide if its driving force exceeds a critical value. Changes in dislocation density, crack tip shielding and residual stresses are obtained. The shielding of a stationary crack tip is found to be small compared with the shielding of a growing crack tip. At steady-state the residual stresses far behind the crack tip are tensile near the crack, decreasing to zero at a certain distance from the crack plane. It is shown that the shielding due to pre-existing dislocations, e.g., for cleavage in α-iron crystals may be considerable. This revised version was published online in July 2006 with corrections to the Cover Date.     

CREEP CRACK GROWTH MICROMECHANISMS AT 405°C IN COLD BENT TUBES OF A C-Mn-Mo STEEL

Abstract— Creep Crack Growth (CCG) tests were performed at 405°C on specimens cut out of the cold bent extrados of five tubes of a C-Mn-Mo steel. Intergranular fracture and grain boundary cavitation was less in the C-Mn-Mo than in the C-Mn steels, in accordance with better CCG resistance of the former material. The dimensions and hardness variation across the crack tip process zone were measured by microhardness profiles performed on metallographic sections of the broken samples. TEM analysis of the dislocation patterns close to the fracture surface confirmed the presence of temperature- and stress-induced plasticity phenomena. A significant enrichment of N at grain boundaries (GB) inside the process zone was detected by Auger spectroscopy; N not only inhibits dislocation motion and stress field relaxation at the crack tip but also causes a decrease in GB cohesion ahead of the crack tip. These results help in understanding the micromechanisms which reduce the creep ductility of C-Mn-Mo and C-Mn cold bent tubes and the role of chemical composition in improving CCG resistance.     

A Griffith crack shielded by a dislocation pile-up

The dislocation free zone at the tip of a mode III shear crack is analyzed. A pile-up of screw dislocations parallel to the crack front, in anti-plane shear, in the stress field of a crack has been solved using a continuous distribution of dislocations. The crack tip remains sharp and is assumed to satisfy Griffith's fracture criteria using the local crack tip stress intensity factor. The dislocation pile-up shield the sharp crack tip from the applied stress intensity factor by simple addition of each dislocation's negative contribution to the applied stress intensity value. The analysis differs substantially from the well known BCS theory in that the local crack tip fracture criteria enters into the dislocation distributions found.     

The brittle-ductile transition in silicon: The influence of pre-existing dislocation arrangements

It has been shown that the value of the brittle-ductile transition temperature in specimens of silicon containing pre-cracks at the surface is strongly dependent on the pre-existing dislocation arrangement close to the crack tip. In particular, removing dislocations from the vicinity of the crack tip has the effect of raising the transition temperature, while introducing more surface dislocations by grinding reduces the transition temperature. In both cases the transition temperature remains sharp, implying that the crack tip sources have to be nucleated in the test before effective shielding occurs. The differences in the transition temperatures reflect the differences in the distances of the dislocations from the crack tips, and in the source lengths of the pre-existing sources which send dislocations to the crack tip.Specimens which are pre-stressed at the brittle-ductile transition to a K value at which crack tip sources are expected to be nucleated exhibit a gradual “soft” transition over a wide range of temperature. This result is in agreement with the predictions of the computer modelling of dislocation emission from crack tips, developed by Hirsch et.al [3,7], and together with the results on the dependence of T_c on dislocation arrangements, provides strong support for their theory of the brittle-ductile transition.     

An energy analysis of triple junction crack nucleation due to the wedging action of grain boundary dislocations

A physical model of wedge crack formation at triple junctions in polycrystalline materials is analyzed in this paper. The origin of the crack formation is the sliding of grain boundaries meeting at triple junctions. Based on the dislocation model of grain boundaries, the sliding is attributed to the gliding of grain boundary dislocations (GBDs). Consequently, the resulting crack formation can be analyzed theoretically in terms of the energetics of the piling up of interacting GBDs. The model permits the determination of crack stability or instability as well as the length of the stable crack. Results are obtained in this paper for polycrystalline ice and aluminium.     

Dislocation Behaviour near the Crack Tip in Bulk Fe-3% Si Single Crystal

An investigation has been made of the disloca-tion distribution and dislocation free zone near thecrack tip in bulk Fe-3% Si single crystal duringdeformation in SEM.It has been found that anumber of dislocations were emitted from the cracktip during deformation.After that,the dislocationsmoved rapidly away from the crack tip,which indi-cated that they were strongly repelled by the stressfield at the crack tip.Between the crack tip and theplastic zone there is a region of dislocation-free,which is referred to as dislocation-free zone (DFZ).The length of DFZs is roughly estimated 100μm which is much longer than that found in thinfoil specimen.The variation of dislocation densityas a function of the distance from the crack tip wasmeasured,which showed that the dislocations areinversely piled up in the plastic zone.The length ofDFZs increased with both the length of pre-crackand the amplitude of applied stress.     

Fracture toughness of layered structures: Embrittlement due to confinement of plasticity

The fracture toughness of a layered composite material is analyzed employing a combined two dimensional dislocation dynamics (DD)-cohesive zone (CZ) model. The fracture mechanism of an elastic-plastic (ductile) material sandwiched within purely elastic layers approaches ideally brittle behaviour with decreasing layer thickness. We investigate the influence of different constitutive parameters concerning dislocation plasticity as well as the effect of cohesive strength of the ductile material on the scaling of fracture toughness with layer thickness. For a constant layer thickness, the results of the numerical model are consistent with the expectation that fracture toughness decreases with increasing yield strength, but increases with the cohesive strength of the material. The scaling behaviour of the fracture toughness with layer thickness depends on these material parameters, but also on the dislocation microstructure in the vicinity of the crack tip. Strain localization due to easy dislocation generation right at the crack tip improves toughness in thin layers and leads to a jump-like increase of fracture toughness with layer thickness. However, the fracture toughness for films that are thick enough to exhibit bulk behaviour proves to be higher when the distribution of dislocations is more homogeneous, because in this case the crack grows in a stable fashion over some distance.     

Inclined pileup of screw dislocations at the crack tip with a dislocation-free zone

A single pileup of screw dislocations extending from the crack tip along an inclined direction has been observed in experiments. It is often associated with dislocation emission mechanisms at the crack tip. This linear pileup is a microplastic slipline emanating from the crack tip. A region near the crack tip is often free from dislocations because of a finite resistance value for the crack tip to emit dislocations. The mathematical problem is solved in this paper by applying the extended Wiener-Hopf method. The condition of finite stress at the end of the plastic zone, the crack opening displacement, and the stress distribution along the slipline are obtained in analytical expressions. Numerical values are calculated and the results can be used to discuss brittle versus ductile fracture for metals as treated in previous studies. A method to approximately calculate the corresponding results for edge dislocations is suggested.     

Elastic interaction between screw dislocations and the internal crack near a free surface

The elastic interaction between screw dislocation and the internal crack near a free surface has been investigated. The stress intensity factor at the crack tip, crack extension force, the image force on the dislocation are affected by the free surface. The number and nature of dislocations, m, inside the crack also play an important role in fracture. In order to understand the plastic zone, the zero-force points of dislocation along the x-axis are involved. The dislocation emitted from the right-hand crack tip is enhanced by positive m and reduced by negative m. On the other hand, if the internal crack is closer to the free surface, a dislocation generated from the right-hand crack tip is easier for negative m and more difficult for positive m. However, the role of m on the dislocation emission for the left-hand crack tip is opposite to that for the right-hand crack tip. Finally, three special cases can be obtained from our results. (1) The interaction between a dislocation and a surface crack; (2) the interaction between a dislocation and an internal crack; (3) the interaction between two dislocations.     

Transient analysis of multiple parallel cracks under anti-plane dynamic loading

The problem of a homogeneous linear elastic body containing multiple non-collinear cracks under anti-plane dynamic loading is considered in this work. The cracks are simulated by distributions of dislocations and an integral equation relating tractions on the crack planes and the dislocation densities is derived. The integral equation in the Laplace transform domain is solved by the Gaussian–Chebyshev integration quadrature. The dynamic stress intensity factor associated with each crack tip is calculated by a numerical inverse Laplace scheme. Numerical results are given for one crack and two or three parallel cracks under normal incidence of a plane horizontally shear stress wave.