*HEADING
BENCHMARK PROBLEM: MODE 0 SPHERE USING ACIN3D3, ACIN3D4, ACIN3D6, ACIN3D8.
*PREPRINT,MODEL=YES
*part,name=inf3
*NODE,NSET=POLE3
      1,           0.,           1.,0.,           0.,           1.,0.
     31,           0.,          -1.,0.,           0.,          -1.,0.
     32,           0.,      1.07958,0.,           0.,      1.07958,0.
     62,           0.,     -1.07958,0.,           0.,     -1.07958,0.
     63,           0.,      1.15915,0.,           0.,      1.15915,0.
     93,           0.,     -1.15915,0.,           0.,     -1.15915,0.
*NODE,NSET=INNER3
      2,     0.104528,     0.994522,0.,     0.104528,     0.994522,0.
      3,     0.207912,     0.978148,0.,     0.207912,     0.978148,0.
      4,     0.309017,     0.951057,0.,     0.309017,     0.951057,0.
      5,     0.406737,     0.913545,0.,     0.406737,     0.913545,0.
      6,          0.5,     0.866025,0.,          0.5,     0.866025,0.
      7,     0.587785,     0.809017,0.,     0.587785,     0.809017,0.
      8,     0.669131,     0.743145,0.,     0.669131,     0.743145,0.
      9,     0.743145,     0.669131,0.,     0.743145,     0.669131,0.
     10,     0.809017,     0.587785,0.,     0.809017,     0.587785,0.
     11,     0.866025,          0.5,0.,     0.866025,          0.5,0.
     12,     0.913545,     0.406737,0.,     0.913545,     0.406737,0.
     13,     0.951057,     0.309017,0.,     0.951057,     0.309017,0.
     14,     0.978148,     0.207912,0.,     0.978148,     0.207912,0.
     15,     0.994522,     0.104528,0.,     0.994522,     0.104528,0.
     16,           1.,     0.0,     0.,           1.,     0.0,     0.
     17,     0.994522,    -0.104528,0.,     0.994522,    -0.104528,0.
     18,     0.978148,    -0.207912,0.,     0.978148,    -0.207912,0.
     19,     0.951057,    -0.309017,0.,     0.951057,    -0.309017,0.
     20,     0.913545,    -0.406737,0.,     0.913545,    -0.406737,0.
     21,     0.866025,         -0.5,0.,     0.866025,         -0.5,0.
     22,     0.809017,    -0.587785,0.,     0.809017,    -0.587785,0.
     23,     0.743145,    -0.669131,0.,     0.743145,    -0.669131,0.
     24,     0.669131,    -0.743145,0.,     0.669131,    -0.743145,0.
     25,     0.587785,    -0.809017,0.,     0.587785,    -0.809017,0.
     26,          0.5,    -0.866025,0.,          0.5,    -0.866025,0.
     27,     0.406737,    -0.913545,0.,     0.406737,    -0.913545,0.
     28,     0.309017,    -0.951057,0.,     0.309017,    -0.951057,0.
     29,     0.207912,    -0.978148,0.,     0.207912,    -0.978148,0.
     30,     0.104528,    -0.994522,0.,     0.104528,    -0.994522,0.
*NCOPY,OLD SET=INNER3,CHANGE NUMBER=100,NEW SET=INNER3,SHIFT
0.,0.,0.
0.,0.,0.,   0.,1.,0., 22.0
*NCOPY,OLD SET=INNER3,CHANGE NUMBER=200,NEW SET=INNER3,SHIFT
0.,0.,0.
0.,0.,0.,   0.,1.,0., 11.0
*ELEMENT,TYPE=ACIN3d3,ELSET=ACIN
1, 3,103,1 
42, 129,29,31
2, 3,105,103
3, 3,5,105
*ELGEN,ELSET=ACIN
3,13,2,2
2,13,2,2
*ELEMENT,TYPE=ASI4,ELSET=COUPLING
1001, 3,103,1,3
**1002,3,5,105,103
1002, 3,105,103,103
1003, 3,5,105,105
1042, 129,29,31,129
*ELGEN,ELSET=COUPLING
1002,13,2,2
1003,13,2,2
*SOLID SECTION,MATERIAL=AETHER,ELSET=ACIN
1.0,   0.0,0.0,0.0
*INTERFACE,ELSET=COUPLING
*NSET,NSET=RECRD3
13,15,115,113,
*TRANSFORM,NSET=INNER3,TYPE=S
0.,0.,0.,  0.,0.,1.
*end part
****************************************************************
*part,name=inf4
*NODE,NSET=POLE4
      1,           0.,           1.,0.,           0.,           1.,0.
     31,           0.,          -1.,0.,           0.,          -1.,0.
     32,           0.,      1.07958,0.,           0.,      1.07958,0.
     62,           0.,     -1.07958,0.,           0.,     -1.07958,0.
     63,           0.,      1.15915,0.,           0.,      1.15915,0.
     93,           0.,     -1.15915,0.,           0.,     -1.15915,0.
