*HEADING
[ACIN3DN, edge-surface impedances, transient]
*preprint,model=yes,history=yes,echo=yes
********************************
*part,name=plates
*NODE, NSET=first
*********************************
*NCOPY,CHANGENUMBER=100,OLDSET=first,NEWSET=nall,SHIFT
*NCOPY,CHANGENUMBER=200,OLDSET=first,NEWSET=nall,SHIFT
*NCOPY,CHANGENUMBER=300,OLDSET=first,NEWSET=nall,SHIFT
*NCOPY,CHANGENUMBER=400,OLDSET=first,NEWSET=nall,SHIFT
*NSET,NSET=ALLNODES
FIRST,NALL
*nset,nset=qa_test,gen
*element,type=m3d4,elset=mass
*element,type=m3d3,elset=mass
*element,type=m3d8,elset=mass
*element,type=m3d8,elset=mass
*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
*end part
********************************
*part,name=aco
*NODE, NSET=first
*********************************
*NCOPY,CHANGENUMBER=100,OLDSET=first,NEWSET=nall,SHIFT
*NCOPY,CHANGENUMBER=200,OLDSET=first,NEWSET=nall,SHIFT
*NCOPY,CHANGENUMBER=300,OLDSET=first,NEWSET=nall,SHIFT
*NSET,NSET=ALLNODES
FIRST,NALL
*ncopy,changenumber=500, oldset=allnodes,newset=nall2,shift
*ncopy,changenumber=1000,oldset=allnodes,newset=nall3,shift
*NSET,NSET=ALLNODES
FIRST,NALL,NALL2,NALL3
*nset,nset=qa_test
*Nset,nset=cloadset
*************
*element,type=ac3d8, elset=inf
*element,type=acin3d4,elset=inf
*surface,name=acinE3d4
1, E1
*surface,name=acinS3d4
11,s2
*************
*element,type=ac3d6, elset=inf
*element,type=acin3d3,elset=inf
*surface,name=acinE3d3
2,E3
*surface,name=acinS3d3
12,s2
*************
*element,type=ac3d15,elset=inf
*element,type=acin3d6,elset=inf
*surface,name=acine3d6
3,E1
*surface,name=acins3d6
13,s2
15,s2
**************
*element,type=ac3d20, elset=inf
*element,type=acin3d8,elset=inf
*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
*end instance
*instance,name=withEdgeZ,part=aco
*end instance
*instance,name=PnoneEdgeZ,part=plates
*end instance
*instance,name=noneEdgeZ,part=aco
*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
*MATERIAL,NAME=AIR320
*DENSITY
1.25,
*ACOUSTIC MEDIUM,BULK MODULUS
*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
5.e-2, 50.e-1
*simpedance
assy1.withEdgeZ.acine3d4,Half_Zmat
assy1.withEdgeZ.acine3d3,Half_Zmat
assy1.withEdgeZ.acine3d6,Half_Zmat
assy1.withEdgeZ.acine3d8,Half_Zmat
*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