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.     

Ductile–brittle transitions in the fracture of plastically-deforming,adhesively-bonded structures. Part I: Experimental studies

Rate effects for adhesively-bonded joints in steel sheets failing by mode-I fracture and plastic deformation were examined. Three types of test geometries were used to provide a range of crack velocities between 0.1 and 5000 mm/s: a DCB geometry under displacement control, a wedge geometry under displacement control, and a wedge geometry loaded under impact conditions. Two fracture modes were observed: quasi-static crack growth and dynamic crack growth. The quasi-static crack growth was associated with a toughened mode of failure; the dynamic crack growth was associated with a more brittle mode of failure. The experiments indicated that the fracture parameters for the quasi-static crack growth were rate independent, and that quasi-static crack growth could occur even at the highest crack velocities. Effects of rate appeared to be limited to the ease with which a transition to dynamic fracture could be triggered. This transition appeared to be stochastic in nature, it did not appear to be associated with the attainment of any critical value for crack velocity or loading rate. While the mode-I quasi-static fracture behavior appeared to be rate independent, an increase in the tendency for dynamic fracture to be triggered as the crack velocity increased did have the effect of decreasing the average energy dissipated during fracture at higher loading rates.     

Influence of temperature on a carbon–fibre epoxy composite subjected to static and fatigue loading under mode-I delamination

In this paper, interlaminar crack initiation and propagation under mode-I with static and fatigue loading of a composite material are experimentally assessed for different test temperatures. The material under study is made of a 3501-6 epoxy matrix reinforced with AS4 unidirectional carbon fibres, with a symmetric laminate configuration [0°]_16/S. In the experimental programme, DCB specimens were tested under static and fatigue loading. Based on the results obtained from static tests, fatigue tests were programmed to analyse the mode-I fatigue behaviour, so the necessary number of cycles was calculated for initiation and propagation of the crack at the different temperatures. G–N curves were determined under fatigue loading, N being the number of cycles at which delamination begins for a given energy release rate. G_ICmax–a, a–N and da/dN–a curves were also determined for different G_cr rates (90%, 85%, 75%, etc.) and different test temperatures: 90 °C, 50 °C, 20 °C, 0 °C, ?30 °C and ?60 °C.     

Energy release rates in rubber during dynamic crack propagation

The theoretical understanding of the fracture mechanics of rubber is not as well developed as for other engineering materials, such as metals. The present study is intended to further the understanding of the dissipative processes that take place in rubber in the vicinity of a propagating crack tip. This dissipation contributes significantly to the total fracture toughness of the rubber and is therefore of great interest from a fracture mechanics point of view. To study this, a computational framework for analysing high-speed crack growth in a biaxially stretched rubber under plane stress is therefore formulated. The main purpose is to investigate the energy release rates required for crack propagation under different modes of biaxial stretching. The results show, that inertia comes into play when the crack speed exceeds about 50 m/s. The total work of fracture by far exceeds the surface energy consumed at the very crack tip, and the difference must be attributed to dissipative damage processes in the vicinity of the crack tip. The size of this damage/dissipation zone is expected to be a few millimetres.     

Lattice orientation effects on void growth and coalescence in fcc single crystals

Void growth and coalescence in fcc single crystals were studied using crystal plasticity under uniaxial and biaxial loading conditions and various orientations of the crystalline lattice. A 2D plane strain unit cell with one and two cylindrical voids was employed using three-dimensional 12 potentially active slip systems. The results were compared to five representative orientations of the tensile axis on the stereographic triangle. For uniaxial tension conditions, the void volume fraction increase under the applied load is strongly dependent on the crystallographic orientation with respect to the tensile axis. For some orientations of the tensile axis, such as [1 0 0] or [1 1 0], the voids exhibited a growth rate twice as fast compared with other orientations ([1 0 0], [2 1 1]). Void growth and coalescence simulations under uniaxial loading indicated that during deformation along some orientations with asymmetry of the slip systems, the voids experienced rotation and shape distortion, due mainly to lattice reorientation. Coalescence effects are shown to diminish the influence of lattice orientation on the void volume fraction increase, but noteworthy differences are still present. Under biaxial loading conditions, practically all differences in the void volume fraction for different orientations of the tensile axes during void growth vanish. These results lead to the conclusion that at microstructural length scales in regions under intense biaxiality/triaxiality conditions, such as crack tip or notched regions, the plastic anisotropy due to the initial lattice orientation has only a minor role in influencing the void growth rate. In such situations, void growth and coalescence are mainly determined by the stress triaxiality, the magnitude of accumulated strain, and the spatial localization of such plastic strains.     

