Crack propagation in graphene

The crack initiation and growth mechanisms in an 2D graphene lattice structure are studied based on molecular dynamics simulations. Crack growth in an initial edge crack model in the arm-chair and the zig-zag lattice configurations of graphene are considered. Influence of the time steps on the post yielding behaviour of graphene is studied. Based on the results, a time step of 0.1 fs is recommended for consistent and accurate simulation of crack propagation. Effect of temperature on the crack propagation in graphene is also studied, considering adiabatic and isothermal conditions. Total energy and stress fields are analyzed. A systematic study of the bond stretching and bond reorientation phenomena is performed, which shows that the crack propagates after significant bond elongation and rotation in graphene. Variation of the crack speed with the change in crack length is estimated. (C) 2015 AIP Publishing LLC.

Statistical characterization of soft tissue inelastic parameters and application to the material failure of the aortic ach

Caracterización de los procesos de disipación mecánica basándose en la microestructura de los tejidos blandos. We present a continuous damage model with regularized softening (smeared crack models) for fiber reinforced soft tissues. Material parameters of the continuous model derive from the mesoscopic scale. In the mesoscopic scale continuum is considered as a collagenous fibrilreinforced composite. We want to study the continnumlevel response as a function of the nanoscale properties of the collagen and the adherent forces between the tropocollagen molecules.

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.

Prediction of load-deflection behaviour using fracture energy GF: A comparative study of two fracture models for concrete

Load-deflection curves for a notched beam under three-point load are determined using the Fictitious Crack Model (FCM) and Blunt Crack Model (BCM). Two values of fracture energy GF are used in this analysis: (i) GF obtained from the size effect law and (ii) GF obtained independently of the size effect. The predicted load-deflection diagrams are compared with the experimental ones obtained for the beams tested by Jenq and Shah. In addition, the values of maximum load (Pmax) obtained by the analyses are compared with the experimental ones for beams tested by Jenq and Shah and by Bažant and Pfeiffer. The results indicate that the descending portion of the load-deflection curve is very sensitive to the GF value used.

Composite material and piezoelectric coefficient uncertainty effects on structural health monitoring using feedback control gains as damage indicators

This article addresses uncertainty effect on the health monitoring of a smart structure using control gain shifts as damage indicators. A finite element model of the smart composite plate with surface-bonded piezoelectric sensors and actuators is formulated using first-order shear deformation theory and a matrix crack model is integrated into the finite element model. A constant gain velocity/position feedback control algorithm is used to provide active damping to the structure. Numerical results show that the response of the structure is changed due to matrix cracks and this change can be compensated by actively tuning the feedback controller. This change in control gain can be used as a damage indicator for structural health monitoring. Monte Carlo simulation is conducted to study the effect of material uncertainty on the damage indicator by considering composite material properties and piezoelectric coefficients as independent random variables. It is found that the change in position feedback control gain is a robust damage indicator.

Modeling of concrete cracking-A hybrid technique of using displacement discontinuity element method and direct boundary element method

This paper presents a new approach by making use of a hybrid method of using the displacement discontinuity element method and direct boundary element method to model concrete cracking by incorporating fictitious crack model. Fracture mechanics approach is followed using the Hillerborg's fictitious crack model. A boundary element based substructure method and a hybrid technique of using displacement discontinuity element method and direct boundary element method are compared in this paper. In order to represent the process zone ahead of the crack, closing forces are assumed to act in such a way that they obey a linear normal stress-crack opening displacement law. Plain concrete beams with and without initial crack under three-point loading were analyzed by both the methods. The numerical results obtained were shown to agree well with the results from existing finite element method. The model is capable of reproducing the whole range of load-deflection response including strain-softening and snap-back behavior as illustrated in the numerical examples. (C) 2011 Elsevier Ltd. All rights reserved.

Bilinear tension softening diagrams of concrete mixes corresponding to their size-independent specific fracture energy

The study of the fracture behaviour of concrete structures using the fictitious crack model requires two fracture properties of the concrete mix, namely, the size-independent specific fracture energy G(F). and the corresponding tension softening relation sigma(w) between the residual stress carrying capacity sigma and the crack opening w in the fracture process zone ahead of a real crack. In this paper, bi-linear tension softening diagrams of three different concrete mixes, ranging in compressive strength from 57 to 122 MPa whose size-independent specific fracture energy has been previously determined, have been constructed in an inverse manner based on the concept of a non-linear hinge from the load-crack mouth opening plots of notched three-point bend beams. (C) 2013 Elsevier Ltd. All rights reserved.

