*HEADING 
     Gas structure interaction with structure modeled using 3D elements.

This test verifies the coupling between cavity fluid pressure and deformation 
of the structure enclosing the fluid cavity.
The deformation of the structure should change the fluid pressure according to
the fluid behavior definition and the pressure should get transferred to the 
structure.  
Fluid cavity pressure x Area of any face of the structure  = 
        Sum of RF in a direction normal to the selected face on all the nodes 
        constituting that face. 
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16, UNIVERSAL GAS CONSTANT=8.314
**=============================================================================
**Cavity I: Structure enclosing cavity modeled using S4R, SFM3D4R, M3D4R 
**          elements
**=============================================================================
*NODE, NSET=shellSurfMemQuad
11, 0, 0, 0
12, 2, 0, 0
13, 2, 2, 0
14, 0, 2, 0
15, 0, 0, 2
16, 2, 0, 2
17, 2, 2, 2
18, 0, 2, 2
*ELEMENT, TYPE=SFM3D4R, ELSET=quadSurf
11, 13, 17, 18, 14
12, 11, 12, 13, 14
*ELEMENT, TYPE=M3D4R, ELSET=quadMem
13, 12, 16, 17, 13
14, 11, 14, 18, 15
*ELEMENT, TYPE=S4R, ELSET=quadShell
15, 18, 17, 16, 15
16, 11, 15, 16, 12
*NSET, NSET=shellSurfMemQuadFixBC
11, 12, 15, 16
*NSET, NSET=shellSurfMemQuadDispBC
13, 14, 17, 18
**
**Cavity ref node
**
*NODE, NSET=shellSurfMemQuad_Ref
19999, 1, 1, 1
*ELSET, ELSET=shellSurfMemQuad
11, 12, 13, 14, 15, 16
*SURFACE, NAME=shellSurfMemQuad_Surf, TYPE=ELEMENT
shellSurfMemQuad, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=shellSurfMemQuad_Cav, REF NODE=shellSurfMemQuad_Ref, 
BEHAVIOR=AIR, SURFACE=shellSurfMemQuad_Surf, AMBIENT PRESSURE=1.0E5

**============================================================================
**Cavity II: Structure enclosing cavity modeled using S3R, SFM3D3, M3D3 
**           elements
**============================================================================
*SYSTEM
4, 0, 0
*NODE, NSET=shellSurfMemTri
21, 0, 0, 0
22, 2, 0, 0
23, 2, 2, 0
24, 0, 2, 0
25, 0, 0, 2
26, 2, 0, 2
27, 2, 2, 2
28, 0, 2, 2
29, 1, 1, 0
210, 2, 1, 1
211, 1, 1, 2
212, 0, 1, 1
213, 1, 2, 1
214, 1, 0, 1
**
**Cavity ref node
**
*NODE, NSET=shellSurfMemTri_Ref
29999, 1, 1, 1
*ELEMENT, TYPE=SFM3D3, ELSET=triSurf
21, 24, 23, 213
22, 23, 27, 213
23, 27, 28, 213
24, 28, 24, 213
25, 23, 24, 29
26, 24, 21, 29
27, 21, 22, 29
28, 22, 23, 29
*ELEMENT, TYPE=M3D3, ELSET=triMem
29, 23, 22, 210
210, 22, 26, 210
211, 26, 27, 210
212, 27, 23, 210
213, 21, 24, 212
214, 24, 28, 212
215, 28, 25, 212
216, 25, 21, 212
*ELEMENT, TYPE=S3R, ELSET=triShell
217, 28, 27, 211
218, 27, 26, 211
219, 26, 25, 211
220, 25, 28, 211
221, 22, 21, 214
222, 21, 25, 214
223, 25, 26, 214
224, 26, 22, 214
*NSET, NSET=shellSurfMemTriFixBC
21, 22, 25, 26, 29, 210, 211, 212, 214
*NSET, NSET=shellSurfMemTriDispBC
23, 24, 27, 28, 213
*ELSET, ELSET=shellSurfMemTri, GENERATE
21, 29, 1
210, 224, 1
*SURFACE, NAME=shellSurfMemTri_Surf, TYPE=ELEMENT
shellSurfMemTri, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=shellSurfMemTri_Cav, REF NODE=shellSurfMemTri_Ref, 
BEHAVIOR=AIR, SURFACE=shellSurfMemTri_Surf, AMBIENT PRESSURE=1.0E5

