Modelle und Prototypen für die Medizin

Der Anwendungsbereich für Modelle und Prototypen hat sich mittlerweile auch auf unterschiedlichste medizinische Fragestellungen ausgedehnt. Der vorliegende Beitrag zeigt an verschiedenen Praxisbeispielen und Verfahren Möglichkeiten eines in der Technik etablierten Verfahrens für eine erweiterte Anwendung auf.

Generic Finite Element Model for Mechanically Consistent Scaling of Composite Beam-to-column Joint with Welded Connection

To study the behaviour of beam-to-column composite connection more sophisticated finite element models is required, since component model has some severe limitations. In this research a generic finite element model for composite beam-to-column joint with welded connections is developed using current state of the art local modelling. Applying mechanically consistent scaling method, it can provide the constitutive relationship for a plane rectangular macro element with beam-type boundaries. Then, this defined macro element, which preserves local behaviour and allows for the transfer of five independent states between local and global models, can be implemented in high-accuracy frame analysis with the possibility of limit state checks. In order that macro element for scaling method can be used in practical manner, a generic geometry program as a new idea proposed in this study is also developed for this finite element model. With generic programming a set of global geometric variables can be input to generate a specific instance of the connection without much effort. The proposed finite element model generated by this generic programming is validated against testing results from University of Kaiserslautern. Finally, two illustrative examples for applying this macro element approach are presented. In the first example how to obtain the constitutive relationships of macro element is demonstrated. With certain assumptions for typical composite frame the constitutive relationships can be represented by bilinear laws for the macro bending and shear states that are then coupled by a two-dimensional surface law with yield and failure surfaces. In second example a scaling concept that combines sophisticated local models with a frame analysis using a macro element approach is presented as a practical application of this numerical model.

The effect of elastic modulus and equilibrium vapor adsorption on capillary forces

Die Kapillarkraft entsteht durch die Bildung eines Meniskus zwischen zwei Festkörpen. In dieser Doktorarbeit wurden die Auswirkungen von elastischer Verformung und Flϋssigkeitadsorption auf die Kapillarkraft sowohl theoretisch als auch experimentell untersucht. Unter Verwendung eines Rasterkraftmikroskops wurde die Kapillarkraft zwischen eines Siliziumoxid Kolloids von 2 µm Radius und eine weiche Oberfläche wie n.a. Polydimethylsiloxan oder Polyisopren, unter normalen Umgebungsbedingungen sowie in variierende Ethanoldampfdrϋcken gemessen. Diese Ergebnisse wurden mit den Kapillarkräften verglichen, die auf einem harten Substrat (Silizium-Wafer) unter denselben Bedingungen gemessen wurden. Wir beobachteten eine monotone Abnahme der Kapillarkraft mit zunehmendem Ethanoldampfdruck (P) fϋr P/Psat > 0,2, wobei Psat der Sättigungsdampfdruck ist.rnUm die experimentellen Ergebnisse zu erklären, wurde ein zuvor entwickeltes analytisches Modell (Soft Matter 2010, 6, 3930) erweitert, um die Ethanoladsorption zu berϋcksichtigen. Dieses neue analytische Modell zeigte zwei verschiedene Abhängigkeiten der Kapillarkraft von P/Psat auf harten und weichen Oberflächen. Fϋr die harte Oberfläche des Siliziumwafers wird die Abhängigkeit der Kapillarkraft vom Dampfdruck vom Verhältnis der Dicke der adsorbierten Ethanolschicht zum Meniskusradius bestimmt. Auf weichen Polymeroberflächen hingegen hängt die Kapillarkraft von der Oberflächenverformung und des Laplace-Drucks innerhalb des Meniskus ab. Eine Abnahme der Kapillarkraft mit zunehmendem Ethanoldampfdruck hat demnach eine Abnahme des Laplace-Drucks mit zunehmendem Meniskusradius zur folge. rnDie analytischen Berechnungen, fϋr die eine Hertzsche Kontakt-deformation angenommen wurde, wurden mit Finit Element Methode Simulationen verglichen, welche die reale Deformation des elastischen Substrats in der Nähe des Meniskuses explizit berϋcksichtigen. Diese zusätzliche nach oben gerichtete oberflächenverformung im Bereich des Meniskus fϋhrt zu einer weiteren Erhöhung der Kapillarkraft, insbesondere fϋr weiche Oberflächen mit Elastizitätsmodulen < 100 MPa.rn

