options.dparam[sn.SICONOS_DPARAM_TOL] = 1e-3
with MechanicsHdf5Runner(mode='w', io_filename=fn) as io:
ch = chute.create_chute(io,
box_height=box_height,
box_length=box_length,
box_width=box_width,
plane_thickness=plane_thickness,
scale=1,
trans=[-0.9, -1.8, -1])
rcs = rocas.create_rocas(io,
n_layer=n_layer,
n_row=n_row,
n_col=n_col,
x_shift=1.8,
roca_size=0.1,
top=3,
rate=0.2,
density=density)
io.add_Newton_impact_friction_nsl('contact', mu=1.0, e=0.01)
run_options = MechanicsHdf5Runner_run_options()
run_options['t0'] = 0
run_options['T'] = T
run_options['h'] = hstep
run_options['theta'] = 1.0
bullet_options = SiconosBulletOptions()
bullet_options.perturbationIterations = 0.
collision_group1=2,
collision_group2=0)
ch = chute_con_vibradores.create_chute(io,
box_height=box_height,
box_length=box_length,
box_width=box_width,
plane_thickness=plane_thickness,
scale=1,
trans=[-0.6, -1.8, -1])
rcs = rocas.create_rocas(io,
n_layer=20,
n_row=2,
n_col=5,
x_shift=1.5,
roca_size=0.13,
top=3,
rate=0.2,
density=density,
distribution=(dist, 0.12))
hstep = 1e-3
with MechanicsHdf5Runner(mode='r+', io_filename=fn) as io:
io.run(
with_timer=False,
t0=0,
T=45, #100*hstep,
h=hstep,
multipoints_iterations=True,
theta=0.50001,
Newton_max_iter=1,
cube_size = 0.1
with MechanicsHdf5Runner(mode='w', io_filename=fn) as io:
ch = chute.create_chute(io,
box_height=box_height,
box_length=box_length,
box_width=box_width,
plane_thickness=plane_thickness,
scale=1,
trans=[-0.6, -1.8, -1])
rcs = rocas.create_rocas(io,
n_layer=200,
n_row=2,
n_col=16,
x_shift=2.0,
roca_size=0.1,
top=3,
rate=0.2,
density=density)
io.add_Newton_impact_friction_nsl('contact', mu=1.0, e=0.01)
with MechanicsHdf5Runner(mode='r+', io_filename=fn) as io:
io.run(t0=0,
T=4,
h=1e-4,
multipoints_iterations=True,
theta=1.0,
Newton_max_iter=1,
solver=Numerics.SICONOS_FRICTION_3D_NSGS,