基于圆弧底试件的动态裂纹扩展及止裂规律研究

为了研究脆性材料的动态裂纹扩展及止裂规律，设计了一种带圆弧形底边的梯形开口边裂纹（trapezoidal opening crack with arc bottom，TOCAB）构型的试件。在落锤冲击设备加载下，对圆心角为0°、60°、90°和120°的TOCAB试件进行了冲击实验，并采用裂纹扩展计(crack propagation gauge，CPG)监测裂纹起裂和扩展时间，从而获得裂纹扩展速度。采用有限差分软件AUTODYN对落锤冲击设备和试件进行数值模拟，研究了裂纹的动态扩展过程及止裂规律。还基于实验和数值方法，计算了裂纹的临界动态应力强度因子。实验和数值结果均表明:3种弧度的TOCAB试件都可以实现运动裂纹止裂，该构型可用于研究动态裂纹止裂问题；数值计算的裂纹扩展路径与实验结果基本一致，验证了数值模型的有效性；裂纹起裂和止裂时刻的临界动态应力强度因子大于裂纹动态扩展过程中的临界动态应力强度因子。     

含预制裂纹L形梁柱试件动态断裂过程

针对含预制裂纹L形梁柱试件,为研究预制裂纹动态扩展的力学特征及其对梁柱试件破坏模式的影响,采用数字动态焦散线实验系统,对距节点核心区不同距离l处含有预制裂纹的试件进行落锤冲击实验,得到预制裂纹的扩展轨迹、速度、动态应力强度因子的变化规律。结果表明,l值增大,扩展裂纹在梁下边缘的贯通点与预制裂纹的夹角逐渐增大,曲裂程度变大。裂纹扩展速度随着l的增大振荡性增强,裂纹扩展平均速度逐渐降低。l值为2 mm时,裂尖表现为Ⅰ型断裂,l值增大,裂尖受到剪应力作用增强,Ⅰ型动态应力强度因子减小,Ⅱ型动态应力强度因子增大,断裂逐渐转变为Ⅰ-Ⅱ复合型。     

Ⅰ型裂纹中低速冲击荷载下起裂韧度测试新方法

为了探寻更加合理的构型试件来研究纯Ⅰ型裂纹在冲击荷载下的起裂及扩展行为, 提出一种新构型试件, 即双倾斜底边中心裂纹试件(double inclined bottom central cracked, DIBCC)。借助于中低速落锤式冲击实验装置进行冲击实验, 通过应力波来使试件内预制裂纹起裂并扩展, 同时利用应变片测试系统监测裂纹起裂时刻, 并采用AUTODYN有限差分软件对实验过程进行数值模拟, 最后计算裂纹的动态应力强度因子, 利用实验测得的起裂时刻, 确定试件的起裂韧度。结果表明:(1)在反射拉伸波作用下, 预制裂纹两侧会产生垂直于裂纹面向外的位移, 使预制裂纹扩张, 从而使裂纹起裂。(2)数值模拟结果与实验结果在裂纹扩展路径上具有一致性, 说明本文中提出的DIBCC构型试件有效, 可以用来测试裂纹在冲击载荷下的断裂韧度。     

旋转防喷器壳体疲劳裂纹扩展研究

受井内高压和岩屑颗粒反复碰撞，旋转防喷器壳体会出现裂纹，裂纹的进一步扩展会导致壳体发生破裂。本研究基于Paris公式研究旋转防喷器壳体的疲劳裂纹扩展寿命，通过疲劳裂纹扩展试验确定Paris公式的材料常数，通过有限元模拟对CT试样进行疲劳裂纹扩展研究，分析初始裂纹尺寸、载荷大小的影响。结果显示：实验结果与有限元仿真结果能较好吻合；随着初始裂纹尺寸的增加，裂纹前缘应力强度因子逐渐增大，初始裂纹开始扩展，深度亦会随应力幅和应力比的增大而增大，而裂纹扩展寿命随着裂纹长度和载荷增大而减小，随着应力比的增大而增大。该研究为旋转防喷器安全使用提供一定的理论依据和研究方法。     

冲击作用下含预制裂纹梁柱试件的动态断裂

利用焦散线实验系统,进行了冲击加载下含预制裂纹梁柱试件的断裂实验,研究了梁柱试件的梁柱节点、梁上和柱端裂纹的扩展轨迹、扩展速度和应力强度因子的变化规律。实验结果表明:受冲击后,试件首先在梁柱节点处开裂,并在裂纹扩展过程中发生明显的曲裂运动,证明梁柱节点处最容易受到破坏。预制裂纹条数越多,梁柱节点处开裂越晚,说明裂纹条数少,能量可以在裂纹尖端积聚得更集中、更快。含柱端预制裂纹的试件,2条裂纹的开裂相隔时间要长于含梁上预制裂纹的试件;同时,试件的第2条裂纹优先在梁上裂纹处开裂,说明固端支座比简支梁断裂需要更多的能量。随着预制裂纹的增多,梁上裂纹在扩展过程中的曲裂现象减弱, 由于部分能量在柱端裂纹处积聚,用于推动梁上裂纹扩展的能量相应地减少。并且由于柱端裂纹的存在,梁上裂纹受到的拉应力分量减小,导致裂纹尖端受到弯矩变小,影响了裂纹的曲裂运动。     

