*HEADING
[ACIN3DN, edge-surface impedances]
*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
*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
*Nset,nset=cloadset
1,101,201,301,
*element,type=acin3d4,elset=inf
1,  1,2,3,4
*element,type=acin3d3,elset=inf
2, 101,102,103
*surface,name=acinE3d4
1,E1
*surface,name=acinS3d4
1,SPOS
*element,type=acin3d3,elset=inf
2, 101,102,103
*surface,name=acinE3d3
2,E3
*surface,name=acinS3d3
2,SPOS
*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
3,SPOS
5,SPOS
*element,type=acin3d8,elset=inf
4, 301,302,303,304,  305,306,307,308
*surface,name=acine3d8
4,E1
*surface,name=acins3d8
4,SPOS
***************************************
*SOLID SECTION,MATERIAL=AIR320,ELSET=INF
1.0, .5, .5,-1.0e2
*endpart
***************************************
*assembly,name=assy1
*instance,name=PwithEdgeZ,part=plates
0.0,0.0,0.0
*end instance
*instance,name=withEdgeZ,part=aco
0.0,0.0,0.0
*end instance
*instance,name=PnoneEdgeZ,part=plates
2.0,0.0,0.0
*end instance
*instance,name=noneEdgeZ,part=aco
2.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
*************************************
*STEP,inc=1000
*STEADYSTATEDYNAMICS,direct,frequency=linear
1.,200.,20
*simpedance,property=Half_Zmat
assy1.withEdgeZ.acine3d4,
assy1.withEdgeZ.acine3d3,
assy1.withEdgeZ.acine3d6,
assy1.withEdgeZ.acine3d8,
*BOUNDARY,type=acceleration
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=1
*NODE OUTPUT
POR,
*OUTPUT,FIELD,FREQUENCY=99999
*ELEMENT OUTPUT
*OUTPUT,HISTORY,FREQUENCY=99999
*END STEP