Mortar contact formulation using smooth active set strategy applied to 3d crack propagation

I Athanasiadis, Zahur Ullah, Lukasz Kaczmarczyk, Chris J. Pearce

Research output: Contribution to conferenceOther

Abstract

This work focuses on 3D modelling of multi-body frictionless contact of brittle materials by means of the mortar method under the assumption of small displacements and infinitesimal deformations. The aim is to understand the influence of contact on quasistatic fracture propagation in nuclear graphite bricks. Surfaces of two bodies that are candidate to come to contact are denoted as master and slave surfaces. No self-contact is considered and body domains are discretised by tetrahedral elements. When two triangular faces of tetrahedral elements are candidate to come to contact, with one face belonging to the master and the other to the slave surface, they are used to create a special triangular prism element. Only small displacements are taken into account, hence prisms are created only at the beginning of an analysis. These elements are used to perform integration of displacements and Lagrange multipliers shape functions and mapping of the displacements shape functions from the master to the slave surface according to the mortar method [1]. At each step of the analysis, prism triangular faces can either be in contact or form a gap denoting an active or passive state, respectively. The state of a prism is determined by evaluation of the complementary function proposed by [2] that, in the present work, is modified in order to yield a smooth Newton method algorithm. This modification involves two steps, expression of max function encompassed in the complementary function by means of an absolute value and then substitution of the C0 absolute function with a similar C1 function. Furthermore, the focus of this work is the contact between nuclear graphite bricks, that are also subjected to fracture, where crack propagation is solved in the context of configurational mechanics for infinitesimal deformations [3]. The model will be used to investigate the influence of contact when a slice of a nuclear graphite brick, that has an initial crack, is loaded through graphite keys. The keys are subjected to displacement boundary conditions that transfer tractions to the brick resulting to propagation of the crack. Results will be compared to numerical results where keys and brick form of a continuous body and experimental results presented in [3].

REFERENCES

[1] Ullah, Z., Kaczmarczyk, L, and Pearce, C.J., Three-dimensional mortar contact formulation: an efficient and accurate numerical implementation. In: 25th UKACM Conference on Computational Mechanics, University of Birmingham Birmingham, United Kingdom 4 pp.

[2] Popp, A. and Gee, M.W. and Wall, W.A. O.C. Zienkiewicz, O.C. and R.L. Taylor, A finite deformation mortar contact formulation using a primal-dual active set strategy. Int. J. Numer. Methods Eng., Vol. 79, pp. 1354–1391, 2009.

[3] Kaczmarczyk, L, Ullah, Z. and Pearce, C.J., Energy consistent framework for continuously evolving 3D crack propagation. Comput. Meth. Appl. Mech. Eng., Vol. 324, pp. 54–73, 2017.

Conference

Conference6th European Conference on Computational Mechanics (ECCM 6)
7th European Conference on Computational Fluid Dynamics (ECFD 7)
Abbreviated titleECCM 6 and ECFD 7
CountryUnited Kingdom
CityGlasgow
Period11/06/1815/06/18
Internet address

Fingerprint

Mortar
Crack propagation
Brick
Prisms
Graphite
Computational mechanics
Cracks
Lagrange multipliers
Newton-Raphson method
Brittleness
Mechanics
Substitution reactions
Boundary conditions

Keywords

  • Contact
  • Mortar Method
  • Active Set
  • Small Strains
  • Small Displacements
  • Brittle Fracture

Cite this

