Analytical solution for the effective elastic properties of rocks with the tilted penny-shaped cracks in the transversely isotropic background

Seismic prediction of cracks is of great significance in many disciplines,for which the rock physics model is indispensable.However,up to now,multitudinous analytical models focus primarily on the cracked rock with the isotropic background,while the explicit model for the cracked rock with the anisotropic background is rarely investigated in spite of such case being often encountered in the earth.Hence,we first studied dependences of the crack opening displacement tensors on the crack dip angle in the coordinate systems formed by symmetry planes of the crack and the background anisotropy,respectively,by forty groups of numerical experiments.Based on the conclusion from the experiments,the analytical solution was derived for the effective elastic properties of the rock with the inclined penny-shaped cracks in the transversely isotropic background.Further,we comprehensively analyzed,according to the developed model,effects of the crack dip angle,background anisotropy,filling fluid and crack density on the effective elastic properties of the cracked rock.The analysis results indicate that the dip angle and background anisotropy can significantly either enhance or weaken the anisotropy degrees of the P-and SH-wave velocities,whereas they have relatively small effects on the SV-wave velocity anisotropy.Moreover,the filling fluid can increase the stiffness coefficients related to the compressional modulus by reducing crack compliance parameters,while its effects on shear coefficients depend on the crack dip angle.The increasing crack density reduces velocities of the dry rock,and decreasing rates of the velocities are affected by the crack dip angle.By comparing with exact numerical results and experimental data,it was demonstrated that the proposed model can achieve high-precision estimations of stiffness coefficients.Moreover,the assumption of the weakly anisotropic background results in the consistency between the proposed model and Hudson's published theory for the orthorhombic rock.

Mechanical properties of AISI 1045 ceramic coated materials by nano indentation and crack opening displacement method

Abstract： An effective approach was conducted for estimating fracture toughness using the crack opening displacement （COD） method for plasma enhanced chemical vapor deposition （PECVD） coating materials. For this evaluation, an elastoplastic analysis was used to estimate critical COD values for single edge notched bending （SENB） specimens. The relationship between fracture toughness （Kic） and critical COD for SENB specimens was obtained. Microstructure of the interface between AleO3-TiO2 composite ceramic coatings and AISI 1045 steel substrates was studied by using scanning electron microscope （SEM）. Chemical compositions were clarified by energy-dispersive X-ray spectroscopy （EDS）. The results show that the interface between of Al203-TiO2 and substrate has mechanical combining. The nanohardness of the coatings can reach 1 200 GPa examined by nanoindentation. The Klc was calculated according to this relationship from critical COD. The bending process produces a significant relationship of COD independent of the axial force applied. Fractographic analysis was conducted to determine the crack length. From the physical analysis of nanoindentation curves, the elastic modulus of 1045/AI2O3-TiO2 is 180 GPa for the 50 μm film. The highest value of fracture toughness for 1045/A1203-TiO2-250 μm is 348 MPa·mv2.

Critical crack tip opening displacement of different strength concrete

Critical crack tip opening displacement (CTODc) of concrete using experimental and analytical evaluation with seven different compressive strengths ranging from 30 up to 150 MPa was studied based on two types of fracture tests:three-point bending (TPB) and wedge splitting (WS).In the tests,the values of CTODc were experimentally recorded using a novel technique,in which fiber Bragg grating (FBG) sensors were used,and two traditional techniques,in which strain gauges and clip gauges were deployed.The values of CTODc of tested concrete were also predicted using two existing analytical formulae proposed by JENQ & SHAH and XU,respectively.It is found that the values of CTODc obtained by both experimental measurements and analytical formulae exhibit a negligible variation as the compressive strength of concrete increases,and the test geometry adopted has little impact on the value of CTODc.Regarding the experimental measurement of CTODc,the clip gauge method generally leads to a larger value of CTODc and shows a more significant scatter as compared with the other two methods,while the strain gauge method leads to a slightly lower CTODc as compared with the FBG sensor method.The analytical formula proposed by JENQ and SHAH is found to generally lead to an overestimation,while the analytical formula proposed by XU shows a good accuracy.

Failure mechanisms and destruction characteristics of cemented coal gangue backfill under compression effect of non-uniform load