Mechanisms of thermal shock failure for ultra-high temperature ceramic

The materials considered in our analysis were ZrB2 ceramic matrix composites. Effect of two different additives (graphite and AlN) on thermal shock stability for the materials was measured by water quench test. It showed that it may provide more stable thermal shock properties with additives of graphite. It was explained by different thermal properties and crack resistance of the two materials in detail. Surface oxidation was one of main reasons for strength degradation of ceramic with additives of graphite after quenched in water, and surface crack was one of main reasons for strength degradation of ceramic with additives of AlN after quenched in water. It was presented that it was a potential method for improving thermal shock stability of ZrB2 ceramic matrix composites by introducing proper quantities of graphite.

Caracterização geológico-geofísica do meio aqüífero fissural : uma contribuição aos modelos de fluxo e armazenamento de água subterrânea

Northeastern Brazil is mainly formed by crystalline terrains (around 60% in area). Moreover, this region presents a semi-arid climate so that it is periodically subject to drought seasons. Furthermore, ground water quality extracted fromwells usually presents poor quality because of their high salinity contents. Nevertheless, ground water is still a very important source of water for human and animal consumption in this region. Well sitting in hard rocks terrains in Northeastern Brazil offers a mean success index of aboul 60%, given that a successful siting is defined by a well producing at least 0.5 m³/h. This low index reveals lack of knowledga about the true conditions of storage and percolation of ground water in crystalline rocks. Two models for structures storing and producing ground water in crystalline rocks in Northeastem Brazil have been proposed in the literature. The first model,tradnionally used for well sitting since the sixties are controlled by faults or fractures zones. This model is commonly referred, in Brazilian hydrogeological literature, as the "creek-crack" model (riacho-fenda in Portuguese). Sites appearing to present dense drainage network are preferred for water well siting - particularly at points where the drainages cross-cul each other. Field follow up work is usually based only on geological criteria. The second model is the "eluvio-alluvial through" (calha eluvio-aluvionar in Portuguese); it is also described in the literature but it is not yet incorporated in well sitting practice. This model is based on the hypothesis that reclilinear drainages can also be controlled by the folietion of the rock. Eventually, depending upon the degree of weathering, a through-shaped structure filled with sediments (alluvium and regolith) can be developed which can store and water can be produced from. Using severalfield case studies, this Thesis presents a thorough analysis ofthe two above cited models and proposes a new model. The analysis is based on an integrated methodological approach using geophysics and structural geology. Both land (Resitiviy and Ground Penetrating Radar- GPR) and aerogeophysical (magnetics and frequency domain eletromagnetics) surveys were used. Slructural analysis emphasized neolectonic aspects; in general, itwas found that fractures in the E-W direction are relatively open, as compared to fracturas inthe N-S direction, probably because E-W fractures were opened by the neotectonic stress regime in Northeastern Brazil, which is controlled by E-W compression and N-S extension. The riacho-fenda model is valid where drainages are controlled by fractures. The degree of fracturing and associated weathering dictale the hydrogeological potential of the structure. Field work in structural analogues reveals that subvertical fractures show consistent directions both in outcrop and aerophotograph scales. Geophysical surveys reveal subvertical conductive anomalies associated to the fracture network controlling the drainage; one of the borders of the conductive anomaly usually coincide wih the drainage. An aspect of particular importance to the validation of fracture control are the possible presence of relalively deep conductive anomalies wihoul continuation or propagalion to the surface. The conductive nature of lhe anomaly is due to the presence of wealhered rock and sedirnenls (alluvium and/or regolilh) storing ground waler which occur associated to the fracture network. Magnetic surveys are not very sensisnive to these structures.lf soil or covering sedirnents are resislive (> 100 Ohm.m), GPR can ba used to image precisely lhe fracture network. A major limialion of riacho-fenda model, revealed by GPR images, is associated to the fact thal subhorizontal fractures do play a very important role in connecting the fracture network, besides connect shallow recharge zones to relalively deep subvertical frecture zones. Iffractures play just a secondary control on the drainage, however, r/acho-fenda model may have a very limiled validny; in these cases, large portions oflhe drainage do nol coincide wilh frectures and mosl oflhewells localed in lhe drainage surrounding would resull dry. Usually, a secondary conlrol on lhe drainage by Ihefraclure networkcan be revealed only wilh detailed geophysical survey. The calha elClv1o-aluvlonarmodel is valid where drainages are conlrolled by folialion. The degree 01 wealhering 01 lhe lolialion planes dictales lhe hydrogeological polenlial 01 lhe slruclure. Outcrop analysis reveals Ihal lolialion and drainage direclions are parallel and Ihal no Iraclures, orfraclures wilh diflerent directions 01 lhe drainage direclion occur. Geophysical surveys reveal conduclive anomalies in a slab lorm associaled 10 lhe Ihrough 01 lhe wealhered rock and sedimenls (alluvium and/or regolith). Magnelic surveys can ofler a very good conlrol on lolialion direclion. An importanl aspect 10 validale lolialion conlrol are lhe presence 01 conductive anomalies showing shallow and deep portions area which are linked. Illhere is an exlensive soil cover, r/acho-fenda and calha eIClv1o-aluv/onar conlrols can be easily misinlerpreled in lhe absence 01 geophysical conlrol. Certainly, Ihis lacl could explain at leasl a part of lhe failure index in well sitting. The model wealhering sack (bolsllo de Intempertsmo in Portuguese) is proposed to explain cases where a very inlensive wealhering occur over lhe crystalline rock so Ihal a secondary inlerslilial porosity is crealed. The waler is Ihen stored in lhe porous of lhe regolilh in a similar mannerlo sedimentary rocks. A possible example ofthis model was delecled by using land geophysical survey where a relalivelyvery deep isolaled conduclive anomaly, in a slab form, was delected. Iflhis structure does store ground waler, certainly Ihere must be a link 01 lhe deep slructure wilh lhe surface in orderlo provide walerfeeding. This model mighl explain anomalous waler yields as greal as 50 m³/h Ihalsomelimescan occur in crystalline rocks in Northeaslern Brazil