Tensões ósseas peri-implantares em implantes colocados justa e infra crestalmente

Este trabalho divide-se em duas partes distintas: uma longa e detalhada revisão bibliográfica acerca das temáticas anatomia peri-implantar, espaço biológico, osso alveolar, osteointegração, cone Morse e platform-switching e FEA (Finit Element Analisys) ; e um estudo sobre tensões peri-implantares em implantes do tipo cone Morse colocados infra e justa crestalmente. Foi possível concluir com este estudo laboratorial que os implantes colocados justacrestalmente apresentam melhores resultados biomecanicamente, ou seja, apresentam um menor volume de osso em tensão. Materiais e métodos: Foi realizada uma pesquisa bibliográfica na PubMed e Medline explorando os seguintes items: “osteointegração”, “saucerização”, “platform switching”, “cone Morse”, “osso alveolar”, “anatomina peri-implantar”, “espaço biológico”, “osteoclastos”, “osteoblastos”, “remodelação óssea”, “colocação de implantes justacrestalmente”, “colocação de implantes infra-crestalmente” e “análise de FEA”. Na bibliografia encontrada com as temáticas supra-citadas foi feita uma cuidadosa selecção de acordo com aquilo a que este trabalho se propunha. Simultaneamente, um modelo 3D de dois implantes, um de conexão externa hexagonal e outro de conexão interna do tipo cone Morse, exactamente iguais com exceção da já referida conexão, de 10mm de comprimento e 4mm de diâmetro, foram inseridos num bloco ósseo obtido através de uma CT e sujeitos a uma força axial de 150N e uma força oblíqua de 150N a 45º, tendo sido avaliados por uma análise de elementos finitos.

Finite-element discretization of 3D energy-transport equations for semiconductors

In this thesis a mathematical model was derived that describes the charge and energy transport in semiconductor devices like transistors. Moreover, numerical simulations of these physical processes are performed. In order to accomplish this, methods of theoretical physics, functional analysis, numerical mathematics and computer programming are applied. After an introduction to the status quo of semiconductor device simulation methods and a brief review of historical facts up to now, the attention is shifted to the construction of a model, which serves as the basis of the subsequent derivations in the thesis. Thereby the starting point is an important equation of the theory of dilute gases. From this equation the model equations are derived and specified by means of a series expansion method. This is done in a multi-stage derivation process, which is mainly taken from a scientific paper and which does not constitute the focus of this thesis. In the following phase we specify the mathematical setting and make precise the model assumptions. Thereby we make use of methods of functional analysis. Since the equations we deal with are coupled, we are concerned with a nonstandard problem. In contrary, the theory of scalar elliptic equations is established meanwhile. Subsequently, we are preoccupied with the numerical discretization of the equations. A special finite-element method is used for the discretization. This special approach has to be done in order to make the numerical results appropriate for practical application. By a series of transformations from the discrete model we derive a system of algebraic equations that are eligible for numerical evaluation. Using self-made computer programs we solve the equations to get approximate solutions. These programs are based on new and specialized iteration procedures that are developed and thoroughly tested within the frame of this research work. Due to their importance and their novel status, they are explained and demonstrated in detail. We compare these new iterations with a standard method that is complemented by a feature to fit in the current context. A further innovation is the computation of solutions in three-dimensional domains, which are still rare. Special attention is paid to applicability of the 3D simulation tools. The programs are designed to have justifiable working complexity. The simulation results of some models of contemporary semiconductor devices are shown and detailed comments on the results are given. Eventually, we make a prospect on future development and enhancements of the models and of the algorithms that we used.