Backfill mining is one of the most important technical means for controlling strata movement and reducing surface subsidence and environmental damage during exploitation of underground coal resources. Ensuring the stability of the backfill bodies is the primary prerequisite for maintaining the safety of the backfilling working face, and the loading characteristics of backfill are closely related to the deformation and subsidence of the roof. Elastic thin plate model was used to explore the non-uniform subsidence law of the roof, and then the non-uniform distribution characteristics of backfill bodies’ load were revealed. Through a self-developed non-uniform loading device combined with acoustic emission (AE) and digital image correlation (DIC) monitoring technology, the synergistic dynamic evolution law of the bearing capacity, apparent crack, and internal fracture of cemented coal gangue backfills (CCGBs) under loads with different degrees of non-uniformity was deeply explored. The results showed that: 1) The uniaxial compressive strength (UCS) of CCGB increased and then decreased with an increase in the degree of non-uniformity of load (DNL). About 40% of DNL was the inflection point of DNL-UCS curve and when DNL exceeded 40%, the strength decreased in a cliff-like manner;2) A positive correlation was observed between the AE ringing count and UCS during the loading process of the specimen, which was manifested by a higher AE ringing count of the high-strength specimen. 3) Shear cracks gradually increased and failure mode of specimens gradually changed from “X” type dominated by tension cracks to inverted “Y” type dominated by shear cracks with an increase in DNL, and the crack opening displacement at the peak stress decreased and then increased. The crack opening displacement at 40% of the DNL was the smallest. This was consistent with the judgment of crack size based on the AE b-value, i. e., it showed the typical characteristics of “small b-value-large crack and large b-value-small crack”. The research results are of significance for preventing the instability and failure of backfill.

CLOSED-FORM WEIGHT FUNCTION FOR EDGE CRACK PROBLEMS

An accurate crack surface opening displacement formula for edge cracks in finite bodies has been derived,and was subsequently used to develop wide-range closed-form weight functions for edge cracks. The analytical weight function is general and can be easily used to determine accurate stress intensity factors and crack surface opening displacements for arbitrarily loaded edge cracks with high efficiency.Examples have been given for problems of edge crack(s)with a number of typical load conditions for illustration.

Strip-coalesced interior zone model for two unequal collinear cracks weakening piezoelectric media

In this paper, a mathematical strip-saturation model is proposed for a poled transversely isotropic piezoelectric plate weakened by two impermeable unequal-collinear hairline straight cracks. Remotely applied in-plane unidirectional electromechanical loads open the cracks in mode-I such that the saturation zone developed at the interior tips of cracks gets coalesced. The developed saturation zones are arrested by distributing over their rims in-plane normal cohesive electrical displacement. The problem is solved using the Stroh formalism and the complex variable technique. The expressions are derived for the stress intensity factors （SIFs）, the lengths of the saturation zones developed, the crack opening displacement （COD）, and the energy release rate. An illustrative numerical case study is presented for the poled PZT-5H ceramic to investigate the effect of prescribed electromechanical loads on parameters affecting crack arrest. Also, the effect of different lengths of cracks on the SIFs and the local energy release rate （LERR） has been studied. The results obtained are graphically presented and analyzed.

Modeling of Fatigue Crack Growth Closure Considering the Integrative Effect of Cyclic Stress Ratio,Specimen Thickness and Poisson's Ratio

