Experimental and numerical investigation of fracture mechanics of bitumen beams

This work presents an experimental and numerical investigation to characterise the fracture properties of pure bitumen (the binder in asphalt paving materials). The paper is divided into two parts. The first part describes an experimental study of fracture characterisation parameters of pure bitumen as determined by three-point bend tests. The second part deals with modelling of fracture and failure of bitumen by Finite Element analysis. Fracture mechanics parameters, stress intensity factor, KIC, fracture energy, GIC, and J-integral, JIC, are used for evaluation of bitumen's fracture properties. The material constitutive model developed by Ossa et al. [40,41] which was implemented into a FE code by Costanzi [18] is combined with cohesive zone models (CZM) to simulate the fracture behaviour of pure bitumen. Experimental and numerical results are presented in the form of failure mechanism maps where ductile, brittle and brittle-ductile transition regimes of fracture behaviour are classified. The FE predictions of fracture behaviour match well with experimental results. © 2012 Elsevier Ltd.

First principles interatomic potential for tungsten based on Gaussian process regression

An accurate description of atomic interactions, such as that provided by first principles quantum mechanics, is fundamental to realistic prediction of the properties that govern plasticity, fracture or crack propagation in metals. However, the computational complexity associated with modern schemes explicitly based on quantum mechanics limits their applications to systems of a few hundreds of atoms at most. This thesis investigates the application of the Gaussian Approximation Potential (GAP) scheme to atomistic modelling of tungsten - a bcc transition metal which exhibits a brittle-to-ductile transition and whose plasticity behaviour is controlled by the properties of $\frac{1}{2} \langle 111 \rangle$ screw dislocations. We apply Gaussian process regression to interpolate the quantum-mechanical (QM) potential energy surface from a set of points in atomic configuration space. Our training data is based on QM information that is computed directly using density functional theory (DFT). To perform the fitting, we represent atomic environments using a set of rotationally, permutationally and reflection invariant parameters which act as the independent variables in our equations of non-parametric, non-linear regression. We develop a protocol for generating GAP models capable of describing lattice defects in metals by building a series of interatomic potentials for tungsten. We then demonstrate that a GAP potential based on a Smooth Overlap of Atomic Positions (SOAP) covariance function provides a description of the $\frac{1}{2} \langle 111 \rangle$ screw dislocation that is in agreement with the DFT model. We use this potential to simulate the mobility of $\frac{1}{2} \langle 111 \rangle$ screw dislocations by computing the Peierls barrier and model dislocation-vacancy interactions to QM accuracy in a system containing more than 100,000 atoms.

The effect of hole size upon the strength of metallic and polymeric foams

Tensile and compressive tests have been performed on centre-hole panels, made from three types of metallic foams and two polymeric foams. In compression, the foams fail in a ductile, notch-insensitive manner, in support of a "net section strength" criterion. In tension, a ductile-brittle transition is observed for some of the foams at sufficiently large specimen sizes: for a small hole diameter the net section strength criterion is obeyed, whereas for a large hole a local stress criterion applies and the net section strength is reduced. For a number of the foams, the panel size was not sufficiently large to observe this ductile-brittle switch in behaviour. The predictions of a cohesive zone model are compared with the measured strengths and are found to be in good agreement. © 2001 Elsevier Science Ltd. All rights reserved.

Modelling and Fragility Analysis of Non-Ductile Reinforced Concrete Buildings in Low-to-Moderate Seismic Zones

Reinforced concrete buildings in low-to-moderate seismic zones are often designed only for gravity loads in accordance with the non-seismic detailing provisions. Deficient detailing of columns and beam-column joints can lead to unpredictable brittle failures even under moderate earthquakes. Therefore, a reliable estimate of structural response is required for the seismic evaluation of these structures. For this purpose, analytical models for both interior and exterior slab-beam-column subassemblages and for a 1/3 scale model frame were implemented into the nonlinear finite element platform OpenSees. Comparison between the analytical results and experimental data available in the literature is carried out using nonlinear pushover analyses and nonlinear time history analysis for the subassemblages and the model frame, respectively. Furthermore, the seismic fragility assessment of reinforced concrete buildings is performed on a set of non-ductile frames using nonlinear time history analyses. The fragility curves, which are developed for various damage states for the maximum interstory drift ratio are characterized in terms of peak ground acceleration and spectral acceleration using a suite of ground motions representative of the seismic hazard in the region.

Fracture mechanics of idealised bituminous mixes

© 2014 Taylor & Francis. The durability of asphalt pavements is strongly impaired by cracks, caused primarily by traffic loads and environmental effects. In this work, fracture behaviour of idealised asphalt mixes is investigated. Experiments on idealised asphalt mixes under pure-tension mode (mode I cracking) were performed and fracture parameters were evaluated. In these three-point bend fracture tests, the test variables were temperature and load rate. The test data were stored in an asphalt materials database and special-purpose tools were implemented to analyse and handle the laboratory data automatically. Fracture mechanism maps were constructed, showing the conditions associated with ductile, brittle and ductile-brittle transition regimes of behaviour. The mechanism maps show the failure response of the material in terms of the stress intensity factor, strain energy release rate and J-integral as a function of the temperature-compensated crack mouth opening strain rate. Fracture behaviour of asphalt mix specimens was simulated by cohesive zone model in conjunction with a novel material constitutive model for asphalt mixes. The finite element model agrees well with the experimental results and provides insights into fracture response of the notched asphalt mix beam specimens.

Mode I crack growth resistance of metallic foams

