Analysis of reflection characteristics of a normal incidence plane wave on resonant sound absorbers: A finite element approach

Resonant sound absorbers are used widely as anechoic coatings in underwater applications. In this paper a finite element scheme based on the Galerkin technique is used to analyze the reflection characteristics of the resonant absorber when insonified by a normal incidence plane wave. A waveguide theory coupled with an impedance matching condition in the fluid is used to model the problem. It is shown in this paper that the fluid medium encompassing the absorber can be modeled as an elastic medium with equivalent Lamé constants. Quarter symmetry conditions within the periodic unit cell are exploited. The finite element results are compared with analytical results, and with results published elsewhere in the literature. It is shown in the process that meshing of the fluid domain can be obviated if the transmission coefficients or reflection coefficients only are desired as is often the case. Finally, some design curves for thin resonant absorbers with water closure are presented in this paper.

Finite element simulation of BAW propagation in inhomogeneous plate due to piezoelectric actuation

A set of finite elements (FEs) is formulated to analyze wave propagation through inhomogeneous material when subjected to mechanical, thermal loading or piezo-electric actuation. Elastic, thermal and electrical properties of the materials axe allowed to vary in length and thickness direction. The elements can act both as sensors and actuators. These elements are used to model wave propagation in functionally graded materials (FGM) and the effect of inhomogeneity in the wave is demonstrated. Further, a surface acoustic wave (SAW) device is modeled and wave propagation due to piezo-electric actuation from interdigital transducers (IDTs) is studied.

Wave propagation analysis in laminated composite plates with transverse cracks using the wavelet spectral finite element method

This paper presents a newly developed wavelet spectral finite element (WFSE) model to analyze wave propagation in anisotropic composite laminate with a transverse surface crack penetrating part-through the thickness. The WSFE formulation of the composite laminate, which is based on the first-order shear deformation theory, produces accurate and computationally efficient results for high frequency wave motion. Transverse crack is modeled in wavenumber-frequency domain by introducing bending flexibility of the plate along crack edge. Results for tone burst and impulse excitations show excellent agreement with conventional finite element analysis in Abaqus (R). Problems with multiple cracks are modeled by assembling a number of spectral elements with cracks in frequency-wavenumber domain. Results show partial reflection of the excited wave due to crack at time instances consistent with crack locations. (C) 2014 Elsevier B.V. All rights reserved.

A New Finite Element Method for Strain Gradient Theories and Applications to Fracture Analyses

A new compatible finite element method for strain gradient theories is presented. In the new finite element method, pure displacement derivatives are taken as the fundamental variables. The new numerical method is successfully used to analyze the simple strain gradient problems – the fundamental fracture problems. Through comparing the numerical solutions with the existed exact solutions, the effectiveness of the new finite element method is tested and confirmed. Additionally, an application of the Zienkiewicz–Taylor C1 finite element method to the strain gradient problem is discussed. By using the new finite element method, plane-strain mode I and mode II crack tip fields are calculated based on a constitutive law which is a simple generalization of the conventional J2 deformation plasticity theory to include strain gradient effects. Three new constitutive parameters enter to characterize the scale over which strain gradient effects become important. During the analysis the general compressible version of Fleck–Hutchinson strain gradient plasticity is adopted. Crack tip solutions, the traction distributions along the plane ahead of the crack tip are calculated. The solutions display the considerable elevation of traction within the zone near the crack tip.

Finite Element Solutions for Plane Strain Mode I Crack with Strain Grading Effect

