Interacting mechanism and initiation prediction of multiple cracks

The maximum Mode Ⅰ and Mode Ⅱ stress intensity factors(SIFs), KI,kmax(θ) and KII,kmax(θ)(0°<θ<360°), of inclined parallel multi-crack varying with relative positions(including horizontal and vertical spacings) are calculated by the complex function and integration method to analyze their interacting mechanism and determine the strengthening and weakening zone of SIFs. The multi-crack initiation criterion is established based on the ratio of maximum tension-shear SIF to predict crack initiation angle, load, and mechanism. The results show that multi-crack always initiates in Mode Ⅰ and the vertical spacing is better not to be times of half crack-length for crack-arrest, which is in good agreement with test results of the red-sandstone cube specimens with three parallel cracks under uniaxial compression. This can prove the validity of the multi-crack initiation criterion.

Multi‑Crack Interaction and Influence on the Spherical Pressure Hull for a Deep‑Sea Manned Submersible

This study investigates the interaction and influence of surface cracks on the spherical pressure hull of a deep-sea manned submersible.The finite element model of the spherical hull is established,and a semi-elliptical surface crack is inserted in the welding toe of the spherical hull as the main crack.Considering the combined effect of external uniform pressure and welding residual stress at the weld toe,the stress intensity factor(SIF)is obtained based on the M-integral method.Inserting disturbing cracks at different positions on the spherical hull surface,the interaction and influence between multi-cracks are revealed by numerical calculation.The results show that the existence of the disturbing crack has a great influence on the stress intensity factor of the main crack,and the influence is different with the different location of disturbing crack.The study of the interaction of multiple cracks under different interference factors and the influence of disturbing cracks on the main crack can provide some reference for future engineering applications.

DYNAMIC FRACTURE MECHANICS ANALYSIS FOR COMPOSITE MATERIAL WITH MATERIAL NON-HOMOGENEITY IN THICKNESS DIRECTION

The problem considered here is the response of a non-homogeneouscomposite material containing some cracks subjected to dynamicloading. It is assumed that the composite material is or- thotropicand all material properties depend only on the coordinate y (alongthe thickness direcion). In the analysis, the elastic region isdivided into a number of plies of infinite length. The materialproperties are taken to be constants for each ply. By utilizing theLaplace transform and Fourier transform technique, the generalsolutions for plies are derived.

Virtual multi-dimensional internal bonds model and its application in simulation of rock mass

VMIB (virtual multi-dimensional internal bonds) is a multiscale mechanical model developed from the VIB (virtual internal bond) theory. In VIB theory,the solid mate-rial is considered to consist of random-distributed material particles in microscale. These particles are connected with normal bonds. The macro constitutive relation is derived from the cohesive law between particles. However,in VMIB,the micro particles are connected with both normal and shear bonds. The macro constitutive relation is derived in terms of bond stiffness coefficients. It has been theoretically certified that there exists a corresponding relationship between the two bond stiffness coefficients and the two macro material constants,i.e. the Young’s modulus and Poisson ratio. This corresponding relationship suggests that it should be necessary and sufficient to simultaneously account for the normal and shear interactions between particles. Due to the fact that the fracture criterion is directly incorporated into the constitutive relation,both VIB and VMIB present many advantages in simulating fractures of materials. In the damage model of rock mass,a damage tensor is usually defined to describe the distribution of cracks. The damage value in one direction determines the relative stiffness of rock mass in this direction. In VMIB solid,the relative distribution density of micro bonds in one direction determines the relative macro stiffness of the material in this direction. The effects of the damage value and the relative distribution density of bonds are consistent. To simulate the failure behavior of rock mass with VMIB,the presented paper sets up a quantitative relationship between the damage tensor and the rela-tive distribution density of bonds. Comparison of the theoretical and the experi-mental results shows that VMIB model can represent the effect of distributed cracks on rock mass with this relationship. The presented work provides a founda-tion for further simulating fracture behavior of rock mass with VMIB model,and an alternative approach for modeling other multi-cracked body.

Triangular element partition method with consideration of crack tip

In fracture simulation,how to model the pre-existing cracks and simulate their propagation without remeshing is an important topic.The newly developed triangular element partition method(TEPM)provides an efficient approach to this problem.It firstly meshes the cracked body regardless of the geometry integrity of the interesting object with triangular elements.After the meshing procedure is completed,some elements are intersected by cracks.For the element intersected by a crack,the TEPM takes the element partition technique to incorporate the discontinuity into the numerical model without any interpolation enrichment.By this approach,the TEPM can simulate fracture without mesh modification.In the TEPM,all the cracked elements are treated as the usual partitioned elements in which the crack runs through.The virtual node pairs(the intersection points of crack faces and elements)at the opposite faces of the crack move independently.Their displacements are respectively determined by their neighbor real nodes(nodes formatted in the original mesh scheme)at the same side of the crack.However,among these cracked elements,the element containing a crack tip,referred to as the crack tip element thereafter,behaves differently from those cut through by the crack.Its influence on the singular field at the vicinity of the fracture tip becomes increasingly significant with the element size increasing.In the crack tip element,the virtual node pair at the crack tip move consistently before fracture occurs while the virtual node pair separate and each virtual node moves independently after the fracture propagates.Accordingly,the crack tip element is automatically transformed into the usual partitioned element.In the present paper,the crack tip element is introduced into the TEPM to account for the effect of the crack tip.Validation examples indicate that the present method is almost free from the element size effect.It can reach the same precision as the conventional finite element method under the same meshing scheme.But the TEPM is much more efficient and convenient than the conventional finite element method because the TEPM avoids the troubles that the conventional finite element method suffers,e.g.,the meshing problem of cracked body,modification of mesh scheme,etc.Though the extended finite element method can also avoid these troubles,it introduces extra degrees of freedom due to node interpolation enrichment.Due to the simplicity of the present TEPM,it is believed that its perspective should be highly inspiring.

Non-Fourier heat conduction analysis of a 2-D plate with inner cracks at arbitrary direction angles

This paper investigates the problem of a 2-D plate containing internal cracks at arbitrary angles to the heating surface under thermal shock.The finite difference method is applied to solve the temperature distribution and a difference scheme of the oblique crack is developed.For the first time,a non-Fourier heat transfer analysis method of a 2-D plate with inner cracks at arbitrary direction angles is established.Specifically,for a strip with an inner crack paralleling to the surface,numerical results are compared with the analytical solutions by rotating the boundary,and perfect agreement is obtained.Numerical experiments are further presented to investigate the influence of crack orientation angle and multi-crack distribution on non-Fourier heat conduction.