Dynamical Elastic Moduli of the Ti-13Nb-13Zr biomaterial alloy by mechanical spectroscopy

Dynamical Elastic Moduli of the Ti-13Nb-13Zr biomaterial alloy were obtained using the mechanical spectroscopy technique. The sample with heat treatment at 1170K for 30 minutes and water quenched with subsequent aging treatment at 670 K for 3 hours (TNZ + WQ + 670 K/3 h), was characterized through decay of free oscillations of the sample in the flexural vibration mode. The spectra of anelastic relaxation (internal friction and frequency) in the temperature range from 300 K to 625 K not revealed the presence of relaxation process. As shown in the literature, the hcp structure usually does not exhibit any relaxation due to the symmetry of the sites in the crystalline lattice, but if there is some relaxation, this only occurs in special cases such as low concentration of zirconium or saturation of the stoichiometric ratio of oxygen for zirconium. Dynamical elastic modulus obtained for TNZ + WQ + 670 K/3 h alloy was 87 GPa at room temperature, which is higher than the value for Ti-13Nb-13Zr alloy (64 GPa) of the literature. This increment may be related to the change of the proportion of α and β phases. Besides that, the presence of precipitates in the alloy after aging treatment hardens the material and reduces its ductility.

Parametric analysis of effective material properties of thickness-shear piezoelectric macro-fibre composites

A previous study on the characterization of effective material properties of a d15 thickness-shear piezoelectric Macro-Fibre Composite (MFC) made of seven layers (Kapton, Acrylic, Electrode, Piezoceramic Fibre and Epoxy Composite, Electrode, Acrylic, Kapton) using a finite element homogenization method has shown that the packaging reduces significantly the shear stiffness of the piezoceramic material and, thus, leads to significantly smaller effective electromechanical coupling coefficient k15 and piezoelectric stress constant e15 when compared to the piezoceramic fibre properties. Therefore, the main objective of this work is to perform a parametric analysis in which the effect of the variations of fibre volume fraction, Epoxy elastic modulus, electrode thickness and active layer thickness on the MFC effective material properties is evaluated. Results indicate that an effective d15 MFC should use relatively thick fibres having relatively high shear modulus and relatively stiff epoxy filler. On the other hand, the electrode thickness does not affect significantly the MFC performance.

Elastic Propagation in random media: applications to the imaging of volcano structures

High-frequency seismograms contain features that reflect the random inhomogeneities of the earth. In this work I use an imaging method to locate the high contrast small- scale heterogeneity respect to the background earth medium. This method was first introduced by Nishigami (1991) and than applied to different volcanic and tectonically active areas (Nishigami, 1997, Nishigami, 2000, Nishigami, 2006). The scattering imaging method is applied to two volcanic areas: Campi Flegrei and Mt. Vesuvius. Volcanic and seismological active areas are often characterized by complex velocity structures, due to the presence of rocks with different elastic properties. I introduce some modifications to the original method in order to make it suitable for small and highly complex media. In particular, for very complex media the single scattering approximation assumed by Nishigami (1991) is not applicable as the mean free path becomes short. The multiple scattering or diffusive approximation become closer to the reality. In this thesis, differently from the ordinary Nishigami’s method (Nishigami, 1991), I use the mean of the recorded coda envelope as reference curve and calculate the variations from this average envelope. In this way I implicitly do not assume any particular scattering regime for the "average" scattered radiation, whereas I consider the variations as due to waves that are singularly scattered from the strongest heterogeneities. The imaging method is applied to a relatively small area (20 x 20 km), this choice being justified by the small length of the analyzed codas of the low magnitude earthquakes. I apply the unmodified Nishigami’s method to the volcanic area of Campi Flegrei and compare the results with the other tomographies done in the same area. The scattering images, obtained with frequency waves around 18 Hz, show the presence of high scatterers in correspondence with the submerged caldera rim in the southern part of the Pozzuoli bay. Strong scattering is also found below the Solfatara crater, characterized by the presence of densely fractured, fluid-filled rocks and by a strong thermal anomaly. The modified Nishigami’s technique is applied to the Mt. Vesuvius area. Results show a low scattering area just below the central cone and a high scattering area around it. The high scattering zone seems to be due to the contrast between the high rigidity body located beneath the crater and the low rigidity materials located around it. The central low scattering area overlaps the hydrothermal reservoirs located below the central cone. An interpretation of the results in terms of geological properties of the medium is also supplied, aiming to find a correspondence of the scattering properties and the geological nature of the material. A complementary result reported in this thesis is that the strong heterogeneity of the volcanic medium create a phenomenon called "coda localization". It has been verified that the shape of the seismograms recorded from the stations located at the top of the volcanic edifice of Mt. Vesuvius is different from the shape of the seismograms recorded at the bottom. This behavior is justified by the consideration that the coda energy is not uniformly distributed within a region surrounding the source for great lapse time.

