表面裂纹疲劳扩展的数值模拟(Ⅱ)

根据Paris疲劳裂纹扩展规律,对拉伸和纯弯曲疲劳载荷下表面裂纹扩展进行了数值模拟。数值模型中,用三次样条函数曲线拟合裂纹尖端,在裂纹扩展增量计算中考虑了裂纹闭合影响。裂纹形状演化的模拟结果与Newman和Raju经验公式预测结果进行了比较,表明了所采用的数值模拟方法的实用性。研究发现,裂纹闭合对疲劳裂纹扩展过程中的裂纹形状演化以及裂纹尖端的应力强度因子(SIF)分布都有明显影响。同裂纹形状演化一样,疲劳裂纹扩展过程中裂纹尖端的SIF分布表现出明显的特征。最后,建议了一个简单函数来统一描述表面裂纹尖端的SIF分布。     

Analysis of vertical cracking phenomena in tensile-strained epitaxial film on a substrate: Part I. Mathematical formulation

This paper presents an analysis of a single vertical crack and periodically distributed vertical cracks in an epitaxial film on a semi-infinite substrate where the cracks penetrate into the substrate. The film and substrate materials have different anisotropic elastic constants, necessitating Stroh formalism in the analysis. The misfit strain due to the lattice mismatch between the film and the substrate serves as the driving force for crack formation. The solution for a dislocation in an anisotropic trimaterial is used as a Green function, so that the cracks are modeled as the continuous distributions of dislocations to yield the singular integral equations of Cauchy-type. The Gauss–Chebyshev quadrature formula is adopted to solve the singular integral equations numerically. Energy arguments provide the critical condition for crack formation, at which the cracks are energetically favorable configurations, in terms of the ratio of the penetration depth into the substrate to the film thickness, the ratio of the spacing of the periodic cracks to the film thickness, and the generalized Dundurs parameters between the film and substrate materials.     

巴西圆盘泥岩试件裂纹扩展及应变演化实验研究

为了研究准静态加载条件下岩石试件巴西劈裂裂纹扩展规律,采用MTS试验机进行准静态加载,同时用高速摄像机记录裂纹扩展过程。采用白光数字散斑处理软件对摄像机记录的照片进行处理,得到试件裂纹扩展过程中应变场的演化情况。通过实验和分析可以看出,由于端部效应及加载方式的原因,因此裂纹起裂点在底部加载部位;泥岩试件表面裂纹的平均扩展速度为252m/s;岩石的非均质性即内部微缺陷、微裂纹使得泥岩试样的开裂并不是沿着中心直径方向,而是偏离一定的角度,初始偏离角度约为17°。裂纹扩展过程可以划分为三个阶段:泥岩试件宏观变形阶段(宏观无裂纹)、宏观裂纹稳定扩展阶段、宏观裂纹动态张裂阶段。同时,在裂纹扩展过程中,表面第一主应变场、水平位移场等变化明显,在开裂部位第一主应变最大。通过对圆盘泥岩试件裂纹扩展实验研究,可为研究岩石破裂及其演化规律提供依据。     

爆炸荷载作用下两平行裂纹对扩展中裂纹的影响规律

用实验和数值模拟方法,研究在爆炸载荷下岩体内部一对平行裂纹对扩展主裂纹的影响规律。实验中,采用带有中心装药孔及预制裂纹的砂岩圆盘试件,利用由示波器、超动态应变仪及裂纹扩展计所组成的测试系统,监测主裂纹扩展速度和扩展距离;数值模拟中,采用了AUTODYN软件进行,模拟了主裂纹及两平行裂纹的扩展规律,对岩石材料,采用线性状态方程及最大拉应力失效准则,并在两平行裂纹间设置相应的观测点记录应力曲线。通过实验与数值模拟分析,得到:爆炸载荷下,紧随冲击波后的稀疏波经过两平行裂纹面反射后变成压缩波,并在两平行裂纹间产生垂直于主裂纹扩展方向的压应力,对裂纹的扩展有压制、止裂作用;而且,这种压应力的大小与两平行裂纹的间距有关,进而导致了不同的止裂效果,影响裂纹的扩展速度及最终扩展长度。     

