Numerical simulation study on shear resistance of anchorage joints considering tensile-shear fracture criterion of 2G-NPR bolt

2G-NPR bolt (the 2nd generation Negative Poisson’s Ratio bolt) is a new type of bolt with high strength, high toughness and no yield platform. It has signifcant efects on improving the shear strength of jointed rock mass and controlling the stability of surrounding rock. To achieve an accurate simulation of bolted joint shear tests, we have studied a numerical simulation method that takes into account the 2G-NPR bolt's tensile–shear fracture criterion. Firstly, the indoor experimental study on the tensile–shear mechanical properties of 2G-NPR bolt is carried out to explore its mechanical properties under diferent tensile–shear angles, and the fracture criterion of 2G-NPR bolt considering the tensile–shear angle is established. Then, a three-dimensional numerical simulation method considering the tensile–shear mechanical constitutive and fracture criterion of 2G-NPR bolt, the elastoplastic mechanical behavior of surrounding rock and the damage and deterioration of grouting body is proposed. The feasibility and accuracy of the method are verifed by comparing with the indoor shear test results of 2G-NPR bolt anchorage joints. Finally, based on the numerical simulation results, the deformation and stress of the bolt, the distribution of the plastic zone of the rock mass, the stress distribution and the damage of the grouting body are analyzed in detail. The research results can provide a good reference value for the practical engineering application and shear mechanical performance analysis of 2G-NPR bolt.

A DUCTILE FRACTURE CRITERION BASED ON THREEDIMENSIONAL VOID MODEL

A fracture criterion derived from a microscopic point of view is proposed and has proved to be effective in the analysis of uniaxial tension. On the one hand, a method of predicting a ductile fracture is proposed using a three-dimensional void model and the assumption of velocity discontinuity. The relationship between the void volume fraction and the critical strain to fracture, calculated with the help of the new model, shows the same tendency as that obtained from the modified Thomason model. On the other hand, the mechanical and metallographic analyses of the uniaxial tension experiment are performed using four kinds of carbon steel. The relationship between the void volume fraction and the critical strain to fracture, calculated from the new model, agrees better with the result obtained from the experiment, rather than that calculated by the modified Thomason model, which confirms the validity of the ductile fracture criterion based on the three-dimensional void model.

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.

Comparison of Modified Mohr-Coulomb Model and Bai-Wierzbicki Model for Constructing 3D Ductile Fracture Envelope of AA6063

Ductile fracture of metal often occurs in the plastic forming process of parts.The establishment of ductile fracture criterion can effectively guide the selection of process parameters and avoid ductile fracture of parts during machining.The 3D ductile fracture envelope of AA6063-T6 was developed to predict and prevent its fracture.Smooth round bar tension tests were performed to characterize the flow stress,and a series of experiments were conducted to characterize the ductile fracture firstly,such as notched round bar tension tests,compression tests and torsion tests.These tests cover a wide range of stress triaxiality(ST)and Lode parameter(LP)to calibrate the ductile fracture criterion.Plasticity modeling was performed,and the predicted results were compared with corresponding experimental data to verify the plasticity model after these experiments.Then the relationship between ductile fracture strain and ST with LP was constructed using the modified Mohr-Coulomb(MMC)model and Bai-Wierzbicki(BW)model to develop the 3D ductile fracture envelope.Finally,two ductile damage models were proposed based on the 3D fracture envelope of AA6063.Through the comparison of the two models,it was found that BW model had better fitting effect,and the sum of squares of residual error of BW model was 0.9901.The two models had relatively large errors in predicting the fracture strain of SRB tensile test and torsion test,but both of the predicting error of both two models were within the acceptable range of 15%.In the process of finite element simulation,the evolution process of ductile fracture can be well simulated by the two models.However,BW model can predict the location of fracture more accurately than MMC model.

Theoretical and experimental study of 3-D initial fracture and its significance to faulting

The experimental results of 3-D fracture under compression are introduced in brief and the theory of stress criterion of 3-D fracture is studied. Methods to imitate initial fractures are developed. It is pointed that there are important defects in the extreme value （EV） method ever proposed by Palaniswamy and Knauss. The major defect lies in that only two Euler angles （2EA） are considered, but another one is neglected. If the variation of all the three Euler angles （3EA） are considered, one can get better result which is consistent with the observation of faulting that extends on curved surfaces but not on planes. The method of evaluating maximal normal stress direction vector （NSDV） is proposed and further proved to be equivalent to the 3EA method. It is proved that the NSDV method can be further optimized to the method of composition of the first principal differential plane （CFPDP）. The results from CFPDP method can fit the curved surfaces of initial growth observed in the experiments of 3-D fracture. The CFPDP method can also be used to interpret the 3-D fractures of the slipping section between the asperities in the buried fault plane that is modeled as ellipse crack. The results of 3-D fracture can be applied to interpreting the related problems of faulting including the mechanism of a lot of shatter rocks with different dimensions, the cause of earthquakes occurred at the edge of plate under low shear stress, and the mechanism of anisotropy caused by the extensive dilatancy anisotropy （EDA） cracks.

