interface cohesive elements for fracture modeling

29
Crack modelling using zero-thickness interface elements Nguyen Vinh Phu The University of Adelaide

Upload: nguyen-vinh-phu

Post on 14-Jul-2015

219 views

Category:

Engineering


6 download

TRANSCRIPT

Page 1: Interface cohesive elements for fracture modeling

Crack modelling using zero-thickness interface elements

Nguyen Vinh Phu !

The University of Adelaide

Page 2: Interface cohesive elements for fracture modeling

IntroductionThere exists problems in which the crack path is known in advance...

inter-granular cracking in polycrystals

delaminated composites

Page 3: Interface cohesive elements for fracture modeling

Introduction (cont.)

• easy to implement (2D, 3D)

• available in major FE packages: ABAQUS, LS-DYNA

...then the use of interface elements is recommended

Alternatives: XFEM/GFEM

Page 4: Interface cohesive elements for fracture modeling

Cohesive crack model

separation

Constitutive equations

deformation

Governing equations (strong form)

Page 5: Interface cohesive elements for fracture modeling

Cohesive crack modelWeak form

new term

(skipped for static problems)

where

Page 6: Interface cohesive elements for fracture modeling

Discrete equationsupper face

lower face

Solid elements

Discrete equation

N1 = 0.5(1� ⇠)

N2 = 0.5(1 + ⇠)

Page 7: Interface cohesive elements for fracture modeling

Discrete equations (cont.)assembled

f ext = f int + f coh

Static problems

Linearization (Newton-Raphson)transformation matrix

Page 8: Interface cohesive elements for fracture modeling
Page 9: Interface cohesive elements for fracture modeling

Common interface elements

2D

3D

Solid elements Q4/T3 Q8/T6

H8

Page 10: Interface cohesive elements for fracture modeling

Numerical integrationIt has been observed numerically that integrating the internal force and stiffness matrix of interface elements using the standard Gauss rule led to oscillatory response [de Borst, IJNME, 1993].

Newton-Cotes integration scheme for interface elements

Page 11: Interface cohesive elements for fracture modeling

Cohesive lawsMode I Bilinear cohesive law (traction-separation law)

k

tensile strength

fracture energy

crack initiation

crack propagation

GIc =1

2tc�f =

1

2k�c�f

�f =2GIc

k�c�c =

tck

t = (1� d)k�

d =

8><

>:

�f

� if loading

max

�f

max

(�f

� �c

) + �c

otherwise

=� � �c�f � �c

elastic stiffness

Page 12: Interface cohesive elements for fracture modeling

• How to generate interface element meshes?

• Solution control - load control: - displacement control - path following control (arc-length methods)

Implementation aspects

F

u

cannot pass the snap-through

NO

cannot pass the snap-back

Page 13: Interface cohesive elements for fracture modeling

Mesh generation

• read a Gmsh mesh file

• double nodes along a path defined by the user

• modify the solid elements involved and

• generate interface elements

A C++ code was written to

Page 14: Interface cohesive elements for fracture modeling

Some application examples

Page 15: Interface cohesive elements for fracture modeling

Path following method

Newton-Raphson

f ext = �g

where

uI = K�1r, uII = K�1g

arc-length/constraint function�(u,�)

correction

Riks 1972 load factor�

reference load vector

Page 16: Interface cohesive elements for fracture modeling

Energy control

V =1

2

Z

⌦✏T� =

1

2

Z

⌦aTBT� =

1

2aTf int =

1

2�aTg

V � �aTgEnergy release rate

Gutierrez 2004 equilibriumf int =

Z

⌦BT�✏ = Ba

Arc-length function

forward Euler

predefined amount of energy to be released [Nm]

G > 0

Page 17: Interface cohesive elements for fracture modeling

Assumption: secant unloading!!!

�E

straight line

Page 18: Interface cohesive elements for fracture modeling

imagine what if there are 2 cracks???

Indirect displacement control [de Borst 1986]

SEN beam

Indirect displacement control

local quantity!!!�(||uA � uB || ,�l)

Page 19: Interface cohesive elements for fracture modeling

Energy control

- fast to evaluate - global quantity

Advantages of energy control

SEN beam

multiple cracksFP van der Mer

EFM, 2008

Page 20: Interface cohesive elements for fracture modeling

Solution procedure

load control

energy control

FP van der Mer EFM, 2008

Page 21: Interface cohesive elements for fracture modeling

f int + f coh

Page 22: Interface cohesive elements for fracture modeling

Numerical examples

• Simple tests (to debug code)

• Material interface debonding

• Multi-delamination of a composite DCB

• Delamination of the DCB

Page 23: Interface cohesive elements for fracture modeling

Debonding of a material interface

continuum: T3 elementsinterface: linear 4-node

Cohesive law: Xu-Needleman

Page 24: Interface cohesive elements for fracture modeling

Simple test

plane strain

Page 25: Interface cohesive elements for fracture modeling

Multi-delamination

Page 26: Interface cohesive elements for fracture modeling

Multi-delamination

Page 27: Interface cohesive elements for fracture modeling

Matrix kinking

interface elements everywhere except for the hard inclusion

“Discontinuous Galerkin/extrinsic cohesive zone modeling: Implementation caveats and applications in

computational fracture mechanics”, VP Nguyen, Engineering Fracture Mechanics, 2014.

Page 28: Interface cohesive elements for fracture modeling

Final remarks• Basic of interface elements was introduced

• Energy control path-following method was presented

When crack path is not known...

interface elements inserted a priori at every element edges!!!

Page 29: Interface cohesive elements for fracture modeling

Things to explore• New cohesive laws

• New (better) interface element formulations (current element technology does not allow industrial applications to be realized.)

• Mesh topology such that mesh bias can be avoided.