Experimental study on cracking behaviour of moulded gypsum containing two non-parallel overlapping flaws under uniaxial compression

Failure of rock mass that is subjected to compres-sive loads occurs from initiation, propagation, and linkage of new cracks from preexisting fissures. Our research inves-tigates the cracking behaviour and coalescence process in a brittle material with two non-parallel overlapping flaws using a high-speed camera. The coalescence tensile crack and tensile wing cracks were the first cracks to occur from the pre-existing flaws. The initiation stresses of the primary cracks at the two tips of each flaw were simultaneous and decreased with reduced flaw inclination angle. The following types of coalescence cracks were identified between the flaws: pri-mary tensile coalescence crack, tensile crack linkage, shear crack linkage, mixed tensile-shear crack, and indirect crack coalescence. Coalescence through tensile linkage occurred mostly at pre-peak stress. In contrast, coalescence through shear or mixed tensile-shear cracks occurred at higher stress. Overall, this study indicates that the geometry of preexisting flaws affect crack initiation and coalescence behaviour.     

Experimental investigation on strength and failure behavior of pre-cracked marble under conventional triaxial compression

Fractures in natural rocks have an important effect on the strength and failure behavior of rock mass, which are often evaluated in rock engineering practice. The theoretical evaluation of mechanical behavior of fractured rock mass has no satisfactory answer due to the role of confining pressure and crack geometry. Therefore, in this paper, conventional triaxial compression experiments were carried out to study the strength and failure behavior of marble samples with two pre-existing closed cracks in non-overlapping geometry. Based on the experimental results of a number of triaxial compression tests, the effect of crack coalescence on the axial supporting capacity and deformation property were investigated with different confining pressures. The results show that intact samples and flawed samples (marble with pre-existing cracks) have different deformation properties after peak stress, which change from brittleness to plasticity and ductility with the increase of confining pressure. The peak strength and failure mode are found depending not only on the geometry of flaw, but also on the confining pressure. The strength of flawed samples shows distinct non-linear behavior, which is in a better agreement with non-linear Hoek–Brown criterion than linear Mohr–Coulomb criterion. For a kind of rock that has been evaluated as a Hoek–Brown material, a new evaluation criterion is put forward by adopting optimal approximation polynomial theory, which can be used to confirm more precisely the strength parameters (cohesion and internal friction angle) of flawed samples. For intact samples, the marble leads to typical shear failure mode with a single fracture surface under different confining pressures, while for flawed samples, under uniaxial compression and a lower confining pressure (σ₃ = 10 MPa), tests for coarse and medium marble (the coarse and medium refer to the grain size) exhibit three basic failure modes, i.e., tensile mode, shear mode, and mixed mode (tensile and shear). Shear mode is associated with lower strength behavior. However, under higher confining pressures (σ₃ = 30 MPa), for coarse marble, the axial supporting capacity is not related to the geometry of flaw. The friction among crystal grains determines the strength behavior of coarse marble. For medium marble, the failure mode and deformation behavior are dependent on the crack coalescence in the sample. The present research provides increased understanding of the fundamental nature of rock failure under conventional triaxial compression.     

An experimental and numerical study of fracture coalescence in pre-cracked specimens under uniaxial compression

This study presents crack initiation, propagation and coalescence at or near pre-existing open cracks or flaws in a specimen under uniaxial compression. The flaw geometry in the specimen was a combination of a horizontal flaw and an inclined flaw underneath. This flaw geometry is different from those reported in the previous studies, where a pair of parallel flaws was used. Three materials were used, PMMA (Poly Methyl MethAcrylate), Diastone (types of molded gypsum), and Hwangdeung granite. Crack initiation and propagation showed similar and different patterns depending on the material. In PMMA, tensile cracks initiated at the flaw tips and propagated to the tip of the other flaw in the bridge area. The cracks then coalesced at a point of the inclined flaw, which is affected by the flaw inclination angle. For Diastone and Hwangdeung granite, tensile cracks were observed followed by the initiation of shear cracks. Coalescence occurred mainly through the tensile cracks or tensile and shear cracks. Crack coalescence was classified according to the crack coalescence types of parallel flaws for overlapping flaw geometry in the past works. In addition, crack initiation and coalescence stresses in the double-flawed specimens were analyzed and compared with those in the single-flawed specimen. Numerical simulations using PFC^2D (Particle Flow Code in two dimensions) based on the DEM (Discrete Element Method) were carried out and showed a good agreement with the experimental results in the coalescence characteristics in Hwangdeung granite. These experimental and numerical results are expected to improve the understanding of the characteristics of cracking and crack coalescence and can be used to analyze the stability of rock and rock structures, such as the excavated underground openings or slopes, tunneling construction, where pre-existing cracks or fractures play a crucial role in the overall integrity of such structures.     

