*HEADING
 KOYNA DAM: SEISMIC ANALYSIS 
 Units - n, m, sec
***********************************************
** Step 3: Earthquake 
**         (without hydrodynamic interactions)
**
** Requires results from koyna_std_to_xpl, 
** FORTRAN subroutine koyna_hp_vdload, and 
** amplitude curves koyna_haccel.inp and 
** koyna_vaccel.inp
**
** Execution command:
** abaqus job=koyna_xpl oldjob=koyna_std_to_xpl 
**        user=koyna_hp_vdload double
***********************************************
*IMPORT, STEP=2, INCREMENT=1, STATE=YES, UPDATE=NO
 DAM
*MATERIAL, NAME=CONCRETE
*ELASTIC
 3.1027E+10, 0.2
*DENSITY
 2643.0
*DAMPING, BETA=0.00323
*CONCRETE DAMAGED PLASTICITY
 36.31
*CONCRETE COMPRESSION HARDENING
 13.0E+6, 0.000
 24.1E+6, 0.001
*CONCRETE TENSION STIFFENING, TYPE=DISPLACEMENT
 2.9E+6         ,0
 1.94393E+6     ,0.000066185
 1.30305E+6     ,0.00012286
 0.873463E+6    ,0.000173427
 0.5855E+6      ,0.00022019
 0.392472E+6    ,0.000264718
 0.263082E+6    ,0.000308088
 0.176349E+6    ,0.00035105
 0.11821E+6     ,0.000394138
 0.0792388E+6   ,0.000437744
 0.0531154E+6   ,0.000482165
*CONCRETE TENSION DAMAGE, TYPE=DISPLACEMENT
 0          ,0
 0.381217   ,0.000066185
 0.617107   ,0.00012286
 0.763072   ,0.000173427
 0.853393   ,0.00022019
 0.909282   ,0.000264718
 0.943865   ,0.000308088
 0.965265   ,0.00035105
 0.978506   ,0.000394138
 0.9867     ,0.000437744
 0.99177    ,0.000482165
*AMPLITUDE, NAME=HAMP, INPUT=koyna_haccel.inp
*AMPLITUDE, NAME=VAMP, INPUT=koyna_vaccel.inp
*STEP, NLGEOM=YES
 STEP 3 - EARTHQUAKE
*DYNAMIC, EXPLICIT
 , 10.0
*DLOAD
 DAM, GRAV, 9.81, 0, -1
 WDAM, P4NU
*BOUNDARY,TYPE=ACCELERATION,AMPLITUDE=HAMP
 NBASE,1,1,9.81
*BOUNDARY,TYPE=ACCELERATION,AMPLITUDE=VAMP
 NBASE,2,2,9.81
**
*OUTPUT, FIELD, VAR=PRESELECT, NUMBER INTERVAL=40
*ELEMENT OUTPUT
 S,PE,LE,PEEQ,PEEQT,DAMAGEC,DAMAGET,SDEG
*OUTPUT,HISTORY, VAR=PRESELECT
*ELEMENT OUTPUT, ELSET=EOUT
 SP1
*NODE OUTPUT, NSET=NOUT
 U1,U2
*END STEP