**============================================================================
**Cavity III: Structure enclosing cavity modeled using C3D8R, SC8R elements
**============================================================================
*SYSTEM
8, 0, 0
*NODE, NSET=brickShellQuad
31, 0, 0, 0
32, 2, 0, 0
33, 2, 2, 0
34, 0, 2, 0
35, 0, 0, 2
36, 2, 0, 2
37, 2, 2, 2
38, 0, 2, 2
**
**Cavity ref node
**
*NODE, NSET=brickShellQuad_Ref
39999, 1, 1, 1
*SYSTEM
7.9, 0, 0
*NODE, NSET=brickShellQuad
311, 0, 0, 0
314, 0, 2, 0
315, 0, 0, 2
318, 0, 2, 2
*SYSTEM
8.1, 0, 0
*NODE, NSET=brickShellQuad
312, 2, 0, 0
313, 2, 2, 0
316, 2, 0, 2
317, 2, 2, 2
*SYSTEM
8, 0.1, 0
*NODE, NSET=brickShellQuad
323, 2, 2, 0
324, 0, 2, 0
327, 2, 2, 2
328, 0, 2, 2
*SYSTEM
8, -0.1, 0
*NODE, NSET=brickShellQuad
321, 0, 0, 0
322, 2, 0, 0
325, 0, 0, 2
326, 2, 0, 2
*SYSTEM
8, 0, 0.1
*NODE, NSET=brickShellQuad
335, 0, 0, 2
336, 2, 0, 2
337, 2, 2, 2
338, 0, 2, 2
*SYSTEM
8, 0, -0.1
*NODE, NSET=brickShellQuad
331, 0, 0, 0
332, 2, 0, 0
333, 2, 2, 0
334, 0, 2, 0
*ELEMENT, TYPE=C3D8R, ELSET=brickQuad
31, 31, 35, 38, 34, 311, 315, 318, 314 
32, 31, 34, 33, 32, 331, 334, 333, 332
33, 32, 33, 37, 36, 312, 313, 317, 316
34, 35, 36, 37, 38, 335, 336, 337, 338
*ELEMENT, TYPE=SC8R, ELSET=continuumshellQuad
35, 33, 34, 38, 37, 323, 324, 328, 327
36, 31, 32, 36, 35, 321, 322, 326, 325
*NSET, NSET=brickShellQuadDispBC
33, 34, 38, 37, 323, 324, 328, 327
*NSET, NSET=brickShellQuadFixBC
31, 32, 35, 36
311, 312, 313, 314, 315, 316, 317, 318
321, 322, 325, 326
331, 332, 333, 334, 335, 336, 337, 338
*SURFACE, NAME=brickQuad_Surf, TYPE=ELEMENT
brickQuad, S1
*SURFACE, NAME=continuumshellQuad_Surf, TYPE=ELEMENT
continuumshellQuad, S1
*SURFACE, NAME=brickShellQuad_Surf, COMBINE=UNION
brickQuad_Surf, continuumshellQuad_Surf
**
** Fluid cavity definition
**
*FLUID CAVITY, NAME=brickShellQuad_Cav, REF NODE=brickShellQuad_Ref, 
BEHAVIOR=AIR, SURFACE=brickShellQuad_Surf, AMBIENT PRESSURE=1.0E5

