Cracked infinite cylinder with two rigid inclusions under axisymmetric tension

This paper considers the problem of an axisymmetric infinite cylinder with a ring shaped crack at z = 0 and two ring-shaped rigid inclusions with negligible thickness at z = ±L. The cylinder is under the action of uniformly distributed axial tension applied at infinity and its lateral surface is free of traction. It is assumed that the material of the cylinder is linearly elastic and isotropic. Crack surfaces are free and the constant displacements are continuous along the rigid inclusions while the stresses have jumps. Formulation of the mixed boundary value problem under consideration is reduced to three singular integral equations in terms of the derivative of the crack surface displacement and the stress jumps on the rigid inclusions. These equations, together with the single-valuedness condition for the displacements around the crack and the equilibrium equations along the inclusions, are converted to a system of linear algebraic equations, which is solved numerically. Stress intensity factors are calculated and presented in graphical form.     

A unified definition for stress intensity factors of interface corners and cracks

Based upon linear fracture mechanics, it is well known that the singular order of stresses near the crack tip in homogeneous materials is a constant value −1/2, which is nothing to do with the material properties. For the interface cracks between two dissimilar materials, the near tip stresses are oscillatory due to the order of singularity being −1/2 ± iε and −1/2. The oscillation index ε is a constant related to the elastic properties of both materials. While for the general interface corners, their singular orders depend on the corner angle as well as the elastic properties of the materials. Owing to the difference of the singular orders of homogeneous cracks, interface cracks and interface corners, their associated stress intensity factors are usually defined separately and even not compatibly. Since homogenous cracks and interface cracks are just special cases of interface corners, in order to build a direct connection among them a unified definition for their stress intensity factors is proposed in this paper. Based upon the analytical solutions obtained previously for the multibonded anisotropic wedges, the near tip solutions for the general interface corners have been divided into five different categories depending on whether the singular order is distinct or repeated, real or complex. To provide a stable and efficient computing approach for the general mixed-mode stress intensity factors, the path-independent H-integral based on reciprocal theorem of Betti and Rayleigh is established in this paper. The complementary solutions needed for calculation of H-integral are also provided in this paper. To illustrate our results, several different kinds of examples are shown such as cracks in homogenous isotropic or anisotropic materials, central or edge notches in isotropic materials, interface cracks and interface corners between two dissimilar materials.     

Constant moving crack in a magnetoelectroelastic material under anti-plane shear loading

Analytical solutions for an anti-plane Griffith moving crack inside an infinite magnetoelectroelastic medium under the conditions of permeable crack faces are formulated using integral transform method. The far-field anti-plane mechanical shear and in-plane electrical and magnetic loadings are applied to the magnetoelectroelastic material. Expressions for stresses, electric displacements and magnetic inductions in the vicinity of the crack tip are derived. Field intensity factors for magnetoelectroelastic material are obtained. The stresses, electric displacements and magnetic inductions at the crack tip show inverse square root singularities. The moving speed of the crack have influence on the dynamic electric displacement intensity factor (DEDIF) and the dynamic magnetic induction intensity factor (DMIIF), while the dynamic stress intensity factor (DSIF) does not depend on the velocity of the moving crack. When the crack is moving at very lower or very higher speeds, the crack will propagate along its original plane; while in the range of M_c1 < M < M_c2, the propagation of the crack possibly brings about the branch phenomena in magnetoelectroelastic media.     

Image decomposition method for the analysis of a mixed dislocation in a general multilayer

The elastic solutions for a mixed dislocation in a general multilayer with N dissimilar anisotropic layers are obtained via a generalized image decomposition method. The original problem is decomposed into N homogeneous subproblems with strategically placed continuously distributed image (virtual) dislocations which satisfy the consistency conditions for degenerate N − M (M < N) layer problems. The image dislocations are used to satisfy the interface or free surface conditions, and represent the unknowns of the problem. The resulting singular Cauchy integral equations are transformed into non-singular Fredholm integral equations of the second kind using certain H- and I-integral transforms. The Fredholm integral equations are then solved via the classical Nyström method. The general decomposition and the elimination of all singular integrals yield an exact formulation of the problem; the approximation arises only in the Nyström method. The dislocation mixity and the number of layers dissimilar in thickness and elastic anisotropy can be handled without difficulty, constrained only by the number of linear algebraic equations in the Nyström method for large N. For the numerical study, image forces on a dislocation in two- and three-layer systems are calculated. The accuracy of the results is verified by checking the boundary conditions and by comparison with previous results. The dependence of the image force on the dislocation position and mixity, and on the layer thicknesses and elastic anisotropies, is also illustrated via numerical investigations.     

