<div>Underground structures must be able to support static overburden loads, as well as to accommodate additional deformations imposed by seismic motions. Progress has been made in the last few years in understanding the soil-structure interaction mechanisms and the stress and displacement transfer from the ground to the structure during a seismic event. It seems well established that, for most tunnels, the most critical demand to the structure is caused by shear waves traveling perpendicular to the tunnel axis. Those waves cause distortions of the cross section (ovaling for a circular tunnel, and racking for a rectangular tunnel) that result in axial forces (thrusts) and bending moments. While all this has been well-studied for structures placed in linear-elastic ground, there is little information regarding the behavior of buried structures placed in nonlinear ground, especially under undrained conditions, i.e., when excess pore pressures generate and accumulate during the earthquake.</div><div><br></div><div><div>Two-dimensional dynamic numerical analyses are conducted to assess the seismic response of deep circular tunnels located far from the seismic source, under drained or undrained loading conditions. It is assumed that the liner remains elastic and that plane strain conditions apply. </div><div> A new cyclic nonlinear elastoplastic constitutive model is developed and verified, to simulate the nonlinear behavior and excess pore pressures accumulation with cycles of loading in the ground. The results of the numerical analyses show negligible effect of input frequencies on the normalized distortions of a tunnel for input frequencies smaller than 5 Hz (the distortions of the tunnel are normalized with respect to those of the free field); that is, for ratios between the wavelength of the seismic input and the tunnel opening larger than about eight to ten. The results also show that undrained conditions, compared with drained conditions, tend to reduce deformations for flexible liners and increase them for stiffer tunnels, when no accumulation of pore pressures with cycles of loading is assumed. However, when pore pressures increase with the number of cycles, the differences in distortions between drained and undrained loading are reduced, i.e., the normalized distortions increase for flexible and decrease for stiff tunnels, compared to those with drained conditions. </div></div><div><br></div><div><div>Undrained loading produces larger thrust in the liner than drained loading for stiff tunnels with flexibility ratio F ≤ 2.0.</div><div>For more flexible tunnels with F > 2.0, the behavior is the opposite, i.e., smaller axial forces are obtained for undrained loading than for drained loading. Including excess pore pressure accumulation does not introduce significant changes in the axial forces of the liner, irrespective of the flexibility of the tunnel, compared to those obtained from undrained loading without pore pressure accumulation.</div><div>The drainage loading condition (drained or undrained) or the magnitude of the free-field excess pore pressures during undrained loading do not affect the normalized bending moments for flexible tunnels, with F ≥ 2. For stiffer tunnels, with F < 2, the normalized bending moments increase from drained to undrained loading, and with the free field excess pore pressures.</div></div><div><br></div><div><div>It is found that the tunnel’s response is determined by the load on the liner, or by the distortions of the cross section, depending on the flexibility ratio. For stiff structures, with F ≤ 2.0, important axial forces and bending moments are produced in the structure, with larger magnitudes for the undrained case; while the distortions of the cross section are very small. When the tunnel becomes more flexible, the loading on the liner decreases, but the distortions of the cross section start to be important. For flexible structures with initial F ≥ 10 (for the cases investigated), the performance is largely determined by the distortions of the cross section, while the axial forces and bending moments are almost negligible. Such distortions are drastically affected by the drainage loading condition and by the magnitude of pore pressures in the free field. </div></div><div><br></div>

