ttwo-grain-collision.sh - granular - granular dynamics simulation
 (HTM) git clone git://src.adamsgaard.dk/granular
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) README
 (DIR) LICENSE
       ---
       ttwo-grain-collision.sh (874B)
       ---
            1 #!/bin/sh
            2 set -e
            3 
            4 id=two-grain-collision
            5 
            6 rm -f ${id}.grains.*.{tsv,png} ${id}.mp4
            7 
            8 (granulargrain -R -u 0.1; granulargrain -f -x 1.2) | granular -e 4.0 -I 0.1 ${id}
            9 
           10 for f in ${id}.grains.*.tsv; do
           11         #granular2img -f '$50' -l 'number of contacts' -t png < "$f" > "${f%.tsv}.png"
           12         #granular2img -f '$5' -l 'x velocity [m/s]' -t png < "$f" > "${f%.tsv}.png"
           13         granular2img -f '$14' -l 'force_x [N]' -t png < "$f" > "${f%.tsv}.png"
           14 done
           15 
           16 ffmpeg -y -framerate 5 -i ${id}.grains.%05d.png \
           17         -c:v libx264 -r 30 -pix_fmt yuv420p ${id}.mp4
           18 xdg-open ${id}.mp4 
           19 
           20 > "${id}.energy.tsv"
           21 for f in ${id}.grains.*.tsv; do
           22         granularenergy < "$f" >> "${id}.energy.tsv"
           23 done
           24 gnuplot -e "set term png;\
           25             set xlabel 'time step';\
           26             set ylabel 'Energy [J]';\
           27                         plot '-' u 0:1 w lp t 'Total energy'" \
           28                          < "${id}.energy.tsv" > "${id}.energy.png"
           29 xdg-open "${id}.energy.png"