Theoretical simulation studies of pulsed field magnetisation of (RE)BCO bulk superconductors

The magnetisation of bulk high temperature superconductors (HTS), such as RE-Ba-Cu-O [(RE)BCO, where RE is a rare earth element or Y], by a practical technique is essential for their application in high field, permanent magnet-like devices. Research to-date into the pulsed field magnetisation (PFM) of these materials, however, has been limited generally to experimental techniques, with relatively little progress in the development of theoretical models. This is because not only is a multi-physics approach needed to take account of the heating of the samples but also the high electric fields generated are well above the regime in which there are reliable experimental results. This paper describes a framework of theoretical simulation using the finite element method (FEM) that is applicable to both single- and multi-pulse magnetisation processes of (RE)BCO bulk superconductors. The model incorporates the heat equation and provides a convenient way of determining the distribution of trapped field, current density and temperature change within a bulk superconductor at each stage of the magnetisation process. An example of the single-pulse magnetisation of a (RE)BCO bulk is described. Potentially, the model may serve as a cost-effective tool for the optimisation of the bulk geometry and the magnetisation profile in multi-pulse magnetisation processes. © 2010 IOP Publishing Ltd.

Elastodynamic erosion of thermal barrier coatings

The finite element method is used to analyze the elastodynamic response of a columnar thermal barrier coating due to normal impact and oblique impact by an erosive particle. An assessment is made of the erosion by crack growth from preexisting flaws at the edge of each column: it is demonstrated that particle impacts can be sufficiently severe to give rise to columnar cracking. First, the transient stress state induced by the normal impact of a circular cylinder or a sphere is calculated in order to assess whether a 2D calculation adequately captures the more realistic 3D behavior. It is found that the transient stress states for the plane strain and axisymmetric models are similar. The sensitivity of response to particle diameter and to impact velocity is determined for both the cylinder and the sphere. Second, the transient stress state is explored for 2D oblique impact by a circular cylindrical particle and by an angular cylindrical particle. The sensitivity of transient tensile stress within the columns to particle shape (circular and angular), impact angle, impact location, orientation of the angular particle, and to the level of friction is explored in turn. The paper concludes with an evaluation of the effect of inclining the thermal barrier coating columns upon their erosion resistance. © 2011 The American Ceramic Society.

The calculation of acoustic shielding of engine noise by the Silent Aircraft Airframe

The Silent Aircraft airframe has a flying wing design with a large wing planform and a propulsion system embedded in the rear of the airframe with intake on the upper surface of the wing. In the present paper, boundary element calculations are presented to evaluate acoustic shielding at low frequencies. Besides the three-dimensional geometry of the Silent Aircraft airframe, a few two-dimensional problems are considered that provide some physical insight into the shielding calculations. Mean flow refraction effects due to forward flight motion are accounted for by a simple time transformation that decouples the mean-flow and acoustic-field calculations. It is shown that significant amount of shielding can be obtained in the shadow region where there is no direct line of sight between the source and observer. The boundary element solutions are restricted to low frequencies. We have used a simple physically-based model to extend the solution to higher frequencies. Based on this model, using a monopole acoustic source, we predict at least an 18 dBA reduction in the overall sound pressure level of forward-propagating fan noise because of shielding.

Modelling of equipment loaded panels for robust control investigations

This paper develops a modelling technique for equipment load panels which directly produces (adequate) models of the underlying dynamics on which to base robust controller design/evaluations. This technique is based on the use of the Lagrange's equations of motion and the resulting models are verified against those produced by a finite Element Method Model.

A cohesive approach to thin-shell fracture and fragmentation

We develop a finite-element method for the simulation of dynamic fracture and fragmentation of thin-shells. The shell is spatially discretized with subdivision shell elements and the fracture along the element edges is modeled with a cohesive law. In order to follow the propagation and branching of cracks, subdivision shell elements are pre-fractured ab initio and the crack opening is constrained prior to crack nucleation. This approach allows for shell fracture in an in-plane tearing mode, a shearing mode, or a bending of hinge mode. The good performance of the method is demonstrated through the simulation of petalling failure experiments in aluminum plates. © 2005 Elsevier B.V. All rights reserved.

