Dynamic solid mechanics using finite volume methods

A procedure for evaluating the dynamic structural response of elastic solid domains is presented. A prerequisite for the analysis of dynamic fluid–structure interaction is the use of a consistent set of finite volume (FV) methods on a single unstructured mesh. This paper describes a three-dimensional (3D) FV, vertex-based method for dynamic solid mechanics. A novel Newmark predictor–corrector implicit scheme was developed to provide time accurate solutions and the scheme was evaluated on a 3D cantilever problem. By employing a small amount of viscous damping, very accurate predictions of the fundamental natural frequency were obtained with respect to both the amplitude and period of oscillation. This scheme has been implemented into the multi-physics modelling software framework, PHYSICA, for later application to full dynamic fluid structure interaction.

Dynamic mesh partitioning and load-balancing for parallel computational mechanics codes

In this Chapter we discuss the load-balancing issues arising in parallel mesh based computational mechanics codes for which the processor loading changes during the run. We briefly touch on geometric repartitioning ideas and then focus on different ways of using a graph both to solve the load-balancing problem and the optimisation problem, both locally and globally. We also briefly discuss whether repartitioning is always valid. Sample illustrative results are presented and we conclude that repartitioning is an attractive option if the load changes are not too dramatic and that there is a certain trade-off between partition quality and volume of data that the underlying application needs to migrate.

A vertex-based finite volume method applied to non-linear material problems in computational solid mechanics

A vertex-based finite volume (FV) method is presented for the computational solution of quasi-static solid mechanics problems involving material non-linearity and infinitesimal strains. The problems are analysed numerically with fully unstructured meshes that consist of a variety of two- and threedimensional element types. A detailed comparison between the vertex-based FV and the standard Galerkin FE methods is provided with regard to discretization, solution accuracy and computational efficiency. For some problem classes a direct equivalence of the two methods is demonstrated, both theoretically and numerically. However, for other problems some interesting advantages and disadvantages of the FV formulation over the Galerkin FE method are highlighted.

Analysis and modeling of heaping behavior of granular mixtures within a computational mechanics framework

This paper presents a continuum model of the flow of granular material during filling of a silo, using a viscoplastic constitutive relation based on the Drucker-Prager plasticity yield function. The performed simulations demonstrate the ability of the model to realistically represent complex features of granular flows during filling processes, such as heap formation and non-zero inclination angle of the bulk material-air interface. In addition, micro-mechanical parametrizations which account for particle size segregation are incorporated into the model. It is found that numerical predictions of segregation phenomena during filling of a binary granular mixture agree well with experimental results. Further numerical tests indicate the capability of the model to cope successfully with complex operations involving granular mixtures.

Predicting the behaviour of micro-electronic displays using computational mechanics

Micro-electronic displays are indispensible devices used in high performance applications such as aerospace, medical, marine and industrial sectors.These devices provide an interface to real time mission critical devices and therefore require good optical visual performance and high reliability, all this within varied and challenging environments.

Fracture mechanics analysis of cracks in solder joint intermetallic compounds

The trend towards miniaturization of electronic products leads to the need for very small sized solder joints. Therefore, there is a higher reliability risk that too large a fraction of solder joints will transform into Intermetallic Compounds (IMCs) at the solder interface. In this paper, fracture mechanics study of the IMC layer for SnPb and Pb-free solder joints was carried out using finite element numerical computer modelling method. It is assumed that only one crack is present in the IMC layer. Linear Elastic Fracture Mechanics (LEFM) approach is used for parametric study of the Stress Intensity Factors (SIF, KI and KII), at the predefined crack in the IMC layer of solder butt joint tensile sample. Contrary to intuition, it is revealed that a thicker IMC layer in fact increases the reliability of solder joint for a cracked IMC. Value of KI and KII are found to decrease with the location of the crack further away from the solder interfaces while other parameters are constant. Solder thickness and strain rate were also found to have a significant influence on the SIF values. It has been found that soft solder matrix generates non-uniform plastic deformation across the solder-IMC interface near the crack tip that is responsible to obtain higher KI and KII.

The Mechanics of QSR 2000

Assessment of the quality of the marine environment forms an important part of the new 1992 OSPAR Convention for the Protection of the Marine Environment of the North-East Atlantic that was ratified and entered into force on 25 March 1998. In the ministerial statement at the signing of the Convention it was agreed that the first assessment (Quality Status Report, QSR) for all Convention waters should be produced for the year 2000. To oversee this charge a new Environmental Assessment and Monitoring Committee (ASMO) was established and a junior group under this committee, to implement necessary actions, the Assessment Co-ordination Group (ACG). Because of the wide geographical diversity and varying levels of information available in different parts of the Convention area it was decided to produce five regional reports for: I The Arctic; II The North Sea; III The Celtic seas; IV The Bay of Biscay and Iberian Coast; V The Wider Atlantic, which will be synthesised in a holistic QSR for the year 2000. The report for the North Sea will largely be an update of QSR 1993 and forms the third cycle of a developing management system for the North Sea. This paper will present the procedures that have been adopted to implement the QSRs, and outlines the guidelines that have been developed for their structure, format, design and publication.

Fracture Mechanics Analysis of the Dentine-Luting Cement Interface.

The objectives of this study were to determine the fracture toughness of adhesive interfaces between dentine and clinically relevant, thin layers of dental luting cements. Cements tested included a conventional glass-ionomer, F (Fuji I), a resin-modified glass-ionomer, FP (Fuji Plus) and a compomer cement, D (DyractCem). Ten miniature short-bar chevron notch specimens were manufactured for each cement, each comprising a 40 µm thick chevron of lute, between two 1.5 mm thick blocks of bovine dentine, encased in resin composite. The interfacial KIC results (MN/m3/2) were median (range): F; 0.152 (0.14-0.16), FP; 0.306 (0.27-0.37), D; 0.351 (0.31-0.37). Non-parametric statistical analysis showed that the fracture toughness of F was significantly lower (p

An orbital period of 0.94days for the hot-Jupiter planet WASP-18b

The `hot Jupiters' that abound in lists of known extrasolar planets are thought to have formed far from their host stars, but migrate inwards through interactions with the proto-planetary disk from which they were born, or by an alternative mechanism such as planet-planet scattering. The hot Jupiters closest to their parent stars, at orbital distances of only ~0.02 astronomical units, have strong tidal interactions, and systems such as OGLE-TR-56 have been suggested as tests of tidal dissipation theory. Here we report the discovery of planet WASP-18b with an orbital period of 0.94days and a mass of ten Jupiter masses (10MJup), resulting in a tidal interaction an order of magnitude stronger than that of planet OGLE-TR-56b. Under the assumption that the tidal-dissipation parameter Q of the host star is of the order of 106, as measured for Solar System bodies and binary stars and as often applied to extrasolar planets, WASP-18b will be spiralling inwards on a timescale less than a thousandth that of the lifetime of its host star. Therefore either WASP-18 is in a rare, exceptionally short-lived state, or the tidal dissipation in this system (and possibly other hot-Jupiter systems) must be much weaker than in the Solar System.