Hypersonic elastic excitations in soft mesoscopic structures

This dissertation is devoted to the experimental exploration of the propagation of elastic waves in soft mesoscopic structures with submicrometer dimensions. A strong motivation of this work is the large technological relevance and the fundamental importance of the subject. Elastic waves are accompanied by time-dependent fluctuations of local stress and strain fields in the medium. As such, the propagation phase velocities are intimately related to the elastic moduli. Knowledge of the elastic wave propagation directly provides information about the mechanical properties of the probed mesoscopic structures, which are not readily accessible experimentally. On the other hand, elastic waves, when propagating in an inhomogeneous medium with spatial inhomogeneities comparable to their wavelength, exhibit rather rich behavior, including the appearance of novel physical phenomena, such as phononic bandgap formation. So far, the experimental work has been restricted to macroscopic structures, which limit wave propagation below the KHz range. It was anticipated that an experimental approach capable of probing the interplay of the wave propagation with the controlled mesoscopic structures would contribute to deeper insights into the fundamental problem of elastic wave propagation in inhomogeneous systems. The mesoscopic nature of the structures to be studied precludes the use of traditional methods, such as sound transmission, for the study of elastic wave propagation. In this work, an optical method utilizing the inelastic scattering of photons by GHz frequency thermally excited elastic waves, known as Brillouin light scattering spectroscopy (BLS), was employed. Two important classes of soft structures were investigated: thin films and colloidal crystals. For the former, the main interest was the effect of the one-dimensional (1D) confinement on the wave propagation due to the presence of the free-surface or interface of the layer and the utilization of these waves to extract relevant material parameters. For the second system, the primary interest was the interaction of the elastic wave and the strong scattering medium with local resonance units in a three-dimensional (3D) periodic arrangement.

Radiation-induced defects in calcium fluoride and their influence on material properties under 193 nm laser irradiation

Calcium fluoride (CaF2) is one of the key lens materials in deep-ultraviolet microlithography because of its transparency at 193 nm and its nearly perfect optical isotropy. Its physical and chemical properties make it applicable for lens fabrication. The key feature of CaF2 is its extreme laser stability. rnAfter exposing CaF2 to 193 nm laser irradiation at high fluences, a loss in optical performance is observed, which is related to radiation-induced defect structures in the material. The initial rapid damage process is well understood as the formation of radiation-induced point defects, however, after a long irradiation time of up to 2 months, permanent damage of the crystals is observed. Based on experimental results, these permanent radiation-induced defect structures are identified as metallic Ca colloids.rnThe properties of point defects in CaF2 and their stabilization in the crystal bulk are calculated with density functional theory (DFT). Because the stabilization of the point defects and the formation of metallic Ca colloids are diffusion-driven processes, the diffusion coefficients for the vacancy (F center) and the interstitial (H center) in CaF2 are determined with the nudged elastic band method. The optical properties of Ca colloids in CaF2 are obtained from Mie-theory, and their formation energy is determined.rnBased on experimental observations and the theoretical description of radiation-induced point defects and defect structures, a diffusion-based model for laser-induced material damage in CaF2 is proposed, which also includes a mechanism for annealing of laser damage. rn

Measuring the elastic modulus of polymers using the atomic force microscope

