Finite element solution of Navier Stokes equations adapted to a priori error estimates

The paper is based on qualitative properties of the solution of the Navier-Stokes equations for incompressible fluid, and on properties of their finite element solution. In problems with corner-like singularities (e.g. on the well-known L-shaped domain) usually some adaptive strategy is used. In this paper we present an alternative approach. For flow problems on domains with corner singularities we use the a priori error estimates and asymptotic expansion of the solution to derive an algorithm for refining the mesh near the corner. It gives very precise solution in a cheap way. We present some numerical results.

Moderate energy impact analysis combining phenomenological contact law with localised damage and integral equation method

A computational impact analysis methodology has been developed, based on modal analysis and a local contact force-deflection model. The contact law is based on Hertz contact theory while contact stresses are elastic, defines a modified contact theory to take account of local permanent indentation, and considers elastic recovery during unloading. The model was validated experimentally through impact testing of glass-carbon hybrid braided composite panels. Specimens were mounted in a support frame and the contact force was inferred from the deceleration of the impactor, measured by high-speed photography. A Finite Element analysis of the panel and support frame assembly was performed to compute the modal responses. The new contact model performed well in predicting the peak forces and impact durations for moderate energy impacts (15 J), where contact stresses locally exceed the linear elastic limit and damage may be deemed to have occurred. C-scan measurements revealed substantial damage for impact energies in the range of 30-50 J. For this regime the new model predictions might be improved by characterisation of the contact law hysteresis during the unloading phase, and a modification of the elastic vibration response in line with damage levels acquired during the impact. © 2011 Elsevier Ltd. All rights reserved.

A subdivision-based implementation of the hierarchical b-spline finite element method

A novel technique is presented to facilitate the implementation of hierarchical b-splines and their interfacing with conventional finite element implementations. The discrete interpretation of the two-scale relation, as common in subdivision schemes, is used to establish algebraic relations between the basis functions and their coefficients on different levels of the hierarchical b-spline basis. The subdivision projection technique introduced allows us first to compute all element matrices and vectors using a fixed number of same-level basis functions. Their subsequent multiplication with subdivision matrices projects them, during the assembly stage, to the correct levels of the hierarchical b-spline basis. The proposed technique is applied to convergence studies of linear and geometrically nonlinear problems in one, two and three space dimensions. © 2012 Elsevier B.V.

An analytical model for guided wave inspection optimization for prismatic structures of any cross section.

This paper presents an analytical modeling technique for the simulation of long-range ultrasonic guided waves in structures. The model may be used to predict the displacement field in a prismatic structure arising from any excitation arrangement and may therefore be used as a tool to design new inspection systems. It is computationally efficient and relatively simple to implement, yet gives accuracy similar to finite element analysis and semi-analytical finite element analysis methods. The model has many potential applications; one example is the optimization of part-circumferential arrays where access to the full circumference of the pipe is restricted. The model has been successfully validated by comparison with finite element solutions. Experimental validation has also been carried out using an array of piezoelectric transducer elements to measure the displacement field arising from a single transducer element in an 88.9-mm-diameter pipe. Good agreement has been obtained between the two models and the experimental data.

Finite element modelling of the sismic behaviour of water front structures

Water front structures have suffered significant damage in many of the recent earthquakes. These include gravity type quay walls, vertically composite walls, cantilever retaining walls, anchored bulkheads and similar structures. One of the primary causes for the poor performance of these classes of structures is the liquefaction of the foundation soil and in some instances liquefaction of the backfill soil. The liquefaction of the soil in-front of the quay wall tends to cause large lateral displacements and rotation of the wall. Often such gravity walls are placed on rubble mound deposited onto the sea bed.This paper presents finite element analyses of such a problem in which strength degradation of the foundation soil and the backfill material will be modelled using PZ mark III constitutive model. The performance of the wall in terms of its lateral displacement, vertical settlement and/or the rotation suffered by the wall will be presented. In addition, the contours of the horizontal and vertical effective stresses and the excess pore pressure ratio will be presented at different time instants together with hyrdraulic gradients. Immediately after the earthquake, the hydraulic gradients indicate migration of pore water into the region below the wall, suggesting further softening of the foundation soil below the wall.

