application of the x-fem to the fracture of piezoelectric ... · 34 conclusion & future work...

35
1 Application of the X-FEM to the fracture of piezoelectric materials Eric Béchet 1 , Matthias Scherzer 2 , Meinhard Kuna 2 1 LPMM– CNRS UMR 7554 Université de Metz - France 2 IMFD – TU Bergakademie Freiberg - Germany

Upload: others

Post on 24-Jul-2020

4 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

1

Application of the X-FEM to the fracture of piezoelectric

materials

Eric Béchet1, Matthias Scherzer2, Meinhard Kuna2

1 LPMM– CNRS UMR 7554Université de Metz - France

2 IMFD – TU Bergakademie Freiberg - Germany

Page 2: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

2

Outline

● Introduction● Updated X-FEM formulation

– New crack tip enrichment functions

– Updated SIF computation scheme

– Use of specific preconditionner

● Convergence study● Conclusions

Page 3: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

3

Introduction

● Goal– Propose a updated enrichment scheme for a

cracked anisotropic piezoelectric media

– Convergence study of the method● Energy error● SIFS and energy release rate

– Development of a SIF evaluation scheme based on interaction integrals specific to piezoelectric materials

– Numerical crack propagation using empirical laws

Page 4: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

4

Introduction

● Physical model : linear piezoelectric media, electrically impermeable crack

Page 5: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

5

Introduction

● Numerical Model– Xfem field approximation

– No remeshing

– Interaction integrals used to compute the SIFs

Page 6: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

6

X-FEM

● Crack represented by level-sets– Local coodinates at the crack tip

lsn

lst

Page 7: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

7

X-FEM

– Singular functions around crack tip

– Heaviside along crack surface

– Remaining dofs unenriched

● Local partition of unity enrichment

Page 8: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

8

Enrichement functions

● Enrichment functions– Jump across the crack for displacements and

potential :

– Crack tip for in a pure mechanical setting

Page 9: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

9

Enrichement functions

Page 10: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

10

Enrichement functions– Crack tip functions for a

piezoelectrical setting● Must span the

eigenfunction's space at the crack tip for displacements and potential

● Depends on the material characteristics and the orientation

● Depends on the permeability of the crack

Page 11: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

11

Enrichement functions

Page 12: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

12

Enrichement functions

● Modified functions

Page 13: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

13

Enrichement functions

Page 14: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

14

Updated enrichment scheme

“topological” Enrichment

h=1/10 h=1/50

Page 15: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

15

Updated enrichment scheme

“Geometrical” Enrichment

h=1/10 h=1/50

Page 16: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

16

Condition number

● The enrichment may lead to almost-singular matrices difficult to use iterative solvers

Page 17: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

17

Preconditionner

● Orthogonalize each subset of enriched dofs

Page 18: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

18

Preconditionner

● Cholesky decomposition & scaling for node k :

“Assembly” of every submatrix R gives R*

Page 19: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

19

Preconditionner

● Trick for handling non positive definite systems (but blockwise positive definite)– If the matrix A belongs to the electrostatic part:

● a, b and c are negative● we need to take the opposite matrix (which is positive

definite) in order to generate the preconditionner

Page 20: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

20

Preconditionner

Condition number of or

(Geometrical enrichment)

Page 21: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

21

Convergence study● Exact solution – use of complex potentials

– cf. H. Sosa, Plane problems in piezoelectric media with defects, Int. J. Sol. Struct. (1991)

Page 22: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

22

Convergence study

EDx EDy

Sxy

SyyPotential

Sxx

Page 23: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

23

Convergence study

Displacement

Page 24: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

24

Convergence study

● Energy norm with respect to the internal energy

Page 25: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

25

Convergence study● Energy norm

– comparison with standard crack tip enrichment

– Infinite body with embedded crack

– inclined material axes (30°)

– PZT4 orthotropic material

Page 26: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

26

Energy error

● Exact vs finite element fields

Page 27: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

27

Energy error

● The “classical” enrichment gives almost the same results as the specific enrichment, with less computational overhead.

● It is not clear whether different material laws (e.g. “more” anisotropic) lead to different results

Page 28: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

28

SIFs computation

● Contour integrals on Γ −> Domain integrals over V

Page 29: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

29

Interaction integrals

● Same procedure used to compute interaction integrals (no crack loading) :

Page 30: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

30

Interaction integrals

● Relation between G and the K factors– Simpler case of the isotropic elasticity well known

– The Irwin matrix depends on the material orientation and is not explicitely known for piezos.

Page 31: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

31

Interaction integrals

● By using the eigenfunction set, every term in the Irwin matrix can be determined– for instance :

– No need of finite element support because the Irwin matrix is intrinsic (for a given material orientation)

Page 32: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

32

Choice of

The field describes the geometry of the integration domain S.

Two choices of integration domain with regard to h:– Topological

– Geometrical

Page 33: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

33

J error

● Exact J vs computed J-integral

Page 34: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

34

Conclusion & future work

● Application of the X-Fem for piezos● The convergence study shows that the four

classical enrichment functions are enough● Use of equivalement volume integrals to

compute the electromechanical J-integral● Interaction integrals will be used to extract K

factors● Auxiliary fields can also be used to compute

the local Irwin matrix

Page 35: Application of the X-FEM to the fracture of piezoelectric ... · 34 Conclusion & future work Application of the X-Fem for piezos The convergence study shows that the four classical

35

Conclusion & future work

● Systematic investigations of J and the K factors's accuracy

● Propagation laws● Investigation for a electrically permeable

crack● 3D extensions (esp. for the eigenfunctions

needed for K extraction)