Crack mechanical failure in lithium niobate crystal under ion irradiation; novel simulation by extended finite elements

Swift heavy ion irradiation (ions with mass heavier than 15 and energy exceeding MeV/amu) transfer their energy mainly to the electronic system with small momentum transfer per collision. Therefore, they produce linear regions (columnar nano-tracks) around the straight ion trajectory, with marked modifications with respect to the virgin material, e.g., phase transition, amorphization, compaction, changes in physical or chemical properties. In the case of crystalline materials the most distinctive feature of swift heavy ion irradiation is the production of amorphous tracks embedded in the crystal. Lithium niobate is a relevant optical material that presents birefringence due to its anysotropic trigonal structure. The amorphous phase is certainly isotropic. In addition, its refractive index exhibits high contrast with those of the crystalline phase. This allows one to fabricate waveguides by swift ion irradiation with important technological relevance. From the mechanical point of view, the inclusion of an amorphous nano-track (with a density 15% lower than that of the crystal) leads to the generation of important stress/strain fields around the track. Eventually these fields are the origin of crack formation with fatal consequences for the integrity of the samples and the viability of the method for nano-track formation. For certain crystal cuts (X and Y), these fields are clearly anisotropic due to the crystal anisotropy. We have used finite element methods to calculate the stress/strain fields that appear around the ion-generated amorphous nano-tracks for a variety of ion energies and doses. A very remarkable feature for X cut-samples is that the maximum shear stress appears on preferential planes that form +/-45º with respect to the crystallographic planes. This leads to the generation of oriented surface cracks when the dose increases. The growth of the cracks along the anisotropic crystal has been studied by means of novel extended finite element methods, which include cracks as discontinuities. In this way we can study how the length and depth of a crack evolves as function of the ion dose. In this work we will show how the simulations compare with experiments and their application in materials modification by ion irradiation.

Crack mechanical failure in ceramic materials under ion irradiation: case of lithium niobate crystal

Swift heavy ion irradiation (ions with mass heavier than 15 and energy exceeding MeV/amu) transfer their energy mainly to the electronic system with small momentum transfer per collision. Therefore, they produce linear regions (columnar nano-tracks) around the straight ion trajectory, with marked modifications with respect to the virgin material, e.g., phase transition, amorphization, compaction, changes in physical or chemical properties. In the case of crystalline materials the most distinctive feature of swift heavy ion irradiation is the production of amorphous tracks embedded in the crystal. Lithium niobate is a relevant optical material that presents birefringence due to its anysotropic trigonal structure. The amorphous phase is certainly isotropic. In addition, its refractive index exhibits high contrast with those of the crystalline phase. This allows one to fabricate waveguides by swift ion irradiation with important technological relevance. From the mechanical point of view, the inclusion of an amorphous nano-track (with a density 15% lower than that of the crystal) leads to the generation of important stress/strain fields around the track. Eventually these fields are the origin of crack formation with fatal consequences for the integrity of the samples and the viability of the method for nano-track formation. For certain crystal cuts (X and Y), these fields are clearly anisotropic due to the crystal anisotropy. We have used finite element methods to calculate the stress/strain fields that appear around the ion- generated amorphous nano-tracks for a variety of ion energies and doses. A very remarkable feature for X cut-samples is that the maximum shear stress appears on preferential planes that form +/-45º with respect to the crystallographic planes. This leads to the generation of oriented surface cracks when the dose increases. The growth of the cracks along the anisotropic crystal has been studied by means of novel extended finite element methods, which include cracks as discontinuities. In this way we can study how the length and depth of a crack evolves as function of the ion dose. In this work we will show how the simulations compare with experiments and their application in materials modification by ion irradiation.

Microstructure-based numerical modeling of the mechanical behavior of Mg alloys

