Hi Torsten,
On 09/01/11 17:30, Torsten Bronger wrote:
Hallöchen!
I try to model crack propagation in thin film with sfepy. I must admit that I don't understand everything in the docs because my FEM knowledge is mediocre at best, but my plan is the following:
The docs a lacking, so do not hesitate to ask...
I take the "Linear Elasticity" example as a starting point (2.5.1 in the manual)
I substitute my own mesh and regions, and boundary conditions
I solve the problem
I look at the tensions between the nodes. If they exceed a certain threshold, the respective connection breaks up. I change the mesh accordingly.
I repeat 3 to 4 until no connections have to be broken up anymore
Is this a sensible approach?
I am not expert in this field, but your approach seems feasible.
You can make a function similar to Mesh.localize() for item 4. Mesh.localize() is for removing nodes from a mesh, so it is not directly what you need, but it might be an inspiration.
Cheers, r.