Key components of large structures in aeronautics industry are required to be made light and have long enough fatigue lives.It is of vital importance to estimate the fatigue life of these structures accurately.Since the FCG process is affected by various factors,no universal model exists due to the complexity of the mechanisms.Most of the existing models are obtained by fitting the experimental data and could hardly describe the integrative effect of most existing factors simultaneously.In order to account for the integrative effect of specimen parameters,material property and loading conditions on FCG process,a new model named integrative influence factor model（IIF） is proposed based on the plasticity-induced crack closure theory.Accordingly to the predictions of crack opening ratio（γ） and effective stress intensity factor range ratio（U） with different material under various loading conditions,predictions of γ and U by the IIF model are completely identical to the theoretical results from the plane stress state to the plane strain state when Poisson＇s ratio equals 1/3.When Poisson＇s ratio equals 0.3,predictions of γ and U by the IIF model are larger than the predictions by the existing model,and more close to the theoretical results.In addition,it describes the influence of R ratios on γ and U effectively in the whole region from-1.0 to 1.0.Moreover,several sets of test data of FCG rates in 5 kinds of aluminum alloys with various specimen thicknesses under different loading conditions are used to validate the IIF model,most of the test data are situated on the predicted curves or between the two curves that represent the specimen with different thicknesses under the same stress ratio.Some of the test data slightly departure from the predictions by the IIF model due to the surface roughness and errors in measurement.Besides,based on the analysis of the physical rule of crack opening ratios,a relative thickness of specimen is defined to describe the influence of material property,specimen thickness and so forth on FCG characteristics conveniently.In conclusion,the relative thickness of specimen simplifies the expression of FCG characteristic and provides a general parameter to analyze the fatigue characteristics of different materials with various thicknesses under different loading conditions.The IIF model describes the integrative effect of existing influence factors explicitly and quantitatively,and provides a helpful tool for fatigue property estimation of practical component and experiment design.

Solution of stress intensity factors of multiple cracks in plane elasticity with eigen COD formulation of boundary integral equation

The concept of eigen crack opening displacement （COD） can be defined as the COD of a crack in infinite plate under the tractions acting on the crack surface. By introducing this concept, the eigen COD formulation of boundary integral equation is proposed in this paper, together with the solution procedures for multiple crack problems in plane elasticity. With the proposed approach, the multiple crack problems can be solved with the conventional displacement discontinuity boundary integral equations in an iterative fashion with a small size of system matrix as that in the numerical Green’s function （NGF） approach but without the trouble to determine the complementary solutions since the standard boundary element discretization on the crack surface is no longer required with the proposed approach. Some numerical examples computing the stress intensity factors are presented and compared with those in literature to show the accuracy and the effectiveness of the proposed approach.

Digital image correlation based high-speed crack tip locatingmethod and its application

Based on a digital image correlation（DIC）method with the measurements of a high speed crack＇s displacement and strain fields,a technique to accurately and automatically locate its crack tip has been developed.The crack tip is identified by finding the abrupt jump on the opening（or dislocation）curve of a point on the trace of the crack propagation,while the opening is measured through a virtual extensometer technique and the abrupt jump is identified by finding the peak on the curve.The method was verified using a specially designed experiment and applied to measure the critical crack tip opening angle of a rock sample.Because the involvement of analytical models has been avoided and then the good performance could be ensured for low resolution speckle images,this technique is expected to be very useful in the quantitative study of high speed cracks in experiments using high speed cameras.

A Straightforward Direct Traction Boundary Integral Method for Two-Dimensional Crack Problems Simulation of Linear Elastic Materials

This paper presents a direct traction boundary integral equation method(DTBIEM)for two-dimensional crack problems of materials.The traction boundary integral equation was collocated on both the external boundary and either side of the crack surfaces.The displacements and tractions were used as unknowns on the external boundary,while the relative crack opening displacement(RCOD)was chosen as unknowns on either side of crack surfaces to keep the single-domain merit.Only one side of the crack surfaces was concerned and needed to be discretized,thus the proposed method resulted in a smaller system of algebraic equations compared with the dual boundary element method(DBEM).A new set of crack-tip shape functions was constructed to represent the strain field singularity exactly,and the SIFs were evaluated by the extrapolation of the RCOD.Numerical examples for both straight and curved cracks are given to validate the accuracy and efficiency of the presented method.

PERIODICAL INTERFACIAL CRACKS IN ANISOTROPIC ELASTOPLASTIC MEDIA

By using Fourier transformation the boundary problem of periodical interfacial cracks in anisotropic elastoplastic bimaterial was transformed into a set of dual integral equations and then it was further reduced by means of definite integral transformation into a group of singular equations. Closed form of its solution was obtained and three corresponding problems of isotropic bimaterial, of a single anisotropic material and of a bimaterial of isotropy- anisotropy were treated as the specific cases. The plastic zone length of the crack tip and crack openning displacement ( COD) decline as the smaller yield limit of the two bonded materials rises, and they were also determined by crack length and the space between two neighboring cracks . In addition , COD also relates it with moduli of the materials .