Dependência estrutural da capacidade de rotação plástic em vigas de concreto armado

Summary In this work the structural dependence of plastic rotation capacity in RC beams is evaluated using the Finite Element Method. The objective is to achieve a better understanding of the non-linear behavior of reinforced concrete members and perform extensive parameter studies, using a rational model developed by Bigaj [1] to analyze the phenomenon of plastic rotation capacity in reinforced concrete members. It is assumed that only bending failure is relevant due to sufficient member resistance against shear and torsion. The paper begins with the physical and theoretical background of the phenomenon of plastic hinge development in RC structures. Special emphasis is laid on the issue of structural dependence of deformation capacity of plastic hinges in RC members. Member size dependence and influence of properties of construction materials were emphasized as well. The essential components of the Bigajs model for calculating the plastic rotation capacity are discussed. The behaviour of the plastic hinge is analysed taking into account the strain localisation in the damage zones of the hinge region. The Fictitious Crack Model (FCM) and the Compressive Damage Zone Model (CDZ) are adopted in a Fracture Mechanics approach to model the behaviour of concrete in tension and compression, respectively. The approach is implemented in FEMOOP, a FEM in-house solver under development, and applied to evaluate ductility in 2D beams. The models were generated with GiD, a pre-processor and post-processor developed by CIMNE, and analyzed with the capabilities implemented in FEMOOP. © Universitat Politècnica de Catalunya, Barcelona, España 2010.

Mixed Mode Fracture Behaviour of Piezoelectric Materials

