Dissertations / Theses on the topic '3D finite element'

This thesis project provides a method to simulate the internal mechanical properties of the undeformed human lower limb when an external force is applied. The value of the external force is determined by the deformed magnetic resonance imaging (MRI) images, which is converted into displacement instead. The purpose is to predict the area of the lower limb that is at most risk of wounds. The tissues of the human lower limb that are concerned in this study are: skin, fat, muscle, fascia, tibia bone, fibula bone and bone marrow. MRI images taken of an undeformed lower limb from experimental people is used to create a three-dimensional finite element model. During the simulation, the finite element model considers the nonlinear behaviors of individual soft tissues instead of lumping them together. Simulation results are presented as a curve of the external force and deformation of the different types of lower limb tissues and it shows the stress distribution in the lower limb.

This master's thesis assesses the possibility of using the nite element method to solve the electromagneticwave equation in a fusion plasma in 3D. In particular, the frequency is chosen to match that of ioncyclotron resonance heating in the fusion experiment JET. In this work, a brief introduction on fusion isgiven, followed by an explanation of the damping process in a plasma. A projection of the 3D wave eldonto a poloidal plane is compared to the 2D wave eld produced by the code FEMIC for validation ofthe developed 3D code. The comparison was done with good results.The power spectrum and coupling resistance per toroidal mode obtained from the 3D model arealso compared to the corresponding quantities obtained from an analytical slab model. Though somediscrepancies can be seen near the toroidal mode number n = 0 and for higher mode numbers (jnj > 70),the appearance of the power spectra are similar. The dierence near n = 0 is attributed to inducedcurrents in the reactor wall, whereas for higher mode numbers, the dierence is likely due to bad resolution.The induced currents in the wall causes singularities in the chosen model of the coupling resistance. Thisproduces unreliable predictions of the coupling resistance.
Denna masteruppsats utvärderar möjligheten att använda finita elementmetoden till att lösa den elektro-magnetiska vågekvationen i ett fusionsplasma i 3D. Speciellt väljs frekvensen för att matcha frekvensen för uppvärmning genom joncyklotronresonans i fusionsexperimentet JET. I detta arbete ges en översiktlig introduktion till fusion, åtföljd av en förklaring av dämpningsprocessen i ett plasma. En projektion av 3D-vågfältet på ett poloidalt plan jämförs med 2D-vågfältet producerat av 2D-koden FEMIC för att validera den utvecklade 3D-koden. Jämförelsen gjordes med gott resultat.| |Effektspektrumet och kopplingsresistansen per toroidal mod från 3D-modellen jämförs också med motsvarande storheter från en analytisk 1D-modell. Trots att vissa skillnader kan ses nära det toroidala modtalet n = 0 och för högre modtal ( n > 70), är utseendet på effektspektrumen lika. Skillnaden nära n = 0 tillskrivs de inducerade strömmarna i reaktorväggen, medan för högre modtal beror skillnaden troligen på dålig upplösning. De inducerade strömmarna i väggen ger upphov till singulariteter i den valda modellen för kopplingsresistansen. Det resulterar i otillförlitliga värden tansen.

Hybrid stepper motors are being applied to more and more industrial regionsdue to their low cost compared with servo motors and prominent performance.Many industrial applications require accurate and eective methods for predictinga motor’s performance at the design stage. The geometry of the motorsis complicated and the magnetic saturation eect is also serious, giving riseto the diculty of understanding the transient behavior of the motors. Furthermore,the drive circuit and control algorithm are more sophisticated thanthose of traditional AC or DC motors. Lastly, the losses of the motors createthe rising of temperature, while the thermal eect and dynamic performanceaect each other.All these factors can be solved by simulating a hybrid stepper motor witha model combining the eect of electromagnetic field, control algorithm, andmotor loss together. In this thesis, a three-dimension (3D) finite elementmodel is developed in the software Maxwell for studying motor characteristics.The electromagnetic field is analyzed in a static state. The simulatedback electromagnetic force is verified by experiments. The feasibility of fullstepcontrol algorithm is analyzed. The vector control algorithm is applied tothe model through co-simulation of Simulink and Maxwell in Simplorer. The3D model is proved to be unrealistic for co-simulation. In the end, this thesissummarizes the modeling experience and gives recommendations on thetransient simulation of the motor.
Hybridstegsmotorer appliceras i fler ochfler industriapplikationer tack vare deras låga kostnad och förbättrad prestanda jämfört med servomotorer. Många branschapplikationer kräver exakta och effektiva metoder för att förutsäga motorns prestanda redan i konstruktionsstadiet. Motorns geometri är komplicerad och den magnetiska mättnadseffekten är också betydande, vilket försvårar modelleringen. Dessutom är drivkretsen och styralgoritmen mer sofistikerad än den för traditionella växeleller likströmsmotorer. Vidare så resulterar motorns förluster i temperaturökningar vilka påverkar dynamiska.Alla dessa faktorer kan studeras genom att simulera hybrida stegmotorer med en modell som kombinerar effekten av elektromagnetiskt fält, kontrollalgoritm och motorförluster tillsammans. I detta examensarbete utvecklas en tredimensionell finit elementmodell i programvaran Maxwell för att studera motorns elektromagnetiska egenskaper. Det elektromagnetiska fältet analyseras i ett statiskt tillstånd. Den beräknade mot-EMK:n har verifieras genom experiment. Vektorkontrollalgoritmen tillämpas på modellen genom samsimulering i Simulink och Maxwell i Simplorer. Den tredimensionella modellen visade sig vara orealistisk för samsimulering. Till sist summeras uppnådaerfarenheter och rekommendationer för fortsatt arbete ges.

Full 360 degree 3D finite element models are the most complete representation of Oil Country Tubular Goods (OCTG) premium threaded connections. Full 3D models can represent helical threads and boundary conditions required to simulate make-up and service loading. A methodology is developed to create a 360 degree full 3D parametric finite element model with helical threads as an effective design and analysis tool. The approach is demonstrated with the creation of a metal-to-metal seal integral joint full 3D model with manufacturer supplied geometry. The premium connection is decomposed into smaller parts to generate parametric geometric features. A controlled parametric meshing scheme is developed to manage mesh density in contact regions to effectively represent the mechanics in regions of interest while minimizing total element count. The scripted parametric approach allows for efficient geometric and mesh updates. Several methods to reduce and manage model runtimes are presented. An elastic-plastic material model is created with material coupon tensile tests results. Digital Image Correlation (DIC) is used to measure full-field displacement and strain data on the surface of the box. Experimental set up and data processing procedures are discussed. Error metrics are developed to correlate the finite element model results with the DIC experimental data. The DIC make-up experimental results are used to reconcile the finite element model to develop a minimum error make-up model relative to the pin rotation. The friction coefficient is estimated and the make-up torque-turn behavior is verified. The calibrated 3D finite element model is validated with ISO_13769 load series B axial and internal pressure loading experimental DIC data. Metal-to-metal seal metrics of contact pressure and seal length are evaluated.
Master of Science

Thesis (M.S.)--Wichita State University, Dept. of Mechanical Engineering.
"May 2006." Title from PDF title page (viewed on October 29, 2006). Thesis adviser: Behnam Bahr. Includes bibliographic references (leaves 85-89).

This paper deals with the visualization of numerical results as a very convenient method to understand and evaluate a solution which has been calculated as a set of millions of numerical values. One of the central research fields of the Chemnitz SFB 393 is the analysis of parallel numerical algorithms for large systems of linear equations arising from differential equations (e.g. in solid and fluid mechanics). Solving large problems on massively parallel computers makes it more and more impossible to store numerical data from the distributed memory of the parallel computer to the disk for later postprocessing. However, the developer of algorithms is interested in an on-line response of his algorithms. Both visual and numerical response of the running program may be evaluated by the user for a decision how to switch or adjust interactively certain parameters that may influence the solution process. The paper gives a survey of current programmer and user interfaces that are used in our various 2D and 3D parallel finite element programs for the visualization of the solution.