In this paper, a new phenomenological theory with strain gradient effects is proposed to account for the size dependence of plastic deformation at micro- and submicro-length scales. The theory fits within the framework of general couple stress theory and three rotational degrees of freedom omega(i) are introduced in addition to the conventional three translational degrees of freedom mu(i). omega(i) is called micro-rotation and is the sum of material rotation plus the particles' relative rotation. While the new theory is used to analyze the crack tip field or the indentation problems, the stretch gradient is considered through a new hardening law. The key features of the theory are that the rotation gradient influences the material character through the interaction between the Cauchy stresses and the couple stresses; the term of stretch gradient is represented as an internal variable to increase the tangent modulus. In fact the present new strain gradient theory is the combination of the strain gradient theory proposed by Chen and Wang (Int. J. Plast., in press) and the hardening law given by Chen and Wang (Acta Mater. 48 (2000a) 3997). In this paper we focus on the finite element method to investigate material fracture for an elastic-power law hardening solid. With remotely imposed classical K fields, the full field solutions are obtained numerically. It is found that the size of the strain gradient dominance zone is characterized by the intrinsic material length l(1). Outside the strain gradient dominance zone, the computed stress field tends to be a classical plasticity field and then K field. The singularity of stresses ahead of the crack tip is higher than that of the classical field and tends to the square root singularity, which has important consequences for crack growth in materials by decohesion at the atomic scale. (C) 2002 Elsevier Science Ltd. All rights reserved.

Analysis of Si/GaAs Bonding Stresses with the Finite Element Method

In conjunction with ANSYS, we use the finite element method to analyze the bonding stresses of Si/GaAs. We also apply a numerical model to investigate a contour map and the distribution of normal stress,shearing stress,and peeling stress,taking into full consideration the thermal expansion coefficient as a function of temperature. Novel bonding structures are proposed for reducing the effect of thermal stress as compared with conventional structures. Calculations show the validity of this new structure.

A finite element solution scheme for an elastic–viscoplastic soil model

A Newton–Raphson solution scheme with a stress point algorithm is presented for the implementation of an elastic–viscoplastic soilmodel in a finite element program. Viscoplastic strain rates are calculated using the stress and volumetric states of the soil. Sub-incrementsof time are defined for each iterative calculation of elastic–viscoplastic stress changes so that their sum adds up to the time incrementfor the load step. This carefully defined ‘iterative time’ ensures that the correct amount of viscoplastic straining is accumulated overthe applied load step. The algorithms and assumptions required to implement the solution scheme are provided. Verification of the solutionscheme is achieved by using it to analyze typical boundary value problems.

Finite element idealization of a cold-formed steel portal frame

A simple linear beam idealization of a cold-formed steel portal frame is presented in which beam elements are used to idealize the column and rafter members, and rotational spring elements are used to represent the rotational flexibility of the joints. In addition, the beam idealization takes into account the finite connection length of the joints. Deflections predicted using the beam idealization are shown to be comparable to deflections obtained from both a linear finite element shell idealization and full-scale laboratory tests. Using the beam idealization, deflections under rafter load are divided into three components: Deflection due to flexure of the column and rafter members, deflection due to bolt-hole elongation, and deflection due to in-plane bracket deformation. Of these deflection components, the deflection due to bolt-hole elongation is the most significant and cannot, therefore, be ignored. Using the beam idealization, engineers can analyze and design cold-formed steel portal frames, including making appropriate allowances for connection effects, without the need to resort to expensive finite element shell analysis.

An Exponentially Convergent Nonpolynomial Finite Element Method for Time-Harmonic Scattering from Polygons

In recent years nonpolynomial finite element methods have received increasing attention for the efficient solution of wave problems. As with their close cousin the method of particular solutions, high efficiency comes from using solutions to the Helmholtz equation as basis functions. We present and analyze such a method for the scattering of two-dimensional scalar waves from a polygonal domain that achieves exponential convergence purely by increasing the number of basis functions in each element. Key ingredients are the use of basis functions that capture the singularities at corners and the representation of the scattered field towards infinity by a combination of fundamental solutions. The solution is obtained by minimizing a least-squares functional, which we discretize in such a way that a matrix least-squares problem is obtained. We give computable exponential bounds on the rate of convergence of the least-squares functional that are in very good agreement with the observed numerical convergence. Challenging numerical examples, including a nonconvex polygon with several corner singularities, and a cavity domain, are solved to around 10 digits of accuracy with a few seconds of CPU time. The examples are implemented concisely with MPSpack, a MATLAB toolbox for wave computations with nonpolynomial basis functions, developed by the authors. A code example is included.

