Automatic Insertion of Cohesive Elements for Delamination

Composite damage modelling with cohesive elements has initially been limited to the analysis of interface damage or delamination. However, their use is also being extended to the analysis of inplane tensile failure arising from matrix or fibre fracture. These interface elements are typically placed at locations where failure is likely to occur, which infers a certain a priori knowledge of the crack propagation path(s). In the case of a crack jump for example, the location of the jump is usually not obvious, and the simulation would require the placement of cohesive elements at all element faces. A better option, presented here, is to determine the potential location of cohesive elements and insert them during the analysis. The aim of this work is to enable the determination of the crack path, as part of the solution process. A subroutine has been developed and implemented in the commercial finite element package ABAQUS/Standard[1] in order to automatically insert cohesive elements within a pristine model, on the basis of the analysis of the current stress field. Results for the prediction of delamination are presented in this paper.

Improving composite damage modelling through automatic placement of cohesive elements

Numerous experimental studies of damage in composite laminates have shown that intralaminar (in-plane) matrix cracks lead to interlaminar delamination (out-of-plane) at ply interfaces. The smearing of in-plane cracks over a volume, as a consequence of the use of continuum damage mechanics, does not always effectively capture the full extent of the interaction between the two failure mechanisms. A more accurate representation is obtained by adopting a discrete crack approach via the use of cohesive elements, for both in-plane and out-of-plane damage. The difficulty with cohesive elements is that their location must be determined a priori in order to generate the model; while ideally the position of the crack migration, and more generally the propagation path, should be obtained as part of the problem’s solution. With the aim of enhancing current modelling capabilities with truly predictive capabilities, a concept of automatic insertion of interface elements is utilized. The consideration of a simple traction criterion in relation to material strength, evaluated at each node of the model (or of the regions of the model where it is estimated cracks might form), allows for the determination of initial crack location and subsequent propagation by the insertion of cohesive elements during the course of the analysis. Several experimental results are modelled using the commercial package ABAQUS/Standard with an automatic insertion subroutine developed in this work, and the results are presented to demonstrate the capabilities of this technique.

Monte Carlo simulation of complex cohesive fracture in random heterogeneous quasi-brittle materials

A numerical method is developed to simulate complex two-dimensional crack propagation in quasi-brittle materials considering random heterogeneous fracture properties. Potential cracks are represented by pre-inserted cohesive elements with tension and shear softening constitutive laws modelled by spatially varying Weibull random fields. Monte Carlo simulations of a concrete specimen under uni-axial tension were carried out with extensive investigation of the effects of important numerical algorithms and material properties on numerical efficiency and stability, crack propagation processes and load-carrying capacities. It was found that the homogeneous model led to incorrect crack patterns and load–displacement curves with strong mesh-dependence, whereas the heterogeneous model predicted realistic, complicated fracture processes and load-carrying capacity of little mesh-dependence. Increasing the variance of the tensile strength random fields with increased heterogeneity led to reduction in the mean peak load and increase in the standard deviation. The developed method provides a simple but effective tool for assessment of structural reliability and calculation of characteristic material strength for structural design.

A hybrid embedded cohesive element method for predicting matrix cracking in composites

The complex architecture of many fibre-reinforced composites makes the generation of finite element meshes a labour-intensive process. The embedded element method, which allows the matrix and fibre reinforcement to be meshed separately, offers a computationally efficient approach to reduce the time and cost of meshing. In this paper we present a new approach of introducing cohesive elements into the matrix domain to enable the prediction of matrix cracking using the embedded element method. To validate this approach, experiments were carried out using a modified Double Cantilever Beam with ply drops, with the results being compared with model predictions. Crack deflection was observed at the ply drop region, due to the differences in stiffness, strength and toughness at the bi-material interface. The new modelling technique yields accurate predictions of the failure process in composites, including fracture loads and crack deflection path.

CFRP strengthened CFST columns under vehicular impact

