The impact of sand slugs against beams and plates: Coupled discrete particle/finite element simulations

The impact of a slug of dry sand particles against a metallic sandwich beam or circular sandwich plate is analysed in order to aid the design of sandwich panels for shock mitigation. The sand particles interact via a combined linear-spring-and-dashpot law whereas the face sheets and compressible core of the sandwich beam and plate are treated as rate-sensitive, elastic-plastic solids. The majority of the calculations are performed in two dimensions and entail the transverse impact of end-clamped monolithic and sandwich beams, with plane strain conditions imposed. The sand slug is of rectangular shape and comprises a random loose packing of identical, circular cylindrical particles. These calculations reveal that loading due to the sand is primarily inertial in nature with negligible fluid-structure interaction: the momentum transmitted to the beam is approximately equal to that of the incoming sand slug. For a slug of given incoming momentum, the dynamic deflection of the beam increases with decreasing duration of sand-loading until the impulsive limit is attained. Sandwich beams with thick, strong cores significantly outperform monolithic beams of equal areal mass. This performance enhancement is traced to the "sandwich effect" whereby the sandwich beams have a higher bending strength than that of the monolithic beams. Three-dimensional (3D) calculations are also performed such that the sand slug has the shape of a circular cylindrical column of finite height, and contains spherical sand particles. The 3D slug impacts a circular monolithic plate or sandwich plate and we show that sandwich plates with thick strong cores again outperform monolithic plates of equal areal mass. Finally, we demonstrate that impact by sand particles is equivalent to impact by a crushable foam projectile. The calculations on the equivalent projectile are significantly less intensive computationally, yet give predictions to within 5% of the full discrete particle calculations for the monolithic and sandwich beams and plates. These foam projectile calculations suggest that metallic foam projectiles can be used to simulate the loading by sand particles within a laboratory setting. © 2013 Elsevier Ltd.

Finite horizon LQR control with limited controller-system communication

We consider finite-horizon LQR control with limited controller-system communication. Within a time-horizon T , the controller can only communicate with the system d

Cabaret finite-difference schemes for the one-dimensional Euler equations

In the present paper we consider second order compact upwind schemes with a space split time derivative (CABARET) applied to one-dimensional compressible gas flows. As opposed to the conventional approach associated with incorporating adjacent space cells we use information from adjacent time layer to improve the solution accuracy. Taking the first order Roe scheme as the basis we develop a few higher (i.e. second within regions of smooth solutions) order accurate difference schemes. One of them (CABARET3) is formulated in a two-time-layer form, which makes it most simple and robust. Supersonic and subsonic shock-tube tests are used to compare the new schemes with several well-known second-order TVD schemes. In particular, it is shown that CABARET3 is notably more accurate than the standard second-order Roe scheme with MUSCL flux splitting.

Lateral response of monopiles using centrifuge testing and finite element analysis

Offshore wind capacity is expected to grow exponentially over the next decade resulting in the production of a considerable amount of renewable energy. Monopiles are currently the most popular type of foundation for supporting offshore wind turbines in shallow to medium depth waters. In this paper, the load-deformation response of a 3.8 m diameter monopile installed in soft clays when subjected to axial and lateral loading is investigated using centrifuge testing and soil pore-fluid coupled three-dimensional finite element analysis. Monopile deformation is principally assessed in terms of its lateral displacements and bending moments. Its behaviour as a short rigid pile is discussed using concepts such as its rotation at mudline and the pile depth at which pivoting occurs. © 2014 Taylor & Francis Group.

Finite Element studies on indentation size effect using a higher order strain gradient theory

In this work, a Finite Element implementation of a higher order strain gradient theory (due to Fleck and Hutchinson, 2001) has been used within the framework of large deformation elasto-viscoplasticity to study the indentation of metals with indenters of various geometries. Of particular interest is the indentation size effect (ISE) commonly observed in experiments where the hardness of a range of materials is found to be significantly higher at small depths of indentation but reduce to a lower, constant value at larger depths. That the ISE can be explained by strain gradient plasticity is well known but this work aims to qualitatively compare a gamut of experimental observations on this effect with predictions from a higher order strain gradient theory. Results indicate that many of the experimental observations are qualitatively borne out by our simulations. However, areas exist where conflicting experimental results make assessment of numerical predictions difficult. © 2012 Elsevier Ltd. All rights reserved.

Finite element analysis of the ring rolling process with integrated closed-loop control

In FEA of ring rolling processes the tools' motions usually are defined prior to simulation. This procedure neglects the closed-loop control, which is used in industrial processes to control up to eight degrees of freedom (rotations, feed rates, guide rolls) in real time, taking into account the machine's performance limits as well as the process evolution. In order to close this gap in the new simulation approach all motions of the tools are controlled according to sensor values which are calculated within the FE simulation. This procedure leads to more realistic simulation results in comparison to the machine behaviour. © 2012 CIRP.