COORDENAÇÃO DE APERFEIÇOAMENTO DO PESSOAL DE ENSINO SUPERIOR
TECNOLOGIA EM COMPUTAÇÃO GRÁFICA
Este trabalho apresenta uma metodologia para geração de malhas adaptativas de elementos finitos 2D e 3D usando modeladores geométricos com multi-regiões e superfícies paramétricas. A estratégia adaptativa adotada é fundamentada no refinamento independente das curvas, superfícies e sólidos. Inicialmente as curvas são refinadas, no seu espaço paramétrico, usando uma técnica de partição binária da curva (binary-tree). A discretização das curvas é usada como dado de entrada para o refinamento das superfícies. A discretização destas é realizada no seu espaço paramétrico e utiliza uma técnica de avanço de fronteira combinada com uma estrutura de dados do tipo quadtree para gerar uma malha não estruturada de superfície. Essas malhas de superfícies são usadas como dado de entrada para o refinamento dos domínios volumétricos. A discretização volumétrica combina uma estrutura de dados do tipo octree juntamente com a técnica de avanço de fronteira para gerar uma malha sólida não estruturada de elementos tetraédricos. As estruturas de dados auxiliares dos tipos binary-tree, quadtree e octree são utilizadas para armazenar os tamanhos característicos dos elementos gerados no refinamento das curvas, superfícies e regiões volumétricas. Estes tamanhos característicos são definidos pela estimativa de erro numérico associado à malha global do passo anterior do processo adaptativo. A estratégia adaptativa é implementada em dois modeladores: o MTOOL (2D) e o MG (3D), que são responsáveis pela criação de um modelo geométrico, podendo ter, multi-regiões, onde no caso 3D as curvas e superfícies são representadas por NURBS.
This work presents a methodology for adaptive generation of 2D and 3D finite-element meshes using geometric modeling with multi- regions and parametric surfaces. The adaptive strategy adopted in this methodology is based on independent refinements of curves, surfaces and solids. Initially, the model´s curves are refined using a binary-partition algorithm in parametric space. The discratizetion of these curves is used as input for the refinement of adjacent surfaces. Surface discretization is also performed in parametric space and employs a quadtree-based refinement coupled to an advancing-front technique for the generation of an unstructured triangulation. These surface meshes are used as input for the refinement adjacent volumetric domains. Volume discretization combines an octree refinement with an advancing-front technique to generate an unstructural mesh of tetrahedral elements. In all stages of the adaptive strategy, the refinement of curves, surface meshes and solid meshes is based on estimated numerical errors associated to the mesh of the previous step in the adaptive process. In addition, curve and surface refinement takes into account metric distortions between parametric and Cartesian spaces and high curvatures of the model´s geometric entities. The adaptive strategies are implemented in two different modelers: MTOOL (2D) and MG (3D), which are responsible for the creation of a geometric model with multi-regions, where for case 3D the curves and surfaces are represented by NURBS, and for the interactive and automatic finite-element mesh generation associated to surfaces and solid regions. Numerical examples of the simulation of engineering problems are presented in order to validate the methodology proposed in this work.

Le nouvel élément de poutre est une évolution d'un élément de Timoshenko poutre avec un nœud supplémentaire situé à mi-longueur. Ce nœud supplémentaire permet l'introduction de trois composantes supplémentaires de contrainte afin que la loi constitutionnelle 3D complète puisse être utilisée directement. L'élément proposé a été introduit dans un code d'éléments finis dans Matlab et une série d'exemples de linéaires/petites contraintes ont été réalisées et les résultats sont systématiquement comparés avec les valeurs correspondantes des simulations ABAQUS/Standard 3D. Ensuite, la deuxième étape consiste à introduire le comportement orthotrope et à effectuer la validation de déplacements larges / petites contraintes basés sur la formulation Lagrangienne mise à jour. Une série d'analyses numériques est réalisée qui montre que l'élément 3D amélioré fournit une excellente performance numérique. En effet, l'objectif final est d'utiliser les nouveaux éléments de poutre 3D pour modéliser des fils dans une préforme composite textile. A cet effet, la troisième étape consiste à introduire un comportement de contact et à effectuer la validation pour un nouveau contact entre 3D poutres à section rectangulaire. La formulation de contact est dérivée sur la base de formulation de pénalité et de formulation Lagrangian mise à jour utilisant des fonctions de forme physique avec l'effet de cisaillement inclus. Un algorithme de recherche de contact efficace, qui est nécessaire pour déterminer un ensemble actif pour le traitement de contribution de contact, est élaboré. Et une linéarisation constante de la contribution de contact est dérivée et exprimée sous forme de matrice appropriée, qui est facile à utiliser dans l'approximation FEM. Enfin, on présente quelques exemples numériques qui ne sont que des analyses qualitatives du contact et de la vérification de l'exactitude et de l'efficacité de l'élément de 3D poutre proposé
The new beam element is an evolution of a two nodes Timoshenko beam element with an extra node located at mid-length. That extra node allows the introduction of three extra strain components so that full 3D stress/strain constitutive relations can be used directly. The second step is to introduce the orthotropic behavior and carry out validation for large displacements/small strains based on Updated Lagrangian Formulation. A series of numerical analyses are carried out which shows that the enhanced 3D element provides an excellent numerical performance. Indeed, the final goal is to use the new 3D beam elements to model yarns in a textile composite preform. For this purpose, the third step is introducing contact behavior and carrying out validation for new 3D beam to beam contact with rectangular cross section. The contact formulation is derived on the basis of Penalty Formulation and Updated Lagrangian formulation using physical shape functions with shear effect included. An effective contact search algorithm is elaborated. And a consistent linearization of contact contribution is derived and expressed in suitable matrix form, which is easy to use in FEM approximation. Finally, some numerical examples are presented which are only qualitative analysis of contact and checking the correctness and the effectiveness of the proposed 3D beam element

In the present research work, a three-dimensional Cellular Automata Finite Element (CAFE) multi-scale model was developed to simulate, ductile fracture, cleavage and the ductile-brittle transition in a structural steel. For the simulation of the ductile-brittle fracture, at least two Cellular Automata arrays are needed, one to represent the ductile material properties and the other one to account for the brittle fracture process. The cell sizes in both arrays are independent of each other and of the finite element size. The cell sizes in each Cellular Automata array are related to the microstructural process of each fracture mechanism. The finite elements size is chosen to represent the macro strain gradients accurately. The model was implemented through the user define material behavior subroutine VUMAT in the finite element program ABAQUS Explicit Version 5.6. In the CAFE model, the material information is moved from the structural response of finite elements and stored in the appropriated number of Cellular Automata (CA) arrays. In the present CAFE model, the Rousselier ductile damage model was applied to each ductile cell. The critical value of the maximum principal stress was used to assess the failure of each brittle cell. In the brittle CA arrays, four different cleavage fracture nucleation micromechanisms, found experimentally at te.st temperatures down to -196øC in a ferritic-pearlitic Grade A ship plate steel were included in the model. This was done in order to simulate the real microfeatures nucleating cleavage in ferritic steels. In this model, the physical damage parameters of the ductile and brittle parts were calibrated separately. After calibration the CAFE model simulated the experimentally measured distribution of brittle microcracks generated in the notch region of blunt four point double-notch bend tests performed at test temperatures from 25øC to -196øC. The ductile part of the CAFE model was calibrated with the simulation of tensile and impact Charpy tests performed at room temperature. Subsequently the model was applied to simulate the ductile-brittle transition of Grade A ship plate steel. When numerical against experimental data was obtained, the parameters were considered true material model parameters of the steel under analysis.

Thesis: S.M., Massachusetts Institute of Technology, Department of Mechanical Engineering, May, 2020
Cataloged from the official PDF of thesis.
Includes bibliographical references (pages 65-68).
Projection Micro-Stereolithography (P[mu]SLA) is one of the most high-throughput additive manufacturing methods, yet preserving the high-resolution characteristic of light-based polymerization techniques. However, further improvement of fabrication speed and precision is usually hindered by the undesired adhesive forces at the curing interface, which is an inevitable consequence of in situ liquid-to-solid phase transition. To overcome this limitation, a bio-inspired super low adhesive interface has been proposed based on the observation of a slippery water layer on the peristome surface of pitcher plant. This hydrophobic layer provides an effective shield to solid adhesion due to its low adhesive energy, and attracting force between fabricated part and UV curing interface is significantly reduced. The introduction of this new lubrication layer not only remarkably improves the fabrication speed, but also increases the refilling rate of liquid pre-polymer resin.
This ultra-low adhesive interface shows promises for pushing the boundaries of continuous 3D printing into a realm of high-throughput additive manufacturing methods ready for industrial applications. In this thesis, I sought to provide a more comprehensive understanding of the solid-solid interaction at the curing interface of 3D photo-polymerizing systems. The state-of- the-art review of current literature suggested that a surface-based cohesive contact theory from a continuum mechanics perspective was the most appropriate model to establish a connection between interfacial material properties and macroscopic measurement results from experiment. Based on that I analyzed the entire mechanical separation process using finite-element method, and provided a semi-quantitative explanation of the stability of such lubricant-infused nano-cavities against peeling forces.
This research lays the ground for elucidating the physical mechanism behind the general adhesion-separation problem, and framework has been constructed in a more general form to allow for analyzing a wide range of interdisciplinary problems involving the dynamics of anisotropic moving contact lines and the propagation of surface instabilities induced by adhesive contact.
by Huifeng Du.
S.M.
S.M. Massachusetts Institute of Technology, Department of Mechanical Engineering