致密砂岩巷道模型试件动态起裂及止裂全过程分析

为了开展在不同冲击载荷作用下巷道围岩内裂纹的起裂、扩展及止裂等问题，以可调速冲击试验机进行动态加载试验，采用致密青砂岩制作裂纹巷道模型试件，并利用裂纹扩展计分别记录了动态起裂、扩展、止裂等时刻，对裂纹扩展速度的变化规律进行分析；随后采用AUTODYN有限差分法软件进行相应的数值模拟，数值模拟得到的裂纹扩展路径与试验结果基本一致。经过两者对比分析可知:随着冲击载荷作用的增加，裂纹平均扩展速度逐渐增大，随后趋于稳定值；预制裂纹的起裂时间随着冲击速度载荷的增加而逐渐降低，并在稳定值上下波动；随着冲击速度载荷的增加，裂纹扩展路径过程中的止裂时段逐渐变短。     

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

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

正交弹性材料中双周期裂纹反平面问题的封闭解

研究了无限大正交弹性材料中含双周期裂纹的反平面问题,其基本胞元含有三条裂纹,且三条裂纹的中心恰好位于一等腰三角形顶点。运用椭圆函数、保角变换理论、施瓦兹公式获得了该问题应力场的封闭解,并得到了裂纹尖端处的应力强度因子。该问题结果取特殊情形退化对应于单排共线周期裂纹的解答。通过数值算例,分析了双周期裂纹归一化的应力强度因子随双周期裂纹的横向间距和纵向间距之比/b a 分别取10、5、2、1时的变化曲线。结果表明：对于一定的横向间距,应力奇异因子随纵向间距的增大而减小,但随着纵向间距的增大,纵向间距对应力奇异因子的影响变得不明显;对于一定的纵向间距,应力奇异因子随横向间距的减小而减小,但随着横向间距的减小,横向间距对应力奇异因子的影响变得不明显。     

爆炸荷载对不同距离邻近隧洞背爆侧裂纹影响的模拟试验

利用爆炸加载动态焦散线测试系统,模拟了不同起爆距离L时,爆炸荷载作用于邻近隧洞的全过程;通过分析迎爆侧破坏现象和背爆侧预制裂纹扩展的位移、速度以及裂纹尖端动态应力强度因子曲线的变化规律,研究了不同间距并行隧洞爆破开挖对既有隧洞的影响。试验结果表明:随着L增大,隧洞背爆侧裂纹的最终扩展位移呈先增大后减小的变化规律;裂纹扩展速度曲线的振荡幅值以及频率呈降低趋势,Ⅰ型应力强度因子的变化曲线逐渐趋于平缓,且最大峰值均出现延迟现象,最大峰值所对应的时间点由t=50μs~100μs向t100μs转移;起爆距离L从50mm增大到60mm时,裂纹扩展的位移、速度曲线以及应力强度因子曲线出现显著变化;而小于50mm或大于60mm时,它们分别呈相似性变化。研究结果可为现场试验确定邻近隧洞爆破开挖的距离L提供参考。     

基于DIC的爆炸加载下脆性材料裂纹扩展规律的试验研究

通过数字图像相关法（DIC），应用PMMA对爆炸加载条件下脆性材料的裂纹扩展规律进行了试验研究。基于对称性试验模型，实现了裂纹尖端位置和应变场信息的同步记录。以此为基础，通过对比分析获知，主应变场应变值最大点不能作为裂纹尖端的判断依据。并以动态裂纹扩展速度为参量，应用断裂动力学和最小二乘牛顿迭代法，计算出了考虑惯性效应的Ⅰ-Ⅱ混合型裂纹的应力强度因子:K_Ⅰ和K_Ⅱ值会随着裂纹扩展方向改变而发生突变；K_Ⅰ最大值为2.63 MPa·m1/2，最小值为0.89 MPa·m1/2；其整体变化趋势表明，爆炸加载条件下脆性材料裂纹扩展随能量积聚和释放呈循环阶梯式递减发展。     

Effect of plastic zone on the fatigue crack propagation behavior between two fatigue cracks

In this study, the fatigue crack propagation behavior in the stress interaction field between two different fatigue cracks is studied by experiment and finite element analysis. In the experiment, the offset distance between two cracks and the applied stress are varied to create different stress interaction fields. The size of the plastic zone area is used to examine the crack propagation path and rate. Three types of crack propagation in the interaction field were found by experiment, and the crack propagation behavior of two cracks was significantly changed as different stresses were applied. The size of the plastic zone obtained by finite element analysis can be used to explain crack propagation behavior qualitatively.     

