A higher order scheme for twodimensional quasistatic crack growth simulations
(2007) In Computer Methods in Applied Mechanics and Engineering 196(2124). p.25272538 Abstract
 An efficient scheme for the simulation of quasistatic crack growth in twodimensional linearly elastic isotropic specimens is presented. The crack growth is simulated in a stepwise manner where an extension to the already existing crack is added in each step. In a local coordinate system each such extension is represented as a polynomial of some, user specified, degree, n. The coefficients of the polynomial describing an extension are found by requiring that the mode II stress intensity factor is equal to zero at certain points of the extension. If a crack grows from a, preexisting crack so that a kink develops, the leading term describing the crack shape close to the kink will, in a local coordinate system, be proportional to x(3/2). We... (More)
 An efficient scheme for the simulation of quasistatic crack growth in twodimensional linearly elastic isotropic specimens is presented. The crack growth is simulated in a stepwise manner where an extension to the already existing crack is added in each step. In a local coordinate system each such extension is represented as a polynomial of some, user specified, degree, n. The coefficients of the polynomial describing an extension are found by requiring that the mode II stress intensity factor is equal to zero at certain points of the extension. If a crack grows from a, preexisting crack so that a kink develops, the leading term describing the crack shape close to the kink will, in a local coordinate system, be proportional to x(3/2). We therefore allow the crack extensions to contain such a term in addition to the monomial terms. The discontinuity in the crack growth direction at a kink, the kink angle, is determined by requiring that the mode II stress intensity factor should be equal to zero for an infinitesimal extension of the existing crack. To implement the scheme, accurate values of the stress intensity factors and Tstress are needed in each step of the simulation. These fracture parameters are computed using a previously developed integral equation of the second kind. (c) 2007 Elsevier B.V. All rights reserved. (Less)
Please use this url to cite or link to this publication:
https://lup.lub.lu.se/record/629682
 author
 Englund, Jonas ^{LU}
 organization
 publishing date
 2007
 type
 Contribution to journal
 publication status
 published
 subject
 keywords
 fast, integral equation, multipole method, stress intensity factor, crack growth
 in
 Computer Methods in Applied Mechanics and Engineering
 volume
 196
 issue
 2124
 pages
 2527  2538
 publisher
 Elsevier
 external identifiers

 wos:000246126700017
 scopus:33947634583
 ISSN
 00457825
 DOI
 10.1016/j.cma.2007.01.007
 language
 English
 LU publication?
 yes
 additional info
 The information about affiliations in this record was updated in December 2015. The record was previously connected to the following departments: Numerical Analysis (011015004)
 id
 2244dc654b9844feac074ab127de9da9 (old id 629682)
 alternative location
 http://www.maths.lth.se/na/staff/helsing/paperA.pdf
 date added to LUP
 20160401 16:37:31
 date last changed
 20200311 05:38:40
@article{2244dc654b9844feac074ab127de9da9, abstract = {An efficient scheme for the simulation of quasistatic crack growth in twodimensional linearly elastic isotropic specimens is presented. The crack growth is simulated in a stepwise manner where an extension to the already existing crack is added in each step. In a local coordinate system each such extension is represented as a polynomial of some, user specified, degree, n. The coefficients of the polynomial describing an extension are found by requiring that the mode II stress intensity factor is equal to zero at certain points of the extension. If a crack grows from a, preexisting crack so that a kink develops, the leading term describing the crack shape close to the kink will, in a local coordinate system, be proportional to x(3/2). We therefore allow the crack extensions to contain such a term in addition to the monomial terms. The discontinuity in the crack growth direction at a kink, the kink angle, is determined by requiring that the mode II stress intensity factor should be equal to zero for an infinitesimal extension of the existing crack. To implement the scheme, accurate values of the stress intensity factors and Tstress are needed in each step of the simulation. These fracture parameters are computed using a previously developed integral equation of the second kind. (c) 2007 Elsevier B.V. All rights reserved.}, author = {Englund, Jonas}, issn = {00457825}, language = {eng}, number = {2124}, pages = {25272538}, publisher = {Elsevier}, series = {Computer Methods in Applied Mechanics and Engineering}, title = {A higher order scheme for twodimensional quasistatic crack growth simulations}, url = {https://lup.lub.lu.se/search/ws/files/4728289/4254535.pdf}, doi = {10.1016/j.cma.2007.01.007}, volume = {196}, year = {2007}, }