A mathematical framework for finite strain elastoplastic consolidation. Part 1: Balance laws, variational formulation, and linearization

A mathematical formulation for finite strain elasto plastic consolidation of fully saturated soil media is presented. Strong and weak forms of the boundary-value problem are derived using both the material and spatial descriptions. The algorithmic treatment of finite strain elastoplasticity for the solid phase is based on multiplicative decomposition and is coupled with the algorithm for fluid flow via the Kirchhoff pore water pressure. Balance laws are written for the soil-water mixture following the motion of the soil matrix alone. It is shown that the motion of the fluid phase only affects the Jacobian of the solid phase motion, and therefore can be characterized completely by the motion of the soil matrix. Furthermore, it is shown from energy balance consideration that the effective, or intergranular, stress is the appropriate measure of stress for describing the constitutive response of the soil skeleton since it absorbs all the strain energy generated in the saturated soil-water mixture. Finally, it is shown that the mathematical model is amenable to consistent linearization, and that explicit expressions for the consistent tangent operators can be derived for use in numerical solutions such as those based on the finite element method.

Dynamic response of underpasses for high-speed train lines

Underpasses are common in modern railway lines. Wildlife corridors and drainage conduits often fall into this category of partially buried structures. Their dynamic behavior has received far less attention than that of other structures such as bridges, but their large number makes their study an interesting challenge from the viewpoint of safety and cost savings. Here, we present a complete study of a culvert, including on-site measurements and numerical modeling. The studied structure belongs to the high-speed railway line linking Segovia and Valladolid in Spain. The line was opened to traffic in 2004. On-site measurements were performed for the structure by recording the dynamic response at selected points of the structure during the passage of high-speed trains at speeds ranging between 200 and 300 km/h. The measurements provide not only reference values suitable for model fitting, but also a good insight into the main features of the dynamic behavior of this structure. Finite element techniques were used to model the dynamic behavior of the structure and its key features. Special attention is paid to vertical accelerations, the values of which should be limited to avoid track instability according to Eurocode. This study furthers our understanding of the dynamic response of railway underpasses to train loads.

Nonlinear response to ultrasound of encapsulated microbubbles

The acoustic backscatter of encapsulated gas-filled microbubbles immersed in a weak compressible liquid and irradiated by ultrasound fields of moderate to high pressure amplitudes is investigated theoretically. The problem is formulated by considering, for the viscoelastic shell of finite thickness, an isotropic hyperelastic neo-Hookean model for the elastic contribution in addition to a Newtonian viscous component. First and second harmonic scattering cross-sections have been evaluated and the quantitative influence of the driving pressure amplitude on the harmonic resonance frequencies for different initial equilibrium bubble sizes and for different encapsulating physical properties has been determined. Conditions for optimal second harmonic imaging have been also investigated and some regions in the parameters space where the second harmonic intensity is dominant over the fundamental have been identified. Results have been obtained for albumin, lipid and polymer encapsulating shells, respectively.

Experimental and numerical study of bond response in structural concrete with corroded steel bars

Corrosion can affect the bond between reinforcing bars and concrete and hence the transfer of longitudinal stresses. Although a number of experimental studies on bond failure have been conducted in recent years, the findings have diverged rather widely, due primarily to differing test conditions. The present paper reports on an experimental programme consisting of eccentric pull-out tests run on corroded steel bars in specimens subjected to accelerated or natural corrosion. An axisymmetric bi-dimensional FE model with finite deformations initially developed to study bond mechanics with sound steel bars, has been enhanced to consider bond effects in corroded steel bars. The model simulation is compared to some of the experimental results for corroded and sound bars and the findings are analysed.

A numerical study of the influence of microvoids in the transverse mechanical response of unidirectional composites