O trabalho trata da formulação e análise da resposta numérica de um modelo de danificação e plastificação localizadas em zonas previamente definidas ao longo de elementos estruturais de barra, estendendo-se os conceitos dos modelos clássicos de plasticidade concentrada. Inicialmente é feita uma breve revisão bibliográfica sobre modelos de dano e plasticidade, destacando-se alguns que contém propostas voltadas para a simplificação de suas formulações ou então da etapa de implementação numérica adotando-se, neste caso, a técnica de discretização em estratos dos elementos de barra. Em seguida, já no âmbito dos modelos ditos simplificados e estendendo-se a formulação de um modelo para estruturas de barras proposto na bibliografia, deduz-se a matriz de rigidez para o caso de barra em regime elástico com dano. Na dedução, considera-se que os processos de danificação, associados à flexão e restritos às extremidades do elemento de viga, são dependentes entre si e afetam diretamente os coeficientes de transmissão de esforços ao longo do elemento. Nesta etapa, através de dois exemplos numéricos, avalia-se a resposta do modelo quando implementado em um programa de elementos finitos. Apresenta-se, na seqüência, um estudo sobre a viabilidade da aplicação do modelo estudado às estruturas em concreto armado, incluindo-se casos em que a distribuição da armadura longitudinal é assimétrica. Como resultado desse estudo, sugere-se uma modificação da função critério de danificação, originalmente empregada, para melhor representar o processo de dissipação de energia relacionada à danificação. A partir dos confrontos entre respostas experimentais de vigas e pórtico em concreto armado e resultados numéricos, conclui-se pelo bom desempenho do modelo modificado. Finalmente, reúnem-se sugestões para a continuidade dos estudos. / This work is related to formulation and numerical analysis of a damage and plasticity model which considers such effects localised on zones previously defined along the structural beam elements, extending the concepts of the classics lumped plasticity models. First of all, a brief bibliography revision on damage and plasticity based models is given, enlightening the ones which propose some kind of simplification at the level of the formulation or only at the stage of the numerical implementation, for instance using a layered technique in a beam element discretization. Afterwards, on the field of the so called simplified models, a stiffness matrix of an elastic damage beam element is deduced, extending the formulation of a framed structures model proposed in the bibliography. The main characteristic of the new element is that the damage processes is localized on hinges at the extremities of the beam element, being the associated damage variables dependent on each other and reducing progressively the transmission efforts factor along the element. Two basic numerical examples show the performance of the model when implemented in a finite element code. Next, a study on the feasibility of the model to analyse reinforced concrete structures is presented. In order to enlarge the application field, cases where longitudinal reinforcement is asymmetrically distributed in the cross section are also considered. As a result of the study, a change in the damage criteria function originally used is suggested, aiming to improve the valuation of the dissipated energy related to the damage process. The response of the modified resulting model is valuated by a confront between experimental and numerical results of beams and frame reinforced concrete structures. The results show a very good performance of the modified model. Finally, some topics for further research on the theme explored in this work are suggested.

Tunnelling in urban areas causes short-term and long-term ground movements under existing buildings. Finite element analysis provides a useful option for assessing the likely extent of damage induced in these buildings. Although finite element analysis is suggested to be used in the later phases of the building damage assessment procedures employed in practice, only the effect of short-term ground movements is typically considered and there are no detailed guidelines on the specification and complexity of the modelling. This thesis addresses the tunnel-soil-building interaction problem and the effect of long-term consolidation, as well as demonstrating the application of 3D finite element analysis with appropriate simplifications for practical assessment purposes. Finite element models are developed to quantify the effect of shallow tunnelling on an example masonry building founded on strip footings, considering both single- and twin-tunnel scenarios in a typical London soil profile. Total stress and effective stress analyses are adopted with specific modelling procedures to focus on the short-term and long-term response respectively. The analyses use a non-linear model for the masonry, and allow slippage and gapping at the soil-footing interface. Two advanced constitutive models for the soil (the extended Mohr-Coulomb model and the modified two-surface kinematic hardening model) are implemented with customized stress update schemes. The finite element results present the interaction between the soil and the building by comparing with the greenfield ground response. The horizontal coupling between the foundation and the ground is shown to be relatively weak. The dominant deformation mode of the building varies with the tunnel configuration (i.e. single or twin tunnels) and the tunnel eccentricity. Strain localization is found around the explicitly modelled window and door openings. The long-term consolidation is sensitive to the permeability of the tunnel lining. The building response to long-term ground movements is further affected by the tunnel-tunnel interaction in the case of twin-tunnel configuration. Performing 3D analysis of a single facade and foundation provides useful damage predictions, without the need to model a complete building. The proposed result processing methods such as characteristic strain and damage bar chart are practical tools for assessment. The study highlights some limitations of the elastic beam assessment method, which is often adopted in the early phase of the damage assessment process.