Engineering critical assessment on welding flaws of the X65 UOE girth weld based on CRACKWISE

Focused on the girth weld of a φ762 x 28.6 mm （X65） offshore UOE pipe, the fracture toughness properties at different weld joint positions were tested. Meanwhile,the largest size of different cracks under service condition was calculated by using CRACKWISE. The service life of the UOE pipe with postulated crack-like flaws was calculated by considering various fatigue factors, such as vibrations caused by ocean current and fluctuations of inner pressure. The assessment process and results can be used to direct the repair of weld flaws laid in pipes or to assess the reliability of in-serving pipes.

Numerical Comparison Research on the Solution of Stress Intensity Factors of Multiple Crack Problems

A newly developed approach without crack surface discretization for modeling 2D solids with large number of cracks in linear elastic fracture mechanics is proposed with the eigen crack opening displacement (COD) boundary integral equations in this paper. The eigen COD is defined as a crack in an infinite domain under fictitious traction acting on the crack surface. Respect to the computational accuracies and efficiencies, the multiple crack problems in finite and infinite plates are solved and compared numerically using three different kinds of boundary integral equations (BIEs): 1) the dual BIEs require crack surface discretization;2) the BIEs with numerical Green’s functions (NGF) without crack surface discretization, but have to solve a complementary matrix;3) the eigen crack opening displacement (COD) BIEs in the present paper. With the concept of eigen COD, the multiple crack problems can be solved by using a conventional displacement discontinuity boundary integral equation in an iterative fashion with a small size of system matrix as that in the NGF approach, but without troubles to determine the complementary matrix. Solution of the stress intensity factors of multiple crack problems is solved and compared in some numerical examples using the above three computational algorithms. Numerical results clearly demonstrate the numerical models of eigen COD BIEs have much higher efficiency, providing a newly numerical technique for multiple crack problems. Not only the accuracy and efficiency of computation can be guaranteed, but also the overall properties and local details can be obtained. In conclusion, the numerical models of eigen COD BIEs realize the simulations for multiple crack problems with large quantity of cracks.

THE ANALYSES OF THE THREE-DIMENSIONAL STRESS STRUCTURE NEAR THE CRACK TIP OF MODE I CT SPECIMENS IN ELASTICPLASTIC STATE(Ⅰ)--THE ANALYSES OF CONSTRAINT PARAMETERS AND FRACTURE PARAMETERS

In the present paper,three dimensional analyses of some general constraint parameters and fracture parameters near the crack tip of Mode I CT specimens in two different thicknesses are carried out by employing ADINA program.The results reveal that the constraints along the thickness direction are obviously separated into two parts:the keeping similar high constraint field(Z_(1))and rapid reducing constraints one(Z_(2)).The two fields are experimentally confiremed to correspond to the smooth region and the shear lip on the fracture face respectively.So the three dimensional stress structure of Mode I specimens can be derived through discussing the two fields respectively.The distribution of the Crack Tip Opening Displacement(CTOD)along the thickness direction and the three dimensional distribution of the void growth ratio(V_(g))near the crack tip are also obtained.The two fracture parameters are in similar trends along the thickness direction,and both of them can reflect the effect of thickness and that of the loading level to a certain degree.

ANTI-PLANE SHEAR WAVES SCATTERING FROM A PARTIALLY DEBONDED MAGNETO-ELECTRO-ELASTIC CIRCULAR CYLINDRICAL INHOMOGENEITY

This article presents an analysis of the scattering of anti-plane shear waves from a single cylindrical inhomogeneity partially bonded to an unbounded magneto-electro-elastic matrix. The magneto-electric permeable boundary conditions are adopted. The crack opening displacement is represented by Chebyshev polynomials and a system of equations is derived and solved for the unknown coefficients. Some examples are calculated and the results are illustrated. The results show that the COD increases when the piezomagnetic coefficient of the inhomogeneity bonded to the piezoelectric matrix becomes larger, and that the COD decreases when the piezomagnetic coefficient of the matrix with the piezoelectric inhomogeneity increases.

Analysis of crack opening stresses for center- and edge-crack tension specimens

