二维孔隙裂隙双重介质逾渗规律研究

在孔隙介质逾渗理论的基础上，将另外一个非常重要的渗透通道——裂隙引入到介质的逾渗研究中，提出了更为普遍的孔隙裂隙双重介质的逾渗研究方法.通过对二维平面孔隙裂隙双重介质的数值计算，得到了孔隙裂隙双重介质三个重要参数：孔隙率，裂隙分形维数、裂隙数量分布初值与逾渗概率的关系，给出了孔隙裂隙双重介质逾渗阈值的数学表达式，揭示了孔隙裂隙双重介质的逾渗规律. 关键词： 孔隙 裂隙 双重介质 逾渗 逾渗阈值     

Numerical simulation of deformation and fracture of a material with a polysilazane-based coating

The paper studies the localization of plastic deformation and fracture in a material with a porous coating. A dynamic boundary value problem in the plane strain formulation is solved. The numerical simulation is performed by the finite difference method. The composite structure corresponds to the experimentally observed one and is specified explicitly in the calculation. A generation procedure of the initial finite-difference grid is developed to describe the coating structure with adjustable porosity and geometry of the substrate-coating interface. Constitutive equations for the steel substrate include an elastic-plastic model of an isotropically hardening material. The ceramic coating is described by a brittle fracture model on the basis of the Huber criterion which accounts for crack nucleation in triaxial tension zones. It is shown that the specific character of deformation and fracture of the studied composite results from the presence of local tensile regions in the vicinity of pores and along the coating-substrate interface, in both tension and compression of the coated material. The interrelation between inhomogeneous plastic flow in the steel substrate and crack propagation in the coating is studied.     

各向异性渗流条件下弹性波的传播特征

研究了弹性波在非均匀裂纹孔隙介质中的传播特性,建立了各向异性喷射流模型.当弹性波通过裂纹孔隙介质时,由于波的扰动及裂纹和孔隙几何结构的不一致,导致在裂纹内部及裂纹与周边孔隙之间同时存在着流体压力梯度.此时的弹性波波动响应中包含着裂纹内连通性特征和背景孔隙渗透率信息.流体的动态流动过程使得介质的等效弹性参数为复数(非完全弹性),并且具有频率依赖性.当弹性波为低频和高频极限时,介质为完全弹性;当处于中间频段时,波有衰减和频率依赖.裂纹孔隙介质的各向异性连通性(渗透率)对应着各向异性特征频率(当渗流长度等于非均匀尺度时的弹性波频率),波的传播受到裂纹内连通性的影响.在一定频段内,随着裂纹厚度的增加,将出现第二峰值,峰值大小同时受到裂纹厚度和半径的影响.     

Non-equilibrium statistical theory about microscopic fatigue cracks of metal in magnetic field

This paper develops the non-equilibrium statistical fatigue damage theory to study the statistical behaviour of micro-crack for metals in magnetic field. The one-dimensional homogeneous crack system is chosen for study. To investigate the effect caused by magnetic field on the statistical distribution of micro-crack in the system, the theoretical analysis on microcrack evolution equation, the average length of micro-crack, density distribution function of micro-crack and fatigue fracture probability have been performed. The derived results relate the changes of some quantities, such as average length, density distribution function and fatigue fracture probability, to the applied magnetic field, the magnetic and mechanical properties of metals. It gives a theoretical explanation on the change of fatigue damage due to magnetic fields observed by experiments, and presents an analytic approach on studying the fatigue damage of metal in magnetic field.     

In situ EBSD during tensile test of aluminum AA3003 sheet

Miniature tensile-test specimens of soft-annealed, weakly textured AA3003 aluminum sheet in 0.9 mm thickness were deformed until fracture inside a scanning electron microscope. Tensile strength measured by the miniature tensile test stage agreed well with the tensile strength by regular tensile testing. Strain over the microscope field of view was determined from changes in positions of constituent particles. Slip lines were visible in secondary electron images already at 0.3% strain; activity from secondary slip systems became apparent at 2% strain. Orientation rotation behavior of the tensile load axis with respect to the crystallographic axes agreed well with previously reported trends for other aluminum alloys. Start of the fracture and tensile crack propagation were documented in secondary electron images. The region of fracture nucleation included and was surrounded by many grains that possessed high Schmid factors at zero strain. Crystal lattice rotation angles in the grains surrounding the initial fracture zone were higher than average while rotations inside the initial fracture zone were lower than average for strains from zero to 31%. The orientation rotation behavior of the tensile load axes of the grains around the fracture zone deviated from the average behavior in this material.     