Numerical study on coalescence of two pre-existing coplanar flaws in rock

Crack propagation and coalescence processes are the fundamental mechanisms leading to progressive failure processes in rock masses, in which parallel non-persistent rock joints are commonly involved. The coalescence behavior of the latter, which are represented as pre-existing coplanar flaws (cracks), is numerically investigated in the present study. By using AUTODYN as the numerical tool, the present study systematically simulates the coalescence of two pre-existing coplanar flaws in rock under compression. The cumulative damage failure criterion is adopted in the numerical models to simulate the cumulative damage process in the crack initiation and propagation. The crack types (shear or tensile) are identified by analyzing the mechanics information associated with the crack initiation and propagation processes. The simulation results, which are generally in a good accordance with physical experimental results, indicate that the ligament length and the flaw inclination angle have a great influence on the coalescence pattern. The coalescence pattern is relatively simple for the flaw arrangements with a short ligament length, which becomes more complicated for those with a long ligament length. The coalescence trajectory is composed of shear cracks only when the flaw inclination angle is small (such as β ? 30°). When the pre-existing flaws are steep (such as β ? 75°), the coalescence trajectory is composed of tensile cracks as well as shear cracks. When the inclination angle is close to the failure angle of the corresponding intact rock material, and the ligament length is not long (such as L ? 2a), the direct shear coalescence is the more favorable coalescence pattern. In the special case that the two pre-existing flaws are vertical, the model will have a direct tensile coalescence pattern when the ligament length is short (L ? a), while the coalescence between the two inner flaw tips is not easy to achieve if the ligament length is long (L ? 2a).     

The degenerate scale problem for the Laplace equation and plane elasticity in a multiply connected region with an outer circular boundary

This paper investigates the degenerate scale problem for the Laplace equation and plane elasticity in a multiply connected region with an outer circular boundary. Inside the boundary, there are many voids with arbitrary configurations. The problem is analyzed with a relevant homogenous BIE (boundary integral equation). It is assumed that all the inner void boundary tractions are equal to zero, and tractions on the outer circular boundary are constant. Therefore, all the integrations in BIE are performed on the outer circular boundary only. By using the relation z * conjg(z) = a * a, or conjg(z) = a * a/z on the circular boundary with radius a, all integrals can be reduced to an integral for complex variable and they can be integrated in closed form. The degenerate scale a = 1 is found in the Laplace equation and in plane elasticity regardless of the void configuration.     

Irreversible thermodynamics of triple junctions during the intergranular void motion under the electromigration forces

A rigorous reformulation of internal entropy production and the rate of entropy flow is developed for multi-component systems consisting of heterophases, interfaces and/or surfaces. The result is a well-posed moving boundary value problem describing the dynamics of curved interfaces and surfaces associated with voids and/or cracks that are intersected by grain boundaries. Extensive computer simulations are performed for void configuration evolution during intergranular motion. In particular we simulate evolution resulting from the action of capillary and electromigration forces in thin film metallic interconnects having a “bamboo” structure, characterized by grain boundaries aligned perpendicular to the free surface of the metallic film interconnects. Analysis of experimental data utilizing previously derived mean time to failure formulas gives consistent values for interface diffusion coefficients and enthalpies of voids. 3.0 × 10⁻⁶ exp(−0.62 eV/kT) m² s⁻¹ is the value obtained for voids that form in the interior of the aluminum interconnects without surface contamination. 6.5 × 10⁻⁶ exp(−0.84 eV/kT) m² s⁻¹ is obtained for those voids that nucleate either at triple junctions or at the grain boundary-technical surface intersections, where the chemical impurities may act as trap centers for hopping vacancies.     

Extending the bounds of ‘steady’ RANS closures: Toward an instability-sensitive Reynolds stress model