First order k-th moment finite element analysis of nonlinear operator equations with stochastic data

We develop and analyze a class of efficient Galerkin approximation methods for uncertainty quantification of nonlinear operator equations. The algorithms are based on sparse Galerkin discretizations of tensorized linearizations at nominal parameters. Specifically, we consider abstract, nonlinear, parametric operator equations J(\alpha ,u)=0 for random input \alpha (\omega ) with almost sure realizations in a neighborhood of a nominal input parameter \alpha _0. Under some structural assumptions on the parameter dependence, we prove existence and uniqueness of a random solution, u(\omega ) = S(\alpha (\omega )). We derive a multilinear, tensorized operator equation for the deterministic computation of k-th order statistical moments of the random solution's fluctuations u(\omega ) - S(\alpha _0). We introduce and analyse sparse tensor Galerkin discretization schemes for the efficient, deterministic computation of the k-th statistical moment equation. We prove a shift theorem for the k-point correlation equation in anisotropic smoothness scales and deduce that sparse tensor Galerkin discretizations of this equation converge in accuracy vs. complexity which equals, up to logarithmic terms, that of the Galerkin discretization of a single instance of the mean field problem. We illustrate the abstract theory for nonstationary diffusion problems in random domains.

Finite element modeling and analysis of thrombosis in middle cerebral artery

Thrombotic stroke, which is caused by blood clot in the cerebral artery, is a major source of increased mortality and morbidity. Considering as efficient and fastest methods, mathematical approaches have gained significant importance for analyzing and understanding the biological events like thrombosis. This paper presents a computational model to analyze the effects of thrombosis using the theory of coupled fluid dynamics-structure interaction. The finite element method is used for the modeling of thrombosis (blood clot) of different stages in the middle cerebral artery with physiological compliance. The developed model is used to investigate the consequences that occur due to the various sizes of clots in the artery in the form of blood flow velocity, blood pressure, and artery wall stress. Such numerical assessment will facilitate better understanding of the biophysical process in case of thrombosis and thus would support medical practitioners to take faster curing steps.

Periodontal ligament influence on the stress distribution in a removable partial denture supported by implant: a finite element analysis

Objective: The non-homogenous aspect of periodontal ligament (PDL) has been examined using finite element analysis (FEA) to better simulate PDL behavior. The aim of this study was to assess, by 2-D FEA, the influence of non-homogenous PDL on the stress distribution when the free-end saddle removable partial denture (RPD) is partially supported by an osseointegrated implant. Material and Methods: Six finite element (FE) models of a partially edentulous mandible were created to represent two types of PDL (non-homogenous and homogenous) and two types of RPD (conventional RPD, supported by tooth and fibromucosa; and modified RPD, supported by tooth and implant [10.00x3.75 mm]). Two additional FE models without RPD were used as control models. The non-homogenous PDL was modeled using beam elements to simulate the crest, horizontal, oblique and apical fibers. The load (50 N) was applied in each cusp simultaneously. Regarding boundary conditions the border of alveolar ridge was fixed along the x axis. The FE software (Ansys 10.0) was used to compute the stress fields, and the von Mises stress criterion (sigma vM) was applied to analyze the results. Results: The peak of sigma vM in non-homogenous PDL was higher than that for the homogenous condition. The benefits of implants were enhanced for the non-homogenous PDL condition, with drastic sigma vM reduction on the posterior half of the alveolar ridge. The implant did not reduce the stress on the support tooth for both PDL conditions. Conclusion: The PDL modeled in the non-homogeneous form increased the benefits of the osseointegrated implant in comparison with the homogeneous condition. Using the non-homogenous PDL, the presence of osseointegrated implant did not reduce the stress on the supporting tooth.

