Analysis of explosion wave interactions and rock breaking effects during dual initiation

In blasting engineering, the location and number of detonation points, to a certain degree, regulate the propagation direction ofthe explosion stress wave and blasting effect. Herein, we examine the explosion wave field and rock breaking effect in terms of shockwave collision, stress change of the blast hole wall in the collision zone, and crack propagation in the collision zone. The produced shockwave on the collision surface has an intensity surpassing the sum of the intensities of the two colliding explosion shock waves. At the collisionlocation, the kinetic energy is transformed into potential energy with a reduction in particle velocity at the wave front and the wavefront pressure increases. The expansion form of the superposed shock wave is dumbbell-shaped, the shock wave velocity in the collisionarea is greater than the radial shock wave velocity, and the average propagation angle of the explosion shock waves is approximately 60°.Accordingly, a fitted relationship between blast hole wall stress and explosion wave propagation angle in the superposition area is plotted.Under the experimental conditions, the superimposed explosion wave stress of the blast hole wall is approximately 1.73 times the singleexplosionwave incident stress. The results of the model test and numerical simulations reveal that large-scale radial fracture cracks weregenerated on the blast hole wall in the superimposed area, and the width of the crack increased. The width of the large-scale radial fracturecracks formed by a strong impact is approximately 5% of the blast hole length. According to the characteristics of blast hole wallcompression, the mean peak pressures of the strongly superimposed area are approximately 1.48 and 1.84 times those of the weakly superimposedand nonsuperimposed areas, respectively.

Fractal characteristics of surface crack evolution in the process of gas-containing coal extrusion

In this paper, simulated experiment device of coal and gas outburst was employed to perform the experiment on gas-containing coal extrusion. In the experiment, coal surface cracks were observed with a high-speed camera and then the images were processed by sketch. Based on the above description, the paper studied the fractal dimension values from different positions of coal surface as well as their changing laws with time. The results show that there is a growing parabola trend of crack dimension value in the process of coal extrusion. Accordingly, we drew the conclusion that extruded coal crack evolution is a process of fractal dimension value increase. On the basis of fractal dimension values taken from different parts of coal masses, a fractal dimension of the contour map was drawn. Thus, it is clear that the contour map involves different crack fractal dimension values from different positions. To be specific, where there are complicated force and violent movement in coal mass, there are higher fractal dimension values, i.e., the further the middle of observation surface is from the exit of coal mass, and the lower the fractal dimension value is. In line with fractal geometry and energy theory of coal and gas outburst, this study presents the relation between fractal dimension and energy in the process of extruding. In conclusion, the evolution of crack fractal dimension value can signify that of energy, which has laid a solid foundation for the quantification research on the mechanism of gas-containing coal extrusion.

Development of X-FEM methodology and study on mixed-mode crack propagation

The extended finite element method（X-FEM） is a novel numerical methodology with a great potential for using in multi-scale computation and multi-phase coupling problems. The algorithm is discussed and a program is developed based on X-FEM for simulating mixed-mode crack propagation. The maximum circumferential stress criterion and interaction integral are deduced. Some numerical results are compared with the experimental data to prove the capability and efficiency of the algorithm and the program. Numerical analyses of sub-interfacial crack growth in bi-materials give a clear description of the effiect on fracture made by interface and loading condition.

SIF-based fracture criterion for interface cracks

The complex stress intensity factor K governing the stress field of an interface crack tip may be split into two parts, i.e.,■ and s^（-iε）, so that K = ■ s^（-iε）, s is a characteristic length and ε is the oscillatory index. ■ has the same dimension as the classical stress intensity factor and characterizes the interface crack tip field. That means a criterion for interface cracks may be formulated directly with■, as Irwin（ASME J. Appl. Mech. 24：361–364, 1957） did in 1957 for the classical fracture mechanics. Then, for an interface crack,it is demonstrated that the quasi Mode I and Mode II tip fields can be defined and distinguished from the coupled mode tip fields. Built upon SIF-based fracture criteria for quasi Mode I and Mode II, the stress intensity factor（SIF）-based fracture criterion for mixed mode interface cracks is proposed and validated against existing experimental results.

Crack buckling in soft gels under compression