Mediterranean-type climate in the South Aegean (Eastern Mediterranean) during the Late Miocene: Evidence from isotope and element proxies

The present-day climate in the Mediterranean region is characterized by mild, wet winters and hot, dry summers. There is contradictory evidence as to whether the present-day conditions (“Mediterranean climate”) already existed in the Late Miocene. This thesis presents seasonally-resolved isotope and element proxy data obtained from Late Miocene reef corals from Crete (Southern Aegean, Eastern Mediterranean) in order to illustrate climate conditions in the Mediterranean region during this time. There was a transition from greenhouse to icehouse conditions without a Greenland ice sheet during the Late Miocene. Since the Greenland ice sheet is predicted to melt fully within the next millennia, Late Miocene climate mechanisms can be considered as useful analogues in evaluating models of Northern Hemispheric climate conditions in the future. So far, high resolution chemical proxy data on Late Miocene environments are limited. In order to enlarge the proxy database for this time span, coral genus Tarbellastraea was evaluated as a new proxy archive, and proved reliable based on consistent oxygen isotope records of Tarbellastraea and the established paleoenvironmental archive of coral genus Porites. In combination with lithostratigraphic data, global 87Sr/86Sr seawater chronostratigraphy was used to constrain the numerical age of the coral sites, assuming the Mediterranean Sea to be equilibrated with global open ocean water. 87Sr/86Sr ratios of Tarbellastraea and Porites from eight stratigraphically different sampling sites were measured by thermal ionization mass spectrometry. The ratios range from 0.708900 to 0.708958 corresponding to ages of 10 to 7 Ma (Tortonian to Early Messinian). Spectral analyses of multi-decadal time-series yield interannual δ18O variability with periods of ~2 and ~5 years, similar to that of modern records, indicating that pressure field systems comparable to those controlling the seasonality of present-day Mediterranean climate existed, at least intermittently, already during the Late Miocene. In addition to sea surface temperature (SST), δ18O composition of coral aragonite is controlled by other parameters such as local seawater composition which as a result of precipitation and evaporation, influences sea surface salinity (SSS). The Sr/Ca ratio is considered to be independent of salinity, and was used, therefore, as an additional proxy to estimate seasonality in SST. Major and trace element concentrations in coral aragonite determined by laser ablation inductively coupled plasma mass spectrometry yield significant variations along a transect perpendicular to coral growth increments, and record varying environmental conditions. The comparison between the average SST seasonality of 7°C and 9°C, derived from average annual δ18O (1.1‰) and Sr/Ca (0.579 mmol/mol) amplitudes, respectively, indicates that the δ18O-derived SST seasonality is biased by seawater composition, reducing the δ18O amplitude by 0.3‰. This value is equivalent to a seasonal SSS variation of 1‰, as observed under present-day Aegean Sea conditions. Concentration patterns of non-lattice bound major and trace elements, related to trapped particles within the coral skeleton, reflect seasonal input of suspended load into the reef environment. δ18O, Sr/Ca and non-lattice bound element proxy records, as well as geochemical compositions of the trapped particles, provide evidence for intense precipitation in the Eastern Mediterranean during winters. Winter rain caused freshwater discharge and transport of weathering products from the hinterland into the reef environment. There is a trend in coral δ18O data to more positive mean δ18O values (–2.7‰ to –1.7‰) coupled with decreased seasonal δ18O amplitudes (1.1‰ to 0.7‰) from 10 to 7 Ma. This relationship is most easily explained in terms of more positive summer δ18O. Since coral diversity and annual growth rates indicate more or less constant average SST for the Mediterranean from the Tortonian to the Early Messinian, more positive mean and summer δ18O indicate increasing aridity during the Late Miocene, and more pronounced during summers. The analytical results implicate that winter rainfall and summer drought, the main characteristics of the present-day Mediterranean climate, were already present in the Mediterranean region during the Late Miocene. Some models have argued that the Mediterranean climate did not exist in this region prior to the Pliocene. However, the data presented here show that conditions comparable to those of the present-day existed either intermittently or permanently since at least about 10 Ma.

