single_cubic_bv.py

This is an example of how to use the netdem library.

import netdem
 
# units: mm, 1000 kg, N, s
sim = netdem.Simulation()
sim.domain_manager.SetBound(-12, -12, -12, 12, 12, 12)
sim.domain_manager.SetCellSpacing(1.0, 1.0, 1.0)
# sim.dem_solver.contact_solver_factory.settings.gjk_use_erosion = True
 
cnt_model = netdem.LinearSpring(1.4e4, 0.7e4, 0.3, 0.0)
cnt_model_prt = sim.scene.InsertContactModel(cnt_model)
sim.scene.SetNumberOfMaterials(1)
sim.scene.SetCollisionModel(0, 0, cnt_model_prt)
 
bond_model = netdem.ParallelBond(2.8e6, 1.4e6, 3.5e4, 1.75e4, 0.01)
bond_model.label = "bond_model"
bond_model_ptr = sim.scene.InsertContactModel(bond_model)
 
stl_model = netdem.STLReader.ReadFile("data/box.stl")
stl_model.RemoveDuplicateVertices()
stl_model.RemoveUnreferencedVertices()
stl_model.ReorientFacets()
stl_model.SetSize(6.2035)
 
bonded_voronois = netdem.BondedVoronois()
bonded_voronois.SetBondModel(bond_model_ptr)
bonded_voronois.InitFromSTL(stl_model, 1000)
bonded_voronois.SaveAsVTK("tmp/out/shape/bonded_voronois_box.vtk")
sim.scene.InsertParticle(bonded_voronois)
 
sim.scene.SetGravity(0, 0, -9.81e-3)
for p in sim.scene.particle_list:
    p.SetDensity(2650)
 
wall_box = netdem.WallBoxPlane(10, 10, 5.0, 0, 0, 0)
wall_box.ImportToScene(sim.scene)
 
# grav = netdem.Gravity()
# grav.Init(sim)
# sim.modifier_manager.Insert(grav)
# sim.modifier_manager.Enable(grav.label)
 
data_dumper = netdem.DataDumper()
data_dumper.Init(sim)
data_dumper.SetRootPath("tmp/out/")
data_dumper.SetSaveByCycles(100)
data_dumper.dump_contact_info = True
data_dumper.dump_wall_info = True
data_dumper.SaveShapeInfoAsSTL()
sim.modifier_manager.Insert(data_dumper)
sim.modifier_manager.Enable(data_dumper.label)
 
# add deformation drived loading for top wall
disp_control = netdem.WallMotionControl()
disp_control.Init(sim)
disp_control.SetWall([5])
disp_control.SetVelocity(0, 0, -0.25)
sim.modifier_manager.Insert(disp_control)
sim.modifier_manager.Enable(disp_control.label)
 
netdem.omp_set_num_threads(1)
sim.Run(2.0)