**============================================================================
**Cavity IV: Structure enclosing cavity modeled using C3D6, SC6R elements
**============================================================================
*SYSTEM
12, 0, 0
*NODE, NSET=brickShellTri
411, 0, 0, 0
412, 2, 0, 0
413, 2, 2, 0
414, 0, 2, 0
415, 0, 0, 2
416, 2, 0, 2
417, 2, 2, 2
418, 0, 2, 2
419, 1, 1, 0
420, 2, 1, 1
421, 1, 1, 2
422, 0, 1, 1
423, 1, 2, 1
424, 1, 0, 1
**
**Cavity ref node
**
*NODE, NSET=brickShellTri_Ref
49999, 1, 1, 1
*SYSTEM
11.9, 0, 0
*NODE, NSET=brickShellTri
4111, 0, 0, 0
4114, 0, 2, 0
4115, 0, 0, 2
4118, 0, 2, 2
4122, 0, 1, 1
*SYSTEM
12.1, 0, 0
*NODE, NSET=brickShellTri
4112, 2, 0, 0
4113, 2, 2, 0
4116, 2, 0, 2
4117, 2, 2, 2
4120, 2, 1, 1
*SYSTEM
12, 0.1, 0
*NODE, NSET=brickShellTri
4213, 2, 2, 0
4214, 0, 2, 0
4217, 2, 2, 2
4218, 0, 2, 2
4223, 1, 2, 1
*SYSTEM
12, -0.1, 0
*NODE, NSET=brickShellTri
4211, 0, 0, 0
4212, 2, 0, 0
4215, 0, 0, 2
4216, 2, 0, 2
4224, 1, 0, 1
*SYSTEM
12, 0, 0.1
*NODE, NSET=brickShellTri
4315, 0, 0, 2
4316, 2, 0, 2
4317, 2, 2, 2
4318, 0, 2, 2
4321, 1, 1, 2
*SYSTEM
12, 0, -0.1
*NODE, NSET=brickShellTri
4311, 0, 0, 0
4312, 2, 0, 0
4313, 2, 2, 0
4314, 0, 2, 0
4319, 1, 1, 0
*ELEMENT, TYPE=C3D6, ELSET=brickTri
41, 414, 411, 422, 4114, 4111, 4122
42, 411, 415, 422, 4111, 4115, 4122
43, 415, 418, 422, 4115, 4118, 4122
44, 418, 414, 422, 4118, 4114, 4122
45, 412, 413, 420, 4112, 4113, 4120
46, 413, 417, 420, 4113, 4117, 4120
47, 417, 416, 420, 4117, 4116, 4120
48, 416, 412, 420, 4116, 4112, 4120
49, 411, 414, 419, 4311, 4314, 4319
410, 414, 413, 419, 4314, 4313, 4319
411, 413, 412, 419, 4313, 4312, 4319
412, 412, 411, 419, 4312, 4311, 4319
413, 415, 416, 421, 4315, 4316, 4321
414, 416, 417, 421, 4316, 4317, 4321 
415, 417, 418, 421, 4317, 4318, 4321 
416, 418, 415, 421, 4318, 4315, 4321
*ELEMENT, TYPE=SC6R, ELSET=continuumShellTri
417, 413, 414, 423, 4213, 4214, 4223
418, 414, 418, 423, 4214, 4218, 4223
419, 418, 417, 423, 4218, 4217, 4223
420, 417, 413, 423, 4217, 4213, 4223
421, 412, 416, 424, 4212, 4216, 4224
422, 416, 415, 424, 4216, 4215, 4224
423, 415, 411, 424, 4215, 4211, 4224
424, 411, 412, 424, 4211, 4212, 4224
*NSET, NSET=brickShellTriDispBC
413, 414, 417, 418, 423, 4213, 4214, 4217, 4218, 4223
*NSET, NSET=brickShellTriFixBC
411, 412, 415, 416, 419, 420, 421, 422, 424,
4111, 4112, 4113, 4114, 4115, 4116, 4117, 4118,
4211, 4212, 4215, 4216, 4311, 4312, 4313, 4314, 
4315, 4316, 4317, 4318, 4122, 4120, 4224, 4321, 
4319
*SURFACE, NAME=brickTri_Surf, TYPE=ELEMENT
brickTri, S1
*SURFACE, NAME=continuumShellTri_Surf, TYPE=ELEMENT
continuumShellTri, S1
*SURFACE, NAME=brickShellTri_Surf, COMBINE=UNION
brickTri_Surf, continuumShellTri_Surf
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=brickShellTri_Cav, REF NODE=brickShellTri_Ref, 
BEHAVIOR=AIR, SURFACE=brickShellTri_Surf, AMBIENT PRESSURE=1.0E5