Testing a new method of nanoindentation using the atomic force microscope (AFM) was the purpose of this research. Nanoindentation is a useful technique to study the properties of materials on the sub-micron scale. The AFM has been used as a nanoindenter previously; however several parameters needed to obtain accurate results, including tip radius and cantilever sensitivity, can be difficult to determine. To solve this problem, a new method to determine the elastic modulus of a material using the atomic force microscope (AFM) has been proposed by Tang et al. This method models the cantilever and the sample as two springs in a series. The ratio of the cantilever spring constant (k) to diameter of the tip (2a) is treated in the model as one parameter (α=k/2a). The value of a, along with the cantilever sensitivity, are determined on two reference samples with known mechanical properties and then used to find the elastic modulus of an unknown sample. To determine the reliability and accuracy of this technique, it was tested on several polymers. Traditional depth-sensing nanoindentation was preformed for comparison. The elastic modulus values from the AFM were shown to be statistically similar to the nanoindenter results for three of the five samples tested.

An anisotropic elastic-viscoplastic damage model for bone tissue

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain D(κ) , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.

High-Pressure Behavior and Phase Stability of Al5BO9, a Mullite-Type Ceramic Material

Phase stability, elastic behavior, and pressure-induced structural evolution of synthetic boron-mullite Al5BO9 (a = 5.6780(7), b = 15.035(6), and c =7.698(3) Å, space group Cmc21, Z = 4) were investigated up to 25.6(1) GPa by in situ single-crystal synchrotron X-ray diffraction with a diamond anvil cell (DAC) under hydrostatic conditions. No evidence of phase transition was observed up to 21.7(1) GPa. At 25.6(1) GPa, the refined unit-cell parameters deviated significantly from the compressional trend, and the diffraction peaks appeared broader than at lower pressure. At 26.7(1) GPa, the diffraction pattern was not indexable, suggesting amorphization of the material or a phase transition to a high-pressure polymorph. Fitting the P–V data up to 21.7(1) GPa with a second-order Birch–Murnaghan Equation-of-State, we obtained a bulk modulus KT0 = 164(1) GPa. The axial compressibilities, here described as linearized bulk moduli, are as follows: KT0(a) = 244(9), KT0(b) = 120(4), and KT0(c) = 166(11) GPa (KT0(a):KT0(b):KT0(c) = 2.03:1:1.38). The structure refinements allowed a description of the main deformation mechanisms in response to the applied pressure. The stiffer crystallographic direction appears to be controlled by the infinite chains of edge-sharing octahedra running along [100], making the structure less compressible along the a-axis than along the b- and c-axis.

Ion Beam irradiation of copper nitride: electronic vs elastic-collision mechanism

Copper nitride is a metastable material which results very attractive because of their potential to be used in functional device. Cu3 N easily decomposes into Cu and N2 by annealing [1] or irradiation (electron, ions, laser) [2, 3]. Previous studies carried out in N-rich Cu3 N films irradiated with Cu at 42MeV evidence a very efficient sputtering of N whose yield (5×10 3 atom/ion), for a film with a thickness of just 100 nm, suggest that the origin of the sputtering has an electronic nature. This N depletion was observed to be responsible for new phase formation ( Cu2 O) and pure Cu [4]

Discontinuous solutions and boundary value problems for non-linearly elastic fibre-reinforced materials

