*HEADING
BENCHMARK PROBLEM: MODE 0 SPHERE USING ACIN3D3, ACIN3D4, ACIN3D6, ACIN3D8.
*PREPRINT,MODEL=YES
*part,name=inf3
*NODE,NSET=POLE3
*NODE,NSET=INNER3
*NCOPY,OLD SET=INNER3,CHANGE NUMBER=100,NEW SET=INNER3,SHIFT
*NCOPY,OLD SET=INNER3,CHANGE NUMBER=200,NEW SET=INNER3,SHIFT
*ELEMENT,TYPE=ACIN3d3,ELSET=ACIN
*ELGEN,ELSET=ACIN
*ELEMENT,TYPE=ASI4,ELSET=COUPLING
**1002,3,5,105,103
*ELGEN,ELSET=COUPLING
*SOLID SECTION,MATERIAL=AETHER,ELSET=ACIN
1.0,   0.0,0.0,0.0
*INTERFACE,ELSET=COUPLING
*NSET,NSET=RECRD3
*TRANSFORM,NSET=INNER3,TYPE=S
*end part
****************************************************************
*part,name=inf4
*NODE,NSET=POLE4
*NODE,NSET=INNER4
*NCOPY,OLD SET=INNER4,CHANGE NUMBER=100,NEW SET=INNER4,SHIFT
*NCOPY,OLD SET=INNER4,CHANGE NUMBER=200,NEW SET=INNER4,SHIFT
*ELEMENT, TYPE=ACIN3D4,ELSET=ACIN
*ELGEN, ELSET=ACIN
*ELEMENT,TYPE=ASI4,ELSET=COUPLING
*ELGEN,ELSET=COUPLING
*SOLID SECTION,MATERIAL=AETHER,ELSET=ACIN
1.0,   0.0,0.0,0.0
*INTERFACE,ELSET=COUPLING
*NSET,NSET=RECRD4
*TRANSFORM,NSET=INNER4,TYPE=S
*end part
****************************************************************
*part,name=inf6
*NODE,NSET=POLE6
*NODE,NSET=INNER6
*NCOPY,OLD SET=INNER6,CHANGE NUMBER=100,NEW SET=INNER6,SHIFT
*NCOPY,OLD SET=INNER6,CHANGE NUMBER=200,NEW SET=INNER6,SHIFT
*ELEMENT,TYPE=ACIN3d6,ELSET=ACIN
*ELGEN,ELSET=ACIN
*ELEMENT,TYPE=ASI8,ELSET=COUPLING
*ELGEN,ELSET=COUPLING
*SOLID SECTION,MATERIAL=AETHER,ELSET=ACIN
1.0,   0.0,0.0,0.0
*INTERFACE,ELSET=COUPLING
*NSET,NSET=RECRD6
*TRANSFORM,NSET=INNER6,TYPE=S
*end part
****************************************************************
*part,name=inf8
*NODE,NSET=POLE8
*NODE,NSET=INNER8
*NCOPY,OLD SET=INNER8,CHANGE NUMBER=100,NEW SET=INNER8,SHIFT
*NCOPY,OLD SET=INNER8,CHANGE NUMBER=200,NEW SET=INNER8,SHIFT
*ELEMENT, TYPE=ACIN3D8,ELSET=ACIN
*ELGEN, ELSET=ACIN
*ELEMENT,TYPE=ASI8,ELSET=COUPLING
*ELGEN,ELSET=COUPLING
*SOLID SECTION,MATERIAL=AETHER,ELSET=ACIN
1.0,   0.0,0.0,0.0
*INTERFACE,ELSET=COUPLING
**
*NSET,NSET=RECRD8
*TRANSFORM,NSET=INNER8,TYPE=S
*endpart
***********************************************
*assembly,name=whole
*instance,part=inf3,name=i_inf3
*endinstance
*instance,part=inf4,name=i_inf4
*endinstance
*instance,part=inf6,name=i_inf6
*endinstance
*instance,part=inf8,name=i_inf8
*endinstance
*nset,nset=recrd
i_inf3.recrd3,i_inf4.recrd4,i_inf6.recrd6,i_inf8.recrd8
*nset,nset=inner
i_inf3.inner3,i_inf4.inner4,i_inf6.inner6,i_inf8.inner8
*nset,nset=pole
i_inf3.pole3,i_inf4.inner4,i_inf6.pole6,i_inf8.pole8
*end assembly
***********************************************
**
** FLUID MATERIAL PROPERTIES; ARTIFICIAL MEDIUM FOR CONVENIENCE
** USING RHO = 1, BULKMOD = (2*PI)^2, WE HAVE C = 2*PI AND K = F.
**
*MATERIAL,NAME=AETHER
*DENSITY
1.000000,
*ACOUSTIC MEDIUM,BULK MODULUS
**
*STEP
*STEADY STATE DYNAMIC,DIRECT
**
** SURFACE RADIAL ACCELERATION OF 1/(RHO C) MAGNITUDE
**
*BOUNDARY,TYPE=VELOCITY
INNER,1,1,0.159154943
INNER,2,3,0.0
POLE,1,1,0.159154943
POLE,2,3,0.0
*NODE PRINT,NSET=RECRD
POR,PPOR
*OUTPUT,FIELD
*NODE OUTPUT
POR,u
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=RECRD
POR,
*END STEP