*NODE,NSET=INNER4
      2,     0.104528,     0.994522,0.,     0.104528,     0.994522,0.
      3,     0.207912,     0.978148,0.,     0.207912,     0.978148,0.
      4,     0.309017,     0.951057,0.,     0.309017,     0.951057,0.
      5,     0.406737,     0.913545,0.,     0.406737,     0.913545,0.
      6,          0.5,     0.866025,0.,          0.5,     0.866025,0.
      7,     0.587785,     0.809017,0.,     0.587785,     0.809017,0.
      8,     0.669131,     0.743145,0.,     0.669131,     0.743145,0.
      9,     0.743145,     0.669131,0.,     0.743145,     0.669131,0.
     10,     0.809017,     0.587785,0.,     0.809017,     0.587785,0.
     11,     0.866025,          0.5,0.,     0.866025,          0.5,0.
     12,     0.913545,     0.406737,0.,     0.913545,     0.406737,0.
     13,     0.951057,     0.309017,0.,     0.951057,     0.309017,0.
     14,     0.978148,     0.207912,0.,     0.978148,     0.207912,0.
     15,     0.994522,     0.104528,0.,     0.994522,     0.104528,0.
     16,           1.,     0.0,     0.,           1.,     0.0,     0.
     17,     0.994522,    -0.104528,0.,     0.994522,    -0.104528,0.
     18,     0.978148,    -0.207912,0.,     0.978148,    -0.207912,0.
     19,     0.951057,    -0.309017,0.,     0.951057,    -0.309017,0.
     20,     0.913545,    -0.406737,0.,     0.913545,    -0.406737,0.
     21,     0.866025,         -0.5,0.,     0.866025,         -0.5,0.
     22,     0.809017,    -0.587785,0.,     0.809017,    -0.587785,0.
     23,     0.743145,    -0.669131,0.,     0.743145,    -0.669131,0.
     24,     0.669131,    -0.743145,0.,     0.669131,    -0.743145,0.
     25,     0.587785,    -0.809017,0.,     0.587785,    -0.809017,0.
     26,          0.5,    -0.866025,0.,          0.5,    -0.866025,0.
     27,     0.406737,    -0.913545,0.,     0.406737,    -0.913545,0.
     28,     0.309017,    -0.951057,0.,     0.309017,    -0.951057,0.
     29,     0.207912,    -0.978148,0.,     0.207912,    -0.978148,0.
     30,     0.104528,    -0.994522,0.,     0.104528,    -0.994522,0.
*NCOPY,OLD SET=INNER4,CHANGE NUMBER=100,NEW SET=INNER4,SHIFT
0.,0.,0.
0.,0.,0.,   0.,1.,0., 11.0
*NCOPY,OLD SET=INNER4,CHANGE NUMBER=200,NEW SET=INNER4,SHIFT
0.,0.,0.
0.,0.,0.,   0.,1.,0., 5.5
*ELEMENT, TYPE=ACIN3D4,ELSET=ACIN
3,   3,5,105,103
*ELGEN, ELSET=ACIN
3,13,2,1
*ELEMENT,TYPE=ASI4,ELSET=COUPLING
1003, 3,5,105,103
*ELGEN,ELSET=COUPLING
1003,13,2,1
*SOLID SECTION,MATERIAL=AETHER,ELSET=ACIN
1.0,   0.0,0.0,0.0
*INTERFACE,ELSET=COUPLING
*NSET,NSET=RECRD4
3,5,11,16,21,26
*TRANSFORM,NSET=INNER4,TYPE=S
0.,0.,0.,  0.,0.,1.
*end part
****************************************************************
*part,name=inf6
*NODE,NSET=POLE6
      1,           0.,           1.,0.,           0.,           1.,0.
     31,           0.,          -1.,0.,           0.,          -1.,0.
     32,           0.,      1.07958,0.,           0.,      1.07958,0.
     62,           0.,     -1.07958,0.,           0.,     -1.07958,0.
     63,           0.,      1.15915,0.,           0.,      1.15915,0.
     93,           0.,     -1.15915,0.,           0.,     -1.15915,0.
