The post-processed Galerkin method applied to non-linear shell vibrations

We present the results of a computational study of the post-processed Galerkin methods put forward by Garcia-Archilla et al. applied to the non-linear von Karman equations governing the dynamic response of a thin cylindrical panel periodically forced by a transverse point load. We spatially discretize the shell using finite differences to produce a large system of ordinary differential equations (ODEs). By analogy with spectral non-linear Galerkin methods we split this large system into a 'slowly' contracting subsystem and a 'quickly' contracting subsystem. We then compare the accuracy and efficiency of (i) ignoring the dynamics of the 'quick' system (analogous to a traditional spectral Galerkin truncation and sometimes referred to as 'subspace dynamics' in the finite element community when applied to numerical eigenvectors), (ii) slaving the dynamics of the quick system to the slow system during numerical integration (analogous to a non-linear Galerkin method), and (iii) ignoring the influence of the dynamics of the quick system on the evolution of the slow system until we require some output, when we 'lift' the variables from the slow system to the quick using the same slaving rule as in (ii). This corresponds to the post-processing of Garcia-Archilla et al. We find that method (iii) produces essentially the same accuracy as method (ii) but requires only the computational power of method (i) and is thus more efficient than either. In contrast with spectral methods, this type of finite-difference technique can be applied to irregularly shaped domains. We feel that post-processing of this form is a valuable method that can be implemented in computational schemes for a wide variety of partial differential equations (PDEs) of practical importance.

An Arbitrary Lagrangian Eulerian Method for Three-Phase Flows with Triple Junction Points

We present a moving mesh method suitable for solving two-dimensional and axisymmetric three-liquid flows with triple junction points. This method employs a body-fitted unstructured mesh where the interfaces between liquids are lines of the mesh system, and the triple junction points (if exist) are mesh nodes. To enhance the accuracy and the efficiency of the method, the mesh is constantly adapted to the evolution of the interfaces by refining and coarsening the mesh locally; dynamic boundary conditions on interfaces, in particular the triple points, are therefore incorporated naturally and accurately in a Finite- Element formulation. In order to allow pressure discontinuity across interfaces, double-values of pressure are necessary for interface nodes and triple-values of pressure on triple junction points. The resulting non-linear system of mass and momentum conservation is then solved by an Uzawa method, with the zero resultant condition on triple points reinforced at each time step. The method is used to investigate the rising of a liquid drop with an attached bubble in a lighter liquid.

A coupled hydro-mechanical analysis for prediction of hydraulic fracture propagation in saturated porous media using EFG mesh-less method

The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system of equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on the penalty method. In order to model discrete fractures, the so-called diffraction method is used.Examples are presented and the results are compared to some closed-form solutions and FEM approximations in order to demonstrate the validity of the developed model and its capabilities. The model is able to take the anisotropy and inhomogeneity of the material into account. The applicability of the model is examined by simulating hydraulic fracture initiation and propagation process from a borehole by injection of fluid. The maximum tensile strength criterion and Mohr-Coulomb shear criterion are used for modelling tensile and shear fracture, respectively. The model successfully simulates the leak-off of fluid from the fracture into the surrounding material. The results indicate the importance of pore fluid pressure in the initiation and propagation pattern of fracture in saturated soils. © 2013 Elsevier Ltd.

HEAVY VEHICLE VIBRATION - A CASE STUDY.

A method is presented for the digital simulation of multiple degrees-of-freedom lumped parameter vibrating systems with arbitrary constitutive elements in an inertial frame of reference. The geometry of the system is treated independently of the constitutive elements and as a result nonlinear (time domain) or linearised (frequency domain) calculations may be performed using a single input description. The method is used to simulate a 3-axle rigid heavy commercial vehicle for harsh vibrating conditions. Some of the assumptions to which the calculations are sensitive are examined. Agreement between the response of a 3-dimensional whole vehicle model and measurements on the test vehicle is satisfactory.

A comparison of methods for the coupled analysis of floating structures