Piezoelectrics present an interactive electromechanical behaviour that, especially in recent years, has generated much interest since it renders these materials adapt for use in a variety of electronic and industrial applications like sensors, actuators, transducers, smart structures. Both mechanical and electric loads are generally applied on these devices and can cause high concentrations of stress, particularly in proximity of defects or inhomogeneities, such as flaws, cavities or included particles. A thorough understanding of their fracture behaviour is crucial in order to improve their performances and avoid unexpected failures. Therefore, a considerable number of research works have addressed this topic in the last decades. Most of the theoretical studies on this subject find their analytical background in the complex variable formulation of plane anisotropic elasticity. This theoretical approach bases its main origins in the pioneering works of Muskelishvili and Lekhnitskii who obtained the solution of the elastic problem in terms of independent analytic functions of complex variables. In the present work, the expressions of stresses and elastic and electric displacements are obtained as functions of complex potentials through an analytical formulation which is the application to the piezoelectric static case of an approach introduced for orthotropic materials to solve elastodynamics problems. This method can be considered an alternative to other formalisms currently used, like the Stroh’s formalism. The equilibrium equations are reduced to a first order system involving a six-dimensional vector field. After that, a similarity transformation is induced to reach three independent Cauchy-Riemann systems, so justifying the introduction of the complex variable notation. Closed form expressions of near tip stress and displacement fields are therefore obtained. In the theoretical study of cracked piezoelectric bodies, the issue of assigning consistent electric boundary conditions on the crack faces is of central importance and has been addressed by many researchers. Three different boundary conditions are commonly accepted in literature: the permeable, the impermeable and the semipermeable (“exact”) crack model. This thesis takes into considerations all the three models, comparing the results obtained and analysing the effects of the boundary condition choice on the solution. The influence of load biaxiality and of the application of a remote electric field has been studied, pointing out that both can affect to a various extent the stress fields and the angle of initial crack extension, especially when non-singular terms are retained in the expressions of the electro-elastic solution. Furthermore, two different fracture criteria are applied to the piezoelectric case, and their outcomes are compared and discussed. The work is organized as follows: Chapter 1 briefly introduces the fundamental concepts of Fracture Mechanics. Chapter 2 describes plane elasticity formalisms for an anisotropic continuum (Eshelby-Read-Shockley and Stroh) and introduces for the simplified orthotropic case the alternative formalism we want to propose. Chapter 3 outlines the Linear Theory of Piezoelectricity, its basic relations and electro-elastic equations. Chapter 4 introduces the proposed method for obtaining the expressions of stresses and elastic and electric displacements, given as functions of complex potentials. The solution is obtained in close form and non-singular terms are retained as well. Chapter 5 presents several numerical applications aimed at estimating the effect of load biaxiality, electric field, considered permittivity of the crack. Through the application of fracture criteria the influence of the above listed conditions on the response of the system and in particular on the direction of crack branching is thoroughly discussed.

Numerical analysis of mixe-mode fracture of brickwork masonry with embedded discontinuity elements

One of the common pathologies of brickwork masonry structural elements and walls is the cracking associated with the differential settlements and/or excessive deflections of the slabs along the life of the structure. The scarce capacity of the masonry in order to accompany the structural elements that surround it, such as floors, beams or foundations, in their movements makes the brickwork masonry to be an element that frequently presents this kind of problem. This problem is a fracture problem, where the wall is cracked under mixed mode fracture: tensile and shear stresses combination, under static loading. Consequently, it is necessary to advance in the simulation and prediction of brickwork masonry mechanical behaviour under tensile and shear loading. The quasi-brittle behaviour of the brickwork masonry can be studied using the cohesive crack model whose application to other quasibrittle materials like concrete has traditionally provided very satisfactory results.

Experimental and numerical analysis of reinforced concrete elements subjected to blast loading

