*HEADING
 Test of fluid cavity and fluid exchange
Type of fluid exchange behavior: Bulk Viscosity
Cv=10000; Ch=100
*NODE
*NODE,NSET=FREF
*NGEN,LINE=C
*ELEMENT,TYPE=SAX1,ELSET=CAN1
*ELGEN,ELSET=CAN1
*ELEMENT,TYPE=SAX1,ELSET=CAN2
*ELGEN,ELSET=CAN2
*ELSET, ELSET=CAV1SURF_DEF, GENERATE
*SURFACE, TYPE=ELEMENT, NAME=CAV1SURF
CAV1SURF_DEF, SPOS
*ELSET, ELSET=CAV2SURF_DEF, GENERATE
*SURFACE, TYPE=ELEMENT, NAME=CAV2SURF
CAV2SURF_DEF, SPOS
*ELEMENT,TYPE=RAX2,ELSET=IMPACTOR
*ELEMENT,TYPE=RAX2,ELSET=END
*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
*CAPACITY, TYPE=TABULAR
*FLUID EXCHANGE, NAME=TUBE, PROPERTY=PROP_DEF
*FLUID EXCHANGE PROPERTY, NAME=PROP_DEF, TYPE=BULK VISCOSITY
*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
*BIAXIAL TEST DATA
*PLANAR TEST DATA
*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
*OUTPUT,FIELD,VARIABLE=PRESELECT,NUMBER INTERVAL=4
*OUTPUT,HISTORY,TIME INTERVAL=0.00064
*NODE OUTPUT,NSET=FREF
PCAV,CVOL
**
** 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,CVOL
*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