ABSTRACT 3-D SOIL STRUCTURE INTERACTION ANALYSES OF THREE IDENTICAL BUILDINGS IN SAKARYA CITY AFTER 17 AUGUST 1999 KOCAELI EARTHQUAKE &Uuml / nal,Orhan M.S., Department of Civil Engineering, Supervisor: Assist. Prof. Dr Kemal &Ouml / nder &Ccedil / etin October 2003, 116 Pages The aim of this study is to analyze the soil structure interaction of three identical buildings on &ordf / ahinler Street of Sakarya city which had no damage to heavy damage after the Kocaeli (1999) earthquake. For the purpose of 3-D dynamic nonlinear analysis of the soil site and the overlying structures, Flac3D software was chosen as the numerical modeling framework. Soil properties were determined by using the results of available site investigation studies. A three dimensional mesh was created to represent the topographic and geometric constraints of the problem. Linearly elastic perfectly plastic constitutive model was implemented to model the soil behavior. The results of 3-D dynamic numerical analyses in the forms of acceleration, displacement, strain, stress and pore pressure were presented. The higher acceleration, strain and stress levels calculated under the collapsed building can be attributed as the major cause of poor performance of the structure.

This thesis presents a whole-embryo finite element model of neurulation -- the first of its kind. An advanced, multiscale finite element approach is used to capture the mechanical interactions that occur across cellular, tissue and whole-embryo scales. Cell-based simulations are used to construct a system of constitutive equations for embryonic tissue fabric evolution under different scenarios including bulk deformation, cell annealing, mitosis, and Lamellipodia effect. Experimental data are used to determine the parameters in these equations. Techniques for obtaining images of live embryos, serial sections of fixed embryo fabric parameters, and material properties of embryonic tissues are used. Also a spatial-temporal correlation system is introduced to organize and correlate the data and to construct the finite element model. Biological experiments have been conducted to verify the validity of this constitutive model. A full functional finite element analysis package has been written and is used to conduct computational simulations. A simplified contact algorithm is introduced to address the element permeability issue. Computational simulations of different cases have been conducted to investigate possible causes of neural tube defects. Defect cases including neural plate defect, non-neural epidermis defect, apical constriction defect, and convergent extension defect are compared with the case of normal embryonic development. Corresponding biological experiments are included to support these defect cases. A case with biomechanical feedbacks on non-neural epidermis is also discussed in detail with biological experiments and computational simulations. Its comparison with the normal case indicates that the introduction of biomechanical feedbacks can yield more realistic simulation results.

This thesis presents a whole-embryo finite element model of neurulation -- the first of its kind. An advanced, multiscale finite element approach is used to capture the mechanical interactions that occur across cellular, tissue and whole-embryo scales. Cell-based simulations are used to construct a system of constitutive equations for embryonic tissue fabric evolution under different scenarios including bulk deformation, cell annealing, mitosis, and Lamellipodia effect. Experimental data are used to determine the parameters in these equations. Techniques for obtaining images of live embryos, serial sections of fixed embryo fabric parameters, and material properties of embryonic tissues are used. Also a spatial-temporal correlation system is introduced to organize and correlate the data and to construct the finite element model. Biological experiments have been conducted to verify the validity of this constitutive model. A full functional finite element analysis package has been written and is used to conduct computational simulations. A simplified contact algorithm is introduced to address the element permeability issue. Computational simulations of different cases have been conducted to investigate possible causes of neural tube defects. Defect cases including neural plate defect, non-neural epidermis defect, apical constriction defect, and convergent extension defect are compared with the case of normal embryonic development. Corresponding biological experiments are included to support these defect cases. A case with biomechanical feedbacks on non-neural epidermis is also discussed in detail with biological experiments and computational simulations. Its comparison with the normal case indicates that the introduction of biomechanical feedbacks can yield more realistic simulation results.