3D modeling of high-T c superconductors by finite element software

A three-dimensional (3D) numerical model is proposed to solve the electromagnetic problems involving transport current and background field of a high-T c superconducting (HTS) system. The model is characterized by the E-J power law and H-formulation, and is successfully implemented using finite element software. We first discuss the model in detail, including the mesh methods, boundary conditions and computing time. To validate the 3D model, we calculate the ac loss and trapped field solution for a bulk material and compare the results with the previously verified 2D solutions and an analytical solution. We then apply our model to test some typical problems such as superconducting bulk array and twisted conductors, which cannot be tackled by the 2D models. The new 3D model could be a powerful tool for researchers and engineers to investigate problems with a greater level of complicity.

Deformation and fracture of impulsively loaded sandwich panels

Light metal sandwich panel structures with cellular cores have attracted interest for multifunctional applications which exploit their high bend strength and impact energy absorption. This concept has been explored here using a model 6061-T6 aluminum alloy system fabricated by friction stir weld joining extruded sandwich panels with a triangular corrugated core. Micro-hardness and miniature tensile coupon testing revealed that friction stir welding reduced the strength and ductility in the welds and a narrow heat affected zone on either side of the weld by approximately 30%. Square, edge clamped sandwich panels and solid plates of equal mass per unit area were subjected to localized impulsive loading by the impact of explosively accelerated, water saturated, sand shells. The hydrodynamic load and impulse applied by the sand were gradually increased by reducing the stand-off distance between the test charge and panel surfaces. The sandwich panels suffered global bending and stretching, and localized core crushing. As the pressure applied by the sand increased, face sheet fracture by a combination of tensile stretching and shear-off occurred first at the two clamped edges of the panels that were parallel with the corrugation and weld direction. The plane of these fractures always lay within the heat affected zone of the longitudinal welds. For the most intensively loaded panels additional cracks occurred at the other clamped boundaries and in the center of the panel. To investigate the dynamic deformation and fracture processes, a particle-based method has been used to simulate the impulsive loading of the panels. This has been combined with a finite element analysis utilizing a modified Johnson-Cook constitutive relation and a Cockcroft-Latham fracture criterion that accounted for local variation in material properties. The fully coupled simulation approach enabled the relationships between the soil-explosive test charge design, panel geometry, spatially varying material properties and the panel's deformation and dynamic failure responses to be explored. This comprehensive study reveals the existence of a strong instability in the loading that results from changes in sand particle reflection during dynamic evolution of the panel's surface topology. Significant fluid-structure interaction effects are also discovered at the sample sides and corners due to changes of the sand reflection angle by the edge clamping system. © 2012 Elsevier Ltd. All rights reserved.

A finite element model of magnetization of superconducting bulks using a solid-state flux pump

Superconductors have a bright future; they are able to carry very high current densities, switch rapidly in electronic circuits, detect extremely small perturbations in magnetic fields, and sustain very high magnetic fields. Of most interest to large-scale electrical engineering applications are the ability to carry large currents and to provide large magnetic fields. There are many projects that use the first property, and these have concentrated on power generation, transmission, and utilization; however, there are relatively few, which are currently exploiting the ability to sustain high magnetic fields. The main reason for this is that high field wound magnets can and have been made from both BSCCO and YBCO, but currently, their cost is much higher than the alternative provided by low-Tc materials such as Nb3Sn and NbTi. An alternative form of the material is the bulk form, which can be magnetized to high fields. This paper explains the mechanism, which allows superconductors to be magnetized without the need for high field magnets to perform magnetization. A finite-element model is presented, which is based on the E-J current law. Results from this model show how magnetization of the superconductor builds up cycle upon cycle when a traveling magnetic wave is induced above the superconductor. © 2011 IEEE.