A bond model for DEM simulation of cementitious materials and deformable structures

There is an increasing use of the discrete element method (DEM) to study cemented (e.g. concrete and rocks) and sintered particulate materials. The chief advantage of the DEM over continuum based techniques is that it does not make assumptions about how cracking and fragmentation initiate and propagate, since the DEM system is naturally discontinuous. The ability for the DEM to produce a realistic representation of a cemented granular material depends largely on the implementation of an inter-particle bonded contact model. This paper presents a new bonded contact model based on the Timoshenko beam theory which considers axial, shear and bending behaviour of the bond. The bond model was first verified by simulating both the bending and dynamic response of a simply supported beam. The loading response of a concrete cylinder was then investigated and compared with the Eurocode equation prediction. The results show significant potential for the new model to produce satisfactory predictions for cementitious materials. A unique feature of this model is that it can also be used to accurately represent many deformable structures such as frames and shells, so that both particles and structures or deformable boundaries can be described in the same DEM framework. 

Physical rehabilitation interventions for adult intensive care patients across the continuum of recovery (i.e. within the intensive care unit (ICU) and following ICU and hospital discharge): an overview of systematic reviews

Background
Patients admitted to the intensive care unit with critical illness often experience significant physical impairments, which typically persist for many years following resolution of the original illness. Physical rehabilitation interventions that enhance restoration of physical function have been evaluated across the continuum of recovery following critical illness including within the intensive care unit, following discharge to the ward and beyond hospital discharge. Multiple systematic reviews have been published appraising the expanding evidence investigating these physical rehabilitation interventions, although there appears to be variability in review methodology and quality. We aim to conduct an overview of existing systematic reviews of physical rehabilitation interventions for adult intensive care patients across the continuum of recovery.

Methods/design
This protocol has been developed according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses Protocol (PRISMA-P) guidelines. We will search the Cochrane Systematic Review Database, Database of Abstracts of Reviews of Effectiveness, Cochrane Central Register of Controlled Trials, MEDLINE, Excerpta Medica Database and Cumulative Index to Nursing and Allied Health Literature databases. We will include systematic reviews of randomised controlled trials of adult patients, admitted to the intensive care unit and who have received physical rehabilitation interventions at any time point during their recovery. Data extraction will include systematic review aims and rationale, study types, populations, interventions, comparators, outcomes and quality appraisal method. Primary outcomes of interest will focus on findings reflecting recovery of physical function. Quality of reporting and methodological quality will be appraised using the PRISMA checklist and the Assessment of Multiple Systematic Reviews tool.

Discussion
We anticipate the findings from this novel overview of systematic reviews will contribute to the synthesis and interpretation of existing evidence regarding physical rehabilitation interventions and physical recovery in post-critical illness patients across the continuum of recovery.

Improving composite damage modelling through automatic placement of cohesive elements

Numerous experimental studies of damage in composite laminates have shown that intralaminar (in-plane) matrix cracks lead to interlaminar delamination (out-of-plane) at ply interfaces. The smearing of in-plane cracks over a volume, as a consequence of the use of continuum damage mechanics, does not always effectively capture the full extent of the interaction between the two failure mechanisms. A more accurate representation is obtained by adopting a discrete crack approach via the use of cohesive elements, for both in-plane and out-of-plane damage. The difficulty with cohesive elements is that their location must be determined a priori in order to generate the model; while ideally the position of the crack migration, and more generally the propagation path, should be obtained as part of the problem’s solution. With the aim of enhancing current modelling capabilities with truly predictive capabilities, a concept of automatic insertion of interface elements is utilized. The consideration of a simple traction criterion in relation to material strength, evaluated at each node of the model (or of the regions of the model where it is estimated cracks might form), allows for the determination of initial crack location and subsequent propagation by the insertion of cohesive elements during the course of the analysis. Several experimental results are modelled using the commercial package ABAQUS/Standard with an automatic insertion subroutine developed in this work, and the results are presented to demonstrate the capabilities of this technique.

Numerical analysis of complex failure mechanisms in composite panels