Dentro de los materiales estructurales, el magnesio y sus aleaciones están siendo el foco de una de profunda investigación. Esta investigación está dirigida a comprender la relación existente entre la microestructura de las aleaciones de Mg y su comportamiento mecánico. El objetivo es optimizar las aleaciones actuales de magnesio a partir de su microestructura y diseñar nuevas aleaciones. Sin embargo, el efecto de los factores microestructurales (como la forma, el tamaño, la orientación de los precipitados y la morfología de los granos) en el comportamiento mecánico de estas aleaciones está todavía por descubrir. Para conocer mejor de la relación entre la microestructura y el comportamiento mecánico, es necesaria la combinación de técnicas avanzadas de caracterización experimental como de simulación numérica, a diferentes longitudes de escala. En lo que respecta a las técnicas de simulación numérica, la homogeneización policristalina es una herramienta muy útil para predecir la respuesta macroscópica a partir de la microestructura de un policristal (caracterizada por el tamaño, la forma y la distribución de orientaciones de los granos) y el comportamiento del monocristal. La descripción de la microestructura se lleva a cabo mediante modernas técnicas de caracterización (difracción de rayos X, difracción de electrones retrodispersados, así como con microscopia óptica y electrónica). Sin embargo, el comportamiento del cristal sigue siendo difícil de medir, especialmente en aleaciones de Mg, donde es muy complicado conocer el valor de los parámetros que controlan el comportamiento mecánico de los diferentes modos de deslizamiento y maclado. En la presente tesis se ha desarrollado una estrategia de homogeneización computacional para predecir el comportamiento de aleaciones de magnesio. El comportamiento de los policristales ha sido obtenido mediante la simulación por elementos finitos de un volumen representativo (RVE) de la microestructura, considerando la distribución real de formas y orientaciones de los granos. El comportamiento del cristal se ha simulado mediante un modelo de plasticidad cristalina que tiene en cuenta los diferentes mecanismos físicos de deformación, como el deslizamiento y el maclado. Finalmente, la obtención de los parámetros que controlan el comportamiento del cristal (tensiones críticas resueltas (CRSS) así como las tasas de endurecimiento para todos los modos de maclado y deslizamiento) se ha resuelto mediante la implementación de una metodología de optimización inversa, una de las principales aportaciones originales de este trabajo. La metodología inversa pretende, por medio del algoritmo de optimización de Levenberg-Marquardt, obtener el conjunto de parámetros que definen el comportamiento del monocristal y que mejor ajustan a un conjunto de ensayos macroscópicos independientes. Además de la implementación de la técnica, se han estudiado tanto la objetividad del metodología como la unicidad de la solución en función de la información experimental. La estrategia de optimización inversa se usó inicialmente para obtener el comportamiento cristalino de la aleación AZ31 de Mg, obtenida por laminado. Esta aleación tiene una marcada textura basal y una gran anisotropía plástica. El comportamiento de cada grano incluyó cuatro mecanismos de deformación diferentes: deslizamiento en los planos basal, prismático, piramidal hc+ai, junto con el maclado en tracción. La validez de los parámetros resultantes se validó mediante la capacidad del modelo policristalino para predecir ensayos macroscópicos independientes en diferentes direcciones. En segundo lugar se estudió mediante la misma estrategia, la influencia del contenido de Neodimio (Nd) en las propiedades de una aleación de Mg-Mn-Nd, obtenida por extrusión. Se encontró que la adición de Nd produce una progresiva isotropización del comportamiento macroscópico. El modelo mostró que este incremento de la isotropía macroscópica era debido tanto a la aleatoriedad de la textura inicial como al incremento de la isotropía del comportamiento del cristal, con valores similares de las CRSSs de los diferentes modos de deformación. Finalmente, el modelo se empleó para analizar el efecto de la temperatura en el comportamiento del cristal de la aleación de Mg-Mn-Nd. La introducción en el modelo de los efectos non-Schmid sobre el modo de deslizamiento piramidal hc+ai permitió capturar el comportamiento mecánico a temperaturas superiores a 150_C. Esta es la primera vez, de acuerdo con el conocimiento del autor, que los efectos non-Schmid han sido observados en una aleación de Magnesio. The study of Magnesium and its alloys is a hot research topic in structural materials. In particular, special attention is being paid in understanding the relationship between microstructure and mechanical behavior in order to optimize the current alloy microstructures and guide the design of new alloys. However, the particular effect of several microstructural factors (precipitate shape, size and orientation, grain morphology distribution, etc.) in the mechanical performance of a Mg alloy is still under study. The combination of advanced characterization techniques and modeling at several length scales is necessary to improve the understanding of the relation microstructure and mechanical behavior. Respect to the simulation techniques, polycrystalline homogenization is a very useful tool to predict the macroscopic response from polycrystalline microstructure (grain size, shape and orientation distributions) and crystal behavior. The microstructure description is fully covered with modern characterization techniques (X-ray diffraction, EBSD, optical and electronic microscopy). However, the mechanical behaviour of single crystals is not well-known, especially in Mg alloys where the correct parameterization of the mechanical behavior of the different slip/twin modes is a very difficult task. A computational homogenization framework for predicting the behavior of Magnesium alloys has been developed in this thesis. The polycrystalline behavior was obtained by means of the finite element simulation of a representative volume element (RVE) of the microstructure including the actual grain shape and orientation distributions. The crystal behavior for the grains was accounted for a crystal plasticity model which took into account the physical deformation mechanisms, e.g. slip and twinning. Finally, the problem of the parametrization of the crystal behavior (critical resolved shear stresses (CRSS) and strain hardening rates of all the slip and twinning modes) was obtained by the development of an inverse optimization methodology, one of the main original contributions of this thesis. The inverse methodology aims at finding, by means of the Levenberg-Marquardt optimization algorithm, the set of parameters defining crystal behavior that best fit a set of independent macroscopic tests. The objectivity of the method and the uniqueness of solution as function of the input information has been numerically studied. The inverse optimization strategy was first used to obtain the crystal behavior of a rolled polycrystalline AZ31 Mg alloy that showed a marked basal texture and a strong plastic anisotropy. Four different deformation mechanisms: basal, prismatic and pyramidal hc+ai slip, together with tensile twinning were included to characterize the single crystal behavior. The validity of the resulting parameters was proved by the ability of the polycrystalline model to predict independent macroscopic tests on different directions. Secondly, the influence of Neodymium (Nd) content on an extruded polycrystalline Mg-Mn-Nd alloy was studied using the same homogenization and optimization framework. The effect of Nd addition was a progressive isotropization of the macroscopic behavior. The model showed that this increase in the macroscopic isotropy was due to a randomization of the initial texture and also to an increase of the crystal behavior isotropy (similar values of the CRSSs of the different modes). Finally, the model was used to analyze the effect of temperature on the crystal behaviour of a Mg-Mn-Nd alloy. The introduction in the model of non-Schmid effects on the pyramidal hc+ai slip allowed to capture the inverse strength differential that appeared, between the tension and compression, above 150_C. This is the first time, to the author's knowledge, that non-Schmid effects have been reported for Mg alloys.