Bending of Timoshenko beam with effect of crack gap based on equivalent spring model

Considering the effect of crack gap,the bending deformation of the Timoshenko beam with switching cracks is studied.To represent a crack with gap as a nonlinear unidirectional rotational spring,the equivalent flexural rigidity of the cracked beam is derived with the generalized Dirac delta function.A closed-form general solution is obtained for bending of a Timoshenko beam with an arbitrary number of switching cracks.Three examples of bending of the Timoshenko beam are presented.The influence of the beam’s slenderness ratio,the crack’s depth,and the external load on the crack state and bending performances of the cracked beam is analyzed.It is revealed that a cusp exists on the deflection curve,and a jump on the rotation angle curve occurs at a crack location.The relation between the beam’s deflection and load is bilinear,each part corresponding to an open or closed state of crack,respectively.When the crack is open,flexibility of the cracked beam decreases with the increase of the beam’s slenderness ratio and the decrease of the crack depth.The results are useful in identifying non-destructive cracks on a beam.     

基于内聚力模型的白点萌生机理研究

基于内聚力理论、氢压理论和应力作用下氢富集理论,在以往孔洞锻合研究的基础上,建立了白点萌生扩展的有限元分析模型。研究了初始氢浓度、裂纹长度以及多裂纹耦合作用等多种因素对白点萌生的影响,确定了不同条件下白点萌生的临界氢浓度及其变化趋势。结果表明,氢对白点萌生具有重要的影响,在白点萌生过程中,高浓度的氢聚集在裂纹端部,极大地降低了钢的内聚力强度;对于单裂纹,随着裂纹长度的增加,白点萌生的临界氢浓度逐渐减小并趋于稳定值;当多个密集小裂纹同时存在时,裂纹间的耦合效应会导致小裂纹有贯穿形成大裂纹的趋势,而远端裂纹的耦合效应则相对较小;对于多裂纹,白点萌生临界氢浓度具有随裂纹个数的增加而线性递减的趋势。     

Nonlinear behaviour of interacting dielectric cracks in piezoelectric materials

Existing studies on the fracture of cracked piezoelectric materials have been limited mostly to the electrically impermeable and permeable crack models, which represent the limiting cases of the physical boundary condition along the crack surfaces. This paper presents a study on the electromechanical behaviour of interacting dielectric cracks in piezoelectric materials. The cracks are filled with dielectric media and, as the result, the electric boundary condition along the crack surfaces is governed by the opening displacement of the cracks. The formulation of this nonlinear problem is based on simulating the cracks using distributed dislocations and solving the resulting nonlinear singular integral equations. Multiple deformation modes are observed. A solution technique is developed to determine the desired deformation mode of the interacting cracks. Numerical results are given to show the effect of the interaction between parallel cracks. Attention is paid to the transition between permeable and impermeable models with increasing crack opening.     

An experimental and analytical treatment of matrix cracking in cross-ply laminates

The development of damage in cross-ply Hercules AS4/3502 graphite/epoxy laminates has been investigated. Specific endeavors were to identify the mechanisms for initiation and growth of matrix cracks and to determine the effect of matrix cracking on the stiffness loss in cross-ply laminates. Two types of matrix cracks were identified. These include both straight and curved cracks. The experimental study of matrix crack damage revealed that the curved cracks formed after the straight cracks and followed a repeatable pattern of location and orientation relative to the straight cracks. Therefore, it was postulated that the growth mechanism for curved cracks is driven by the stress state resulting from the formation of the straight cracks. This phenomenon was analytically investigated by a finite-element model of straight cracks in a cross-ply laminate. The finite-element results provide supporting evidence for the postulated growth mechanism. The experimental study also revealed that the number of curved cracks increased with the number of consecutive 90-deg plies. Finally, experimental results show as much as 10-percent degradation in axial stiffness due to matrix cracking in cross-ply graphite/epoxy laminates.     