Accurate determination of crack opening stress is of central importance to fatigue crack growth analysis and life prediction based on the crack-closure model. This paper studies the crack opening behavior for center- and edge-crack tension specimens. It is found that the crack opening stress is affected by the crack tip element. By taking the crack tip element into account, a modified crack opening stress equation is given for the center-crack tension specimen. Crack surface displace- ment equations for an edge crack in a semi-infinite plate under remote uniform tension and partially distributed pressure are derived by using the weight function method. Based on these displacements, a crack opening stress equation for an edge crack in a semi-infinite plate under uniform tension has been developed. The study shows that the crack opening stress is geometry-dependent, and the weight function method provides an effective and reliable tool to deal with such geometry depen- dence.

An Improved Test Method and Numerical Analysis for Crack Opening Resistance of FRC Round Determinate Panels

This paper presents a numerical analysis of a material test program investigating the crack opening properties of fibre reinforced concrete （FRC） round determinate panels （RDP）. The objective of this research is to set up a modified RDP test method to improve the current ASTM C1550 test method for FRC and fibre reinforced shotcrete （FRS） composite. By this test method, light small-diameter panels are applied and more crack information can be obtained by a new rotation angle measuring technique. It is shown that this modified test method can be used to effectively evaluate the crack opening resistance of FRC. The finite element analysis was then performed to clarify the crack propaga- tions and failure mode of FRC RDP panels. It helps establish a reasonable theoretical method to predict the structural response of RDP. combining with this modified testing technioue.

Simulating hydraulic fracturing processes in laboratory-scale geological media using three-dimensional TOUGH-RBSN

In this context, recent developments in the coupled three-dimensional(3 D) hydro-mechanical(HM)simulation tool TOUGH-RBSN are presented. This tool is used to model hydraulic fracture in geological media, as observed in laboratory-scale tests. The TOUGH-RBSN simulator is based on the effective linking of two numerical methods: TOUGH2, a finite volume method for simulating mass transport within a permeable medium; and a lattice model based on the rigid-body-spring network(RBSN) concept. The method relies on a Voronoi-based discretization technique that can represent fracture development within a permeable rock matrix. The simulator provides two-way coupling of HM processes, including fluid pressure-induced fracture and fracture-assisted flow. We first present the basic capabilities of the modeling approach using two example applications, i.e. permeability evolution under compression deformation, and analyses of a static fracturing simulation. Thereafter, the model is used to simulate laboratory tests of hydraulic fracturing in granite. In most respects, the simulation results meet expectations with respect to permeability evolution and fracturing patterns. It can be seen that the evolution of injection pressure associated with the simulated fracture developments is strongly affected by fluid viscosity.

An information-based rough set approach to critical engineering factor identification

In order to analyze the main critical engineering factors, an information-based rough set approach that considers conditional information entropy as a measurement of information has been developed. An algorithm for continuous attribute discretization based on conditional information entropy and an algorithm for rule extraction considering the supports of rules are proposed. The initial decision system is established by collecting enough monitoring data. Then, the continuous attributes are discretized, and the condition attributes are reduced. Finally, the rules that indicate the action law of the main factors are extracted and the results are explained. By applying this approach to a crack in an arch gravity dam, it can be concluded that the water level and the temperature are the main factors affecting the crack opening, and there is a negative correlation between the crack opening and the temperature. This conclusion corresponds with the observation that cracks in most concrete dams are influenced mainly by water level and temperature, and the influence of temperature is more evident.

STUDY ON THE RELATIONSHIP BETWEEN J-INTEGRAL AND COD

J-integral arid crack opening displacement 8 （COD） were important parameters for characterization of fractures in engineering materials. The relationship between J-integral arid COD has been investigated for a long time and was generally represented as J = m σy, δ, where σy is the yield strength arid m is a function of specimen geometry arid material properties. To determine the value of m, extensive studies and experiments have been performed. The method that used the fracture-surface topography analysis （ FRASTA ） for determining J-integral from fracture surfaces of materials was introduced. On the basis of the relationship between COD arid fracture surface average profile, the relationship between J-integral arid COD was deduced arid compared with the generally used equation. The method was experimentally confirmed to be able to provide a new way to determine the relationship between J-integral and COD.