Search results

In this work, the cohesive zone model (CZM) developed by some of the authors to simulate the propagation of fatigue defects in two dimensions is extended in order to simulate the propagation of defects in 3D. The paper aims to discuss this issue.

The procedure has been implemented in the finite element (FE) solver (Abaqus) by programming the appropriate software-embedded subroutines. Part of the procedure is devoted to the calculation of the rate of energy release per unit, G, necessary to know the growth of the defect.

The model was tested on different joint geometries, with different load conditions (pure mode I, mode II pure, mixed mode I/II) and the results of the analysis were compared with analytical solutions or virtual crack closure technique (VCCT).

The possibility to simulate the growth of a crack without any re-meshing requirements and the relatively easy possibility to manipulate the constitutive law of the cohesive elements makes the CZM attractive also for the fatigue crack growth simulation. However, differently from VCCT, three-dimensional fatigue de-bonding/delamination with CZM is not yet state-of-art in FE softwares.

The purpose of this paper is to present a 2D combined finite‐discrete element method (FEM‐DEM) to model the multi‐fracture of beam structures and an application of the method to an ice‐structure interaction problem.

In the method, elastic beams and their fracture are modelled according to FEM by using nonlinear Timoshenko beam elements and cohesive crack model. Additionally, the beam elements are used to tie the discrete elements together. The contact forces between the colliding beams are calculated by using the DEM.

Three numerical examples are given to verify the method. Further, the method is applied to model the failure process of a floating ice beam against an inclined structure. Based on the comparison of the experiments and the simulation, a good agreement between the results is observed.

In the context of combined FEM‐DEM, the two novel features presented in this paper are: the use of Timoshenko finite element beams with damping to calculate internal forces and to combine the discrete elements; and the bending failure by the cohesive crack approach while simultaneously keeping track of the position of the neutral axis of the beam.

In certain cases, traction–separation laws do not reflect the behaviour sufficiently so that thin volumetric elements, Internal Thickness Extrapolation formulations, bulk material projections or various other approaches are applied. All of them have disadvantages in the formulation or practical application.

Damage within thin layers is often modelled using at cohesive zone elements (CZE). The constitutive behaviour of cohesive zone elements is usually described by traction–seperation laws (TSLs) that consider the (traction separation) relation in normal opening and tangential shearing direction. Here, the deformation (separation) as well as the reaction (traction) are vectorial quantities.

In this contribution, a CZE is presented that includes damage from membrane modes.

Membrane mode-related damaging effects that can be seen in physical tests that could not be simulated with standard CZEs are well captured by membrane mode–enhanced cohesive zone elements.

The aim of this paper is to present a novel data transfer technique to simulate, by G/XFEM, a cohesive crack propagation coupled with a smeared damage model. The efficiency of this technique is evaluated in terms of processing time, number of Newton–Raphson iterations and accuracy of structural response.

The cohesive crack is represented by the G/XFEM enrichment strategy. The elements crossed by the crack are divided into triangular cells. The smeared crack model is used to describe the material behavior. In the nonlinear solution of the problem, state variables associated with the original numerical integration points need to be transferred to new points created with the triangular subdivision. A nonlocal strategy is tailored to transfer the scalar and tensor variables of the constitutive model. The performance of this technique is numerically evaluated.

When compared with standard Gauss quadrature integration scheme, the proposed strategy may deliver a slightly superior computational efficiency in terms of processing time. The weighting function parameter used in the nonlocal transfer strategy plays an important role. The equilibrium state in the interactive-incremental solution process is not severely penalized and is readily recovered. The advantages of such proposed technique tend to be even more pronounced in more complex and finer meshes.

This work presents a novel data transfer technique based on the ideas of the nonlocal formulation of the state variables and specially tailored to the simulation of cohesive crack propagation in materials governed by the smeared crack constitutive model.

The purpose of this paper is to propose a novel combined finite-discrete element method (FDEM), based on the cohesive zone model, for simulating rockslide problems at the laboratory scale.

The combined FDEM is realized using ABAQUS/Explicit. The rock mass is represented as a collection of elastic bulk elements glued by cohesive elements with zero thickness. To reproduce the tensile and shear micro-fractures in rock material, the Mohr-Coulomb model with tension cut-off is employed as the damage initiation criterion of cohesive elements. Three simulated laboratory tests are considered to verify the capability of combined FDEM in reproducing the mechanical behavior of rock masses. Three slope models with different joint inclinations are taken to illustrate the application of the combined FDEM to rockslide simulation.

The results show that the joint inclination is an important factor for inducing the progressive failure behavior. With a low joint inclination, the slope failure process is observed to be a collapse mode. As the joint inclination becomes higher, the failure mode changes to sliding and the steady time of rock blocks is shortened. Moreover, the runout distance and post-failure slope angle decrease as the joint inclination increases.

These studies indicate that the combined FDEM performed within ABAQUS can simulate slope stability problems for research purposes and is useful for studying the slope failure mechanism comprehensively.

Since failure of laminated composites by delaminations is common, the purpose of this paper is to present a numerical procedure to check the stability of delaminations in fiber metal laminate (Glare), with different possible damage configurations, under uni-axial tension. Deformation behavior of the laminate is also examined. Influence of the type and the extent of damage, represented by varying sizes and number of delaminations, on delamination driving force and laminate deformation is found.

