An analysis of ductile rupture modes at a crack tip

An elastic-Viscoplastic model of a ductile, porous solid is used to study the influence of the nucleation and growth of micro-voids in the material near the tip of a crack. Conditions of small scale yielding are assumed, and the numerical analyses of the stress and strain fields are based on finite strain theory, so that crack tip blunting is fully accounted for. An array of large inclusions or inclusion colonies, with a relatively low strength, results in large voids near the crack tip at a rather early stage, whereas small second phase particles in the matrix material between the inclusions require large strains before cavities nucleate. Various distributions of the large inclusions, and various critical strains for nucleation of the small scale voids between the inclusions, are considered. Localization of plastic flow plays an important role in determining the failure path between the crack tip and the nearest larger void, and the path is strongly sensitive to the distribution of the large inclusions. Values of the J-integral and the crack opening displacement at fracture initiation are estimated, together with values of the tearing modulus during crack growth, and these values are related to experimental results.     

混合层中柱状粒子取向分布函数的研究

数值模拟了平面流固混合层中柱状粒子在平面取向的分布函数，提出了描述流场中某点直粒子取向占优程度的参数Dψ的分析，说明在混合层中心区域粒子的取向排列是有序的，在混合层相邻两个大涡之间的区域，是粒子最可几角占优程度最高的地方，而且最可几取向就是沿流场的拉伸方向，说明流场的拉伸作用及方向对粒子的取向排列具有重要影响。与其它模型相比，粒子取向的分布函数模型能定量反映流场大涡对粒子取向控制的强弱程度。     

Finite element analysis of void growth near a blunting crack tip

Large deformation finite element analysis has been used to study the near crack tip growth of long cylindrical holes aligned parallel to the plane of a mode I plane strain crack. The near crack tip stress and deformation fields are analyzed. The results show that the holes are pulled towards the crack tip and change their shape to approximately elliptical with the major axis radial to the crack. They also grow faster directly ahead of the crack than at an angle to the crack plane. Several crack-hole coalescence criteria are discussed and estimates for the conditions for fracture initiation are given and compared with experimental results. The range of estimates now available from finite element calculations coincides quite well with the range of experimental data for materials containing inclusions which are only loosely bonded to the matrix.     

Two mechanisms of ductile fracture: void by void growth versus multiple void interaction

Two distinct mechanisms of crack initiation and advance by void growth have been identified in the literature on the mechanics of ductile fracture. One is the interaction a single void with the crack tip characterizing initiation and the subsequent void by void advance of the tip. This mechanism is represented by the early model of Rice and Johnson and the subsequent more detailed numerical computations of McMeeking and coworkers on a single void interacting with a crack tip. The second mechanism involves the simultaneous interaction of multiple voids on the plane ahead of the crack tip both during initiation and in subsequent crack growth. This mechanism is revealed by models with an embedded fracture process zone, such as those developed by Tvergaard and Hutchinson. While both mechanisms are based on void nucleation, growth and coalescence, the inferences from them with regard to crack growth initiation and growth are quantitatively different. The present paper provides a formulation and numerical analysis of a two-dimensional plane strain model with multiple discrete voids located ahead of a pre-existing crack tip. At initial void volume fractions that are sufficiently low, initiation and growth is approximately represented by the void by void mechanism. At somewhat higher initial void volume fractions, a transition in behavior occurs whereby many voids ahead of the tip grow at comparable rates and their interaction determines initiation toughness and crack growth resistance. The study demonstrates that improvements to be expected in fracture toughness by reducing the population of second phase particles responsible for nucleating voids cannot be understood in terms of trends of one mechanism alone. The transition from one mechanism to the other must be taken into account.     

爆炸聚能作用下混凝土试件劈裂的高速3D DIC实验

