Mesh6_3D

Refined 3D mesh in zones.

Static Scene

Interactive Scene

Static Scene

Interactive Scene

/home/docs/checkouts/readthedocs.org/user_builds/easyfea/envs/v1.5.5/lib/python3.11/site-packages/EasyFEA/utilities/PyVista.py:969: PyVistaDeprecationWarning:
../../../../envs/v1.5.5/lib/python3.11/site-packages/EasyFEA/utilities/PyVista.py:969: Arguments 'pointa', 'pointb', 'center' must be passed as keyword arguments to function 'CircularArc'.
From version 0.50, passing these as positional arguments will result in a TypeError.
  arc1 = pv.CircularArc(geom.pt1.coord, geom.pt3.coord, geom.center.coord)
/home/docs/checkouts/readthedocs.org/user_builds/easyfea/envs/v1.5.5/lib/python3.11/site-packages/EasyFEA/utilities/PyVista.py:970: PyVistaDeprecationWarning:
../../../../envs/v1.5.5/lib/python3.11/site-packages/EasyFEA/utilities/PyVista.py:970: Arguments 'pointa', 'pointb', 'center' must be passed as keyword arguments to function 'CircularArc'.
From version 0.50, passing these as positional arguments will result in a TypeError.
  arc2 = pv.CircularArc(

from EasyFEA import Display, ElemType, PyVista
from EasyFEA.Geoms import Circle, Domain

if __name__ == "__main__":
    Display.Clear()

    L = 1
    meshSize = L / 4

    contour = Domain((0, 0), (L, L), meshSize)
    circle = Circle((L / 2, L / 2), L / 3, meshSize)
    inclusions = [circle]

    refine1 = Domain((0, L), (L, L * 0.8), meshSize / 8)
    refine2 = Circle(circle.center, L / 2, meshSize / 8)
    refine3 = Circle((0, 0), L / 2, meshSize / 8)
    refineGeoms = [refine1, refine2, refine3]

    PyVista.Plot_Geoms([contour, circle, *refineGeoms]).show()

    mesh = contour.Mesh_Extrude(
        inclusions, [0, 0, -L], [5], ElemType.PRISM15, refineGeoms=refineGeoms
    )
    PyVista.Plot_Mesh(mesh).show()