*HEADING
EXAMPLE PROBLEM: MODE 0 SPHERE USING BGT BOUNDARY CONDITION AND AC3D8R
***************************
*PARAMETER
f=1.0
T=4.0/f
T4=T/4.0
omega=2*3.14159*f
t0=-0.25/f
*PREPRINT,ECHO=YES,HISTORY=YES,MODEL=YES
***************************
*NODE,NSET=POLE
*NODE,NSET=INNER
*NODE,NSET=MID
*NODE,NSET=OUTER
*NSET, NSET=OUT
*NCOPY,OLD SET=INNER,CHANGE NUMBER=100,NEW SET=INNER,SHIFT
*NCOPY,OLD SET=MID,CHANGE NUMBER=100,SHIFT
*NCOPY,OLD SET=OUTER,CHANGE NUMBER=100,SHIFT
*ELEMENT, TYPE=AC3D8R,ELSET=LAYER_1
*ELGEN, ELSET=LAYER_1
*ELEMENT, TYPE=AC3D8R,ELSET=LAYER_2
*ELGEN, ELSET=LAYER_2
*ELSET, ELSET=EALL
LAYER_1, LAYER_2
*NCOPY, CHANGE NUMBER=1000, OLD SET=INNER, SHIFT, NEW SET=SHELL
*ELEMENT, TYPE=S4R, ELSET=SHELL
*ELGEN,ELSET=SHELL
*SURFACE,TYPE=ELEMENT,NAME=SLAVE
LAYER_1, S1
*SURFACE,TYPE=ELEMENT,NAME=MASTER
SHELL,SPOS
*TIE,NAME=COUPLING
SLAVE, MASTER
**
** 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
*SOLID SECTION,MATERIAL=AETHER,ELSET=EALL
**
*SHELL SECTION,MATERIAL=MAT1, ELSET=SHELL
0.01
*MATERIAL,NAME=MAT1
*ELASTIC
207E0,0.3
*DENSITY
7800,
*NSET,NSET=RECRD
*TRANSFORM,NSET=SHELL,TYPE=S
*AMPLITUDE,NAME=COSN,DEF=PERIODIC
1,<omega>,<t0>,0.0
1,0,0,0
*STEP
*DYNAMIC,EXPLICIT
,<T>
*IMPEDANCE,nonreflecting=spherical
LAYER_2, I2,1.15915
**
*BOUNDARY,TYPE=VELOCITY,AMP=COSN
SHELL, 1,1,0.159154943
*BOUNDARY
SHELL, 2,3,0.0
SHELL, 4,6,0.0
*OUTPUT,FIELD,NUMBER=5
*NODE OUTPUT,NSET=SHELL
V,
*NODE OUTPUT
POR
*OUTPUT,HISTORY,FREQ=1
*NODE OUTPUT, NSET=OUT
POR
**
** Output requests for qa testing
**
*NSET,NSET=QA_TEST1
SHELL,
*NSET,NSET=QA_TEST2,ELSET=EALL
*OUTPUT,FIELD,NUM=1
*NODE OUTPUT,NSET=QA_TEST1
U,
*NODE OUTPUT,NSET=QA_TEST2
POR,
*OUTPUT,HISTORY,VAR=PRESELECT,FREQ=1
*END STEP