Experimental study of the blockage boundary for dense-phase pneumatic conveying of powders through a horizontal slit

The estimation of the blockage boundary for pneumatic conveying through a slit is of significant importance. In this paper, we investigate the characteristics for blockage of powder (48 μm average diameter) through a horizontal slit (1.6 m × 0.05 m × 0.002 m). The results show that the required critical solid mass flow rate increases as the superficial air velocity increases superficial air velocity. The solid loading ratio and superficial air velocity displayed a decreasing power law relationship. This finding agrees with existing theory and experimental results. However, a minimum inlet solid loading ratio exists. When the air velocity is greater than the corresponding air velocity of the minimum solid loading ratio, the solid loading ratio exhibits an increasing trend in power law. We also found that when the inlet conveying pressure increased, the critical solid mass flow rate required for blockage, the inlet solid loading ratio, and the minimum inlet solid loading ratio increased.     

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.     

Atomization of glycerin with a twin-fluid swirl nozzle

Atomization of liquids with high viscosity is always a challenge, especially when small diameter droplets and high liquid flow rates are simultaneously required. In the present research, the performance of a Venturi–vortex twin-fluid swirl nozzle is examined, attending to its capabilities to generate droplets with diameters below 20 µm when atomizing pure glycerin at room temperature. In this nozzle, air is injected tangentially in a central convergent section, and discharges suctioning the liquid fed to a coaxial chamber, here using a gear pump. The resulting spray is visualized and analyzed. Droplet size distributions are measured with a laser diffractometer. As expected, droplet diameter increases with liquid flow rate, and quickly diminishes when air flow rate is increased. Sauter mean diameters (SMD) below 15 µm can be obtained even when atomizing pure glycerin. However, these values are obtained for relatively low glycerin flow rates (∼5 l/h), and with rather wide distributions. For 10 l/h and an air-to-liquid mass flow rate ratio (ALR) of 13.7 more than 26% of the glycerin volume is atomized in droplets smaller than 20 µm. Liquid ligaments are observed near the nozzle exit, but they tend to break up while moving downstream.     

Dynamic fracture of concrete – compact tension specimen

The behavior of concrete structures is strongly influenced by the loading rate. Compared to quasi-static loading concrete loaded by impact loading acts in a different way. First, there is a strain-rate influence on strength, stiffness, and ductility, and, second, there are inertia forces activated. Both influences are clearly demonstrated in experiments. Moreover, for concrete structures, which exhibit damage and fracture phenomena, the failure mode and cracking pattern depend on loading rate. In general, there is a tendency that with the increase of loading rate the failure mode changes from mode-I to mixed mode. Furthermore, theoretical and experimental investigations indicate that after the crack reaches critical speed of propagation there is crack branching. The present paper focuses on 3D finite-element study of the crack propagation of the concrete compact tension specimen. The rate sensitive microplane model is used as a constitutive law for concrete. The strain-rate influence is captured by the activation energy theory. Inertia forces are implicitly accounted for through dynamic finite element analysis. The results of the study show that the fracture of the specimen strongly depends on the loading rate. For relatively low loading rates there is a single crack due to the mode-I fracture. However, with the increase of loading rate crack branching is observed. Up to certain threshold (critical) loading rate the maximal crack velocity increases with increase of loading rate, however, for higher loading rates maximal velocity of the crack propagation becomes independent of the loading rate. The critical crack velocity at the onset of crack branching is found to be approximately 500 m/s.     

On the blast resistance of laminated glass

