$
$*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$ ----------- Periodic Pendulum Motion ---------------
$
$*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$*@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
$
$*********************************     set basic parameters
$
&dimen -dimen feet kips
&device  -pri device
&title Sample Of a Simple Pendulum
$
$*********************************     read in the model
$
inmodel
$
$********************************     set the condition
$
&instate -condition test 0 0 90
$
$********************************     edit the model
$
medit
   ~pin  fix
   connector a ~pin    *two
end
$
$*******************************    time for simulation
$
&env test  -time 20 .025
$
$*******************************    time domain
$*******************************    newmark method
$
TDOM -newm  1/4 1/2 -no_cap
$
$*******************************    post processing
$
prcpost
  &describe interest -assoc @
  illustrate
$     pic -events 0 5.5 0.1
  end
  traj
    report loca velo -even 0 800 .10
   end
  traject
    report loca velo -even 0 800 .5
    &channel gra_dev -file pend.ans/pend_mot.eps
    p 1 9  -t_left 'Roll Angle (Deg)' -t_x 'Time (Seconds)'
    &channel gra_dev -file pend.ans/pend_ace.eps
    p 1 21 -t_left 'Roll Angular Acceleration (Deg/Sec^2)' \
           -t_x 'Time (Seconds)'
    end
 conforce -even 0 800 .25
    rep
    end
end
$
$******************************     all done
$
&fini