*HEADING
THREE - D VERIFICATION PROBLEM FOR SURFACE CONTACT IN ACOUSTICS.
TRANSIENT
*PREPRINT,MODEL=YES,HISTORY=YES
*NODE, Nset=ODB_NSET
1
41,,10.
1001,1.
1041,1.,10.
2001,,,1.
2041,,10.,1.
3001,1.,,1.
3041,1.,10.,1.
4001,,,.5
4041,,10.,.5
5001,1.,,.5
5041,1.,10.,.5
6001,0.5,,0.5
6041,0.5,10,0.5
*NGEN,NSET=NALL1
1,41
2001,2041
*NGEN,NSET=NALL2
1001,1041
3001,3041
*NFILL,NSET=NALL
NALL1,NALL2,2,500
*NGEN,NSET=NALL3
4001,4041
5001,5041
6001,6041
*ELEMENT,TYPE=AC3D15,ELSET=ALL_ELEMENTS
    1,    1, 2001, 1001,    5, 2005, 1005,
         4001, 6001,  501, 4005, 6005,  505,
            3, 2003, 1003,
    2, 1001, 2001, 3001, 1005, 2005, 3005,
         6001, 2501, 5001, 6005, 2505, 5005,
         1003, 2003, 3003,
*ELGEN,ELSET=ALL_ELEMENTS
1,10,4,2
2,10,4,2
*ELSET,ELSET=LINER5
1,2
*NGEN,NSET=ALL_NODES
1,41
1001,1041
2001,2041
3001,3041
*NCOPY,OLDSET=ALL_NODES,NEWSET=ALL_NODES,SHIFT,CHANGENUMBER=10000
0.0,10.0,0.0
1.0,0.0,0.0,0.0,1.0,0.0,0.0
*ELEMENT,TYPE=AC3D4,ELSET=ALL_ELEMENTS
20001, 10001, 10003, 12003, 11001
20002, 10003, 12003, 11001, 11003
20003, 12003, 11001, 11003, 13003
20004, 10001, 12003, 12001, 11001
20005, 12003, 12001, 11001, 13003
20006, 12001, 11001, 13003, 13001
*ELGEN,ELSET=ALL_ELEMENTS
20001,20,2,6
20002,20,2,6
20003,20,2,6
20004,20,2,6
20005,20,2,6
20006,20,2,6
*elset,elset=liner4
20004
*elset,elset=liner6
20006
*MATERIAL,NAME=AIR
*DENSITY
1.293
*ACOUSTIC MEDIUM,BULK MODULUS
1.42176e5
*SOLID SECTION,MATERIAL=AIR,ELSET=ALL_ELEMENTS
*NSET,NSET=OUTPUT_NODES,GEN
1,41
10001,10041
*NSET,NSET=DRIVEN_NODES
10041,11041,12041,13041
***********************************************
*Elset,Elset=ODB_ESET
ALL_ELEMENTS,
*Nset, Nset=ODB_NSET,GENERATE
1,14000 
************************************* 
*AMPLITUDE,NAME=SINE,DEFINITION=PERIODIC 
1,628.0,0.0,0.0 
0.0,1.0,0.0,0.0
***********************************************
*SURFACE ,NAME=RIGHT
20,S2
19,S2
*SURFACE ,NAME=LEFT
LINER4,S4
LINER6,S2
*TIE,NAME=LR
LEFT,RIGHT
*************************************
*STEP,INC=100000
*DYNAMIC,DIRECT
5.e-4,0.4e-1,
*BOUNDARY,AMPLITUDE=SINE
DRIVEN_NODES,8,,1.
*IMPEDANCE
LINER5,I1
*OUTPUT,FIELD
*NODE OUTPUT
POR,
*NODE OUTPUT,NSET=ODB_NSET
POR,
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=ODB_NSET
POR,
*ENDSTEP