Influence of flaw inclination angle and loading condition on crack initiation and propagation

With reference to the experimental observation of crack initiation and propagation from pre-existing flaws in rock specimens under compression, the influences of pre-existing flaw inclination angle on the cracking processes were analyzed by means of finite element method (FEM) and non-linear dynamics method. FEM analysis on the stress field distribution induced by the presence of a pre-existing flaw provided better understanding for the influence of flaw inclination angle on the initiation position and initiation angle of the potential cracks. Numerical analysis based on the non-linear dynamics method was performed to simulate the cracking processes. The resultant crack types, crack initiation sequences and the overall crack pattern were different under different loading conditions. Under a relatively low loading rate or a small magnitude of maximum loading pressure, tensile cracks would tend to initiate prior to shear cracks. In contrast, under a relatively high loading rate and a large magnitude of maximum loading pressure, shear cracks would tend to initiate prior to tensile cracks instead.     

Numerical study of failure behaviour of pre-cracked rock specimens under conventional triaxial compression

Macroscopic pre-existing flaws play an important role in evaluating the strength and the failure modes of a heterogeneous rock mass. Crack initiation, propagation and coalescence from macroscopic pre-existing flaws are considered in a 3-D numerical model (RFPA3D) to investigate their effects on the underlying failure modes of rock. A feature of the code RFPA3D is that it can numerically simulate the evolution of cracks in three-dimensional space, as well as the heterogeneity of the rock mass. Three types of flaw geometries were evaluated numerically against experimental results: Type A for intact specimen, and Types B and C for flawed cylindrical specimens with different macroscopic pre-existing flaws, respectively. The effect of confining pressure on the fracture evolution was also considered. Numerical results showed that both the ligament angle and the flaw angle of two pre-existing cracks can affect the uniaxial compressive strength of the specimen and the mechanism of fracture evolution. In addition, both the uniaxial compressive strength and the accumulated acoustic emission increase with increasing heterogeneity.     

The effects of the rock bridge ligament angle and the confinement on crack coalescence in rock bridges: An experimental study and discrete element method

The present article investigates the influences of the rock bridge ligament angle, β, and the confinement on crack coalescence patterns by conducting laboratory and numerical tests on rock-like specimens. Laboratory tests show that no coalescence in the rock bridge occurred for low β. With an increase of β, tensile-shear coalescence and tensile coalescences subsequently occurred. In addition, the increase in the confinement first promoted shear coalescence and then restrained crack coalescence for low β, whereas the tensile coalescence was restrained by the increase in confinement for high β. The numerical results corroborate the laboratory tests in the coalescence patterns. In addition, the numerical study shows that tensile and shear cracks subsequently initiated near crack tips because of the concentrated tensile and shear stresses, respectively. Regarding the influence of β on crack coalescence, tensile or shear stress failed to concentrate in rock bridges for low β. Therefore, the cracks failed to coalesce, whereas with the increase in β, tensile and shear stress concentrations occurred in the bridge and led to either tensile shear or tensile coalescence. Regarding the influence of confinement on crack coalescence, the increase in confinement restrained the tensile stress concentrations and further hindered tensile crack coalescence in rock bridges for high values of β.     

Periodic group crack problems in an infinite plate

This paper investigates periodic group crack problems in an infinite plate. The periodic group crack is composed of infinite groups with numbering from j = −∞, …, −2, −1, 0, 1, 2, …, to j = ∞, and the groups are placed periodically. The same loading condition and the same geometry are assumed for cracks in all groups. A singular integral equation is used to solve the problems. The singular integral equation is formulated on cracks of the 0th group (or the central group) with the collection of influences from the infinite groups. The influences of many neighboring groups to the central group are evaluated exactly. Meantime, the influences of many remote groups to the central group can be summed up into one term approximately. The stress intensity factors at crack tips can be evaluated from the solution of the singular integral equation. It is found from some sample problems that the obtained results are very accurate. Finally, several numerical examples are presented and interaction among the group cracks is addressed.     