The incapability of the conventional Unsteady RANS (Reynolds–Averaged Navier Stokes) models to adequately capture turbulence unsteadiness presents the prime motivation of the present work, which focuses on formulating an instability-sensitive, eddy-resolving turbulence model on the Second-Moment Closure level. The model scheme adopted, functioning as a ‘sub-scale’ model in the Unsteady RANS framework, represents a differential near-wall Reynolds stress model formulated in conjunction with the scale-supplying equation governing the homogeneous part of the inverse turbulent time scale ω_h (ω_h = ɛ_h/k). The latter equation was straightforwardly obtained from the model equation describing the dynamics of the homogeneous part of the total viscous dissipation rate ɛ, defined as ɛ_h = ɛ − 0.5ν∂²k/(∂x_j∂x_j) (Jakirlic and Hanjalic, 2002), by applying the derivation rules to the expression for ω_h. The model capability to account for vortex length and time scales variability was enabled through an additional term in the corresponding length-scale determining equation, providing a selective enhancement of its production, pertinent particularly to the highly unsteady separated shear layer region, modeled in terms of the von Karman length scale (comprising the second derivative of the velocity field) in line with the SAS (Scale-Adaptive Simulation) proposal (Menter and Egorov, 2010). The present model formulation, termed as SRANS model (Sensitized RANS), does not comprise any parameter depending explicitly on grid spacing. The predictive capabilities of the newly proposed length-scale determining model equation, solved in conjunction with Jakirlic and Hanjalic’s (2002) Reynolds stress model equation, are presently demonstrated by computing the flow configurations of increasing complexity featured by boundary layer separation from sharp-edged and continuous curved surfaces: backward-facing step flow, flow over a wall-mounted fence, flow over smoothly contoured periodically arranged hills and flow in a 3-D diffuser. The model performances are also assessed in capturing the natural decay of the homogeneous isotropic turbulence and the near-wall Reynolds stress anisotropy in a plane channel. In most cases considered the fluctuating velocity field was obtained starting from steady RANS results.     

Computational and experimental characterization of a liquid jet plunging into a quiescent pool at shallow inclination

A circular water jet (Re = 1.6 × 10⁵; We = 8.8 × 10³) plunging at shallow angles (θ ≈ 12.5°) into a quiescent pool is investigated computationally and experimentally. A surprising finding from the work is that cavities, of the order of jet diameter, are formed periodically in the impact location, even though the impinging flow is smooth and completely devoid of such a periodicity. Computational prediction of these frequencies was compared with experimental findings, yielding excellent agreement. The region in the vicinity of the impact is characterized by strong churning due to splashing and formation of air cavities. Measured velocity profiles indicate a concentration of momentum beneath the free surface slightly beyond the impact location (X/D_j ≈ 14), with a subsequent shift towards the free surface further downstream of this point (X/D_j ≈ 30). This shift is due primarily to the action of buoyancy on the cavity/bubble population. Comparisons of the mean velocity profile between simulations and experiments are performed, yielding good agreement, with the exception of the relatively small churning flow region. Further downstream (X/D_j ≳ 40), the flow develops mostly due to diffusion and the location of peak velocity coincides with the free surface. In this region, the free surface acts as an adiabatic boundary and restricts momentum diffusion, causing the peak velocity to occur at the free surface.     

Integral-based constitutive equation for rubber at high strain rates

High-speed experiments were conducted to characterize the deformation and failure of Styrene Butadiene Rubber at impact rates. Dynamic tensile stress–strain curves of uniaxial strip specimens and force–extension curves of thin sheets were obtained from a Charpy tensile impact apparatus. Results from the uniaxial tension tests indicated that although the rubber became stiffer with increasing strain rates, the stress–strain curves remained virtually the same above 280 s⁻¹. Above this critical strain rate, strength, fracture strain and toughness decreased with increasing strain rates. When strain rates were below 180 s⁻¹, the initial modulus, tensile strength and breaking extension increased as the strain rate increased. Between strain rates of 180 and 280 s⁻¹, the initial modulus and tensile strength increased with increasing strain rates but the extension at break decreased with increasing strain rates. A hyper-viscoelastic constitutive relation of integral form was used to describe the rate-dependent material behavior of the rubber. Two characteristic relaxation times, 5 ms and 0.25 ms, were needed to fit the proposed constitutive equation to the data. The proposed constitutive equation was implemented in ABAQUS Explicit via a user-defined subroutine and used to predict the dynamic response of the rubber sheets in the experiments. Numerical predictions for the transient deformation and failure of the rubber sheet were within 10% of experimental results.     

Thermo-elastic crack branching in general anisotropic media

Thermal fields may exist in addition to mechanical loading, for example, due to short term exposure to fire. In this paper, the branching of cracks in the presence of combined thermal and mechanical loads is investigated for general anisotropic media by employing the theory of Stroh’s dislocation formalism, extended to thermo-elasticity in matrix notation. A general solution to the thermo-elastic crack problem for an anisotropic material under arbitrary loading is obtained in a compact form. Green’s functions are also presented for a thermal dislocation (heat vortex) and a conventional dislocation (or, referred as mechanical dislocation), which are formulated considering the cuts located at an arbitrary angle with respect to the x₁ axis of the coordinate system (x₁, x₂, x₃). Using the derived compact expressions, the interaction between the crack and the dislocation is studied and a closed form solution for this interaction is obtained. The branching portion of the thermo-elastic crack is modelled as a continuous distribution of dislocations. This problem is then converted into a set of singular integral equations. Numerical results are presented to illustrate the possible effects of thermal loading on the propagation of the branched crack.     

