SIMULATING LARGE DEFORMATION OF INTRANEURAL GANGLION CYST USING FINITE ELEMENT METHOD

Intraneural Ganglion Cyst is disorder observed in the nerve injury, it is still unknown and very difficult to predict its propagation in the human body so many times it is referred as an unsolved history. The treatments for this disorder are to remove the cystic substance from the nerve by a surgery. However these treatments may result in neuropathic pain and recurrence of the cyst. The articular theory proposed by Spinner et al., (Spinner et al. 2003) considers the neurological deficit in Common Peroneal Nerve (CPN) branch of the sciatic nerve and adds that in addition to the treatment, ligation of articular branch results into foolproof eradication of the deficit. Mechanical modeling of the affected nerve cross section will reinforce the articular theory (Spinner et al. 2003). As the cyst propagates, it compresses the neighboring fascicles and the nerve cross section appears like a signet ring. Hence, in order to mechanically model the affected nerve cross section; computational methods capable of modeling excessively large deformations are required. Traditional FEM produces distorted elements while modeling such deformations, resulting into inaccuracies and premature termination of the analysis. The methods described in research report have the capability to simulate large deformation. The results obtained from this research shows significant deformation as compared to the deformation observed in the conventional finite element models. The report elaborates the neurological deficit followed by detail explanation of the Smoothed Particle Hydrodynamic approach. Finally, the results show the large deformation in stages and also the successful implementation of the SPH method for the large deformation of the biological organ like the Intra-neural ganglion cyst.

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

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

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

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

Finite element modeling of tube deformation during cold pilgering

A three-dimensional finite element model of cold pilgering of stainless steel tubes is developed in this paper. The objective is to use the model to increase the understanding of forces and deformations in the process. The focus is on the influence of vertical displacements of the roll stand and axial displacements of the mandrel and tube. Therefore, the rigid tools and the tube are supported with elastic springs. Additionally, the influences of friction coefficients in the tube/mandrel and tube/roll interfaces are examined. A sensitivity study is performed to investigate the influences of these parameters on the strain path and the roll separation force. The results show the importance of accounting for the displacements of the tube and rigid tools on the roll separation force and the accumulative plastic strain.

Introducing the Logarithmic finite element method: a geometrically exact planar Bernoulli beam element

We propose a novel finite element formulation that significantly reduces the number of degrees of freedom necessary to obtain reasonably accurate approximations of the low-frequency component of the deformation in boundary-value problems. In contrast to the standard Ritz–Galerkin approach, the shape functions are defined on a Lie algebra—the logarithmic space—of the deformation function. We construct a deformation function based on an interpolation of transformations at the nodes of the finite element. In the case of the geometrically exact planar Bernoulli beam element presented in this work, these transformation functions at the nodes are given as rotations. However, due to an intrinsic coupling between rotational and translational components of the deformation function, the formulation provides for a good approximation of the deflection of the beam, as well as of the resultant forces and moments. As both the translational and the rotational components of the deformation function are defined on the logarithmic space, we propose to refer to the novel approach as the “Logarithmic finite element method”, or “LogFE” method.

A finite element framework for the interplay between delamination and buckling of rubber-like bi-material systems and stretchable electronics

In this study, a finite element (FE) framework for the analysis of the interplay between buckling and delamination of thin layers bonded to soft substrates is proposed. The current framework incorporates the following modeling features: (i) geometrically nonlinear solid shell elements, (ii) geometrically nonlinear cohesive interface elements, and (iii) hyperelastic material constitutive response for the bodies that compose the system. A fully implicit Newton–Raphson solution strategy is adopted to deal with the complex simultaneous presence of geometrical and material nonlinearities through the derivation of the consistent FE formulation. Applications to a rubber-like bi-material system under finite bending and to patterned stiff islands resting on soft substrate for stretchable solar cells subjected to tensile loading are proposed. The results obtained are in good agreement with benchmark results available in the literature, confirming the accuracy and the capabilities of the proposed numerical method for the analysis of complex three-dimensional fracture mechanics problems under finite deformations.

A Trans-dimensional inversion algorithm to model deformation sources with unconstrained shape in finite element domains

Ground deformation provides valuable insights on subsurface processes with pattens reflecting the characteristics of the source at depth. In active volcanic sites displacements can be observed in unrest phases; therefore, a correct interpretation is essential to assess the hazard potential. Inverse modeling is employed to obtain quantitative estimates of parameters describing the source. However, despite the robustness of the available approaches, a realistic imaging of these reservoirs is still challenging. While analytical models return quick but simplistic results, assuming an isotropic and elastic crust, more sophisticated numerical models, accounting for the effects of topographic loads, crust inelasticity and structural discontinuities, require much higher computational effort and information about the crust rheology may be challenging to infer. All these approaches are based on a-priori source shape constraints, influencing the solution reliability. In this thesis, we present a new approach aimed at overcoming the aforementioned limitations, modeling sources free of a-priori shape constraints with the advantages of FEM simulations, but with a cost-efficient procedure. The source is represented as an assembly of elementary units, consisting in cubic elements of a regular FE mesh loaded with a unitary stress tensors. The surface response due to each of the six stress tensor components is computed and linearly combined to obtain the total displacement field. In this way, the source can assume potentially any shape. Our tests prove the equivalence of the deformation fields due to our assembly and that of corresponding cavities with uniform boundary pressure. Our ability to simulate pressurized cavities in a continuum domain permits to pre-compute surface responses, avoiding remeshing. A Bayesian trans-dimensional inversion algorithm implementing this strategy is developed. 3D Voronoi cells are used to sample the model domain, selecting the elementary units contributing to the source solution and those remaining inactive as part of the crust.