Ventricle fibrillation (VF) is a condition in which the heart's lower chambers show an unsynchronized and chaotic motion which prevents the heart from pumping blood and oxygen to the body. VF is considered a sudden cardiac arrest and it is responsible for 300,000 sudden deaths in the USA yearly. The most effective way of reversing this life threatening condition is to apply an electrical shock directly to the heart using an Implantable Cardioverter-Defibrillator (ICD). The main issue in using ICDs is the placement of the defibrillating electrodes so that the current can be optimally channeled through the cardiac muscle, particularly in the left ventricular myocardium. According to the critical mass hypothesis, defibrillation will be successful when 75% of the myocardium tissue is halted by the defibrillation shock. The defibrillation threshold (DFT) or the minimum effective voltage required for successful results is suggested to be related to the myocardial voltage gradient (VG) distribution, but it has not been quantified. Moreover, the goal is to keep the DFT as low as possible to try to maximize the success of defibrillation, minimize the chance of myocardium damage and cardiac arrhythmias caused by high-intensity shocks, and also potentially reduce the battery size and as well as prolong the device's useful lifespan. Various numerical techniques have been used to model the heart to solve the governing equations required to obtain the myocardium VG distribution during electrical defibrillation. The Finite Element Method (FEM) has been of particular interest since it can handle the irregular domains, material inhomogeneities, and complex boundary conditions of problems in bioelectricity. In this thesis, a finite element model of the heart tissue is constructed in order to study and optimize the defibrillation mechanism. The modelling process starts with a surface reconstruction based on radial basis function interpolation to generate the triangular surface me
La fibrillation ventriculaire (FV) est un état dans lequel la cavité inférieure du coeur montre une motion asynchrone et chaotique, empêchant le coeur de pomper le sang et l'oxygène au corps. La FV est considérée comme un arrêt cardiaque soudain, responsable de la mort subite de 300,000 personnes chaque année aux Etats-Unis. Afin d'inverser cette condition mortelle, le recours le plus efficace est la délivrance d'un choc électrique directement au niveau du coeur à l'aide d'un Défibrillateur Cardioverteur Implantable (DCI). Le principal problème de l'utilisation des DCIs est le placement des électrodes défibrillateurs pour permettre au courant d'être conduit optimalment à travers du muscle cardiaque, en particulier, le myocarde ventriculaire gauche. Selon l'hypothèse de la masse critique, la défibrillation sera réussi quand 75% du tissu de myocarde est inactivé par le choc de défibrillation. Le seuil de défibrillation ou la tension efficace minimale exigée pour donner des résultats réussis est suggéré d'être liés à la distribution du gradient de la tension myocardique, toutefois, il n'a pas été mesurée. En outre, le but est de maintenir le seuil de défibrillation aussi bas que possible pour essayer de maximiser le succès de la défibrillation, réduire au minimum le risque des dommages de myocarde et des arythmies cardiaques provoqués par des chocs à haute intensité, et aussi permettre de réduire la taille de la batterie ainsi que prolonger la durée de la vie utile du dispositif. De diverses techniques numériques ont été utilisées pour modeler le coeur afin de résoudre les équations régissant requises pour obtenir la distribution de gradient de la tension myocardique pendant la défibrillation. La méthode des éléments finis (FEM) a été d'intérêt particulier car elle peut gérer les domaines irréguliers, les hétérogénéité de matériel, et les conditions aux limites complexe de problèmes dans la bioélec

Often bone fractures are joined by inserting metal plates and screws to hold the fragmented bone under compression. However, after the fractured bone is healed removing the screws leaves holes in the bone which takes months to fill up and heal completely. The goal of this research is to investigate those voids specifically in a finite element model of a femur. The holes were found to experience high stress that can easily lead to crack propagations during everyday activities. Finite element models of femurs were modeled after two common fracture fixation systems, specifically just after the plates, rods and screws are removed. To observe the stress levels bones are likely to experience, common mechanical tests that are relevant to or associated with common daily activities were performed. While the 3-point bending tests did not yield significant results, the compression and torsion tests produced high stress areas near the screw holes. In certain cases, the von Mises’ stress reached 3.66 x 106 N/mm2. Our finite element modeling seeks to establish groundwork for future explorations on the holes created by fracture fixation hardware. In the future, this work will lead to redesigning of fixation systems with reduced stress concentration around the holes. Therefore, the initiation of new cracks around these holes will be limited during everyday activity.

The influence of the surface excavation and loading on the existing tunnels has been investigated using a Finite Element Method program, Plaxis 3D Tunnel. A parametric study has been carried out where the parameters were the depth of the surface excavation, the eccentricity of the excavation with respect to the tunnel axis, the height of the embankment fill and the stiffness of the soil. It is found that, excavations over the existing tunnels have a negative effect on the tunnel lining capacity since the unloading due the excavation reduces the normal forces and increases the bending moments. On the other hand, it is found that surface loading within the limits considered in this study is not critical in terms of the tunnel stability due to the increase in normal forces and decrease in bending moments.

The J-integral has been used over the past twenty years to predict crack initiation and crack growth in elastic-plastic materials. Fundamental to the J-integral gaining acceptance as a suitable fracture parameter for elastic-plastic solids were the papers of Rice and Rosengren, and Hutchinson. In 1968 Rice and Rosengren, and Hutchinson both presented theoretical models for J scaling of the stresses and strains around a stationary crack tip, under plane stress and plane constraint. However, the theories of Rice and Rosengren, and Hutchinson are unable to take into account any thickness effect, and other workers have investigated the limitations on the validity of J controlled fracture, culminating in empirical size requirements for test specimens and structural components. Druce and co-workers at U.K.A.E.A., Harwell carried out a comprehensive range of tests to investigate the size requirements for a C-Mn steel and found significant variations from currently accepted theory. Druce found, that while the crack initiation of J was fairly independent of specimen thickness, the slope of the J crack growth resistance curve decreased dramatically with increasing thickness. Thus there is much less potential for stable crack growth in thicker sections of the same material. This clearly contradicts the simplistic engineering approach, where an increased factor of structural safety may be obtained by increasing the dimensions of sections, which make up the structure. The main purpose of this thesis is to describe a finite analysis of the specimens tested by Druce, in order to explain the discrepancies from accepted theory. In addition, the numerical results were also used to verify theoretical models of crack tip stress and strain fields, and the assumptions used in J controlled crack growth.

Thesis (M.S.)--West Virginia University, 1999.
Title from document title page. Document formatted into pages; contains ix, 166 p. : ill. (some col.). Vita. Includes abstract. Includes bibliographical references (p. 149-158).

Made available in DSpace on 2014-06-11T19:28:58Z (GMT). No. of bitstreams: 0 Previous issue date: 2008-08-21Bitstream added on 2014-06-13T19:58:24Z : No. of bitstreams: 1 almeida_eo_me_araca.pdf: 3755451 bytes, checksum: 666f30ecb05cff86a47b4d7dfe771edb (MD5)
A barra pré-fabricada é uma opção de tratamento vantajosa por minimizar o tempo de confecção da prótese implantossuportada do tipo protocolo de carga imediata com níveis de adaptação e passividade previsíveis. No entanto, para a sua aplicação, o arco do paciente deve adaptar-se anatomicamente ao formato da barra pré-fabricada, e a qualidade óssea deve ser compatível com a indicação da carga imediata. Além disso, a literatura não apresenta dados numéricos sobre a influência de diferentes formatos horizontais do arco mandibular na distribuição das tensões no osso, assim como não é homogênea a definição do critério de análise a ser utilizado quando materiais frágeis, como o osso, são estudados utilizando o método dos elementos finitos. Devido a isso, foi realizada uma revisão da literatura com o objetivo de relacionar as propriedades dos materiais estudados com o critério de análise adotado. Para isto, foram levantados 2061 artigos através do indexador Medline database, durante os anos de 2004 a 2008. Observou-se que tanto o critério das tensões de von Mises como a máxima tensão principal foram utilizados para o estudo dos materiais considerados frágeis ou não. Para o estudo da influência de diferentes formatos horizontais do arco mandibular e diferentes tipos de osso na distribuição das tensões no osso cortical e medular na simulação de uma prótese fixa implantossuportada do tipo protocolo mandibular confeccionada através do sistema pré-fabricado Neopronto (Neodent, Implante Osseointegrável), quatro modelos (M) representativos...
Prefabricated bar is a profitable treatment option for minimizing the preparation time of fixed protocol-type prosthesis in immediate load with adaptation levels and predictable passivity. Nevertheless, in order to be applied, the patient's arch must be anatomically adapted to the format of the prefabricated bar, and bone quality must be compatible to immediate load indication. However, literature does not show any numerical data about the influence of different horizontal formats of the mandibular arch on bone stress distribution, as well as it is not homogeneous which analysis criterion is to be used when friable materials, such as the bone, are studied by making use of the finite element method. Due to this, a literature review was elaborated with the purpose to relate the materials' properties through the adopted analysis criterion. Therefore, 2061 articles were raised through the Medline Database Indexer from 2004 to 2008. It was noticed that not only Von Mises's equivalent stress criterion but also maximum principal stress were used for the study of materials considered as friable or not. For the study of the influence of different horizontal formats of mandibular archs and different types of bone in stress distribution on the cortical and medular bone in simulating a fixed implant-supported protocol-type mandibular prosthesis prepared through the Neopronto prefabricated system (Neodent, Osseointegratable Implant), four representative models (M) of a totally endentulate mandibular arch restored by the Neopronto prefabricated system (Neodent, Osseointegratable Implant) over four interforaminal implants were prepared through the SolidWorks 2007 program with a bone range (I to IV) (M.I - M.II - M.III - M.IV), considering the horizontal dimension as of regular size... (Complete abstract click electronic access below)