Recent interest in designing soft gels with high fracture toughness has called for simple and robust methods to test fracture behavior. The conventional method of applying tension to a gel sample suffers from a difficulty of sample gripping. In this paper, we study a possible fracture mechanism of soft gels under uni-axial compression. We show that the surfaces of a pre-existing crack, oriented parallel to the loading axis, can buckle at a critical compressive stress. This buckling instability can open the crack surfaces and cre- ate highly concentrated stress fields near the crack tip, which can lead to crack growth. We show that the onset of crack buckling can be deduced by a dimensional argument com- bined with an analysis to determine the critical compression needed to induce surface instabilities of an elastic half space. The critical compression for buckling was verified for a neo- Hookean material model using finite element simulations.

CRACK TIP STRUCTURE AND DUCTILE-BRITTLE TRANSITION IN BULK SINGLE CRYSTALS

Based on the crack tip structure a new model of ductile -brittle transition was proposed. Using this new model we calculated the dependence of the transition temperature-strain rate over a wide range of strain rate. Finally the significance of this new model is discussed in detail.

COMPUTER SIMULATED ANALYSES ON DEFORMATION AND FRACTURE OF NON-CRACKED AND CRACKED SPECIMENS

A unified damage and fracture model,the combinatory work density model,which is suitable for ei- ther non-cracked body or cracked body has been suggested.In the present paper,the deformation and fracture of the two kinds of tensile spceimen and TPB specimen made of 40Cr steel have been simulated by using the new mod- el together with the large èlastic-plastic deformation finite element method.The results give a good picture of the whole deformation and fracture processes of the specimens in experiments;especially,the results on the TPB specimen can be used to obtain the relationship between load and displacement at the loading point P-Δ,and between crack ex- tension and displacement at the loading point Δα-Δ,the resistance curve J_R-Δa and the fracture toughness J_(IC).All the results are in remarkable agreement with those obtained by experiments.Therefore the model suggested here can be used to simulate crack initiation and propagation in non-cracked body and fracture initiation and crack stable propa- gation in cracked body.

Correlative Factors Research of Tibial Shaft Fracture with Ipsilateral Posterior Malleolar Crack Fractures

Background:Missed diagnosis rate of spiral tibia shaft fracture with posterior malleolus crack fracture(PMCF)is high in the clinical.However,the mechanism and related factors of fracture are still unclear.Moreover,PMCF has been observed in other types of tibial shaft fractures.Objective:To explore the correlative factors of tibial shaft fracture with ipsilateral PMCF,decrease the rate of clinical missed diagnosis,strengthen the effective fixation of PMCF,and reduce the incidence of traumatic arthritis.Methods:From September 2014 to May 2019,we collected 137 tibiofibular fracture.Only 68 cases involved in ankle joint CT examination and were retrospectively analyzed.The patients were divided into posterior malleolus group(30 cases)and non-posterior malleolus group(38 cases)according to whether come up PMCF or not.The posterior malleolus group contained 24 males and 6 females,27-77(47.57±11.79)years old,the non-posterior malleolus group contained 23 males and 15 females,18-85(48.71±13.84)years old.The gender,age,location,fibula fracture and tibial shaft fracture classification were observed for univariate and multivariate analysis.Results:The probability of PMCF was higher with right tibial shaft fracture(OR=3.6995%CI:1.13-12.08 P<0.05);the probability of PMCF following distal fibular fracture was higher than that without fibular fracture(OR=11.3695%CI:1.72-75.05 P<0.05);the probability of PMCF with type A tibial shaft fracture was higher than type C(OR=4.8295%CI:1.19-19.58 P<0.05).Conclusion:Right type A tibial shaft fracture accompanied by distal fibular fracture are very important factors related to PMCF,which needs highly attention to avoid clinical missed diagnosis.

Effect of thermal aging on the fatigue crack growth behavior of cast duplex stainless steels

The effect of thermal aging on the fatigue crack growth（FCG） behavior of Z3CN20？09M cast duplex stainless steel with low ferrite content was investigated in this study. The crack surfaces and crack growth paths were analyzed to clarify the FCG mechanisms. The microstructure and micromechanical properties before and after thermal aging were also studied. Spinodal decomposition in the aged ferrite phase led to an increase in the hardness and a decrease in the plastic deformation capacity, whereas the hardness and plastic deformation capacity of the austenite phase were almost unchanged after thermal aging. The aged material exhibited a better FCG resistance than the unaged material in the near-threshold regime because of the increased roughness-induced crack closure associated with the tortuous crack path and rougher fracture surface; however, the tendency was reversed in the Paris regime because of the cleavage fracture in the aged ferrite phases.