The effect of porosity on the transverse mechanical properties of unidirectional fiber-reinforced composites is studied by means of computational micromechanics. The composite behavior is simulated by the finite element analysis of a representative volume element of the composite microstructure in which the random distribution of fibers and the voids are explicitly included. Two types of voids – interfiber voids and matrix voids – were included in the microstructure and the actual damage mechanisms in the composite, namely matrix and interface failure, were accounted for. It was found that porosity (in the range 1–5%) led to a large reduction in the transverse strength and the influence of both types of voids in the onset and propagation of damage throughout the microstructure was studied under transverse tension and compression. Finally, the failure locus of the composite lamina under transverse tension/compression and out-of-plane shear was obtained by means of computational micromechanics and compared with the predictions of Puck’s model and with experimental data available in the literature. The results show that the strength of composites is significantly reduced by the presence of voids

A methodology to calibrate structural finite element models for reinforced concrete structures subject to blast loads

Civil buildings are not specifically designed to support blast loads, but it is important to take into account these potential scenarios because of their catastrophic effects, on persons and structures. A practical way to consider explosions on reinforced concrete structures is necessary. With this objective we propose a methodology to evaluate blast loads on large concrete buildings, using LS-DYNA code for calculation, with Lagrangian finite elements and explicit time integration. The methodology has three steps. First, individual structural elements of the building like columns and slabs are studied, using continuum 3D elements models subjected to blast loads. In these models reinforced concrete is represented with high precision, using advanced material models such as CSCM_CONCRETE model, and segregated rebars constrained within the continuum mesh. Regrettably this approach cannot be used for large structures because of its excessive computational cost. Second, models based on structural elements are developed, using shells and beam elements. In these models concrete is represented using CONCRETE_EC2 model and segregated rebars with offset formulation, being calibrated with continuum elements models from step one to obtain the same structural response: displacement, velocity, acceleration, damage and erosion. Third, models basedon structural elements are used to develop large models of complete buildings. They are used to study the global response of buildings subjected to blast loads and progressive collapse. This article carries out different techniques needed to calibrate properly the models based on structural elements, using shells and beam elements, in order to provide results of sufficient accuracy that can be used with moderate computational cost.

A gradient-enhanced continuum damage model with application to fibre-reinforced tissues at finite strains

A non-local gradient-based damage formulation within a geometrically non-linear set- ting is presented. The hyperelastic constitutive response at local material point level is governed by a strain energy function which is additively composed by an isotropic neo-Hookean matrix and by an anisotropic fibre-reinforced material based on the model proposed by T. Gasser, R. Ogden, and G. Holzapfel.

Forced Response of Mistuned Bladed Disks: Quantitative Validation of the Asymptotic Description

The effect of small mistuning in the forced response of a bladed disk is analyzed using a recently introduced methodology: the asymptotic mistuning model. The asymptotic mistuning model is an extremely reduced, simplified model that is derived directly from the full formulation of the mistuned bladed disk using a consistent perturbative procedure based on the relative smallness of the mistuning distortion. A detailed description of the derivation of the asymptotic mistuning model for a realistic bladed disk configuration is presented. The asymptotic mistuning model results for several different mistuning patterns and forcing conditions are compared with those from a high-resolution finite element model. The asymptotic mistuning model produces quantitatively accurate results, and, probably more relevant, it gives precise information about the factors (tuned modes and components of the mistuning pattern) that actually play a role in the vibrational forced response of mistuned bladed disks.

Influence of soil stiffness on the dynamic response of an offshore wind turbine

