Visualization of triply periodic minimal surfaces (TPMS)#

importing any supporting modules#

[1]:
import numpy as np
import pyvista as pv
import ipywidgets as widgets
import damask

pv.set_jupyter_backend('pythreejs')

defining variable aspects of the resulting geometry#

[2]:
surface = widgets.Dropdown(
    options=list(damask.Grid._minimal_surface),
    description='Minimal Surface:',
    width=400,
    disabled=False,
    )
grid = widgets.IntSlider(
    value=64,
    min=16,
    max=256,
    step=2,
    description='Grid:',
    disabled=False,
    continuous_update=False,
    orientation='horizontal',
    readout=True,
    readout_format='d'
    )
periods = widgets.IntSlider(
    value=1,
    min=1,
    max=4,
    step=1,
    description='Periods:',
    disabled=False,
    continuous_update=False,
    orientation='horizontal',
    readout=True,
    readout_format='d'
    )
threshold = widgets.FloatSlider(
    value=0.0,
    min=-1,
    max=1,
    step=0.05,
    description='Threshold:',
    disabled=False,
    continuous_update=False,
    orientation='horizontal',
    readout=True,
    readout_format='.2f'
    )

problem setup#

Select which triply periodic minimal surface to visualize at what grid resolution and with how many periods in the filled volume

[3]:
display(surface,grid,periods,threshold)

generate resulting TPMS#

[4]:
tpms = damask.Grid.from_minimal_surface([grid.value]*3,np.ones(3),
                                        surface.value,
                                        threshold=threshold.value,
                                        periods=periods.value)\
                  .vicinity_offset()
[5]:
# transform into pyvista mesh
mesh = pv.UniformGrid(dims=tpms.cells+1,
                      spacing=tpms.size/tpms.cells)
mesh['material'] = tpms.material.flatten(order='F')

visualize TPMS#

[6]:
# generate an interactive visualization of the interface
pl = pv.Plotter()
pl.set_background('white')
pl.add_mesh_threshold(mesh)
pl.show()

../../_images/documentation_examples_triply_periodic_minimal_surfaces_11_0.png
[7]:
# export VTK dataset
mesh.save(f'TPMS_{surface.value}_{threshold.value}_{grid.value}.vti')