[Getfem-users] Crack SIF calculation

[Oaks, Aaron Jameson]

Greetings,

I've got a crack set up in an elastic deformation system as in demo_crack.py and it looks to be giving me the right answer for the given crack, so I'm looking at getting the crack to propagate. I've been reading over the [XFEM] reference by Moes et al. [1] that propagates the crack by calculating stress intensity factors from interaction integrals, but I'm not sure how to implement these calculations in GetFEM. I saw some work on this was done in contrib/crack_plate/crack_bilaplacian_sif.cc, but it is done in the low level assembly (which I'm not yet familiar with), and contains lines like:

    cerr << __PRETTY_FUNCTION__ << " IS TO BE DONE\n";

    /* too lazy to do it now */

So I'm not sure how much of it is correct, or is compatible with the setup in demo_crack.py.

Have any had success performing these interaction integral SIF calculations using some form of this contrib/ code, or any other implementation within GetFEM?

Thank you for your time.

Aaron Oaks

[1] http://dx.doi.org/10.1002/(SICI)1097-0207(19990910)46:1<131::AID-NME726>3.0.CO;2-J