**============================================================================
**Cavity V: Structure enclosing cavity modeled using C3D4 elements
**============================================================================
*SYSTEM
16, 0, 0
*NODE, NSET=tet
51, 0, 0, 0
52, 2, 0, 0
53, 2, 2, 0
54, 0, 2, 0
55, 0, 0, 2
56, 2, 0, 2
57, 2, 2, 2
58, 0, 2, 2
599, 1, 1, 0
5100, 2, 1, 1
5101, 1, 1, 2
5102, 0, 1, 1
5103, 1, 2, 1
5104, 1, 0, 1
5999, 1, 1, -1
51000, 3, 1, 1
51001, 1, 1, 3
51002, -1, 1, 1
51003, 1, 3, 1
51004, 1, -1, 1
**
**Cavity ref node
**
*NODE, NSET=tet_Ref
59999, 1, 1, 1
*ELEMENT, TYPE=C3D4, ELSET=tet
51, 52, 51, 599, 5999
52, 51, 54, 599, 5999
53, 54, 53, 599, 5999
54, 53, 52, 599, 5999
55, 52, 53, 5100, 51000
56, 53, 57, 5100, 51000
57, 57, 56, 5100, 51000
58, 56, 52, 5100, 51000
59, 56, 57, 5101, 51001
510, 57, 58, 5101, 51001
511, 58, 55, 5101, 51001
512, 55, 56, 5101, 51001
513, 54, 51, 5102, 51002
514, 51, 55, 5102, 51002
515, 55, 58, 5102, 51002
516, 58, 54, 5102, 51002
517, 53, 54, 5103, 51003
518, 54, 58, 5103, 51003
519, 58, 57, 5103, 51003
520, 57, 53, 5103, 51003
521, 51, 52, 5104, 51004
522, 52, 56, 5104, 51004
523, 56, 55, 5104, 51004
524, 55, 51, 5104, 51004
*NSET, NSET=tetFixBC
51, 52, 55, 56, 599, 5100, 5101, 5102, 5104, 5999, 51000, 51001, 51002, 51004
*NSET, NSET=tetDispBC
53, 54, 57, 58, 5103, 51003
*SURFACE, NAME=tet_Surf, TYPE=ELEMENT
tet, S1
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=tet_Cav, REF NODE=tet_Ref, BEHAVIOR=AIR, SURFACE=tet_Surf,
AMBIENT PRESSURE=1.0E5

**============================================================================
**Cavity VI: Structure enclosing cavity modeled using R3D3, R3D4 and M3D4R 
**          elements
**============================================================================
*SYSTEM
20, 0, 0
*NODE, NSET=rigid
61, 0, 0, 0
62, 2, 0, 0
63, 2, 2, 0
64, 0, 2, 0
65, 0, 0, 2
66, 2, 0, 2
67, 2, 2, 2
68, 0, 2, 2
69, 1, 0, 1
*NODE, NSET=rigidQuad_Ref
699, 1, 3, 1
*NODE, NSET=rigidTri_Ref
6100, 1, -1, 1
**
**Cavity ref node
**
*NODE, NSET=rigid_Ref
69999, 1, 1, 1
*ELEMENT, TYPE=M3D4R, ELSET=rigidMem
61, 61, 62, 63, 64
62, 62, 66, 67, 63
63, 68, 67, 66, 65
64, 61, 64, 68, 65
*ELEMENT, TYPE=R3D4, ELSET=rigidQuad
65, 63, 67, 68, 64
*ELEMENT, TYPE=R3D3, ELSET=rigidTri
66, 62, 61, 69
67, 61, 65, 69
68, 65, 66, 69
69, 66, 62, 69
*SURFACE, NAME=rigid_surf, TYPE=ELEMENT
rigidMem, SPOS
rigidQuad, SPOS
rigidTri, SPOS
**
**Fluid cavity definition
**
*FLUID CAVITY, NAME=rigid_Cav, REF NODE=rigid_Ref, BEHAVIOR=AIR, 
SURFACE=rigid_Surf, AMBIENT PRESSURE=1E5