单轴压缩下页岩类脆性材料细观损伤演化的离散元分析

了解广泛存在的类似页岩的脆性材料各向异性对工程安全具有重要意义。本研究将页岩视为粘结颗粒材料,基于离散单元方法研究了横观各向同性脆性页岩的损伤演化。再现了不同层理角的页岩试样的破坏模式,并对比了实验和数值模拟的抗压强度和弹性模量。引入微裂纹的概念,通过定义裂纹密度函数,系统地研究了单轴压缩条件下,页岩层理角对细观结构的影响。此外基于平均配位数建立了配位数变化与细观损伤的联系,并根据配位数的变化与裂纹数量将加载过程分为三个阶段,分析了不同阶段配位数与裂纹数量的对应变化关系。研究表明,页岩的裂纹密度随着层理角的增加而增加,而试样的平均配位数在加载过程中先上升后剧烈下降,颗粒集合体在单轴压缩条件下的应力应变及裂纹数量曲线与平均配位数曲线有良好的一致性。该研究揭示了横观各向同性脆性岩石的破坏过程和内在机理,将为页岩类脆性材料的工程应用提供理论指导。     

双周期裂纹压电材料反平面剪切问题断裂力学分析

研究压电材料双周期裂纹反平面剪切与平面电场作用的问题．运用复变函数方法，获得了该问题严格的闭合解，并由此给出了裂纹尖端应力强度因子和电位移强度因子的精确公式．数值算例显示了裂纹分布特征对材料断裂行为的重要影响．叠间小裂纹能够对主裂纹的应力和电位移场起着屏蔽作用，相反行间小裂纹却起着放大作用，至于钻石形分布裂纹的影响规律则更为复杂．对于某些特殊情形给予了解答并导出一系列有意义的结果。     

Environmentally assisted initiation and growth of multiple surface cracks

The initial stages of stress corrosion on an amorphous polymer is investigated. This is done by exposing stressed specimens of polycarbonate to an acetone and water solution. The surface develops two distinct features of degradation that appear on different length scales when subjected to tensile stress. Small pits form on the surface and make it rough. These pits are in the order of micrometers, and are found to be randomly distributed. They occur even without load and seem to slightly increase in number with increasing stress. In the millimeter domain, visible to the bare eye, surface cracks are formed transverse to the direction of loading. The occurrence of cracks is seen to have a positive stress threshold value, exceeding which, a linear increase of number of cracks with stress is found. The manners in which the cracks grow and coalesce on the surface are examined. It is seen that they do not meet crack tip to crack tip. Instead, they avoid each other initially and coalesce crack tip to crack side. The results are discussed in the light of mechanical considerations. A stress analysis for a few configurations of meeting cracks supports the experimental observations. With assumptions of stress corrosion crack growth and coalescence, a simulation of cracks growing from randomly distributed initiation sites is performed. Similar crack patterns as obtained in the experiments are found.     

Electrically driven cracks in piezoelectric ceramics: experiments and fracture mechanics analysis

Piezoelectric systems like multilayer actuators are susceptible to damage by crack propagation induced by strain incompatibilities. These can arise under electric fields for example between the electroded and external regions. Such incompatibilities have been realised in thin rectangular model specimens from PZT-piezoelectric ceramics with top and bottom electrodes only close to one edge. Under an electric field, controlled crack propagation has been observed in situ in an optical microscope. The crack paths are reproducible with very high accuracy. Small electrode widths lead to straight cracks with two transitions between stable and unstable crack growth regions, while large electrode widths result in curved cracks with four transitions. Fracture mechanics analysis is able to explain the different crack paths. An iteration method is developed to simulate the curved crack propagation also for strong curvature of the crack paths using the finite element method. The computed crack contours exhibit excellent quantitative agreement with the experiment with respect to their shape, the stages of stable and unstable crack propagation and the transitions between them. Finally, also the crack length as a function of the electric field can be predicted.     