This study presents a new approach for nonlinear multi-scale constitutive models using artificial neural networks (ANNs). Three ANN classes are proposed to characterize the nonlinear multi-axial stress-strain behavior of metallic, polymeric, and fiber reinforced polymeric (FRP) materials, respectively. Load-displacement responses from nanoindentation of metallic and polymeric materials are used to train new generation of dimensionless ANN models with different micro-structural properties as additional variables to the load-deflection. The proposed ANN models are effective in inverse-problems set to back-calculate in-situ material parameters from given overall nanoindentation test data with/without time-dependent material behavior. Towards that goal, nanoindentation tests have been performed for silicon (Si) substrate with/without a copper (Cu) film. Nanoindentation creep test data, available in the literature for Polycarbonate substrate, are used in these inverse problems. The predicted properties from the ANN models can also be used to calibrate classical constitutive parameters. The third class of ANN models is used to generate the effective multi-axial stress-strain behavior of FRP composites under plane-stress conditions. The training data are obtained from coupon tests performed in this study using off-axis tension/compression and pure shear tests for pultruded FRP E-glass/polyester composite systems. It is shown that the trained nonlinear ANN model can be directly coupled with finite-element (FE) formulation as a material model at the Gaussian integration points of each layered-shell element. This FE-ANN modeling approach is applied to simulate an FRP plate with an open-hole and compared with experimental results. Micromechanical nonlinear ANN models with damage formulation are also formulated and trained using simulated FE modeling of the periodic microstructure. These new multi-scale ANN constitutive models are effective and can be extended by including more material variables to capture complex material behavior, such as softening due to micro-structural damage or failure.

[Truncated abstract] An important field of research in orthopaedic biomechanics is the elucidation and mathematical modelling of the mechanical response of cartilaginous tissues. Such research has applications in the understanding of joint function and degenerative processes, as well as in surgical planning and simulation, and engineering of tissue replacements. In the case of surgical and tissue engineering applications especially, patient-specific mechanical properties are highly desirable. Unfortunately, obtaining such information would generally involve destructive mechanical testing of patient tissue, thus rendering the tissue functionally unusable. Development of a laser scanning confocal arthroscope (LSCA) within our School will soon allow non-invasive extraction of 3D microstructural images of cartilaginous tissues in vivo. It is also envisaged that, linked to a suitably formulated constitutive formulation, such information could allow estimation of tissue mechanical response without physical biopsy. This thesis describes the development of techniques to potentially allow non-invasive patient-specific estimation of tissue mechanical response based on confocal arthroscopy data. A microstructural constitutive model is developed which is capable of directly incorporating LSCA-derived patient-specific structural information. A fibre composite type homogenisation approach is used as the basis for the model. ... The result is a series of orientation tensors describing the 3D orientation of linear features in the image stack. The developed analysis techniques are used to estimate fibre volume fraction and orientation distribution for each of the meniscal specimens. The developed constitutive model and image-derived structural parameters are finally used to estimate the reaction force history of two meniscal cartilage specimens subjected to partially confined compression. The predictions are made on the basis of the specimens? individual structural condition as assessed by confocal microscopy and involve no tuning of material parameters. Although the model does not reproduce all features of the experimental curves, as an unfitted estimate of mechanical response the prediction is quite accurate. In light of the obtained results it is judged that more general non-invasive estimation of tissue mechanical properties is possible using the developed framework. The likely limitations and potential areas of improvement are discussed.

