*HEADING
Test of fluid cavity and fluid exchange
Type of fluid exchange behavior: Bulk Viscosity
Cv=10000; Ch=100
*Preprint, Contact=Yes, Model=Yes, History=Yes
*NODE,NSET=ALL
*NODE,NSET=FREF
*NGEN,LINE=C,NSET=ALL
*ELEMENT,TYPE=M3D4R,ELSET=CAN1
*ELGEN,ELSET=CAN1
*ELEMENT,TYPE=M3D4R,ELSET=CAN2
*ELGEN,ELSET=CAN2
*ELEMENT,TYPE=R3D3,ELSET=IMPACTOR
*ELEMENT,TYPE=R3D3,ELSET=END
*NSET,NSET=MID
*BOUNDARY
79991,1
79991,3,6
5001,XSYMM
5019,XSYMM
6001,XSYMM
6019,XSYMM
19991,1,6
17001,XSYMM
17001,3,4
15001,XSYMM
15019,XSYMM
16001,XSYMM
16019,XSYMM
11001,1,6
ALL,3
MID,1,3
*ELEMENT,TYPE=SFM3D3,ELSET=FLU1
*ELEMENT,TYPE=SFM3D3,ELSET=FLU2
*ELSET, ELSET=CAV1SURF_DEF
*SURFACE, TYPE=ELEMENT, NAME=CAV1SURF
CAV1SURF_DEF, SPOS
*ELSET, ELSET=CAV2SURF_DEF
*SURFACE, TYPE=ELEMENT, NAME=CAV2SURF
CAV2SURF_DEF, SPOS
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16, UNIVERSAL GAS CONSTANT=8.314
*FLUID CAVITY, NAME=CAV1, REF NODE=991, SURFACE=CAV1SURF, AMBIENT=50000., BEHAVIOR=GAS
*FLUID CAVITY, NAME=CAV2, REF NODE=992, SURFACE=CAV2SURF, AMBIENT=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
*ELEMENT,TYPE=MASS,ELSET=MASS
*MASS,ELSET=MASS
1.E-10
*MEMBRANE SECTION,ELSET=CAN1,MATERIAL=ST
0.05,
*MEMBRANE SECTION,ELSET=CAN2,MATERIAL=ST
0.05,
*SURFACE, TYPE=ELEMENT, NAME=FLU1SURF
FLU1,
*SURFACE, TYPE=ELEMENT, NAME=FLU2SURF
FLU2,
*SURFACE SECTION, ELSET=FLU1
*SURFACE SECTION, ELSET=FLU2
*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
991, 8273.6
992, 8273.6
*INITIAL CONDITIONS,TYPE=STRESS
CAN1,82736.,41368.
CAN2,82736.,41368.
*RESTART,WRITE,NUM=2,TIMEMARKS=NO
*RIGID BODY,ELSET=IMPACTOR,REF NODE=79991
*RIGID BODY,ELSET=END,REF NODE=19991
*************************************************************
**General contact definition
*Contact
*Contact Inclusions, All Element Based
*Contact Exclusions
,FLU1SURF
,FLU2SURF
*************************************************************
*STEP
*DYNAMIC,EXPLICIT
,0.64
*BOUNDARY,TYPE=VELOCITY
79991,2,,-4.
*FLUID EXCHANGE INTERACTION
TUBE
*NSET, NSET=RFCHECK
*OUTPUT, FIELD, NUMBER INTERVAL=1
*ELEMENT OUTPUT
S,SP,LE,LEP,ERV
*NODE OUTPUT
PCAV,CVOL
U,
*CONTACT OUTPUT
CSTRESS, CFORCE
*OUTPUT,HISTORY,TIME INTERVAL=0.00064
*NODE OUTPUT,NSET=FREF
PCAV,CVOL
**
** Output requests for qa testing
**
*NSET,NSET=QA_TEST1
FREF,
*OUTPUT,FIELD,NUM=1
*NODE OUTPUT,NSET=QA_TEST1
U,CVOL
*OUTPUT,HISTORY,VAR=PRESELECT,TIME INTERVAL=0.64
*END STEP