A finite element head injury model. Volume II: computer program documentation. Final report.

National Highway Traffic Safety Administration, Washington, D.C.

Rapid solution of finite element equations on locally refined grids by multi-level methods /

"UILU-ENG 80 1712."

The generalized element method /

"UILU-ENG 78 1738."

Analysis of underground openings in rock by finite element methods : final report /

"This research was supported by the Advanced Research Projects Agency of the Department of Defense and was monitored by the Bureau of Mines ..."

Implementation of a bond model, including dilation, for reinforced materials in a finite element analysis /

Cover title.

Finite element analysis of fault bend influence on stick-slip instability along an intra-plate fault

Earthquakes have been recognized as resulting from stick-slip frictional instabilities along the faults between deformable rocks. A three-dimensional finite-element code for modeling the nonlinear frictional contact behaviors between deformable bodies with the node-to-point contact element strategy has been developed and applied here to investigate the fault geometry influence on the nucleation and development process of the stick-slip instability along an intra-plate fault through a typical fault bend model, which has a pre-cut fault that is artificially bent by an angle of 5.6degrees at the fault center. The numerical results demonstrate that the geometry of the fault significantly affects nucleation, termination and restart of the stick-slip instability along the intra-plate fault, and all these instability phenomena can be well simulated using the current finite-element algorithm.

Acoustic Propagation over Periodic and Quasi-Period Rough Surfaces Proceedings

Numerical techniques such as the Boundary Element Method, Finite Element Method and Finite Difference Time Domain have been used widely to investigate plane and curved wave-front scattering by rough surfaces. For certain shapes of roughness elements (cylinders, semi-cylinders and ellipsoids) there are semi-analytical alternatives. Here, we present a theory for multiple scattering by cylinders on a hard surface to investigate effects due to different roughness shape, the effects of vacancies and variation of roughness element size on the excess attenuation due to a periodically rough surfaces.

Three-dimensional finite element thermal analysis of dental tissues irradiated with Er,Cr:YSGG laser

In the present study, a finite element model of a half-sectioned molar tooth was developed in order to understand the thermal behavior of dental hard tissues (both enamel and dentin) under laser irradiation. The model was validated by comparing it with an in vitro experiment where a sound molar tooth was irradiated by an Er,Cr:YSGG pulsed laser. The numerical tooth model was conceived to simulate the in vitro experiment, reproducing the dimensions and physical conditions of the typical molar sound tooth, considering laser energy absorption and calculating the heat transfer through the dental tissues in three dimensions. The numerical assay considered the same three laser energy densities at the same wavelength (2.79 mu m) used in the experiment. A thermographic camera was used to perform the in vitro experiment, in which an Er, Cr: YSGG laser (2.79 mu m) was used to irradiate tooth samples and the infrared images obtained were stored and analyzed. The temperature increments in both the finite element model and the in vitro experiment were compared. The distribution of temperature inside the tooth versus time plotted for two critical points showed a relatively good agreement between the results of the experiment and model. The three dimensional model allows one to understand how the heat propagates through the dentin and enamel and to relate the amount of energy applied, width of the laser pulses, and temperature inside the tooth. (C) 2008 American Institute of Physics. [DOI: 10.1063/1.2953526]

An exact conserving algorithm for nonlinear dynamics with rotational DOFs and general hyperelasticity. Part 1: Rods

A fully conserving algorithm is developed in this paper for the integration of the equations of motion in nonlinear rod dynamics. The starting point is a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, which results in an extremely simple update of the rotational variables. The weak form is constructed with a non-orthogonal projection corresponding to the application of the virtual power theorem. Together with an appropriate time-collocation, it ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that nonlinear hyperelastic materials (and not only materials with quadratic potentials) are permitted without any prejudice on the conservation properties. Spatial discretization is performed via the finite element method and the performance of the scheme is assessed by means of several numerical simulations.

Identification of Elastic, Dielectric, and Piezoelectric Constants in Piezoceramic Disks

Three-dimensional modeling of piezoelectric devices requires a precise knowledge of piezoelectric material parameters. The commonly used piezoelectric materials belong to the 6mm symmetry class, which have ten independent constants. In this work, a methodology to obtain precise material constants over a wide frequency band through finite element analysis of a piezoceramic disk is presented. Given an experimental electrical impedance curve and a first estimate for the piezoelectric material properties, the objective is to find the material properties that minimize the difference between the electrical impedance calculated by the finite element method and that obtained experimentally by an electrical impedance analyzer. The methodology consists of four basic steps: experimental measurement, identification of vibration modes and their sensitivity to material constants, a preliminary identification algorithm, and final refinement of the material constants using an optimization algorithm. The application of the methodology is exemplified using a hard lead zirconate titanate piezoceramic. The same methodology is applied to a soft piezoceramic. The errors in the identification of each parameter are statistically estimated in both cases, and are less than 0.6% for elastic constants, and less than 6.3% for dielectric and piezoelectric constants.