Effect Of Platform Connection And Abutment Material On Stress Distribution In Single Anterior Implant-supported Restorations: A Nonlinear 3-dimensional Finite Element Analysis.

Although various abutment connections and materials have recently been introduced, insufficient data exist regarding the effect of stress distribution on their mechanical performance. The purpose of this study was to investigate the effect of different abutment materials and platform connections on stress distribution in single anterior implant-supported restorations with the finite element method. Nine experimental groups were modeled from the combination of 3 platform connections (external hexagon, internal hexagon, and Morse tapered) and 3 abutment materials (titanium, zirconia, and hybrid) as follows: external hexagon-titanium, external hexagon-zirconia, external hexagon-hybrid, internal hexagon-titanium, internal hexagon-zirconia, internal hexagon-hybrid, Morse tapered-titanium, Morse tapered-zirconia, and Morse tapered-hybrid. Finite element models consisted of a 4×13-mm implant, anatomic abutment, and lithium disilicate central incisor crown cemented over the abutment. The 49 N occlusal loading was applied in 6 steps to simulate the incisal guidance. Equivalent von Mises stress (σvM) was used for both the qualitative and quantitative evaluation of the implant and abutment in all the groups and the maximum (σmax) and minimum (σmin) principal stresses for the numerical comparison of the zirconia parts. The highest abutment σvM occurred in the Morse-tapered groups and the lowest in the external hexagon-hybrid, internal hexagon-titanium, and internal hexagon-hybrid groups. The σmax and σmin values were lower in the hybrid groups than in the zirconia groups. The stress distribution concentrated in the abutment-implant interface in all the groups, regardless of the platform connection or abutment material. The platform connection influenced the stress on abutments more than the abutment material. The stress values for implants were similar among different platform connections, but greater stress concentrations were observed in internal connections.

Photoelastic And Finite Element Analyses Of Occlusal Loads In Mandibular Body.

This study proposed to evaluate the mandibular biomechanics in the posterior dentition based on experimental and computational analyses. The analyses were performed on a model of human mandible, which was modeled by epoxy resin for photoelastic analysis and by computer-aided design for finite element analysis. To standardize the evaluation, specific areas were determined at the lateral surface of mandibular body. The photoelastic analysis was configured through a vertical load on the first upper molar and fixed support at the ramus of mandible. The same configuration was used in the computer simulation. Force magnitudes of 50, 100, 150, and 200 N were applied to evaluate the bone stress. The stress results presented similar distribution in both analyses, with the more intense stress being at retromolar area and oblique line and alveolar process at molar level. This study presented the similarity of results in the experimental and computational analyses and, thus, showed the high importance of morphology biomechanical characterization at posterior dentition.

Finite-element Analysis Of 3 Situations Of Trauma In The Human Edentulous Mandible.

