$ file ortx72   
$ file f108  from KIZs pr9_1_111  Changes in lower case
diag 8,15
SOL 107 $
INCLUDE 'SSSALTERDIR:ortx7a.alt'
CEND
TITLE= RECT. VOLUME ( 1X1X1)- 3 SIDES FIXED 	ortx72
SUBTITLE= FLUID-STRUCTURE INTERACTION,4TH SIDE ON BEAM ON ELAST. FOUND.
LABEL= RESIDUAL ONLY- complex eigenvalues, scaled eigenvector alter
$echo= both
CMETHOD = 10
$PARAM,USETPRT,0
param,fluidmp,-1
param,structmp,-1
param,panelmp,-1
param,grdpnt,0
$
set 99=100 thru 224
$pressure(punch,sort1)=99
set 9=41 thru 65
$disp=9
DISP=ALL
STRESS = ALL
SPC =1313
BEGIN BULK
EIGC,   10,     HESS,   ,       ,       ,       ,       10
spc,	1313,	441,	123456 $ get the ones autospc won't
spc,	1313,	201001,	0
tabdmp1	1000
	0.0	10.0	1000.0	10.0	endt
$
$ MODAL REDUCTION ON BOTH
$
$ PRODUCE OUTPUT FORTRAN FILES
$PARAM,POST,-2
$
$ DYNAMIC LOADING
LSEQ	2000	101	1101
RLOAD1	1100 	101			1001
SLOAD 	1101	120	125.
SLOAD	1101	124	125.
SLOAD	1101	220	125.
SLOAD	1101	224	125.
TABLED4	1001	0.0	1.0	0.0	1.E6
	0.0	0.0	1.0	0.0	ENDT
$ FRQUENCY RANGE
FREQ1	200	0.00001	1.0	4
$ METHOD CARD FOR EIGEN VALUE
PARAM,NEWSEQ,-1
$EIGR	30	GIV				50
$EIGR	40	GIV				5
EIGRL	30   	-0.001		50
EIGRL	20   	-0.001		20
$ 
$THE STRUCTURAL POINTS
GRID   	41	   	0.0	0.0	0.0
GRID   	441	   	0.00	0.0	0.0
GRID	45	  	1.0	0.0	0.0		
$
GRID	61	   	0.0	0.20	0.0		
GRID	65	  	1.0	0.20	0.0		
$
$THE STRUCTURAL ELEMENTS
CQUAD4	100	4444	41	45	65	61
$ STRUCTURAL PROPERTIES
PSHELL	4444	77	.05	77	1.0
MAT1	77	100.	   	.333	1.000
$
$ FLUID POINTS
$ GRIDS 200 THRU 204 AND 100 THRU 104 ARE CONTACT GRIDS. KEEP THEM IN THE
$ RESIDUAL
GRID	200	0       0.0	0.2	0.0	-1
GRID      201001       0  350.60  2931.3  714.33      -1
GRID	204	  	1.0	0.2	.0	-1
GRID	220	  	0.0	0.2	1.0	-1
GRID	224		1.0	0.2	1.0	-1
GRID	100	   	0.0	0.0	0.0	-1
GRID	104	  	1.0	0.0	.0	-1
GRID	120	  	0.0	0.0	1.0	-1
GRID	124		 1.0	0.0	1.0	-1
$
$ FLUID ELEMENTS
CHEXA	1000	15	100	104	124	120	200	204
	224	220
$
$ FLUID PROPERTIES
PSOLID	15	25 				1	PFLUID
MAT10	25    	1.	0.1
$
$ THE ELASTIC SPRINGS OF THE FOUNDATION
$ SPRINGS FROM THE QUAD4'S TO THE GROUND
$ KIZ CORRECTED TO THE VALUE IN THE REFERENCE
$ I.E. .1 IS CHANGED TO .125
$ AT CORNERS HALF SPRINGS (THE K IS/UNIT LENGTH SO AT CORNERS ITS HALF)
CELAS2	1007	.0625	61	3
CELAS2	1008	.0625	65	3
CELAS2	1009	.0625	41	3
CELAS2	1010	.0625	45	3
$
$ SPC ALL THE STRUCTRUAL POINTS NOT TO ROTATE ABOUT X OR Z
SPC1	1313	1246	41	THRU	65
$ SPC THE END STRUCT. POINTS NOT TO ROATE ABOUT Y (STATED IN THE PROBLEM)
SPC1	1313	5	41	61	65	45
$ define structural panels
panel	firstone111
set1	111	41  	45	65	61
set1	222	61	65	41	45
set1	333	41	41	45	45
set1	444	41	45	61	65
set1	200	100	104	200	204
acmodl	ident	all	200	111
$all		111
param, debug, yes $ print eigenvector scaling vector
ENDDATA