Energy consumption for creep fracture of metallic materials

The energy conservation law is applied to formulate the ductile and brittle creep fracture criterion for metallic materials. The criterion contains a summary of heat and latent energies. Assuming that the heat energy is given out so it has no effect on the fracture process, the ductile creep fracture criterion is simplified. To take into account the evaluation of the damage state of materials the compressibility condition is introduced and the brittle creep fracture law is formulated.

Effect of Fracture Criteria on Taylor Impact Fracture

Taylor impact test on aluminum alloy 7A04-T6 presented in a previous study has been numerically evaluated using the finite element code ABAQUS/Explicit in this paper. The influence of fracture criterion in numerical simulations of the deformation and fracture behavior of Taylor rod has been studied. The following fracture criterions, included a modified version of Johnson-Cook, the Cockcroft-Latham(C-L), the constant fracture strain, the maximum shear stress and the maximum principle stress fracture models, have been used in this paper. Model constants for each criterion are calibrated from material tests. The modified version of Johnson-Cook fracture criterion with the stress triaxiality cut-off idea used can give good prediction to the Taylor impact fracture behavior. However, the C-L fracture criterion with only a constant required for calibrating by a simple test is found to give reasonable predictions. Unfortunately, the other three criteria are not able to repeat the experimentally obtained fracture behavior. The result indicates that the stress triaxiality cut-off idea is necessary to predict the Taylor impact fracture.

ON PROBLEMS OF FRACTURE OF MATERIALS IN COMPRESSION ALONG TWO INTERNAL PARALLEL CRACKS

The fracture of materials under the action of compressive forces, directed along cracks which are parallel in plane can,lot be described within the framework of the linear fracture mechanics. The criteria of fracture of the Griffith-Irvin or COC type, used in classical linear fracture mechanics, are not applicable in this problem, since these forces have no influence on stress intensity coefficients and on values of cracks opening([1, 2]). The problems of such a class may be described only by using new approaches. One of possible approaches is presented by the first author, which involves using linearized relations, derived from exact non-linear equations of deformable solid body mechanics([3, 4, 5]). It should be remarked here that this approach has been widely used in problems of deformable bodies stability. As a criterion of the initiation of fracture the criterion of local instability near defects of the crack type is used. In these cases the process of loss of stability initiates the fracture process.

ANALYSIS OF A CRACK IN A FUNCTIONALLY GRADED STRIP WITH A POWER FORM SHEAR MODULUS

The plane strain problem of a crack in a functionally graded strip with a power form shear modulus is studied. The governing equation in terms of Airy＇s stress function is solved exactly by means of Fourier transform. The mixed boundary problem is then reduced to a system of singular integral equations and is solved numerically to obtain the stress intensity factor at crack-tip. The maximum circumferential stress criterion and the strain energy density criterion are both employed to predict the direction of crack initiation. Numerical examples are given to show the influence of the material gradation models and the crack sizes on the mode-I and mode-II stress intensity factors. The dependence of the critical kink-angle on the crack size is examined and it is found that the crack kink-angle decreases with the increase of the normalized crack length, indicating that a longer crack tends to follow the original crack-line while it is much easier for a shorter crack to deviate from the original crack-line.

RANDOM-FUZZY MODEL FOR THE DUCTILE/BRITTLE TRANSITION

Based on a large experimental sample, the influence of stress triaxiality and temperature have been studied on the growth of micro voids and the ductile/brittle transition (DBT) of 40Cr steel. Ductile and brittle fractures have been investigated simultaneously. The experiments show that the ductile fracture is controlled by the micro void mechanism,and the critical void growth ratio can be used as the criterion of ductile fracture. The brittle fracture is modeled by an embedded penny crack, and a micro stress intensity with a characteristic length can be used as the criterion of the brittle fracture. The DBT is the result of the competition of the two mechanisms. Which exists in the fracture of all specimens simultaneously. The distributions of model parameters were measured experimentally. A random model and a random-fuzzy model for DBT were presented. The comparison between the calculated and experimental results shows that the random-fuzzy model can model the DBT satisfactorily.

Strength evaluation of ultra-high strength steels and spot weld of automotive bodies