A three-dimensional continuum damage mechanics-based material model has been implemented in an implicit Finite Element code to simulate the progressive degradation of advanced composite materials. The damage model uses seven damage variables assigned to tensile, compressive and non-linear shear damage at a laminae level. The objectivity of the numerical discretization is assured using a smeared formulation. The material model was benchmarked against experimental uniaxial coupon tests and it is shown to reproduce key aspects observable during failure, such as the inclined fracture plane in matrix compression and the shear band in a ±45° tension specimen.

Numerical considerations for the implementation of a comprehensive composite damage model

A comprehensive continuum damage mechanics model [1] had been developed to capture the detailed
behaviour of a composite structure under a crushing load. This paper explores some of the difficulties
encountered in the implementation of this model and their mitigation. The use of reduced integration
element and a strain softening model both negatively affect the accuracy and stability of the
simulation. Damage localisation effects demanded an accurate measure of characteristic length. A
robust algorithm for determining the characteristic length was implemented. Testing showed that this
algorithm produced marked improvements over the use of the default characteristic length provided
by Abaqus. Zero-energy or hourglass modes, in reduced integration elements, led to reduced
resistance to bending. This was compounded by the strain softening model, which led to the formation
of elements with little resistance to deformation that could invert if left unchecked. It was shown,
through benchmark testing, that by deleting elements with excess distortions and controlling the mesh
using inbuilt distortion/hourglass controls, these issues can be alleviated. These techniques
contributed significantly to the viability and usability of the damage model.

Validation of a 3D damage model for predicting the response of composite structures under crushing loads

A 3D intralaminar continuum damage mechanics based material model, combining damage mode interaction and material nonlinearity, was developed to predict the damage response of composite structures undergoing crush loading. This model captures the structural response without the need for calibration of experimentally determined material parameters. When used in the design of energy absorbing composite structures, it can reduce the dependence on physical testing. This paper validates this model against experimental data obtained from the literature and in-house testing. Results show that the model can predict the force response of the crushed composite structures with good accuracy. The simulated energy absorption in each test case was within 12% of the experimental value. Post-crush deformation and the damage morphologies, such as ply splitting, splaying and breakage, were also accurately reproduced. This study establishes the capability of this damage model for predicting the responses of composite structures under crushing loads.

M Dwarf Flare Continuum Variations on One-second Timescales: Calibrating and Modeling of ULTRACAM Flare Color Indices

We present a large data set of high-cadence dMe flare light curves obtained with custom continuum filters on the triple-beam, high-speed camera system ULTRACAM. The measurements provide constraints for models of the near-ultraviolet (NUV) and optical continuum spectral evolution on timescales of ≈1 s. We provide a robust interpretation of the flare emission in the ULTRACAM filters using simultaneously obtained low-resolution spectra during two moderate-sized flares in the dM4.5e star YZ CMi. By avoiding the spectral complexity within the broadband Johnson filters, the ULTRACAM filters are shown to characterize bona fide continuum emission in the NUV, blue, and red wavelength regimes. The NUV/blue flux ratio in flares is equivalent to a Balmer jump ratio, and the blue/red flux ratio provides an estimate for the color temperature of the optical continuum emission. We present a new “color-color” relationship for these continuum flux ratios at the peaks of the flares. Using the RADYN and RH codes, we interpret the ULTRACAM filter emission using the dominant emission processes from a radiative-hydrodynamic flare model with a high nonthermal electron beam flux, which explains a hot, T ≈ 104 K, color temperature at blue-to-red optical wavelengths and a small Balmer jump ratio as observed in moderate-sized and large flares alike. We also discuss the high time resolution, high signal-to-noise continuum color variations observed in YZ CMi during a giant flare, which increased the NUV flux from this star by over a factor of 100. Based on observations obtained with the Apache Point Observatory 3.5 m telescope, which is owned and operated by the Astrophysical Research Consortium, based on observations made with the William Herschel Telescope operated on the island of La Palma by the Isaac Newton Group in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofsica de Canarias, and observations, and based on observations made with the ESO Telescopes at the La Silla Paranal Observatory under programme ID 085.D-0501(A).