Crack tip blunting and cleavage under dynamic conditions

In structural materials with both brittle and ductile phases, cracks often initiate within the brittle phase and propagate dynamically towards the ductile phase. The macroscale, quasistatic toughness of the material thus depends on the outcome of this microscale, dynamic process. Indeed, dynamics has been hypothesized to suppress dislocation emission, which may explain the occurrence of brittle transgranular fracture in mild steels at low temperatures (Lin et al., 1987). Here, crack tip blunting and cleavage under dynamic conditions are explored using continuum mechanics and molecular dynamics simulations. The focus is on two questions: (1) whether dynamics can affect the energy barriers for dislocation emission and cleavage, and (2) what happens in the dynamic “overloaded” situation, in which both processes are energetically possible. In either case, dynamics may shift the balance between brittle cleavage and ductile blunting, thereby affecting the intrinsic ductility of the material. To explore these effects in simulation, a novel interatomic potential is used for which the intrinsic ductility is tunable, and a novel simulation technique is employed, termed as a “dynamic cleavage test”, in which cracks can be run dynamically at a prescribed energy release rate into a material. Both theory and simulation reveal, however, that the intrinsic ductility of a material is unaffected by dynamics. The energy barrier to dislocation emission appears to be identical in quasi-static and dynamic conditions, and, in the overloaded situation, ductile crack tip behavior ultimately prevails since a single emission event can blunt and arrest the crack, preventing further cleavage. Thus, dynamics cannot embrittle a ductile material, and the origin of brittle failure in certain alloys (e.g., mild steels) appears unrelated to dynamic effects at the crack tip.     

Simulating a crack by a notch in the solution of dynamic problems

In the article the method of dynamic photoelasticity is used to investigate the formation of a state of stress at the apexes of cracks and notches under the action of stress waves. Edge and closed cracks in the case of the action of longitudinal and surface waves are discussed. The investigations show that a crack may be simulated by a notch only under definite conditions.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 173–178, July–August, 1973.The author thanks Yu. N. Rabotnov for his advice and for his observations during evaluation of the work.     

不同节理角度和地应力条件下岩石双孔爆破的数值模拟

在采矿工程、地下交通工程和水利水电工程等岩体开挖工程中,爆破法仍是一种主要的破岩方法. 实际天然岩体中存在的裂隙、节理等不连续面和所处的地应力场环境,都会对爆破荷载传播过程和最终的破碎效果产生重要影响. 本文把岩石爆破视为爆炸应力波动态作用和爆生气体压力准静态作用两个独立的先后作用过程,同时考虑初始应力场的静态作用,建立岩石爆破的力学模型. 基于渗流方程描述爆生气体在爆生裂纹中的传播过程,进而基于流固耦合理论实现爆生气体准静态压力对裂纹尖端的力学作用. 爆炸应力波主要在孔周引起压碎区和径向微裂纹区,随后爆生气体压力楔入径向裂纹尖端,促使裂纹进一步扩展,形成径向主裂纹. 数值模型可以再现孔周压碎区、径向微裂纹区、径向主裂纹区萌生、扩展的完整演化过程. 针对不同节理角度和地应力条件下岩石双孔爆破过程的数值模拟结果表明,初始地应力场的压应力作用不利于爆生裂纹的萌生与扩展,但节理的存在对裂纹的扩展具有明显的导向和促进作用,有利于爆生裂纹沿节理面方向的扩展.      

Numerical computation of central crack growth in an active particle of electrodes influenced by multiple factors

