*HEADING
[ACIN3DN, edge-surface impedances, transient]
*preprint,model=yes,history=yes,echo=yes
********************************
*part,name=plates
*NODE, NSET=first
1, 0., 0., 0.,0,0,1
2, 1., 0., 0.,0,0,1
3, 1., 1., 0.,0,0,1
4, 0., 1., 0.,0,0,1
5, .5, 0., 0.,0,0,1
6, 1., .5, 0.,0,0,1
7, .5, 1., 0.,0,0,1
8, 0., .5, 0.,0,0,1
9, .5, .5, 0.,0,0,1
*********************************
*NCOPY,CHANGENUMBER=100,OLDSET=first,NEWSET=nall,SHIFT
0,2,0,
1,0,0, 0,1,0, 0
*NCOPY,CHANGENUMBER=200,OLDSET=first,NEWSET=nall,SHIFT
0,4,0
1,0,0, 0,1,0, 0
*NCOPY,CHANGENUMBER=300,OLDSET=first,NEWSET=nall,SHIFT
0,6,0,
1,0,0, 0,1,0, 0
*NCOPY,CHANGENUMBER=400,OLDSET=first,NEWSET=nall,SHIFT
0,6,0,
1,0,0, 0,1,0, 0
*NSET,NSET=ALLNODES
FIRST,NALL
*nset,nset=qa_test,gen
1,70000,1
*element,type=m3d4,elset=mass
201,1,2,3,4
*element,type=m3d3,elset=mass
202,101,102,103
*element,type=m3d8,elset=mass
203,201,202,203,204,  205,206,207,208
*element,type=m3d8,elset=mass
204,301,302,303,304,  305,306,307,308
*surface,name=sm3d4
201,spos
*surface,name=sm3d3
202,spos
*surface,name=sm3d6
203,spos
*surface,name=sm3d8
204,spos
*MEMBRANESECTION,MATERIAL=METAL,ELSET=MASS
0.05,
*end part
********************************
*part,name=aco
*NODE, NSET=first
1, 0., 0., 0.,0,0,1
2, 1., 0., 0.,0,0,1
3, 1., 1., 0.,0,0,1
4, 0., 1., 0.,0,0,1
5, .5, 0., 0.,0,0,1
6, 1., .5, 0.,0,0,1
7, .5, 1., 0.,0,0,1
8, 0., .5, 0.,0,0,1
9, .5, .5, 0.,0,0,1
*********************************
*NCOPY,CHANGENUMBER=100,OLDSET=first,NEWSET=nall,SHIFT
0,2,0,
1,0,0, 0,1,0, 0
*NCOPY,CHANGENUMBER=200,OLDSET=first,NEWSET=nall,SHIFT
0,4,0
1,0,0, 0,1,0, 0
*NCOPY,CHANGENUMBER=300,OLDSET=first,NEWSET=nall,SHIFT
0,6,0,
1,0,0, 0,1,0, 0
*NSET,NSET=ALLNODES
FIRST,NALL
*ncopy,changenumber=500, oldset=allnodes,newset=nall2,shift
0,0,-1
1,0,0, 0,1,0, 0
*ncopy,changenumber=1000,oldset=allnodes,newset=nall3,shift
0,0,-.5
1,0,0, 0,1,0, 1
*NSET,NSET=ALLNODES
FIRST,NALL,NALL2,NALL3
*nset,nset=qa_test
1,101,201,301
*Nset,nset=cloadset
1,101,201,301,
*************
*element,type=ac3d8, elset=inf
11, 1,4,3,2,  501,504,503,502
*element,type=acin3d4,elset=inf
1,  1,2,3,4
*surface,name=acinE3d4
1, E1
*surface,name=acinS3d4
11,s2
*************
*element,type=ac3d6, elset=inf
12, 101,103,102, 601,603,602
*element,type=acin3d3,elset=inf
2, 101,102,103
*surface,name=acinE3d3
2,E3
*surface,name=acinS3d3
12,s2
*************
*element,type=ac3d15,elset=inf
13, 201,203,202, 701,703,702, 209,206,205, 
    709,706,705, 1201,1203,1202
