Is it possible to apply pressure-load into cracks-geometry?

[oradzhabov]

Here is one of basic examples
https://easyfea.readthedocs.io/en/latest/examples/Meshes/Mesh5_2D.html

Where I can see that after cracks crack1, crack2, crack3, crack4, crack5 have been added to mesh, in result we see opened cracks.

Model has external boundary conditions.
Is it possible to apply pressure-load into these crack-lines?

I want to simulate air-pressure, which applied inside in crack-surface and increases crack-propagation.

I've tried create smth similar like in example above:

crack_lines = []
for poly in cracks_xy:
    pts = [Point(p[0], p[1], isOpen=True) for p in poly]
    crack_lines.append(Points(pts, meshSize=meshSize, isHollow=False, isOpen=True))

...
mesh = out_contour.Mesh_2D([], ElemType.TRI3, cracks=crack_lines, refineGeoms=refine, additionalSurfaces=holes)
...
crack_nodes_arr = [mesh.Nodes_Points(ll.points) for ll in crack_lines]
...
for i in range(len(crack_nodes_arr)):
sim.add_pressureLoad(crack_nodes_arr[i], pressure)
...
but it crashed in "mesh.Get_normals(nodes)".

I've fixed simple bug locally to check further behavior, but when I've checked
Display.Plot_BoundaryConditions(sim)

it showed that my boundary conditions did not applied properly.
It shows crack-lines, but BC covered only part of mesh-vertices.
And I'm not sure that pressure will push in both sides of mesh from the crack-line.

Thanks you.

[oradzhabov]

After some investigation, I found that "crack_lines.append(Points(pts, meshSize=meshSize, isHollow=False, isOpen=True))" not the best,
Better to use

        for i in range(len(poly)-1):
            ll = Line(Point(poly[i][0], poly[i][1]),
                      Point(poly[i+1][0], poly[i+1][1]), meshSize=meshSize, isOpen=False)
            crack_lines.append(ll)

because Line() will find all mesh-points, instead of Points()

[matnoel]

Hi @oradzhabov,

I made a few minor updates to the EasyFEA/FEM/_mesher module. You can install the latest version using the following command:

pip install git+https://github.com/matnoel/EasyFEA.git@dev

To test the changes, you can use this example file: Mesh5_2D.py.

If you are interested in running a 3D simulation, feel free to let me know—I’d be happy to assist!

Best regards,
Matthieu