tshear.c - slidergrid - grid of elastic sliders on a frictional surface
 (HTM) git clone git://src.adamsgaard.dk/slidergrid
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) README
 (DIR) LICENSE
       ---
       tshear.c (1261B)
       ---
            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 = "shear";
           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 = nx*ny*nz;
           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 = 1.0;
           33         sim.sliders[i].moment_of_inertia = 1.0;
           34         sim.sliders[i].bond_parallel_kv_stiffness = 1.0e5;
           35         //sim.sliders[i].bond_parallel_kv_viscosity = 1.0e2;
           36         sim.sliders[i].bond_shear_kv_stiffness = 1.0e5;
           37     }
           38 
           39     sim.sliders[0].vel.y = 0.1;
           40 
           41     // set temporal parameters
           42     sim.time_end = 200.0;
           43     sim.file_interval = 0.1;
           44     //sim.time_end = 1.0;
           45     //sim.file_interval = 1.0e-3;
           46 
           47     return sim;
           48 }