This dissertation developed a method that can accurately and efficiently capture the response of a structure by rigorous combination of a reduced-dimensional beam finite element model with a model based on full two-dimensional (2D) or three-dimensional (3D) finite elements. As a proof of concept, a joint 2D-beam approach is studied for planar-inplane deformation of strip-beams. This approach is developed for obtaining understanding needed to do the joint 3D-beam model. A Matlab code is developed to solve achieve this 2D-beam approach. For joint 2D-beam approach, the static response of a basic 2D-beam model is studied. The whole beam structure is divided into two parts. The root part where the boundary condition is applied is constructed as a 2D model. The free end part is constructed as a beam model. To assemble the two different dimensional model, a transformation matrix is used to achieve deflection continuity or load continuity at the interface. After the transformation matrix from deflection continuity or from load continuity is obtained, the 2D part and the beam part can be assembled together and solved as one linear system. For a joint 3D-beam approach, the static and dynamic response of a basic 3D-beam model is studied. A Fortran program is developed to achieve this 3D-beam approach. For the uniform beam constrained at the root end, similar to the joint 2D-beam analysis, the whole beam structure is divided into two parts. The root part where the boundary condition is applied is constructed as a 3D model. The free end part is constructed as a beam model. To assemble the two different dimensional models, the approach of load continuity at the interface is used to combine the 3D model with beam model. The load continuity at the interface is achieved by stress recovery using the variational-asymptotic method. The beam properties and warping functions required for stress recovery are obtained from VABS constitutive analysis. After the transformation matrix from load continuity is obtained, the 3D part and the beam part can be assembled together and solved as one linear system. For a non-uniform beam example, the whole structure is divided into several parts, where the root end and the non-uniform parts are constructed as 3D models and the uniform parts are constructed as beams. At all the interfaces, the load continuity is used to connect 3D model with beam model. Stress recovery using the variational-asymptotic method is used to achieve the load continuity at all interfaces. For each interface, there is a transformation matrix from load continuity. After we have all the transformation matrices, the 3D parts and the beam parts are assembled together and solved as one linear system.

La meccanica del contatto fra ruota ferroviaria e rotaia, è una delle più importanti aree di studio nell’Ingegneria Ferroviaria. Ad oggi, un vasto numero di formulazioni analitiche sono state proposte dai ricercatori, con lo scopo di valutare i parametri di contatto tra ruota e rotaia ed arrivare ad una descrizione affidabile delle forze che agiscono nell’area di contatto. Comunque, solo alcuni dei metodi disponibili permettono di considerare, nelle loro elaborazioni, la reale geometria di contatto tra ruota e rotaia o la non linearità delle proprietà dei materiali impiegati. Un’alternativa ai metodi analitici per descrivere la fisica del contatto tra ruota e rotaia e di fare uso di tecniche computazionali numeriche quale l’Analisi agli Elementi Finiti (FEA – Finite Element Analysis). Il vantaggio di questo metodo è dato dalla possibilità di poter modellare complesse geometrie che permettono di simulare più accuratamente il contatto tra ruota e rotaia e determinare l’ampiezza degli stress e la loro distribuzione, oltre alla dimensione e forma dell’area di contatto. Per questa ragione, un modello agli Elementi Finiti tridimensionale della ruota ferroviaria e della rotaia è stato creato utilizzando Ansys Parametric Design Language di ANSYS per studiare il problema di contatto normale. Per verificare l’accuratezza dei risultati forniti dal modello, questo è stato validato in confronto alla teoria di Hertz sul contatto elastico. Questa teoria rappresenta la base sulla quale la maggior parte dei modelli computazionali vengono sviluppati. Il principali scopi di questa tesi sono quindi di studiare i problemi del contatto tra ruota e rotaia e di valutare l’influenza dei parametri operativi quali il coefficiente di frizione, lo spostamento laterale della sala montata e l’inclinazione della rotaia sul piano orizzontale, sulla fisica del contatto.

The main purpose of this thesis is to develop a 3D beam element in order to model wind turbine wings made of composite materials. The proposed element is partly based on the formulation of the classical beam element of constant cross-section without shear deformation (Euler-Bernoulli) and including Saint-Venant torsional effects for isotropic materials, similarly to the one presented in Batoz & Dhatt (1990, pp.147-190). The main novelty consists in the addition of the coupling between axial and bending with torsional effects that may arise when using composite materials. PreComp, a free access code developed by the National Renewable Energy Laboratory (NREL) to provide structural properties for composite blades, is used to obtain the section properties for the beam element. Its performance is assessed, showing its inaccuracy especially when calculating torsional related constants when webs are present in the cross-section. Shell models of constant cross-section cantilever blades are developed to assess the performance of the beam elements, including or not coupling terms. Natural frequencies and displacements under static loads are compared for different study cases of increasing complexity. For fiber-reinforced materials, elements with coupling terms show good agreement with the shell model, especially for the dynamic problem. Elements without coupling terms are unable to capture the dynamic behavior, as these terms seem to have a higher effect on the results when compared to the static case (especially the FT term).

PONTIFÍCIA UNIVERSIDADE CATÓLICA DO RIO DE JANEIRO
COORDENAÇÃO DE APERFEIÇOAMENTO DO PESSOAL DE ENSINO SUPERIOR
FUNDAÇÃO DE APOIO À PESQUISA DO ESTADO DO RIO DE JANEIRO
PROGRAMA DE SUPORTE À PÓS-GRADUAÇÃO DE INSTS. DE ENSINO
PROGRAMA DE EXCELENCIA ACADEMICA
BOLSA NOTA 10
A grande maioria das obras geotécnicas, tais como algumas fundações, aterros e escavações, apresentam uma configuração geométrica tipicamente tridimensional. No entanto, em função da complexidade da obtenção de soluções considerando os aspectos inerentes à condição de deformação e tensão tridimensionais, as análises destas obras têm sido comumente realizadas adotando-se as aproximações de deformação plana e/ou axissimétrica. Esta dissertação apresenta um modelo computacional baseado na formulação em deslocamento do método dos elementos finitos para a simulação de problemas mecânicos de equilíbrio estático de obras geotécnicas em condição tridimensional. São ressaltados os problemas de capacidade de carga de fundações superficiais e a simulação da construção de aterros. O solo é considerado como um material não linear elasto-plástico e os modelos não associados Mohr-Coulomb modificado e Lade-Kim são adotados para representar sua relação tensão-deformação resistência. As estratégias de solução não linear, tanto a nível global quanto em nível de ponto de Gauss, são apresentadas e discutidas. Os exemplos de capacidade de carga das fundações superficiais são comparados com os resultados da teoria do equilíbrio limite e da análise limite indicando bons resultados. O efeito da tensão intermediária na trajetória de tensão é apresentado e discutido nos exemplos de simulação da construção de aterros.
The vast majority of geotechnical works, such as some foundations, embankments and excavations, show a typical three-dimensional geometric configuration. However, depending on the complexity of obtaining solutions considering the aspects inherent to the three-dimensional deformation and stress, the analysis of these works have been performed by adopting common approaches and plane strain or axisymmetric. This dissertation presents a computational model based on displacement formulation of finite element method for simulation of mechanical problems of static equilibrium condition of geotechnical works in three dimensions. It highlights the problems of load capacity of shallow foundations and the simulation of embankment construction. The soil is considered as material non-linear elastic-plastic models and non-associated Mohr- Coulomb and modified Lade-Kim are adopted to represent their relationship stress-strain-resistance. The nonlinear solution strategies, both overall and in the Gauss point level, are presented and discussed. Examples of load capacity of shallow foundations are compared with the results of the theory of limit equilibrium and limit analysis showing good results. The effect of tension in the intermediate stress path is presented and discussed in the simulation examples of the construction of embankments.

