*HEADING
 Test of fluid cavity and fluid exchange
Type of fluid exchange behavior: Bulk Viscosity
Cv=10000; Ch=100
*NODE
5001,0.0,0.
5010,1.,1.
5019,0.,2.
6001,0.0,2.125
6010,1.0,3.125
6019,0.,4.125
7001,-0.01,4.13
7002,2.,4.13
1001,0.,-0.0001
1002,2.,-0.0001
79991,-0.01,4.13
19991,0.,-0.0001
*NODE,NSET=FREF
1,0.,1.
2,0.,3.125
*NGEN,LINE=C
5001,5010,1,,0.,1.
5010,5019,1,,0.,1.
6001,6010,1,,0.,3.125
6010,6019,1,,0.,3.125
*ELEMENT,TYPE=SAX1,ELSET=CAN1
5001,5001,5002
*ELGEN,ELSET=CAN1
5001,18,1,1
*ELEMENT,TYPE=SAX1,ELSET=CAN2
6001,6001,6002
*ELGEN,ELSET=CAN2
6001,18,1,1
*ELSET, ELSET=CAV1SURF_DEF, GENERATE
 5001,  5018,     1
*SURFACE, TYPE=ELEMENT, NAME=CAV1SURF
CAV1SURF_DEF, SPOS
*ELSET, ELSET=CAV2SURF_DEF, GENERATE
 6001,  6018,     1
*SURFACE, TYPE=ELEMENT, NAME=CAV2SURF
CAV2SURF_DEF, SPOS
*ELEMENT,TYPE=RAX2,ELSET=IMPACTOR
7001,7001,7002
*ELEMENT,TYPE=RAX2,ELSET=END
1001,1001,1002
*BOUNDARY
79991,1
79991,3,6
5001,XSYMM
5019,XSYMM
6001,XSYMM
6019,XSYMM
19991,1,6
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16, UNIVERSAL GAS CONSTANT=8.314
*FLUID CAVITY, NAME=CAV1, REF NODE=1, SURFACE=CAV1SURF, AMBIENT PRESSURE=50000., BEHAVIOR=GAS

*FLUID CAVITY, NAME=CAV2, REF NODE=2, SURFACE=CAV2SURF, AMBIENT PRESSURE=50000., BEHAVIOR=GAS

*FLUID BEHAVIOR, NAME=GAS
*MOLECULAR WEIGHT
31.178
*CAPACITY, TYPE=TABULAR
28.110
*FLUID EXCHANGE, NAME=TUBE, PROPERTY=PROP_DEF
1, 2
*FLUID EXCHANGE PROPERTY, NAME=PROP_DEF, TYPE=BULK VISCOSITY
10000., 100.
*SHELL SECTION,ELSET=CAN1,MATERIAL=ST,SECTION INTEGRATION=GAUSS
0.05,
*SHELL SECTION,ELSET=CAN2,MATERIAL=ST,SECTION INTEGRATION=GAUSS
0.05,
*MATERIAL,NAME=ST
*DENSITY
1000., 
*HYPERELASTIC,TEST DATA INPUT,OGDEN,N=3
*UNIAXIAL TEST DATA
  1.5506E5,  0.1338
  2.4367E5,  0.2675
  3.1013E5,  0.3567
  4.2089E5,  0.6242
  5.3165E5,  0.8917
  5.9810E5,  1.1592
  6.8671E5,  1.4268
  8.8608E5,  2.0510
 10.6329E5,  2.5860
 12.4051E5,  3.0318
 16.1709E5,  3.7898
 19.9367E5,  4.3694
 23.4810E5,  4.8153
 27.4684E5,  5.1720
 31.0127E5,  5.4395
 34.5570E5,  5.7070
 38.3228E5,  5.9299
 42.0886E5,  6.0637
 45.6329E5,  6.1975
 49.3987E5,  6.3312
 53.1646E5,  6.4650
 56.9304E5,  6.5541
 64.2405E5,  6.6433
*BIAXIAL TEST DATA
  0.9384E5,  0.0200
  1.5900E5,  0.0600
  2.4087E5,  0.1100
  2.6220E5,  0.1400
  3.3240E5,  0.2000
  4.4278E5,  0.3100
  5.1830E5,  0.4200
  6.6024E5,  0.6800
  7.7794E5,  0.9400
  9.7857E5,  1.4900
 12.6351E5,  2.0300
 14.6804E5,  2.4300
 17.4000E5,  2.7500
 20.1058E5,  3.0700
 22.4502E5,  3.2600
 24.6530E5,  3.4500
*PLANAR TEST DATA
  0.6000E5,  0.0690
  1.6000E5,  0.1034
  2.4000E5,  0.1724
  3.3600E5,  0.2828
  4.2000E5,  0.4276
  6.0000E5,  0.8483
  7.8000E5,  1.3862
  9.6000E5,  2.0000
 11.1200E5,  2.4897
 12.9600E5,  3.0345
 14.8800E5,  3.4483
 16.5800E5,  3.7793
 18.2000E5,  4.0621
*INITIAL CONDITIONS,TYPE=FLUID PRESSURE
1, 8273.6
2, 8273.6
*INITIAL CONDITIONS,TYPE=STRESS
CAN1,82736.,82736.
CAN2,82736.,82736.
*RESTART,WRITE,NUM=2,TIMEMARKS=NO
*SURFACE,TYPE=ELEMENT,NAME=CANN1
CAN1,SNEG
*SURFACE,TYPE=ELEMENT,NAME=CANN2
CAN2,SNEG
*SURFACE,TYPE=ELEMENT,NAME=FIX
END,SPOS
*SURFACE,TYPE=ELEMENT,NAME=RIGID
IMPACTOR,SNEG
*RIGID BODY,ELSET=IMPACTOR,REF NODE=79991
*RIGID BODY,ELSET=END,REF NODE=19991
*STEP
*DYNAMIC,EXPLICIT
,0.64
*BOUNDARY,TYPE=VELOCITY
79991,2,,-4.
*CONTACT PAIR
RIGID,CANN2
*CONTACT PAIR
CANN1,CANN2
*CONTACT PAIR
FIX,CANN1
*FLUID EXCHANGE ACTIVATION
TUBE
*NSET, NSET=RFCHECK
79991, 19991
*OUTPUT,FIELD,VARIABLE=PRESELECT,NUMBER INTERVAL=4
*OUTPUT,HISTORY,TIME INTERVAL=0.00064
*NODE OUTPUT,NSET=FREF
PCAV,CVOL,CSAREA
**
** Output requests for qa testing
**
*NSET,NSET=QA_TEST1
FREF,
*NSET,NSET=QA_TEST2
RFCHECK,
*OUTPUT,FIELD,NUM=1
*NODE OUTPUT,NSET=QA_TEST1
U,
*NODE OUTPUT,NSET=QA_TEST2
RF,
*OUTPUT,HISTORY,VAR=PRESELECT,TIME INTERVAL=0.16
*FILE OUTPUT,NUMBER INTERVAL=1
*NODE FILE, NSET=FREF
U,CVOL
*NODE FILE, NSET=RFCHECK
RF, 
*END STEP