Maxillofacial trauma resulting from falls in elderly patients is a major social and health care concern. Most of these traumatic events involve mandibular fractures. The aim of this study was to analyze stress distributions from traumatic loads applied on the symphyseal, parasymphyseal, and mandibular body regions in the elderly edentulous mandible using finite-element analysis (FEA). Computerized tomographic analysis of an edentulous macerated human mandible of a patient approximately 65 years old was performed. The bone structure was converted into a 3-dimensional stereolithographic model, which was used to construct the computer-aided design (CAD) geometry for FEA. The mechanical properties of cortical and cancellous bone were characterized as isotropic and elastic structures, respectively, in the CAD model. The condyles were constrained to prevent free movement in the x-, y-, and z-axes during simulation. This enabled the simulation to include the presence of masticatory muscles during trauma. Three different simulations were performed. Loads of 700 N were applied perpendicular to the surface of the cortical bone in the symphyseal, parasymphyseal, and mandibular body regions. The simulation results were evaluated according to equivalent von Mises stress distributions. Traumatic load at the symphyseal region generated low stress levels in the mental region and high stress levels in the mandibular neck. Traumatic load at the parasymphyseal region concentrated the resulting stress close to the mental foramen. Traumatic load in the mandibular body generated extensive stress in the mandibular body, angle, and ramus. FEA enabled precise mapping of the stress distribution in a human elderly edentulous mandible (neck and mandibular angle) in response to 3 different traumatic load conditions. This knowledge can help guide emergency responders as they evaluate patients after a traumatic event.

Keratoconus prediction using a finite element model of the cornea with local biomechanical properties

PURPOSE: The ability to predict and understand which biomechanical properties of the cornea are responsible for the stability or progression of keratoconus may be an important clinical and surgical tool for the eye-care professional. We have developed a finite element model of the cornea, that tries to predicts keratoconus-like behavior and its evolution based on material properties of the corneal tissue. METHODS: Corneal material properties were modeled using bibliographic data and corneal topography was based on literature values from a schematic eye model. Commercial software was used to simulate mechanical and surface properties when the cornea was subject to different local parameters, such as elasticity. RESULTS: The simulation has shown that, depending on the corneal initial surface shape, changes in local material properties and also different intraocular pressures values induce a localized protuberance and increase in curvature when compared to the remaining portion of the cornea. CONCLUSIONS: This technique provides a quantitative and accurate approach to the problem of understanding the biomechanical nature of keratoconus. The implemented model has shown that changes in local material properties of the cornea and intraocular pressure are intrinsically related to keratoconus pathology and its shape/curvature.

Stress distribution in the cervical region of an upper central incisor in a 3D finite element model

The aim of this study was to evaluate the stress distribution in the cervical region of a sound upper central incisor in two clinical situations, standard and maximum masticatory forces, by means of a 3D model with the highest possible level of fidelity to the anatomic dimensions. Two models with 331,887 linear tetrahedral elements that represent a sound upper central incisor with periodontal ligament, cortical and trabecular bones were loaded at 45º in relation to the tooth's long axis. All structures were considered to be homogeneous and isotropic, with the exception of the enamel (anisotropic). A standard masticatory force (100 N) was simulated on one of the models, while on the other one a maximum masticatory force was simulated (235.9 N). The software used were: PATRAN for pre- and post-processing and Nastran for processing. In the cementoenamel junction area, tensile forces reached 14.7 MPa in the 100 N model, and 40.2 MPa in the 235.9 N model, exceeding the enamel's tensile strength (16.7 MPa). The fact that the stress concentration in the amelodentinal junction exceeded the enamel's tensile strength under simulated conditions of maximum masticatory force suggests the possibility of the occurrence of non-carious cervical lesions such as abfractions.

Rare earth element and yttrium geochemistry applied to the genetic study of cryolite ore at the Pitinga Mine (Amazon, Brazil)

This work aims at the geochemical study of Pitinga cryolite mineralization through REE and Y analyses in disseminated and massive cryolite ore deposits, as well as in fluorite occurrences. REE signatures in fluorite and cryolite are similar to those in the Madeira albite granite. The highest ΣREE values are found in magmatic cryolite (677 to 1345 ppm); ΣREE is lower in massive cryolite. Average values for the different cryolite types are 10.3 ppm, 6.66 ppm and 8.38 ppm (for nucleated, caramel and white types, respectively). Disseminated fluorite displays higher ΣREE values (1708 and 1526ppm) than fluorite in late veins(34.81ppm). Yttrium concentration is higher in disseminated fluorite and in magmatic cryolite. The evolution of several parameters (REEtotal, LREE/HREE, Y) was followed throughout successive stages of evolution in albite granites and associated mineralization. At the end of the process, late cryolite was formed with low REEtotal content. REE data indicate that the MCD was formed by, and the disseminated ore enriched by (additional formation of hydrothermal disseminated cryolite), hydrothermal fluids, residual from albite granite. The presence of tetrads is poorly defined, although nucleated, caramel and white cryolite types show evidence for tetrad effect.