Micromechanics of dentin/adhesive interface in function of dentin depth: 3d finite element analysis

Objectives: The aim of this study was to analyze the stress distribution on dentin/adhesive interface (d/a) through a 3-D finite element analysis (FEA) varying the number and diameter of the dentin tubules orifice according to dentin depth, keeping hybrid layer (HL) thickness and TAǴs length constant. Materials and Methods: 3 models were built through the SolidWorks software: SD - specimen simulating superficial dentin (41 x 41 x 82 μm), with a 3 μm thick HL, a 17 μm length Tag, and 8 tubules with a 0.9 μm diameter restored with composite resin. MD - similar to M1 with 12 tubules with a 1.2 μm diameter, simulating medium dentin. DD - similar to M1 with 16 tubules with a 2.5 μm diameter, simulating deep dentin. Other two models were built in order to keep the diameter constant in 2.5 μm: MS - similar to SD with 8 tubules; and MM - similar to MD with 12 tubules. The boundary condition was applied to the base surface of each specimen. Tensile load (0.03N) was performed on the composite resin top surface. Stress field (maximum principal stress in tension - σMAX) was performed using Ansys Wokbench 10.0. Results: The peak of σMAX (MPa) were similar between SD (110) and MD (106), and higher for DD (134). The stress distribution pathway was similar for all models, starting from peritubular dentin to adhesive layer, intertubular dentin and hybrid layer. The peak of σMAX (MPa) for those structures was, respectively: 134 (DD), 56.9 (SD), 45.5 (DD), and 36.7 (MD). Conclusions: The number of dentin tubules had no influence in the σMAX at the dentin/adhesive interface. Peritubular and intertubular dentin showed higher stress with the bigger dentin tubules orifice condition. The σMAX in the hybrid layer and adhesive layer were going down from superficial dentin to deeper dentin. In a failure scenario, the hybrid layer in contact with peritubular dentin and adhesive layer is the first region for breaking the adhesion. © 2011 Nova Science Publishers, Inc.

Stress distribution on dentin-cement-post interface varying root canal and glass fiber post diameters. A three-dimensional finite element analysis based on micro-C data

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Influence of Tapered and External Hexagon Connections on Bone Stresses Around Tilted Dental Implants: Three-Dimensional Finite Element Method With Statistical Analysis

Background: The purpose of this study is to analyze the tension distribution on bone tissue around implants with different angulations (0 degrees, 17 degrees, and 30 degrees) and connections (external hexagon and tapered) through the use of three-dimensional finite element and statistical analyses.Methods: Twelve different configurations of three-dimensional finite element models, including three inclinations of the implants (0 degrees, 17 degrees, and 30 degrees), two connections (an external hexagon and a tapered), and two load applications (axial and oblique), were simulated. The maximum principal stress values for cortical bone were measured at the mesial, distal, buccal, and lingual regions around the implant for each analyzed situation, totaling 48 groups. Loads of 200 and 100 N were applied at the occlusal surface in the axial and oblique directions, respectively. Maximum principal stress values were measured at the bone crest and statistically analyzed using analysis of variance. Stress patterns in the bone tissue around the implant were analyzed qualitatively.Results: The results demonstrated that under the oblique loading process, the external hexagon connection showed significantly higher stress concentrations in the bone tissue (P < 0.05) compared with the tapered connection. Moreover, the buccal and mesial regions of the cortical bone concentrated significantly higher stress (P < 0.005) to the external hexagon implant type. Under the oblique loading direction, the increased external hexagon implant angulation induced a significantly higher stress concentration (P = 0.045).Conclusions: The study results show that: 1) the oblique load was more damaging to bone tissue, mainly when associated with external hexagon implants; and 2) there was a higher stress concentration on the buccal region in comparison to all other regions under oblique load.