$ Simple Example to illustrate use of ridgyro alter
$ Free Vibrations with Rotation
$
$ =============================================================
$
sol 107
diag 8
include 'ridgyroa.v2001'
cend
$
title=ridgyro Example
subtitle=Free Vibrations with Rotation
echo=sort
$
disp=all
cmethod=10
$
begin bulk
$
param, grdpnt, 0
$
eigc, 10, hess,,,,, 10
$
$ Reference: 500 Rotations/Minute
param, rpm, 500.
$
dti, rgyro, 0
$          Rec #, Grid, Iaa , Axis, Factor
dti, rgyro, 1   ,    3, 200.,    1,    1.
$ 
$ Remarks:
$ - Axis: 1=x-axis, 2=y-axis, 3=z-axis
$ - Iaa is the moment of inertia with respect to the rotation
$   axis.
$ - The speed of rotations is Factor*RPM, given in rotations per
$   minute.
$
$include 'model.bdf'
$
$ The model consists of a shaft of length 1m which is rigidly 
$ supported at both ends. In the middle of the shaft there
$ is a rigid rotor.
$
$
$         z ^   
$           |           III Rotor
$                       III/
$           |           III           |
$           |-------------------------|  --> x
$           |1     2    I3I    4     5|
$                       III
$                       III
$
$ Problem Data:
$
$    Shaft:   Circular Cross Section with a radius of 0.01m
$             Material is steel with a Young's modulus
$             of 2.E11N/m and a mass density of 7850kg/m^3.
$
$    Rotor:   Mass:  m  =10kg
$             Polar Moment of Inertia: 
$                    Ixx=200kgm^2
$             Diametral Moments of Inertia:
$                    Iyy=Izz=100kgm^2
$
$ Units: N, kg, m
$
$ =============================================================
$
$ Data of the Shaft
$
grid, 1,,  0., 0., 0.,, 123456
=,   *4,, *1.,==
grid, 2,, 0.25
=,   *1,, *0.25,==
=(1)
$
cbeam, 1, 1, 1, 2, 0., 1.
=,    *1, =, *1, *1, ==
=(2)
$
pbeaml, 1, 1,, rod
      , 0.01
mat1, 1, 2.E11,, 0.1, 7850.
$
$ Data of the Rotor
$
conm2, 100, 3,, 10.,
     , 200.,, 100.,,, 100.
enddata