Efficient parametric uncertainty analysis within the hybrid Finite Element/Statistical Energy Analysis method

This paper is concerned with the development of efficient algorithms for propagating parametric uncertainty within the context of the hybrid Finite Element/Statistical Energy Analysis (FE/SEA) approach to the analysis of complex vibro-acoustic systems. This approach models the system as a combination of SEA subsystems and FE components; it is assumed that the FE components have fully deterministic properties, while the SEA subsystems have a high degree of randomness. The method has been recently generalised by allowing the FE components to possess parametric uncertainty, leading to two ensembles of uncertainty: a non-parametric one (SEA subsystems) and a parametric one (FE components). The SEA subsystems ensemble is dealt with analytically, while the effect of the additional FE components ensemble can be dealt with by Monte Carlo Simulations. However, this approach can be computationally intensive when applied to complex engineering systems having many uncertain parameters. Two different strategies are proposed: (i) the combination of the hybrid FE/SEA method with the First Order Reliability Method which allows the probability of the non-parametric ensemble average of a response variable exceeding a barrier to be calculated and (ii) the combination of the hybrid FE/SEA method with Laplace's method which allows the evaluation of the probability of a response variable exceeding a limit value. The proposed approaches are illustrated using two built-up plate systems with uncertain properties and the results are validated against direct integration, Monte Carlo simulations of the FE and of the hybrid FE/SEA models. © 2013 Elsevier Ltd.

A fixed-grid b-spline finite element technique for fluid-structure interaction

We present a fixed-grid finite element technique for fluid-structure interaction problems involving incompressible viscous flows and thin structures. The flow equations are discretised with isoparametric b-spline basis functions defined on a logically Cartesian grid. In addition, the previously proposed subdivision-stabilisation technique is used to ensure inf-sup stability. The beam equations are discretised with b-splines and the shell equations with subdivision basis functions, both leading to a rotation-free formulation. The interface conditions between the fluid and the structure are enforced with the Nitsche technique. The resulting coupled system of equations is solved with a Dirichlet-Robin partitioning scheme, and the fluid equations are solved with a pressure-correction method. Auxiliary techniques employed for improving numerical robustness include the level-set based implicit representation of the structure interface on the fluid grid, a cut-cell integration algorithm based on marching tetrahedra and the conservative data transfer between the fluid and structure discretisations. A number of verification and validation examples, primarily motivated by animal locomotion in air or water, demonstrate the robustness and efficiency of our approach. © 2013 John Wiley & Sons, Ltd.

Multi-hit ballistic damage characterisation of 304 stainless steel plates with finite elements

Multi-impact of projectiles on thin 304 stainless steel plates is investigated to assess the degradation of ballistic performance, and to characterise the inherent mechanisms. Assessment of ballistic degradation is by means of a double-impact of rigid spheres at the same site on a circular clamped plate. The limiting velocity of the second impact, will be altered by the velocity of the antecedent impact. Finite element analyses were used to elucidate experimental results and understand the underlying mechanisms that give rise to the performance degradation. The effect of strength and ductility on the single and multi-impact performance was also considered. The model captured the experimental results with excellent agreement. Moreover, the material parameters used within the model were exclusively obtained from published works with no fitting or calibration required. An attempt is made to quantify the elevation of the ballistic limit of thin plates by the dynamic mechanism of travelling hinges. Key conclusions: The multi-hit performance scales linearly with the single-hit performance; and strength is a significantly greater effector of increased ballistic limit than ductility, even at the expense of toughness. © 2014 Elsevier Ltd.

An important pitfall of pseudo-static finite element analysis

Finite Element (FE) pseudo-static analysis can provide a good compromise between simplified methods of dynamic analysis and time domain analysis. The pseudo-static FE approach can accurately model the in situ, stresses prior to seismic loading (when it follows a static analysis simulating the construction sequence) is relatively simple and not as computationally expensive as the time domain approach. However this method should be used with caution as the results can be sensitive to the choice of the mesh dimensions. In this paper two simple examples of pseudo-static finite element analysis are examined parametrically, a homogeneous slope and a cantilever retaining wall, exploring the sensitivity of the pseudo-static analysis results on the adopted mesh size. The mesh dependence was found to be more pronounced for problems with high critical seismic coefficients values (e.g. gentle slopes or small walls), as in these cases a generalised layer failure mechanism is developed simultaneously with the slope or wall mechanism. In general the mesh width was found not to affect notably the predicted value of critical seismic coefficient but to have a major impact on the predicted movements. © 2012 Elsevier Ltd.