Delaminated Glare is modeled by finite element method. Interface cohesive elements are used to model the delaminations. Finite element results provide the deflection/deformation characteristics of the laminate. Driving forces of delaminations are estimated by J integrals that are numerically obtained over cyclic paths near delamination tips. Laminates with different types of delaminations are also fabricated and externally delaminated for measurement of their interlaminar fracture toughness. The delamination is considered to be stable if its driving force is less than corresponding interlaminar fracture toughness of the laminate.

Delaminations are found to be stable in laminates with lower number of delaminations and unstable in laminates with higher number of delaminations. Increase in size of delaminations increases the deformations but reduces the delamination driving force whereas increase in number of delaminations increases both deformations and driving forces. The trends change in case of laminates with symmetrical damage. Shape of delamination is also found to influence the deformations and driving forces. The finite element model is validated.

There is scope for validating the numerical results reported in the paper by theoretical models.

Checking the stability of delaminations and their effect on deformation behavior of the laminate helps is assessment of safety and remaining life of the laminate. If failure is predicted, preemptive action is taken by using repair patch ups at identified critical locations in order to avoid failures in service conditions.

The paper offers the following benefits: use of cohesive zone method that is readily possible in finite element procedures and is relatively simple, fast and reasonably accurate is demonstrated; suitability of using J integrals over paths crossing non-homogeneous and property mismatched material layers is tested; and influence of the type and the extent of damage in the laminate on its deformation behavior and delamination driving forces is found. This type of work has not been reported so far.

The purpose of this paper is to describe finite element modelling for fracture and fatigue behaviour of zirconia toughened alumina microstructures.

A two‐dimensional finite element model is developed with an actual Al₂O₃‐10 vol% ZrO₂ microstructure. A bilinear, time‐independent cohesive zone law is implemented for describing fracture behaviour of grain boundaries. Simulation conditions are similar to those found at contact between a head and a cup of hip prosthesis. Residual stresses arisen from the mismatch of thermal coefficient between grains are determined. Then, effects of a micro‐void and contact stress magnitude are investigated with models containing residual stresses. For the purpose of simulating fatigue behaviour, cyclic loadings are applied to the models.

Results show that crack density is gradually increased with increasing magnitude of contact stress or number of fatigue cycles. It is also identified that a micro‐void brings about the increase of crack density rate.

This paper is the first step for predicting the lifetime of ceramic implants. The social implications would appear in the next few years about health issues.

This proposed finite element method allows describing fracture and fatigue behaviours of alumina‐zirconia microstructures for hip prosthesis, provided that a microstructure image is available.

This study aims to accurately simulate the tilting separation process of mask projection stereolithography (MPSL) and verify the tilting theory.

The finite element separation models of MPSL 3D printing process were established. The established models simulated both tilting and pulling-up separation process by changing the constraints and boundary conditions. The bilinear cohesive curves were used to define the separation interface. The stress distribution of the cured part and FEP film at different times during the whole separation process was extracted. Different orientations of pulling-up and tilting were also compared for stress distribution. The stress change was analyzed for the center and edge points of the upper surface of cured part.

The results showed that the stress increased with the separation speed, and the stress at the edge position of exposure area was greater than the internal position. The tilting traction stress distribution was affected by the exposure area function and the velocity distribution. Alternation of the exposure area function changed the cohesive stiffness. The non-coincidence of the calculated traction stress with the input bilinear cohesive curve reflected the influence of the material properties and the separation methods. The high-speed side of tilting had fast separation and high traction stress.

This study proposes a technical method for simulation tilting separation and verified the tilting theory. The cohesive zone model was proved applicable to the tilting traction stress calculation.

The purpose of this paper is to assess a numerical tool to simulate and predict the onset and the propagation of the delaminations in a composite structure.

The approach to the work is done through the cohesive zone model technique applied to the finite element method.

Double cantilever beam, end notched flexure and mixed mode bending tests have been performed and correlated to benchmark cases, in order to validate the procedure. Numerical test campaign on specimens of the skirts with delaminations has been performed to analyze the behaviour under compressive load and the buckling.

This tool is applied to the study of the behaviour of some components in carbon/epoxy composite of a space structure in which one or more delaminations are eventually present following impact damage or manufacturing process. The components in particular are the booster's skirts of a small class launcher, subjected to a compressive load.

The parameters of cohesive elements have to be chosen correctly in the simulation of composite delamination by finite element method: such as interface strength, interface stiffness and shape of cohesive law. The purpose of this work is to investigate their influence on the accuracy of the results obtained. A three-dimensional cohesive-zone model has been established using Ls-dyna to simulate Double-Cantilever-Beam mode I (DCB) and Edge-Notched-Flexure mode II (ENF) tests. The influence of these parameters of cohesive element on the maximum load and the slope of load-displacement curve have been discussed by comparing experimental and numerical results. Four traction-separation laws: bilinear, linear-parabolic, exponential and trapezoidal were considered associated with a large variation of the interface strength and of the initial interface stiffness.