Coplanar propagation paths of 3D cracks in infinite bodies loaded in shear

Bower and Ortiz, recently followed by Lazarus, developed a powerful method, based on a theoretical work of Rice, for numerical simulation of planar propagation paths of mode 1 cracks in infinite isotropic elastic bodies. The efficiency of this method arose from the need for the sole 1D meshing of the crack front. This paper presents an extension of Rice’s theoretical work and the associated numerical scheme to mixed-mode (2 + 3) shear loadings. Propagation is supposed to be channeled along some weak planar layer and to remain therefore coplanar, as in the case of a geological fault for instance. The capabilities of the method are illustrated by computing the propagation paths of cracks with various initial contours (circular, elliptic, rectangular, heart-shaped) in both fatigue and brittle fracture. The crack quickly reaches a stable, almost elliptic shape in all cases. An approximate but accurate analytic formula for the ratio of the axes of this stable shape is derived.     

Experimental investigation of horizontal forced-vibration effect on air-water two-phase flow

In order to investigate the potential seismic vibrations effect on two-phase flow in an annular channel, experimental tests with air-water two-phase flow under horizontal vibrations were carried out. A low-speed eccentric-cam vibration module capable of operating at motor speed of 45–1200 rpm (f = 0.75–20 Hz) was attached to an annular channel, which was scaled down from a prototypic BWR fuel sub-channel with inner and outer diameters of 19.1 mm and 38.1 mm, respectively. The two-phase flow was operated in the ranges of 〈j_f〉 = 0.25–1.00 m/s and 〈j_g〉 = 0.03–1.46 m/s with 27 flow conditions, and the vibration amplitudes controlled by cam eccentricity (E) were designed for the range of 0.8–22.2 mm. Ring-type impedance void meters were utilized to detect the area-averaged time-averaged void fraction under stationary and vibration conditions. A systematic experimental database was built and analyzed with effective maps in terms of flow conditions (〈j_g〉-〈j_f〉) and vibration conditions (E-f and f-a), and the potential effects were expressed by regions on the maps. In the 〈j_g〉-〈j_f〉 maps, the void fraction was found to potentially decrease under vibrations in bubbly flow regime and relatively lower liquid flow conditions, which may be explained by the increase of distribution parameter. Whereas and the void fraction may increase at the region closed to bubbly-to-slug transition boundary under vibrations, which may be explained by the changes of drift velocity due to flow regime change from bubbly to slug flows. No significant change in void fraction was found in slug flow regime under the present test conditions.     

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.     

Emission of organic carbon,elemental carbon and water-soluble ions from crop straw burning under flaming and smoldering conditions

Emissions from major agricultural residues were measured using a self-designed combustion system. Emission factors (EFs) of organic carbon (OC), elemental carbon (EC), and water-soluble ions (WSIs) (K⁺, NH₄⁺, Na⁺, Mg²⁺, Ca²⁺, Cl⁻, NO₃⁻, SO₄^2–) in smoke from wheat and rice straw were measured under flaming and smoldering conditions. The OC₁/TC (total carbon) was highest (45.8% flaming, 57.7% smoldering) among carbon fractions. The mean EFs for OC (EF_OC) and EC (EF_EC) were 9.2 ± 3.9 and 2.2 ± 0.7 g/kg for wheat straw and 6.4 ± 1.9 and 1.1 ± 0.3 g/kg for rice straw under flaming conditions, while they were 40.8 ± 5.6 and 5.8 ± 1.0 g/kg and 37.6 ± 6.3 and 5.0 ± 1.4 g/kg under smoldering conditions, respectively. Higher EC ratios were observed in particulate matter (PM) mass under flaming conditions. The OC and EC for the two combustion patterns were significantly correlated (p < 0.01, R = 0.95 for wheat straw; p < 0.01, R = 0.97 for rice straw), and a higher positive correlation between OC₃ and EC was observed under both combustion conditions. WSIs emitted from flaming smoke were dominated by Cl⁻ and K⁺, which contributed 3.4% and 2.4% of the PM mass for rice straw and 2.2% and 1.0% for wheat straw, respectively. The EFs of Cl⁻ and K⁺ were 0.73 ± 0.16 and 0.51 ± 0.14 g/kg for wheat straw and 0.25 ± 0.15 and 0.12 ± 0.05 g/kg for rice straw under flaming conditions, while they were 0.42 ± 0.28 and 0.12 ± 0.06 g/kg and 0.30 ± 0.27 and 0.05 ± 0.03 g/kg under smoldering conditions, respectively. Na⁺, Mg²⁺, and NH₄⁺ were vital components in PM, comprising from 0.8% (smoldering) to 3.1% (flaming) of the mass. Strong correlations of Cl⁻ with K⁺, NH₄⁺, and Na⁺ ions were observed in rice straw and the calculated diagnostic ratios of OC/EC, K⁺/Na⁺ and Cl⁻/Na⁺ could be useful to distinguishing crop straw burning from other sources of atmospheric pollution.     