A Dugdale-type cohesive zone model is used to predict the mode I crack growth resistance (R-curve) of metallic foams, with the fracture process characterized by an idealized traction-separation law that relates the crack surface traction to crack opening displacement. A quadratic yield function, involving the von Mises effective stress and mean stress, is used to account for the plastic compressibility of metallic foams. Finite element calculations are performed for the crack growth resistance under small scale yielding and small scale bridging in plane strain, with K-field boundary conditions. The following effects upon the fracture process are quantified: material hardening, bridging strength, T-stress (the non-singular stress acting parallel to the crack plane), and the shape of yield surface. To study the failure behaviour and notch sensitivity of metallic foams in the presence of large scale yielding, a study is made for panels embedded with either a centre-crack or an open hole and subjected to tensile stressing. For the centre-cracked panel, a transition crack size is predicted for which the fracture response switches from net section yielding to elastic-brittle fracture. Likewise, for a panel containing a centre-hole, a transition hole diameter exists for which the fracture response switches from net section yielding to a local maximum stress criterion at the edge of the hole.

On the transition from displacement to mixing ventilation with a localized heat source

We investigate the steady state natural ventilation of a room heated at the base and consisting of two vents at different levels. We explore how the air flow rate and internal temperature relative to the exterior vary as a function of the vent areas, position of the vents and heat load in order to establish appropriate ventilation strategies for a room. When the room is heated by a distributed source, the room becomes well mixed and the steady state ventilation rate depends on the heating rate, the area of the vents and the distance between the lower and upper level vents. However, when the room is heated by a localised source the room becomes stratified. If the effective ventilation area is sufficiently large, then the interface separating the two layers lies above the inlet vent and the lower layer is comprised of ambient fluid. In this case the upper layer is warmer than in the well mixed case and the ventilation rate is smaller. However, if the effective area for ventilation is sufficiently small, then the interface separating the two layers lies below the inlet vent and the lower layer is comprised of warm fluid which originates as the cold incoming fluid mixes during descent from the vent through the upper layer. In this case both the ventilation rate and the upper layer temperature are the same as in the case of a distributed heat load. As the vertical separation between lower and upper level vents decreases, then the temperature difference between the layers falls to zero and the room becomes approximately well mixed. These findings suggest how the appropriate ventilation strategy for a room can be varied depending on the exterior temperature, with mixing ventilation more suitable for winter conditions and displacement ventilation for warmer external temperatures.

Experimental investigation of a control method for SI-HCCI-SI transition in a multi-cylinder gasoline engine

In HCCI engines, the Air/Fuel Ratio (AFR) and Residual Gas Fraction (RGF) are difficult to control during the SI-HCCI-SI transition, and this may result in incomplete combustion and/or high pressure raise rates. As a result, there may be undesirably high engine load fluctuations. The objectives of this work are to further understand this process and develop control methods to minimize these load fluctuations. This paper presents data on instantaneous AFR and RGF measurements, both taken by novel experimental techniques. The data provides an insight into the cyclic AFR and RGF fluctuations during the switch. These results suggest that the relatively slow change in the intake Manifold Air Pressure (MAP) and actuation time of the Variable Valve Timing (VVT) are the main causes of undesired AFR and RGF fluctuations, and hence an unacceptable Net IMEP (NIMEP) fluctuation. We also found large cylinder-to-cylinder AFR variations during the transition. Therefore, besides throttle opening control and VVT shifting, cyclic and individual cylinder fuel injection control is necessary to achieve a smooth transition. The control method was developed and implemented in a test engine, and the result was a considerably reduced NIMEP fluctuation during the mode switch. The instantaneous AFR and RGF measurements could furthermore be adopted to develop more sophisticated control methods for SI-HCCI-SI transitions. © 2010 SAE International.

Modal scattering at an impedance transition in a lined flow duct

An explicit Wiener-Hopf solution is derived to describe the scattering of duct modes at a hard-soft wall impedance transition in a circular duct with uniform mean flow. Specifically, we have a circular duct r = 1, - ∞ < x < ∞ with mean flow Mach number M > 0 and a hard wall along x < 0 and a wall of impedance Z along x > 0. A minimum edge condition at x = 0 requires a continuous wall streamline r = 1 + h(x, t), no more singular than h = Ο(x1/2) for x ↓ 0. A mode, incident from x < 0, scatters at x = 0 into a series of reflected modes and a series of transmitted modes. Of particular interest is the role of a possible instability along the lined wall in combination with the edge singularity. If one of the "upstream" running modes is to be interpreted as a downstream-running instability, we have an extra degree of freedom in the Wiener-Hopf analysis that can be resolved by application of some form of Kutta condition at x = 0, for example a more stringent edge condition where h = Ο(x3/2) at the downstream side. The question of the instability requires an investigation of the modes in the complex frequency plane and therefore depends on the chosen impedance model, since Z = Z (ω) is essentially frequency dependent. The usual causality condition by Briggs and Bers appears to be not applicable here because it requires a temporal growth rate bounded for all real axial wave numbers. The alternative Crighton-Leppington criterion, however, is applicable and confirms that the suspected mode is usually unstable. In general, the effect of this Kutta condition is significant, but it is particularly large for the plane wave at low frequencies and should therefore be easily measurable. For ω → 0, the modulus fends to |R001| → (1 + M)/(1 -M) without and to 1 with Kutta condition, while the end correction tends to ∞ without and to a finite value with Kutta condition. This is exactly the same behaviour as found for reflection at a pipe exit with flow, irrespective if this is uniform or jet flow.