*HEADING
 KOYNA DAM: SEISMIC ANALYSIS
 (without hydrodynamic interactions)
 Units - n, m, sec
*************************************************
** Step 1: Gravity load
** Step 2: Hydrostatic pressure load
** Step 3: Earthquake
**         (without hydrodynamic interactions)
**
** Requires amplitude curves koyna_haccel.inp
** and koyna_vaccel.inp
*************************************************
*PREPRINT, MODEL=YES
*NODE
*NGEN, NSET=NBASE
*NGEN, NSET=NMID
*NGEN, NSET=NWL
*NGEN, NSET=NTOP
*NFILL, BIAS=1.05
 NBASE, NMID, 24, 100
*NFILL, BIAS=0.92
 NMID, NWL, 12, 100
*NFILL
 NWL, NTOP, 2, 100
*NSET, NSET=NOUT
*ELSET, ELSET=EOUT
*ELEMENT, TYPE=CPS4R
*ELGEN, ELSET=DAM
*ELSET, ELSET=WDAM, GENERATE
*SOLID SECTION, ELSET=DAM, MATERIAL=CONCRETE
 1.0,
*MATERIAL, NAME=CONCRETE
*ELASTIC
 3.1027E+10, 0.2
*DENSITY
 2643.0
*DAMPING, BETA=0.00323
*CONCRETE DAMAGED PLASTICITY
*CONCRETE COMPRESSION HARDENING
*CONCRETE TENSION STIFFENING, TYPE=DISPLACEMENT
*CONCRETE TENSION DAMAGE, TYPE=DISPLACEMENT
*BOUNDARY
 NBASE, 1,2
*AMPLITUDE, NAME=HAMP, INPUT=koyna_haccel.inp
*AMPLITUDE, NAME=VAMP, INPUT=koyna_vaccel.inp
**
*STEP, NLGEOM, UNSYMM=YES
 STEP 1 - GRAVITY LOAD
*STATIC
 1.0E-10, 1.0E-10
*DLOAD
 DAM, GRAV, 9.81, 0, -1
*OUTPUT, FIELD, VARIABLE=PRESELECT, FREQ=10
*ELEMENT OUTPUT
 S,PE,LE,PEEQ,PEEQT,DAMAGEC,DAMAGET,SDEG
*OUTPUT, HISTORY, VARIABLE=PRESELECT, FREQ=1
*ELEMENT OUTPUT, ELSET=EOUT
 SP1
*NODE OUTPUT, NSET=NOUT
 U
*END STEP
**
*STEP, NLGEOM, UNSYMM=YES
 STEP 2 - HYDROSTATIC LOAD
*STATIC
 1.0E-10, 1.0E-10
*DLOAD
 WDAM, HP4, 900067.6, 91.75, 0
*END STEP
**
*STEP, INC=2000, NLGEOM, UNSYMM=YES
 STEP 3 - EARTHQUAKE
*DYNAMIC, HAFTOL=1.0E7
 0.02, 10.0, 1E-15,0.02
*BOUNDARY,TYPE=ACCELERATION,AMPLITUDE=HAMP
 NBASE,1,1,9.81
*BOUNDARY,TYPE=ACCELERATION,AMPLITUDE=VAMP
 NBASE,2,2,9.81
*CONTROLS, ANALYSIS=DISCONTINUOUS
*END STEP