$ file ir21  inertia relief by vecplot rigid body modes 
diag 8
SOL 101
Include 'ir2a.v707' $ use add instead of merge for lmt coeff.
CEND
TITLE = inertia relief with vecplot r.b. modes constraints		ir21
Subtitle = all basic coords, 6 rb modes, no s.e.s 
Echo = both
DISPL=ALL
Oload = all
Elforce = all
Gpfo = all
Spc = 1 $
Subcase 1
Load = 1
Subcase 2
Load = 2
Spcf = all
BEGIN BULK
$ new data for alter
Grid, 99999999 $ make it last in sequence
Spcoff,	99999999,	123456  $ protect from autospc 
$ end of alter-required data
Param, tiny, 0.0 $ print small spcforces
$ data required for old method
param, grdpnt, 3
param, inrel, -2 $ request autosup using vecplot 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
PBAR,   12      21      1.0     2.0     1.0     1.25
MAT1,   21      1.0+4,	,	.3 $ tune so first mode about 5 hz.	
$	eid	gid	cord $ tied to cid
spc,	1,	1,	1 $ remove fore and aft mode
$	gid	cp	x	y	z	cd
GRID,	1,       ,       ,       ,       ,	
GRID,	2,       ,       0.5,	,	,		
GRID,	3,       ,       1.,	,	,	  
Grid,	4,	,	1.5,	,	,	 
Grid,	5,	,	2.0,	,	,	  
CONM2,  1    	1,	,	1.0
CONM2,  2    	2,	,	1.0
CONM2,  3    	3,	,	1.0 
CONM2,  4,	4,	,	1.0 
CONM2,  5    	5,	,	1.0
param, usetprt, 0 $ review sets
force,	1,	1,	,	1.	0.	0.,	1.
force,	2,	1,	,	1.	0.	1.
force,	2,	3,	,	1.	0.	-3.
force,	2,	2,	,	1.	0.	-2.
force,	2,	4,	,	1.	0.	4.
CONM2,  14,    	1,		
,	2.0 $ define torsion modes
CONM2,  24,    	2,		
,	2.0
CONM2,  340,    	3,		
,	2.0
CONM2,  44,    	4,		
,	2.0
CONM2,  54,    	5,		
,	2.0
conm2,	356,	3,	 $ make order r1, r2, r3
,	,	,	10.0,	,	,	12.0
ENDDATA