**============================================================================
**Fluid behavior definition
**============================================================================
*FLUID BEHAVIOR, NAME=AIR
*MOLECULAR WEIGHT
0.0289
*CAPACITY, TYPE=POLYNOMIAL
28.110, 1.967E-3, 4.802E-6, -1.966E-9, 0.0
**============================================================================
**Boundary conditions
**============================================================================
*BOUNDARY
shellSurfMemQuadFixBC, 1, 6
shellSurfMemQuadDispBC, 1, 3
17, 4, 6
18, 4, 6
shellSurfMemTriFixBC, 1, 6
shellSurfMemTriDispBC, 1, 3
27, 4, 6
28, 4, 6
brickShellQuadFixBC, 1, 3
brickShellQuadDispBC, 1, 3
brickShellTriFixBC, 1, 3
brickShellTriDispBC, 1, 3
tetDispBC, 1, 3
tetFixBC, 1, 3
*BOUNDARY
rigidQuad_Ref, 1, 6
rigidTri_Ref, 1, 6
**============================================================================
**Initial conditions
**============================================================================
*INITIAL CONDITIONS, TYPE=TEMPERATURE
shellSurfMemQuad_Ref, 20
shellSurfMemTri_Ref, 20
brickShellQuad_Ref, 20
brickShellTri_Ref, 20
tet_Ref, 20
rigid_Ref, 20
*INITIAL CONDITIONS, TYPE=FLUID PRESSURE
shellSurfMemQuad_Ref, 4.0E5
shellSurfMemTri_Ref, 4.0E5
brickShellQuad_Ref, 4.0E5
brickShellTri_Ref, 4.0E5
tet_Ref, 4.0E5
rigid_Ref, 4.0E5
**============================================================================
**Section assignment
**============================================================================
*MEMBRANE SECTION, ELSET=quadMem, MATERIAL=MAT
0.1
*MEMBRANE SECTION, ELSET=triMem, MATERIAL=MAT
0.1
*MEMBRANE SECTION, ELSET=rigidMem, MATERIAL=MAT
0.1
*RIGID BODY, REF NODE=rigidQuad_Ref, ELSET=rigidQuad
*RIGID BODY, REF NODE=rigidTri_Ref, ELSET=rigidTri
*SHELL GENERAL SECTION, ELSET=quadShell, MATERIAL=MAT, OFFSET=SPOS
0.1
*SHELL GENERAL SECTION, ELSET=triShell, MATERIAL=MAT, OFFSET=SPOS
0.1
*SHELL GENERAL SECTION, ELSET=continuumshellQuad, MATERIAL=MAT, OFFSET=SPOS
0.1
*SHELL GENERAL SECTION, ELSET=continuumshellTri, MATERIAL=MAT, OFFSET=SPOS
0.1
*SOLID SECTION, ELSET=brickQuad, MATERIAL=MAT
*SOLID SECTION, ELSET=brickTri, MATERIAL=MAT
*SOLID SECTION, ELSET=tet, MATERIAL=MAT
*SURFACE SECTION, ELSET=quadSurf
*SURFACE SECTION, ELSET=triSurf
**============================================================================
**Material defintion
**============================================================================
*MATERIAL, NAME=MAT
*DENSITY
5000
*ELASTIC
2E-10, 0.1
**=============================================================================
**Create node sets
**=============================================================================
*NSET, NSET=qa_test_shellSurfMemQuad_Ref
shellSurfMemQuad_Ref
*NSET, NSET=qa_test_shellSurfMemTri_Ref
shellSurfMemTri_Ref
*NSET, NSET=qa_test_brickShellQuad_Ref
brickShellQuad_Ref
*NSET, NSET=qa_test_brickShellTri_Ref
brickShellTri_Ref
*NSET, NSET=qa_test_tet_Ref
tet_Ref
*NSET, NSET=qa_test_rigid_Ref
rigid_Ref
*NSET, NSET=qa_test_shellSurfMemQuadRF, ELSET=shellSurfMemQuad
*NSET, NSET=qa_test_shellSurfMemTriRF, ELSET=shellSurfMemTri
*NSET, NSET=qa_test_brickShellQuadRF
31, 32, 33, 34, 35, 36, 37, 38
*NSET, NSET=qa_test_brickShellTriRF
411, 412, 413, 414, 415, 416, 417, 418, 419, 420, 421, 422, 423, 424
*NSET, NSET=qa_test_tetRF
52, 53, 56, 57
*NSET, NSET=qa_test_rigidRF
rigidQuad_Ref, rigidTri_Ref
**============================================================================
** Amplitude definition
**============================================================================
*AMPLITUDE, NAME=RAMP
0, 0, 1, 1
**=============================================================================
**Step 1: Apply displacement BC of 0.50 to all the nodes of the face with
**        normal in +2 direction.
**        CVOL (at the end of step) = 10
**        PCAV (at the end of step) = 3.2E5
**        Area of the face with normal in +2 direction = 2*2 = 4
**        Area of face with normal in +1 or +3 direction = 2*(2+0.5) = 5
**        Sum of RF2 on the nodes of the face with normal in +2 direction =   
**                                                                      PCAV*4
**        Sum of RF1 on the nodes of the face with normal in +1 direction =   
**                                                                      PCAV*5
**        Sum of RF3 on the nodes of the face with normal in +3 direction =   
**                                                                      PCAV*5
**=============================================================================
*STEP
*DYNAMIC, EXPLICIT, DIRECT USER CONTROL
0.001, 1
*BOUNDARY, OP=MOD, AMPLITUDE=RAMP
shellSurfMemQuadDispBC, 2, 2, 0.50
shellSurfMemTriDispBC, 2, 2, 0.50
brickShellQuadDispBC, 2, 2, 0.50
brickShellTriDispBC, 2, 2, 0.50
tetDispBC, 2, 2, 0.50
rigidQuad_Ref, 2, 2, 0.50
*OUTPUT, HISTORY, TIME INTERVAL=0.25
*NODE OUTPUT, NSET=qa_test_shellSurfMemQuad_Ref
PCAV, CVOL
*NODE OUTPUT, NSET=qa_test_shellSurfMemTri_Ref
PCAV, CVOL
*NODE OUTPUT, NSET=qa_test_brickShellQuad_Ref
PCAV, CVOL
*NODE OUTPUT, NSET=qa_test_brickShellTri_Ref
PCAV, CVOL
*NODE OUTPUT, NSET=qa_test_tet_Ref
PCAV, CVOL
*NODE OUTPUT, NSET=qa_test_rigid_Ref
PCAV, CVOL
*NODE OUTPUT, NSET=qa_test_shellSurfMemQuadRF
RF1, RF2, RF3
*NODE OUTPUT, NSET=qa_test_shellSurfMemTriRF
RF1, RF2, RF3
*NODE OUTPUT, NSET=qa_test_brickShellQuadRF
RF1, RF2, RF3
*NODE OUTPUT, NSET=qa_test_brickShellTriRF
RF1, RF2, RF3
*NODE OUTPUT, NSET=qa_test_tetRF
RF1, RF2, RF3
*NODE OUTPUT, NSET=qa_test_rigidRF
RF2
*NODE OUTPUT, NSET=shellSurfMemQuadDispBC
U2
*NODE OUTPUT, NSET=shellSurfMemTriDispBC
U2
*NODE OUTPUT, NSET=brickShellQuadDispBC
U2
*NODE OUTPUT, NSET=brickShellTriDispBC
U2
*NODE OUTPUT, NSET=tetDispBC
U2
*NODE OUTPUT, NSET=rigidQuad_Ref
U2
*END STEP
**=============================================================================
**Step 2: Apply displacement BC of -0.50 to all the nodes of the face 
**        with normal in +2 direction.
**        All the parameters e.g. CVOL, PCAV, RF's at the end of this step 
**        will be the same as at the beginning of Step 1. 
**=============================================================================
*STEP
*DYNAMIC, EXPLICIT, DIRECT USER CONTROL
0.001, 1
*BOUNDARY, OP=MOD, AMPLITUDE=RAMP
shellSurfMemQuadDispBC, 2, 2, -0.50
shellSurfMemTriDispBC, 2, 2, -0.50
brickShellQuadDispBC, 2, 2, -0.50
brickShellTriDispBC, 2, 2, -0.50
tetDispBC, 2, 2, -0.50
rigidQuad_Ref, 2, 2, -0.50
*END STEP