Assessing model credibility of CT-Based finite element model to predict bone strength

Osteoporosis is one of the major causes of mortality among the elderly. Nowadays, areal bone mineral density (aBMD) is used as diagnostic criteria for osteoporosis; however, this is a moderate predictor of the femur fracture risk and does not capture the effect of some anatomical and physiological properties on the bone strength estimation. Data from past research suggest that most fragility femur fractures occur in patients with aBMD values outside the pathological range. Subject-specific finite element models derived from computed tomography data are considered better tools to non-invasively assess hip fracture risk. In particular, the Bologna Biomechanical Computed Tomography (BBCT) is an In Silico methodology that uses a subject specific FE model to predict bone strength. Different studies demonstrated that the modeling pipeline can increase predictive accuracy of osteoporosis detection and assess the efficacy of new antiresorptive drugs. However, one critical aspect that must be properly addressed before using the technology in the clinical practice, is the assessment of the model credibility. The aim of this study was to define and perform verification and uncertainty quantification analyses on the BBCT methodology following the risk-based credibility assessment framework recently proposed in the VV-40 standard. The analyses focused on the main verification tests used in computational solid mechanics: force and moment equilibrium check, mesh convergence analyses, mesh quality metrics study, evaluation of the uncertainties associated to the definition of the boundary conditions and material properties mapping. Results of these analyses showed that the FE model is correctly implemented and solved. The operation that mostly affect the model results is the material properties mapping step. This work represents an important step that, together with the ongoing clinical validation activities, will contribute to demonstrate the credibility of the BBCT methodology.

Root cause analysis applied to a finite element model's refinement of a negative stiffness structure

Negative Stiffness Structures are mechanical systems that require a decrease in the applied force to generate an increase in displacement. They are structures that possess special characteristics such as snap-through and bi-stability. All of these features make them particularly suitable for different applications, such as shock-absorption, vibration isolation and damping. From this point of view, they have risen awareness of their characteristics and, in order to match them to the application needed, a numerical simulation is of great interest. In this regard, this thesis is a continuation of previous studies in a circular negative stiffness structure and aims at refine the numerical model by presenting a new solution. To that end, an investigation procedure is needed. Amongst all of the methods available, root cause analysis was the chosen one to perform the investigation since it provides a clear view of the problem under analysis and a categorization of all the causes behind it. As a result of the cause-effect analysis, the main causes that have influence on the numerical results were obtained. Once all of the causes were listed, solutions to them were proposed and it led to a new numerical model. The numerical model proposed was of nonlinear type of analysis with hexagonal elements and a hyperelastic material model. The results were analyzed through force-displacement curves, allowing for the visualization of the structure’s energy recovery. When compared to the results obtained from the experimental part, it is evident that the trend is similar and the negative stiffness behaviour is present.

Comparison between one-dimensional uncoupled and convection-conduction conjugated approaches in finned surface heat transfer

This work studies the forced convection problem in internal flow between concentric annular ducts, with radial fins at the internal tube surface. The finned surface heat transfer is analyzed by two different approaches. In the first one, it is assumed one-dimensional heat conduction along the internal tube wall and fins, with the convection heat transfer coefficient being a known parameter, determined by an uncoupled solution. In the other way, named conjugated approach, the mathematical model (continuity, momentum, energy and K-epsilon equations) applied to tube annuli problem was numerically solved using finite element technique in a coupled formulation. At first time, a comparison was made between results obtained for the conjugated problem and experimental data, showing good agreement. Then, the temperature profiles under these two approaches were compared to each other to analyze the validity of the one-dimensional classical formulation that has been utilized in the heat exchanger design.

Análise de treliças planas viscoelásticas com deformações finitas e pórticos planos viscoelásticos com vigas mistas

Este trabalho compõe-se de duas partes. A primeira parte propõe-se a apresentar um estudo e um programa computacional para a análise não linear geométrica de treliças planas com propriedades: viscoelásticas. Na segunda parte, tem-se o estudo e um programa sobre pórticos planos com propriedades viscoelásticas, usando o modelo reológico standard e o dado pelo CEB. Leva-se em consideração o efeito de temperatura e retração nesta análise. Estende-se o trabalho sobre pórtico para o estudo sobre vigas mistas, levando em consideração a mudança da linha neutra. A formulação está baseada no método dos elementos finitos para grandes deformações, particularizada para treliça e pórtico. É feita a descrição de ambos os programas e rodados diversos exemplos.

Avaliação de métodos eletromagnéticos aplicando campos polarizados e focalizados