为探寻非钻孔条件下露天爆破大块二次破碎形态的控制方法, 应用线性聚能射流对圆柱混凝土模型试件进行侵彻实验, 使用High-speed 3D DIC(高速三维数字图像相关方法)方法分析试件劈裂发展过程的全场三维形变特征。研究结果表明, 数据分析区内劈裂裂纹扩展速度在4个区间内呈阶梯式变化趋势, 峰值速度为235.52 m/s, 平均速度为140.89 m/s;线性聚能射流侵彻对劈裂裂纹扩展有明显导向作用, 应力集中作用使得劈裂裂纹围绕线性射流侵彻对称轴扩展, 扩展方向变化幅度较小;在劈裂裂纹扩展速度突变的3个时刻, 劈裂裂纹路径产生了3处明显拐点, 在拐点处伴随有支裂纹的产生, 支裂纹的扩展距离均未超过5 cm;主应变集中带形状及分布位置决定了裂纹扩展路径及趋势, 拉应变集中先于裂纹出现, 试件呈现准静态劈裂形态, 劈裂面平整度较高。     

Two illustrations of rapid crack growth in a pre-stressed compressible neo-Hookean material

An unbounded isotropic compressible neo-Hookean solid is initially in equilibrium under uniform tensile (possibly large) pre-stress. In one case, plane strain conditions generate slit crack growth at a constant sub-critical rate; in the other, axial symmetry produces penny-shaped crack growth. The procedure of superposing infinitesimal deformations upon those that are large is carried out in terms of tractable exact full-field solutions.These solutions are examined apart from a specific fracture mechanics model, nevertheless, they show that pre-stress induces, in addition to the expected anisotropy, a critical value above which a negative Poisson effect occurs. It is also found that dilatational, rotational and Rayleigh wave speeds decrease, and that the decrease is greater for the plane strain state associated with slit crack growth than for the axially symmetric state of the penny-shaped crack.Dynamic stress intensity factors are also extracted, and found to fall below those for a linear isotropic solid at the same pre-stress and crack growth rate. Moreover, the range of growth rates for sub-critical crack propagation is also decreased.     

晶体取向对单晶体断裂特征影响的模拟分析

结合Ni基单晶合金制三种不同晶体取向的紧凑拉伸试样试验,本文利用考虑有限变形和品格转动效应的晶体滑移有限元程序对单晶体三维断裂特征进行了模拟计算分析,详细考察了裂纹尖端三维应力场特征和断裂特征,结果表明:晶体取向对裂纹尖端应力场有较大影响,但应力沿试样厚度方向明显分成两个部分,在试样心部,应力沿厚度方向变化不大,在试样外表面则明显变化。裂纹尖端张开位移(CTDD)沿厚度方向类似分成两个部分。垂直于滑移面的应力分量致单晶体的准解理断裂,即裂纹的起裂和扩展途径均与该应力分量有关。     

花键轴冷挤压成形后的三维光塑性模拟研究

本文利用光塑性方法,以聚碳酸酯（Ｐｏｌｙｃａｒｂｏｎａｔｅ）为材料研究了矩形花键轴冷挤压成形过程,获得了花键挤压成形的三维塑性应变场分布,提出其变形的应变分布特征,并对典型截面的应变进行了计算分析;文中还通过对聚碳酸酯材料的实验,得到了使花键齿形充满的最佳坯料尺寸,为实际工艺生产优选出合理的工艺参数。     

APPLICATION OF THE DIGITAL MOIRE METHOD IN FRACTURE ANALYSIS OF A CRACKED RUBBER SHEET

This paper deals with the mode I crack problem of a cracked rubber sheet under plane stress condition using the delicate digital moire technique. Through the four step phase-shifting method of automated fringe analysis, the displacement fields in the Cartesian coordinate system are given. By the coordinate-transform equation, the radial and circular displacement distributions in the polar coordinate system are obtained. The displacement isoline distributions and displacement vector distributions near the crack tip are discussed. The strain isoline distributions near the crack tip are also analyzed in this paper. Finally, the distribution rules for the mechanical fields near the crack tip are discussed with the sector division method.     

Size effects on void growth in single crystals with distributed voids

The effect of void size on void growth in single crystals with uniformly distributed cylindrical voids is studied numerically using a finite deformation strain gradient crystal plasticity theory with an intrinsic length parameter. A plane strain cell model is analyzed for a single crystal with three in-plane slip systems. It is observed that small voids allow much larger overall stress levels than larger voids for all the stress triaxialities considered. The amount of void growth is found to be suppressed for smaller voids at low stress triaxialities. Significant differences are observed in the distribution of slips and on the shape of the deformed voids for different void sizes. Furthermore, the orientation of the crystalline lattice is found to have a pronounced effect on the results, especially for the smaller void sizes.     