La industria de la energía eólica marina ha crecido de forma significativa durante los últimos 15 años, y se espera que siga creciendo durante los siguientes. La construcción de torres en aguas cada vez más profundas y el aumento en potencia y tamaño de las turbinas han creado la necesidad de diseñar estructuras de soporte cada vez más fiables y optimizadas, lo que requiere un profundo conocimiento de su comportamiento. Este trabajo se centra en la respuesta dinámica de una turbina marina con cimentación tipo monopilote y sobre la que actúa la fuerza del viento. Se han realizado cálculos con distintas propiedades del suelo para cubrir un rango de rigideces que va desde una arena muy suelta a una muy densa. De este modo se ha analizado la influencia que tiene la rigidez del suelo en el comportamiento de la estructura. Se han llevado a cabo análisis estáticos y dinámicos en un modelo de elementos finitos implementado en Abaqus. El desplazamiento en la cabeza de la torre y la tensión en su base se han obtenido en función de la rigidez del suelo, y con ellos se ha calculado la amplificación dinámica producida cuando la frecuencia natural del sistema suelo‐cimentación torre se aproxima a la frecuencia de la carga. Dos diferentes enfoques a la hora de modelizar el suelo se han comparado: uno utilizando elementos continuos y otro utilizando muelles elásticos no lineales. Por último, un análisis de fiabilidad se ha llevado a cabo con un modelo analítico para calcular la probabilidad de resonancia del sistema, en el que se han considerado las propiedades de rigidez del suelo como variables aleatorias. Offshore wind energy industry has experienced a significant growth over the past 15 years, and it is expected to continue its growth in the coming years. The expansion to increasingly deep waters and the rise in power and size of the turbines have led to a need for more reliable and optimized support designs, which requires an extensive knowledge of the behaviour of these structures. This work focuses on the dynamic response of an offshore wind turbine founded on a monopile and subjected to wind loading. Different soil properties have been considered in order to cover the range of stiffness from a very loose to a very dense sand. In this way, the influence of stiffness on the structure behaviour has been assessed. Static and dynamic analyses have been carried out by means of a finite element model implemented in Abaqus. Head displacement and stress at the tower base have been obtained as functions of soil stiffness, and they have been used to calculate the dynamic amplification that is produced when the natural frequency of the system soil‐foundation‐tower approaches the load frequency. Two different approaches of soil modelling have been compared: soil modelled as a continuum and soil simulated with non linear elastic springs. Finally, a reliability analysis to assess the probability of resonance has been performed with an analytical model, in which soil stiffness properties are considered as stochastic variables.

A study of the mechanical behaviour of the human aortic artery by means of non-linear finite element models

En la presente tesis desarrollamos una estrategia para la simulación numérica del comportamiento mecánico de la aorta humana usando modelos de elementos finitos no lineales. Prestamos especial atención a tres aspectos claves relacionados con la biomecánica de los tejidos blandos. Primero, el análisis del comportamiento anisótropo característico de los tejidos blandos debido a las familias de fibras de colágeno. Segundo, el análisis del ablandamiento presentado por los vasos sanguíneos cuando estos soportan cargas fuera del rango de funcionamiento fisiológico. Y finalmente, la inclusión de las tensiones residuales en las simulaciones en concordancia con el experimento de apertura de ángulo. El análisis del daño se aborda mediante dos aproximaciones diferentes. En la primera aproximación se presenta una formulación de daño local con regularización. Esta formulación tiene dos ingredientes principales. Por una parte, usa los principios de la teoría de la fisura difusa para garantizar la objetividad de los resultados con diferentes mallas. Por otra parte, usa el modelo bidimensional de Hodge-Petruska para describir el comportamiento mesoscópico de los fibriles. Partiendo de este modelo mesoscópico, las propiedades macroscópicas de las fibras de colágeno son obtenidas a través de un proceso de homogenización. En la segunda aproximación se presenta un modelo de daño no-local enriquecido con el gradiente de la variable de daño. El modelo se construye a partir del enriquecimiento de la función de energía con un término que contiene el gradiente material de la variable de daño no-local. La inclusión de este término asegura una regularización implícita de la implementación por elementos finitos, dando lugar a resultados de las simulaciones que no dependen de la malla. La aplicabilidad de este último modelo a problemas de biomecánica se estudia por medio de una simulación de un procedimiento quirúrgico típico conocido como angioplastia de balón. In the present thesis we develop a framework for the numerical simulation of the mechanical behaviour of the human aorta using non-linear finite element models. Special attention is paid to three key aspects related to the biomechanics of soft tissues. First, the modelling of the characteristic anisotropic behaviour of the softue due to the collagen fibre families. Secondly, the modelling of damage-related softening that blood vessels exhibit when subjected to loads beyond their physiological range. And finally, the inclusion of the residual stresses in the simulations in accordance with the opening-angle experiment The modelling of damage is addressed with two major and different approaches. In the first approach a continuum local damage formulation with regularisation is presented. This formulation has two principal ingredients. On the one hand, it makes use of the principles of the smeared crack theory to avoid the mesh size dependence of the structural response in softening. On the other hand, it uses a Hodge-Petruska bidimensional model to describe the fibrils as staggered arrays of tropocollagen molecules, and from this mesoscopic model the macroscopic material properties of the collagen fibres are obtained using an homogenisation process. In the second approach a non-local gradient-enhanced damage formulation is introduced. The model is built around the enhancement of the free energy function by means of a term that contains the referential gradient of the non-local damage variable. The inclusion of this term ensures an implicit regularisation of the finite element implementation, yielding mesh-objective results of the simulations. The applicability of the later model to biomechanically-related problems is studied by means of the simulation of a typical surgical procedure, namely, the balloon angioplasty.

