|
From: | Oaks, Aaron Jameson |
Subject: | [Getfem-users] Crack SIF calculation |
Date: | Wed, 29 Jul 2015 20:53:57 +0000 |
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
|
[Prev in Thread] | Current Thread | [Next in Thread] |