Fracture micromechanisms and mechanical behavior of YBCO bulk superconductors at 77 and 300 K

In this study two YBa2Cu3O7−δ bulk superconductors were evaluated, with the aim of analyzing the influence of the processing method (TSMG and Bridgman) and the test temperature on their mechanical behavior. The relationship between their mechanical properties and fracture micromechanisms has also been studied. Both materials were tested at room and at service temperature. TPB tests were carried out to determine their mechanical behavior, strength and toughness. Moreover, one of the two materials, characterized by transversal microstructural anisotropy, was tested in two directions. Hardness of both materials at nano and micro scale was studied. The results show that the mechanical behavior of the materials is controlled by the defects and cracks that have been introduced during the processing of the materials. A good degree of agreement was found between the experimental crack defects detected by means of SEM and those gathered from the fracture mechanical analysis of the experimental data

Mechanical characterization of GdBCO/Ag and YBCO single grains fabricated by top-seeded melt growth at 77 and 300 K

YBaCuO and GdBaCuO + 15 wt% Ag large, single-grain, bulk superconductors have been fabricated via the top-seeded, melt-growth (TSMG) process using a generic NdBCO seed. The mechanical behavior of both materials has been investigated by means of three-point bending (TPB) and transversal tensile tests at 77 and 300 K. The strength, fracture toughness and hardness of the samples were studied for two directions of applied load to obtain comprehensive information about the effect of microstructural anisotropy on the macroscopic and microscopic mechanical properties of these technologically important materials. Splitting (Brazilian) tests were carried out on as-melt-processed cylindrical samples following a standard oxygenation process and with the load applied parallel to the growth-facet lines characteristic of the TSMG process. In addition, the elastic modulus of each material was measured by three different techniques and related to the microstructure of each sample using optical microscopy. The results show that both the mechanical properties and the elastic modulus of both YBCO and GdBCP/Ag are improved at 77 K. However, the GdBCO/Ag samples are less anisotropic and exhibit better mechanical behavior due to the presence of silver particles in the bulk, superconducting matrix. The splitting tensile strength was determined at 77 K and both materials were found to exhibit similar behavior, independently of their differences in microstructure.

Improvement of the mechanical properties of 2124 Al/SiCp MMC plate by optimisation of the solution treatment

