H∞ sampled-data synthesis and related numerical issues

In this paper we present a new, compact derivation of state-space formulae for the so-called discretisation-based solution of the H∞ sampled-data control problem. Our approach is based on the established technique of continuous time-lifting, which is used to isometrically map the continuous-time, linear, periodically time-varying, sampled-data problem to a discretetime, linear, time-invariant problem. State-space formulae are derived for the equivalent, discrete-time problem by solving a set of two-point, boundary-value problems. The formulae accommodate a direct feed-through term from the disturbance inputs to the controlled outputs of the original plant and are simple, requiring the computation of only a single matrix exponential. It is also shown that the resultant formulae can be easily re-structured to give a numerically robust algorithm for computing the state-space matrices. © 1997 Elsevier Science Ltd. All rights reserved.

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.

Robust multiplexed MPC for distributed multi-agent systems

This paper extends the authors' earlier work which adapted robust multiplexed MPC for application to distributed control of multi-agent systems with non-interacting dynamics and coupled constraint sets in the presence of persistent unknown, but bounded disturbances. Specifically, we propose exploiting the single agent update nature of the multiplexed approach, and fix the update sequence to enable input move-blocking and increased discretisation rates. This permits a higher rate of individual policy update to be achieved, whilst incurring no additional computational cost in the corresponding optimal control problems to be solved. A disturbance feedback policy is included between updates to facilitate finding feasible solutions. The new formulation inherits the property of rapid response to disturbances from multiplexing the control and numerical results show that fixing the update sequence does not incur any loss in performance. © 2011 IFAC.

An integrated model of rolling tyre vibrations on a smooth road

We describe new results on the vibrations of rolling tyres, aimed at noise prediction for tyres of given design on a smooth road surface. This new approach incorporates our existing models, of smooth road-tyre interaction and belt vibration but includes additional features that are required for real tyre patterns. To this end, the model allows variable tread block size and grooves along the belt circumference; the density and angle of these grooves may also vary laterally. The key innovation is to treat the tyre belt as a laterally stacked series of rings, each of which is equipped with a set of viscoelastic springs around its circumference. It is shown how to use this construction to mimic the details of actual tyre patterns and, in conjunction with existing models, predict belt vibrations. The construction is applied to develop a ring discretisation for a real tyre that shows strong lateral variations. It is shown that the vibration amplitude is concentrated on a set of parallel lines in frequency-wavenumber space and that the tread pattern dictates the occurrence and spacing of these lines. Linkage to a boundary element calculation then allows quantification of the influence of tread parameters on radiated noise. Keywords: Vibration, tread pattern, tyre noise. Copyright © (2011) by the Institute of Noise Control Engineering.

Development of realistic liquid equations of state for multi-phase CFD modelling

Several equations of state (EOS) have been incorporated into a novel algorithm to solve a system of multi-phase equations in which all phases are assumed to be compressible to varying degrees. The EOSs are used to both supply functional relationships to couple the conservative variables to the primitive variables and to calculate accurately thermodynamic quantities of interest, such as the speed of sound. Each EOS has a defined balance of accuracy, robustness and computational speed; selection of an appropriate EOS is generally problem-dependent. This work employs an AUSM+-up method for accurate discretisation of the convective flux terms with modified low-Mach number dissipation for added robustness of the solver. In this paper we show a newly-developed time-marching formulation for temporal discretisation of the governing equations with incorporated time-dependent source terms, as well as considering the system of eigenvalues that render the governing equations hyperbolic.

A comparison of two models for the vibration response of piled foundations to inertial and underground-railway-induced loadings

The vibration response of piled foundations due to ground-borne vibration produced by an underground railway is a largely-neglected area in the field of structural dynamics. However, this continues to be an important aspect of research as it is expected that the presence of piled foundations can have a significant influence on the propagation and transmission of the wavefield produced by the underground railway. This paper presents a comparison of two methods that can be employed in calculating the vibration response of a piled foundation: an efficient semi-analytical model, and a Boundary Element model. The semi-analytical model uses a column or an Euler beam to model the pile, and the soil is modelled as a linear, elastic continuum that has the geometry of a thick-walled cylinder with an infinite outer radius and an inner radius equal to the radius of the pile. The boundary element model uses a constant-element BEM formulation for the halfspace, and a rectangular discretisation of the circular pile-soil interface. The piles are modelled as Timoshenko beams. Pile-soil-pile interactions are inherently accounted for in the BEM equations, whereas in the semi-analytical model these are quantified using the superposition of interaction factors. Both models use the method of joining subsystems to incorporate the incident wavefield generated by the underground railway into the pile model. Results are computed for a single pile subject to an inertial loading, pile-soil-pile interactions, and a pile group subjected to excitation from an underground railway. The two models are compared in terms of accuracy, computation time, versatility and applicability, and guidelines for future vibration prediction models involving piled foundations are proposed.

Development of accurate, robust liquid equations of state for multi-phase CFD simulations with a modified AUSM+-up scheme

Eight equations of state (EOS) have been evaluated for the simulation of compressible liquid water properties, based on empirical correlations, the principle of corresponding states and thermodynamic relations. The IAPWS-IF97 EOS for water was employed as the reference case. These EOSs were coupled to a modified AUSM+-up convective flux solver to determine flow profiles for three test cases of differing flow conditions. The impact of the non-viscous interaction term discretisation scheme, interfacial pressure method and selection of low-Mach number diffusion were also compared. It was shown that a consistent discretisation scheme using the AUSM+-up solver for both the convective flux and the non-viscous interfacial term demonstrated both robustness and accuracy whilst facilitating a computationally cheaper solution than discretisation of the interfacial term independently by a central scheme. The simple empirical correlations gave excellent results in comparison to the reference IAPWS-IF97 EOS and were recommended for developmental work involving water as a cheaper and more accurate EOS than the more commonly used stiffened-gas model. The correlations based on the principles of corresponding-states and the modified Peng-Robinson cubic EOS also demonstrated a high degree of accuracy, which is promising for future work with generic fluids. Further work will encompass extension of the solver to multiple dimensions and to account for other source terms such as surface tension, along with the incorporation of phase changes. © 2013.

Modelling the vibration of tyre sidewalls

The dynamical behaviour of the sidewall has an important influence on tyre vibration characteristics. Nonetheless, it remains crudely represented in many existing models. The current work considers a geometrically accurate, two-dimensional, sidewall description, with a view to identifying potential shortcomings in the approximate formulations and identifying the physical characteristics that must be accounted for. First, the mean stress state under pressurisation and centrifugal loading is investigated. Finite-Element calculations show that, while the loaded sidewall shape remains close to a toroid, its in-plane tensions differ appreciably from the associated analytical solution. This is largely due to the inability of the anisotropic sidewall material to sustain significant azimuthal stress. An approximate analysis, based on the meridional tension alone, is therefore developed, and shown to yield accurate predictions. In conjunction with a set of formulae for the 'engineering constants' of the sidewall material, the approximate solutions provide a straightforward and efficient means of determining the base state for the vibration analysis. The latter is implemented via a 'waveguide' discretisation of a variational formulation. Its results show that, while the full geometrical description is necessary for a complete and reliable characterisation of the sidewall's vibrational properties, a one-dimensional approximation will often be satisfactory in practice. Meridional thickness variations only become important at higher frequencies (above 500 Hz for the example considered here), and rotational inertia effects appear to be minor at practical vehicle speeds. © 2013 Elsevier Ltd. All rights reserved.