kin2.fbl

# run the simulations and generate the images
ulin Kinematic coupling
sys ccx solve2 > solve-k.log
read solve2.frd new
read solve2.inp nom

frame
view elem
rot -z
rot u 90
rot l 45
rot u 45
tra l 30
plot fv all
# movie counter
valu counter 0
# number of modes
valu num 10
# increment for modes
valu incr 1
# loop over movies
while counter < num
  valu counter + counter incr
  # convert counter value to integer
  valu counter int counter
  # build output file name from prefix, counter and extension
  valu target & Refs2/kin counter
  valu target & target .gif
  # make sure that the frequency is correctly displayed
  ds counter e 4
  # write 11 frames from the output of the subsequent command
  movi frames auto
  ds counter a 4
  # store the movie under the target name target
  sys mv movie.gif target
endwhile
ulin Axial tension in part 2
rot x
rot c -90
rot u 20
ds 12 e 7
view disp
scal d 100
frame
plot fv all
hcpy png Refs2/se1-kin
ulin Vertical load at the end of part 2
rot y
rot c 180
ds 14 e 7
scal d
scal d 5
hcpy png Refs2/se2-kin
sys dat2txt.py solve2
sys gnuplot frequency.plt
sys mv f.png Refs2/f-kin.png