Thesis (MScEng)--Stellenbosch University, 2011.
ENGLISH ABSTRACT: Several transverse flux and longitudinal flux linear generator topologies exist for freepiston Stirling engine applications. In this thesis the transverse flux permanent magnet linear generators are investigated together with a back-to-back converter which can deliver the electrical energy from the linear generator to the electrical network. The transverse flux permanent magnet linear generator is geometrically optimised with the aim to maximise the power-to-weight ratio while maintaining preset power and efficiency levels. An optimised 3 kW linear generator is built and the measured results correlate to the simulation results. A close-loop current control scheme is introduced to control the current of the rectifier, which is part of the back to back converter. The transverse flux permanent magnet linear generator is connected to the input of the rectifier which has the ability to force a specific current from the generator. The measured results of the rectifier correlate to the results of the simulations that were done. The current control present some complications and it is suggested that another control scheme is used. A close-loop voltage control scheme is introduced for the control of the DC bus voltage. The DC bus is connected between the rectifier and the inverter, which is the other part of the back-to-back converter. A close-loop current control scheme is introduced to control the inverter current that flows from the inverter to the electrical network. The measured results of the inverter and the DC bus correlate to the results of the simulations that were done. The results of the system, including the generator, rectifier and inverter, tested as a unit is presented and discussed.
AFRIKAANSE OPSOMMING: Verskeie tranverse vloed en longitudinale vloed lineˆere generator topologie¨e bestaan vir vrysuier Stirling enjin toepassings. In hierdie tesis word ’n transverse vloed permanente magneet lineˆere generator ondersoek saam met ’n omsetter. Die omsetter dra die elektriese energie van die generator oor aan die elektriese netwerk. Die transverse vloed permanente magneet lineˆere generator word geometries geoptimeer met die doel om die drywing-tot-gewig verhouding te maksimiseer terwyl vasgestelde drywing en effektiwiteit vlakke behou word. ’n Geoptimeerde 3kW lineˆere generator prototipe is vervaardig en die gemete resultate is geverifieer met die simulasie resultate. ’n Geslote lus stroombeheer strategie word voorgestel om die stroom te beheer van die gelykrigter, wat deel is van die omsetter. Die transverse vloed permanente magneet lineˆere generator word aan die gelykrigter, wat die vermo¨e het om ’n spesifieke stroom uit die generator te forseer, se intree verbind. Die gemete resultate van die gelykrigter wat gebou is stem goed ooreen met die van die simulasies wat gedoen is. Die stroombeheer hou komplikasies in wat bespreek word. Dus word die gebruik van ’n alternatiewe stroombeheer voorgestel. ’n Geslote lus spannings beheer strategie¨e word voorgestel om die gs. busspanning te beheer. Die gs. bus is gekonnekteer tussen die gelykrigter en die wisselrigter, wat ook deel uitmaak van die omsetter. ’n Geslote lus stroom beheer word voorgestel om die stroom te beheer wat vanaf die wisselrigter na die elektriese netwerk toe vloei. Die gemete resultate van die wisselrigter en die gs. bus stem goed ooreen met die van die simulasies wat gedoen is. Die resultate van die hele stelsel, wat die generator, gelykrigter en die wisselrigter insluit, wat as ’n eenheid getoets is word weergegee en bespreek.

Doctor of Philosophy
Department of Mechanical and Nuclear Engineering
Youqi Wang
Textile fabrics and textile composite materials demonstrate exceptional mechanical properties, including high stiffness, high strength to weight ratio, damage tolerance, chemical resistance, high temperature tolerance and low thermal expansion. Recent advances in weaving techniques have caused various textile fabrics to gain applications in high performance products, such as aircrafts frames, aircrafts engine blades, ballistic panels, helmets, aerospace components, racing car bodies, net-shape joints and blood vessels. Fabric mechanical properties are determined by fabric internal architectures and fabric micro-geometries are determined by the textile manufacturing process. As the need for high performance textile materials increases, textile preforms with improved thickness and more complex structures are designed and manufactured. Therefore, the study of textile fabrics requires a reliable and efficient CAD/CAM tool that models fabric micro-geometry through computer simulation and links the manufacturing process with fabric micro-geometry, mechanical properties and weavability. Dynamic Weaving Process Simulation is developed to simulate the entire textile process. It employs the digital element approach to simulate weaving actions, reed motion, boundary tension and fiber-to-fiber contact and friction. Dynamic Weaving Process Simulation models a Jacquard loom machine, in which the weaving process primarily consists of four steps: weft insertion, beating up, weaving and taking up. Dynamic Weaving Process Simulation simulates these steps according to the underlying loom kinematics and kinetics. First, a weft yarn moves to the fell position under displacement constraints, followed by a beating-up action performed by reed elements. Warp yarns then change positions according to the yarn interlacing pattern defined by a weaving matrix, and taking-up action is simulated to collect woven fabric for continuous weaving process simulation. A Jacquard loom machine individually controls each warp yarn for maximum flexibility of warp motion, managed by the weaving matrix in simulation. Constant boundary tension is implemented to simulate the spring at each warp end. In addition, process simulation adopts re-mesh function to store woven fabric and add new weft yarns for continuous weaving simulation. Dynamic Weaving Process Simulation fully models loom kinetics and kinematics involved in the weaving process. However, the step-by-step simulation of the 3D weaving process requires additional calculation time and computer resource. In order to promote simulation efficiency, enable finer yarn discretization and improve accuracy of fabric micro geometry, parallel computing is implemented in this research and efficiency promotion is presented in this dissertation. The Dynamic Weaving Process Simulation model links fabric micro-geometry with the manufacturing process, allowing determination of weavability of specific weaving pattern and process design. Effects of various weaving process parameters on fabric micro-geometry, fabric mechanical properties and weavability can be investigated with the simulation method.

Essential modern transportation systems together with the high demand for sustainable pavements under applied vehicular loading have led to a great deal of research worldwide of concrete pavements. Despite progressive knowledge of concrete pavement behaviour under applied loads, concrete pavements are still subject to deterioration due to crack initiation and propagation, indicating the need for further research. Cracks can be related to fatigue of concrete or erosion of materials in sub-layers. Transverse joints in concrete pavements are the locations where most pavement distress appears, leading to deterioration of the riding quality and featuring high maintenance cost. The state of stresses in the concrete surrounding dowel bars, in dowel jointed concrete pavements, are a major factor that contribute to transverse joint distress. A three-dimensional (3D) finite element model is developed in this study for analysing a dowel-jointed concrete pavement. The effects of different pavement and joint related parameters on the load transfer characteristics of a joint have been evaluated using the 3D finite element model. The numerical results from FE modelling are validated with classical analytical solutions of shear and moment along the dowel. Five loading cases are applied in the model to replicate realistic vehicular loadings approaching and leaving the joint. Group action of the dowel bar system has also been examined.
Thesis (PhD Doctorate)
Doctor of Philosophy (PhD)
Griffith School of Engineering
Science, Environment, Engineering and Technology
Full Text

This thesis is aimed on three-dimensional finite element time domain analysis. Vector finite elements are used for analysis of cavity resonator made from PEC and filled with vacuum. Results are compared with classical “nodal” finite elements and with published literature.

Lors des accidents de la route, les lésions de l’abdomen mettent des vies en jeu et nécessitent un traitement médical long et coûteux. Il est important d’offrir une meilleure prévention des lésions abdominales traumatiques. Des modèles numériques de corps humain ont été développés afin de comprendre les mécanismes lésionnels. L’objectif de cette thèse est de contribuer au développement d’un modèle numérique de foie humain pour la prédiction de lésions hépatiques en cas de choc.Dans un premier temps, une caractérisation morphologique du foie a été réalisée sur 78 scanners et a mis en évidence l’existence de 4 morphotypes. Ensuite, des expérimentations ont mis en évidence l’influence de la pressurisation des vaisseaux sur la déformation de la capsule et ont permis d’identifier la déformation à rupture locale lors d’un choc.Enfin, les modèles numériques des 4 morphotypes identifiés ont été construits pour simuler les essais de décélération réalisés dans la partie expérimentale
In road accidents, lesions of the abdomen are life-threatening, and require a long and expensive medical treatment. It is important to offer a better prevention of traumatic abdominal injuries. Numerical models of the human body have been developed to understand the mechanisms of injury. The aim of this thesis is to contribute to the development of a numerical model of human liver for the prediction of liver lesions in case of shock.First, a morphological characterization of the liver was performed on 78 healthy patient scanners and highlight the existence of 4 morphotypes.Then, experiments were carried out to highlight the influence of this pressurization of the vessels on the strain state of the Glisson capsule and to identify the local ultimate strain during an impact.Finally, the numerical models of the 4 identified morphotypes were constructed and used to simulate the deceleration tests performed in the experimental part

A new three-dimensional finite solid element (3D FE) model for dowel-type wood connections was developed using the concept of a beam on a nonlinear wood foundation, which addresses the intricate wood crushing behaviour under the connector in a dowel type connection. In order to implement the concept of wood foundation with solid elements, a 3D FE wood foundation model was defined within a prescribed foundation zone surrounding the dowel. Based on anisotropic plasticity material theory, the material model for the foundation zone was developed using effective foundation material constants that were defined from dowel-embedment test data. New 3D FE single nail connection models were developed that incorporated the wood foundation model. The 3D wood foundation model was justified and validated using dowel-embedment test data with a range of dowel diameters, from a 2.5-mm nail to a25.4-mm bolt. The connection models provided successful results in simulating the characteristics of load-slip behaviour that were experimentally observed. Based on the success of the single nail connection models, several applications of the3D FE connection models were investigated including statistical wood material models, bolted connection models and a multiple nail connection model. Throughout the application studies, discussion of the benefits and limitations of the new model approach using the 3D FE wood foundation are presented. Also, future areas of study are proposed in order to improve the 3D FE dowel-type wood connections models.