Squeeze flow of a power-law fluid between non-parallel plates

Squeeze flow in the gap between non-parallel circular plates of radius R is discussed. The test material is assumed to be a power-law fluid, with a no-slip boundary condition at the plates. If the mean separation between the plates is h, and the angle of inclination between the plates is ? ? h/R, the force on the plates is perturbed only at O(?²) and is increased by less than 10% if ? < 0.35h/R. A torque O(?) tends to return the plates to a parallel configuration.     

Experimental quantification of the plastic blunting process for stage II fatigue crack growth in one-phase metallic materials

The plastic blunting process during stage II fatigue crack growth was studied in pure polycrystalline Ni to investigate effects of strain localization and inelastic behavior on the kinematics of crack advance. Correlations were obtained between strain fields ahead of a fatigue crack, crack advance per cycle and crack growth kinetics. Strain fields were quantified using a combination of in situ loading experiments, scanning electron microscopy and digital image correlation for 8 < ΔK < 20 MPa m^1/2 and a fixed load ratio of 0.1. Results indicate that strain localized along a dominant deformation band, which was usually crystallographic and carried mostly pure shear for large loads and was of mixed character for lower loads. Instances of double deformation bands were observed, with bands acting either in a simultaneous or alternating fashion. It was found that the area integral of the opening strain for values larger than a given threshold, an “integrated” strain, had a power-law relationship with ΔK, with the exponent approximately equal to the Paris exponent (m). Therefore, the crack growth rate was proportional to the integrated strain. An analysis based on this correlation and the presence of dominant shear bands indicated that the integrated strain is related to the accumulated displacement in the band. This, in turn, is proportional to the product of the cyclic plastic zone radius and the average shear strain ahead of the tip, which represents a basic length scale for plastic blunting. Assumptions on the load dependence of these quantities, based on their observed spatial variation, allowed estimating $m=2\left(1+\frac{1}{1+n^{\prime }}\right)$, where n′ is the cyclic hardening exponent (0 < n′ < 1). This gives 3 < m < 4, which accounts for about 50% of the observed values of m between 1.5 and 6 for a wide variety of metallic materials.     

Experimental study of parallel and inclined turbulent wall jets

An experimental study of the flow field in a two-dimensional wall jet has been conducted. All measurements were carried out using hot-wire anemometry. The experimental facility has a rectangular slot nozzle of high aspect ratio l/b = 100 (where l and b are the length and height slot, respectively). Mean velocities and Reynolds stresses were determined with three nozzle Reynolds numbers (Re = 1 × 10⁴, 2 × 10⁴ and 3 × 10⁴) and four different inclination angles between the wall and the flow velocity at the nozzle (β = 0°, 10°, 20° and 30°). Results indicate that all wall jets are self-preserving in the developed region. Normal to the wall two regions can be identified: one similar to a plane free jet and the other similar to a boundary layer. Downstream the interaction between these two regions creates a mixed or third region. The logarithmic region increases with the distance from the nozzle and with the Reynolds number. For the inclined wall jet, the spreading rate expressed in terms of jet half-width or maximum velocity decay with respect to the streamwise distance, asymptotes to a linear law. The streamwise locations where the jet becomes self-similar are farther from the exit than in parallel wall jet. The slope of both half-width and maximum velocity decay in the developed region are affected by both wall jet inclination angle and nozzle exit Reynolds number.     

Experimental investigation of a developing two-phase bubbly flow in horizontal pipe

Experimental results for various water and air superficial velocities in developing adiabatic horizontal two-phase pipe flow are presented. Flow pattern maps derived from videos exhibit a new boundary line in intermittent regime. This transition from water dominant to water–gas coordinated regimes corresponds to a new transition criterion C_T = 2, derived from a generalized representation with the dimensionless coordinates of Taitel and Dukler.Velocity, turbulent kinetic energy and dissipation rate, void fraction and bubble size radial profiles measured at 40 pipe diameters for J_L = 4.42 m/s by hot film velocimetry and optical probes confirm this transition: the gas influence is not continuous but strongly increases beyond J_G = 0.06 m/s. The maximum dissipation rate, derived from spectra, is increased in two-phase flow by a factor 5 with respect to the single phase case.The axial evolution of the bubble intercept length histograms also reveal the flow organization in horizontal layers, driven by buoyancy effects. Bubble coalescence is attested by a maximum bubble intercept evolving from 2.5 to 4.5 mm along the pipe. Turbulence generated by the bubbles is also manifest by the 4-fold increase of the maximum turbulent dissipation rate along the pipe.     