Static mechanical properties of 2124 Al/SiCp MMC have been measured as a function of solution temperature and time. An optimum solution treatment has been established which produces significant improvements in static mechanical properties and fatigue crack growth resistance over conventional solution treatments. Increasing the solution treatment parameters up to the optimum values improves the mechanical properties because of intermetallic dissolution, improved solute and GPB zone strengthening and increased matrix dislocation density. Increasing the solution treatment parameters beyond the optimum values results in a rapid reduction in mechanical properties due to the formation of gas porosity and surface blisters. The optimum solution treatment improves tensile properties in the transverse orientation to a greater extent than in the longitudinal orientation and this results in reduced anisotropy. © 1996 Elsevier Science Limited.

Effect of anisotropy on elastic moduli measured by nanoindentation in human tibial cortical bone

Many biological materials are known to be anisotropic. In particular, microstructural components of biological materials may grow in a preferred direction, giving rise to anisotropy in the microstructure. Nanoindentation has been shown to be an effective technique for determining the mechanical properties of microstructures as small as a few microns. However, the effects of anisotropy on the properties measured by nanoindentation have not been fully addressed. This study presents a method to account for the effects of anisotropy on elastic properties measured by nanoindentation. This method is used to correlate elastic properties determined from earlier nanoindentation experiments and from earlier ultrasonic velocity measurements in human tibial cortical bone. Also presented is a procedure to determine anisotropic elastic moduli from indentation measurements in multiple directions. © 2001 John Wiley & Sons, Inc. J Biomed Mater Res.

The impact of mechanical compression on cortical microtubules in Arabidopsis: a quantitative pipeline

Exogenous mechanical perturbations on living tissues are commonly used to investigate whether cell effectors can respond to mechanical cues. However, in most of these experiments, the applied mechanical stress and/or the biological response are described only qualitatively. We developed a quantitative pipeline based on microindentation and image analysis to investigate the impact of a controlled and prolonged compression on microtubule behaviour in the Arabidopsis shoot apical meristem, using microtubule fluorescent marker lines. We found that a compressive stress, in the order of magnitude of turgor pressure, induced apparent microtubule bundling. Importantly, that response could be reversed several hours after the release of compression. Next, we tested the contribution of microtubule severing to compression-induced bundling: microtubule bundling seemed less pronounced in the katanin mutant, in which microtubule severing is dramatically reduced. Conversely, some microtubule bundles could still be observed 16 hours after the release of compression in the spiral2 mutant, in which severing rate is instead increased. To quantify the impact of mechanical stress on anisotropy and orientation of microtubule arrays, we used the nematic tensor based FibrilTool ImageJ/Fiji plugin. To assess the degree of apparent bundling of the network, we developed several methods, some of which were borrowed from geostatistics. The final microtubule bundling response could notably be related to tissue growth velocity that was recorded by the indenter during compression. Because both input and output are quantified, this pipeline is an initial step towards correlating more precisely the cytoskeleton response to mechanical stress in living tissues.

The mechanical strength of phosphates under friction-induced cross-linking

Purpose: In the present study, we consider mechanical properties of phosphate glasses under high temperatureinduced and under friction-induced cross-linking, which enhance the modulus of elasticity. Design/methodology/approach: Two nanomechanical properties are evaluated, the first parameter is the modulus of elasticity (E) (or Young's modulus) and the second parameter is the hardness (H). Zinc meta-, pyro - and orthophosphates were recognized as amorphous-colloidal nanoparticles were synthesized under laboratory conditions and showed antiwear properties in engine oil. Findings: Young's modulus of the phosphate glasses formed under high temperature was in the 60-89 GPa range. For phosphate tribofilm formed under friction hardness and the Young's modulus were in the range of 2-10 GPa and 40-215 GPa, respectively. The degree of cross-linking during friction is provided by internal pressure of about 600 MPa and temperature close to 1000°C enhancing mechanical properties by factor of 3 (see Fig 1). Research limitations/implications: The addition of iron or aluminum ions to phosphate glasses under high temperature - and friction-induced amorphization of zinc metaphosphate and pyrophosphate tends to provide more cross-linking and mechanically stronger structures. Iron and aluminum (FeO4 or AlO4 units), incorporated into phosphate structure as network formers, contribute to the anion network bonding by converting the P=O bonds into bridging oxygen. Future work should consider on development of new of materials prepared by solgel processes, eg., zinc (II)-silicic acid. Originality/value: This paper analyses the friction pressure-induced and temperature–induced the two factors lead phosphate tribofilm glasses to chemically advanced glass structures, which may enhance the wear inhibition. Adding the coordinating ions alters the pressure at which cross-linking occurs and increases the antiwear properties of the surface material significantly.