Finite Element Analysis (FEA) is proposed to simulate the connection response of practical in-service conditions and test the performance of Oil Country Tubular Goods (OCTG) premium threaded connections. A plug-in is developed in Abaqus/CAE for creating the 360-degree full 3D parametric finite element model with helical threads as an effective design and analysis tool. All size, position and material data of the model are parameterized. The premium connection plug-in accepts input from the Graphical User Interface (GUI) for further modification. Each premium connection component is programed as a collection of single-purpose independent functions organized as an independent module in order to allow users to modify specific function behavior conveniently. A main program is designed as an Abaqus kernel plug-in to achieve all the functions by calling these independent functions, making the plug-in flexible. Each single script file is not too long to jeopardize readability. The GUI of the plug-in is designed with proper layout arrangement and illustrations to make the plug-in user-friendly and easy to use. The premium connection FE model is used in a virtual test to validate the model against the ISO 13679 test protocol and is used to develop the seal metrics for points on the ISO 13679 sealability envelope. The plug-in can be used to develop and evaluate the design envelope of the premium connection.
Master of Science

In recent decades, railway transport infrastructures have been regaining their importance due to their efficiency and environmentally friendly technologies. This has led to increasing train speeds, higher axle loads and more frequent train usage. These improved service provisions have however brought new challenges to traditional railway track engineering, especially to track geotechnical dynamics. These challenges demanded for a better understanding of the track dynamics. Due to the large cost and available load conditions limitation, experimental investigation is not always the best choice for the dynamic effect study of railway track structure. Comparatively speaking, an accurate mathematical modeling and numerical solution of the dynamic interaction of the track structural components reveals distinct advantage for understanding the response behavior of the track structure. The purpose of this thesis is to study the influence of design parameters on dynamic response of the railway track structure by implementing Finite Element Method (FEM). According to the complexity, different railway track systems have been simulated, including: Beam on discrete support model, Discretely support track including ballast mass model and Rail on sleeper on continuum model. The rail and sleeper have been modeled by Euler-Bernoulli beam element. Spring and dashpot has been used for the simulation of railpads and the connection between the sleeper and ballast ground. Track components have been studied separately and comparisons have been made between different models. The finite element analysis is divided into three categories: eigenvalue analysis, dynamic analysis and general static analysis. The eigenfrequencies and corresponding vibration modes were extracted from all the models. The main part of the finite element modeling involves the steady-state dynamic analysis, in which receptance functions were obtained and used as the criterion for evaluating the dynamic properties of track components. Dynamic explicit analysis has been used for the simulation of a moving load, and the train speed effect has been studied. The displacement of the trackbed has been evaluated and compared to the measurement taken in Sweden in the static analysis.

Made available in DSpace on 2014-06-11T19:35:03Z (GMT). No. of bitstreams: 0 Previous issue date: 2009-09-17Bitstream added on 2014-06-13T19:44:52Z : No. of bitstreams: 1 gomes_ea_dr_araca.pdf: 12533208 bytes, checksum: 62db616b630b0f9a9807e6b24297d9e8 (MD5)
Objetivo: A análise pelo método dos elementos finitos tridimensional foi utilizada para avaliar o efeito dos desajustes angular e vertical simulado em próteses fixas implantossuportadas parafusadas de 3 elementos, na resposta biomecânica do tecido ósseo peri-implantar, implantes e componentes protéticos. Material e método: Foram confeccionados 4 modelos tridimensionais, representativos de secção mandibular direita contendo um implante na região de segundo pré-molar (2oPM) e outro na região de segundo molar (2oM) ferulizados por uma prótese metalocerâmica implantossuportada parafusada de 3 elementos, diferenciados pelo tipo de desajuste protético, gerando os grupos: GC (grupo controle) - prótese totalmente adaptada aos implantes, DAU (desajuste angular unilateral) - prótese com desajuste angular unilateral de 100μm na região mesial do 2oM, DVU (desajuste vertical unilateral) - prótese com desajuste vertical unilateral de 100μm na região mesial do 2oM e DVT (desajuste vertical total) - prótese com desajuste vertical de 100μm em toda a plataforma de assentamento da infraestrutura no 2oM. Com a utilização do programa de elementos finitos Ansys® uma força vertical total de 400N foi distribuída em 12 pontos de parada cêntrica, sendo que cada molar recebeu um carregamento total de 150N e o 2oPM de 100N, gerando mapas de tensões. Resultados: Verificou-se que a distribuição e os valores de tensões para o tecido ósseo peri-implantar foram similares para todos os grupos. Os grupos com desajustes alteraram o padrão de distribuição e elevaram a magnitude das tensões em relação ao grupo GC. No grupo DAU os maiores valores de tensões foram observados no corpo do implante e parafuso de retenção. Já nos grupos DVU e DVT elevados valores de tensões foram encontrados na plataforma de assentamento da infraestrutura...
Purpose: The three-dimensional finite element analysis was used to evaluate the effect of vertical and angular misfit in 3-piece implant-supported and screw-retained fixed prostheses on biomechanical response in peri-implant bone tissue, implants and prosthetic components. Material and method: Four three-dimensional models were fabricated to represent a right posterior mandibular section with one implant in the region of second premolar (2°PM) and other in the region of second molar (2°M). The implants were splinted by a 3- piece implant-supported metal ceramic prosthesis and differed according to the type of misfit as represented by the groups: CG (control group) - prosthesis with complete fit to the implants, UAM (unilateral angular misfit) – prosthesis presenting unilateral angular misfit of 100μm in the mesial region of the 2°M, UVM (unilateral vertical misfit) – prosthesis presenting unilateral vertical misfit of 100μm in the mesial region of the 2°M, and TVM (total vertical misfit) – prosthesis presenting total vertical misfit of 100μm in the platform of the framework in the 2°M. A vertical load of 400N was distributed and applied on 12 centric points by the software Ansys®. So, vertical loading of 150N was applied on each molar while vertical load of 100N was applied on the 2°PM. Stress maps were obtained by processing of the finite element software. Results: The stress values and distribution on peri-implant bone tissue were similar for all groups. The groups presenting misfit exhibited different distribution pattern and increased stress magnitude in comparison to the group CG. The highest stress values in group UAM were observed in the implant body and retention screw. The groups UVM and TVM exhibited high stress values in the platform of the framework... (Complete abstract click electronic access below)

Peat deposits are distributed mostly along the coast of West Kalimantan, Indonesia, particularly in and around the provincial capital of Pontianak. There are many problems with constructing over peat soil as the existence of this type of soil always generates geotechnical engineering problems for regional development. For light construction on peat land, several variations of Indonesian traditional floating wood foundations, commonly called tiang tongkat foundations, are still being used today. A tiang tongkat foundation of any dimension is constructed over different fields. The foundation was modelled as three-dimensional linear elastic and the Pontianak soft organic soil was modelled as undrained Soft-Soil-Creep Model. All of the 324 models were made to be used for simulation by means of the Plaxis 3-Dimensional Foundation Program. The purpose of this analysis is to predict the load-settlement behaviour and the capacity of traditional foundations. The comparison between field tests and numerical analysis and analytical solutions are also demonstrated.