Stabilized Petrov-Galerkin methods for the convection-diffusion-reaction and the Helmholtz equations

We present two new stabilized high-resolution numerical methods for the convection–diffusion–reaction (CDR) and the Helmholtz equations respectively. The work embarks upon a priori analysis of some consistency recovery procedures for some stabilization methods belonging to the Petrov–Galerkin framework. It was found that the use of some standard practices (e.g. M-Matrices theory) for the design of essentially non-oscillatory numerical methods is not feasible when consistency recovery methods are employed. Hence, with respect to convective stabilization, such recovery methods are not preferred. Next, we present the design of a high-resolution Petrov–Galerkin (HRPG) method for the 1D CDR problem. The problem is studied from a fresh point of view, including practical implications on the formulation of the maximum principle, M-Matrices theory, monotonicity and total variation diminishing (TVD) finite volume schemes. The current method is next in line to earlier methods that may be viewed as an upwinding plus a discontinuity-capturing operator. Finally, some remarks are made on the extension of the HRPG method to multidimensions. Next, we present a new numerical scheme for the Helmholtz equation resulting in quasi-exact solutions. The focus is on the approximation of the solution to the Helmholtz equation in the interior of the domain using compact stencils. Piecewise linear/bilinear polynomial interpolation are considered on a structured mesh/grid. The only a priori requirement is to provide a mesh/grid resolution of at least eight elements per wavelength. No stabilization parameters are involved in the definition of the scheme. The scheme consists of taking the average of the equation stencils obtained by the standard Galerkin finite element method and the classical finite difference method. Dispersion analysis in 1D and 2D illustrate the quasi-exact properties of this scheme. Finally, some remarks are made on the extension of the scheme to unstructured meshes by designing a method within the Petrov–Galerkin framework.

Deep-hole drilling process, monitoring system and sensors modifications

 In the drilling processes and especially deep-hole drilling process, the monitoring system and having control on mechanical parameters (e.g. Force, Torque,Vibration and Acoustic emission) are essential. The main focus of this thesis work is to study the characteristics of deep-hole drilling process, and optimize the monitoring system for controlling the process. The vibration is considered as a major defect area of the deep-hole drilling process which often leads to breakage of the drill, therefore by vibration analysis and optimizing the workpiecefixture, this area is studied by finite element method and the suggestions are explained. By study on a present monitoring system, and searching on the new sensor products, the modifications and recommendations are suggested for optimize the present monitoring system for excellent performance in deep-hole drilling process research and measurements.

Consistent Dirichlet Boundary Conditions for Numerical Solution of Moving Boundary Problems

We consider the imposition of Dirichlet boundary conditions in the finite element modelling of moving boundary problems in one and two dimensions for which the total mass is prescribed. A modification of the standard linear finite element test space allows the boundary conditions to be imposed strongly whilst simultaneously conserving a discrete mass. The validity of the technique is assessed for a specific moving mesh finite element method, although the approach is more general. Numerical comparisons are carried out for mass-conserving solutions of the porous medium equation with Dirichlet boundary conditions and for a moving boundary problem with a source term and time-varying mass.

Sistemas dinâmicos com amortecedores ativos controlados por atuadores piezelétricos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Análise eletromagnética de cabos OPGW utilizando métodos numérico e analítico

Este trabalho apresenta o estudo eletromagnético de cabos OPGW (Optical Ground Wire) os quais têm dupla função: de pára-raios para linhas de transmissão de alta tensão e de canal de comunicação através de fibras ópticas embutidas na estrutura do cabo. Descargas atmosféricas ou curtos-circuitos podem comprometer a integridade do cabo, devido ao aquecimento nas regiões onde há maior concentração de corrente. Para a análise deste problema foram feitos cálculos eletromagnéticos relacionando-os aos efeitos térmicos no cabo. Nesta análise foram consideradas três diferentes geometrias: o modelo de cabo real, o modelo de cabo com camadas homogêneas e o modelo de cabo com uma camada modificada; esta modificação está relacionada à forma geométrica dos fios da armação do cabo. As ferramentas utilizadas em tal estudo foram o software comercial FEMLAB Multiphysics, baseado no método dos elementos finitos, e um método analítico desenvolvido a partir das equações de Maxwell no domínio da freqüência, que foi implementado utilizando o software MATLAB. Os principais resultados deste trabalho são gráficos de distribuição de densidade de corrente na seção reta do cabo para diferentes freqüências, estudo do efeito pelicular e do efeito de proximidade entre os condutores do cabo.