Athanasiadis, I., Ullah, Z., Kaczmarczyk, L., & Pearce, C. J. (2018). Mortar contact formulation using smooth active set strategy applied to 3d crack propagation. 6th European Conference on Computational Mechanics (ECCM 6)
7th European Conference on Computational Fluid Dynamics (ECFD 7) , Glasgow, United Kingdom.
Athanasiadis, I ; Ullah, Zahur ; Kaczmarczyk, Lukasz ; Pearce, Chris J. / Mortar contact formulation using smooth active set strategy applied to 3d crack propagation. 6th European Conference on Computational Mechanics (ECCM 6)
7th European Conference on Computational Fluid Dynamics (ECFD 7) , Glasgow, United Kingdom.
@conference{6d21b540c43c49a4b094f69b42715d27,
title = "Mortar contact formulation using smooth active set strategy applied to 3d crack propagation",
abstract = "This work focuses on 3D modelling of multi-body frictionless contact of brittle materials by means of the mortar method under the assumption of small displacements and infinitesimal deformations. The aim is to understand the influence of contact on quasistatic fracture propagation in nuclear graphite bricks. Surfaces of two bodies that are candidate to come to contact are denoted as master and slave surfaces. No self-contact is considered and body domains are discretised by tetrahedral elements. When two triangular faces of tetrahedral elements are candidate to come to contact, with one face belonging to the master and the other to the slave surface, they are used to create a special triangular prism element. Only small displacements are taken into account, hence prisms are created only at the beginning of an analysis. These elements are used to perform integration of displacements and Lagrange multipliers shape functions and mapping of the displacements shape functions from the master to the slave surface according to the mortar method [1]. At each step of the analysis, prism triangular faces can either be in contact or form a gap denoting an active or passive state, respectively. The state of a prism is determined by evaluation of the complementary function proposed by [2] that, in the present work, is modified in order to yield a smooth Newton method algorithm. This modification involves two steps, expression of max function encompassed in the complementary function by means of an absolute value and then substitution of the C0 absolute function with a similar C1 function. Furthermore, the focus of this work is the contact between nuclear graphite bricks, that are also subjected to fracture, where crack propagation is solved in the context of configurational mechanics for infinitesimal deformations [3]. The model will be used to investigate the influence of contact when a slice of a nuclear graphite brick, that has an initial crack, is loaded through graphite keys. The keys are subjected to displacement boundary conditions that transfer tractions to the brick resulting to propagation of the crack. Results will be compared to numerical results where keys and brick form of a continuous body and experimental results presented in [3].REFERENCES [1] Ullah, Z., Kaczmarczyk, L, and Pearce, C.J., Three-dimensional mortar contact formulation: an efficient and accurate numerical implementation. In: 25th UKACM Conference on Computational Mechanics, University of Birmingham Birmingham, United Kingdom 4 pp. [2] Popp, A. and Gee, M.W. and Wall, W.A. O.C. Zienkiewicz, O.C. and R.L. Taylor, A finite deformation mortar contact formulation using a primal-dual active set strategy. Int. J. Numer. Methods Eng., Vol. 79, pp. 1354–1391, 2009. [3] Kaczmarczyk, L, Ullah, Z. and Pearce, C.J., Energy consistent framework for continuously evolving 3D crack propagation. Comput. Meth. Appl. Mech. Eng., Vol. 324, pp. 54–73, 2017.",
keywords = "Contact, Mortar Method, Active Set, Small Strains, Small Displacements, Brittle Fracture",
author = "I Athanasiadis and Zahur Ullah and Lukasz Kaczmarczyk and Pearce, {Chris J.}",
year = "2018",
month = "6",
day = "11",
language = "English",
note = "6th European Conference on Computational Mechanics (ECCM 6)<br/>7th European Conference on Computational Fluid Dynamics (ECFD 7) , ECCM 6 and ECFD 7 ; Conference date: 11-06-2018 Through 15-06-2018",
url = "http://www.eccm-ecfd2018.org/frontal/default.asp",

}

Athanasiadis, I, Ullah, Z, Kaczmarczyk, L & Pearce, CJ 2018, 'Mortar contact formulation using smooth active set strategy applied to 3d crack propagation' 6th European Conference on Computational Mechanics (ECCM 6)
7th European Conference on Computational Fluid Dynamics (ECFD 7) , Glasgow, United Kingdom, 11/06/18 - 15/06/18, .

Mortar contact formulation using smooth active set strategy applied to 3d crack propagation. / Athanasiadis, I; Ullah, Zahur; Kaczmarczyk, Lukasz; Pearce, Chris J.

2018. 6th European Conference on Computational Mechanics (ECCM 6)
7th European Conference on Computational Fluid Dynamics (ECFD 7) , Glasgow, United Kingdom.

Research output: Contribution to conferenceOther

TY - CONF

T1 - Mortar contact formulation using smooth active set strategy applied to 3d crack propagation

AU - Athanasiadis, I

AU - Ullah, Zahur

AU - Kaczmarczyk, Lukasz

AU - Pearce, Chris J.

PY - 2018/6/11

Y1 - 2018/6/11

