A finite-element analysis of viscoelastically damped sandwich plates

A finite element analysis associated with an asymptotic solution method for the harmonic flexural vibration of viscoelastically damped unsymmetrical sandwich plates is given. The element formulation is based on generalization of the discrete Kirchhoff theory (DKT) element formulation. The results obtained with the first order approximation of the asymptotic solution presented here are the same as those obtained by means of the modal strain energy (MSE) method. By taking more terms of the asymptotic solution, with successive calculations and use of the Padé approximants method, accuracy can be improved. The finite element computation has been verified by comparison with an analytical exact solution for rectangular plates with simply supported edges. Results for the same plates with clamped edges are also presented.

A Finite Element Analysis of a Crack Penetrating or Deflecting into an Interface in a Thin Laminate

The crack tip driving force of a crack growing from a pre-crack that is perpendicular to and terminating at an interface between two materials is investigated using a linear fracture mechanics theory. The analysis is performed both for a crack penetrating the interface, growing straight ahead, and for a crack deflecting into the interface. The results from finite element calculations are compared with asymptotic solutions for infinitesimally small crack extensions. The solution is found to be accurate even for fairly large amounts of crack growth. Further, by comparing the crack tip driving force of the deflected crack with that of the penetrating crack, it is shown how to control the path of the crack by choosing the adhesion of the interface relative to the material toughness.

Finite element analysis on a circular wedge prism with 400 mm diameter

The paper comprehensively analyzes the distortions of a circular wedge prism with 400 mm diameter in a scanner by method of optical-mechanical-thermal integrating analysis. The structure and intensity of the prism assembly is verified and checked, and the surface deformations of the prism under gravity load, as well as the thermo-elastic distortions of the prism, are analyzed in detail and evaluated, which is finally contrasted with the measured values of Zygo Mark interferometer. The results show: the maximal distortion of the prism assembly is 10 nm magnitude and the maximal stress is 0.441 Mpa, which has much tolerance to the precision requirement of structure and the admissible stress of material; the influence of heat effect on the surface deformations of prism is proved to be far greater than the influence of gravity load, so some strict temperature-controlled measures are to be considered when the scanner is used. (c) 2006 Elsevier GmbH. All rights reserved.

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.

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.

Stress-strain response of hydrate-bearing sands: Numerical study using discrete element method simulations

Gas hydrate is a crystalline solid found within marine and subpermafrost sediments. While the presence of hydrates can have a profound effect on sediment properties, the stress-strain behavior of hydrate-bearing sediments is poorly understood due to inherent limitations in laboratory testing. In this study, we use numerical simulations to improve our understanding of the mechanical behavior of hydrate-bearing sands. The hydrate mass is simulated as either small randomly distributed bonded grains or as "ripened hydrate" forming patchy saturation, whereby sediment clusters with 100% pore-filled hydrate saturation are distributed within a hydrate-free sediment. Simulation results reveal that reduced sand porosity and higher hydrate saturation cause an increase in stiffness, strength, and dilative tendency, and the critical state line shifts toward higher void ratio and higher shear strength. In particular, the critical state friction angle increases in sands with patchy saturation, while the apparent cohesion is affected the most when the hydrate mass is distributed in pores. Sediments with patchy hydrate distribution exhibit a slightly lower strength than sediments with randomly distributed hydrate. Finally, hydrate dissociation under drained conditions leads to volume contraction and/or stress relaxation, and pronounced shear strains can develop if the hydrate-bearing sand is subjected to deviatoric loading during dissociation.

Implementation of a PI phase angle controller for finite element analysis of the BDFM

Time-stepping finite element analysis of the BDFM for a specific load condition is shown to be a challenging problem because the excitation required cannot be predetermined and the BDFM is not open loops stable for all operating conditions. A simulation approach using feedback control to set the torque and stabilise the BDFM is presented together with implementation details. The performance of the simulation approach is demonstrated with an example and computed results are compared with measurements.

Linear multiscale analysis and finite element validation of stretching and bending dominated lattice materials

The paper presents a multiscale procedure for the linear analysis of components made of lattice materials. The method allows the analysis of both pin-jointed and rigid-jointed microtruss materials with arbitrary topology of the unit cell. At the macroscopic level, the procedure enables to determine the lattice stiffness, while at the microscopic level the internal forces in the lattice elements are expressed in terms of the macroscopic strain applied to the lattice component. A numeric validation of the method is described. The procedure is completely automated and can be easily used within an optimization framework to find the optimal geometric parameters of a given lattice material. © 2011 Elsevier Ltd. All rights reserved.

A fixed-grid b-spline finite element technique for fluid-structure interaction

We present a fixed-grid finite element technique for fluid-structure interaction problems involving incompressible viscous flows and thin structures. The flow equations are discretised with isoparametric b-spline basis functions defined on a logically Cartesian grid. In addition, the previously proposed subdivision-stabilisation technique is used to ensure inf-sup stability. The beam equations are discretised with b-splines and the shell equations with subdivision basis functions, both leading to a rotation-free formulation. The interface conditions between the fluid and the structure are enforced with the Nitsche technique. The resulting coupled system of equations is solved with a Dirichlet-Robin partitioning scheme, and the fluid equations are solved with a pressure-correction method. Auxiliary techniques employed for improving numerical robustness include the level-set based implicit representation of the structure interface on the fluid grid, a cut-cell integration algorithm based on marching tetrahedra and the conservative data transfer between the fluid and structure discretisations. A number of verification and validation examples, primarily motivated by animal locomotion in air or water, demonstrate the robustness and efficiency of our approach. © 2013 John Wiley & Sons, Ltd.