The elastic-plastic indentation response of a columnar thermal barrier coating

Thermal barrier coatings with a columnar microstructure are prone to erosion damage by a mechanism of surface cracking upon impact by small foreign particles. In order to explore this erosion mechanism, the elastic indentation and the elastic-plastic indentation responses of a columnar thermal barrier coating to a spherical indenter were determined by the finite element method and by analytical models. It was shown that the indentation response is intermediate between that of a homogeneous half-space and that given by an elastic-plastic mattress model (with the columns behaving as independent non-linear springs). The sensitivity of the indentation behaviour to geometry and to the material parameters was explored: the diameter of the columns, the gap width between columns, the coefficient of Coulomb friction between columns and the layer height of the thermal barrier coating. The calculations revealed that the level of induced tensile stress is sufficient to lead to cracking of the columns at a depth of about the column radius. It was also demonstrated that the underlying soft bond coat can undergo plastic indentation when the coating comprises parallel columns, but this is less likely for the more realistic case of a random arrangement of tapered columns. © 2009 Elsevier B.V.

Mechanical constraint and loading on ferroelectric memory capacitors

The influence of mechanical constraint imposed by device geometry upon the switching response of a ferroelectric thin film memory capacitor is investigated. The memory capacitor was represented by two-dimensional ferroelectric islands of different aspect ratio, mechanically constrained by surrounding materials. Its ferroelectric non-linear behaviour was modeled by a crystal plasticity constitutive law and calculated using the finite element method. The switching response of the device, in terms of remnant charge storage, was determined as a function of geometry and constraint. The switching response under applied in-plane tensile stress and hydrostatic pressure was also studied experimentally. Our results showed that (1) the capacitor's aspect ratio could significantly affect the clamping behaviour and thus the remnant polarization, (2) it was possible to maximise the switching charge through the optimisation of the device geometry, and (3) it is possible to find a critical switching stress at zero electric field and a critical coercive field at zero residual stress. © 2009 Materials Research Society.

Effects of soil fabric on behaviors of granular soils: Microscopic modeling

The effects of initial soil fabric on behaviors of granular soils are investigated by using Distinct Element Method (DEM) numerical simulation. Soil specimens are represented by an assembly of non-uniform sized spheres with different initial contact normal distributions. Isotropically consolidated triaxial compression loading and extension unloading in both undrained and drained conditions are simulated for vertically- and horizontally-sheared specimens. The numerical simulation results are compared qualitatively with the published experimental data and the effects of initial soil fabric on resulting soil behaviors are discussed, including the effects of specimen reconstitution methods, effects of large preshearing, and anisotropic characteristics in undrained and drained conditions. The effects of initial soil fabric and mode of shearing on the quasi-steady state line are also investigated. The numerical simulation results can systematically explain that the observed experimental behaviors of granular soils are due principally to their conditions of the initial soil fabric. This outcome provides insights into the observed phenomena in microscopic view. © 2011 Elsevier Ltd.

Vibration of fluid loaded conical shells.

An analytical model is presented to describe the vibration of a truncated conical shell with fluid loading in the low frequency range. The solution for the dynamic response of the shell is presented in the form of a power series. Fluid loading is taken into account by dividing the shell into narrow strips which are considered to be locally cylindrical. Analytical results are presented for different boundary conditions and have been compared with the computational results from a boundary element model. Limitations of the model to the low frequency range are discussed.

Collapse mechanism maps for the hollow pyramidal core of a sandwich panel under transverse shear