Unconstrained Finite Element for Geometrical Nonlinear Dynamics of Shells

This paper presents a positional FEM formulation to deal with geometrical nonlinear dynamics of shells. The main objective is to develop a new FEM methodology based on the minimum potential energy theorem written regarding nodal positions and generalized unconstrained vectors not displacements and rotations. These characteristics are the novelty of the present work and avoid the use of large rotation approximations. A nondimensional auxiliary coordinate system is created, and the change of configuration function is written following two independent mappings from which the strain energy function is derived. This methodology is called positional and, as far as the authors' knowledge goes, is a new procedure to approximated geometrical nonlinear structures. In this paper a proof for the linear and angular momentum conservation property of the Newmark beta algorithm is provided for total Lagrangian description. The proposed shell element is locking free for elastic stress-strain relations due to the presence of linear strain variation along the shell thickness. The curved, high-order element together with an implicit procedure to solve nonlinear equations guarantees precision in calculations. The momentum conserving, the locking free behavior, and the frame invariance of the adopted mapping are numerically confirmed by examples. Copyright (C) 2009 H. B. Coda and R. R. Paccola.

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]

Alpha element abundances and gradients in the Milky Way bulge from FLAMES-GIRAFFE spectra of 650 K giants

Aims. We present the analysis of the [alpha/Fe] abundance ratios for a large number of stars at several locations in the Milky Way bulge with the aim of constraining its formation scenario. Methods. We obtained FLAMES-GIRAFFE spectra (R = 22 500) at the ESO Very Large Telescope for 650 bulge red giant branch (RGB) stars and performed spectral synthesis to measure Mg, Ca, Ti, and Si abundances. This sample is composed of 474 giant stars observed in 3 fields along the minor axis of the Galactic bulge and at latitudes b = -4 degrees, b = -6 degrees, b = -12 degrees. Another 176 stars belong to a field containing the globular cluster NGC 6553, located at b = -3 degrees and 5 degrees away from the other three fields along the major axis. Stellar parameters and metallicities for these stars were presented in Zoccali et al. (2008, A&A, 486, 177). We have also re-derived stellar parameters and abundances for the sample of thick and thin disk red giants analyzed in Alves-Brito et al. (2010, A&A, 513, A35). Therefore using a homogeneous abundance database for the bulge, thick and thin disk, we have performed a differential analysis minimizing systematic errors, to compare the formation scenarios of these Galactic components. Results. Our results confirm, with large number statistics, the chemical similarity between the Galactic bulge and thick disk, which are both enhanced in alpha elements when compared to the thin disk. In the same context, we analyze [alpha/Fe] vs. [Fe/H] trends across different bulge regions. The most metal rich stars, showing low [alpha/Fe] ratios at b = -4 degrees disappear at higher Galactic latitudes in agreement with the observed metallicity gradient in the bulge. Metal-poor stars ([Fe/H] < -0.2) show a remarkable homogeneity at different bulge locations. Conclusions. We have obtained further constrains for the formation scenario of the Galactic bulge. A metal-poor component chemically indistinguishable from the thick disk hints for a fast and early formation for both the bulge and the thick disk. Such a component shows no variation, neither in abundances nor kinematics, among different bulge regions. A metal-rich component showing low [alpha/Fe] similar to those of the thin disk disappears at larger latitudes. This allows us to trace a component formed through fast early mergers (classical bulge) and a disk/bar component formed on a more extended timescale.