Mixed mode (I and II) crack tip fields in bulk metallic glasses

Stationary crack tip fields in bulk metallic glasses under mixed mode (I and II) loading are studied through detailed finite element simulations assuming plane strain, small scale yielding conditions. The influence of internal friction or pressure sensitivity on the plastic zones, notch deformation, stress and plastic strain fields is examined for different mode mixities. Under mixed mode loading, the notch deforms into a shape such that one part of its surface sharpens while the other part blunts. Increase in mode II component of loading dramatically enhances the normalized plastic zone size, lowers the stresses but significantly elevates the plastic strain levels near the notch tip. Higher internal friction reduces the peak tangential stress but increases the plastic strain and stretching near the blunted part of the notch. The simulated shear bands are straight and extend over a long distance ahead of the notch tip under mode II dominant loading. The possible variations of fracture toughness with mode mixity corresponding to failure by brittle micro-cracking and ductile shear banding are predicted employing two simple fracture criteria. The salient results from finite element simulations are validated by comparison with those from mixed mode (I and II) fracture experiments on a Zr-based bulk metallic glass.     

Higher order fields for damaged nonlinear antiplane shear notch, crack and inclusion problems

Damaged nonlinear antiplane shear problems with a variety of singularities are studied analytically. A deformation plasticity theory coupled with damage is employed in analysis. The effect of microscopic damage is considered in terms of continuum damage mechanics approach. An exact solution for the general damaged nonlinear singular antiplane shear problem is derived in the stress plane by means of a hodograph transformation, then corresponding higher order asymptotic solutions are obtained by reversing the stress plane solution to the physical plane. As example, traction free sharp notch and crack, rigid sharp wedge and flat inclusion, and mixed boundary sharp notch problems are investigated, respectively. Consequently, higher order fields are obtained, in which analytical expressions of the dominant and second order singularity exponents and angular distribution functions of the near tip fields are derived. Effects of the damage and hardening exponents of materials and the geometric angle of notch/wedge on the near tip quantities are discussed in detail. It is found that damage leads to a weaker dominant singularity of stress, but to little stronger singularities of the dominant and second order terms of strain compared to that for undamaged material. It is also seen that damage has important effect on the angular distribution functions of the near tip stress and strain fields. As special cases, higher order analytical solutions of the crack and rigid flat inclusion tip fields are obtained, respectively, by reducing the notch/wedge tip solutions. Effects of damage and hardening exponents on the dominant and second order terms in the solutions of the crack and inclusion tip fields are discussed.     

Embedded optical fiber Bragg grating sensor in a nonuniform strain field: Measurements and simulations

This paper investigates the use of embedded optical fiber Bragg gratings to measure strain near a stress concentration within a solid structure. Due to the nature of a stress concentration (i.e., the strong nonuniformity of the strain field), the assumption that the grating spectrum in reflection remains a single peak with a constant bandwidth is not valid. Compact tension specimens including a controlled notch shape are fabricated, and optical fiber Bragg gratings with different gage lengths are embedded near the notch tip. The form of the spectra in transmission varies between gages that are at different distances from the notch tip under given loading conditions. This variation is shown to be due to the difference in the distribution of strain along the gage length. By using the strain field measured using electronic speckle pattern interferometry on the specimen surface and a discretized model of the grating, the spectra in transmission are then calculated analytically. For a known strain distribution, it is then shown that one can determine the magnitude of the applied force on the specimen. Thus, by considering the nonuniformity of the strain field, the optical fiber Bragg gage functions well as an embedded strain gage near the stress concentration.     

Modeling of fatigue crack growth from a notch

