Search results

Simple discrete element models using PFC^2D models with bonded assemblies of particles were used to numerically simulate direct tension and block punching tests on thin spray‐on tunnel liner materials to gain insight about the liner support mechanisms. PFC^2D input parameters were calibrated such that the rupture load and elongation at rupture were similar to the laboratory test data. The calibrated model of the liner material was then used to simulate a liner around a highly stressed tunnel in rock where stresses caused extensive fracturing near the top of the tunnel. The effect of the liner was analysed by modelling the tunnel with and without the liner and showed that the liner had minimal impact on fracture propagation in the rock because of the liner's highly deformable nature. However, the liner was able to retain the fractured rock in place.

This paper discusses the issues involved in the development of combined finite/discrete element methods; both from a fundamental theoretical viewpoint and some related algorithmic considerations essential for the efficient numerical solution of large scale industrial problems. The finite element representation of the solid region is combined with progressive fracturing, which leads to the formation of discrete elements, which may be composed of one or more deformable finite elements. The applicability of the approach is demonstrated by the solution of a range of examples relevant to various industrial sections.

The purpose of this paper is to propose a novel combined finite-discrete element method (FDEM), based on the cohesive zone model, for simulating rockslide problems at the laboratory scale.

The combined FDEM is realized using ABAQUS/Explicit. The rock mass is represented as a collection of elastic bulk elements glued by cohesive elements with zero thickness. To reproduce the tensile and shear micro-fractures in rock material, the Mohr-Coulomb model with tension cut-off is employed as the damage initiation criterion of cohesive elements. Three simulated laboratory tests are considered to verify the capability of combined FDEM in reproducing the mechanical behavior of rock masses. Three slope models with different joint inclinations are taken to illustrate the application of the combined FDEM to rockslide simulation.

The results show that the joint inclination is an important factor for inducing the progressive failure behavior. With a low joint inclination, the slope failure process is observed to be a collapse mode. As the joint inclination becomes higher, the failure mode changes to sliding and the steady time of rock blocks is shortened. Moreover, the runout distance and post-failure slope angle decrease as the joint inclination increases.

These studies indicate that the combined FDEM performed within ABAQUS can simulate slope stability problems for research purposes and is useful for studying the slope failure mechanism comprehensively.

Describes a probabilistic methodology for fracture assessments of flawed structures constructed of ferritic steels using the research code WSTRESS. The probabilistic formulation for cleavage fracture implements a multiaxial form of the weakest link model which couples the macroscopic fracture behavior with a micromechanics model based on the statistics of microcracks. The Weibull stress, σ_w, emerges as a suitable near‐tip parameter to provide a connection between the microregime of failure and remote loading (J). WSTRESS builds on an iterative procedure to incorporate a 3‐D finite element description of the crack‐tip stress field and measured values of fracture toughness to calibrate the Weibull modulus, m, and the scale parameter, σ_u. Specific features of the code include statistical inference of Weibull parameters based on uncensored and censored models (with maximum likelihood method), construction of confidence intervals, several definitions for the near‐tip fracture process zone and other general facilities such as spatial integration of the stress field (to incorporate the random orientation of microcracks) and stochastic simulation of fracture data using the Monte Carlo method. The code also includes a convenient free‐form command language and a seamless interface with finite element results files stored in Patran binary or ASCII format.

The purpose of this paper is to present a 2D combined finite‐discrete element method (FEM‐DEM) to model the multi‐fracture of beam structures and an application of the method to an ice‐structure interaction problem.

In the method, elastic beams and their fracture are modelled according to FEM by using nonlinear Timoshenko beam elements and cohesive crack model. Additionally, the beam elements are used to tie the discrete elements together. The contact forces between the colliding beams are calculated by using the DEM.

Three numerical examples are given to verify the method. Further, the method is applied to model the failure process of a floating ice beam against an inclined structure. Based on the comparison of the experiments and the simulation, a good agreement between the results is observed.

In the context of combined FEM‐DEM, the two novel features presented in this paper are: the use of Timoshenko finite element beams with damping to calculate internal forces and to combine the discrete elements; and the bending failure by the cohesive crack approach while simultaneously keeping track of the position of the neutral axis of the beam.

To contribute for a reliable estimation of the compressive strength of unreinforced masonry from the properties of the constituents (units and mortar).

Sophisticated non‐linear continuum models, based on damage, plasticity, cracking or other formulation, are today standard in several finite element programs. The adequacy of such models to provide reliable estimates of masonry compressive strength, from the properties of the constituents, remains unresolved. The authors have shown recently that continuum models might significantly overestimate the prediction of the compressive strength. Hence, an alternative phenomenological approach developed in a discrete framework is proposed, based on attributing to masonry components a fictitious micro‐structure composed of linear elastic particles separated by non‐linear interface elements. The model is discussed in detail and a comparison with experimental results and numerical results using a standard continuum model is provided.

Clear advantages in terms of compressive strength and peak strain prediction were found using the particle model when compared with standard continuum models. Moreover, compressive and tensile strength values provided by the model were found to be particle size‐ and particle distortion‐independent for practical purposes. It is also noted that size‐dependent responses were obtained and that shear parameters rather than tensile parameters were found to play a major role at the meso‐level of the phenomenological model.

This paper provides further insight into the compressive behaviour of quasi‐brittle materials, with an emphasis on the strength prediction of masonry composites. Reliable prediction of masonry strength is of great use in the civil engineering field, allowing one to reduce experimental testing in expensive wallets and to avoid the usage of conservative empirical formulae.

This paper aims to present a multi-scale stochastic damage model (SDM) for concrete and apply it to the stochastic response analysis of reinforced concrete shear wall structures.