Resonadores piezoeléctricos de cizalladura: estudio de la formación de biofilms y caracterización de fluidos magneto‐reológicos

La presente tesis revisa y analiza algunos aspectos fundamentales relativos al comportamiento de los sensores basados en resonadores piezoeléctricos TSM (Thickness Shear Mode), así como la aplicación de los mismos al estudio y caracterización de dos medios viscoelásticos de gran interés: los fluidos magnetoreológicos y los biofilms microbianos. El funcionamiento de estos sensores está basado en la medida de sus propiedades resonantes, las cuales varían al entrar en contacto con el material que se quiere analizar. Se ha realizado un análisis multifrecuencial, trabajando en varios modos de resonancia del transductor, en algunas aplicaciones incluso de forma simultánea (excitación pulsada). Se han revisado fenómenos como la presencia de microcontactos en la superficie del sensor y la resonancia de capas viscoelásticas de espesor finito, que pueden afectar a los sensores de cuarzo de manera contraria a lo que predice la teoría convencional (Sauerbrey y Kanazawa), pudiéndonos llevar a incrementos positivos de la frecuencia de resonancia. Además, se ha estudiado el efecto de una deposición no uniforme sobre el resonador piezoeléctrico. Para ello se han medido deposiciones de poliuretano, modelándose la respuesta del resonador con estas deposiciones mediante FEM. El modelo numérico permite estudiar el comportamiento del resonador al modificar distintas variables geométricas (espesor, superficie, no uniformidad y zona de deposición) de la capa depositada. Se ha demostrado que para espesores de entre un cuarto y media longitud de onda aproximadamente, una capa viscoelástica no uniforme sobre la superficie del sensor, amplifica el incremento positivo del desplazamiento de la frecuencia de resonancia en relación con una capa uniforme. Se ha analizado también el patrón geométrico de la sensibilidad del sensor, siendo también no uniforme sobre su superficie. Se han aplicado sensores TSM para estudiar los cambios viscoelásticos que se producen en varios fluidos magneto-reológicos (FMR) al aplicarles distintos esfuerzos de cizalla controlados por un reómetro. Se ha podido ver que existe una relación directa entre diversos parámetros reológicos obtenidos con el reómetro (fuerza normal, G’, G’’, velocidad de deformación, esfuerzo de cizalla…) y los parámetros acústicos, caracterizándose los FMR tanto en ausencia de campo magnético, como con campo magnético aplicado a distintas intensidades. Se han estudiado las ventajas que aporta esta técnica de medida sobre la técnica basada en un reómetro comercial, destacando que se consigue caracterizar con mayor detalle algunos aspectos relevantes del fluido como son la deposición de partículas (estabilidad del fluido), el proceso de ruptura de las estructuras formadas en los FMR tanto en presencia como en ausencia de campo magnético y la rigidez de los microcontactos que aparecen entre partículas y superficies. También se han utilizado sensores de cuarzo para monitorear en tiempo real la formación de biofilms de Staphylococcus epidermidis y Eschericia coli sobre los propios resonadores de cristal de cuarzo sin ningún tipo de recubrimiento, realizándose ensayos con cepas que presentan distinta capacidad de producir biofilm. Se mostró que, una vez que se ha producido una primera adhesión homogénea de las bacterias al sustrato, podemos considerar el biofilm como una capa semi-infinita, de la cual el sensor de cuarzo refleja las propiedades viscoelásticas de la región inmediatamente contigua al resonador, no siendo sensible a lo que sucede en estratos superiores del biofilm. Los experimentos han permitido caracterizar el módulo de rigidez complejo de los biofilms a varias frecuencias, mostrándose que el parámetro característico que indica la adhesión de un biofilm tanto en el caso de S. epidermidis como de E. coli, es el incremento de G’ (relacionado con la elasticidad o rigidez de la capa), el cual viene ligado a un incremento de la frecuencia de resonancia del sensor. ABSTRACT This thesis reviews and analyzes some key aspects of the behavior of sensors based on piezoelectric resonators TSM (Thickness Shear Mode) and their applications to the study and characterization in two viscoelastic media of great interest: magnetorheological fluids and microbial biofilms. The operation of these sensors is based on the analysis of their resonant properties that vary in contact with the material to be analyzed. We have made a multi-frequency analysis, working in several modes of resonance of the transducer, in some applications even simultaneously (by impulse excitation). We reviewed some phenomena as the presence of micro-contacts on the sensor surface and the resonance of viscoelastic layers of finite thickness, which can affect quartz sensors contrary to the conventional theory predictions (Sauerbrey and Kanazawa), leading to positive resonant frequency shifts. In addition, we studied the effect of non-uniform deposition on the piezoelectric resonator. Polyurethane stools have been measured, being the resonator response to these depositions modeled by FEM. The numerical model allows studying the behavior of the resonator when different geometric variables (thickness, surface non-uniformity and deposition zone) of the deposited layer are modified. It has been shown that for thicknesses between a quarter and a half of a wavelength approximately, non-uniform deposits on the sensor surface amplify the positive increase of the resonance frequency displacement compared to a uniform layer. The geometric pattern of the sensor sensitivity was also analyzed, being also non-uniform over its surface. TSM sensors have been applied to study the viscoelastic changes occurring in various magneto-rheological fluids (FMR) when subjected to different controlled shear stresses driven by a rheometer. It has been seen that there is a direct relationship between various rheological parameters obtained with the rheometer (normal force, G', G'', stress, shear rate ...) and the acoustic parameters, being the FMR characterized both in the absence of magnetic field, and when the magnetic field was applied at different intensities. We have studied the advantages of this technique over the characterization methods based on commercial rheometers, noting that TSM sensors are more sensitive to some relevant aspects of the fluid as the deposition of particles (fluid stability), the breaking process of the structures formed in the FMR both in the presence and absence of magnetic field, and the rigidity of the micro-contacts appearing between particles and surfaces. TSM sensors have also been used to monitor in real time the formation of biofilms of Staphylococcus epidermidis and Escherichia coli on the quartz crystal resonators themselves without any coating, performing tests with strains having different ability to produce biofilm. It was shown that, once a first homogeneous adhesion of bacteria was produced on the substrate, the biofilm can be considered as a semi-infinite layer and the quartz sensor reflects only the viscoelastic properties of the region immediately adjacent to the resonator, not being sensitive to what is happening in upper layers of the biofilm. The experiments allow the evaluation of the biofilm complex stiffness module at various frequencies, showing that the characteristic parameter that indicates the adhesion of a biofilm for the case of both S. epidermidis and E. coli, is an increased G' (related to the elasticity or stiffness of the layer), which is linked to an increase in the resonance frequency of the sensor.

Endothelin-induced conversion of embryonic heart muscle cells into impulse-conducting Purkinje fibers

A regular heart beat is dependent on a specialized network of pacemaking and conductive cells. There has been a longstanding controversy regarding the developmental origin of these cardiac tissues which also manifest neural-like properties. Recently, we have shown conclusively that during chicken embryogenesis, impulse-conducting Purkinje cells are recruited from myocytes in spatial association with developing coronary arteries. Here, we report that cultured embryonic myocytes convert to a Purkinje cell phenotype after exposure to the vascular cytokine, endothelin. This inductive response declined gradually during development. These results yield further evidence for a role of arteriogenesis in the induction of impulse-conducting Purkinje cells within the heart muscle lineage and also may provide a basis for tissue engineering of cardiac pacemaking and conductive cells.

Coupled Biomechanical Response of the Cornea Assessed by Non-Contact Tonometry. A Simulation Study

The mechanical response of the cornea subjected to a non-contact air-jet tonometry diagnostic test represents an interplay between its geometry, the corneal material behavior and the loading. The objective is to study this interplay to better understand and interpret the results obtained with a non-contact tonometry test. A patient-specific finite element model of a healthy eye, accounting for the load free configuration, was used. The corneal tissue was modeled as an anisotropic hyperelastic material with two preferential directions. Three different sets of parameters within the human experimental range obtained from inflation tests were considered. The influence of the IOP was studied by considering four pressure levels (10–28 mmHg) whereas the influence of corneal thickness was studied by inducing a uniform variation (300–600 microns). A Computer Fluid Dynamics (CFD) air-jet simulation determined pressure loading exerted on the anterior corneal surface. The maximum apex displacement showed a linear variation with IOP for all materials examined. On the contrary, the maximum apex displacement followed a cubic relation with corneal thickness. In addition, a significant sensitivity of the apical displacement to the corneal stiffness was also obtained. Explanation to this behavior was found in the fact that the cornea experiences bending when subjected to an air-puff loading, causing the anterior surface to work in compression whereas the posterior surface works in tension. Hence, collagen fibers located at the anterior surface do not contribute to load bearing. Non-contact tonometry devices give useful information that could be misleading since the corneal deformation is the result of the interaction between the mechanical properties, IOP, and geometry. Therefore, a non-contact tonometry test is not sufficient to evaluate their individual contribution and a complete in-vivo characterization would require more than one test to independently determine the membrane and bending corneal behavior.

NONSAP : a structural analysis program for static and dynamic response of nonlinear systems /

Photocopy of original: Berkeley : Structural Engineering Laboratory, University of California, 1974.

Adaptation is not required to explain the long-term response of axons to molecular gradients

It has been suggested that growth cones navigating through the developing nervous system might display adaptation, so that their response to gradient signals is conserved over wide variations in ligand concentration. Recently however, a new chemotaxis assay that allows the effect of gradient parameters on axonal trajectories to be finely varied has revealed a decline in gradient sensitivity on either side of an optimal concentration. We show that this behavior can be quantitatively reproduced with a computational model of axonal chemotaxis that does not employ explicit adaptation. Two crucial components of this model required to reproduce the observed sensitivity are spatial and temporal averaging. These can be interpreted as corresponding, respectively, to the spatial spread of signaling effects downstream from receptor binding, and to the finite time over which these signaling effects decay. For spatial averaging, the model predicts that an effective range of roughly one-third of the extent of the growth cone is optimal for detecting small gradient signals. For temporal decay, a timescale of about 3 minutes is required for the model to reproduce the experimentally observed sensitivity.