The dynamic analysis of a deepwater floating platform and the associated mooring/riser system should ideally be fully coupled to ensure a reliable response prediction. It is generally held that a time domain analysis is the only means of capturing the various coupling and nonlinear effects accurately. However, in recent work it has been found that for an ultra-deepwater floating system (2000m water depth), the highly efficient frequency domain approach can provide highly accurate response predictions. One reason for this is the accuracy of the drag linearization procedure over both first and second order motions, another reason is the minimal geometric nonlinearity displayed by the mooring lines in deepwater. In this paper, the aim is to develop an efficient analysis method for intermediate water depths, where both mooring/vessel coupling and geometric nonlinearity are of importance. It is found that the standard frequency domain approach is not so accurate for this case and two alternative methods are investigated. In the first, an enhanced frequency domain approach is adopted, in which line nonlinearities are linearized in a systematic way. In the second, a hybrid approach is adopted in which the low frequency motion is solved in the time domain while the high frequency motion is solved in the frequency domain; the two analyses are coupled by the fact that (i) the low frequency motion affects the mooring line geometry for the high frequency motion, and (ii) the high frequency motion affects the drag forces which damp the low frequency motion. The accuracy and efficiency of each of the methods are systematically compared. Copyright © 2007 by ASME.

Monitoring tunnel deformation induced by close-proximity bored tunneling using distributed optical fiber strain measurements

The paper describes the use of optical fiber Brillouin Optical Time Domain Reflectometry (BOTDR) to monitor the strain distribution in an existing tunnel while a twin tunnel was bored at close-proximity. The twin circular bored tunnels between Serangoon and Bartley stations on the new Circle Line Stage 3 subway in Singapore were constructed at close-proximity to avoid underpinning the foundations of adjacent buildings. The minimum clear separation of the two tunnels is 2.3m (0.4 times the tunnel diameter). The Outer Tunnel was constructed first, followed by the Inner Tunnel, with the earth-pressure balance tunnel boring machines maintained at a minimum of 100m apart. In this trial application of BOTDR, the strain distribution along the Outer Tunnel was measured, in order to monitor its deformation due to the boring of the Inner Tunnel at close-proximity. The aim of the trial application was to determine the practicality of this monitoring method for future use in 'live' tunnels. This paper compares the measurements obtained from optical fiber BOTDR with conventional methods of tunnel monitoring and describes preliminary installation and workmanship guidelines derived from lessons learnt during this trial. © 2007 ASCE.

Fibre optics monitoring of clay cuttings and embankments along London's ring motorway

Although a wide range of techniques exist for slope monitoring, the task of monitoring slopes is sometimes complicated by the extensive nature and unpredictability of slope movements. The Brillouin optical time-domain reflectometer (BOTDR) is a distributed optical fiber strain measurement technology utilising Brillouin scattering. This method measures continuous strain along a standard optical fibre over a distance up to 10 km and hence has potential to detect deformations and diagnose problems along large sections of slopes and embankments. This paper reports the demonstration of BOTDR method for monitoring surface ground movements of clay cuttings and embankments along London's ring M25 motorway. A field trial investigating varying methods of onsite fibre optic installations was conducted. The surrounding ground was artificially moved by excavating a 3 m deep trench perpendicular to the instrumented sections. Results obtained from onsite installations after slope movement demonstrate a half-pipe covered fibre optic installed on wide (200mm) Tensar ™SS20 geogrid gives the most consistent recorded strain change profile. Initial conclusions suggest this method best represents induced ground motion at the surface and hence is recommended for implementation in future sitework. Copyright ASCE 2008.

Efficient minimum manoeuvre time optimisation of an oversteering vehicle at constant forward speed

A receding horizon steering controller is presented, capable of pushing an oversteering nonlinear vehicle model to its handling limit while travelling at constant forward speed. The controller is able to optimise the vehicle path, using a computationally efficient and robust technique, so that the vehicle progression along a track is maximised as a function of time. The resultant method forms part of the solution to the motor racing objective of minimising lap time. © 2011 AACC American Automatic Control Council.

The use of spatial impulse responses to characterise flexible forming processes with mobile tools

A novel test method for the characterisation of flexible forming processes is proposed and applied to four flexible forming processes: Incremental Sheet Forming (ISF), conventional spinning, the English wheel and power hammer. The proposed method is developed in analogy with time-domain control engineering, where a system is characterised by its impulse response. The spatial impulse response is used to characterise the change in workpiece deformation created by a process, but has also been applied with a strain spectrogram, as a novel way to characterise a process and the physical effect it has on the workpiece. Physical and numerical trials to study the effects of process and material parameters on spatial impulse response lead to three main conclusions. Incremental sheet forming is particularly sensitive to process parameters. The English wheel and power hammer are strongly similar and largely insensitive to both process and material parameters. Spinning develops in two stages and is sensitive to most process parameters, but insensitive to prior deformation. Finally, the proposed method could be applied to modelling, classification of existing and novel processes, product-process matching and closed-loop control of flexible forming processes. © 2012 Elsevier B.V.

