Research article Topical Sections

Cohesive delamination and frictional contact on joining surface via XFEM

  • In the present paper, the complex mechanical behaviour of the surfaces joining two different bodies is analysed by a cohesive-frictional interface constitutive model. The kinematical behaviour is characterized by the presence of discontinuous displacement fields, that take place at the internal connecting surfaces, both in the fully cohesive phase and in the delamination one. Generally, in order to catch discontinuous displacement fields, internal connecting surfaces (adhesive layers) are modelled by means of interface elements, which connect, node by node, the meshes of the joined bodies, requiring the mesh to be conforming to the geometry of the single bodies and to the relevant connecting surface. In the present paper, the Extended Finite Element Method (XFEM) is employed to model, both from the geometrical and from the kinematical point of view, the whole domain, including the connected bodies and the joining surface. The joining surface is not discretized by specific finite elements, but it is defined as an internal discontinuity surface, whose spatial position inside the mesh is analytically defined. The proposed approach is developed for two-dimensional composite domains, formed by two or more material portions joined together by means of a zero thickness adhesive layer. The numerical results obtained with the proposed approach are compared with the results of the classical interface finite element approach. Some examples of delamination and frictional contact are proposed with linear, circular and curvilinear adhesive layer.

    Citation: Francesco Parrinello, Giuseppe Marannano. Cohesive delamination and frictional contact on joining surface via XFEM[J]. AIMS Materials Science, 2018, 5(1): 127-144. doi: 10.3934/matersci.2018.1.127

    Related Papers:

    [1] Nikolay A. Voronin . Specialties of deformation and damage of the topocomposite on a ductile substrate during instrumental indentation. AIMS Materials Science, 2020, 7(4): 453-467. doi: 10.3934/matersci.2020.4.453
    [2] Shashi Bahl, Tarunpreet Singh, Virinder Kumar, Shankar Sehgal, Ashok Kumar Bagha . A systematic review on recent progress in advanced joining techniques of the lightweight materials. AIMS Materials Science, 2021, 8(1): 62-81. doi: 10.3934/matersci.2021005
    [3] Rafael Resende Lucas, Rita de Cássia Mendonça Sales-Contini, Luis Felipe Barbosa Marques, Jonas Frank Reis, Ana Beatriz Ramos Moreira Abrahão, Edson Cocchieri Botelho, Rogério Pinto Mota . Characterization of the hybrid joint between AA2024-T3 alloy and thermoplastic composite obtained by oxy-fuel welding (OFW). AIMS Materials Science, 2024, 11(3): 585-601. doi: 10.3934/matersci.2024029
    [4] Clara Argerich, Ruben Ibáñez, Angel León, Anaïs Barasinski, Emmanuelle Abisset-Chavanne, Francisco Chinesta . Tape surface characterization and classification in automated tape placement processability: Modeling and numerical analysis. AIMS Materials Science, 2018, 5(5): 870-888. doi: 10.3934/matersci.2018.5.870
    [5] Sarah Trinschek, Karin John, Uwe Thiele . From a thin film model for passive suspensions towards the description of osmotic biofilm spreading. AIMS Materials Science, 2016, 3(3): 1138-1159. doi: 10.3934/matersci.2016.3.1138
    [6] Nikolay A. Voronin . Analysis of the mechanisms of deformation of topocomposites by modeling of the indentation load-displacement curves. AIMS Materials Science, 2019, 6(3): 397-405. doi: 10.3934/matersci.2019.3.397
    [7] Yuan Guo, Xiong (Bill) Yu . Characterizing the surface charge of clay minerals with Atomic Force Microscope (AFM). AIMS Materials Science, 2017, 4(3): 582-593. doi: 10.3934/matersci.2017.3.582
    [8] Claude Stolz . Moving surfaces and interfaces : application to damage, fracture and wearcontact. AIMS Materials Science, 2016, 3(3): 881-907. doi: 10.3934/matersci.2016.3.881
    [9] Mica Grujicic, Jennifer S. Snipes, S. Ramaswami . Polyurea/Fused-silica interfacial decohesion induced by impinging tensile stress-waves. AIMS Materials Science, 2016, 3(2): 486-507. doi: 10.3934/matersci.2016.2.486
    [10] Mohamed Samy El-Feky, Passant Youssef, Ahmed Maher El-Tair, Sara Ibrahim, Mohamed Serag . Effect of nano silica addition on enhancing the performance of cement composites reinforced with nano cellulose fibers. AIMS Materials Science, 2019, 6(6): 864-883. doi: 10.3934/matersci.2019.6.864
  • In the present paper, the complex mechanical behaviour of the surfaces joining two different bodies is analysed by a cohesive-frictional interface constitutive model. The kinematical behaviour is characterized by the presence of discontinuous displacement fields, that take place at the internal connecting surfaces, both in the fully cohesive phase and in the delamination one. Generally, in order to catch discontinuous displacement fields, internal connecting surfaces (adhesive layers) are modelled by means of interface elements, which connect, node by node, the meshes of the joined bodies, requiring the mesh to be conforming to the geometry of the single bodies and to the relevant connecting surface. In the present paper, the Extended Finite Element Method (XFEM) is employed to model, both from the geometrical and from the kinematical point of view, the whole domain, including the connected bodies and the joining surface. The joining surface is not discretized by specific finite elements, but it is defined as an internal discontinuity surface, whose spatial position inside the mesh is analytically defined. The proposed approach is developed for two-dimensional composite domains, formed by two or more material portions joined together by means of a zero thickness adhesive layer. The numerical results obtained with the proposed approach are compared with the results of the classical interface finite element approach. Some examples of delamination and frictional contact are proposed with linear, circular and curvilinear adhesive layer.


    1. Introduction

    A problem of significant interest and importance in solid mechanics is the analysis of delamination phenomena in joined solids. In this sense, composite laminates, materials with inclusions, masonry structures with mortar beds, may be considered as typical joined bodies. The structural response of joined solids strongly depends on the mechanical behaviour of the internal adhesive layers, where several microstructural failure mechanisms lead to progressive degradation of the mechanical properties, up to complete delamination (macrocrack). Moreover, under compressive loading conditions, macrocraks may close and produce frictional effects.

    Adhesive layers are generally modelled as zero thickness surfaces by means of interface mechanical device, whose internal constitutive law is defined as a relation between traction and jump-displacement components.

    Interface constitutive modelling may be formulated at mesoscale level, reproducing the average mechanical behaviour of an adhesive material portion, known as representative surface element (RSE) whose size is sufficiently large with respect to the microstructure characteristic dimensions.

    In the literature, interface constitutive laws have been initially used in the pioneering cohesive-zone models of Dugdale and Barenblatt [1,2], to describe the progressive separation process in the fracture process zone (FPZ). These cohesive-zone models are defined in the framework of fracture mechanics, possibly coupled with damage mechanics or with plasticity theory. Interface models in the analysis of composite delamination are developed in the damage mechanics framework [3,4,5,6,7,8,9,10,11], where plastic and frictional phenomena are neglected. Alternative formulations [12,13,14] are developed in the framework of non-associative softening-plasticity theory, allow to catch the residual frictional behaviour of the fully cracked material.

    Cohesive damaging with residual frictional effects has been initially proposed in [15] for composite laminates and improved in [16,17]. In [18] and, by the same author in [19,20,21], an original mesoscale interpretation of the damage parameter allows to model the transition from the pure cohesive damaging behaviour of the sound material to the pure frictional behaviour of the fully cracked one. In [20], the cohesive-frictional constitutive model has been successfully used for the analysis of frictional effects in Mode Ⅱ delamination tests of composite material, focusing the attention on experimental tests, which were been performed with an increasing cycling loading path. In the present paper the theoretical and computational issues, related to the use of a cohesive-frictional model combined with the kinematic formulation of XFEM, are analysed in detail.

    Most of the proposed constitutive models assume the initial behaviour of interface as elastic, by a linear relationship between traction and jump displacement. The normal and tangential elastic parameters are defined in function of the elastic behaviour of the adhesive medium keeping joined the two bodies. For problems of joined solids without any adhesive layer, where bodies are directly bonded one each other, the initial elastic behaviour is assumed as an approximate penalty approach to numerically simulate the mechanical constrain between the bodies.

    From the kinematical point of view, the joining surface is modelled as jump displacement between the edges of the joined bodies to describe opening and relative sliding phenomena in the delamination phase. Jump displacements can be active also in the interface sound condition and considered as elastic deformation, which can produce material penetration in compressive state, or elastic delamination in tensile state. In order to bound such effects, a high initial elastic stiffness has to be assumed, but it can produce numerical problems such as spurious traction oscillations, as shown in [22].

    In the framework of FEM resolution strategy, adhesive layers are generally modelled by means of interface elements, which are defined in a zero thickness domain (surface element for 3D mechanical problems or linear element for 2D mechanical problems) connecting, node by node, the meshes of the joined bodies. Interface finite elements are of simple implementation but, as a drawback, they require two different meshes for the connected bodies with coincident nodes on the connecting surface. Generally, the domain discretization problem by use of interface elements can be quit burdensome.

    An alternative formulation is the extended finite element method (XFEM), which allows to define discontinuous displacement fields inside a generic domain, without the need of the mesh to be conforming to the discontinuity surface. The standard displacement based formulation is enriched by a discontinuous displacement field, which is governed by additional degrees of freedom. The spatial position of displacements discontinuity is implicitly defined by the level set method and it can be arbitrarily positioned inside the mesh.

    The XFEM has been developed by Belytschko el al. [23,24] to model crack formation and propagation, voids and inhomogeneity. Several XFEM schemes have been developed and proposed in literature: for two dimensional linear elastic fracture mechanics (LEFM) in [23], with the step function enrichment in [24], for three dimensional crack propagation in [25], for material discontinuity problems in [26,27], with second-order XFEM formulation in [27], for cohesive crack models in [28,29,30]. The XFEM formulation has also been proposed for the modelling of physically existing discontinuity surface. In [31], XFEM and level set method are analysed for the modelling of interface material failure with a damage separation law for the interface, neglecting frictional effects. In [32], the XFEM formulation has been employed to model the frictional effects on the contact surface between two adjacent bodies, but neither decohesion nor delamination have been considered. In [33], the material microstructure has been modelled by XFEM formulation, including voids and material interface, in a finite element multiscale resolution strategy. In [34], material interface has been modelled in a higher order XFEM framework, including curved discontinuity surface, for which attention is focused on the integration procedure and convergence of numerical solutions.

    The main target of the present paper is the modelling, in a unified framework, of several mechanical problems such as: frictional contact, delamination, decohesion and smooth transition from the initial cohesive behaviour up to the frictional residual on. Such a complex mechanical behaviour can be modelled, for a joining or discontinuity surface, in the kinematic framework of XFEM, for which the interface can be geometrically placed arbitrarily in a unique mesh, without the requirement of defining two or more meshes with coincident nodes at the joining surface.

    The implemented constitutive model has been developed by the same author in [19]. The solution proposed in the present paper is limited to two-dimensional problems for which the connection surface is geometrically represented by a line. The finite element formulation is based on the triangular six nodes eXtended element published in [30], which allows to model displacement fields with a strong discontinuity at a generic line inside the mesh.


    2. Extended finite element problem

    Consider two solid domains Ωa and Ωb perfectly bonded and in mutual contact on the common internal surface Γc, with applied external tractions t on the free boundary Γt and with imposed displacement ˉu on Γu, as represented in Figure 1a in a two-dimensional case. The hypotheses of homogeneous material and small strain are assumed. In Figure 1b the overall domain Ω=ΩaΩb is discretized in triangular finite elements and the mesh does not conform the internal connecting surface Γc.

    Figure 1. a) domain and boundary of the joined solids; b) mesh not conforming to internal surface Γc.

    The kinematic formulation developed in [30] allows to model a discontinuous displacement field for both the partially cracked finite element ("crack tip element") and for the completely cracked finite element. Moreover both a first-order scheme, with a three node triangular element and a second-order scheme, with a six node triangular element, are developed.

    The displacement field of the whole domain Ω is defined as the additive combination of a continuous part uc and a discontinuous one ud :

    u(x)=uc(x)+ud(x) (1)

    where the continuous part is defined by standard shape functions

    uc(x)=iINi(x)ui (2)

    and I is the set of all nodes of the meshed domain and ui is the displacement vector of the i-th node. With reference to domain Ωenr, defined by the set of elements entirely intersected by the connecting surface Γc, and the relevant set on nodes IenrI, the discontinuous part of the displacement field is written in the following form

    ud(x)=iIenrNi(x)Ψi(x)ai (3)

    where Ψi(x) is the enrichment function and ai contains the relevant degree of freedoms.

    The spatial position of the connecting surface Γc is defined by the level set method, that is

    Γc={xΩ:ϕ(x)=0} (4)

    where function ϕ(x) assumes different signs in the two domain Ωa and Ωb and whose iso-zero curve coincides with surface Γc, as shown in Figure 1a. A local reference with tangential and normal axes (t,n) is defined at a generic point of the connecting surface Γc and the normal axis n is oriented towards the positive values of level set function.

    For elements completely intersected by the surface Γc, the enrichment function is defined as

    Ψi(x)=12sign(ϕ(x))12sign(ϕi),  iIenr,  xΩenr (5)

    where ϕi=ϕ(xi) is the level set value at the i-th node.

    Due to the specific form of the enrichment functions Ψi(x), the enrichment displacement field vanishes outside the enriched elements, as shown in Figure 2 for a three node triangular element.

    Figure 2. Enrichment shape functions for a three node triangular element crosses by the connecting surface from side to side.

    In the finite element framework the interface deformation is defined as a jump displacement between the adjacent edges of the two connected bodies, lying on the surface Γc. The jump displacement across the connecting surface Γc is defined as the difference of displacement between the two adjacent points on the positive and on the negative part of Γc, that is

    [[u(x)]]=u(x+)u(x),    xΓc (6)

    where x+ and x are points respectively on the positive side of Γc and on its negative side, with reference to the sign of the level-set function. By substitution of Eqs 1, 2, 3 and 5 in Eq 6, it can be observed that jump displacement is function of the enriching parameters ai, that is

    [[u(x)]]=iIenrNi(x)ai,    xΓc (7)

    In the present paper, the initial behaviour of the interface is assumed elastic and the jump displacements can be different than zero also in the undamaged state of the connecting surface. The enriching degrees of freedom ai with iIenr, which belong to nodes contained in the domain Ωenr of elements crossed by the interface, are always unconstrained and some numerical problems, such as spurious traction oscillations, can emerge (see [22]).

    In [31], two different formulations are considered respectively for the sound interface and for the cracked one. In the former case, the finite element model is enhanced with a weak discontinuity field, in order to model the kinematical behaviour at the initial state, when the two bodies are perfectly bonded with zero jump displacements and with independent deformation states on the two element parts. Finally, when in a finite element the internal interface attains the delamination condition, the element is enriched with a strong discontinuity field by releasing the relevant degrees of freedom. Nevertheless, in such formulation numerical problem of stress singularity or high stress concentration are observed.


    2.1. Position of surface discontinuity

    The position inside the mesh of the connecting surface is defined in the input file by means of linear or quadratic parametric equations, or circle equations. Inside a single finite element the connecting surface is linearised and the following intersection schemes are considered:

    a) side to side, when connecting surface intersects two element sides (Figure 2);

    b) node to side, when connecting surface contains one node and intersects its opposite side;

    c) node to node, when connecting surface contains two vertex nodes;

    d) one node, when connecting surface contains just one vertex node, without cutting the element internal domain.

    In (a) and (b) schemes, the finite element is completely cut, whereas in scheme (c) the interface coincides with a finite element side.


    3. Cohesive-frictional constitutive model for interface

    The interface constitutive model proposed in [19] is adopted for the numerical analysis; the model is developed in two-dimensional space, with reference to a zero thickness linear interface and it is described in the present section.

    Interface constitutive models typically link the displacement discontinuity [[u(x)]] and the adhesion forces t(x). As described in previous section, the kinematics is defined by the displacement discontinuity across the interface, and it can be decomposed in two components: the normal one (Mode Ⅰ) and the tangential one (Mode Ⅱ). The interface equilibrium condition imposes continuity of traction t across interface, therefore it is assumed t=t+=t, where t+ and t are the traction vectors, respectively, on the positive edge and on the negative one. Statics and kinematics of an infinitesimal element dS of the connecting surface Γc are represented in Figures 3a and b.

    Figure 3. a) Statics and b) kinematics of an infinitesimal element dS of the connecting surface Γc.

    The proposed constitutive model is developed in the damage mechanics framework [35] considering the Representative Surface Element (RSE), whose size is assumed sufficiently large compared to the adhesive material inhomogeneities (Figure 4), and average damage in the RSE is defined as

    ω=dScdS (8)

    where dSc measures the area of the micro-cracked RSE portion and dS measures the RSE area.

    The damage variable ω measures the specific areal crack density of the infinitesimal surface element dS and it allows us to define, at the material point, the extensions of the two complementary fractions: the micro-cracked one ωdS and of the sound one (1ω)dS (see Figure 4). Two different constitutive laws can be formulated for the two fractions and the damage variable governs the progressive transition of the adhesive constitutive model from the initial elastic-cohesive behaviour to the residual frictional one. In particular, a cohesive linear elastic law is assumed for the sound fraction and a non-associative elastic-plastic constitutive relation, with Coulomb yield surface, is assumed for the cracked fraction. The compressive behaviour of both the sound fraction and the damaged one, is modelled as linearly elastic.

    Figure 4. Damage meso-scale level interpretation.

    The constitutive model formulation is developed considering the index s for the sound fraction variables and index c for the micro-cracked fraction ones. For the sound fraction the elastic deformation is define as

    δes=[[u]]δps (9)

    where δps measures the plastic contribution. For the cracked fraction, which is represented in Figure 5a, the variable δec measures the asperities elastic deformation, represented in Figure 5b, and it is defined by the following relation

    Figure 5. Meso-scale level representation of frictional phenomena: a) cracked fraction; b) micro-cracked fraction elastic deformation; c) micro-cracked fraction elasto-plastic deformation; d) interface edges detachment.

    where δpc measures the plastic component (sliding and dilatancy represented in Figure 5c) and δdc is the detachment jump displacement (Figure 5d).

    The vector δpc represents the jump displacement caused by the irreversible sliding of asperities among each others and it develops when the limit of frictional strength is attained. Moreover, since frictional behaviour is active only for compressive normal traction, positive normal traction produces loss of contact between the crack edges, which is measured by the detachment displacement δdc.

    The constitutive model is developed in a thermodynamic framework following the classic Coleman and Noll procedures (for more detail see [35,36,37]). The Helmholtz free energy density (for unit interface length) is defined as

    Ψ=Ψel(δes,δec,ω)+Ψin(ξ)==12(1ω)δeTsKsδes+12ωδeTcKcδec+Ψin(ξ) (11)

    where Ψel and Ψin are respectively the elastic part and the internal one; Ks and Kc are diagonal elastic matrices, respectively, of the sound fraction and of the micro-cracked one; ξ is a scalar state internal variable, which governs damage hardening. The second thermodynamic principle, in the form of Clausius-Duhem inequality, imposes a pointwise positiveness of rate dissipation energy density given as the difference between the work done by external forces and the variation of the free energy, for any change of the state variables. This condition is mathematically written as

    D=tT[[˙u]]˙Ψ0, (12)

    where rate of Helmholtz free energy is defined as

    ˙Ψ=Ψelω˙ω+Ψelδes˙δes+Ψelδec˙δec+Ψinξ˙ξ==(Y˙ωtTs˙δestTc˙δecχ˙ξ) (13)

    ts and tc are conjugated variable of the elastic deformations δes and δec, that is

    ts=Ψelδes=(1ω)Ks (14)
    tc=Ψelδec=ωKcδec, (15)

    Y is the energy release rate and χ is the static hardening-like variable

    Y=Ψelω=12δeTsKsδes12δeTcKcδec (16)
    χ=Ψinξ. (17)

    For an elastic loading step, where ˙ω=0, ˙ξ=0, ˙δps=0 and ˙δpc=0, rate elastic deformations are ˙δes=˙δec=[[˙u]], dissipation has to be null (D=0) and after substitution Eq 12 reads

    D=(ttstc)T[[˙u]]=0. (18)

    Since Eq 18 has to be verified for any elastic loading step [[˙u]], the following internal equilibrium condition is obtained

    t=ts+tc (19)

    and, after substitution of Eqs 14 and 15, the overall interface elastic traction-displacement relation is obtained

    t=(1ω)Ksδes+ωKcδec. (20)

    Equation 20 is assumed to be valid also for inelastic loading steps, in which dissipation is

    D=tTs˙δps+tTc˙δpc+Y˙ωχ˙ξ0. (21)

    where detachment δdc does not produce any dissipation.

    Dissipation positiveness is ensured by means of suitable activation criteria, which govern the evolution of plastic and damaging phenomena. In particular, damage evolution is governed by the following yield function

    ϕd(Y,χ)=Yχ(ξ)Y00, (22)

    with Y0 assumed as initial yielding threshold. The relevant flow rules with the loading-unloading conditions read

    ˙ω=ϕdY˙λd=˙λd,˙ξ=ϕdχ˙λd=˙λd,˙λd0,ϕd˙λd=0,˙ϕd˙λd=0. (23)

    The damage hardening-like law is

    χ(ξ):=12KsNˉu2e[(ˉufˉuf(1ξ)+ˉueξ)21], (24)

    where ˉue and ˉuf are jumps displacement limit values, respectively, at the elastic threshold and at the unitary damage condition in a pure tensile state and KsN is the normal elastic stiffness of the interface sound fraction. Hardening law in Eq 24 produce linear softening in the traction-jump displacement space and it produces the following simple Mode Ⅰ fracture energy

    GI=12KsNˉueˉuf. (25)

    Plastic phenomena in the sound fraction are neglected (δps=0), whereas, in micro-cracked fraction, plastic phenomena evolution is governed by means of the classical Coulomb yield function, generally used for frictional materials,

    ϕp=|tcT|+αtcN0. (26)

    under the hypothesis of non-associative plasticity, by means of the following plastic potential

    Ωp=|tcT|+βtcN, (27)

    where α and β, with αβ, are respectively the frictional coefficient and the dilatancy one. The plastic flow rules and loading/unloading conditions are

    ˙δpcT=ΩptcT˙λp=sgn(tcT)˙λp,˙δpcN=ΩptcN˙λp=β˙λp,˙λp0,ϕp˙λp=0,˙ϕp˙λp=0. (28)

    By substitution of Eqs 28 and 23 in Eq 21, dissipation assumes the following form

    D=(|tcT|+βtcN)˙λp+Y0˙λd0 (29)

    and positiveness is ensured for any damaging or plastic loading step. The cohesive-frictional interface constitutive model has been implemented in the open source finite element code FEAP [38]. The interface kinematics is implicitly defined inside the triangular six nodes extended finite element. The enrichment degrees of freedom are active for the nodes of all the finite elements crossed the joining surface ΓC; on the contrary, the other enrichment degrees of freedom are restrained. Three Gauss integration points are considered for each finite element, where the two components of the jump displacement are computed through Eq 7 and passed to the subroutine of the interface constitutive model, which give back the tractions components and the constitutive stiffness matrix. Details on the implementation of both user defined finite elements and user defined constitutive models can be found in [38].


    4. Numerical simulations

    The proposed interface constitutive model and the triangular extended finite element have been implemented in the open source finite element code FEAP [38]. The numerical results of delamination phenomenon in joined bodies are presented.

    The first simulation is the two-dimensional analysis of delamination between a square matrix the the circular inclusion, represented in Figure 6 with dimensions and boundary conditions. Matrix and inclusion are considered as elastic isotropic with the following constitutive parameters: EM=40000 N/mm2, νM=0.2, EI=400000 N/mm2, νI=0.2. With reference to the proposed interface constitutive model, the elastic parameters of the sound fractions are: KsN=400000 N/mm,KsT=347830 N/mm; the elastic parameters of the cracked fractions are: KcN=400000 N/mm,KcT=6956.6 N/mm; the jump displacement limit values are ˉue=106 mm,ˉuf=5×105 mm; frictional and dilatancy coefficients are α=0.364,β=0.176.

    Figure 6. Matrix with circular inclusion: a) numerical simulation scheme; b) mesh with interface elements; c) mesh with XFEM elements.

    The domains of the two joined solids are numerically modelled by the not conforming mesh of XFEM triangular elements represented in Figure 6c; they are also modelled by quadrilateral isoparametric elements and interface elements in the mesh represented in Figure 6b, which conforms to the joining surface between matrix and inclusion. The numerical responses carried out by the two approaches, XFEM and classical interface elements, are compared in terms of force vs displacement in Figure 7 and the results are numerically coincident. Figures 8a and b plot the maps of vertical normal stress carried out, respectively, at the reference loading steps a and b represented in Figure 7.

    Figure 7. Comparison of the numerical responses obtained by XFEM and by interface element of the square matrix with circular inclusion.
    Figure 8. Vertical normal stress maps for the square matrix with circular inclusion at the reference loading steps a) and b).

    The second proposed simulation is the two-dimensional sliding test between two elastic solids in mutual contact along a curvilinear surface, without any cohesive strength. The configuration test is represented in Figure 9a together with the adopted mesh, which is defined for the whole domain and it does not conform the contact surface. The hypothesis of small displacement is assumed. In the considered sliding test, the lower elastic solid is constrained along the lower boundary side and the upper elastic solid is subjected to a constant vertical force Fy=1000 N and to a monotonically increasing horizontal displacement at the upper boundary side, as shown in Figure 9a. Rotation of the upper boundary side is restrained. The elastic parameters of the two solids are E=40000 N/mm2, ν=0.2; for the contacts surface, frictional and dilatancy coefficients respectively are α=tan20, β=0. The map of vertical stresses, obtained by the numerical simulation for the imposed horizontal displacement ux=3.0 mm, is plotted in Figure 9b. The results of the numerical simulation are plotted in Figure 10 in terms of horizontal displacement vs horizontal force. The second numerical simulation shows that the proposed formulation allows to model problems of frictional contact with a very simple mesh generation, where a quite complex contact surface is defined in the input file as an isoparametric quadratic function.

    Figure 9. Two-dimensional sliding test of two solid in mutual contact along a curvilinear surface: a) test configuration with the XFEM model; b) Vertical stress in the frictional sliding test at the loading condition corresponding to the imposed displacement ux=3.0 mm.
    Figure 10. Numerical response in terms of horizontal displacement vs horizontal force.

    The last numerical simulation reproduces the experimental results of a three points bending test performed on a carbon fibres-epoxy composite beam, whose dimensions are represented in Figure 11. The composite panel has been manufactured by hand lay-up and vacuum-bag techniques by a sequence of 12 unidirectional plies of carbon fibres. A Mylar foil was inserted at the mid-plane of the composite panels during the lay-up process in order to produce a starter delamination crack. Four specimens have been tested under displacement control and simultaneously the crack extension has been monitored by an high resolution camera. The experimental test results are plotted in Figure 12 in terms of imposed displacement vs vertical reaction. The fracture energy values obtained for the four experimental tests are reported in Table 1 with the relevant mean value. The three points bending test has been numerically simulated by using a two-dimensional mesh of 9×110 six nodes extended triangular elements in plane strain condition. The composite material is assumed as elastic isotropic with the following constitutive parameters: E=80580 N/mm2 and ν=0.23. The interface constitutive parameters are: KsN=KsT=90000 N/mm, KcN=90000 N/mm,KcT=3000 N/mm; the jump displacement limit values are ˉue=6.9×104 mm,ˉuf=0.014593 mm; frictional and dilatancy coefficients respectively are α=tan35, β=0. The adopted constitutive parameters produces the fracture energy value GII=12KsNˉueˉuf=0.4531 N/mm. The results of the numerical simulation are compared to the experimental results in Figure 12 in terms of imposed displacement vs load P, showing a good agreement between the two results. The numerical simulation produces: the same stiffness in the elastic branch; the same values of load and imposed displacement at the delamination condition; a response close to the experimental data also in the post-delamination branch. Moreover, the set of parameters chosen for the cohesive law produce a mode Ⅱ fracture energy very close to the values obtained by the experimental data. The discretization of the pre-cracked specimen, by XFEM formulation, has been obtained by the generation the mesh of the overall rectangular domain and by defining, in the input file, the line of the pre-existing crack, with unitary damage, and the line of possible delamination. The discretization of the same problem with standard interface element would be extremely much more onerous.

    Figure 11. Specimen size for the three points bending test.
    Figure 12. Specimen size for the three points bending test.
    Table 1. Fracture energy values of the four three points bending tests.
    Specimen number Fracture energy GII [N/mm]
    1 0.4706
    2 0.4737
    3 0.4198
    4 0.4481
    Mean value 0.4531
     | Show Table
    DownLoad: CSV

    Finally, the map of horizontal normal stress at the maximum loading condition is plotted in Figure 13 for the whole specimen together with a zoom around the crack tip.

    Figure 13. Map of σx stress component of the whole specimen at the maximum loading condition and a zoom around the crack tip.

    The first numerical simulation proposed in the present paper has been discretized either by the standard interface elements and by the extended finite elements and the results of the two numerical simulations are practically coincident, showing the effectiveness of the implemented XFEM formulation.

    Whereas interface finite elements are of very simple implementation in a finite element code, the most difficult in their use is the discretization phase, which requires a separate mesh for each connected body, with coincident nodes of the discretized bodies on the joining surface. The domain discretization phase, by use of interface elements, can be quite burdensome. Some alternative formulations of non-matching interface elements have been proposed in literature [39,40], which allow the discretization of the two connected bodies with independent meshes. The two meshes can be defined with non-coincident nodes on the joining surface ΓC, with an easier discretization phase.

    On the contrary, the principal advantage of the XFEM formulation in the modelling of different bodies connected by joining surface is the discretization. The whole domain can be discretized with a unique mesh and the position of joining surface can be analytically defined, with much lower time needed for such operation, with respect to the interface formulation. Moreover, the XFEM formulation allows also to model extrinsic cohesive law, with initial rigid behaviour. Disadvantages of the XFEM formulation is the implementation, which can be very onerous, and the mesh density of the joining surface, which cannot be easily controlled for each element.


    5. Conclusions

    The present paper proposes a unified approach for the modelling of friction and delamination phenomena on discontinuity surfaces, discretized in an overall mesh based on an XFEM formulation. The proposed approach allows us to model both the delamination phenomenon with crack propagation of two joined solids and the frictional contact phenomenon between two adjacent solids. The transition, from the initially cohesive behaviour with delamination phenomena, to the final behaviour with residual frictional strength is also considered.

    The discontinuity surfaces has been arbitrarily defined inside the mesh in parametric formulation, a unique mesh can be defined for whole domain, including connected bodies and discontinuity surface. On the contrary, in a standard approach, separated meshes have to be defined for the connected or adjacent bodies and specific interface finite elements for the discontinuity surface. Moreover, the meshes of the adjacent or connected bodies have to be defined with spatially coincident nodes on the discontinuity surface.


    6. Acknowledgments

    Funding: The financial support of the Italian Ministry for University and Research (MIUR), under the grant PRIN-2015, project No.2015LYYXA8, "Multiscale mechanical models for the design and optimization of microstructured smart materials and metamaterials" is gratefully acknowledged.


    7. Conflict of interest

    The authors declare that they have no conflict of interest.


    [1] Dugdale D (1960) Yielding of steel sheets containing slits. J Mech Phys Solids 8: 100–104. doi: 10.1016/0022-5096(60)90013-2
    [2] Barenblatt G (1962) The mathematical theory of equilibrium cracks in brittle fracture. Adv Appl Mech 7: 55–129. doi: 10.1016/S0065-2156(08)70121-2
    [3] Allix O, Blanchard L (2006) Mesomodeling of delamination: towards industrial applications. Compos Sci Technol 66: 731–744. doi: 10.1016/j.compscitech.2004.12.023
    [4] Allix O, Ladevéze P (1992) Interlaminar interface modeling for the prediction of delamination. Compos Struct 22: 235–242. doi: 10.1016/0263-8223(92)90060-P
    [5] Borino G, Fratini L, Parrinello F (2009) Mode I failure modeling of friction stir welding joints. Int J Adv Manuf Tech 41: 498–503. doi: 10.1007/s00170-008-1498-1
    [6] Corigliano A (1993) Formulation, identification and use of interface models in the numerical analysis of composite delamination. Int J Solids Struct 30: 2779–2811. doi: 10.1016/0020-7683(93)90154-Y
    [7] Mi Y, Crisfield MA, Davies GAO, et al. (1998) Progressive delamination using Interface elements. J Compos Mater 32: 1246–1272. doi: 10.1177/002199839803201401
    [8] Alfano G, Crisfield MA (2001) Finite element interface models for the delamination analysis of laminated composites: mechanical and computational issues. Int J Numer Meth Eng 50: 1701–1736. doi: 10.1002/nme.93
    [9] Qiu Y, Crisfield MA, Alfano G (2001) An interface element formulation for the simulation of delamination with buckling. Eng Fract Mech 68: 1755–1776. doi: 10.1016/S0013-7944(01)00052-2
    [10] Zou Z, Reid SR, Li S (2003) A continuum damage model for delamination in laminated composites. J Mech Phys Solids 51: 333–356. doi: 10.1016/S0022-5096(02)00075-3
    [11] Point N, Sacco E (1996) A delamination model for laminated composites. Int J Solids Struct 33: 483–509. doi: 10.1016/0020-7683(95)00043-A
    [12] Mortara G, Boulon M, Ghionna VN (2002) A 2-D constitutive model for cyclic interface behaviour. Int J Numer Anal Met 26: 1071–1096. doi: 10.1002/nag.236
    [13] Carol I, López CM, Roa O (2001) Micromechanical analysis of quasi-brittle materials using fracture-based interface elements. Int J Numer Meth Eng 52: 193–215. doi: 10.1002/nme.277
    [14] Cocchetti G, Maier G, Shen XP (2002) Piecewise linear models for interfaces and mixed mode cohesive cracks. CMES-Comp Model Eng 3: 279–298.
    [15] Tvergaard V (1990) Effect of fiber debonding in a whisker-reinforced metal. Mater Sci Eng 125: 203–213. doi: 10.1016/0921-5093(90)90170-8
    [16] Gambarotta L (2004) Friction-damage coupled model for brittle materials. Eng Fract Mech 71: 829–836. doi: 10.1016/S0013-7944(03)00020-1
    [17] Gambarotta L, Logomarsino S (1997) Damage models for the seismic response of brick masonry shear walls. Part I: the mortar joint model and its application. Earthq Eng Struct D 26: 423–439.
    [18] Alfano G, Sacco E (2006) Combining interface damage and friction in a cohesive-zone model. Int J Numer Meth Eng 68: 542–582. doi: 10.1002/nme.1728
    [19] Parrinello F, Failla B, Borino G (2009) Cohesive-frictional interface constitutive model. Int J Solids Struct 46: 2680–2692. doi: 10.1016/j.ijsolstr.2009.02.016
    [20] Parrinello F, Marannano G, Borino G, et al. (2013) Frictional effect in mode II delamination: Experimental test and numerical simulation. Eng Fract Mech 110: 258–269. doi: 10.1016/j.engfracmech.2013.08.005
    [21] Parrinello F, Marannano G, Borino G (2016) A thermodynamically consistent cohesive-frictional interface model for mixed mode delamination. Eng Fract Mech 153: 61–79. doi: 10.1016/j.engfracmech.2015.12.001
    [22] Simone A (2004) Partition of unity-based discontinuous elements for interface phenomena: computational issues. Int J Numer Meth Bio 20: 465–478.
    [23] Belytschko T, Black T (1999) Elastic crack growth in finite elements with minimal remeshing. Int J Numer Meth Eng 45: 601–620. doi: 10.1002/(SICI)1097-0207(19990620)45:5<601::AID-NME598>3.0.CO;2-S
    [24] Moës N, Dolbow J, Belytschko T (1999) A finite element method for crack growth without remeshing. Int J Numer Meth Eng 46: 131–150. doi: 10.1002/(SICI)1097-0207(19990910)46:1<131::AID-NME726>3.0.CO;2-J
    [25] Sukumar N, Moës N, Moran N, et al. (2000) Extended finite element method for threedimensional crack modelling. Int J Numer Meth Eng 48: 1549–1570. doi: 10.1002/1097-0207(20000820)48:11<1549::AID-NME955>3.0.CO;2-A
    [26] Belytschko T, Moës N, Usui S, et al. (2001) Arbitrary discontinuities in finite elements. Int J Numer Meth Eng 50: 993–1013. doi: 10.1002/1097-0207(20010210)50:4<993::AID-NME164>3.0.CO;2-M
    [27] Sukumar N, Chopp DL, Moës N, et al. (2000) Modelling holes and inclusions by level sets in the extended finite element method. Comput Method Appl M 190: 6183–6200.
    [28] Belytschko T, Moës N (2002) Extended finite element method for cohesive crack growth. Eng Fract Mech 69: 813–883. doi: 10.1016/S0013-7944(01)00128-X
    [29] Wells GN, Sluys LJ (2001) A new method for modelling cohesive cracks using finite elements. Int J Numer Meth Eng 50: 2667–2682. doi: 10.1002/nme.143
    [30] Zi G, Belytschko T (2003) New crack-tip elements for XFEM and applications to cohesive cracks. Int J Numer Meth Eng 57: 2221–2240. doi: 10.1002/nme.849
    [31] Hettich T, Ramm E (2006) Interface material failure modeled by the extended finite-element method and level sets. Comput Method Appl M 195: 4753–4767. doi: 10.1016/j.cma.2005.09.022
    [32] Khoei AR, Nikbakht M (2007) An enriched finite element algorithm for numerical computation of contact friction problems. Int J Mech Sci 49: 183–199. doi: 10.1016/j.ijmecsci.2006.08.014
    [33] Moës N, Cloirec M, Cartraud P, et al. (2003) A computational approach to handle complex microstructure geometries. Comput Method Appl M 192: 3163–3177. doi: 10.1016/S0045-7825(03)00346-3
    [34] Drau K, Chevaugeon N, Moës N (2010) Studied X-FEM enrichment to handle material interfaces with higher order finite element. Comput Method Appl M 199: 1922–1936. doi: 10.1016/j.cma.2010.01.021
    [35] Lemaitre J, Chaboche JL (1990) Mechanics of solids materials, Cambridge University Press.
    [36] Coleman BD, Noll W (1963) The thermodynamics of elastic materials with heat conduction and viscosity. Arch Ration Mech An 13: 167–178. doi: 10.1007/BF01262690
    [37] Coleman B (1971) Thermodynamics of Materials with Memory, CISM, Springer.
    [38] Zienkiewicz OC, Taylor RL (2000) The finite element method: solid mechanics, Butterworth-Heinemann.
    [39] Paggi M, Wriggers P (2016) Node-to-segment and node-to-surface interface finite elements for fracture mechanics. Comput Method Appl M 300: 540–560. doi: 10.1016/j.cma.2015.11.023
    [40] Nguyen VP, Nguyen CT, Bordas S, et al. (2016) Modelling interfacial cracking with non-matching cohesive interface elements. Comput Mech 58: 731–746. doi: 10.1007/s00466-016-1314-y
  • This article has been cited by:

    1. F. Parrinello, G. Borino, Non associative damage interface model for mixed mode delamination and frictional contact, 2019, 76, 09977538, 108, 10.1016/j.euromechsol.2019.03.012
    2. Francesco Parrinello, Vincenzo Gulizzi, Ivano Benedetti, A Model for Low-Cycle Fatigue in Micro-Structured Materials, 2019, 827, 1662-9795, 134, 10.4028/www.scientific.net/KEM.827.134
    3. Francesco Parrinello, Hybrid equilibrium element with interelement interface for the analysis of delamination and crack propagation problems, 2021, 122, 0029-5981, 190, 10.1002/nme.6531
    4. Francesco Parrinello, Guido Borino, 2020, Chapter 35, 978-3-030-41056-8, 419, 10.1007/978-3-030-41057-5_35
    5. Francesco Parrinello, Guido Borino, An extrinsic interface developed in an equilibrium based finite element formulation, 2019, 18, 24523216, 616, 10.1016/j.prostr.2019.08.207
    6. Francesco Parrinello, Hybrid Equilibrium Finite Element Formulation for Cohesive Crack Propagation, 2019, 827, 1662-9795, 104, 10.4028/www.scientific.net/KEM.827.104
    7. F. Parrinello, I. Benedetti, A coupled plasticity-damage cohesive-frictional interface for low-cycle fatigue analysis, 2022, 224, 00207403, 107298, 10.1016/j.ijmecsci.2022.107298
    8. Aazim Shafi Lone, G. A. Harmain, Azher Jameel, Modeling of contact interfaces by penalty based enriched finite element method, 2022, 1537-6494, 1, 10.1080/15376494.2022.2034075
    9. András Szekrényes, Differential quadrature solution for composite flat plates with delamination using higher-order layerwise models, 2022, 248, 00207683, 111621, 10.1016/j.ijsolstr.2022.111621
    10. András Szekrényes, Application of differential quadrature method to delaminated first-order shear deformable composite plates, 2021, 166, 02638231, 108028, 10.1016/j.tws.2021.108028
    11. Guido Borino, Francesco Parrinello, A symmetric tangent stiffness approach to cohesive mechanical interfaces in large displacements, 2022, 23, 1550-2287, 551, 10.1080/15502287.2022.2032481
  • Reader Comments
  • © 2018 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(5939) PDF downloads(906) Cited by(11)

Article outline

Figures and Tables

Figures(13)  /  Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog