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.

Finite element analysis of compressor disc heat transfer using the transient conduction method with sparse boundary collocation points

Surface temperature measurements from two discs of a gas turbine compressor rig are used as boundary conditions for the transient conduction solution (inverse heat transfer analysis). The disc geometry is complex, and so the finite element method is used. There are often large radial temperature gradients on the discs, and the equations are therefore solved taking into account the dependence of thermal conductivity on temperature. The solution technique also makes use of a multigrid algorithm to reduce the solution time. This is particularly important since a large amount of data must be analyzed to obtain correlations of the heat transfer. The finite element grid is also used for a network analysis to calculate the radiant heat transfer in the cavity formed between the two compressor discs. The work discussed here proved particularly challenging as the disc temperatures were only measured at four different radial locations. Four methods of surface temperature interpolation are examined, together with their effect on the local heat fluxes. It is found that the choice of interpolation method depends on the available number of data points. Bessel interpolation gives the best results for four data points, whereas cubic splines are preferred when there are considerably more data points. The results from the analysis of the compressor rig data show that the heat transfer near the disc inner radius appears to be influenced by the central throughflow. However, for larger radii, the heat transfer from the discs and peripheral shroud is found to be consistent with that of a buoyancy-induced flow.

Optimal control with stochastic PDE constraints and uncertain controls

The optimal control of problems that are constrained by partial differential equations with uncertainties and with uncertain controls is addressed. The Lagrangian that defines the problem is postulated in terms of stochastic functions, with the control function possibly decomposed into an unknown deterministic component and a known zero-mean stochastic component. The extra freedom provided by the stochastic dimension in defining cost functionals is explored, demonstrating the scope for controlling statistical aspects of the system response. One-shot stochastic finite element methods are used to find approximate solutions to control problems. It is shown that applying the stochastic collocation finite element method to the formulated problem leads to a coupling between stochastic collocation points when a deterministic optimal control is considered or when moments are included in the cost functional, thereby forgoing the primary advantage of the collocation method over the stochastic Galerkin method for the considered problem. The application of the presented methods is demonstrated through a number of numerical examples. The presented framework is sufficiently general to also consider a class of inverse problems, and numerical examples of this type are also presented. © 2011 Elsevier B.V.

Shear-flexible subdivision shells

We present a new shell model and an accompanying discretisation scheme that is suitable for thin and thick shells. The deformed configuration of the shell is parameterised using the mid-surface position vector and an additional shear vector for describing the out-of-plane shear deformations. In the limit of vanishing thickness, the shear vector is identically zero and the Kirchhoff-Love model is recovered. Importantly, there are no compatibility constraints to be satisfied by the shape functions used for discretising the mid-surface and the shear vector. The mid-surface has to be interpolated with smooth C 1-continuous shape functions, whereas the shear vector can be interpolated with C 0-continuous shape functions. In the present paper, the mid-surface as well as the shear vector are interpolated with smooth subdivision shape functions. The resulting finite elements are suitable for thin and thick shells and do not exhibit shear locking. The good performance of the proposed formulation is demonstrated with a number of linear and geometrically non-linear plate and shell examples. © 2012 John Wiley & Sons, Ltd.

Predicting the acoustics of squeak and rattle

This paper discusses the development of a computationally efficient numerical method for predicting the acoustics of rattle events upfront in the design cycle. The method combines Finite Elements, Boundary Elements and SEA and enables the loudness of a large number of rattle events to be efficiently predicted across a broad frequency range. A low frequency random vibro-acoustic model is used in conjunction with various closed form analytical expressions in order to quickly predict impact probabilities and locations. An existing method has been extended to estimate the statistics of the contact forces across a broad frequency range. Finally, broadband acoustic radiation is predicted using standard low, mid and high frequency vibro-acoustic methods and used to estimate impact loudness. The approach is discussed and a number of validation examples are presented.

Wave propagation under vertically excited surface foundation

Recently developed equipment allows measurement of the shear modulus of soil in situ as a function of level of strain. In these field experiments, the excitation is applied on the ground surface using large scale shakers, and the response of the soil deposit is recorded through embedded receivers. The focus of this paper is on the simulation of signals which would be recorded at the receiver locations in idealized conditions to provide guidelines on the interpretation of field measurements. Discrete and finite element methods are employed to model one dimensional and three dimensional geometries, respectively, under various lateral boundary conditions. When the first times of arrival are detected by receivers under the vertical impulse, they coincide with the arrival of the P wave, related to the constrained modulus of the material, regardless of lateral boundary conditions. If one considers, on the other hand, phase differences between the motions at two receivers the picture is far more complicated and one would obtain propagation velocities, function of frequency and depth, which do not correspond to either the constrained modulus or Young's modulus. It is thus necessary to apply some care when interpreting the data from field tests based on vertical steady state vibrations. The use of inverse analysis can be considered as a way of extracting the shear modulus of soil from the field test measurements. © 2008 ASCE.

What is the weakest topology in which feedback stability is robust?

Mathematical theorems in control theory are only of interest in so far as their assumptions relate to practical situations. The space of systems with transfer functions in ℋ∞, for example, has many advantages mathematically, but includes large classes of non-physical systems, and one must be careful in drawing inferences from results in that setting. Similarly, the graph topology has long been known to be the weakest, or coarsest, topology in which (1) feedback stability is a robust property (i.e. preserved in small neighbourhoods) and (2) the map from open-to-closed-loop transfer functions is continuous. However, it is not known whether continuity is a necessary part of this statement, or only required for the existing proofs. It is entirely possible that the answer depends on the underlying classes of systems used. The class of systems we concern ourselves with here is the set of systems that can be approximated, in the graph topology, by real rational transfer function matrices. That is, lumped parameter models, or those distributed systems for which it makes sense to use finite element methods. This is precisely the set of systems that have continuous frequency responses in the extended complex plane. For this class, we show that there is indeed a weaker topology; in which feedback stability is robust but for which the maps from open-to-closed-loop transfer functions are not necessarily continuous. © 2013 Copyright Taylor and Francis Group, LLC.

Unified form language: A domain-specific language for weak formulations of partial differential equations

We present the Unified Form Language (UFL), which is a domain-specific language for representing weak formulations of partial differential equations with a view to numerical approximation. Features of UFL include support for variational forms and functionals, automatic differentiation of forms and expressions, arbitrary function space hierarchies formultifield problems, general differential operators and flexible tensor algebra. With these features, UFL has been used to effortlessly express finite element methods for complex systems of partial differential equations in near-mathematical notation, resulting in compact, intuitive and readable programs. We present in this work the language and its construction. An implementation of UFL is freely available as an open-source software library. The library generates abstract syntax tree representations of variational problems, which are used by other software libraries to generate concrete low-level implementations. Some application examples are presented and libraries that support UFL are highlighted. © 2014 ACM.