Discrete Element Modeling of Asphalt Concrete Cracking Using a User-def ined Three-dimensional Micromechanical Approach

We established a user-defined micromechanical model using discrete element method （DEM） to investigate the cracking behavior of asphalt concrete （AC）. Using the ＂Fish＂ language provided in the particle flow code in 3-Demensions （PFC3D）, the air voids and mastics in asphalt concrete were realistically built as two distinct phases. With the irregular shape of individual aggregate particles modeled using a clump of spheres of different sizes, the three-dimensional （3D） discrete element model was able to account for aggregate gradation and fraction. Laboratory uniaxial complex modulus test and indirect tensile strength test were performed to obtain input material parameters for the numerical simulation. A set of the indirect tensile test were simulated to study the cracking behavior of AC at two levels of temperature, i e, -10 ℃ and 15 ℃. The predicted results of the numerical simulation were compared with laboratory experimental measurements. Results show that the 3D DEM model is able to predict accurately the fracture pattern of different asphalt mixtures. Based on the DEM model, the effects of air void content and aggregate volumetric fraction on the cracking behavior of asphalt concrete were evaluated.

The interaction between a screw dislocation and a wedge-shaped crack in one-dimensional hexagonal piezoelectric quasicrystals

Based on the fundamental equations of piezoelasticity of quasicrystal material,we investigated the interaction between a screw dislocation and a wedge-shaped crack in the piezoelectricity of one-dimensional hexagonal quasicrystals.Explicit analytical solutions are obtained for stress and electric displacement intensity factors of the crack,as well as the force on dislocation.The derivation is based on the conformal mapping method and the perturbation technique.The influences of the wedge angle and dislocation location on the image force are also discussed.The results obtained in this paper can be fully reduced to some special cases already available or deriving new ones.

THE CONSTITUTIVE RELATION OF CRACK-WEAKENED ROCK MASSES UNDER AXIAL-DIMENSIONAL UNLOADING

An accurate and efficient numerical method for solving the crack-crack interaction problem is presented. The method is mainly by means of the dislocation model, stress superposition principle and Chebyshev polynomial expansion of the pseudo-traction. This method can be applied to compute the stress intensity factors of multiple kinked cracks and multiple rows of periodic cracks as well as the overall strains of rock masses containing multiple kinked cracks under complex loads. Many complex computational examples are given. The dependence of the crack-crack interaction on the crack configuration, the geometrical and physical parameters, and loads pattern, is investigated. By comparison with numerical results under confining pressure unloading, it is shown that the crack-crack interaction under axial-dimensional unloading is weaker than those under confining pressure unloading. Numerical results for single faults and crossed faults show that the single faults are more unstable than the crossed faults. It is found from numerical results for different crack lengths and different crack spacing that the interaction among kinked cracks decreases with an increase in length of the kinked cracks and the crack spacing under axial-dimensional unloading.

THREE-DIMENSIONAL ELLIPTIC CRACK UNDER IMPACT LOADING

The dynamic stress intensity factor of a three-dimensionalelliptic crack under impact loading is determined with the finiteelement method. The computation results can take into account theinfluence of time and the ratio of the wave speeds on the stressintensity factor. The present method is suitable not only forthree-dimensional dynamic crack, but also for three-dimensionaldynamic contact.

Exact analytic solutions for an elliptic hole with asymmetric collinear cracks in a one-dimensional hexagonal quasi-crystal

Using the complex variable function method and the technique of conformal mapping, the anti-plane shear problem of an elliptic hole with asymmetric colfinear cracks in a one-dimensional hexagonal quasi-crystal is solved, and the exact analytic solutions of the stress intensity factors （SIFs） for mode Ⅲ problem are obtained. Under the limiting conditions, the present results reduce to the Griffith crack and many new results obtained as well, such as the circular hole with asymmetric collinear cracks, the elliptic hole with a straight crack, the mode T crack, the cross crack and so on. As far as the phonon field is concerned, these results, which play an important role in many practical and theoretical applications, are shown to be in good agreement with the classical results.

Interactions of penny-shaped cracks in three-dimensional solids