The dimensional tolerances are among the most important manufacturing criteria in the forging of aerofoil blades for aero-engine applications. These are usually classified in relation to aerofoil shape, thickness and twist along cross-sections and are. affected by many factors such as preform shape, temperature, die elasticity and component springback. The interaction of all these factors makes design for net-shape forging extremely difficult. The focus was on the final shape of the forged aerofoil blades, which was affected by temperature, die elasticity and springback. A thermaIly coupled 3D thermo-elasto-visco-plastic analysis was employed to . investigate the entire forging cycle including forging, removal of dies and cooling of the aerofoil section, so that the e{fect of forging and post forging conditions such as temperature can be taken into account. A die-shape compensation strategy is presented in this· thesis. A di~ cleaning technology has been dev.eloped to represent the dies as a series of sampled points. These points were subsequently utilized as spline control points in the definition of a 3D spline surface. Methods are presented to evaluate the net shape error and to use this error data iri order to redefine the die . surface for a successive iteration. An iterative procedure was applied and a significant reduction in deviation from nominal shape was achieved The various components in this study have been developed in parallel to aIlow the successful integration of CAD, CAE, FE simulation and die shape modification in the optimisation of aerofoil blade dies to achieve net shape forging.

Statistical data suggests that increased use of evidence-based medical therapies has largely contributed to the decrease in American death rate caused by heart disease. And my studies are about two applications of magnetic resonance imaging (MRI) as a non-invasive approach in evidence-based health care research. In my first study, the achievement of a pulmonary valve replacement surgery was assessed on a patient with tetralogy of Fallot (TOF). In order to evaluate the remodeling of right ventricle, two biventricular finite element models were built up for pre-surgical images and post-surgical images. In my second study, 3D Lagrangian strains and torsion in the left ventricle of ten rats were investigated using Displacement ENcoding with Stimulated Echoes (DENSE) cardiac magnetic resonance (CMR) images. Tools written in MATLAB were developed for 2D contouring, 3D modeling, strain and torsion computations, and statistical comparison across subjects.

As rail transportation is significantly more virtuous than airplanes or cars in terms of greenhousegases emissions, its development is being encouraged in several European countries, includingSweden. In addition, the development of railway lines on which trains can travel at higher speeds ismade in Sweden with the integration of existing infrastructure. On railway bridges, an increased trainspeed potentially leads to an increase in vibrations during passage, for which the structure may not bedesigned. It is therefore essential to know the dynamic properties of the structures used.Several studies highlight the influence of friction phenomena in sliding bearings on the dynamicproperties of bridges equipped with them. This Master Thesis is based on previous works that led tothe development of a finite element modelling the friction mechanisms that occur in these bearings.The friction occurring between a PTFE sliding plate and a steel surface is thus modelled using the Bouc-Wen model, a model for hysteresis phenomena. The finite element was developed as a Fortransubroutine, which can be integrated into the finite element calculation software Abaqus as a "userdefinedelement". It allows friction to be modelled along the longitudinal direction of the bridge onlyand can therefore only be used in two-dimensional models. The user-defined element is also based ona model that takes into account the influence of contact pressure and sliding velocity on the steel-PTFEcoefficient of friction. As several studies indicate, contact temperature can also have a significantinfluence on the value of the coefficient of friction but is not taken into account in the current model.In this project, the previously developed finite element was therefore generalized to account forfriction in both directions of the sliding plate by the means of a two-dimensional generalization of theBouc-Wen model. Based on experimental data available in scientific literature, the model forcalculating the coefficient of friction was also extended to take into account the influence of thecontact temperature. In addition, a model to update the contact temperature based on the theory ofsurface heating of semi-infinite bodies has been incorporated. Finally, this thesis presents theintegration of this updated finite element on three-dimensional models of the Banafjäl Bridge, locatedin northern Sweden. Simulations to estimate the fundamental frequencies and resonance modes ofthe structure as well as the temperature increase that can occur in a bearing during the passage of atrain were carried out on this model.

The finite element method is a powerful analysis tool which has facilitated a better understanding of the behaviour of reinforced concrete structures. Its use in the research field is widespread and complements experimental tests and the development of new analytical models. Its application in practice engineering has permitted to deal with complex elements. However, the general structural engineer is still reluctant to consider finite element modelling for his work as he finds most of these models excessively sophisticated for his needs and knowledge. In particular, complexity of many finite element tools usually derives from the adoption of advanced concrete constitutive models. Implementation of more simple models based on engineering practice could facilitate its use by less experienced finite element users. In structural engineering practice finite element analysis can be of great usefulness to deal with those more problematic elements and/or where the application of traditional analysis methods presents limitations. This includes the so-called D-regions with a 3D behaviour. The strut-and-tie method and the stress field method are consistent and rational tools for the analysis and design of D-regions, but while their application to 2D elements is well covered in literature, its extension to 3D is problematic. This generally explains why excessively conservative assumptions are still common in the design of these elements. Refinement of current analytical and design approaches or the use of finite element analysis could lead to more rational solutions which in turn will reduce material requirements and costs. A 3D nonlinear finite element-based tool was developed in this thesis oriented towards the analysis and design of 3D D-regions by less experienced finite element users. Regarding material modelling, an orthotropic concrete model was adopted to permit the use of uniaxial stress-strain relationships. Only one single parameter, the uniaxial compressive strength of concrete, needs to be defined. Additionally, several aid functions were implemented, among which the following can be highlighted: a comprehensive, embedded reinforcement model to facilitate the introduction of complex rebar geometries; special support and load elements permitting an integrated and simple treatment of the boundary conditions imposed by them; and a simple design algorithm for the automatic determination of the required rebar areas. Three examples of applications to representative 3D D-regions are presented to show the capabilities of the tool. In particular, the analyses of fourteen four-pile caps, three socket base column-to-foundations connections and one anchorage block are described in the third part of the thesis. Results prove that realistic response predictions can be obtained considering relatively simple constitutive models. The capacity of the tool to configure consistent stress field models depending on the reinforcement arrangement is also demonstrated. The generation of rational reinforcement configurations by applying the implemented design algorithm is also shown. A strut-and-tie-based method for the analysis and design of four-pile caps with rectangular geometries is proposed in the fourth part. The method is based on a refined 3D strut-and-tie model and the consideration of three potential modes of failure: exceeding the reinforcement strength, crushing of the diagonal strut at the base of the column with narrowing of the strut and splitting of the diagonal strut due to transverse cracking. The main innovation is that the strut inclination is not fixed as in current strut-and-tie-based design procedures, but determined by maximizing the pile cap strength. The method accounts for strength softening of cracked concrete, compatibility constraints and reinforcement details. Its application to 162 specimens of literature led to very good predictions of the ultimate strength and, to a lesser extent, of the mode of failure.
El método de los elementos finitos es una potente herramienta de análisis que ha facilitado un mejor conocimiento del comportamiento de las estructuras de hormigón armado. Su uso en el ámbito de la investigación está ampliamente extendido. Su aplicación en la práctica ingenieril ha permitido la resolución de elementos complejos. Sin embargo, el ingeniero estructural común todavía es reticente a usar la modelización por elementos finitos ya que considera que la mayoría de estos modelos son excesivamente sofisticados para sus necesidades. La complejidad de muchas herramientas de elementos finitos suele derivarse de la adopción de modelos constitutivos de hormigón avanzados. La implementación de modelos más sencillos podría facilitar su uso por usuarios menos experimentados. En la práctica ingenieril el análisis con elementos finitos puede ser de gran utilidad para tratar aquellos elementos más problemáticos y/o donde la aplicación de los métodos de análisis tradicionales presenta limitaciones. Esto incluye las llamadas regiones D con comportamiento 3D. El método de bielas y tirantes y el método de campos de tensiones son herramientas racionales para el análisis y dimensionamiento de regiones D, pero su extensión a 3D es problemática. Este hecho explica por qué se adoptan todavía hipótesis excesivamente conservadoras en el dimensionamiento de estos elementos. La propuesta de métodos analíticos y de diseño más adecuados o la modelización con elementos finitos podría conducir a soluciones más racionales, lo que a su vez reduciría las necesidades de material y los costes. Como parte de esta tesis se ha desarrollado una herramienta de cálculo no lineal basada en el método de los elementos finitos orientada al análisis y dimensionamiento de regiones D tridimensionales por usuarios con menos experiencia en la modelización con elementos finitos. Se ha adoptado un modelo ortotrópico para el hormigón para permitir el uso de relaciones uniaxiales de tensión-deformación. Sólo es necesario definir un único parámetro, la resistencia a compresión uniaxial del hormigón. Adicionalmente, se han implementado varias funciones de ayuda, entre las que destacan: un modelo de armadura embebida para facilitar la introducción de geometrías de armado complejas; elementos especiales de apoyo y de carga que permiten un tratamiento integral de las condiciones de contorno; y un algoritmo de diseño para la determinación automática del área de armado necesaria. Se presentan tres ejemplos de aplicación a regiones D 3D representativas para mostrar las capacidades de la herramienta. En concreto, en la tercera parte del documento se describen los análisis de catorce encepados, tres cálices de cimentación y un bloque de anclaje. Los resultados muestran que se pueden obtener predicciones bastante realistas considerando modelos constitutivos relativamente sencillos. También se demuestra la capacidad de la herramienta para configurar modelos de campo de tensiones consistentes dependiendo de la configuración de armado. Además se muestra la capacidad del algoritmo de diseño para configurar disposiciones de armado racionales. En la cuarta parte se propone un método para el análisis y dimensionamiento de encepados sobre cuatro pilotes con geometría rectangular. El método se basa en un modelo 3D de bielas y tirantes refinado y la consideración de tres modos de fallo posibles: rotura del acero, aplastamiento de la biela diagonal en la base de la columna con estrechamiento de la misma y splitting de la biela diagonal debido a la fisuración transversal. La principal novedad es que el ángulo de la biela no se fija como en otros modelos, sino que se determina mediante la maximización de la resistencia del encepado. El método considera el debilitamiento de la resistencia del hormigón fisurado, condiciones de compatibilidad de deformaciones y detalles de armado. Su aplicación a 162 especímenes dio luga
El mètode dels elements finits és una potent eina d'anàlisi que ha facilitat un millor coneixement del comportament de les estructures de formigó armat. El seu ús en l'àmbit de la investigació està àmpliament estès. La seua aplicació en la pràctica enginyeril ha permès la resolució d'elements més complexos. No obstant això, l'enginyer estructural comú encara és reticent a fer servir la modelització per elements finits ja que considera que la majoria d'aquests models són excessivament sofisticats per a les seues necessitats i el seu conèixement. En concret, la complexitat de moltes eines d'elements finits sol derivar-se de l'adopció de models constitutius avançats de formigó. La implementació de models més senzills basats en la pràctica enginyeril podria facilitar el seu ús per a usuaris menys experimentats en la modelització amb elements finits. A la pràctica enginyeril l'anàlisi amb elements finits pot ser de gran utilitat per a tractar aquells elements més problemàtics i/o on l'aplicació dels mètodes d'anàlisi tradicionals presenta limitacions. Això inclou les anomenades regions D amb comportament 3D. El mètode de bieles i tirants i el mètode de camps de tensions són eines racionals per a l'anàlisi i dimensionament de regions D, però la seua extensió a 3D és problemàtica. Aquest fet explica per què s'adopten encara hipòtesis excessivament conservadores en el dimensionament d'aquests elements. La proposta de mètodes analítics i de disseny més adequats o la modelització amb elements finits podria conduir a solucions més racionals, amb el que també es reduirien les necessitats de material i els costos. Com a part d'aquesta tesi s'ha desenvolupat una eina de càlcul no lineal basada en el mètode dels elements finits orientada a l'anàlisi i dimensionament de regions D tridimensionals per a usuaris amb menys experiència en la modelització amb elements finits. S'ha adoptat un model ortotròpic per al formigó per permetre l'ús de relacions uniaxials de tensió-deformació. Només cal definir un únic paràmetre, la resistència a compressió uniaxial del formigó. Addicionalment, s'han implementat diverses funcions d'ajuda, entre les quals destaquen: un model d'armadura embeguda per facilitar la introducció de geometries d'armat complexes; elements especials de suport i de càrrega que permeten un tractament integral i senzill de les condicions de contorn; i un algoritme de disseny per a la determinació automàtica de l'àrea d'armat necessari. Es presenten tres exemples d'aplicació a regions D 3D representatives per mostrar les capacitats de l'eina. En particular, en la tercera part del document es descriuen les anàlisis de catorze enceps sobre quatre pilons, 3 calzes de fonamentació i un bloc d'ancoratge. Els resultats mostren que es poden obtenir prediccions prou realistes considerant models constitutius relativament senzills. També es demostra la capacitat de l'eina per configurar models de camp de tensions consistents depenent de la configuració d'armat. A més es mostra la capacitat de l'algoritme de disseny per configurar disposicions d'armat racionals. En la quarta part es proposa un mètode per a l'anàlisi i dimensionament d'enceps sobre quatre pilons amb geometria rectangular. El mètode es basa en un model 3D de bieles i tirants refinat i la consideració de tres modes de fallada possibles: trencament de l'acer, aixafament de la biela diagonal a la base de la columna amb estrenyiment de la mateixa i splitting de la biela diagonal per causa de la fissuració transversal. La principal novetat és que l'angle de la biela no es fixa com en els models actuals de bieles i tirants, sinó que es determina mitjan\c{c}ant la maximització de la resistència de l'encep. El mètode proposat considera el debilitament de la resistència del formigó fissurat, condicions de compatibilitat de deformacions i detalls d'armat. La seua aplicació a 162 espècimens de la liter
Meléndez Gimeno, C. (2017). A finite element-based approach for the analysis and design of 3D reinforced concrete elements and its applications to D-regions [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/86193
TESIS

In the course of thermal treatments, materials are usually subjected to continuous heating and cooling cycles during which microstructural evolution and mechanical interactions occur simultaneously at different length and time scales. Modeling of these processes necessitates dealing with inherent complexities such as large material property variations, complex couplings and boundary conditions, coupled heat and mass transfer mechanisms and phase transformations. In this study, a mathematical framework based on finite element method (FEM) capable of predicting temperature history, evolution of phases and internal stresses during heat treatment of metals and alloys was developed. The model was integrated into the commercial FEA software MSC.Marc®
by user subroutines. The accuracy of the model was verified by simulating the quenching of eccentrically drilled steel cylinders. Simulation results were justified via SEM observations and XRD residual stress measurements. According to the results, the model can effectively predict the trends in the distribution of microstructure and residual stresses with a remarkable accuracy.

The monopiles with typical diameters of 4-6m are mostly used as foundations for the offshore wind turbine structures to resist the vertical and horizontal loads. The response of the laterally loaded pile depends upon the soil and pile behavior and also on the interaction between them. Various methods for the development of p-y curves exist and the accuracy of such methods depends upon the data from which it was developed. The accuracy of the analysis of the pile depends on the accuracy at which the p-y curve represents the soil response to the lateral pile deflection because most of the commonly used p-y curve criteria are based on a very limited number of tests. The p-y curve recommended by American Petroleum Institute (API) code which is widely used by geotechnical engineers for the analysis of laterally loaded piles; however the method was developed for the slender piles with diameters up to approximately 2m. Hence several authors have claimed that p-y curves should be revised to be applicable for large diameter rigid pile (Augustesen, Sorensen, & Ibsen, 2010), (Achmus, Abdel-Rahman, & Kuo, 2008).The thesis presents the finite element model (FEM) analysis of monopiles in marine clay with the diameter from 1m to 6m by means of the Plaxis 3D Foundation. To avoid the influence of the pile flexibility and pile rotation, lateral translation of rigid piles is analyzed.The rigid pile and soil with hardening soil parameter is modeled in FEM Plaxis 3D Foundation. Piles are loaded with lateral loads to obtain the soil response on pile and pile deflection at different stages which in turn is used for the development of p-y curve and are compared with the API p-y curves. The result does not show any diameter effect in the initial stiffness but the initial stiffness of soil is very low in comparison to that given by the p-y curve (API code). Low stiffness compared to API p-y curve may be due to the selected parameters of the soil model.As the result shows difference in p-y curve from FEM 3D plaxis foundation and API code, some field test results for larger diameter will be fruitful to decide which method is valid for the analysis of the larger diameter monopiles.

In the context of interior noise reduction, the present work aims at proposing Finite Element (FE) solution strategies for interior structural-acoustic applications including 3D modelling of homogeneous and isotropic poroelastic materials, under timeharmonic excitations, and in the low frequency range. A model based on the Biot-Allard theory is used for the poroelastic materials, which is known to be very costly in terms of computational resources. Reduced models offer the possibility to enhance the resolution of such complex problems. However, their applicability to porous materials remained to be demonstrated.First, this thesis presents FE resolutions of poro-elasto-acoustic coupled problems using modal-based approaches both for the acoustic and porous domains. The original modal approach proposed for porous media, together with a dedicated mode selection and truncation procedure, are validated on 1D to 3D applications.In a second part, modal-reduced models are combined with a Padé approximants reconstruction scheme in order to further improve the efficiency.A concluding chapter presents a comparison and a combination of the proposed methods on a 3D academic application, showing promising performances. Conclusions are then drawn to provide indications for future research and tests to be conducted in order to further enhance the methodologies proposed in this thesis.
Dans le contexte de lutte contre les nuisances sonores, cette thèse porte sur le développement de méthodes de résolution efficaces par éléments finis, pour des problèmes de vibroacoustique interne avec interfaces dissipatives, dans le domaine des basses fréquences. L’étude se limite à l’utilisation de solutions passives telles que l’intégration de matériaux poreux homogènes et isotropes, modélisés par une approche fondée sur la théorie de Biot-Allard. Ces modèles étant coûteux en terme de résolution, un des objectifs de cette thèse est de proposer une approche modale pour la réduction du problème poroélastique, bien que l’adéquation d’une telle approche avec le comportement dynamique des matériaux poreux soit à démontrer.Dans un premier temps, la résolution de problèmes couplés élasto-poro-acoustiques par sous-structuration dynamique des domaines acoustiques et poreux est établie. L’approche modale originale proposée pour les milieux poroélastiques, ainsi qu’une procédure de sélection des modes significatifs, sont validées sur des exemples 1D à 3D.Une deuxième partie présente une méthode combinant l’utilisation des modèles réduits précédemment établis avec une procédure d’approximation de solution par approximants de Padé. Il est montré qu’une telle combinaison offre la possibilité d’accroître les performances de la résolution (allocation mémoire et ressources en temps de calcul).Un chapitre dédié aux applications permet d’évaluer et comparer les approches sur un problème académique 3D, mettant en valeur leurs performances encourageantes. Afin d’améliorer les méthodes établies dans cette thèse, des perspectives à ces travaux de recherche sont apportées en conclusion.

QC 20120224

FP6 Marie-Curie Smart Structures
FP7 Marie-Curie Mid-Frequency

In the context of interior noise reduction, the present work aims at proposing Finite Element (FE) solution strategies for interior structural-acoustic applications including 3D modelling of homogeneous and isotropic poroelastic materials, under timeharmonic excitations, and in the low frequency range. A model based on the Biot-Allard theory is used for the poroelastic materials, which is known to be very costly in terms of computational resources. Reduced models offer the possibility to enhance the resolution of such complex problems. However, their applicability to porous materials remained to be demonstrated.First, this thesis presents FE resolutions of poro-elasto-acoustic coupled problems using modal-based approaches both for the acoustic and porous domains. The original modal approach proposed for porous media, together with a dedicated mode selection and truncation procedure, are validated on 1D to 3D applications.In a second part, modal-reduced models are combined with a Padé approximants reconstruction scheme in order to further improve the efficiency.A concluding chapter presents a comparison and a combination of the proposed methods on a 3D academic application, showing promising performances. Conclusions are then drawn to provide indications for future research and tests to be conducted in order to further enhance the methodologies proposed in this thesis.