Experimental study on cracking behaviour of moulded gypsum containing two non-parallel overlapping flaws under uniaxial compression

Failure of rock mass that is subjected to compres-sive loads occurs from initiation, propagation, and linkage of new cracks from preexisting fissures. Our research inves-tigates the cracking behaviour and coalescence process in a brittle material with two non-parallel overlapping flaws using a high-speed camera. The coalescence tensile crack and tensile wing cracks were the first cracks to occur from the pre-existing flaws. The initiation stresses of the primary cracks at the two tips of each flaw were simultaneous and decreased with reduced flaw inclination angle. The following types of coalescence cracks were identified between the flaws: pri-mary tensile coalescence crack, tensile crack linkage, shear crack linkage, mixed tensile-shear crack, and indirect crack coalescence. Coalescence through tensile linkage occurred mostly at pre-peak stress. In contrast, coalescence through shear or mixed tensile-shear cracks occurred at higher stress. Overall, this study indicates that the geometry of preexisting flaws affect crack initiation and coalescence behaviour.     

Microscopic analysis of crack propagation for multiple cracks, inclusions and voids

The elastic crack interaction with internal defects, such as microcracks, voids and rigid inclusions, is investigated in this study for the purpose of analyzing crack propagation. The elastic stress field is obtained using linear theory of elasticity for isotropic materials. The cracks are modeled as pile-ups of edge dislocations resulting into a coupled set of integral equations, whose kernels are those of a dislocation in a medium with or without an inclusion or void. The numerical solution of these equations gives the stress intensity factors and the complete stress field in the given domain. The solution is valid for a general solid, however the propagation analysis is valid mostly for brittle materials. Among different propagation models the ones based on maximum circumferential stress and minimum strain energy density theories, are employed. A special emphasis is given to the estimation of the crack propagation direction that defines the direction of crack branching or kinking. Once a propagation direction is determined, an improved model dealing with kinked cracks must be employed to follow the propagation behavior.     

残余应力对压力容器破坏模式的影响分析

破前漏(简称LBB)是压力容器、核电站设备结构设计与评价中的一个重要准则.表面裂纹准静态扩展的几何形貌变化规律的预测是破前漏(LBB)评判十分重要的课题之一.本文对特定焊接残余应力场加载作用下,含三维表面裂纹的压力容器模型,用有限元软件(ABAQUS)进行了表面裂纹准静态扩展模拟计算,得到在此残余应力场作用下应力强度因子沿裂纹前缘的分布规律.结合外载引起的应力强度因子,就可以判别裂纹的扩展形貌,从而判断结构是否满足LBB要求.     

Brittle fracture dynamics with arbitrary paths I. Kinking of a dynamic crack in general antiplane loading

We present a new method for determining the elasto-dynamic stress fields associated with the propagation of anti-plane kinked or branched cracks. Our approach allows the exact calculation of the corresponding dynamic stress intensity factors. The latter are very important quantities in dynamic brittle fracture mechanics, since they determine the crack path and eventual branching instabilities. As a first illustration, we consider a semi-infinite anti-plane straight crack, initially propagating at a given time-dependent velocity, that changes instantaneously both its direction and its speed of propagation. We will give the explicit dependence of the stress intensity factor just after kinking as a function of the stress intensity factor just before kinking, the kinking angle and the instantaneous velocity of the crack tip.     

Basic solution of two parallel non-symmetric permeable cracks in piezoelectric materials

The behavior of two parallel non-symmetric cracks in piezoelectric materials subjected to the anti-plane shear loading was studied by the Schmidt method for the permeable crack electric boundary conditions. Through the Fourier transform, the present problem can be solved with two pairs of dual integral equations ip which the unknown variables are the jumps of displacements across crack surfaces. To solve the dual integral equations, the jumps of displacements across crack surfaces were directly expanded in a series of Jacobi polynomials. Finally, the relations between electric displacement intensity factors and stress intensity factors at crack tips can be obtained. Numerical examples are provided to show the effect of the distance between two cracks upon stress and electric displacement intensity factors at crack tips. Contrary to the impermeable crack surface condition solution, it is found that electric displacement intensity factors for the permeable crack surface conditions are much smaller than those for the impermeable crack surface conditions. At the same time, it can be found that the crack shielding effect is also present in the piezoelectric materials.     

The simulation of dynamic crack propagation using the cohesive segments method

The cohesive segments method is a finite element framework that allows for the simulation of the nucleation, growth and coalescence of multiple cracks in solids. In this framework, cracks are introduced as jumps in the displacement field by employing the partition of unity property of finite element shape functions. The magnitude of these jumps are governed by cohesive constitutive relations. In this paper, the cohesive segments method is extended for the simulation of fast crack propagation in brittle solids. The performance of the method is demonstrated in several examples involving crack growth in linear elastic solids under plane stress conditions: tensile loading of a block; shear loading of a block and crack growth along and near a bi-material interface.