Equilibrium of a pressurized plastic fluid in a wellbore crack

This paper investigates equilibrium of a pressurized plastic fluid invading a tensile wellbore crack in a linear elastic, permeable rock. The crack is initially filled by pore fluid at ambient pressure, that is immiscibly displaced by the plastic fluid invading from the wellbore. The plastic fluid comes to rest to form a “plug” within the elastically deformed crack when the limit equilibrium between the shear stresses generated at the fracture walls and the pressure drop between the wellbore wall and the crack tip is reached. The model assumes that the leak-off of the plastic fluid into the rock is negligible, while the displaced pore fluid in the crack tip region is freely exchanged with the surrounding permeable rock to maintain the ambient pressure level. When the crack length ? is small or large compared to the wellbore radius R, the problem reduces to that of a pressurized edge or Griffith’s crack, respectively, subjected to a uniform far-field confining stress. In these two end-member cases, the normalized solution for the net pressure distribution, the plug length, and the stress intensity factor at the crack tip is obtained as a function of two numbers – the normalized net fluid pressure at the crack inlet and at the crack tip (partial plugs only) – that embody the solution’s dependence on the wellbore and the far field loading, the fluid yield strength, and the rock modulus. In the general case of an intermediate crack length (? ～ R), the normalized solution is a function of two additional parameters, the length-to-radius ratio and a normalized measure of the far field stress anisotropy, respectively, which accurate approximation is devised from an end-member solution using a rescaling argument. The equilibrium plug solutions are used to evaluate the breakdown pressure, the critical wellbore pressure at which the crack propagation condition is first met, and to analyze the stability of the ensuing crack propagation.     

Some exact solutions for rotating flows of a generalized Burgers’ fluid in cylindrical domains

The velocity field and the adequate shear stress corresponding to the flow of a generalized Burgers’ fluid model, between two infinite co-axial cylinders, are determined by means of Laplace and finite Hankel transforms. The motion is due to the inner cylinder that applies a time dependent torsional shear to the fluid. The solutions that have been obtained, presented in series form in terms of usual Bessel functions J₁( ? ), J₂( ? ), Y₁( ? ) and Y₂( ? ), satisfy all imposed initial and boundary conditions. Moreover, the corresponding solutions for Burgers’, Oldroyd-B, Maxwell, second grade, Newtonian fluids and large-time transient solutions for generalized Burgers’ fluid are also obtained as special cases of the present general solutions. The effect of various parameters on large-time and transient solutions of generalized Burgers’ fluid is also discussed. Furthermore, for small values of the material parameters, λ₂ and λ₄ or λ₁, λ₂, λ₃ and λ₄, the general solutions corresponding to generalized Burgers’ fluids are going to those for Oldroyd-B and Newtonian fluids, respectively. Finally, the influence of the pertinent parameters on the fluid motion, as well as a comparison between models, is shown by graphical illustrations.     

The influence of pipe length on thermal statistics computed from DNS of turbulent heat transfer

We present results from direct numerical simulation of turbulent heat transfer in pipe flow at a bulk flow Reynolds number of 5000 and Prandtl numbers ranging from 0.025 to 2.0 in order to examine the effect of streamwise pipe length (πδ ≡ πD/2 ? L ? 12πδ) on the convergence of thermal turbulence statistics. Various lower and higher order thermal statistics such as mean temperature, rms of fluctuating temperature, turbulent heat fluxes, two-point auto and cross-correlations, skewness and flatness were computed and it is found that the value of L required for convergence of the statistics depends on the Prandtl number: larger Prandtl numbers requires comparatively shorter pipe length for convergence of most of the thermal statistics.     

Numerical analysis of turbulent convection heat transfer from an array of perforated fins

