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...

1. I take the "Linear Elasticity" example as a starting point (2.5.1 in the manual)

2. I substitute my own mesh and regions, and boundary conditions

3. I solve the problem

4. I look at the tensions between the nodes. If they exceed a certain threshold, the respective connection breaks up. I change the mesh accordingly.

5. 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.