The interaction of arbitrarily distributed penny-shaped cracks in three-dimensional solids is analyzed in this paper. Using oblate spheroidal coordinates and displacement functions, an analytic method is devel- oped in which the opening and the sliding displacements on each crack surface are taken as the basic unknown functions. The basic unknown functions can be expanded in series of Legendre polynomials with unknown coefficients. Based on superposition technique, a set of governing equations for the unknown coefficients are formulated from the traction free conditions on each crack surface. The boundary collocation procedure and the average method for crack-surface tractions are used for solving the governing equations. The solution can be obtained for quite closely located cracks. Numerical examples are given for several crack problems. By comparing the present results with other existing results, one can conclude that the present method provides a direct and efficient approach to deal with three-dimensional solids containing multiple cracks.

Three dimensional K-Tz stress fields around the embedded center elliptical crack front in elastic plates

Through detailed three-dimensional （3D） finite element （FE） calculations, the out-of-plane constraints Tz along embedded center-elliptical cracks in mode I elastic plates are studied. The distributions of Tz are obtained near the crack front with aspect ratios （a/c） of 0.2, 0.4, 0.5, 0.6, 0.8 and 1.0. Tz decreases from an approximate value of Poisson ratio v at the crack tip to zero with increasing normalized radial distances （r/a） in the normal plane of the crack front line, and increases gradually when the elliptical parameter angle φ changes from 0° to 90°at the same r/a. With a/c rising to 1.0, Tz is getting nearly independent of φ and is only related to r/a. Based on the present FE calculations for Tz, empirical formulas for Tz are obtained to describe the 3D distribution of Tz for embedded center-elliptical cracks using the least squares method in the range of 0.2 ≤ a/c ≤ 1.0. These Tz results together with the corresponding stress intensity factor K are well suitable for the analysis of the 3D embedded centerelliptical crack from field, and a two-parameter K-Tz principle is proposed.

A CLOSED FORM SOLUTION OF STRESS INTENSITY FACTORS FOR THREE DIMENSIONAL FINITE BODIES WITH ECCENTRIC CRACKS

In this paper, a new analytical-engineering method of closed form solution about stress intensity factors for three dimensional finite bodies with eccentric cracks is derived by means of energy release rate method. The results of stress intensity factors can be obtained. The results provided ir this method are in nice agreement with those of the famous alternating method by which only special cases can be solved.

A modified maximum tangential tensile stress criterion for three-dimensional crack propagation

The three-dimensional （3D） crack propagation is a hot issue in rock mechanics. To properly simulate 3D crack propagation, a modified maximum tangential tensile stress criterion is proposed. In this modified criterion, it is supposed that cracks propagate only at crack front in the principal normal plane. The tangential tensile stress at crack front in the principal normal plane in local coordinates is employed to determine crack propagation, which is calculated through coordinate transformation from global to local coordinates. New cracks will propagate when the maximum tangential tensile stress at crack front in the principal normal plane reaches the tensile strength of rock-like materials. Compared with the previous crack propagation criteria, the modified crack propagation criterion is helpful in calculating 3D crack stress intensity factor, and can overcome the limitations of propagation step determined by individual experiences in previous studies. Finally, the 3D crack propagation process is traced by element-free Galerkin method. The numerical results agree well with the experimental ones for a frozen resin sample with prefabricated 3D cracks.

A STUDY ON STRESS INTENSITY FACTORS AND SINGULAR STRESS FIELDS OF THREE-DIMENSIONAL INTERFACE CRACK

By using the finite-part integral concepts and limit technique,the hypersingular inte- grodifferential equations ofthree-dimensional(3D)planar interface crack were obtained; then thedominant-part analysis of 2D hypersingular integral was further usedto investigate the stress fields near the crack front theoretically,and the accurate formulae were obtained for the singular stressfields and the complex stress intensity factors.

基于Deepcrack网络的混凝土裂缝检测方法