An important pitfall of pseudo-static finite element analysis

Finite Element (FE) pseudo-static analysis can provide a good compromise between simplified methods of dynamic analysis and time domain analysis. The pseudo-static FE approach can accurately model the in situ, stresses prior to seismic loading (when it follows a static analysis simulating the construction sequence) is relatively simple and not as computationally expensive as the time domain approach. However this method should be used with caution as the results can be sensitive to the choice of the mesh dimensions. In this paper two simple examples of pseudo-static finite element analysis are examined parametrically, a homogeneous slope and a cantilever retaining wall, exploring the sensitivity of the pseudo-static analysis results on the adopted mesh size. The mesh dependence was found to be more pronounced for problems with high critical seismic coefficients values (e.g. gentle slopes or small walls), as in these cases a generalised layer failure mechanism is developed simultaneously with the slope or wall mechanism. In general the mesh width was found not to affect notably the predicted value of critical seismic coefficient but to have a major impact on the predicted movements. © 2012 Elsevier Ltd.

Enhancements of the simulation method on the edge effect in resin transfer molding processes

In resin transfer molding processes, small clearances exist between the fiber preform and the mold edges, which result in a preferential resin flow in the edge channel and then disrupt the flow patterns during the mold filling stage. A mathematical model including the effect of cavity thickness on resin flow was developed for flow behavior involving the interface between an edge channel and a porous medium. According to the mathematical analysis of momentum equations in a fully developed rectangular duct and formulations of the equivalent edge permeability, comparing with three-dimensional Navier-Stokes equations, the governing equations were modified in the edge channel. The volume of fluid (VOF) method was applied to track the flow front. A simple case is numerically simulated using the modified governing equations. The effects of edge channel width and cavity thickness on flow front and inlet pressure are analyzed, and the evolution characteristics of simulated results are in agreement with the experimental results. (c) 2007 Elsevier B.V. All rights reserved

Evaluation of distortions in laser welded shipbuilding parts using local-global finite element approach

In recognition of the differences of scale between the welding pool and the heat affected zone along the welding line on one hand, and the overall size of the components being welded on the other, a local-global finite element approach was developed for the evaluation of distortions in laser welded shipbuilding parts. The approach involves the tandem use of a 'local' and a 'global' step. The local step involves a three-dimensional finite element model for the simulation of the laser welding process using the Sysweld finite element code, which takes into account thermal, metallurgical, and mechanical aspects. The simulation of the laser welding process was performed using a non-linear heat transfer analysis, based on a keyhole formation model, and a coupled transient thermomechanical analysis, which takes into account metallurgical transformations using the temperature dependent material properties and the continuous cooling transformation diagram. The size and shape of the keyhole used in the local finite element analysis was evaluated using a keyhole formation model and the Physica finite volume code. The global step involves the transfer of residual plastic strains and the stiffness of the weld obtained from the local model to the global analysis, which then provides the predicted distortions for the whole part. This newly developed methodology was applied to the evaluation of global distortions due to laser welding of stiffeners on a shipbuilding part. The approach has been proved reliable in comparison with experiments and of practical industrial use in terms of computing time and storage.

Simple control of the polarisation in uniform hybrid waveguide-planar leaky-wave antennas

A simple and original mechanism to control the polarisation of uniform hybrid waveguide-planar leaky-wave antennas is proposed. The operation is based on introducing simple modifications of the planar dimensions of the structure cross-section, which is shown to control the horizontal and vertical components of the radiated fields. The proposed antenna dispenses with the need for periodic elements, commonly used in flexible polarised leaky-wave antennas, and therefore significantly reduces the design complexity. Parametric curves have been obtained to assist in the simple and efficient design of the proposed antenna. The novel mechanism is illustrated by means of several antenna prototypes operating at 5.7 GHz, producing linear, elliptical and circular polarisations. Commercial three-dimensional Finite Element Method has been used for the simulations, and the results are validated with experimental testing.[br].

Evaluation of cement stresses in finite element analyses of cemented orthopaedic implants

Stress analysis of the cement fixation of orthopaedic implants to bone is frequently? carried out using finite element analysis. However, the stress distribution in the cement laver is usually intricate, and it is difficult to report it in a way that facilitates comparison of implants for pre-clinical testing. To study this problem, and make recommendations for stress reporting, a finite element analysis of a hip prosthesis implanted into a synthetic composite femur is developed. Three cases are analyzed: a fully bonded implant, a debonded implant, and a debonded implant where the cement is removed distal to the stein tip. In addition to peak stresses, and contour and vector plots, a stressed volume and probability-of-failure analysis is reported. It is predicted that the peak stress is highest for the debonded stem, and that removal of the distal cement more than halves this peak stress. This would suggest that omission of the distal cement is good for polished prostheses (as practiced for the Exeter design). However; if the percentage of cement stressed above a certain threshold (say 3 MPa) is considered, then the removal of distal cement is shown to be disadvantageous because a higher volume of cement is stressed to above the threshold. Vector plots clearly demonstrate the different load transfer for bonded and debonded prostheses: A bonded stein generates maximum tensile stresses in the longitudinal direction, whereas a debonded stem generates most tensile stresses in the hoop direction, except near the tip where tensile longitudinal stresses occur due to subsidence of the stein. Removal of the cement distal to the tip allows greater subsidence but alleviates these large stresses at the tip, albeit at the expense of increased hoop stresses throughout the mantle. It is concluded that a thorough analysis of cemented implants should not report peak stress, which can be misleading, but rather stressed volume, and that vector plots should be reported if a precise analysis of the load transfer mechanism is required.