En este trabajo se han analizado varios problemas en el contexto de la elasticidad no lineal basándose en modelos constitutivos representativos. En particular, se han analizado problemas relacionados con el fenómeno de perdida de estabilidad asociada con condiciones de contorno en el caso de material reforzados con fibras. Cada problema se ha formulado y se ha analizado por separado en diferentes capítulos. En primer lugar se ha mostrado el análisis del gradiente de deformación discontinuo para un material transversalmente isótropo, en particular, el modelo del material considerado consiste de una base neo-Hookeana isótropa incrustada con fibras de refuerzo direccional caracterizadas con un solo parámetro. La solución de este problema se vincula con instabilidades que dan lugar al mecanismo de fallo conocido como banda de cortante. La perdida de elipticidad de las ecuaciones diferenciales de equilibrio es una condición necesaria para que aparezca este tipo de soluciones y por tanto las inestabilidades asociadas. En segundo lugar se ha analizado una deformación combinada de extensión, inación y torsión de un tubo cilíndrico grueso donde se ha encontrado que la deformación citada anteriormente puede ser controlada solo para determinadas direcciones de las fibras refuerzo. Para entender el comportamiento elástico del tubo considerado se ha ilustrado numéricamente los resultados obtenidos para las direcciones admisibles de las fibras de refuerzo bajo la deformación considerada. En tercer lugar se ha estudiado el caso de un tubo cilíndrico grueso reforzado con dos familias de fibras sometido a cortante en la dirección azimutal para un modelo de refuerzo especial. En este problema se ha encontrado que las inestabilidades que aparecen en el material considerado están asociadas con lo que se llama soluciones múltiples de la ecuación diferencial de equilibrio. Se ha encontrado que el fenómeno de instabilidad ocurre en un estado de deformación previo al estado de deformación donde se pierde la elipticidad de la ecuación diferencial de equilibrio. También se ha demostrado que la condición de perdida de elipticidad y ^W=2 = 0 (la segunda derivada de la función de energía con respecto a la deformación) son dos condiciones necesarias para la existencia de soluciones múltiples. Finalmente, se ha analizado detalladamente en el contexto de elipticidad un problema de un tubo cilíndrico grueso sometido a una deformación combinada en las direcciones helicoidal, axial y radial para distintas geotermias de las fibras de refuerzo . In the present work four main problems have been addressed within the framework of non-linear elasticity based on representative constitutive models. Namely, problems related to the loss of stability phenomena associated with boundary value problems for fibre-reinforced materials. Each of the considered problems is formulated and analysed separately in different chapters. We first start with the analysis of discontinuous deformation gradients for a transversely isotropic material under plane deformation. In particular, the material model is an augmented neo-Hookean base with a simple unidirectional reinforcement characterised by a single parameter. The solution of this problem is related to material instabilities and it is associated with a shear band-type failure mode. The loss of ellipticity of the governing differential equations is a necessary condition for the existence of these material instabilities. The second problem involves a detailed analysis of the combined non-linear extension, inflation and torsion of a thick-walled circular cylindrical tube where it has been found that the aforementioned deformation is controllable only for certain preferred directions of transverse isotropy. Numerical results have been illustrated to understand the elastic behaviour of the tube for the admissible preferred directions under the considered deformation. The third problem deals with the analysis of a doubly fibre-reinforced thickwalled circular cylindrical tube undergoing pure azimuthal shear for a special class of the reinforcing model where multiple non-smooth solutions emerge. The associated instability phenomena are found to occur prior to the point where the nominal stress tensor changes monotonicity in a particular direction. It has been also shown that the loss of ellipticity condition that arises from the equilibrium equation and ^W=2 = 0 (the second derivative of the strain-energy function with respect to the deformation) are equivalent necessary conditions for the emergence of multiple solutions for the considered material. Finally, a detailed analysis in the basis of the loss of ellipticity of the governing differential equations for a combined helical, axial and radial elastic deformations of a fibre-reinforced circular cylindrical tube is carried out.

Optimización de la densidad de energía en vigas de material compuesto (PRF) sometidas a flexión pura