El hormigón es uno de los materiales de construcción más empleados en la actualidad debido a sus buenas prestaciones mecánicas, moldeabilidad y economía de obtención, entre otras ventajas. Es bien sabido que tiene una buena resistencia a compresión y una baja resistencia a tracción, por lo que se arma con barras de acero para formar el hormigón armado, material que se ha convertido por méritos propios en la solución constructiva más importante de nuestra época. A pesar de ser un material profusamente utilizado, hay aspectos del comportamiento del hormigón que todavía no son completamente conocidos, como es el caso de su respuesta ante los efectos de una explosión. Este es un campo de especial relevancia, debido a que los eventos, tanto intencionados como accidentales, en los que una estructura se ve sometida a una explosión son, por desgracia, relativamente frecuentes. La solicitación de una estructura ante una explosión se produce por el impacto sobre la misma de la onda de presión generada en la detonación. La aplicación de esta carga sobre la estructura es muy rápida y de muy corta duración. Este tipo de acciones se denominan cargas impulsivas, y pueden ser hasta cuatro órdenes de magnitud más rápidas que las cargas dinámicas impuestas por un terremoto. En consecuencia, no es de extrañar que sus efectos sobre las estructuras y sus materiales sean muy distintos que las que producen las cargas habitualmente consideradas en ingeniería. En la presente tesis doctoral se profundiza en el conocimiento del comportamiento material del hormigón sometido a explosiones. Para ello, es crucial contar con resultados experimentales de estructuras de hormigón sometidas a explosiones. Este tipo de resultados es difícil de encontrar en la literatura científica, ya que estos ensayos han sido tradicionalmente llevados a cabo en el ámbito militar y los resultados obtenidos no son de dominio público. Por otra parte, en las campañas experimentales con explosiones llevadas a cabo por instituciones civiles el elevado coste de acceso a explosivos y a campos de prueba adecuados no permite la realización de ensayos con un elevado número de muestras. Por este motivo, la dispersión experimental no es habitualmente controlada. Sin embargo, en elementos de hormigón armado sometidos a explosiones, la dispersión experimental es muy acusada, en primer lugar, por la propia heterogeneidad del hormigón, y en segundo, por la dificultad inherente a la realización de ensayos con explosiones, por motivos tales como dificultades en las condiciones de contorno, variabilidad del explosivo, o incluso cambios en las condiciones atmosféricas. Para paliar estos inconvenientes, en esta tesis doctoral se ha diseñado un novedoso dispositivo que permite ensayar hasta cuatro losas de hormigón bajo la misma detonación, lo que además de proporcionar un número de muestras estadísticamente representativo, supone un importante ahorro de costes. Con este dispositivo se han ensayado 28 losas de hormigón, tanto armadas como en masa, de dos dosificaciones distintas. Pero además de contar con datos experimentales, también es importante disponer de herramientas de cálculo para el análisis y diseño de estructuras sometidas a explosiones. Aunque existen diversos métodos analíticos, hoy por hoy las técnicas de simulación numérica suponen la alternativa más avanzada y versátil para el cálculo de elementos estructurales sometidos a cargas impulsivas. Sin embargo, para obtener resultados fiables es crucial contar con modelos constitutivos de material que tengan en cuenta los parámetros que gobiernan el comportamiento para el caso de carga en estudio. En este sentido, cabe destacar que la mayoría de los modelos constitutivos desarrollados para el hormigón a altas velocidades de deformación proceden del ámbito balístico, donde dominan las grandes tensiones de compresión en el entorno local de la zona afectada por el impacto. En el caso de los elementos de hormigón sometidos a explosiones, las tensiones de compresión son mucho más moderadas, siendo las tensiones de tracción generalmente las causantes de la rotura del material. En esta tesis doctoral se analiza la validez de algunos de los modelos disponibles, confirmando que los parámetros que gobiernan el fallo de las losas de hormigón armado ante explosiones son la resistencia a tracción y su ablandamiento tras rotura. En base a los resultados anteriores se ha desarrollado un modelo constitutivo para el hormigón ante altas velocidades de deformación, que sólo tiene en cuenta la rotura por tracción. Este modelo parte del de fisura cohesiva embebida con discontinuidad fuerte, desarrollado por Planas y Sancho, que ha demostrado su capacidad en la predicción de la rotura a tracción de elementos de hormigón en masa. El modelo ha sido modificado para su implementación en el programa comercial de integración explícita LS-DYNA, utilizando elementos finitos hexaédricos e incorporando la dependencia de la velocidad de deformación para permitir su utilización en el ámbito dinámico. El modelo es estrictamente local y no requiere de remallado ni conocer previamente la trayectoria de la fisura. Este modelo constitutivo ha sido utilizado para simular dos campañas experimentales, probando la hipótesis de que el fallo de elementos de hormigón ante explosiones está gobernado por el comportamiento a tracción, siendo de especial relevancia el ablandamiento del hormigón. Concrete is nowadays one of the most widely used building materials because of its good mechanical properties, moldability and production economy, among other advantages. As it is known, it has high compressive and low tensile strengths and for this reason it is reinforced with steel bars to form reinforced concrete, a material that has become the most important constructive solution of our time. Despite being such a widely used material, there are some aspects of concrete performance that are not yet fully understood, as it is the case of its response to the effects of an explosion. This is a topic of particular relevance because the events, both intentional and accidental, in which a structure is subjected to an explosion are, unfortunately, relatively common. The loading of a structure due to an explosive event occurs due to the impact of the pressure shock wave generated in the detonation. The application of this load on the structure is very fast and of very short duration. Such actions are called impulsive loads, and can be up to four orders of magnitude faster than the dynamic loads imposed by an earthquake. Consequently, it is not surprising that their effects on structures and materials are very different than those that cause the loads usually considered in engineering. This thesis broadens the knowledge about the material behavior of concrete subjected to explosions. To that end, it is crucial to have experimental results of concrete structures subjected to explosions. These types of results are difficult to find in the scientific literature, as these tests have traditionally been carried out by armies of different countries and the results obtained are classified. Moreover, in experimental campaigns with explosives conducted by civil institutions the high cost of accessing explosives and the lack of proper test fields does not allow for the testing of a large number of samples. For this reason, the experimental scatter is usually not controlled. However, in reinforced concrete elements subjected to explosions the experimental dispersion is very pronounced. First, due to the heterogeneity of concrete, and secondly, because of the difficulty inherent to testing with explosions, for reasons such as difficulties in the boundary conditions, variability of the explosive, or even atmospheric changes. To overcome these drawbacks, in this thesis we have designed a novel device that allows for testing up to four concrete slabs under the same detonation, which apart from providing a statistically representative number of samples, represents a significant saving in costs. A number of 28 slabs were tested using this device. The slabs were both reinforced and plain concrete, and two different concrete mixes were used. Besides having experimental data, it is also important to have computational tools for the analysis and design of structures subjected to explosions. Despite the existence of several analytical methods, numerical simulation techniques nowadays represent the most advanced and versatile alternative for the assessment of structural elements subjected to impulsive loading. However, to obtain reliable results it is crucial to have material constitutive models that take into account the parameters that govern the behavior for the load case under study. In this regard it is noteworthy that most of the developed constitutive models for concrete at high strain rates arise from the ballistic field, dominated by large compressive stresses in the local environment of the area affected by the impact. In the case of concrete elements subjected to an explosion, the compressive stresses are much more moderate, while tensile stresses usually cause material failure. This thesis discusses the validity of some of the available models, confirming that the parameters governing the failure of reinforced concrete slabs subjected to blast are the tensile strength and softening behaviour after failure. Based on these results we have developed a constitutive model for concrete at high strain rates, which only takes into account the ultimate tensile strength. This model is based on the embedded Cohesive Crack Model with Strong Discontinuity Approach developed by Planas and Sancho, which has proved its ability in predicting the tensile fracture of plain concrete elements. The model has been modified for its implementation in the commercial explicit integration program LS-DYNA, using hexahedral finite elements and incorporating the dependence of the strain rate, to allow for its use in dynamic domain. The model is strictly local and does not require remeshing nor prior knowledge of the crack path. This constitutive model has been used to simulate two experimental campaigns, confirming the hypothesis that the failure of concrete elements subjected to explosions is governed by their tensile response, being of particular relevance the softening behavior of concrete.

