GopherProxy

	ttwist.c - slidergrid - grid of elastic sliders on a frictional surface
	git clone git://src.adamsgaard.dk/slidergrid
	Log
	Files
	Refs
	README
	LICENSE
	---
	ttwist.c (1264B)
	---
	1 #include "../../slidergrid/simulation.h"
	2 #include "../../slidergrid/grid.h"
	3 #include "../../slidergrid/slider.h"
	4
	5 #include <stdio.h>
	6
	7 // test a regular, 2d, orthogonal grid of sliders
	8 simulation setup_simulation()
	9 {
	10 // create empty simulation structure with default values
	11 simulation sim = create_simulation();
	12 sim.id = "twist";
	13
	14 // initialize grid of sliders
	15 //int nx = 20;
	16 int nx = 2;
	17 int ny = 1;
	18 int nz = 1;
	19 sim.N = nxnynz;
	20 sim.sliders = create_regular_slider_grid(nx, ny, nz, 1.0, 1.0, 1.0);
	21
	22 sim.bond_length_limit = 1.5;
	23
	24 // set slider masses and moments of inertia
	25 int i;
	26 for (i=0; i<sim.N; i++) {
	27
	28 // set default values
	29 initialize_slider_values(&sim.sliders[i]);
	30
	31 // set custom values for certain parameters
	32 sim.sliders[i].mass = 10.0;
	33 sim.sliders[i].moment_of_inertia = 1.0e3;
	34 sim.sliders[i].bond_parallel_kv_stiffness = 1.0e4;
	35 sim.sliders[i].bond_shear_kv_stiffness = 1.0e5;
	36 sim.sliders[i].bond_twist_kv_stiffness = 1.0e5;
	37 }
	38
	39 sim.sliders[0].angvel.x = 1.0e1;
	40
	41 // set temporal parameters
	42 //sim.time_end = 100.0;
	43 //sim.file_interval = 1.0;
	44 sim.time_end = 1.0;
	45 sim.file_interval = 1.0e-3;
	46
	47 return sim;
	48 }