The proposed SDM is constructed at two scales, i.e. the macro-scale and the micro-scale. The general framework of the SDM is established on the basis of the continuum damage mechanics (CDM) at the macro-scale, whereas the detailed damage evolution is determined through a parallel fiber buddle model at the micro-scale. The parallel buddle model is made up of micro-elements with stochastic fracture strains, and a one-dimensional random field is assumed for the fracture strain distribution. To represent the random field, a random functional method is adopted to quantify the stochastic damage evolution process with only two variables; thus, the numerical efficiency is greatly enhanced. Meanwhile, the probability density evolution method (PDEM) is introduced for the structural stochastic response analysis.

By combing the SDM and PDEM, the probabilistic analysis of a shear wall structure is performed. The mean value, standard deviation and the probability density function of the shear wall responses, e.g., shear capacity, accumulated energy consumption and damage evolution, are obtained.

It is noted that the proposed method can reflect the influences of randomness from material level to structural level, and is efficient for stochastic response determination of shear wall structures.

Simple but also accurate models are needed to predict the failure response of concrete structures. Simplicity involves modelling assumptions while accuracy involves objectivity of both the experimentally identified model parameters and the numerica results. For concrete‐like heterogeneous and brittle materials, the modelling assumptions idealizing the material as a homogeneous continuum with classical linear or non‐linear behaviour, leads to some problems at the identification stage, namely the size effect phenomena. A continuum damage model, representing the non‐linear behaviour due to microcracking, is proposed here for predictive computations of structural responses. A Weibull based theory is used to determine, in a statistical sense, the value of the initial damage threshold. The essential influence of material heterogeneity on the damage evolution, is accounted for with a bi‐scale approach which is based on the idea of the non‐local continuum with local strain. It has already established that the non‐local approaches yield realistic failure predictions and the numerical results are convergent for subsequent mesh refinements. The applications presented here show the ability of the approach to predict the failure response of concrete structures without being obscured by size effect problems.

In order to help explain experimental findings related to the stabbing- and ballistic-penetration resistance of flexible body-armor, single-yarn pull-out tests, involving specially prepared fabric-type test coupons, are often carried out. The purpose of this paper is to develop a finite-element-based computational framework for the simulation of the single-yarn pull-out test, and applied to the case of Kevlar® KM2 fabric.

Three conditions of the fabric are considered: neat, i.e, as-woven; polyethylene glycol (PEG)-infiltrated; and shear-thickening fluid (STF)-infiltrated. Due to differences in the three conditions of the fabric, the computational framework had to utilize three different finite-element formulations: standard Lagrangian formulation for the neat fabric; combined Eulerian-Lagrangian formulation for the PEG-infiltrated fabric (an Eulerian subdomain had to be used to treat the PEG solvent/dispersant); and combined continuum Lagrangian/discrete-particle formulation for the STF-infiltrated fabric (to account for the interactions of the particles suspended in PEG, which give rise to the STF character of the suspension, with the yarns, the particles had to be treated explicitly).

The results obtained for the single-yarn pull-out virtual tests are compared with the authors’ experimental counterparts, and a reasonably good agreement is obtained, for all three conditions of the fabric.

To the authors’ knowledge, the present work represents the first attempt to simulate single-yarn pull-out tests of Kevlar® KM2 fabric.

Ballistic missile-resistant impact panels have seen fair advancement over the years, especially in military applications. However, high cost, as well as a changing materials landscape, has impressed the need for a deeper understanding of impact mechanism as well as of new permutations in design strategy development. Parameters such as projectile impact angle, panel impact location as well as application of multilayer sandwich panels are not fully explored and characterised. In this work, finite element method simulation methodology is used on a 25 mm by 25 mm plate of 3.5 mm thickness to investigate the above-mentioned parameters and conditions. Solid elements using Johnson–Cook damage material models are developed. Two common impact angles of 90 and 45° at centre and plate-edge locations are investigated for single-layer titanium alloy and carbon steel panels. Subsequently, a bilayer panel comprising of titanium alloy at the impact layer with the same overall plate thickness is investigated for impact at five different impact speed (ranging from 100 ms^-1 to 500 ms^-1). The displacements and von Mises stresses are documented for all cases, and it is shown that angular impact angles bring about greater plastic deformations as well as higher fracture likelihood compared to normal angle impact. Findings also indicate that with an addition of 1 mm thick Ti-6Al-4V front bilayer, the impact resistance of the high carbon steel is significantly improved (up to twice the impact load), especially at higher impact velocities. The study documents the properties of titanium alloy–carbon steel bilayer armoured panel, which shows good promise for its implementation due to its superior performance and its cost-savings potential.

In this work, finite element method simulation methodology is used to investigate the above-mentioned parameters and conditions. Solid elements using Johnson–Cook damage material models are developed. Two common impact angles 90 and 45° at centre and plate-edge locations are investigated for single-layer titanium alloy and carbon steel panels, and, subsequently, a bilayer panel comprising of titanium alloy for the outer layer is investigated for the combination of the same aforementioned materials. Five different impact speed effects are studied.

The effects and trends of displacements and stresses are documented for all cases and shown to indicate angular impact angles bringing about greater plastic deformations as well as higher fracture likelihood compared to normal angle impact. Findings also show that with an addition of 1 mm thick Ti-6Al-4V front bilayer, the impact resistance of the high carbon steel is significantly improved, especially at higher impact velocities.

The study documents the properties of titanium alloy–carbon steel bilayer armoured panel, which shows good promise for its implementation due to its superior performance and its cost-savings potential.