Concrete filled steel tubular (CFST) columns are increasingly used in bridge piers and high-rise buildings due to their excellent axial load bearing capacity. These columns may experience severe damage or failure due to transverse impact of vehicle collisions. In this study, numerical investigation is carried out to evaluate the effect of carbon fibre reinforced polymer (CFRP) strengthening CFST columns under vehicular impact. The CFRP composites damage mechanisms are simulated to account four different failure criteria. The cohesive elements are introduced as interface element to properly simulate the adhesively bonded regime. Simplified vehicle model is also developed to represent real vehicle behaviour. The FE analysis results show that externally bonded CFRP composites improve the impact resistance capacity compared to bare CFST column.

Att skapa en helhet av fragment Om översättningen av kohesiva element i Sjóns Stålnatt

Abstract (Creating a whole from fragments. On the translation of cohesive elements in Sjón's Steelnight): In this paper, I discuss the Swedish translation of the Icelandic writer Sjón’s debut novel Steelnight: a story (in Swedish Stålnatt – en berättelse). Steelnight is a fragmented novel which is located at the border between poetry and prose. At the beginning of the story, the writer introduces two main plots (with subplots). At first glance, it may seem as though the plots do not share much in common. However, a closer analysis of the text reveals that there are striking parallels between the plots which are created partly by the writer’s use of cohesive elements. Thus, despite the fragmented style, the writer creates cohesion in the text by using stylistic means such as repetition of certain nouns and adjectives and by a frequent use of personal and demonstrative pronouns. This use of cohesive elements creates a challenge for the translator.

Numerical analysis of stiffener runout sections

The recent trend of incorporating more composite material in primary aircraft structures has highlighted the vulnerability of stiffened aerostructures to through-thickness stresses, which may lead to delamination and debonding at the skin-stiffener interface, leading to collapse. Stiffener runout regions are particularly susceptible to this problem and cannot be avoided due to the necessity to terminate stiffeners at rib intersections or at cutouts, interrupting the stiffener load path. In this paper, experimental tests relating to two different stiffener runout specimens are presented and the failure modes of both specimens are discussed in detail. A thinner-skinned specimen showed sudden and unstable crack propagation, while a thicker-skinned specimen showed initially unstable but subsequent stable crack growth. Detailed finite element models of the two specimens are developed, and it is shown how such models can explain and predict the behaviour and failure mode of stiffener runouts. The models contain continuum shell elements to model the skin and stiffener, while cohesive elements using a traction-separation law are placed at the skin-stiffener interface to effectively model the debonding which promotes structural failure.

Buckling strength of adhesively-bonded single and double-strap repairs on carbon-epoxy structures

This work reports on an experimental and finite element method (FEM) parametric study of adhesively-bonded single and double-strap repairs on carbon-epoxy structures under buckling unrestrained compression. The influence of the overlap length and patch thickness was evaluated. This loading gains a particular significance from the additional characteristic mechanisms of structures under compression, such as fibres microbuckling, for buckling restrained structures, or global buckling of the assembly, if no transverse restriction exists. The FEM analysis is based on the use of cohesive elements including mixed-mode criteria to simulate a cohesive fracture of the adhesive layer. Trapezoidal laws in pure modes I and II were used to account for the ductility of most structural adhesives. These laws were estimated for the adhesive used from double cantilever beam (DCB) and end-notched flexure (ENF) tests, respectively, using an inverse technique. The pure mode III cohesive law was equalled to the pure mode II one. Compression failure in the laminates was predicted using a stress-based criterion. The accurate FEM predictions open a good prospect for the reduction of the extensive experimentation in the design of carbon-epoxy repairs. Design principles were also established for these repairs under buckling.

A produção textual escrita de alunos de 1ª série do ensino fundamental: sob a influência dos métodos de alfabetização

Com base em estudos que tomam como eixo central a alfabetização e o letramento, esta dissertação compreende uma investigação acerca da influência dos métodos de alfabetização, mais especificamente do método silábico e da proposta de alfabetização com base lingüística (ABL), na produção textual escrita de alunos de 1 série do Ensino Fundamental, alfabetizados ou pelo método silábico ou pela proposta ABL. Para tanto, este estudo partiu da hipótese de que a relação oralidade/escrita estabelecida pelo aluno em suas produções e a utilização de um número restrito de elementos coesivos que estabelecem as relações de sentido entre os enunciados que compõem o texto é, em parte, responsabilidade dos métodos e/ou propostas de alfabetização adotados. O diálogo estabelecido com autores que dão visibilidade à prática de alfabetizar letrando; os estudos sobre letramento, alfabetização, psicogênese da língua escrita e a produção de textos escolares ajudaram na condução desse trabalho, que além da pesquisa bibliográfica se valeu da pesquisa de campo realizada nas escolas N.S.A e L.C pertencentes à Rede Municipal de Ensino de Ananindeua e à Rede Privada de Ensino de Belém, respectivamente. Durante esta pesquisa, foram observados o desenvolvimento do trabalho de quatro professoras com turmas de alfabetização e da produção textual de doze alunos de duas turmas de 1 série da escola N.S.A; também foram aplicados dois questionários a sete professoras de turmas de alfabetização e 1 série das duas escolas citadas. A referida pesquisa tornou possível perceber que a alfabetização já começa a ser vista como um processo que visa não só ao domínio da tecnologia da leitura e da escrita, mas também aos usos sociais dessa tecnologia; mostra também que os alunos alfabetizados pela proposta ABL usam um maior número de elementos coesivos em suas produções que os alunos alfabetizados pelo método silábico e que a relação entre a fala e a escrita estabelecida na escrita de palavras é muito mais forte nas produções textuais escritas dos alunos alfabetizados pela ABL.

High fidelity simulations of failure in fiber-reinforced composites

Los ensayos virtuales de materiales compuestos han aparecido como un nuevo concepto dentro de la industria aeroespacial, y disponen de un vasto potencial para reducir los enormes costes de certificación y desarrollo asociados con las tediosas campañas experimentales, que incluyen un gran número de paneles, subcomponentes y componentes. El objetivo de los ensayos virtuales es sustituir algunos ensayos por simulaciones computacionales con alta fidelidad. Esta tesis es una contribución a la aproximación multiescala desarrollada en el Instituto IMDEA Materiales para predecir el comportamiento mecánico de un laminado de material compuesto dadas las propiedades de la lámina y la intercara. La mecánica de daño continuo (CDM) formula el daño intralaminar a nivel constitutivo de material. El modelo de daño intralaminar se combina con elementos cohesivos para representar daño interlaminar. Se desarrolló e implementó un modelo de daño continuo, y se aplicó a configuraciones simples de ensayos en laminados: impactos de baja y alta velocidad, ensayos de tracción, tests a cortadura. El análisis del método y la correlación con experimentos sugiere que los métodos son razonablemente adecuados para los test de impacto, pero insuficientes para el resto de ensayos. Para superar estas limitaciones de CDM, se ha mejorado la aproximación discreta de elementos finitos enriqueciendo la cinemática para incluir discontinuidades embebidas: el método extendido de los elementos finitos (X-FEM). Se adaptó X-FEM para un esquema explícito de integración temporal. El método es capaz de representar cualitativamente los mecanismos de fallo detallados en laminados. Sin embargo, los resultados muestran inconsistencias en la formulación que producen resultados cuantitativos erróneos. Por último, se ha revisado el método tradicional de X-FEM, y se ha desarrollado un nuevo método para superar sus limitaciones: el método cohesivo X-FEM estable. Las propiedades del nuevo método se estudiaron en detalle, y se concluyó que el método es robusto para implementación en códigos explícitos dinámicos escalables, resultando una nueva herramienta útil para la simulación de daño en composites. Virtual testing of composite materials has emerged as a new concept within the aerospace industry. It presents a very large potential to reduce the large certification costs and the long development times associated with the experimental campaigns, involving the testing of a large number of panels, sub-components and components. The aim of virtual testing is to replace some experimental tests by high-fidelity numerical simulations. This work is a contribution to the multiscale approach developed in Institute IMDEA Materials to predict the mechanical behavior of a composite laminate from the properties of the ply and the interply. Continuum Damage Mechanics (CDM) formulates intraply damage at the the material constitutive level. Intraply CDM is combined with cohesive elements to model interply damage. A CDM model was developed, implemented, and applied to simple mechanical tests of laminates: low and high velocity impact, tension of coupons, and shear deformation. The analysis of the results and the comparison with experiments indicated that the performance was reasonably good for the impact tests, but insuficient in the other cases. To overcome the limitations of CDM, the kinematics of the discrete finite element approximation was enhanced to include mesh embedded discontinuities, the eXtended Finite Element Method (X-FEM). The X-FEM was adapted to an explicit time integration scheme and was able to reproduce qualitatively the physical failure mechanisms in a composite laminate. However, the results revealed an inconsistency in the formulation that leads to erroneous quantitative results. Finally, the traditional X-FEM was reviewed, and a new method was developed to overcome its limitations, the stable cohesive X-FEM. The properties of the new method were studied in detail, and it was demonstrated that the new method was robust and can be implemented in a explicit finite element formulation, providing a new tool for damage simulation in composite materials.

Modelización tridimensional de ensayos de adherencia tipo viga del refuerzo NSM-FRP en hormigón

El uso de refuerzos NSM‐FRP en estructuras de hormigón armado se ha incrementado considerablemente en los últimos años como método de refuerzo estructural. Los ensayos de arrancamiento en viga de los refuerzos NSM‐FRP permiten el estudio del comportamiento de la unión pegada. El principal objetivo del presente trabajo aborda la simulación numérica de este tipo de ensayos, con el propósito de caracterizar correctamente la adherencia entre las barras de NSM‐FRP y el hormigón. En una fase inicial se simuló un modelo bidimensional para conseguir evaluar y verificar el comportamiento de los elementos cohesivos y ver su comportamiento primero ante diferentes modelos de material y segundo ante un modo mixto de fallo, debido a la aplicación simultanea de carga axial y carga cortante. En una segunda fase se creó un modelo tridimensional para estudiar el arrancamiento de una barra de material compuesto insertada en hormigón, creando un modelo de material de hormigón y viendo el comportamiento cualitativo del sistema ante variaciones en los parámetros de los diferentes materiales. En la tercera fase, la más importante del presente trabajo, se abordó la simulación numérica del ensayo de arrancamiento en viga. Se simularon todos los componentes del ensayo y se evaluaron diferentes alternativas para representar la interfase NSM‐FRP ‐ hormigón, usando elementos cohesivos y diferentes distribuciones de los mismos en la interfase. Para conseguir representar lo más fielmente posible las condiciones del ensayo, se diseñó también un controlador PID que permite realizar las simulaciones numéricas mediante un control en desplazamientos, lo cual permite capturar más correctamente el comportamiento de reblandecimiento de la unión pegada. El controlador PID aplica técnicas de ingeniería de control para conseguir calcular a priori la amplitud necesaria del desplazamiento impuesto que provoque una evolución establecida en una variable interna del sistema. La variable usada para correlacionar los ensayos es la diferencia en desplazamientos entre dos puntos y se escoge una evolución lineal de la misma, pero en la tesis también se exponen los resultados de escoger otras posibles variables internas con diferentes evoluciones. Se compararon las simulaciones numéricas con resultados de mediciones experimentales previamente publicadas. Los resultados carga‐deslizamiento obtenidos encajan bien con los datos experimentales. El modelo propuesto es también capaz de predecir el modo de fallo en la interfase NSM‐FRP ‐ hormigón. Finalmente, también se han llevado a cabo estudios paramétricos, para evaluar la influencia de cada parámetro en los resultados. También se realizó un estudio cualitativo de cómo se comporta la unión pegada en cada momento de la simulación, mediante el uso macros y gráficas tridimensionales, para conseguir una mejor visualización y facilitar el análisis de los resultados. ABSTRACT The use of near‐surface mounted FRP reinforcement in reinforced concrete structures has seen a considerable increase in recent years as a strengthening method. Beam pull‐out tests for near‐surface reinforcement allow obtaining the local bond‐slip behavior of a bonded joint. The main objective of the current work deals with the three‐dimensional modeling of this kind of test with the purpose of characterizing suitably the mechanics of bond between FRP rods and concrete. In an initial stage, a two bidimensional in order to evaluate and to verify the behavior of the cohesive elements. Its behavior was evaluated first testing different material models and second testing the behavior when mixed mode failure appears, due to simultaneous axial and shear load. In a second stage a tridimensional model was created in order to study the pull‐out of an inserted beam of composite material in concrete. A concrete material model was created and the influence of each material parameter was studied qualitatively. The third part, the most relevant of the present work, the numerical simulation of the Beam Pull‐Out test was faced. All the parts of the Beam Pull‐Out test were included inthe simulation and different alternatives to represent the FRP bar – concrete interface have been evaluated, using cohesive elements and different distributions of them. In order to reproduce the test conditions more reliably, a PID controller has also been designed to conduct suitably the numerical tests in order to properly capture the softening branch of the load‐slip behaviour. The PID controller applies control techniques to calculate a priori the necessary amplitude of the load in order to achieve a given evolution through the simulation of an internal variable previously chosen. The variable used in order to correlate the simulation with the test results is the difference in displacements between two points and a linear evolution was chosen, but in the thesis the results of choosing other possible internal variables with different evolutions are also shown. The numerical FE simulations were compared with experimental measurements previously published. Load‐slip predictions compare well with the corresponding experimental data. The proposed model is also able to predict the failure mode at the FRP‐concrete interface. Some parametric studies have also been carried out, in order to evaluate the influence of each material parameter in the results. A qualitative study of the behaviour of the joint was also performed, using the results of the numeric simulations and through the use of macros and 3D graphs, the tensional state of each point of the joint can be visualized in each moment of the simulation.

Vertical uplift resistance of two closely spaced horizontal strip anchors embedded in cohesive-frictional weightless medium

The vertical uplift resistance of two closely spaced horizontal strip plate anchors has been investigated by using lower and upper bound theorems of the limit analysis in combination with finite elements and linear optimization. The interference effect on uplift resistance of the two anchors is evaluated in terms of a nondimensional efficiency factor (eta(c)). The variation of eta(c) with changes in the clear spacing (S) between the two anchors has been established for different combinations of embedment ratio (H/B) and angle of internal friction of the soil (phi). An interference of the anchors leads to a continuous reduction in uplift resistance with a decrease in spacing between the anchors. The uplift resistance becomes a minimum when the two anchors are placed next to each other without any gap. The critical spacing (S-cr) between the two anchors required to eliminate the interference effect increases with an increase in the values of both H/B and phi. The value of S-cr was found to lie approximately in the range 0.65B-1.5B with H/B = 1 and 11B-14B with H/B = 7 for phi varying from 0 degrees to 30 degrees.

Solving Axisymmetric Stability Problems by Using Upper Bound Finite Elements, Limit Analysis, and Linear Optimization

A numerical formulation has been proposed for solving an axisymmetric stability problem in geomechanics with upper bound limit analysis, finite elements, and linear optimization. The Drucker-Prager yield criterion is linearized by simulating a sphere with a circumscribed truncated icosahedron. The analysis considers only the velocities and plastic multiplier rates, not the stresses, as the basic unknowns. The formulation is simple to implement, and it has been employed for finding the collapse loads of a circular footing placed over the surface of a cohesive-frictional material. The formulation can be used to solve any general axisymmetric geomechanics stability problem.

A cohesive approach to thin-shell fracture and fragmentation

We develop a finite-element method for the simulation of dynamic fracture and fragmentation of thin-shells. The shell is spatially discretized with subdivision shell elements and the fracture along the element edges is modeled with a cohesive law. In order to follow the propagation and branching of cracks, subdivision shell elements are pre-fractured ab initio and the crack opening is constrained prior to crack nucleation. This approach allows for shell fracture in an in-plane tearing mode, a shearing mode, or a bending of hinge mode. The good performance of the method is demonstrated through the simulation of petalling failure experiments in aluminum plates. © 2005 Elsevier B.V. All rights reserved.