When a fatigue crack is nucleated and propagates into the vicinity of the notch, the crack growth rate is generally higher than that can be expected by using the stress intensity factor concept. The current study attempted to describe the crack growth at notches quantitatively with a detailed consideration of the cyclic plasticity of the material. An elastic–plastic finite element analysis was conducted to obtain the stress and strain histories of the notched component. A single multiaxial fatigue criterion was used to determine the crack initiation from the notch and the subsequent crack growth. Round compact specimens made of 1070 steel were subjected to Mode I cyclic loading with different R-ratios at room temperature. The approach developed was able to quantitatively capture the crack growth behavior near the notch. When the R-ratio was positive, the crack growth near a notch was mainly influenced by the plasticity created by the notch and the resulted fatigue damage during crack initiation. When the R-ratio was negative, the contact of the cracked surfaces during a part of a loading cycle reduced the cyclic plasticity of the material near the crack tip. The combined effect of notch plasticity and possible contact of cracked surface were responsible for the observed crack growth phenomenon near a notch.     

Large strain field near a notch-tip under tension

This paper considers the plane strain notch-tip filed in a rubber-like material under tension. Based on a new constitutive relation, the asymptotic equations of the near tip field are derived and solved. It is shown that the notch tip field is composed of two narrowing sectors and an expanding sector. The dominant stress and strain near the notch tip are found to be in a state of uniaxial tension.     

环形含内钝缺口试件在斜加载作用下弹塑性场的分布特征

运用有限元数值计算方法,对含内缺口的环形试件在受到与缺口对称面成不同角度的载荷作用下,缺口前端附近的弹塑性场分布特征进行了研究,得到了缺口前端应力和应变场以及塑性区的变化情况,此外,选取缺口弧度不同的模型,探讨了缺口顶端钝化程度对弹生场的影响,为了对计算结果的可靠性进行评价,运用高精度云纹干涉实验力学方法进行了验证性试验,试验得到的位移,应变分布情况与计算结果基本吻合。文中得到的结果为钝口试件的材     

Reliability assessment of a bi-material notch: Strain energy density factor approach

Joints of different materials have many applications in structural engineering and microelectronics. In the present contribution the joint is modelled as a bi-material notch. The singular stress field near the notch tip is investigated. Depending on the notch geometry and materials, the stress field can have one or two singularities. It is shown that to study the problem of a crack onset at the notch, both terms have to be taken into account. Criteria for the direction and for crack nucleation are formulated. The approach utilizes the knowledge of the strain energy density factor distribution in a bi-material notch vicinity.     

橡胶楔体与刚性缺口接触大变形分析

利用Ｋｕｏｗｌｅｓ与Ｓｔｅｒｎｂｅｒｇ提出的非线性弹性大变形应变能函数,对橡胶楔体与刚性缺口接触问题进行大变形渐近分析,推导了楔体尖端场的渐近方程,得到楔体尖端附近的应力应变场及应力的奇异性指数与橡胶楔体角度、刚性缺口角度及材料常数有关的表达式;楔尖附近同一半径上应力分量为常数,同时,利用非线性有限元理论编制了大变形有限元程序,考虑楔体尖端与缺口接触边界条件,计算得到了与分析解一致的结论,当缺口角     

三点弯曲试样动态应力强度因子计算研究

利用Hopkinson压杆对三点弯曲试样进行冲击加载,采集了垂直裂纹面距裂尖2mm和与裂纹面成60°距裂尖5mm处的应变信号。根据裂尖附近测试的应变信号计算试样的动态应力强度因子,并与有限元计算结果进行比较,结果表明由于裂尖有一段疲劳裂纹区,通过裂尖附近应变信号来计算动态应力强度因子时,如果裂尖位置确定不准及粘贴应变片位置不够准确对计算结果将带来很大影响。因此利用应变片法计算动态应力强度因子时,为了获得更准确的计算结果,在实验后应对试件裂纹面进行分析测量,重新确定裂尖位置,必要时需对应变片至裂尖距离进行修正后再计算动态应力强度因子值。     

橡胶楔体与刚性缺口接触的大变形分析

利用一种新的橡胶材料应变能函数，对橡胶楔体与刚性缺口接触大变形问题进行了分析。得到了接触尖点附近变形的奇异性特征，给出了奇异性指数与材料常数、橡胶楔体角度、刚性缺口角度之间的关系式。同时编制了大变形有限元程序，计算得到了与理论解一致的结论。     