Atomistic scale fracture behaviours in hierarchically nanotwinned metals

In the present study, a series of large-scale molecular dynamics simulations have been performed to investigate the atomistic scale fracture behaviours along the boundaries of primary twins in Cu with hierarchically nanotwinned structures (HTS), and compare their fracture behaviours with those in monolithic twins. The results indicate that crack propagation along [1?1?2] on the twin plane in monolithic nanotwins is brittle cleavage and fracture, resulting in low crack resistance and fracture toughness. However, the crack resistance along the boundaries of primary twins in HTS is much higher, and a smaller spacing of secondary twins (λ ₂) leads to even higher fracture toughness. With large λ ₂, the crack growth is achieved by void nucleation, growth and coalescence. However, considerable plastic deformation and enhanced fracture toughness in HTS could be achieved by the crack blunting and by the extensive dislocation accommodation ahead of the crack tip when λ ₂ is small.     

Structural and phase changes in metals during formation of the fatigue fracture surface in metals and alloys

The mechanisms of deformation-induced structural and phase changes, which accompany crack nucleation and growth during fatigue fracture of metals and alloys, are considered. It is shown that the crack steady growth phase is determined by the processes of local melting of the material in extension phase, and melt crystallization processes with formation of shrinkage pores in front of the crack tip in the unloading phases and crack narrowing.     

Dislocation Plasticity Effects on Interfacial Fracture

Analyses are reviewed where plastic flow in the vicinity of an interfacial crack is represented in terms of the nucleation and glide of discrete dislocations. Attention is confined to cracks along a metal-ceramic interface, with the ceramic idealized as being rigid. Both monotonic and fatigue loading are considered. The main focus is on the stress and deformation fields near the crack tip predicted by discrete dislocation plasticity, in comparison with those obtained from conventional continuum plasticity theory. The role that discrete dislocation plasticity can play in interpreting interface fracture properties in the presence of plastic flow is discussed.     

Thermal shock fracture of piezoelectric materials

This paper deals with a case study for the piezoelectric materials suddenly exposed to an environmental medium of different temperature. The problem is idealized to a plate containing an edge crack or an embedded crack. The stress and electric displacement histories in an uncracked plate are calculated. These stresses and electric displacements are then added to the crack surface tractions and electric displacements with opposite sign to formulate a mixed boundary value problem. The cracking problem is thus reduced to a singular integral equation of Cauchy type, which is then solved numerically. Both impermeable crack assumption and permeable crack assumption are considered. The results for stress and electric displacement intensity factors are computed as a function of normalized time and crack size. Lower bound solutions are obtained for the maximum thermal shock that the material can sustain without catastrophic failure according to the two distinct criteria: (i) The maximum local tensile stress equals the tensile strength of the medium. (ii) The maximum stress intensity factor for the pre-existing representative crack equals the fracture toughness of the medium. The parameters that control the transient thermal stress and electric displacement are also identified. The method can be used to explore susceptibility to thermal fracture in piezoelectric materials containing pre-cracks.     

Crack growth on the basal plane in single crystal zinc: experiments and computations

Crack propagation on the basal planes in zinc was examined by means of in situ fracture testing of pre-cracked single crystals, with specific attention paid to the fracture mechanism. During quasistatic loading, crack propagation occurred in short bursts of dynamic crack extension followed by periods of arrests, the latter accompanied by plastic deformation and blunting of the crack-tip. In situ observations confirmed nucleation and propagation of microcracks on parallel basal planes and plastic deformation and failure of the linking ligaments. Pre-existing twins in the crack path serve as potent crack arrestors. The crystallographic orientation of the crack growth direction on the basal plane was found to influence both the fracture load as well as the deformation at the crack-tip, producing fracture surfaces of noticeably different appearances. Finite element analysis incorporating crystal plasticity was used to identify dominant slip systems and the stress distribution around the crack-tip in plane stress and plane strain. The computational results are helpful in rationalizing the experimental observations including the mechanism of crack propagation, the orientation dependence of crack-tip plasticity and the fracture surface morphology.     