Las necesidades energéticas actuales requieren el desarrollo de tecnologías eficaces y eficientes en producción, transporte y distribución de energía. Estas necesidades han impulsado nuevos desarrollos en el ámbito energético, entre los cuales se encuentran sistemas de almacenamiento de energía. El avance en ingeniería de materiales permite pensar en la posibilidad de almacenamiento mediante deformación elástica de vigas. Concretamente se parte de un concepto de mecanismo acumulador de energía basado en la deformación elástica de resortes espirales a torsión. Dichos resortes se pueden considerar como elementos vigas sometidos a flexión pura y grandes deflexiones. Esta Tesis de centra en el diseño y optimización de estos elementos con el fin de maximizar la densidad de energía que son capaces de absorber. El proceso de optimización comienza con la identificación del factor crítico del que depende dicho proceso, en este caso de trata de la densidad de energía. Dicho factor depende de la geometría de la sección resistente y del material empleado en su construcción. En los últimos años ha existido un gran desarrollo de los materiales compuestos de tipo polimérico reforzados con fibras (PRF). Estos materiales están sustituyendo gradualmente a otros materiales, como los metales, debido principalmente a su excelente relación entre propiedades mecánicas y peso. Por otro lado, analizando las posibles geometrías para la sección resistente, se observó que la más adecuada es una estructura tipo sándwich. Se implementa así un procedimiento de diseño de vigas sándwich sometidas a flexión pura, con las pieles fabricadas en materiales compuestos tipo PRF y un núcleo que debe garantizar el bajo peso de la estructura. Se desarrolla así un procedimiento sistemático que se puede particularizar dependiendo de los parámetros de entrada de la viga, y que tiene en cuenta y analiza la aparición de todos los posibles modos de fallo posibles. Así mismo se desarrollan una serie de mapas o ábacos de diseño que permiten seleccionar rápidamente las dimensiones preliminares de la viga. Finalmente se llevan a cabo ensayos que permiten, por un lado, validar el concepto del mecanismo acumulador de energía a través del ensayo de un muelle con sección monolítica, y por otro validar los distintos diseños de vigas sándwich propuestos y mostrar el incremento de la densidad de energía con respecto a la alternativa monolítica. Como líneas futuras de investigación se plantean la investigación en nuevos materiales, como la utilización de nanotubos de carbono, y la optimización del mecanismo de absorción de energía; optimizando el mecanismo de absorción a flexión pura e implementando sistemas que permitan acumular energía mediante la deformación elástica debida a esfuerzos de tracción-compresión. ABSTRACT Energy supply requires the development of effective and efficient technologies for the production, transport and distribution of energy. In recent years, many energy storage systems have been developed. Advances in the field of materials engineering has allowed the development of new concepts as the energy storage by elastic deformation of beams. Particularly, in this Thesis an energy storage device based on the elastic deformation of torsional springs has been studied. These springs can be considered as beam elements subjected to pure bending loads and large deflections. This Thesis is focused on the design and optimization of these beam elements in order to maximize its density of stored energy. The optimization process starts with the identification of the critical factors for the elastic energy storage: the density. This factor depends on the geometry of the cross section of the beam and the materials from which it is made. In the last 20 years, major advances in the field of composite materials have been made, particularly in the field of fiber reinforced polymers (FRP). This type of material is substituting gradually metallic materials to their excellent weight-mechanical properties ratio. In the other side, several possible geometries are analyzed for its use in the cross section of the beam; it was concluded that the best option, for maximum energy density, is using a sandwich beam. A design procedure for sandwich beams with skins made up with FRP composites and a light weight core is developed. This procedure can be particularized for different input parameters and it analyzes all the possible failure modes. Abacus and failure mode maps have been developed in order to simplify the design process. Finally several tested was made. Firstly, a prototype of the energy storage system which uses a monolithic composite beam was tested in order to validate the concept of the energy storage by elastic deformation. After that sandwich beam samples are built and tested, validating the design and showing the increase of energy density with respect to the monolithic beam. As futures research lines the following are proposed: research in new materials, as carbon nanotubes; and the optimization of the energy storage mechanism, that means optimizing the pure bending storage mechanism and developing new ones based on traction-compression mechanisms.

Métodos computacionales para visco-hiperelasticidad anisótropa en grandes deformaciones

