*HEADING
EXAMPLE PROBLEM: MODE 0 SPHERE USING BGT BOUNDARY CONDITION AND ACAX4R.
***************************
*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=INNER
*NODE,NSET=MID
*NODE,NSET=OUTER
*NSET, NSET=OUT
*NCOPY, CHANGE NUMBER=100, OLD SET=INNER, SHIFT, NEW SET=BEAMS
*ELEMENT, TYPE=ACAX4R, ELSET=LAYER_1
*ELEMENT, TYPE=ACAX4R, ELSET=LAYER_2
*ELGEN, ELSET=LAYER_1
*ELGEN, ELSET=LAYER_2
*ELSET, ELSET=EALL
LAYER_1, LAYER_2
*ELEMENT, TYPE=B21,ELSET=BEAMS
*ELGEN,ELSET=BEAMS
*SURFACE,TYPE=ELEMENT,NAME=SLAVE
LAYER_1, S1
*SURFACE,TYPE=ELEMENT,NAME=MASTER
BEAMS,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
**
*BEAM SECTION,MATERIAL=MAT1,ELSET=BEAMS,SECTION=RECT
0.1,0.1
*MATERIAL,NAME=MAT1
*ELASTIC
207E0,0.3
*DENSITY
7800,
*NSET,NSET=RECRD
*TRANSFORM,NSET=BEAMS,TYPE=S
*AMPLITUDE,NAME=COSN,DEF=PERIODIC
1,<omega>,<t0>,0.0
1,0,0,0
*STEP
*DYNAMIC,EXPLICIT
,<T>
*IMPEDANCE,non=SPHER
LAYER_2, I3,1.15915
**
*BOUNDARY,TYPE=VELOCITY,AMP=COSN
BEAMS, 1,1,0.159154943
*BOUNDARY
BEAMS, 2,2,0.0
BEAMS, 6,6,0.0
*OUTPUT,FIELD,NUMBER=5
*NODE OUTPUT,NSET=BEAMS
V,
*NODE OUTPUT
POR
*OUTPUT,HISTORY,FREQ=1
*NODE OUTPUT, NSET=OUT
POR
**
** Output requests for qa testing
**
*NSET,NSET=QA_TEST1
BEAMS,
*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