Blast resistant glazing systems typically use laminated glass to reduce the risk of flying glass debris in the event of an explosion. Laminated glass has one or more bonded polymer interlayers to retain glass fragments upon fracture. With good design, the flexibility of the interlayer and the adhesion between layers enable laminated glass to continue to resist blast after the glass layers fracture. This gives protection from significantly higher blast loads when compared to a monolithic pane. Full-scale open-air blast tests were performed on laminated glass containing a polyvinyl butyral (PVB) interlayer. Test windows of size 1.5 m × 1.2 m were secured to robust frames using structural silicone sealant. Blast loads were produced using charge masses of 15 kg and 30 kg (TNT equivalent) at distances of 10–16 m. Deflection and shape measurements of deforming laminated glass were obtained using high-speed digital image correlation. Measurements of loading at the joint, between the laminated glass and the frame, were obtained using strain gauges. The main failure mechanisms observed were the cohesive failure of the bonded silicone joint and delamination between the glass and interlayer at the pane edge. A new finite element model of laminated glass is developed and calibrated using laboratory based tests. Predictions from this model are compared against the experimental results.     

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.     

Experimental pool boiling investigations of FC-72 on silicon with artificial cavities and integrated temperature microsensors

In this experimental study, fluorinert FC-72 is boiled on a silicon chip with artificial cavities and integrated microsensors. The horizontal silicon chip with dimensions of 39.5 × 19 × 0.38 mm is completely immersed in FC-72. The integrated nickel–titanium temperature microsensors on the back of the chip are calibrated individually and exhibit a near-linear increase of electrical resistance with temperature. The applied heat fluxes and the resulting wall superheat at the boiling surface are varied by means of an integral thin-film resistance heater (95% Al, 4% Cu and 1% Si), also on the back of the silicon chip. Artificial cylindrical cavities with a mouth diameter of 10 μm and depths of 40, 80 or 100 μm situated above the microthermometers serve as artificial nucleation sites, due to trapped vapour. Bubble growth rates, frequencies, departure diameters of bubbles and waiting times between bubbles from an isolated cavity for different wall superheats and pressures were obtained by analysing high-speed video images and the simultaneously measured temperature below the artificial cavity.     

Atomization of liquids by disintegrating thin liquid films using gas jets

Liquid atomization is useful in many applications, such as engineering, science, pharmaceutics, medicine, forensics and others. In the present research, an innovative methodology and a new device for atomization of liquids into mists of micron and submicron droplets have been developed. The new liquid-atomization method exploits the physical phenomenon of fragmentation of thin liquid films into fine micron and submicron droplets by gas jets. For several tested prototypes, the direct observations using a high-speed visualization technique have demonstrated that bubbles were generated within a liquid and their shells have been subsequently destroyed by applying a mechanical impulse (pressure of a compressed air) once the bubbles came over the liquid surface. The main characteristics of the generated tap water mists have been experimentally measured by means of the laser diffraction technique under various conditions for each prototype. One of the prototype devices allowed obtaining mists containing 90–99% of droplets smaller than 1 µm, with the minimum arithmetic and Sauter mean droplet diameters of 1.48 µm and 2.66 µm, and the 2.64 ml/min of droplet flow rate for 3.5 bar manometer pressure of atomizing air. The gas to liquid mass ratios (GLR) in the new device are depending on the atomizing tube length and the number of perforated orifices in the tube: more the tube length, hence more the number of perforated orifices, and therefore more liquid droplets will form for the same gas flow rate. The measured GLR values related to 1 m length of the utilized atomizing tube were in the range of 0.65–1.06, and for the specifically utilized atomizing tube of 72 mm length were among 9.07–14.67. The results of this study demonstrate that the developed method of generation of very fine droplet mists has many advantages over the existing techniques and can be perspective for many practical applications.     

Capture of phosphates in surface water by TiO2 nanoparticles under UV irradiation