平面幂律塑性Ⅰ型裂纹尖端应力场全解

在航空航天、船舶、石油管道和核电等领域,服役结构或部件在长期极端条件下运行,不可避免地会产生裂纹,因此,为研究含裂纹结构的准静态断裂行为,必须了解裂纹尖端附近区域的应力应变场特点.对于幂律材料裂纹构元,研究平面应变和平面应力条件下Ⅰ型裂纹尖端应力场的解析分布.基于能量密度等效和量纲分析,推导了能量密度中值点代表性体积单元(representative volume element, RVE)的等效应力解析方程,并定义其为应力因子,进而针对有限平面应变和平面应力紧凑拉伸(compact tension, CT)试样和单边裂纹弯曲(single edge bend, SEB)试样,以应力因子作为应力特征量,并构造用于表征裂尖等效应力等值线的蝶翅轮廓式和扇贝轮廓式三角特殊函数,提出描述幂律塑性条件下平面I型裂纹尖端应力场的半解析模型.该半解析模型形式简单,对CT和SEB试样的裂尖应力场的预测结果与有限元分析的结果比较表明,两者之间均密切吻合,模型公式可直接用于预测Ⅰ型裂纹尖端应力分布,方便于断裂安全评价和理论发展.     

An investigation of hydraulic fracturing. Part I: one-phase flow model

The prediction of the growth of a hydraulic fracture in an oil bearing formation based on the injection rate of fluid is valuable in applications of the waterflood technique in secondary oil recovery. In this paper, the problem of hydraulic fracture growth is studied under the assumption of uniform distribution of pressure in the fracture and unidirectional permeating flow in an infinitely large isothermal linearly elastic porous medium saturated with a one-phase incompressible fluid. The condition of plane strain is imposed in the study. A comparison of the constant fracture toughness criterion based on the asymptotic value for large crack growth with the crack tip ductility criterion for an ideally plastic solid under plane strain and small-scale yielding conditions indicates that the effect of ductility of rock on the crack growth is so small that the steady state value of the energy release rate can be reached within a short period of crack growth. Thus we can employ the constant fracture toughness criterion in our study. The analysis includes the effects of both fracture volume increase and leak-off of fluid from the surface of the fracture. A nonlinear singular integro-differential equation can be formulated for the quasi-static hydraulic fracture growth under a prescribed injection rate. It is solved numerically by a modified fourth order Runge-Kutta method.     

含微裂纹弹性体的应力应变关系

本义建立了考虑裂纹闭合和裂纹表面摩擦影响的含微裂纹弹性体的应力应变关系，给出了柔度张量增量的显式表达式。对于二维平面应力和平面应变状态，给出了等效工程弹性系数。数值计算结果表明，裂纹闭合和裂纹面摩擦对裂纹体的应力应变关系和等效工程弹性系数有重要影响。     

含微裂纹弹性体的应力应变关系

本义建立了考虑裂纹闭合和裂纹表面摩擦影响的含微裂纹弹性体的应力应变关系，给出了柔度张量增量的显式表达式。对于二维平面应力和平面应变状态，给出了等效工程弹性系数。数值计算结果表明，裂纹闭合和裂纹面摩擦对裂纹体的应力应变关系和等效工程弹性系数有重要影响。     

Modeling dynamic brittle behavior of materials with circular flaws or pores

Compressive failure of brittle materials is driven primarily by crack growth from pre-existing flaws in the material. These flaws, such as grain boundaries, pores, preexisting cracks, inclusions and missing grains, are randomly spaced and have a range of possible shapes and sizes. The current work proposes a micromechanics-based model for compressive dynamic failure of brittle materials with circular pore flaws, which incorporates both the number density and the size distribution of flaws. Results show that the distribution of flaw sizes is very important, particularly at moderate strain rate, since analyses based solely on the mean flaw size overpredict strength. Therefore, in order to increase dynamic strength at low to moderate strain rates, it is most effective to control the presence of large flaws. At very high strain rates, however, crack growth is activated even in small flaws and therefore controlling the total number density rather than the size of the flaws is effective for increasing dynamic strength. Finally, the model shows that neglecting very small flaws in the pore population may not have significant effects on the results in many cases, suggesting that the model is a useful tool for identifying a minimum resolution required for experimental characterization of microstructure.     