El comportamiento mecánico de muchos materiales biológicos y poliméricos en grandes deformaciones se puede describir adecuadamente mediante formulaciones isocóricas hiperelásticas y viscoelásticas. Las ecuaciones de comportamiento elástico y viscoelástico y las formulaciones computacionales para materiales incompresibles isótropos en deformaciones finitas están ampliamente desarrolladas en la actualidad. Sin embargo, el desarrollo de modelos anisótropos no lineales y de sus correspondientes formulaciones computacionales sigue siendo un tema de investigación de gran interés. Cuando se consideran grandes deformaciones, existen muchas medidas de deformación disponibles con las que poder formular las ecuaciones de comportamiento. Los modelos en deformaciones cuadráticas facilitan la implementación en códigos de elementos finitos, ya que estas medidas surgen de forma natural en la formulación. No obstante, pueden dificultar la interpretación de los modelos y llevar a resultados pocos realistas. El uso de deformaciones logarítmicas permite el desarrollo de modelos más simples e intuitivos, aunque su formulación computacional debe ser adaptada a las exigencias del programa. Como punto de partida, en esta tesis se demuestra que las deformaciones logarítmicas representan la extensión natural de las deformaciones infinitesimales, tanto axiales como angulares, al campo de las grandes deformaciones. Este hecho permite explicar la simplicidad de las ecuaciones resultantes. Los modelos hiperelásticos predominantes en la actualidad están formulados en invariantes de deformaciones cuadráticas. Estos modelos, ya sean continuos o microestructurales, se caracterizan por tener una forma analítica predefinida. Su expresión definitiva se calcula mediante un ajuste de curvas a datos experimentales. Un modelo que no sigue esta metodología fue desarrollado por Sussman y Bathe. El modelo es sólo válido para isotropía y queda definido por una función de energía interpolada con splines, la cual reproduce los datos experimentales de forma exacta. En esta tesis se presenta su extensión a materiales transversalmente isótropos y ortótropos utilizando deformaciones logarítmicas. Asimismo, se define una nueva propiedad que las funciones de energía anisótropas deben satisfacer para que su convergencia al caso isótropo sea correcta. En visco-hiperelasticidad, aparte de las distintas funciones de energía disponibles, hay dos aproximaciones computational típicas basadas en variables internas. El modelo original de Simó está formulado en tensiones y es válido para materiales anisótropos, aunque sólo es adecuado para pequeñas desviaciones con respecto al equilibrio termodinámico. En cambio, el modelo basado en deformaciones de Reese y Govindjee permite grandes deformaciones no equilibradas pero es, en esencia, isótropo. Las formulaciones anisótropas en este último contexto son microestructurales y emplean el modelo isótropo para cada uno de los constituyentes. En esta tesis se presentan dos formulaciones fenomenológicas viscoelásticas definidas mediante funciones hiperelásticas anisótropas y válidas para grandes desviaciones con respecto al equilibrio termodinámico. El primero de los modelos está basado en la descomposición multiplicativa de Sidoroff y requiere un comportamiento viscoso isótropo. La formulación converge al modelo de Reese y Govindjee en el caso especial de isotropía elástica. El segundo modelo se define a partir de una descomposición multiplicativa inversa. Esta formulación está basada en una descripción co-rotacional del problema, es sustancialmente más compleja y puede dar lugar a tensores constitutivos ligeramente no simétricos. Sin embargo, su rango de aplicación es mucho mayor ya que permite un comportamiento anisótropo tanto elástico como viscoso. Varias simulaciones de elementos finitos muestran la gran versatilidad de estos modelos cuando se combinan con funciones hiperelásticas formadas por splines. ABSTRACT The mechanical behavior of many polymeric and biological materials may be properly modelled be means of isochoric hyperelastic and viscoelastic formulations. These materials may sustain large strains. The viscoelastic computational formulations for isotropic incompressible materials at large strains may be considered well established; for example Ogden’s hyperelastic function and the visco-hyperelastic model of Reese and Govindjee are well known models for isotropy. However, anisotropic models and computational procedures both for hyperelasticity and viscohyperelasticity are still under substantial research. Anisotropic hyperelastic models are typically based on structural invariants obtained from quadratic strain measures. These models may be microstructurallybased or phenomenological continuum formulations, and are characterized by a predefined analytical shape of the stored energy. The actual final expression of the stored energy depends on some material parameters which are obtained from an optimization algorithm, typically the Levenberg-Marquardt algorithm. We present in this work anisotropic spline-based hyperelastic stored energies in which the shape of the stored energy is obtained as part of the procedure and which (exactly in practice) replicates the experimental data. These stored energies are based on invariants obtained from logarithmic strain measures. These strain measures preserve the metric and the physical meaning of the trace and deviator operators and, hence, are interesting and meaningful for anisotropic formulations. Furthermore, the proposed stored energies may be formulated in order to have material-symmetries congruency both from a theoretical and from a numerical point of view, which are new properties that we define in this work. On the other hand, visco-hyperelastic formulations for anisotropic materials are typically based on internal stress-like variables following a procedure used by Sim´o. However, it can be shown that this procedure is not adequate for large deviations from thermodynamic equilibrium. In contrast, a formulation given by Reese and Govindjee is valid for arbitrarily large deviations from thermodynamic equilibrium but not for anisotropic stored energy functions. In this work we present two formulations for visco-hyperelasticity valid for anisotropic stored energies and large deviations from thermodynamic equilibrium. One of the formulations is based on the Sidoroff multiplicative decomposition and converges to the Reese and Govindjee formulation for the case of isotropy. However, the formulation is restricted to isotropy for the viscous component. The second formulation is based on a reversed multiplicative decomposition. This last formulation is substantially more complex and based on a corotational description of the problem. It can also result in a slightly nonsymmetric tangent. However, the formulation allows for anisotropy not only in the equilibrated and non-equilibrated stored energies, but also in the viscous behavior. Some examples show finite element implementation, versatility and interesting characteristics of the models.