N2 - This work focuses on 3D modelling of multi-body frictionless contact of brittle materials by means of the mortar method under the assumption of small displacements and infinitesimal deformations. The aim is to understand the influence of contact on quasistatic fracture propagation in nuclear graphite bricks. Surfaces of two bodies that are candidate to come to contact are denoted as master and slave surfaces. No self-contact is considered and body domains are discretised by tetrahedral elements. When two triangular faces of tetrahedral elements are candidate to come to contact, with one face belonging to the master and the other to the slave surface, they are used to create a special triangular prism element. Only small displacements are taken into account, hence prisms are created only at the beginning of an analysis. These elements are used to perform integration of displacements and Lagrange multipliers shape functions and mapping of the displacements shape functions from the master to the slave surface according to the mortar method [1]. At each step of the analysis, prism triangular faces can either be in contact or form a gap denoting an active or passive state, respectively. The state of a prism is determined by evaluation of the complementary function proposed by [2] that, in the present work, is modified in order to yield a smooth Newton method algorithm. This modification involves two steps, expression of max function encompassed in the complementary function by means of an absolute value and then substitution of the C0 absolute function with a similar C1 function. Furthermore, the focus of this work is the contact between nuclear graphite bricks, that are also subjected to fracture, where crack propagation is solved in the context of configurational mechanics for infinitesimal deformations [3]. The model will be used to investigate the influence of contact when a slice of a nuclear graphite brick, that has an initial crack, is loaded through graphite keys. The keys are subjected to displacement boundary conditions that transfer tractions to the brick resulting to propagation of the crack. Results will be compared to numerical results where keys and brick form of a continuous body and experimental results presented in [3].REFERENCES [1] Ullah, Z., Kaczmarczyk, L, and Pearce, C.J., Three-dimensional mortar contact formulation: an efficient and accurate numerical implementation. In: 25th UKACM Conference on Computational Mechanics, University of Birmingham Birmingham, United Kingdom 4 pp. [2] Popp, A. and Gee, M.W. and Wall, W.A. O.C. Zienkiewicz, O.C. and R.L. Taylor, A finite deformation mortar contact formulation using a primal-dual active set strategy. Int. J. Numer. Methods Eng., Vol. 79, pp. 1354–1391, 2009. [3] Kaczmarczyk, L, Ullah, Z. and Pearce, C.J., Energy consistent framework for continuously evolving 3D crack propagation. Comput. Meth. Appl. Mech. Eng., Vol. 324, pp. 54–73, 2017.

AB - This work focuses on 3D modelling of multi-body frictionless contact of brittle materials by means of the mortar method under the assumption of small displacements and infinitesimal deformations. The aim is to understand the influence of contact on quasistatic fracture propagation in nuclear graphite bricks. Surfaces of two bodies that are candidate to come to contact are denoted as master and slave surfaces. No self-contact is considered and body domains are discretised by tetrahedral elements. When two triangular faces of tetrahedral elements are candidate to come to contact, with one face belonging to the master and the other to the slave surface, they are used to create a special triangular prism element. Only small displacements are taken into account, hence prisms are created only at the beginning of an analysis. These elements are used to perform integration of displacements and Lagrange multipliers shape functions and mapping of the displacements shape functions from the master to the slave surface according to the mortar method [1]. At each step of the analysis, prism triangular faces can either be in contact or form a gap denoting an active or passive state, respectively. The state of a prism is determined by evaluation of the complementary function proposed by [2] that, in the present work, is modified in order to yield a smooth Newton method algorithm. This modification involves two steps, expression of max function encompassed in the complementary function by means of an absolute value and then substitution of the C0 absolute function with a similar C1 function. Furthermore, the focus of this work is the contact between nuclear graphite bricks, that are also subjected to fracture, where crack propagation is solved in the context of configurational mechanics for infinitesimal deformations [3]. The model will be used to investigate the influence of contact when a slice of a nuclear graphite brick, that has an initial crack, is loaded through graphite keys. The keys are subjected to displacement boundary conditions that transfer tractions to the brick resulting to propagation of the crack. Results will be compared to numerical results where keys and brick form of a continuous body and experimental results presented in [3].REFERENCES [1] Ullah, Z., Kaczmarczyk, L, and Pearce, C.J., Three-dimensional mortar contact formulation: an efficient and accurate numerical implementation. In: 25th UKACM Conference on Computational Mechanics, University of Birmingham Birmingham, United Kingdom 4 pp. [2] Popp, A. and Gee, M.W. and Wall, W.A. O.C. Zienkiewicz, O.C. and R.L. Taylor, A finite deformation mortar contact formulation using a primal-dual active set strategy. Int. J. Numer. Methods Eng., Vol. 79, pp. 1354–1391, 2009. [3] Kaczmarczyk, L, Ullah, Z. and Pearce, C.J., Energy consistent framework for continuously evolving 3D crack propagation. Comput. Meth. Appl. Mech. Eng., Vol. 324, pp. 54–73, 2017.

KW - Contact

KW - Mortar Method

KW - Active Set

KW - Small Strains

KW - Small Displacements

KW - Brittle Fracture

M3 - Other

ER -

Athanasiadis I, Ullah Z, Kaczmarczyk L, Pearce CJ. Mortar contact formulation using smooth active set strategy applied to 3d crack propagation. 2018. 6th European Conference on Computational Mechanics (ECCM 6)
7th European Conference on Computational Fluid Dynamics (ECFD 7) , Glasgow, United Kingdom.