*HEADING Test of fluid cavity and fluid exchange Type of fluid exchange behavior: Bulk Viscosity Cv=10000; Ch=100 *NODE,NSET=ALL 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 15001,0.0,0.,1. 15010,1.,1.,1. 15019,0.,2.,1. 16001,0.0,2.125,1. 16010,1.0,3.125,1. 16019,0.,4.125,1. 17001,-0.01,4.13,1. 17002,2.,4.13,1. 11001,0.,-0.0001,1. 11002,2.,-0.0001,1. 1,0.,1. 2,0.,3.125 10001,0.,1.,1. 10002,0.,3.125,1. 79991,-0.01,4.13 19991,0.,-0.0001 *NODE,NSET=FREF 991,0.,1.,0.5 992,0.,3.125,0.5 *NGEN,LINE=C,NSET=ALL 5001,5010,1,,0.,1. 5010,5019,1,,0.,1. 6001,6010,1,,0.,3.125 6010,6019,1,,0.,3.125 15001,15010,1,,0.,1.,1. 15010,15019,1,,0.,1.,1. 16001,16010,1,,0.,3.125,1. 16010,16019,1,,0.,3.125,1. *ELEMENT,TYPE=M3D4R,ELSET=CAN1 5001,5002,5001,15001,15002 *ELGEN,ELSET=CAN1 5001,18,1,1 *ELEMENT,TYPE=M3D4R,ELSET=CAN2 6001,6002,6001,16001,16002 *ELGEN,ELSET=CAN2 6001,18,1,1 *ELEMENT,TYPE=R3D3,ELSET=IMPACTOR 7001,7001,7002,17002 7002,17002,17001,7001 *ELEMENT,TYPE=R3D3,ELSET=END 1001,1001,1002,11002 1002,11002,11001,1001 *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 *ELEMENT,TYPE=SFM3D3,ELSET=FLU1 8021,5001,5002,1 8022,5002,5003,1 8023,5003,5004,1 8024,5004,5005,1 8025,5005,5006,1 8026,5006,5007,1 8027,5007,5008,1 8028,5008,5009,1 8029,5009,5010,1 8030,5010,5011,1 8031,5011,5012,1 8032,5012,5013,1 8033,5013,5014,1 8034,5014,5015,1 8035,5015,5016,1 8036,5016,5017,1 8037,5017,5018,1 8038,5018,5019,1 8041,15019,15018,10001 8042,15018,15017,10001 8043,15017,15016,10001 8044,15016,15015,10001 8045,15015,15014,10001 8046,15014,15013,10001 8047,15013,15012,10001 8048,15012,15011,10001 8049,15011,15010,10001 8050,15010,15009,10001 8051,15009,15008,10001 8052,15008,15007,10001 8053,15007,15006,10001 8054,15006,15005,10001 8055,15005,15004,10001 8056,15004,15003,10001 8057,15003,15002,10001 8058,15002,15001,10001 *ELEMENT,TYPE=SFM3D3,ELSET=FLU2 9021,6001,6002,2 9022,6002,6003,2 9023,6003,6004,2 9024,6004,6005,2 9025,6005,6006,2 9026,6006,6007,2 9027,6007,6008,2 9028,6008,6009,2 9029,6009,6010,2 9030,6010,6011,2 9031,6011,6012,2 9032,6012,6013,2 9033,6013,6014,2 9034,6014,6015,2 9035,6015,6016,2 9036,6016,6017,2 9037,6017,6018,2 9038,6018,6019,2 9040,16019,16018,10002 9041,16018,16017,10002 9042,16017,16016,10002 9043,16016,16015,10002 9044,16015,16014,10002 9045,16014,16013,10002 9046,16013,16012,10002 9047,16012,16011,10002 9048,16011,16010,10002 9049,16010,16009,10002 9050,16009,16008,10002 9051,16008,16007,10002 9052,16007,16006,10002 9053,16006,16005,10002 9054,16005,16004,10002 9055,16004,16003,10002 9056,16003,16002,10002 9057,16002,16001,10002 *ELSET, ELSET=CAV1SURF_DEF 5001, 5002, 5003, 5004, 5005, 5006, 5007, 5008, 5009, 5010, 5011, 5012, 5013, 5014, 5015, 5016 5017, 5018, 8021, 8022, 8023, 8024, 8025, 8026, 8027, 8028, 8029, 8030, 8031, 8032, 8033, 8034 8035, 8036, 8037, 8038, 8041, 8042, 8043, 8044, 8045, 8046, 8047, 8048, 8049, 8050, 8051, 8052 8053, 8054, 8055, 8056, 8057, 8058 *SURFACE, TYPE=ELEMENT, NAME=CAV1SURF CAV1SURF_DEF, SPOS *ELSET, ELSET=CAV2SURF_DEF 6001, 6002, 6003, 6004, 6005, 6006, 6007, 6008, 6009, 6010, 6011, 6012, 6013, 6014, 6015, 6016 6017, 6018, 9021, 9022, 9023, 9024, 9025, 9026, 9027, 9028, 9029, 9030, 9031, 9032, 9033, 9034 9035, 9036, 9037, 9038, 9040, 9041, 9042, 9043, 9044, 9045, 9046, 9047, 9048, 9049, 9050, 9051 9052, 9053, 9054, 9055, 9056, 9057 *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 PRESSURE=50000., BEHAVIOR=GAS *FLUID CAVITY, NAME=CAV2, REF NODE=992, 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 991,992 *FLUID EXCHANGE PROPERTY, NAME=PROP_DEF, TYPE=BULK VISCOSITY 10000., 100. *ELEMENT,TYPE=MASS,ELSET=MASS 100000,1 100001,2 100002,10001 100003,10002 *MASS,ELSET=MASS 1.E-10 *MEMBRANE SECTION,ELSET=CAN1,MATERIAL=ST 0.05, *MEMBRANE SECTION,ELSET=CAN2,MATERIAL=ST 0.05, *SURFACE SECTION, ELSET=FLU1 *SURFACE SECTION, ELSET=FLU2 *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 991, 8273.6 992, 8273.6 *INITIAL CONDITIONS,TYPE=STRESS CAN1,82736.,41368. CAN2,82736.,41368. *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,SNEG *SURFACE,TYPE=ELEMENT,NAME=RIGID IMPACTOR,SPOS *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, OP=ADD, TIMEMARKS=NO, NUMBER INTERVAL=1 *ELEMENT OUTPUT S,SP,LE,LEP,ERV *NODE OUTPUT U, *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, *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