Prediction of tread block forces for a free-rolling tyre in contact with a rough road

This paper describes a new approach to model the forces on a tread block for a free-rolling tyre in contact with a rough road. A theoretical analysis based on realistic tread mechanical properties and road roughness is presented, indicating partial contact between a tread block and a rough road. Hence an asperity-scale indentation model is developed using a semi-empirical formulation, taking into account both the rubber viscoelasticity and the tread block geometry. The model aims to capture the essential details of the contact at the simplest level, to make it suitable as part of a time-domain dynamic analysis of the coupled tyre-road system. The indentation model is found to have a good correlation with the finite element (FE) predictions and is validated against experimental results using a rolling contact rig. When coupled to a deformed tyre belt profile, the indentation model predicts normal and tangential force histories inside the tyre contact patch that show good agreement with FE predictions. © 2012 Elsevier B.V..

Robust F0 estimation based on a multi-microphone periodicity function for distant-talking speech

This work addresses the problem of deriving F0 from distanttalking speech signals acquired by a microphone network. The method here proposed exploits the redundancy across the channels by jointly processing the different signals. To this purpose, a multi-microphone periodicity function is derived from the magnitude spectrum of all the channels. This function allows to estimate F0 reliably, even under reverberant conditions, without the need of any post-processing or smoothing technique. Experiments, conducted on real data, showed that the proposed frequency-domain algorithm is more suitable than other time-domain based ones.

Nonlinear phenomena in thermoacoustics: A comparison between single-mode and multi-mode methods

Nonlinear analysis of thermoacoustic instability is essential for prediction of frequencies and amplitudes of limit cycles. In frequency domain analyses, a quasi-linear transfer function between acoustic velocity and heat release rate perturbations, called the flame describing function (FDF), is obtained from a flame model or experiments. The FDF is a function of the frequency and amplitude of velocity perturbations but only contains the heat release response at the forcing frequency. While the gain and phase of the FDF provide insight into the nonlinear dynamics of the system, the accuracy of its predictions remains to be verified for different types of nonlinearity. In time domain analyses, the governing equations of the fully coupled problem are solved to find the time evolution of the system. One method is to discretize the governing equations using a suitable basis, such as the natural acoustic modes of the system. The number of modes used in the discretization alters the accuracy of the solution. In our previous work we have shown that predictions using the FDF are almost exactly the same as those obtained from the time-domain using only one mode for the discretization. We call this the single-mode method. In this paper we compare results from the single-mode and multi-mode methods, applied to a thermoacoustic system of a premixed flame in a tube. For some cases, the results differ greatly in both amplitude as well as frequency content. This study shows that the contribution from higher and subharmonics to the nonlinear dynamics can be significant and must be considered for an accurate and comprehensive analysis of thermoacoustic systems. Hence multi-mode simulations are necessary, and the single-mode method or the FDF may be insufficient to capture some of the complex nonlinear behaviour in fhermoacoustics.

Shock dynamics of random structures.

Predicting the response of a structure following an impact is of interest in situations where parts of a complex assembly may come into contact. Standard approaches are based on the knowledge of the impulse response function, requiring the knowledge of the modes and the natural frequencies of the structure. In real engineering structures the statistics of higher natural frequencies follows those of the Gaussian Orthogonal Ensemble, this allows the application of random point process theory to get a mean impulse response function by the knowledge of the modal density of the structure. An ensemble averaged time history for both the response and the impact force can be predicted. Once the impact characteristics are known in the time domain, a simple Fourier Transform allows the frequency range of the impact excitation to be calculated. Experimental and numerical results for beams, plates, and cylinders are presented to confirm the validity of the method.

On the time delay in continuous in-cylinder sampling from IC engines

When gas sample is continuously drawn from the cylinder of an internal combustion engine, the sample that appears at the end of the sampling system corresponds to the in-cylinder content sometime ago because of the finite transit time which is a function of the cylinder pressure history. This variable delay causes a dispersion of the sample signal and makes the interpretation of the signal difficult An unsteady flow analysis of a typical sampling system was carried out for selected engine loads and speeds. For typical engine operation, a window in which the delay is approximately constant may be found. This window gets smaller with increase in engine speed, with decrease in load, and with the increase in exit pressure of the sampling system.