Estudio numérico y experimental del efecto de tamaño en modo mixto de los paneles de yeso laminado y lana de roca

Los paneles sándwich de yeso laminado y lana de roca presentan una abundante patología de fisuración debida a flechas excesivas de forjados. Existe, por tanto, la necesidad de avanzar en la simulación y predicción de comportamiento bajo solicitaciones de tracción y cortante de ese tipo de paneles, a pesar de que en las aplicaciones habituales no tienen responsabilidad estructural. El comportamiento de este material puede ser considerado cuasi-frágil, y en base a ello en este trabajo ha sido estudiado haciendo uso de modelos de fisura cohesiva, cuya aplicación a otros materiales cuasifrágiles, como el hormigón, ha aportado resultados muy satisfactorios. En esta comunicación se presenta el trabajo realizado para estudiar el efecto del tamaño del elemento de yeso laminado y lana de roca en su comportamiento mecánico-resistente. Para ello se diseñó una campaña de ensayos en modo mixto sobre probetas de diferente tamaño. Se han realizado ensayos de flexión en tres puntos en modo mixto de unas probetas entalladas, geométricamente similares y de diferente tamaño, obteniéndose las curvas carga-desplazamiento y cargaabertura de la boca de la entalla. Para simular numéricamente el comportamiento en fractura del panel en modo mixto se ha utilizado un modelo de elementos finitos con fisura embebida basado en la fisura cohesiva en el que se introducen como entrada los parámetros obtenidos a partir de la experimentación de trabajos anteriores, obteniéndose un buen ajuste. En función de estos resultados se analiza el efecto del tamaño en los paneles. Sandwich panels of laminated gypsum and rockwool have an abundant pathology of cracking due to excessive slabs deflection. Therefore, it is necessary to progress in the simulation and prediction of behaviour under tensile and shear load of such panels, although in typical applications have no structural responsability. The behaviour of this material may be considered quasi-brittle and, based on this idea, in this work has been studied using a cohesive crack model that has been applied to other quasi-brittle materials, such as concrete, and has provided very satisfactory results. This communication presents the work carried out to study the size effect of the specimen of plasterboard and rockwool in its mechanical and resistant behaviour. The authors designed an experimental campaign under mixed mode composed by testing specimens of different sizes. Assymetrical three-point bending tests have been performed on notched specimens, geometrically similar and of different size, to obtain load-displacement and load-crack moutn opening displacement curves. To numerically simulate the mixed-mode fracture behaviour of the panels we have used a finite element model with embedded crack, based on the cohesive crack model, using as input the experimental parameters obtained in previous work, obtaining a good adjustment. Based on these results we analyze the size effect of the panels