*NODE,NSET=INNER6
      2,     0.104528,     0.994522,0.,     0.104528,     0.994522,0.
      3,     0.207912,     0.978148,0.,     0.207912,     0.978148,0.
      4,     0.309017,     0.951057,0.,     0.309017,     0.951057,0.
      5,     0.406737,     0.913545,0.,     0.406737,     0.913545,0.
      6,          0.5,     0.866025,0.,          0.5,     0.866025,0.
      7,     0.587785,     0.809017,0.,     0.587785,     0.809017,0.
      8,     0.669131,     0.743145,0.,     0.669131,     0.743145,0.
      9,     0.743145,     0.669131,0.,     0.743145,     0.669131,0.
     10,     0.809017,     0.587785,0.,     0.809017,     0.587785,0.
     11,     0.866025,          0.5,0.,     0.866025,          0.5,0.
     12,     0.913545,     0.406737,0.,     0.913545,     0.406737,0.
     13,     0.951057,     0.309017,0.,     0.951057,     0.309017,0.
     14,     0.978148,     0.207912,0.,     0.978148,     0.207912,0.
     15,     0.994522,     0.104528,0.,     0.994522,     0.104528,0.
     16,           1.,     0.0,     0.,           1.,     0.0,     0.
     17,     0.994522,    -0.104528,0.,     0.994522,    -0.104528,0.
     18,     0.978148,    -0.207912,0.,     0.978148,    -0.207912,0.
     19,     0.951057,    -0.309017,0.,     0.951057,    -0.309017,0.
     20,     0.913545,    -0.406737,0.,     0.913545,    -0.406737,0.
     21,     0.866025,         -0.5,0.,     0.866025,         -0.5,0.
     22,     0.809017,    -0.587785,0.,     0.809017,    -0.587785,0.
     23,     0.743145,    -0.669131,0.,     0.743145,    -0.669131,0.
     24,     0.669131,    -0.743145,0.,     0.669131,    -0.743145,0.
     25,     0.587785,    -0.809017,0.,     0.587785,    -0.809017,0.
     26,          0.5,    -0.866025,0.,          0.5,    -0.866025,0.
     27,     0.406737,    -0.913545,0.,     0.406737,    -0.913545,0.
     28,     0.309017,    -0.951057,0.,     0.309017,    -0.951057,0.
     29,     0.207912,    -0.978148,0.,     0.207912,    -0.978148,0.
     30,     0.104528,    -0.994522,0.,     0.104528,    -0.994522,0.
*NCOPY,OLD SET=INNER6,CHANGE NUMBER=100,NEW SET=INNER6,SHIFT
0.,0.,0.
0.,0.,0.,   0.,1.,0., 16.0
*NCOPY,OLD SET=INNER6,CHANGE NUMBER=200,NEW SET=INNER6,SHIFT
0.,0.,0.
0.,0.,0.,   0.,1.,0., 8.0
*ELEMENT,TYPE=ACIN3d6,ELSET=ACIN
1, 3,103,1 ,  203, 102,2
42, 129,29,31,  229,30,130
2, 3,105,103, 204,104,203
3, 3,5,105, 4,205,204
*ELGEN,ELSET=ACIN
3,13,2,2
2,13,2,2
*ELEMENT,TYPE=ASI8,ELSET=COUPLING
1001, 3,103,1,1 ,  203, 102, 1,2
1002, 3,105,103,103, 204,104,103,203
1003, 3,5,105,105, 4,205,105,204
1042, 129,29,31,31,  229,30,31,130
*ELGEN,ELSET=COUPLING
1003,13,2,2
1002,13,2,2
*SOLID SECTION,MATERIAL=AETHER,ELSET=ACIN
1.0,   0.0,0.0,0.0
*INTERFACE,ELSET=COUPLING
*NSET,NSET=RECRD6
213,214,215
*TRANSFORM,NSET=INNER6,TYPE=S
0.,0.,0.,  0.,0.,1.
*end part
****************************************************************
*part,name=inf8
*NODE,NSET=POLE8
      1,           0.,           1.,0.,           0.,           1.,0.
     31,           0.,          -1.,0.,           0.,          -1.,0.
     32,           0.,      1.07958,0.,           0.,      1.07958,0.
     62,           0.,     -1.07958,0.,           0.,     -1.07958,0.
     63,           0.,      1.15915,0.,           0.,      1.15915,0.
     93,           0.,     -1.15915,0.,           0.,     -1.15915,0.