Brittle and ductile fracture of semiconductor nanowires – molecular dynamics simulations

Fracture of silicon and germanium nanowires in tension at room temperature is studied by molecular dynamics simulations using several interatomic potential models. While some potentials predict brittle fracture initiated by crack nucleation from the surface, most potentials predict ductile fracture initiated by dislocation nucleation and slip. A simple parameter based on the ratio between the ideal tensile strength and the ideal shear strength is found to correlate very well with the observed brittle versus ductile behaviours for all the potentials used in this study. This parameter is then computed by ab initio methods, which predict brittle fracture at room temperature. A brittle-to-ductile transition (BDT) is observed in MD simulations at higher temperature. The BDT mechanism in semiconductor nanowires is different from that in the bulk, due to the lack of a pre-existing macrocrack that is always assumed in bulk BDT models.     

Some aspects of strength in sea ice mechanics

The paper provides estimates of the relation between the strength of sea ice and its porosity in the framework of a porous solid model proposed earlier. The influence of the ice capillary structure filled with brine on the ice strength is analyzed, showing that the rise of local stresses in the ice changes its limit equilibrium conditions. Comparison is made to elucidate how the sea ice strength responds to loading length- and crosswise the capillary axes. Patterns by which ice structure elements merge resulting mostly in main cracks or in multiple ordered crack systems under compression are discussed. The tendency toward generation of ordered crack systems or propagation of single main cracks in the structured material is determined by the activity of its structural elements involved in fracture. A method is proposed for estimating the average strength of sea ice cover relative to horizontal compression with regard to partial loss of the bearing capacity of ice layers during deformation.     

Effect of inclusions on heterogeneous crack nucleation in nanocomposites

A two-dimensional theoretical model is proposed for the heterogeneous nucleation of a grain-boundary nanocrack in a nanocomposite consisting of a nanocrystalline matrix and nanoinclusions whose elastic moduli are identical to those of the matrix. The inclusions have the form of rods with a rectangular cross section and undergo dilatation eigenstrain induced by the differences in the lattice parameters and thermal expansion coefficients of the matrix and inclusions. In terms of the model, a mode-I–II nanocrack nucleates at the negative disclination of a biaxial dipole consisting of wedge grain-boundary (or junction) disclinations; then, the nanocrack opens along a grain boundary and reaches an inclusion boundary. Depending on the relative positions and orientations of the initial segment of the nanocrack and the inclusion, the nanocrack can either penetrate into the inclusion or bypass it along the matrix-inclusion interface. The nanocrack nucleation probability increases near an inclusion with negative (compressive) dilatation eigenstrain. A decrease in the inclusion size decreases (increases) the probability of a crack opening along the interface if the dilatation eigenstrain is negative (positive).     

Influence of nanotwin generation near crack twins on the fracture toughness of nanomaterials

A theoretical model of microscopic mechanisms of the nucleation and development of deformation twins in nanocrystalline materials has been developed. Within the model, we have studied the generation of deformation twins near crack tips, which occurs through multiple nanoscopic shears that represent nanoscopic regions of an ideal plastic shear. It has been shown that the nucleation of such nanotwins near crack tips reduces the high local stresses that arise near these tips. Thus, the generation and development of nanotwins near crack tips increases the fracture toughness of brittle nanocrystalline materials and serves as an efficient mechanism of improving the crack resistance of deformed nanocrystalline materials.     

Some aspects of nucleation and evolution of microscopic and mesoscopic cracks and quasi-brittle fracture of homogeneous materials

The dynamics of nucleation of microcracks in the region of a developing macroscopic crack has been studied using scanning electron microscopy during in situ experiment and the acoustic emission method. An explosive-like nucleation of microcracks and the influence of dissipative properties of a material on the size and the rate of redistribution of local stresses of the microcracks have been established. Each of nucleations of microscopic and mesoscopic defects is considered as an act of local testing of the dissipative ability of the system, and the interaction between microcracks is determined by relaxation processes when the ability is sufficient. With deteriorating dissipative properties (with increasing residence time under loading, deformation rate, etc.), an elastic linear interaction becomes possible, and the macroscopic fracture occurs. In the microcrack dynamics, a specific ductile-brittle transition is determined by a time deterioration of the dissipative properties of the material under the action of a mechanical load. Generally, as a large latent energy is stored during deformation, a macroscopic fracture can be caused by any little discrete structural transformation.     