混凝土结构裂缝对建筑安全构成了极大的潜在威胁,裂缝检测对建筑结构的维护具有重要意义,当前基于深度学习的裂缝检测针对提取裂缝细节的能力仍有待提高。因此,该文对Deepcrack网络进行优化,提出了基于金字塔分割注意力和全局上下文的混凝土裂缝检测算法PG-Deepcrack。首先,在编码器中提出双卷积-注意力并行模块,增加金字塔分割注意力分支为卷积层提供更丰富的多尺度裂缝信息;其次,为了捕获长距离依赖关系,并行模块操作后引入全局上下文模块,进一步提升网络对裂缝细节的表达能力;最后,在特征融合阶段利用全维动态卷积和GELU激活函数,对编解码器特征层联级融合,使网络更全面地保留不同尺寸的裂缝信息并提高模型的泛化性能。为验证网络模型的有效性,在Deepcrack数据集上与7个网络模型进行对比试验,所提出的网络表现了最佳性能,IoU达到了72.78%。

CONVOLUTION OF THE IMPACT THREE-DIMENSIONAL ELASTO-DYNAMICS AND DYNAMIC STRESS INTENSITY FACTOR FOR AN ELLIPTIC CRACK

This paper presents a formulation for three-dimensional elasto-dynamics with an elliptic crack based on the Laplace and Fourier transforms and the convolution theorem. The dynamic stress intensity factor for the crack is determined by solving a Fredholm integral equation of the first kind. The results of this paper are very close to those given by the two-dimensional dual integral equation method.

Prediction of Three Dimensional Crack Path Brittle Fracture in Weldment under Dynamic Load

Three dimensional dynamic stress intensity factors are analyzed for a curved crack with a second order perturbation method. The method is extended to obtain an approximate representation of a three dimensional dynamic stress intensity factors at the tip of a curved crack. Due to three dimensional curved crack growth the dynamic energy release rate can be calculated by using the Irwin＇s formula. A three dimensional curved crack in materials with inhomogeneous fracture toughness are considered. Paths of a brittle three dimensional curved crack propagating along a welded joint are predicted via the present method, where the effects of dynamic applied stresses, residual stresses, and material deterioration due to welding are taken into considerations.

A MODE Ⅱ CRACK IN A TWO-DIMENSIONAL OCTAGONAL QUASICRYSTALS

The present study develops the fracture theory for a two-dimensional octagonal quasicrystals. The exact analytic solution of a Mode Ⅱ Griffith crack in the material was obtained by using the Fourier transform and dual integral equations theory, then the displacement and stress fields, stress intensity factor and strain energy release rate were determined, the physical sense of the results relative to phason and the difference between mechanical behaviors of the crack problem in crystal and quasicrystal were figured out. These provide important information for studying the deformation and fracture of the new solid phase.

INTERFACIAL CRACK ANALYSIS IN THREE-DIMENSIONAL TRANSVERSELY ISOTROPIC BI-MATERIALS BY BOUNDARY INTEGRAL EQUATION METHOD

The integral-differential equations for three-dimensional planar interfacial cracks of arbitrary shape in transversely isotropic bimaterials were derived by virtue of the Somigliana identity and the fundamental solutions, in which the displacement discontinuities across the crack faces are the unknowns to be determined. The interface is parallel to both the planes of isotropy. The singular behaviors of displacement and stress near the crack border were analyzed and the stress singularity indexes were obtained by integral equation method. The stress intensity factors were expressed in terms of the displacement discontinuities. In the non-oscillatory case, the hyper-singular boundary integral-differential equations were reduced to hyper-singular boundary integral equations similar to those of homogeneously isotropic materials.

Calculating changes in fractal dimension of surface cracks to quantify how the dynamic loading rate affects rock failure in deep mining

The split-Hopkinson pressure bar(SHPB)and digital image correlation(DIC)techniques are combined to analyze the dynamic compressive failure process of coal samples,and the box fractal dimension is used to quantitatively analyze the dynamic changes in the coal sample cracks under impact load conditions with different loading rates.The experimental results show that the fractal dimension can quantitatively describe the evolution process of coal fractures under dynamic load.During the dynamic compression process,the evolution of the coal sample cracks presents distinct stages.In the crack propagation stage,the fractal dimension increases rapidly with the progress of loading,and in the crack widening stage,the fractal dimension increases slowly with the progress of loading.The initiation of the crack propagation phase of the coal samples gradually occurs more quickly with increasing loading rate;the initial cracks appear earlier.At the same loading time point,when the loading rate is greater,the fractal dimension of the cracks observed in the coal sample is greater.

Relevance between abutment pressure and fractal dimension of crack network induced by mining

