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
_______________________________________________
Getfem-users mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-users