反平面V形切口塑性应力奇异性分析

论文提出了用插值矩阵法计算幂硬化塑性材料反平面V形切口和裂纹尖端区域的应力奇异性.首先在切口和裂纹尖端区域采用自尖端径向度量的渐近位移场假设,将其代入塑性全量理论的基本微分方程后,推导出包含应力奇异性特征指数和特征角函数的非线性常微分方程特征值问题.然后采用插值矩阵法迭代求解导出的控制方程,得到一般的塑性材料反平面V形切口和裂纹的前若干阶应力奇异阶和相应的特征角函数,该法的重要优点是以上求解的特征角函数和它们各阶导函数具有同阶精度,并且一次性地求出前若干阶特征对.同时,插值矩阵法计算量小,易于和其他方法联合使用,这些优点在后续求解尖端区域完全应力场非常优越.论文方法的计算结果与现有结果对照,发现吻合良好,表明了论文方法的有效性.     

Finite deformation analysis of crack tip fields in plastically compressible hardening–softening–hardening solids

Crack tip fields are calculated under plane strain small scale yielding conditions. The material is characterized by a finite strain elastic–viscoplastic constitutive relation with various hardening–softening–hardening hardness functions. Both plastically compressible and plastically incompressible solids are considered. Displacements corresponding to the isotropic linear elastic mode I crack field are prescribed on a remote boundary. The initial crack is taken to be a semi-circular notch and symmetry about the crack plane is imposed. Plastic compressibility is found to give an increased crack opening displacement for a given value of the applied loading. The plastic zone size and shape are found to depend on the plastic compressibility, but not much on whether material softening occurs near the crack tip.On the other hand, the near crack tip stress and deformation fields depend sensitively on whether or not material softening occurs. The combination of plastic compressibility and softening(or softening–hardening) has a particularly strong effect on the near crack tip stress and deformation fields.     

A MOVING CRACK IN A NONHOMOGENEOUS MATERIAL STRIP

This paper considers an anti-plane moving crack in a nonhomogeneous material strip of finite thickness. The shear modulus and the mass density of the strip are considered for a class of functional forms for which the equilibrium equation has analytical solutions. The problem is solved by means of the singular integral equation technique. The stress field near the crack tip is obtained. The results are plotted to show the effect of the material non-homogeneity and crack moving velocity on the crack tip field. Crack bifurcation behaviour is also discussed. The paper points out that use of an appropriate fracture criterion is essential for studying the stability of a moving crack in nonhomogeneous materials. The prediction whether the unstable crack growth will be enhanced or retarded is strongly dependent on the type of the fracture criterion used. Based on the analysis, it seems that the maximum 'anti-plane shear' stress around the crack tip is a suitable failure criterion for moving cracks in nonhomogeneous materials.     

The expression of stress and strain at the tip of three-dimensional notch

ＴＨＥＥＸＰＲＥＳＳＩＯＮＯＦＳＴＲＥＳＳＡＮＤＳＴＲＡＩＮＡＴＴＨＥＴＩＰＯＦＴＨＲＥＥ－ＤＩＭＥＮＳＩＯＮＡＬＮＯＴＣＨＱｉａｎＪｕｎ（钱俊）（ＥａｓｔＣｈｉｎａＩｎｓｔｉｔｕｔｅｏｆＴｅｃｈｍｏｌｏｇｙ）ＮａｎｊｉｎｇＬｏｎｇＹｕ－ｑｉｕ（龙驭...     

Logarithm strain singularity at tip of mode I growing crack in elasticand perfectly-plastic material