Analysis of the electromechanical behavior of ferroelectric ceramics based on a nonlinear finite element model

A nonlinear finite element （FE） model based on domain switching was proposed to study the electromechanical behavior of ferroelectric ceramics. The incremental FE formulation was improved to avoid any calculation instability. The problems of mesh sensitivity and convergence, and the efficiency of the proposed nonlinear FE technique have been assessed to illustrate the versatility and potential accuracy of the said technique. The nonlinear electromechanical behavior, such as the hysteresis loops and butterfly curves, of ferroelectric ceramics subjected to both a uniform electric field and a point electric potential has been studied numerically. The results obtained are in good agreement with those of the corresponding theoretical and experimental analyses. Furthermore, the electromechanical coupling fields near （a） the boundary of a circular hole, （b） the boundary of an elliptic hole and （c） the tip of a crack, have been analyzed using the proposed nonlinear finite element method （FEM）. The proposed nonlinear electromechanically coupled FEM is useful for the analysis of domain switching, deformation and fracture of ferroelectric ceramics.

A New 3D Meso-mechanical Modeling Method of Coral Aggregate Concrete Considering Interface Characteristics

On the basis of the three-dimensional(3D)random aggregate&mortar two-phase mesoscale finite element model,C++programming was used to identify the node position information of the interface between the aggregate and mortar elements.The nodes were discretized at this position and the zero-thickness cohesive elements were inserted.After that,the crack energy release rate fracture criterion based on the fracture mechanics theory was assigned to the failure criterion of the interface transition zone(ITZ)elements.Finally,the three-phase mesomechanical model based on the combined finite discrete element method(FDEM)was constructed.Based on this model,the meso-crack extension and macro-mechanical behaviour of coral aggregate concrete(CAC)under uniaxial compression were successfully simulated.The results demonstrated that the meso-mechanical model based on FDEM has excellent applicability to simulate the compressive properties of CAC.

Effect of shot peening on hydrogen embrittlement of high strength steel

The effect of shot peening（SP） on hydrogen embrittlement of high strength steel was investigated by electrochemical hydrogen charging, slow strain rate tensile tests, and hydrogen permeation tests. Microstructure observation, microhardness, and X-ray diffraction residual stress studies were also conducted on the steel. The results show that the shot peening specimens exhibit a higher resistance to hydrogen embrittlement in comparison with the no shot peening（NSP） specimens under the same hydrogen-charging current density. In addition, SP treatment sharply decreases the apparent hydrogen diffusivity and increases the subsurface hydrogen concentration. These findings are attributed to the changes in microstructure and compressive residual stress in the surface layer by SP. Scanning electron microscope fractographs reveal that the fracture surface of the NSP specimen exhibits the intergranular and quasi-cleavage mixed fracture modes, whereas the SP specimen shows only the quasi-cleavage fractures under the same hydrogen charging conditions, implying that the SP treatment delays the onset of intergranular fracture.

Incubation time fracture criterion for FEM simulations

The paper is discussing problems connected with embedment of the incubation time criterion for brittle fracture into finite element computational schemes. Incubation time fracture criterion is reviewed; practical questions of its numerical implementation are extensively discussed. Several examples of how the incubation time fracture criterion can be used as fracture condition in finite element computations are given. The examples include simulations of dynamic crack propagation and arrest, impact crater formation （i.e. fracture in initially intact media）, spall fracture in plates, propagation of cracks in pipelines. Applicability of the approach to model initiation, development and arrest of dynamic fracture is claimed.

Mechanical properties and failure characteristics of fractured sandstone with grouting and anchorage

Based on uniaxial compression experimental results on fractured sandstone with grouting and anchorage, we studied the strength and deformation properties, the failure model, crack formation and evolution laws of fractured sandstone under different conditions of anchorage. The experimental results show that the strength and elastic modulus of fractured sandstone with different fracture angles are significantly lower than those of intact sandstone. Compared with the fractured samples without anchorage,the peak strength, residual strength, peak and ultimate axial strain of fractured sandstone under different anchorage increase by 64.5–320.0%, 62.8–493.0%, and 31.6–181.4%, respectively. The number of bolts and degree of pre-stress has certain effects on the peak strength and failure model of fractured sandstone. The peak strength of fractured sandstone under different anchorage increases to some extent, and the failure model of fractured sandstone also transforms from tensile failure to tensile–shear mixed failure with the number of bolts. The pre-stress can restrain the formation and evolution process of tensile cracks, delay the failure process of fractured sandstone under anchorage and impel the transformation of failure model from brittle failure to plastic failure.