Numerical investigation is made for three-dimensional fluid flow and convective heat transfer from an array of solid and perforated fins that are mounted on a flat plate. Incompressible air as working fluid is modeled using Navier–Stokes equations and RNG based k ? ? turbulent model is used to predict turbulent flow parameters. Temperature field inside the fins is obtained by solving Fourier’s conduction equation. The conjugate differential equations for both solid and gas phase are solved simultaneously by finite volume procedure using SIMPLE algorithm. Perforations such as small channels with square cross section are arranged streamwise along the fin’s length and their numbers varied from 1 to 3. Flow and heat transfer characteristics are presented for Reynolds numbers from 2 × 10⁴ to 4 × 10⁴ based on the fin length and Prandtl number is taken Pr = 0.71. Numerical computations are validated with experimental studies of the previous investigators and good agreements were observed. Results show that fins with longitudinal pores, have remarkable heat transfer enhancement in addition to the considerable reduction in weight by comparison with solid fins.     

Unsteady gravity-driven slender rivulets of a power-law fluid

Unsteady gravity-driven flow of a thin slender rivulet of a non-Newtonian power-law fluid on a plane inclined at an angle α to the horizontal is considered. Unsteady similarity solutions are obtained for both converging sessile rivulets (when 0 < α < π/2) in the case x < 0 with t < 0, and diverging pendent rivulets (when π/2 < α < π) in the case x > 0 with t > 0, where x denotes a coordinate measured down the plane and t denotes time. Numerical and asymptotic methods are used to show that for each value of the power-law index N there are two physically realisable solutions, with cross-sectional profiles that are ‘single-humped’ and ‘double-humped’, respectively. Each solution predicts that at any time t the rivulet widens or narrows according to |x | ^{(2N+1)/2(N+1)} and thickens or thins according to |x | ^N/(N+1) as it flows down the plane; moreover, at any station x, it widens or narrows according to |t | ^?N/2(N+1) and thickens or thins according to |t | ^?N/(N+1). The length of a truncated rivulet of fixed volume is found to behave according to |t | ^N/(2N+1).     

Effect of gravity on the mechanical properties of lunar regolith tested using a low gravity simulation device

Simulations of the bearing capacity and shear strength of regolith under Earth’s gravity produce different results from those under low gravity. A low-gravity simulation device was developed in this study, and an internal stress model of regolith simulant was established to correct the errors. The model revealed additional force on both shear plane in the shear test and the press plate area in the pressure–sinkage test. The sinkage and shear test results showed that low gravity decreased the deformable index n, frictional modulus k_φ and cohesion c, whereas there were no obvious changes to the cohesive modulus k_c and internal friction angle φ. The sinkage generally increased as the gravity decreased under a consistent normal load larger than 50 N, but when the wheel load was lower than 50 N, the sinkage of the TYII-1 simulant was larger under 1 G than 1/6 G. Gravity had little effect on the shear strength of the regolith. However, the tractive thrust of the TYII-1 simulant was lower under 1/6 G than 1 G. The smaller difference was due to differences in the way the soils responded to changes in the gravity level for the TYII-2 simulant.     

An improved atomistic simulation based mixed-mode cohesive zone law considering non-planar crack growth

A novel and improved atomistic simulation based cohesive zone law characterizing interfacial debonding is developed which explicitly accounts for the non-planarity of the crack propagation. Group of atoms in the simulation constituting cohesive zones which are used to obtain local stress and crack opening displacement data are determined dynamically during the non-planar crack growth as they cannot be determined apriori. The methodology is used to study the debonding of Σ5 (2 1 0)/[0 0 1] symmetric tilt grain boundary interface in a Cu bicrystal under several mixed mode loading conditions. Simulations show that such bicrystalline specimen exhibits three types of energy dissipative mechanisms – shear coupled GB migration (SCM) away from the crack-tips, change in spacial orientation of GB structural units rendering highly disordered grain boundary near the crack tips and brittle intergranular fracture. Which combination of these three deformation mechanism will be active influencing the degree of non-planarity of the crack propagation at various stages of loading depends on the loading mode-mixity. As the ratio of shear component of the loading parallel to the GB plane and normal to the tilt axis with respect to the normal loading increases (thereby increasing the mode-mixity), overall strain-to-failure also increases and SCM tends to become the dominant deformation mechanism. Through this framework, analytical functional forms and parameters describing cohesive laws for both normal and shear traction as a function of the mode-mixity of the loading and crack opening displacement are predicted.