15, 201,204,203, 701,704,703, 208,207,209, 
    708,707,709, 1201,1204,1203
*element,type=acin3d6,elset=inf
3, 201,202,203,   205,206,209
5, 201,203,204,   209,207,208
*surface,name=acine3d6
3,E1
*surface,name=acins3d6
13,s2
15,s2
**************
*element,type=ac3d20, elset=inf
14,301,304,303,302,  801,804,803,802, 308,307,306,305,
   808,807,806,805,  1301,1304,1303,1302
*element,type=acin3d8,elset=inf
4, 301,302,303,304,  305,306,307,308
*surface,name=acine3d8
4,E1
*surface,name=acins3d8
14,s2
***************************************
*SOLID SECTION,MATERIAL=AIR320,ELSET=INF
1.0, 1.0,1.0,-1.0e2
*endpart
***************************************
*assembly,name=assy1
*instance,name=PwithEdgeZ,part=plates
0.0,0.0,-1.0
*end instance
*instance,name=withEdgeZ,part=aco
0.0,0.0,0.0
*end instance
*instance,name=PnoneEdgeZ,part=plates
1.0,0.0,-1.0
*end instance
*instance,name=noneEdgeZ,part=aco
1.0,0.0,0.0
*end instance
*tie,name=quad4
withEdgeZ.acins3d4,PwithEdgeZ.sm3d4
noneEdgeZ.acins3d4,PnoneEdgeZ.sm3d4
withEdgeZ.acins3d3,PwithEdgeZ.sm3d3
noneEdgeZ.acins3d3,PnoneEdgeZ.sm3d3
withEdgeZ.acins3d6,PwithEdgeZ.sm3d6
noneEdgeZ.acins3d6,PnoneEdgeZ.sm3d6
withEdgeZ.acins3d8,PwithEdgeZ.sm3d8
noneEdgeZ.acins3d8,PnoneEdgeZ.sm3d8
*end assembly
*************************************
** MATERIAL IMPEDANCE FOR AIR320 IS 400.0
** THE ADMITTANCE FOR HALF THIS VALUE IS USED
**
*impedance property,name=Half_Zmat
0.0, 0.005
*MATERIAL,NAME=AIR320
*DENSITY
1.25, 
*ACOUSTIC MEDIUM,BULK MODULUS
1.28E5,
*material,name=metal
*density
7800.0,
*elastic
200.0e9,0.3
*amplitude,definition=perio,name=TenPiHz
1, 5.0, 0.0, 0.0
0.0,1.0,0.0,0.0
*************************************
*STEP,inc=1000
*dynamic,direct
5.e-2, 50.e-1
*simpedance,property=Half_Zmat
assy1.withEdgeZ.acine3d4,
assy1.withEdgeZ.acine3d3,
assy1.withEdgeZ.acine3d6,
assy1.withEdgeZ.acine3d8,
*BOUNDARY,type=acceleration,amplitude=tenpihz
ASSY1.PWITHEDGEZ.ALLNODES,3,3,1.0e-4
ASSY1.PNONEEDGEZ.ALLNODES,3,3,1.0e-4
*BOUNDARY
ASSY1.PWITHEDGEZ.ALLNODES,1,2,0.0
ASSY1.PNONEEDGEZ.ALLNODES,1,2,0.0
*OUTPUT,FIELD,FREQUENCY=100
*NODE OUTPUT
POR,
*OUTPUT,FIELD,FREQUENCY=99999
*ELEMENT OUTPUT
*OUTPUT,HISTORY,FREQUENCY=1
*NODE OUTPUT,NSET=ASSY1.WITHEDGEZ.QA_TEST
POR,
*NODE OUTPUT,NSET=ASSY1.NONEEDGEZ.QA_TEST
POR,
*END STEP