Computational Fracture Analysis of an AFM-Specimen under Mixed Mode Loading Conditions

Fracture processes in ship-building structures are in many cases of a 3-D character. A finite element （FE） model of an all fracture mode （AFM） specimen was built for the study of 3-D mixed mode crack fracture behavior including modes Ⅰ,Ⅱ, and Ⅲ. The stress intensity factors （SIFs） were calculated by the modified virtual crack closure integral （MVCCI） method, and the crack initiation angle assessment was based on a recently developed 3-D fracture criterion--the Richard criterion. It was shown that the FE model of the AFM-specimen is applicable for investigations under general mixed mode loading conditions, and the computational results of crack initiation angles are in agreement with some available experimental findings. Thus, the applicability of the FE model of the AFM-specimen for mixed mode loading conditions and the validity of the Richard criterion can be demonstrated.

APPLICATION OF SINTAP TO THE SAFETY ASSESSMENT OF X65 PIPELINE STEEL

The recently developed European flaw assessment procedure, structural integrity assessment procedure（SiNTAP） is applied to assessment for welded joints of the API 5L X65 pipeline steel with an assumed embedded flaw and surface flaw at the weld toe. As one of the basic input data, fracture toughness crack tip opening displacement （CTOD） tests are conducted at 0℃ and performed according to the requirements of the standard of BS7448. For the heat affected zone （HAZ） specimens, the microstructure observation is performed to insure that the tip of the crack is located in the coarse grain zone. The result explains the dispersity of the test values. In structural integrity assessment procedure-fracture assessment diagram（SINTAP FAD） method, the failure curves of welded joints at level 1 and 3 are derived from the tensile test results. The results of the assessment show that all assessment points are located within the failure lines of analysis level 1 and 3. So the welded joint of the pipeline is safe. This study laid the foundation of application of SINTAP to pipeline structure assessment.

Experimental and Numerical Study of Shear Fracture in Brittle Materials with Interference of Initial Double Cracks

A simultaneous experimental and numerical study of shear fracture of concrete-like materials is carried out using Brazilian disc specimens with initial double edge cracks and fourpoint bending beam specimens with double edge-notches.The interference effects of two cracks/notches are investigated through varied ligament angles and crack lengths.It is shown that shear fracturing paths change remarkably with the initial ligament angles and crack lengths.The cracked specimens are numerically simulated by an indirect boundary element method.A comparison between the numerical results and the experimental ones shows good agreement.

Semi-permeable Yoffe-type interfacial crack analysis in MEE composites based on the strip electro-magnetic polarization saturation model

Fracture analysis of a semi-permeable Yoffe-type interfacial crack propagating subsonically in magneto-electro-elastic （MEE） composites is presented based on the strip electro- magnetic polarization saturation （SEMPS） model. The electro-magnetic fields inside the crack are considered under the semi-permeable boundary condition. Nonlinear effects near the interfacial crack tip are represented by different electro-magnetic saturation zones. Utilizing the extended Stroh＇s method, we derive the moving dislocation densities as well as intensity factor and energy release rate for Yoffe-type MEE interracial crack. Numerical re- sults through an iterative approach are presented to show the characteristics of fracturedominant parameters with respect to propagation velocity and boundary condition category. The fracture-dominant parameters under the semi-permeable boundary condition are lower than those under the impermeable one, which implies that the electro-magnetic fields in the crack gap can retard the propagation of MEE interfacial crack.

Domain-switching embedded nonlinear electromechanical finite element method for ferroelectric ceramics

To simulate the nonlinear behavior of ferroelectric structures and devices under non-uniform electromechanical loadings,a domain-switching embedded electromechanical finite element method is developed in this paper.Following continuum assumption,the electromechanical behavior of each representative material point can be obtained by averaging the behavior of the local corresponding microstructure,e.g.42 domains used in this work.A new Double Gibbs free energy criterion for domain-switching is proposed to ensure the convergence and stability of the simulations on ferroelectrics under non-uniform field.Several computational examples are given to demonstrate that this nonlinear finite element method can yield reasonable and stable simulation results which can be used to explain some experimental results and assist the design of ferroelectric devices.