Em investigações geofísicas rasas que empregam os métodos eletromagnéticos indutivos mais avançados, alvos com baixo número de indução (Low Induction Number – LIN) produzem anomalias eletromagnéticas muito baixas e de difícil interpretação. Para suprir esta deficiência, neste trabalho são estudados a aplicabilidade de campos eletromagnéticos polarizados e focalizados – POLFOCEM como fonte primária de indução. Os campos E.M. focalizados e polarizados, vertical e horizontalmente, são obtidos pelas combinações vetoriais de pares de dipolos transmissores e, ocorrem na região central entre eles. A focalização é observada nesta região na profundidade de 0,25 do espaçamento entre esses transmissores – L. Portanto, máximos acoplamentos podem ser obtidos através da seleção da polarização de acordo com a geometria do alvo, ocorrendo um aumento na densidade de fluxo magnético sobre ele e, máximas anomalias produzidas. É utilizada uma metodologia numérica para o cômputo dessas anomalias por meio da técnica dos elementos finitos para solução do problema 2,5-D. Em todos os experimentos numéricos são realizadas comparações qualitativas e quantitativas entre as respostas obtidas pelos sistemas POLFOCEM e convencional, o qual emprega um único dipolo como transmissor (dipolo-dipolo). As anomalias produzidas pelo sistema POLFOCEM, em que os dipolos transmissores são acionados simultaneamente, correspondem à soma das anomalias produzidas por cada um desses dipolos independentes, caracterizando, desta forma, a linearidade dos campos eletromagnéticos. Os experimentos numéricos são realizados para alvos prismáticos bidimensionais com três diferentes inclinações, inseridos num semi-espaço resistivo, e para as freqüências das fontes na faixa das ondas de rádio. As anomalias assimétricas no sistema convencional, que se tornam simétricas no sistema POLFOCEM, apresentam valores menores em amplitude. Contudo, aquelas anomalias tanto assimétricas quanto simétricas que se tornam anti-simétricas apresentam valores maiores. Em decorrência dessas diminuições e aumentos nas amplitudes ocorrem rotações nos diagramas de Argand, no sentido horário e anti-horário para alvos com baixos valores de condutividade, respectivamente. Em experimentos de identificação de presença de dois alvos próximos, o sistema convencional é capaz de identificá-los primeiramente, prevalecendo o seu uso.

Modelagem direta de dados de eletrorresistividade 3-D

Dentre os métodos geofísicos, o da Eletrorresistividade é um dos mais tradicionais, com o seu desenvolvimento ocorrido há mais de 80 anos. Durante esse tempo o seu uso acompanhou o avanço no poder de processamento numérico e mais recentemente, a modelagem e inversão tridimensional tornou-se uma possibilidade para o geofísico. Apresentamos, neste trabalho, a técnica de elementos finitos aplicada ao método da eletrorresistividade 3-D, através do cálculo do potencial secundário. Para o desenvolvimento da metodologia, simulamos o levantamento do método da eletrorresistividade 3-D com os arranjos Dipolo-Dipolo e Schlumberger, visando medir as variações laterais e verticais da resistividade aparente do solo. Estes arranjos consistem na injeção de corrente elétrica na superfície e de medidas de diferenças de potencial elétrico, resultante da interação da corrente elétrica com o solo. Sendo que, as fontes e receptores são localizados de acordo com os arranjos escolhidos para o levantamento. Neste trabalho, as curvas de sondagem e as pseudo-secções de resistividade aparente, são obtidas através da modelagem de eletrorresistividade 3-D, usando malha de elementos finitos regular. Para efeito de validação, os resultados são comparados com a resposta 3-D obtida a partir dos potenciais totais.

An algorithm for graphical computer results in shells

A Finite Element technique to interpolate general data (function values and its derivatives) has been developped. The technique can be considered as a generalized solution of the classical polynomial interpolation, because the condition for the interpolating function to be a polynomial is replaced by a minimizing condition of a given “smoothing” functional. In this way it is possible to find interpolating functions with a given level of continuity according to the class of finite elements used. Examples have been presented in order to assess the accuracy and efficiency of the procedure.

Yield functions for porous materials with cubic symmetry using different definitions of yield

Plastic yield criteria for porous ductile materials are explored numerically using the finite-element technique. The cases of spherical voids arranged in simple cubic, body-centred cubic and face-centred cubic arrays are investigated with void volume fractions ranging from 2 % through to the percolation limit (over 90 %). Arbitrary triaxial macroscopic stress states and two definitions of yield are explored. The numerical data demonstrates that the yield criteria depend linearly on the determinant of the macroscopic stress tensor for the case of simple-cubic and body-centred cubic arrays - in contrast to the famous Gurson-Tvergaard-Needleman (GTN) formula - while there is no such dependence for face-centred cubic arrays within the accuracy of the finite-element discretisation. The data are well fit by a simple extension of the GTN formula which is valid for all void volume fractions, with yield-function convexity constraining the form of the extension in terms of parameters in the original formula. Simple cubic structures are more resistant to shear, while body-centred and face-centred structures are more resistant to hydrostatic pressure. The two yield surfaces corresponding to the two definitions of yield are not related by a simple scaling.