A class of domain decomposition preconditioners for hp-discontinuous Galerkin finite element methods

In this article we address the question of efficiently solving the algebraic linear system of equations arising from the discretization of a symmetric, elliptic boundary value problem using hp-version discontinuous Galerkin finite element methods. In particular, we introduce a class of domain decomposition preconditioners based on the Schwarz framework, and prove bounds on the condition number of the resulting iteration operators. Numerical results confirming the theoretical estimates are also presented.

High-order hp-adaptive discontinuous Galerkin finite element methods for compressible fluid flows

This article is concerned with the construction of general isotropic and anisotropic adaptive strategies, as well as hp-mesh refinement techniques, in combination with dual-weighted-residual a posteriori error indicators for the discontinuous Galerkin finite element discretization of compressible fluid flow problems.

Development of a Coupling Model for Fluid-Structure Interaction using the Mesh-free Finite Element Method and the Lattice Boltzmann Method

In the presented thesis work, the meshfree method with distance fields was coupled with the lattice Boltzmann method to obtain solutions of fluid-structure interaction problems. The thesis work involved development and implementation of numerical algorithms, data structure, and software. Numerical and computational properties of the coupling algorithm combining the meshfree method with distance fields and the lattice Boltzmann method were investigated. Convergence and accuracy of the methodology was validated by analytical solutions. The research was focused on fluid-structure interaction solutions in complex, mesh-resistant domains as both the lattice Boltzmann method and the meshfree method with distance fields are particularly adept in these situations. Furthermore, the fluid solution provided by the lattice Boltzmann method is massively scalable, allowing extensive use of cutting edge parallel computing resources to accelerate this phase of the solution process. The meshfree method with distance fields allows for exact satisfaction of boundary conditions making it possible to exactly capture the effects of the fluid field on the solid structure.

Frequency adjustable MEMS vibration energy harvester

Ambient mechanical vibrations offer an attractive solution for powering the wireless sensor nodes of the emerging "Internet-of-Things". However, the wide-ranging variability of the ambient vibration frequencies pose a significant challenge to the efficient transduction of vibration into usable electrical energy. This work reports the development of a MEMS electromagnetic vibration energy harvester where the resonance frequency of the oscillator can be adjusted or tuned to adapt to the ambient vibrational frequency. Micro-fabricated silicon spring and double layer planar micro-coils along with sintered NdFeB micro-magnets are used to construct the electromagnetic transduction mechanism. Furthermore, another NdFeB magnet is adjustably assembled to induce variable magnetic interaction with the transducing magnet, leading to significant change in the spring stiffness and resonance frequency. Finite element analysis and numerical simulations exhibit substantial frequency tuning range (25% of natural resonance frequency) by appropriate adjustment of the repulsive magnetic interaction between the tuning and transducing magnet pair. This demonstrated method of frequency adjustment or tuning have potential applications in other MEMS vibration energy harvesters and micromechanical oscillators.

Characterisation and investigation of local variations in mechanical behaviour in cast aluminium using gradient solidification, Digital Image Correlation and finite element simulation

Due to design and process-related factors, there are local variations in the microstructure and mechanical behaviour of cast components. This work establishes a Digital Image Correlation (DIC) based method for characterisation and investigation of the effects of such local variations on the behaviour of a high pressure, die cast (HPDC) aluminium alloy. Plastic behaviour is studied using gradient solidified samples and characterisation models for the parameters of the Hollomon equation are developed, based on microstructural refinement. Samples with controlled microstructural variations are produced and the observed DIC strain field is compared with Finite Element Method (FEM) simulation results. The results show that the DIC based method can be applied to characterise local mechanical behaviour with high accuracy. The microstructural variations are observed to cause a redistribution of strain during tensile loading. This redistribution of strain can be predicted in the FEM simulation by incorporating local mechanical behaviour using the developed characterization model. A homogeneous FEM simulation is unable to predict the observed behaviour. The results motivate the application of a previously proposed simulation strategy, which is able to predict and incorporate local variations in mechanical behaviour into FEM simulations already in the design process for cast components.