The capture of orthophosphates and total phosphorus from the Pudong Canal river in the Pudong District of Shanghai by TiO₂ nanoparticles is studied using a rotating photoreactor and the nano-TiO₂ photocatalyst Degussa P25. The effects of UV irradiation intensity in a range of 20–74 mW/cm², the loading of the TiO₂ nanoparticles in a range of 0.05–0.1 g/L, irradiation time up to 4 h, and pH values in a range of 2–10.5 on the capture efficiency are investigated. The results show that the capture of orthophosphates and total P are significantly enhanced by UV irradiation; at a loading of 0.1 g/L and an irradiation intensity above 36 mW/cm², orthophosphates and total phosphorus are rapidly captured by TiO₂ nanoparticles, causing an observed reduction from 0.4 mg/L down to 0.02 mg/L. pH values in a range of 2–10.5 have little effect on the capture efficiency of orthophosphates and total phosphorus.     

Thermo-mechanical behaviour of TRIP 1000 steel sheets subjected to low velocity perforation by conical projectiles at different temperatures

This paper presents and analyzes the behaviour of TRIP 1000 steel sheets subjected to low velocity perforation by conical projectiles. The relevance of this material resides in the potential transformation of retained austenite to martensite during impact loading. This process leads to an increase in strength and ductility of the material. However, this transformation takes place only under certain loading conditions strongly dependent on the initial temperature and deformation rate. In order to study the material behaviour under impact loading, perforation tests have been performed using a drop weight tower. Experiments were carried out at two different initial temperatures T₀ = 213 K and T₀ = 288 K, and within the range of impact velocities 2.5 m/s ? V₀ ? 4.5 m/s. The experimental setup enabled the measuring of impact velocity, residual velocity, load-time history and failure mode. In addition, dry and lubricated contacts between the striker and the plate have been investigated. Finally, by using X-ray diffraction it has been shown that no martensitic transformation takes place during the perforation process. The causes involving the none-appearance of martensite are examined.     

Innovative iteration technique for nonlinear ordinary differential equations of large deflection problem of circular plates

The large deflection problem of a bending circular plate is reduced to a nonlinear ordinary differential equation, and a type of pseudo-linearization is used to obtain the final solution. An innovative iteration technique is suggested. The initial iteration values of some functions for a given loading are adopted from the previous solution for a lower loading. This will significantly extend the range of solution for the non-dimensional loading. Previously, the solution for the non-dimensional loading Q = 100 was obtained. However, in this paper the solution for the non-dimensional loading Q = 1000 is achieved.     

A thermomechanical constitutive model for phase transformations in silicon under pressure and contact loading conditions

Most of the technologically relevant abrasive machining techniques for silicon (Si) such as lapping, sawing and grinding are based on the interaction of the silicon surface with a hard particle or asperity. It has been long established that the governing deformation mechanism for Si under such contact loading conditions is stress induced phase transformation. The present work introduces a novel phenomenological constitutive model for phase transformations of silicon set up in a thermomechanical framework of broad applicability. Taking into account experimental observations as well as first principle and molecular dynamics calculations, it captures both the cd-Si → β-Si transition upon compression and the β-Si → a-Si transition upon rapid decompression, which are most relevant for indenter loading. The model was numerically implemented in analogy to incremental plasticity and successfully applied for finite-element (FE) simulations of nanoindentation.     

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.     

Coupled annealing temperature and layer thickness effect on strengthening mechanisms of Ti/Ni multilayer thin films

A systematic study was performed on mechanical and microstructural properties of Ti/Ni multilayers with layer thickness from 200 nm to 6 nm and annealing temperature from room temperature to 500 °C. Based on the observed hardness evolution, a coupled layer-thickness and annealing-temperature dependent strengthening mechanism map is proposed. For as-deposited films, the deformation behavior follows the traditional trend of dislocation mediated strengthening to grain boundary mediated softening with decreasing layer thickness. For annealed films, grain boundary relaxation is considered to be the initial strengthening mechanism with higher activation temperature required for thicker layers. Under further annealing, solid solution hardening, intermetallic precipitation hardening, and fully intermixed alloy structure continue to strengthen the thin layered films, while recrystallization and grain-growth lead to the eventual softening of thick layered films. For the films with intermediate layer thickness, a strong orientation dependent hardness behavior is exhibited under high temperature annealing due to mechanism switch from grain growth softening to intermetallic precipitation hardening when changing the loading orientation from perpendicular to parallel to the layer interfaces.