*NODE,NSET=INNER8
      2,     0.104528,     0.994522,0.,     0.104528,     0.994522,0.
      3,     0.207912,     0.978148,0.,     0.207912,     0.978148,0.
      4,     0.309017,     0.951057,0.,     0.309017,     0.951057,0.
      5,     0.406737,     0.913545,0.,     0.406737,     0.913545,0.
      6,          0.5,     0.866025,0.,          0.5,     0.866025,0.
      7,     0.587785,     0.809017,0.,     0.587785,     0.809017,0.
      8,     0.669131,     0.743145,0.,     0.669131,     0.743145,0.
      9,     0.743145,     0.669131,0.,     0.743145,     0.669131,0.
     10,     0.809017,     0.587785,0.,     0.809017,     0.587785,0.
     11,     0.866025,          0.5,0.,     0.866025,          0.5,0.
     12,     0.913545,     0.406737,0.,     0.913545,     0.406737,0.
     13,     0.951057,     0.309017,0.,     0.951057,     0.309017,0.
     14,     0.978148,     0.207912,0.,     0.978148,     0.207912,0.
     15,     0.994522,     0.104528,0.,     0.994522,     0.104528,0.
     16,           1.,     0.0,     0.,           1.,     0.0,     0.
     17,     0.994522,    -0.104528,0.,     0.994522,    -0.104528,0.
     18,     0.978148,    -0.207912,0.,     0.978148,    -0.207912,0.
     19,     0.951057,    -0.309017,0.,     0.951057,    -0.309017,0.
     20,     0.913545,    -0.406737,0.,     0.913545,    -0.406737,0.
     21,     0.866025,         -0.5,0.,     0.866025,         -0.5,0.
     22,     0.809017,    -0.587785,0.,     0.809017,    -0.587785,0.
     23,     0.743145,    -0.669131,0.,     0.743145,    -0.669131,0.
     24,     0.669131,    -0.743145,0.,     0.669131,    -0.743145,0.
     25,     0.587785,    -0.809017,0.,     0.587785,    -0.809017,0.
     26,          0.5,    -0.866025,0.,          0.5,    -0.866025,0.
     27,     0.406737,    -0.913545,0.,     0.406737,    -0.913545,0.
     28,     0.309017,    -0.951057,0.,     0.309017,    -0.951057,0.
     29,     0.207912,    -0.978148,0.,     0.207912,    -0.978148,0.
     30,     0.104528,    -0.994522,0.,     0.104528,    -0.994522,0.
*NCOPY,OLD SET=INNER8,CHANGE NUMBER=100,NEW SET=INNER8,SHIFT
0.,0.,0.
0.,0.,0.,   0.,1.,0., 11.0
*NCOPY,OLD SET=INNER8,CHANGE NUMBER=200,NEW SET=INNER8,SHIFT
0.,0.,0.
0.,0.,0.,   0.,1.,0., 5.5
*ELEMENT, TYPE=ACIN3D8,ELSET=ACIN
3, 3,5,105,103, 4,205,104,203
*ELGEN, ELSET=ACIN
3,13,2,1
*ELEMENT,TYPE=ASI8,ELSET=COUPLING
1003, 3,5,105,103, 4,205,104,203
*ELGEN,ELSET=COUPLING
1003,13,2,1
*SOLID SECTION,MATERIAL=AETHER,ELSET=ACIN
1.0,   0.0,0.0,0.0
*INTERFACE,ELSET=COUPLING
**
*NSET,NSET=RECRD8
213,215
*TRANSFORM,NSET=INNER8,TYPE=S
0.,0.,0.,  0.,0.,1.
*endpart
***********************************************
*assembly,name=whole
*instance,part=inf3,name=i_inf3
0,0,0
0,-1,0,  0,1,0,  0
*endinstance
*instance,part=inf4,name=i_inf4
0,0,0
0,-1,0,  0,1,0,  90
*endinstance
*instance,part=inf6,name=i_inf6
0,0,0
0,-1,0,  0,1,0,  180
*endinstance
*instance,part=inf8,name=i_inf8
0,0,0
0,-1,0,  0,1,0,  270
*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
39.4784176, 2
**
*STEP
*STEADY STATE DYNAMIC,DIRECT
00.01, 
0.05, 
0.1, 
0.13, 
0.16, 
0.2, 
0.3, 
0.4, 
0.5, 
0.6, 
0.7, 
0.8, 
0.9, 
1.0, 
1.1, 
1.2, 
1.3, 
1.4, 
1.5, 
1.6, 
1.7, 
1.8, 
1.9, 
2.0, 
2.3, 
2.5, 
3.0, 
4.0, 
5.0, 
6.0, 
7.0, 
8.0, 
9.0, 
10.0, 
12.0, 
15.0, 
20.0, 
**
** 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