Note
Go to the end to download the full example code.
Mesh6_3D#
Refined 3D mesh in zones.


13 from EasyFEA import Display, ElemType, PyVista
14 from EasyFEA.Geoms import Circle, Domain
15
16 if __name__ == "__main__":
17 Display.Clear()
18
19 L = 1
20 meshSize = L / 4
21
22 contour = Domain((0, 0), (L, L), meshSize)
23 circle = Circle((L / 2, L / 2), L / 3, meshSize)
24 inclusions = [circle]
25
26 refine1 = Domain((0, L), (L, L * 0.8), meshSize / 8)
27 refine2 = Circle(circle.center, L / 2, meshSize / 8)
28 refine3 = Circle((0, 0), L / 2, meshSize / 8)
29 refineGeoms = [refine1, refine2, refine3]
30
31 PyVista.Plot_Geoms([contour, circle, *refineGeoms]).show()
32
33 mesh = contour.Mesh_Extrude(
34 inclusions, [0, 0, -L], [5], ElemType.PRISM15, refineGeoms=refineGeoms
35 )
36 PyVista.Plot_Mesh(mesh).show()
Total running time of the script: (0 minutes 1.001 seconds)