Fatigue Design Challenges of Transition Pieces from Decommissioned Platforms for Offshore Wind Energy

All structures are subjected to various loading conditions and combinations. For offshore structures, these loads include permanent loads, hydrostatic pressure, wave, current, and wind loads. Typically, sea environments in different geographical regions are characterized by the 100-year wave height, surface currents, and velocity speeds. The main problems associated with the commonly used, deterministic method is the fact that not all waves have the same period, and that the actual stochastic nature of the marine environment is not taken into account. Offshore steel structure fatigue design is done using the DNVGL-RP-0005:2016 standard which takes precedence over the DNV-RP-C203 standard (2012). Fatigue analysis is necessary for oil and gas producing offshore steel structures which were first constructed in the Gulf of Mexico North Sea (the 1930s) and later in the North Sea (1960s). Fatigue strength is commonly described by S-N curves which have been obtained by laboratory experiments. The rapid development of the Offshore wind industry has caused the exploration into deeper ocean areas and the adoption of new support structural concepts such as full lattice tower systems amongst others. The optimal design of offshore wind support structures including foundation, turbine towers, and transition piece components putting into consideration, economy, safety, and even the environment is a critical challenge. In this study, fatigue design challenges of transition pieces from decommissioned platforms for offshore wind energy are proposed to be discussed. The fatigue resistance of the material and structural components under uniaxial and multiaxial loading is introduced with the new fatigue design rules whilst considering the combination of global and local modeling using finite element analysis software programs.

Biofilm And Saliva Affect The Biomechanical Behavior Of Dental Implants.

Friction coefficient (FC) was quantified between titanium-titanium (Ti-Ti) and titanium-zirconia (Ti-Zr), materials commonly used as abutment and implants, in the presence of a multispecies biofilm (Bf) or salivary pellicle (Pel). Furthermore, FC was used as a parameter to evaluate the biomechanical behavior of a single implant-supported restoration. Interface between Ti-Ti and Ti-Zr without Pel or Bf was used as control (Ctrl). FC was recorded using tribometer and analyzed by two-way Anova and Tukey test (p<0.05). Data were transposed to a finite element model of a dental implant-supported restoration. Models were obtained varying abutment material (Ti and Zr) and FCs recorded (Bf, Pel, and Ctrl). Maximum and shear stress were calculated for bone and equivalent von Misses for prosthetic components. Data were analyzed using two-way ANOVA (p<0.05) and percentage of contribution for each condition (material and FC) was calculated. FC significant differences were observed between Ti-Ti and Ti-Zr for Ctrl and Bf groups, with lower values for Ti-Zr (p<0.05). Within each material group, Ti-Ti differed between all treatments (p<0.05) and for Ti-Zr, only Pel showed higher values compared with Ctrl and Bf (p<0.05). FC contributed to 89.83% (p<0.05) of the stress in the screw, decreasing the stress when the FC was lower. FC resulted in an increase of 59.78% of maximum stress in cortical bone (p=0.05). It can be concluded that the shift of the FC due to the presence of Pel or Bf is able to jeopardize the biomechanical behavior of a single implant-supported restoration.

Influence of specimens' design and manufacturing process on microtensile bond strength to enamel: laboratory and FEA comparison