This work presents a strategy for the coupled poro-elasto-plastic formulation. The Finite Element Method (FEM) is used to solve the differential equations, interpolating displacement and pore pressure fields. This problem is solved fully coupled, based on an only one system of equations. The nonlinear problem is globally solved by the Newton-Raphson procedure, and the Closest Point algorithm is implemented for the returning map in the elasto-plastic models. Based on a computational module that has already been developed (PORO), which is written using C++ language and Object-Oriented Programming (OOP), this work expands this program creating new classes for different elasto-plastic constitutive models. The program is verified by classical examples in the literature such as the poro-elastic column and the problem of Schiffman. Some strategies for optimization the computational cost are presented, which use specialized math libraries (MKL) and code parallelization (OpenMP). / Este trabalho apresenta, discute e implementa a formulação poro-elastoplástica fortemente acoplada. A discretização espacial das equações diferenciais governantes é realizada através do Método dos Elementos Finitos (MEF), com interpolação do campo de deslocamento e da poropressão. O problema poro-mecânico é resolvido de forma totalmente acoplada, com base em um único sistema de equações. O método iterativo de Newton-Rhapson é empregado para a solução global do problema não linear, tendo ainda o algoritmo implícito iterativo Closest Point para a integração local das equações da plasticidade. Baseando-se em um programa computacional pré-existente denominado PORO, escrito na linguagem C++ e que utiliza o paradigma de Programação Orientada a Objetos (POO), faz-se a adaptação desse código através da criação de novas classes para permitir o uso de modelos constitutivos elastoplásticos e lei de fluxo associada no acoplamento poro-mecânico. Para verificação do programa são analisados problemas clássicos da literatura, a exemplo da coluna poro-elástica e o caso de Schiffman. Descrevem-se ainda algumas estratégias de otimização do custo computacional, implementando-se o uso de bibliotecas matemáticas (MKL) e paralelização do código (OpenMP).

O conhecimento do comportamento mecânico do material é vital para a previsão da durabilidade e atendimento aos requisitos iniciais de projeto. Por meio de ensaios experimentais é possível mensurar esse comportamento e utilizar como parâmetro de entrada em simulações numéricas. Nesse ponto a temperatura altera significativamente as propriedades mecânicas estáticas e dinâmicas do material. Testes de fadiga devem ser considerados na caracterização do material mesmo que um alto custo e tempo sejam necessários. Assim o objetivo deste trabalho é avaliar o comportamento termo-mecânico do policloropreno em tração e fadiga. Este estudo se fez necessário para prever a durabilidade em diversas temperaturas em mancais de rolamento de cardans que exigem rigidez e durabilidade. A abordagem proposta consiste em formular equações constitutivas que considerem as constantes dependentes da temperatura para a situação estática (ensaio de tração) e dinâmica (ensaio de fadiga). Foram estudados 4 modelos constitutivos já consagrados na literatura. Corpos de prova normatizados foram utilizados para a avaliação experimental do material tanto em tração quanto em fadiga. Pode-se concluir que esta abordagem conseguiu modelar adequadamente o comportamento mecânico em ambas as situações estudadas. / Mechanical behavior of materials is vital for predicting durability and compliance with initial requirements of the project. Through experimental tests can measure this behavior and use as an input parameter in numerical simulations. At this point temperature significantly changes static and dynamic mechanical properties of the material. Fatigue tests should be considered in material characterization even if a high-cost and time are required. So the aim of this study is evaluate the thermo-mechanical behavior of polychloroprene in tension and fatigue. This study was necessary to predict durability in various temperatures in driveshafts of rolling bearings that requires rigidity and durability. The proposed approach is formulating constitutive equations to consider the dependent constant temperature for static situation (tensile test) and dynamic (fatigue test). Four constitutive models already established in literature were studied. Standardized test specimens were used for the experimental evaluation of the material both in tension and in fatigue. It can be concluded that this approach could adequately model the mechanical behavior in both cases studied.