三相压电复合本构模型中的弧形界面裂纹

深入研究了三相同心圆柱压电复合本构模型中的弧形绝缘界面裂纹问题。采用复势方法获得了该问题的级数形式的解答,并给出了应力、应变、电位移和电场强度等物理量在全场及界面上的分布,同时推导了裂尖处广义强度因子及裂面张开位移和裂面上电势差的表达式。具体计算表明该级数解答收敛迅速,同时显示出第三相混杂区的影响是不能忽略的。由于裂尖处应力奇异性为－1／2,则这种解答不会出现平面应变状态下界面裂纹裂尖处的振荡奇异性,从而不会产生违反物理实际的裂面相互嵌入现象,则该弹性解答也是建立了坚实的物理基础之上。     

Spread of plasticity from stacked stress concentrations

A large elastic solid containing an infinite sequence of slitlike relaxed cracks with a constant distance of vertical separation is considered. The solid is deforming under plane strain shear conditions (mode II). The plastic relaxation around each of the cracks is represented by a suitable distribution of edge dislocations coplanar with the crack itself, the distribution being determined from a singular integral equation. This equation is solved numerically using an expansion of the non-singular part of the kernel in a series of Chebyshev polynomials. Solutions are obtained for the extent of spread of plasticity from each crack and for the associated dislocation distribution as a function of the crack spacing and the applied load. The results are applied to a brief discussion of the fracture process at stress concentrations using the crack opening displacement criterion.     

The influence of particle shape,volume fraction and distribution on post-necking deformation and fracture in uniaxial tension of AA5754 sheet materials

Finite element analysis was performed over a small particle field, edge constraint plane strain post-necking model. The aim is to understand the roles of particle shape, volume fraction and distribution over the post-necking deformation and fracture of AA5754-O sheet materials. For models containing one single particle, the post-necking deformation decreases when the particle varies from circular to elliptical. The inter-particle spacing, the major parameter of distribution to determine whether a pair of particles belongs to a stringer or not, was varied for models with two particles of circular or elliptical shape. The general trend is that the post-necking deformation and fracture strains decrease with decreasing spacing between particles. There is considerable difference in terms of both fracture topographies and strains for models containing 16 particles when distributions varied from random/uniform to stringer distributions. The post-necking deformation and fracture strains monotonically decrease with particle volume fractions for models with 4–64 particles of random or stringer distribution. This indicates that the post-necking behavior for AA5754-O alloys where the matrix material is rather ductile is not solely controlled by a single or pair of particles although they may become initiation places of damage. Multiple damaging sources such as stringers or large particles can act cooperatively and speed up the damaging propagation of the material, and therefore produce small post-necking deformation and early fracture. The center clustering of particles can be beneficial for post-necking behavior and bendability of sheet materials.     

Statistical evolution of small fatigue crack in 1Cr18Ni9Ti weld metal

Statistical evolution of small fatigue crack behavior is analyzed. Three considerations are given. They consist of the total fit, the consistency of statistical parameters with test data, and the practice of commonly used distributions, namely Weibull (two- and three-parameters), normal, lognormal, extreme minimum, extreme maximum, and exponential. The crack growth rate data of 1Cr18Ni9Ti stainless steel-weld metal reveal that fatigue damage is a gradual process involving an initial non-ordered (chaotic) or load history-independent random state to an ordered load history-dependent random state. The seven distributions are examined. Different from existent conclusion, appropriate assumed distribution for characterizing the evolution data is the normal or extreme minimum value model rather than the Weibull two- or three-parameter. The other distributions do not fit the data as they violated two or all three of the mentioned considerations.     

Large deformation field near a crack tip in rubber-like material

This paper deals with strain field near a crack tip in a rubber-like material under plane strain condition. The constitutive relation adopted here is valid for both small and large strain. The asymptotic equations are derived for a shrinking sector and expanding sector. The closed mathematic solution is obtained for the latter while a numerical solution is found for the former. By connecting deformation of the two sectors, the crack tip field character is found.