This study evaluated the effect of specimens' design and manufacturing process on microtensile bond strength, internal stress distributions (Finite Element Analysis - FEA) and specimens' integrity by means of Scanning Electron Microscopy (SEM) and Laser Scanning Confocal Microscopy (LCM). Excite was applied to flat enamel surface and a resin composite build-ups were made incrementally with 1-mm increments of Tetric Ceram. Teeth were cut using a diamond disc or a diamond wire, obtaining 0.8 mm² stick-shaped specimens, or were shaped with a Micro Specimen Former, obtaining dumbbell-shaped specimens (n = 10). Samples were randomly selected for SEM and LCM analysis. Remaining samples underwent microtensile test, and results were analyzed with ANOVA and Tukey test. FEA dumbbell-shaped model resulted in a more homogeneous stress distribution. Nonetheless, they failed under lower bond strengths (21.83 ± 5.44 MPa)c than stick-shaped specimens (sectioned with wire: 42.93 ± 4.77 MPaª; sectioned with disc: 36.62 ± 3.63 MPa b), due to geometric irregularities related to manufacturing process, as noted in microscopic analyzes. It could be concluded that stick-shaped, nontrimmed specimens, sectioned with diamond wire, are preferred for enamel specimens as they can be prepared in a less destructive, easier, and more precise way.

The role of crystalline anisotropy in mechanical property extractions through Berkovich indentation

This work uses crystal plasticity finite element simulations to elucidate the role of elastoplastic anisotropy in instrumented indentation P-h(s) curve measurements in face-centered Cubic (fcc) crystals. It is shown that although the experimental fluctuations in the loading stage of the P-h(s) curves can be attributed to anisotropy, the variability in the unloading stage of the experiments Is much greater than that resulting from anisotropy alone. Moreover, it is found that the conventional procedure used to evaluate the contact variables ruling the unloading P-h(s) curve introduces all uncertainty that approximates to the more fundamental influence of anisotropy. In view of these results, a robust procedure is proposed that uses contact area measurements in addition to the P-h(s) curves to extract homogenized J(2)-Plasticity-equivalent mechanical properties from single crystals.

Chaos-Galerkin solution of stochastic Timoshenko bending problems

This paper presents an accurate and efficient solution for the random transverse and angular displacement fields of uncertain Timoshenko beams. Approximate, numerical solutions are obtained using the Galerkin method and chaos polynomials. The Chaos-Galerkin scheme is constructed by respecting the theoretical conditions for existence and uniqueness of the solution. Numerical results show fast convergence to the exact solution, at excellent accuracies. The developed Chaos-Galerkin scheme accurately approximates the complete cumulative distribution function of the displacement responses. The Chaos-Galerkin scheme developed herein is a theoretically sound and efficient method for the solution of stochastic problems in engineering. (C) 2011 Elsevier Ltd. All rights reserved.

Analytical solutions for geometrically nonlinear trusses

This paper presents an analytical method for analyzing trusses with severe geometrically nonlinear behavior. The main objective is to find analytical solutions for trusses with different axial forces in the bars. The methodology is based on truss kinematics, elastic constitutive laws and equilibrium of nodal forces. The proposed formulation can be applied to hyper elastic materials, such as rubber and elastic foams. A Von Mises truss with two bars made by different materials is analyzed to show the accuracy of this methodology.

A geometrical approach of 3-D FEA for educational purposes applied to electrostatic fields