Cockcrofi-Latham fracture criterion was applied to predict the fracture of high strength steels. A Marciniak-type biaxial stretching test of the four classes of high strength steels was carried out to measure the material damage limit of Cockcrofi-Latham fracture criterion. Furthermore, in order to improve the simulation accuracy, the local anisotropic parameters depending on the plastic strain （strain dependent model of anisotropy ） were measured by digital image correlation method and incorporated into Hill＇ s anisotropic yield condition by authors. To confirm the validity of Cockcrofi-Latham fracture criterion, the uniaxial tensile tests based on JIS No. 5 tensile specimen were performed. The force-displacement history and fracture happening strokes were predicted with high accuracy. Then, Cockcrofi-Latham fracture criterion was applied to predict the failure of four types of spot welded joints. To simulate the local bending and warping deformations around the heataffected zone, the discrete Kirchhoff triangle element was adapted. FEM results for four classes of high strength steels and four types of spot welded joints had a good correlation with experimental ones.

Deformation and Densification Laws of Powder Cold Forging

The deformation and densification laws of preform upsetting and closed-die forging were researched based on experimental results of cold forging of deoxidized Fe powder sintering porous material under different initial conditions such as friction factor, ratio between height and diameter and relative density. The fracture limit criteria＂ for powder cold-forging upsetting and the limit strain curve were achieved. The effect of friction facto,, ratlt, between height and diameter and relative density on fracture strain limitation was emphatically analyzed. The limit process parameter curves for the deformation of upsetting were also confirmed. Laws of deformation, densification and density distribution for closed-die forging of powder perform during cold-forging were further analyzed and discussed with the help of experimental phase analysis. As a result, this experiment established theoretical foundations for the design of preform and die as well as optimization of technological process parameters.

Ductile Fracture Prediction of 316LN Stainless Steel in Hot Deformation Process

A ductile fracture criterion of 316LN stainless steel, combined with the plastic deformation capacity of ma- terial and the stress state dependent damages, was proposed to predict ductile fracture during hot deformation. To the end, tensile tests at high temperatures were first performed to investigate the fracture behavior of 316LN stain- less steel. The experimental results show the variation of the critical fracture strain as a function of temperature and strain rate. Second, the criterion was calibrated by using the upsetting tests and the corresponding numerical simula- tions. Finally, the proposed fracture criterion was validated by the designed tests and the corresponding finite ele- ment （FE） simulation. The results show that the criterion can successfully predict the onset of ductile fracture at ele- vated temperatures.

A UNIVERSAL CRACK EXTENSION CRITERION BASED ON THE EQUIVALENT STRESS GRADIENT:I.THEORY AND NUMERICAL VERIFICATION

In this contribution,the microscopic fracture mechanism and extension criterion for mixed type crack in ductile material under plane mixed mode loading are investigated in details.A universal extension criterion for the mixed type crack,i.e.the crack propagates along the direction of the maximum gradient of equivalent stress,is suggested.This new criterion is used to predict the propagation direction of mixed type crack,showing a good agreement with other theories for different types of mode mixity.Moreover,the numerical verification is also carried out for the case of an edge crack with different mixed mode loadings.Finally,a potential application to three-dimensional fracture in the ductile material induced by holes is also discussed.

A constitutive model coupling damage and material anisotropy for wide stress triaxiality

A constitutive model that can describe the damage evolution of anisotropic metal sheets during the complex forming processes which experience wide stress triaxiality history is essential to accurately predict the deformation and rupture behaviors of the processes.In this study,a modified Lemaitre damage criterion which couples with the anisotropic Barlat 89 yield function is established.The effects of stress triaxiality,Lode parameter and shear stress on damage accumulation are considered in the constitutive model.The model is numerically implemented and applied to fracture prediction in tensile tests with different stress triaxialities and a complex deformation process with wide stress triaxiality history.The good consistency of predictions and experiments indicates that the modified Lemaitre damage model has excellent fracture prediction ability.Finally,the accuracy of the model is analyzed and discussed.

ROCK INITIATION AND PROPAGATION SIMULATION UNDER COMPRESSION-SHEAR LOADING USING CONTINUOUS-DISCONTINUOUS CELLULAR AUTOMATON METHOD

A continuous-discontinuous cellular automaton method is developed for rock initiation and propagation simulations, in which the level set method, discontinuous enrichment shape functions and discontinuous cellular automaton are combined. No renmshing is needed for crack growth analysis, and all calculations are restricted to cells without an assembled global stiffness matrix. The frictional contact theory is employed to construct the contact model of normal pressure and tangential shear on crack surfaces. A discontinuous cellular automaton updating rule suitable for frictional contact of rock is proposed simultaneously with Newton＇s iteration method for nonlinear iteration. Besides, a comprehensive fracturing criterion for brittle rock under compression-shear loading is developed. The accuracy and effectivenesss of the proposed method is proved by numerical simulation.