The finite element method has been used to develop collapse mechanism maps for the shear response of sandwich panels with a stainless steel core comprising hollow struts. The core topology comprises either vertical tubes or inclined tubes in a pyramidal arrangement. The dependence of the elastic and plastic buckling modes upon core geometry is determined, and optimal geometric designs are obtained as a function of core density. For the hollow pyramidal core, strength depends primarily upon the relative density ρ̄ of the core with a weak dependence upon tube slenderness. At ρ̄ below about 3%, the tubes of the pyramidal core buckle plastically and the peak shear strength scales linearly with ρ̄. In contrast, at ρ̄ above 3%, the tubes do not buckle and a stable shear response is observed. The predictions of the current study are in excellent agreement with previous measurements on the shear strength of the hollow pyramidal core, and suggest that this core topology is attractive from the perspectives of both core strength and energy absorption. © 2011 Elsevier Ltd. All rights reserved.

Design methodology for composite structures: A small low air-speed wind turbine blade case study

A small low air-speed wind turbine blade case study is used to demonstrate the effectiveness of a materials and design selection methodology described by Monroy Aceves et al. (2008) [24] for composite structures. The blade structure comprises a shell of uniform thickness and a unidirectional reinforcement. The shell outer geometry is fixed by aerodynamic considerations. A wide range of lay-ups are considered for the shell and reinforcement. Structural analysis is undertaken using the finite element method. Results are incorporated into a database for analysis using material selection software. A graphical selection stage is used to identify the lightest blade meeting appropriate design constraints. The proposed solution satisfies the design requirements and improves on the prototype benchmark by reducing the mass by almost 50%. The flexibility of the selection software in allowing identification of trends in the results and modifications to the selection criteria is demonstrated. Introducing a safety factor of two on the material failure stresses increases the mass by only 11%. The case study demonstrates that the proposed design methodology is useful in preliminary design where a very wide range of cases should be considered using relatively simple analysis. © 2011 Elsevier Ltd.

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.

Dynamic characterisation and longitudinal strength of swaged joints

Swaging is a cold working process involving plastic deformation of the work piece to change its shape. A swaged joint is a connection between two components whereby a swaging tool induces plastic deformation of the components at their junction to effectively bind them together. This is commonly used when welding or other standard joining techniques are not viable. Swaged joints can be found for example, in nuclear fuel assemblies to connect the edges of thin rectangular plates to a supporting structure or frame. The aim of this work is to find a model to describe the vibrational behaviour of a swaged joint and to estimate its strength in resisting a longitudinally applied load. The finite element method and various experimental rigs were used in order to find relationships between the natural frequencies of the plate, the joint stiffness and the force required to shift the plate against the restraining action of the swage connection. It is found that a swaged joint is dynamically equivalent to a simple support with the rotation elastically restrained and a small stiffness is enough to resist an important load. © 2011 Elsevier Ltd. All rights reserved.

Effects of initial soil fabric and mode of shearing on quasi-steady state line for monotonic undrained behavior

The effects of initial soil fabric and mode of shearing on quasi-steady state line in void ratiostress space are studied by employing the Distinct Element Method numerical analysis. The results show that the initial soil fabric and the mode of shearing have a profound effect on the location of the quasi-steady state line. The evolution of the soil fabric during the course of undrained shearing shows that the specimens with different initial soil fabrics reach quasi-steady state at various soil fabric conditions. At quasi-steady state, the soil fabric has a significant adjustment to change its behavior from contractive to dilative. As the stress state approaches the steady state, the soil fabrics of different initial conditions become similar. The numerical analysis results are compared qualitatively with the published experimental data and the effects of specimen reconstitution methods and mode of shearing found in the experimental studies canbe systematically explained by the numerical analysis. © 2009 Taylor & Francis Group.

Fabrication of carbon nanotubes on inter-digitated metal electrode for switchable nanophotonic devices

This paper reports the modeling and characterization of interdigitated rows of carbon nanotube electrodes used to address a liquid crystal media. Finite Element Method modeling of the nanotube arrays was performed to analyze the static electric Fields produced to Find suitable electrode geometry. A device was fabricated based on the simulation results and electro optics characteristics of the device are presented. This Finding has applications in the development of micron and submicron pixels, precise beem steering and nanotube based active back planes.