Determining the elastic-plastic properties of metallic materials through instrumented indentation

En los últimos años, y asociado al desarrollo de la tecnología MEMS, la técnica de indentación instrumentada se ha convertido en un método de ensayo no destructivo ampliamente utilizado para hallar las características elástico-plásticas de recubrimientos y capas delgadas, desde la escala macroscópica a la microscópica. Sin embargo, debido al complejo mecanismo de contacto debajo de la indentación, es urgente proponer un método más simple y conveniente para obtener unos resultados comparables con otras mediciones tradicionales. En este estudio, el objetivo es mejorar el procedimiento analítico para extraer las propiedades elástico-plásticas del material mediante la técnica de indentación instrumentada. La primera parte se centra en la metodología llevada a cabo para medir las propiedades elásticas de los materiales elásticos, presentándose una nueva metodología de indentación, basada en la evolución de la rigidez de contacto y en la curva fuerza-desplazamiento del ensayo de indentación. El método propuesto permite discriminar los valores de indentación experimental que pudieran estar afectados por el redondeo de la punta del indentador. Además, esta técnica parece ser robusta y permite obtener valores fiables del modulo elástico. La segunda parte se centra en el proceso analítico para determinar la curva tensión-deformación a partir del ensayo de indentación, empleando un indentador esférico. Para poder asemejar la curva tension-deformación de indentación con la que se obtendría de un ensayo de tracción, Tabor determinó empíricamente un factor de constricción de la tensión () y un factor de constricción de la deformación (). Sin embargo, la elección del valor de y  necesitan una derivación analítica. Se describió analíticamente una nueva visión de la relación entre los factores de constricción de tensión y la deformación basado en la deducción de la ecuación de Tabor. Un modelo de elementos finitos y un diseño experimental se realizan para evaluar estos factores de constricción. A partir de los resultados obtenidos, las curvas tension-deformación extraidas de los ensayos de indentación esférica, afectadas por los correspondientes factores de constricción de tension y deformación, se ajustaron a la curva nominal tensión-deformación obtenida de ensayos de tracción convencionales. En la última parte, se estudian las propiedades del revestimiento de cermet Inconel 625-Cr3C2 que es depositado en el medio de una aleación de acero mediante un láser. Las propiedades mecánicas de la matriz de cermet son estudiadas mediante la técnica de indentación instrumentada, haciendo uso de las metodologías propuestas en el presente trabajo. In recent years, along with the development of MEMS technology, instrumented indentation, as one type of a non-destructive measurement technique, is widely used to characterize the elastic and plastic properties of metallic materials from the macro to the micro scale. However, due to the complex contact mechanisms under the indentation tip, it is necessary to propose a more convenient and simple method of instrumented indention to obtain comparable results from other conventional measurements. In this study, the aim is to improve the analytical procedure for extracting the elastic plastic properties of metallic materials by instrumented indentation. The first part focuses on the methodology for measuring the elastic properties of metallic materials. An alternative instrumented indentation methodology is presented. Based on the evolution of the contact stiffness and indentation load versus the depth of penetration, the possibility of obtaining the actual elastic modulus of an elastic-plastic bulk material through instrumented sharp indentation tests has been explored. The proposed methodology allows correcting the effect of the rounding of the indenter tip on the experimental indentation data. Additionally, this technique does not seem too sensitive to the pile-up phenomenon and allows obtaining convincing values of the elastic modulus. In the second part, an analytical procedure is proposed to determine the representative stress-strain curve from the spherical indentation. Tabor has determined the stress constraint factor (stress CF), and strain constraint factor (strain CF), empirically but the choice of a value for and is debatable and lacks analytical derivation. A new insight into the relationship between stress and strain constraint factors is analytically described based on the formulation of Tabor’s equation. Finite element model and experimental tests have been carried out to evaluate these constraint factors. From the results, representative stress-strain curves using the proposed strain constraint factor fit better with the nominal stress-strain curve than those using Tabor’s constraint factors. In the last part, the mechanical properties of an Inconel 625-Cr3C2 cermet coating which is deposited onto a medium alloy steel by laser cladding has been studied. The elastic and plastic mechanical properties of the cermet matrix are studied using depth-sensing indentation (DSI) on the micro scale.