A geometrical approach of the finite-element analysis applied to electrostatic fields is presented. This approach is particularly well adapted to teaching Finite Elements in Electrical Engineering courses at undergraduate level. The procedure leads to the same system of algebraic equations as that derived by classical approaches, such as variational principle or weighted residuals for nodal elements with plane symmetry. It is shown that the extension of the original procedure to three dimensions is straightforward, provided the domain be meshed in first-order tetrahedral elements. The element matrices are derived by applying Maxwell`s equations in integral form to suitably chosen surfaces in the finite-element mesh.

Effects of elastic indenter deformation on spherical instrumented indentation tests: the reduced elastic modulus

Although the Hertz theory is not applicable in the analysis of the indentation of elastic-plastic materials, it is common practice to incorporate the concept of indenter/specimen combined modulus to consider indenter deformation. The appropriateness was assessed of the use of reduced modulus to incorporate the effect of indenter deformation in the analysis of the indentation with spherical indenters. The analysis based on finite element simulations considered four values of the ratio of the indented material elastic modulus to that of the diamond indenter, E/E(i) (0, 0.04, 0.19, 0.39), four values of the ratio of the elastic reduced modulus to the initial yield strength, E(r)/Y (0, 10, 20, 100), and two values of the ratio of the indenter radius to maximum total displacement, R/delta(max) (3, 10). Indenter deformation effects are better accounted for by the reduced modulus if the indented material behaves entirely elastically. In this case, identical load-displacement (P - delta) curves are obtained with rigid and elastic spherical indenters for the same elastic reduced modulus. Changes in the ratio E/E(i), from 0 to 0.39, resulted in variations lower than 5% for the load dimensionless functions, lower than 3% in the contact area, A(c), and lower than 5% in the ratio H/E(r). However, deformations of the elastic indenter made the actual radius of contact change, even in the indentation of elastic materials. Even though the load dimensionless functions showed only a little increase with the ratio E/E(i), the hardening coefficient and the yield strength could be slightly overestimated when algorithms based on rigid indenters are used. For the unloading curves, the ratio delta(e)/delta(max), where delta(e) is the point corresponding to zero load of a straight line with slope S from the point (P(max), delta(max)), varied less than 5% with the ratio E/E(i). Similarly, the relationship between reduced modulus and the unloading indentation curve, expressed by Sneddon`s equation, did not reveal the necessity of correction with the ratio E/E(i). The most affected parameter in the indentation curve, as a consequence of the indentation deformation, was the ratio between the residual indentation depth after complete unloading and the maximum indenter displacement, delta(r)/delta(max) (up to 26%), but this variation did not significantly decrease the capability to estimate hardness and elastic modulus based on the ratio of the residual indentation depth to maximum indentation depth, h(r)/h(max). In general, the results confirm the convenience of the use of the reduced modulus in the spherical instrumented indentation tests.

A study of plastic deformation around a defect using the magnetic Barkhausen noise in ASTM 36 steel

The present work presents the measurements of the magnetic Barkhausen noise (MBN) in ASTM 36 steel samples around a pit under plastic deformation. The contour maps obtained from these Barkhausen noise measurements are compared with the finite element analysis of the ideal plastic deformation. Also, a parameter of the Barkhausen signal to detect the plastic deformation around the pit in ASTM 36 steel is obtained. Additionally to that, we propose another MBN parameter to estimate the pit width using the Barkhausen noise. (c) 2007 Elsevier Ltd. All rights reserved.

Dynamic Design of Piezoelectric Laminated Sensors and Actuators using Topology Optimization

Sensors and actuators based on piezoelectric plates have shown increasing demand in the field of smart structures, including the development of actuators for cooling and fluid-pumping applications and transducers for novel energy-harvesting devices. This project involves the development of a topology optimization formulation for dynamic design of piezoelectric laminated plates aiming at piezoelectric sensors, actuators and energy-harvesting applications. It distributes piezoelectric material over a metallic plate in order to achieve a desired dynamic behavior with specified resonance frequencies, modes, and enhanced electromechanical coupling factor (EMCC). The finite element employs a piezoelectric plate based on the MITC formulation, which is reliable, efficient and avoids the shear locking problem. The topology optimization formulation is based on the PEMAP-P model combined with the RAMP model, where the design variables are the pseudo-densities that describe the amount of piezoelectric material at each finite element and its polarization sign. The design problem formulated aims at designing simultaneously an eigenshape, i.e., maximizing and minimizing vibration amplitudes at certain points of the structure in a given eigenmode, while tuning the eigenvalue to a desired value and also maximizing its EMCC, so that the energy conversion is maximized for that mode. The optimization problem is solved by using sequential linear programming. Through this formulation, a design with enhancing energy conversion in the low-frequency spectrum is obtained, by minimizing a set of first eigenvalues, enhancing their corresponding eigenshapes while maximizing their EMCCs, which can be considered an approach to the design of energy-harvesting devices. The implementation of the topology optimization algorithm and some results are presented to illustrate the method.