Experimental and theoretical study of multiscale damage-failure transition in very high cycle fatigue

Multiscale mechanisms of failure of metals (Armco iron, titanium, aluminum) are studied for high cycle and very high cycle fatigue. By correlating with the results of structural studies, a theoretical approach is developed to describe fatigue crack kinetics in damaged material under high cycle and very high cycle fatigue loading conditions. Stages of crack nucleation and propagation are analyzed using the profilometry data from the fracture surface. The scale invariance of fracture surface roughness is established, which allows an explanation of the self-similar nature of fatigue crack kinetics under high cycle and very high cycle fatigue. Variation of elastic-plastic properties of Armco iron under very high cycle fatigue is studied using an acoustic resonance method. It is found that the material density decreases during fatigue damage accumulation, with the minimum of the material density in the bulk of the specimen.     

Influence of the Defect Structure of the Material on the Fracture and Crack Motion

The influence of the structure of the medium on the crack propagation in a Finsler space is examined in the context of the fracture theory. The structure of the medium is taken into account via the connectivity coefficient of the Finsler space and its metric. The deformation of the fibered space caused by the crack motion is analyzed.     

The Uphill Turtle Race; On Short Time Nucleation Probabilities

The short time behavior of nucleation probabilities is studied by representing nucleation as a diffusion process in a potential well with escape over a barrier. If initially all growing nuclei start at the bottom of the well, the first nucleation time on average is larger than the inverse nucleation frequency. Explicit expressions are obtained for the short time probability of first nucleation. For very short times these become independent of the shape of the potential well. They agree well with numerical results from an exact enumeration scheme. For a large number N of growing nuclei the average first nucleation time scales as 1/log N in contrast to the long-time nucleation frequency, which scales as 1/N. For linear potential wells closed form expressions are given for all times.     

Critical conditions for dislocation nucleation at surface steps

Dislocation nucleation at a surface step is analyzed based on a general variational boundary integral formulation of the Peierls–Nabarro dislocation model. By modelling the surface step as part of the surface of a three-dimensional crack, the free surface effect is taken into account by transferring the half space problem into an equivalent one in the infinite medium. The profiles of embryonic dislocations, corresponding to the relative displacements between the two adjacent atomic layers along the slip planes, are then rigorously solved through the variational boundary integral method. The critical conditions for dislocation nucleation are determined by solving the stress-dependent activation energies required to activate the embryonic dislocations from their stable to unstable saddle-point configurations. In particular, the effect of the step geometry, such as the height of the step, the dip angle of the slip plane and the inclined angle of the step surface on dislocation nucleation, is quantitatively ascertained. The results show that the atomic-scale surface step can rapidly reduce the critical stress required for dislocation nucleation from the surface by nearly an order of magnitude. The decrease in critical stress as a function of the height of the step is more significant for slip planes with smaller dip angles and surface steps with smaller inclined angles.     

Dynamic crack nucleation,propagation, and interactions with crystalline secondary phases in aluminum alloys subjected to large deformations

The major objective of this work was to model within a continuum framework the dynamic nucleation and evolution of failure surfaces in aluminum alloys with complex microstructures, using a recently developed compatibility-based fracture criterion for large deformations. Computational analyses were conducted to understand how Mn-bearing dispersoids, Ω and θ′ precipitates affect dynamic fracture processes in an Al–Cu–Mg–Ag alloys (2139-Al). High strain-rate simulations were based on a rate-dependent dislocation-density-based crystalline plasticity formulation and a nonlinear explicit dynamic finite-element approach. Results indicate that the fracture criterion elucidated how dispersoids and precipitates have a dominant role in dynamic crack blunting, branching and arrest. Rationally orientated precipitates result in overall dynamic microstructural strengthening and enhanced uniformity of deformation. These precipitates, however, accelerated unstable crack propagation, and this is amplified in the presence of a pre-crack. In contrast, dispersoids decreased microstructural toughness and ductility, but greatly improved dynamic damage tolerance, especially in the presence of a pre-crack. It can also be predicted that low angle boundaries can change the propagation direction of ductile cracks, and contribute to damage tolerance without crack initiation. Collectively, rationally oriented precipitates and dispersoids can significantly improve the combined dynamic strength, toughness and damage tolerance of crystalline aluminum alloys.