Relationship between static and dynamic elastic modulus of calcarenite heated at different temperatures: the San Julián's stone

The San Julián’s stone is the main material used to build the most important historical buildings in Alicante city (Spain). This paper describes the analysis developed to obtain the relationship between the static and the dynamic modulus of this sedimentary rock heated at different temperatures. The rock specimens have been subjected to heating processes at different temperatures to produce different levels of weathering on 24 specimens. The static and dynamic modulus has been measured for every specimen by means of the ISRM standard and ultrasonic tests, respectively. Finally, two analytic formulas are proposed for the relationship between the static and the dynamic modulus for this stone. The results have been compared with some relationships proposed by different researchers for other types of rock. The expressions presented in this paper can be useful for the analysis, using non-destructive techniques, of the integrity level of historical constructions built with San Julián’s stone affected by fires.

On the initiation of boudinage and folding instabilities in layered elasto-visco-plastic solids - Insights from natural microstructures and numerical simulations

The present understanding of the initiation of boudinage and folding structures is based on viscosity contrasts and stress exponents, considering an intrinsically unstable state of the layer. The criterion of localization is believed to be prescribed by geometry-material interactions, which are often encountered in natural structures. An alternative localization phenomenon has been established for ductile materials, in which instability emerges for critical material parameters and loading rates from homogeneous conditions. In this thesis, conditions are sought under which this type of instability prevails and whether localization in geological materials necessarily requires a trigger by geometric imperfections. The relevance of critical deformation conditions, material parameters and the spatial configuration of instabilities are discussed in a geological context. In order to analyze boudinage geometries, a numerical eigenmode analysis is introduced. This method allows determining natural frequencies and wavelengths of a structure and inducing perturbations on these frequencies. In the subsequent coupled thermo-mechanical simulations, using a grain size evolution and end-member flow laws, localization emerges when material softening through grain size sensitive viscous creep sets in. Pinch-and-swell structures evolve along slip lines through a positive feedback between the matrix response and material bifurcations inside the layer, independent from the mesh-discretization length scale. Since boudinage and folding are considered to express the same general instability, both structures should arise independently of the sign of the loading conditions and for identical material parameters. To this end, the link between material to energy instabilities is approached by means of bifurcation analyses of the field equations and finite element simulations of the coupled system of equations. Boudinage and folding structures develop at the same critical energy threshold, where dissipative work by temperature-sensitive creep overcomes the diffusive capacity of the layer. This finding provides basis for a unified theory for strain localization in layered ductile materials. The numerical simulations are compared to natural pinch-and-swell microstructures, tracing the adaption of grain sizes, textures and creep mechanisms in calcite veins. The switch from dislocation to diffusion creep relates to strain-rate weakening, which is induced by dissipated heat from grain size reduction, and marks the onset of continuous necking. The time-dependent sequence uncovers multiple steady states at different time intervals. Microstructurally and mechanically stable conditions are finally expressed in the pinch-and-swell end members. The major outcome of this study is that boudinage and folding can be described as the same coupled energy-mechanical bifurcation, or as one critical energy attractor. This finding allows the derivation of critical deformation conditions and fundamental material parameters directly from localized structures in the field.