$ file rnormal1  use suport entry for lanczos rb modes 
$ mag  January 13, 2000
diag 8
SOL 111
Include 'rnormala.v2001' $
CEND
TITLE = use suport entry (r-set) when normalizing rb modes
Subtitle = modal frequency response example
Echo = both
Svect=all $ print a-set eigenvectors
Sdisp(sort2)= all $ print modal solutions
Dload = 1
Sdamp=1 $
Freq=1 $
acce(sort2,phase)=ALL
Method = 1
BEGIN BULK
$ required  data-stucture and modes
$.......2.......3.......4.......5.......6.......7.......8.......9.......10.....
$	2	3	4	5	6	7	8	9
CBAR,   12,      12,     1,      2,      ,       ,       1.0   $ I1 PLANE IS 1-3
CBAR    23       12      2       3                       1.0
Cbar,	34,	12,	3,	4,	,	,	1.0
Cbar,	45,	12,	4,	5,	,	,	1.0 $ straight shaft
PBAR,   12      21      1.0     2.0     1.0     1.25
MAT1,   21      1.0+4,	,	.3 $ tune so first mode about 30 hz.	
$	gid	cp	x	y	z	
GRID,	1,       ,       ,       ,       ,	
GRID,	2,       ,       0.5,	,	,		
GRID,	3,       ,       1.,	,	,	 
Grid,	4,	,	1.5,	,	,	
Grid,	5,	,	2.0,	,	,	 
Suport,	3	123456
CONM2,  1    	1,	,	0.2
CONM2,  2    	2,	,	0.2
CONM2,  3    	3,	,	0.2
CONM2,  4,	4,	,	0.2
CONM2,  5    	5,	,	0.2 $ tune for unit inertia along all axes
conm2,	33,	3,	,	0.0
,	1.0, 	,	0.5,	,	,	0.5 
$ put mass on torsion mode, bring other I's up to unity
eigrl,	1,	,	,	12 $ 6 rigid plus 6 flexible
$ required data, frequency reponse
FREQ1,   1     	10.      10.     10 $ from 10 by 10 for 10
TABDMP1, 1       G
,        0.      .03     1.0      0.03	ENDT $ constant damping
RLOAD1,  1,	5,	,	,	1
Darea,	5,	3,	2,	1.0+4 $ lateral force at center of shaft
TABLED1  1
,         .0      1.      1.      1.      ENDT $ constant table
$optional data
param, grdpnt, 3 $ the center point.  Output weight summary here
param, usetprt, 0 $ review sets
ENDDATA