Reanalyzed in detail is the stress and strain distribution near the tip of a Mode I steadily growing crack in an elastic and perfectly-plastic material. The crack tip region is divided into five angular sectors, one of which is singular in character and represents a rapid transition zone that becomes a line of strain discontinuity in the limit as crack tip is approached. It is shown for an incompressible material (ν=0.5) under plane strain that the local strain in all the angular sectors possesses the same logarithm singularity, i.e., In r where r is the radial distance measured from the crack tip. This result also prevails for the compressible material ( v < 0.5) and resolves a long standing controversy concerning the strain singularity in the sector just ahead of the crack tip.     

Finite strain analysis of crack tip fields in incompressible hyperelastic solids loaded in plane stress

A new method is developed to determine the dominant asymptotic stress and deformation fields near the tip of a Mode-I traction free plane stress crack. The analysis is based on the fully nonlinear equilibrium theory of incompressible hyperelastic solids. We show that the dominant singularity of the near tip stress field is governed by the asymptotic solution of a linear second order ordinary differential equation. Our method is applicable to any hyperelastic material with a smooth work function that depends only on the trace of the Cauchy-Green tensor and is particularly useful for materials that exhibit severe strain hardening. We apply this method to study two types of soft materials: generalized neo-Hookean solids and a solid that hardens exponentially. For the generalized neo-Hookean solids, our method is able to resolve a difficulty in the previous work by Geubelle and Knauss (1994a). Our theoretical results are compared with finite element simulations.     

A new version of the Neuber rule accounting for the influence of the notch opening angle for out-of-plane shear loads

The paper deals with nonlinear stress and strain distributions at the root of sharp and rounded notches with different opening angles under antiplane shear loading and small scale yielding. In order to make an easier comparison with the Neuber rule, the material is thought of as obeying the particular nonlinear law used in the past just by Neuber.By solving the linear differential equation resulting from the use of the hodograph transformation, a new relationship linking linear and nonlinear stress and strain concentrations is found. The relationship is written also in terms of the relevant notch stress intensity factors. In contrast with the Neuber rule, this relationship strictly depends on the notch opening angle. Even when the notch opening angle is zero, it does not match the Neuber Rule, but results in an additional factor 2 which is in agreement with Hult and McClintock’s solution when the notch tip radius tends to zero and the notch becomes a crack.     

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

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

Strain and failure analysis of graphite/epoxy plates with cracks

The objective of this investigation was to study the deformation and failure of uniaxially loaded graphite/epoxy plates with cracks and to determine the influence of notch size on failure. The specimens were quasi-isotropic laminates with cracks of various lengths. They were instrumented with strain gages, photoelastic coatings and moiré grids. Strains near the crack tip show two distinct points of rate change at strain levels of 0.002 and 0.006, the latter corresponding to the ultimate strain of the 90-deg plies. Failure near the crack tip takes the form of a damage zone consisting of ply subcracking along fibers, local delamination and fiber breakage. Failure occurs when this damage zone reaches some critical value. Measured maximum strains at failure exceeded twice the ultimate strain of the unnotched laminate. The average stress over a characteristic distance (5 mm) from the crack tip was used as a criterion to describe the influence of crack length on failure. Comparison of results with those from similar specimens with circular holes showed that strength was nearly independent of notch geometry in this case, i.e., specimens with holes and cracks of the same size had nearly the same strength.     

Study on fracture behavior of ferroelectric ceramics under combined electromechanical loading by using a moiré interferometry technique

This paper discusses an in situ observation of fracture behavior around a crack tip in ferroelectric ceramics under combined electromechanical loading by use of a moiré interferometry technique. The deformation field induced by the electric field and the stress concentration near the crack tip in three-points bending experiments was measured. By analysis of the moiré images it is found that under a constant mechanical load, the electric field almost has no effect on the crack extension in the case that the directions of the poling, electric field and crack extension are perpendicular to each other. When the poling direction is parallel to the crack extension direction and perpendicular to the electric field, the strain decreases faster than that calculated by FEM with and without electrical loading as one goes away from the crack tip. In addition, as the electric field intensity increases, the strain near the crack tip increases, and the strain concentration becomes more significant. The project supported by the National Natural Science Foundation of China (10132010, 10025209, 10232023)