Estudio numérico y experimental del efecto de tamaño en modo mixto de los paneles de yeso laminado y lana de roca

Los paneles sándwich de yeso laminado y lana de roca presentan una abundante patología de fisuración debida a flechas excesivas de forjados. Existe, por tanto, la necesidad de avanzar en la simulación y predicción de comportamiento bajo solicitaciones de tracción y cortante de ese tipo de paneles, a pesar de que en las aplicaciones habituales no tienen responsabilidad estructural. El comportamiento de este material puede ser considerado cuasi-frágil, y en base a ello en este trabajo ha sido estudiado haciendo uso de modelos de fisura cohesiva, cuya aplicación a otros materiales cuasifrágiles, como el hormigón, ha aportado resultados muy satisfactorios. En esta comunicación se presenta el trabajo realizado para estudiar el efecto del tamaño del elemento de yeso laminado y lana de roca en su comportamiento mecánico-resistente. Para ello se diseñó una campaña de ensayos en modo mixto sobre probetas de diferente tamaño. Se han realizado ensayos de flexión en tres puntos en modo mixto de unas probetas entalladas, geométricamente similares y de diferente tamaño, obteniéndose las curvas carga-desplazamiento y cargaabertura de la boca de la entalla. Para simular numéricamente el comportamiento en fractura del panel en modo mixto se ha utilizado un modelo de elementos finitos con fisura embebida basado en la fisura cohesiva en el que se introducen como entrada los parámetros obtenidos a partir de la experimentación de trabajos anteriores, obteniéndose un buen ajuste. En función de estos resultados se analiza el efecto del tamaño en los paneles. Sandwich panels of laminated gypsum and rockwool have an abundant pathology of cracking due to excessive slabs deflection. Therefore, it is necessary to progress in the simulation and prediction of behaviour under tensile and shear load of such panels, although in typical applications have no structural responsability. The behaviour of this material may be considered quasi-brittle and, based on this idea, in this work has been studied using a cohesive crack model that has been applied to other quasi-brittle materials, such as concrete, and has provided very satisfactory results. This communication presents the work carried out to study the size effect of the specimen of plasterboard and rockwool in its mechanical and resistant behaviour. The authors designed an experimental campaign under mixed mode composed by testing specimens of different sizes. Assymetrical three-point bending tests have been performed on notched specimens, geometrically similar and of different size, to obtain load-displacement and load-crack moutn opening displacement curves. To numerically simulate the mixed-mode fracture behaviour of the panels we have used a finite element model with embedded crack, based on the cohesive crack model, using as input the experimental parameters obtained in previous work, obtaining a good adjustment. Based on these results we analyze the size effect of the panels.