Based on the geological conditions of coal mining face No.15-14120 at No.8 mine of Pingdingshan coal mining group,the real-time evolution of coal-roof crack network with working face advancing was collected with the help of intrinsically safe borehole video instrument.And according to the geology of this working face,a discrete element model was calculated by UDEC.Combining in situ experimental data with numerical results,the relationship between the fractal dimension of boreholes'wall and the distribution of advanced abutment pressure was studied under the condition of mining advance.The results show that the variation tendency of fractal dimension and the abutment pressure has the same characteristic value.The distance between working face and the peak value of the abutment pressure has a slight increasing trend with the advancing of mining-face.When the working face is set as the original point,the trend of fractal dimension from the far place to the origin can be divided into three phases:constant,steady increasing and constant.And the turning points of these phases are the max-influencing distance(50 m)and peak value(15 m)of abutment pressure.

Characterizing Dessication Cracking of a Remolded Clay (Eutric Vertisol) Using the Fractal Dimension Approach

Fractal dimension fd, was used as one of the parameters to describe dessicationcracking pattern of a remolded Black Cotton soil (Eutric Vertisol). The fractal dimension computed from filtered, thinned and skeletonized binary images of soil cracks using the Fractal3 software provided an insight into temporal variability of fd as well as its relationship with the Crack Intensity Factor (CIF) and Soil Moisture Content (SMC). The results showed that even for single crack, the fd prior to filtering and thinning were higher than after. Cracking patterns were observedfroma chosen soil sample during dessication and the corresponding relationship between fd and CIF compared and monitored. As the critical SMC decreased during drying (45% to 27%), the CIF soil increased (0.023% - 5.75%), so did the fd (1.233 to 1.7193). The fd showed a positive linear correlation with CIF at r2 = 0.247 (P fd with SMC was best described using a polynomial function at r2 = 0.969 (P fd was sensitive to dessication cracking and therefore on SMC changes. Visual observation of dessication cracking showed that CIF increased and attained stability after day 4 while the computed and logarithmic transformed crack area attained stability between days 7 to 10 gradually decreasing to values below 2%. The estimated crack Cover or Brightness of the digitized binary images also gave better approximation of the CIF though this was slightly higher. Our results showed that dessication cracking of the Eutric Vertisol was independent of antecedent critical SMC and was time-constrained. Further soil cracking therefore stopped once maximum CIF was attained and only widening and deepening of pre-existing cracks continued.

炸药模拟裂纹宽度的CT检测与测量分析

炸药作为武器系统的核心部件,其内部结构质量直接关系到武器系统的可靠性,目前CT技术已被用于炸药试件的内部质量检测与评价研究,特别用于炸药试件内部裂纹的定性与定量分析中。然而选择不同的检测工艺参数会造成CT图像质量的变化,直接影响到后续定量表征的精度。通过在炸药试件表面加工不同宽度的标准尺寸线槽,再将其与相同规格的炸药试件对粘,以实现在炸药试件内部预制标准尺寸裂纹的目的。然后,改变检测参数与工艺条件,得到不同检测参数对试件裂纹宽度测量精度的影响规律,分析讨论了CT测量炸药试件(如尺寸小于10 mm的始发药等)裂纹宽度的最佳检测工艺条件。

喷射混凝土围岩复合体受压力学特性及本构模型

为探究喷射混凝土围岩复合体受压力学特性,结合三维数字图像相关技术,针对4种不同喷射混凝土围岩复合体试件开展了单轴受压试验,探究了纤维和岩石种类对复合体破坏模式和受压性能的影响.基于Weibull统计损伤理论,构建了适用于复合体的受压应力应变本构模型.试验结果表明,不同围岩基复合体破坏模式存在显著差异,但微裂缝均先从喷射混凝土中发展.破坏时,围岩部分表现单一裂纹模式,而喷射混凝土部分呈现多裂纹扩展状态.4种复合体峰值强度均高于单一喷射混凝土破坏强度,且砂岩基复合体的峰值强度均高于花岗岩基复合体,平均高出5.55%.同种围岩条件下,纤维增强喷射混凝土复合体强度会存在轻微降低,较普通喷射混凝土复合体平均下降1.95%.此外,所建本构模型计算结果与试验数据吻合较好.