Mechanical degradation, especially fractures in active particles in an electrode, is a major reason why the capacity of lithium-ion batteries fades. This paper proposes a model that couples Li-ion diffusion, stress evolution, and damage mechanics to simulate the growth of central cracks in cathode particles \((\hbox {LiMn}_{2}\hbox {O}_{4})\) by an extended finite element method by considering the influence of multiple factors. The simulation shows that particles are likely to crack at a high discharge rate, when the particle radius is large, or when the initial central crack is longer. It also shows that the maximum principal tensile stress decreases and cracking becomes more difficult when the influence of crack surface diffusion is considered. The fracturing process occurs according to the following stages: no crack growth, stable crack growth, and unstable crack growth. Changing the charge/discharge strategy before unstable crack growth sets in is beneficial to prevent further capacity fading during electrochemical cycling.     

Special rotational deformation as a toughening mechanism in nanocrystalline solids

A theoretical model is suggested which describes the effect of special rotational deformation on crack growth in deformed nanocrystalline ceramics and metals. Within the model, the special rotational deformation (driven by the external stress concentrated near the tip of a mode I crack) occurs in a nanograin through formation of immobile disclinations whose strengths gradually increase during the formation process conducted by grain boundary sliding and diffusion. The special rotational deformation releases, in part, local stresses near the crack tip, thus serving as a toughening mechanism in nanocrystalline materials. The effects of the special rotational deformation on the growth of pre-existent, comparatively large cracks in nanocrystalline metals and ceramics are estimated.     

基于等效弹簧模型的裂纹Euler-Bernoulli梁弯曲变形分析

考虑裂纹的缝隙效应,研究了开闭裂纹Euler-Bernoulli梁的弯曲变形．首先,将裂纹等效为内部旋转弹簧,利用广义函数,给出了考虑裂纹缝隙影响的Euler-Bernoulli梁的等效抗弯刚度,推导了具有任意数目开闭裂纹梁弯曲变形的显式通解．在此基础上,研究了均布载荷作用下上侧单裂纹简支梁以及裂纹处承受集中力和集中力偶共同作用的固支梁的弯曲变形,分析了梁长细比、裂纹深度和位置以及载荷等对裂纹开闭状态和梁弯曲变形的影响。结果表明：梁挠度分布在裂纹处存在尖点,而转角分布存在跳跃;梁挠度与载荷的响应关系一般为双折线形式,分别对应于裂纹的张开和闭合状态;且裂纹张开时,裂纹梁的柔度随着梁长细比的增加和裂纹深度的减小而减小。这些结果对梁裂纹无损检测具有指导意义．     

Combined analysis of fracture under stresses acting along cracks

A joint approach to the study of two non-classical fracture mechanisms, namely fracture of cracked materials with initial (residual) stresses acting along the crack planes and fracture under compression along parallel cracks, is considered in the framework of three-dimensional linearized solid mechanics. Mathematical statements of problems for pre-stressed solids that contain interacting circular cracks are given. Problems for an infinite solid containing two parallel co-axial cracks and for a space with the periodical set of co-axial parallel cracks as well as for a half-space with near-the-surface crack are solved. Several patterns of loading on the crack faces (normal loading, radial shear and torsion) are considered. The effects of initial stresses on stress intensity factors are analyzed for highly elastic materials with some types of elastic potentials. Formulation of fracture criteria accounting effect of initial (residual) stresses is given. Critical parameters of fracture of solids containing interacting cracks under compression along the cracks are calculated. The influence of geometrical parameters of the problems as well as physical and mechanical properties of materials on these critical parameters is analyzed.     

On the dynamics of crack growth in glass

An experimental investigation of crack growth has been performed using plates of glass subjected to impulsive loading by means of small explosive charges. For comparison, some experiments were also performed with two “glass-like” polymers. The formation of the resultant crack patterns in glass was observed by the use of high-speed photography during the period following the detonation. Information has been obtained about the physical features of crack growth with regard to crack velocity, bifurcation, and the interaction between growing cracks and propagating stress waves. The existence of a terminal velocity of crack propagation is considered and, in particular, the degree to which crack bifurcation provides an explanation for this phenomenon is discussed.