Multiple-zone sliding contact with friction on an anisotropic thermoelastic half-space

The paper studies a class of multiple-zone sliding contact problems. This class is general enough to include frictional and thermal effects, and anisotropic response of the indented material. In particular, a rigid die (indenter) slides with Coulomb friction and at constant speed over the surface of a deformable and conducting body in the form of a 2D half-space. The body is assumed to behave as a thermoelastic transversely isotropic material. Thermoelasticity of the Green–Lindsay type is assumed to govern. The solution method is based on integral transforms and singular integral equations. First, an exact transform solution for the auxiliary problem of multiple-zone (integer n > 1) surface tractions is obtained. Then, an asymptotic form for this auxiliary problem is extracted. This form can be inverted analytically, and the result applied to sliding contacts with multiple zones. For illustration, detailed calculations are provided for the case of two (n = 2) contact zones. The solution yields the contact zone width and location in terms of sliding speed, friction, die profile, and also the force exerted. Calculations for the hexagonal material zinc illustrate effects of speed, friction and line of action of the die force on relative contact zone size, location of maximal values for the temperature and the compressive stress, and the maximum temperature for a given maximum stress. Finally, from our general results, a single contact zone solution follows as a simple limit.     

Superiority of PANS compared to LES in predicting a rudimentary landing gear flow with affordable meshes

Results of simulations of the flow around a rudimentary landing gear are presented in the paper. A newly proposed improved Partially-Averaged Navier–Stokes (PANS) method using k − ε − ζ − f turbulence model is used for prediction of the flow. The results are compared with the experimental data but also with the results of two LES simulations performed using the PANS computational grids. PANS simulations predicted the flow in good agreement with the experimental data. LES predicted a non-physical creation of separation over the front wheels that does not exist in the PANS prediction and was not observed in the experimental oil film. PANS simulations showed low sensitivity to the grid refinement. They show clear advantage compared with the LES simulations when the computational grid is inadequate for resolution of the near-wall flow structures.     

Spatiotemporal variations of PM2.5 and PM10 concentrations between 31 Chinese cities and their relationships with SO2, NO2, CO and O3

The variations of mass concentrations of PM_2.5, PM₁₀, SO₂, NO₂, CO, and O₃ in 31 Chinese provincial capital cities were analyzed based on data from 286 monitoring sites obtained between March 22, 2013 and March 31, 2014. By comparing the pollutant concentrations over this length of time, the characteristics of the monthly variations of mass concentrations of air pollutants were determined. We used the Pearson correlation coefficient to establish the relationship between PM_2.5, PM₁₀, and the gas pollutants. The results revealed significant differences in the concentration levels of air pollutants and in the variations between the different cities. The Pearson correlation coefficients between PMs and NO₂ and SO₂ were either high or moderate (PM_2.5 with NO₂: r = 0.256–0.688, mean r = 0.498; PM₁₀ with NO₂: r = 0.169–0.713, mean r = 0.493; PM_2.5 with SO₂: r = 0.232–0.693, mean r = 0.449; PM₁₀ with SO₂: r = 0.131–0.669, mean r = 0.403). The correlation between PMs and CO was diverse (PM_2.5: r = 0.156–0.721, mean r = 0.437; PM₁₀: r = 0.06–0.67, mean r = 0.380). The correlation between PMs and O₃ was either weak or uncorrelated (PM_2.5: r = −0.35 to 0.089, mean r = −0.164; PM₁₀: r = −0.279 to 0.078, mean r = −0.127), except in Haikou (PM_2.5: r = 0.500; PM₁₀: r = 0.509).