*HEADING
EXAMPLE PROBLEM: MODE 0 SPHERE USING BGT BOUNDARY CONDITION AND ACAX8. AUTOMATIC BGT.
*NODE,NSET=INNER
*NODE,NSET=MID
*NODE,NSET=OUTER
*ELEMENT, TYPE=ACAX8,ELSET=EALL
*ELGEN, ELSET=EALL
*SURFACE,NAME=RADIATING
EALL,S3
*ELEMENT,TYPE=ASI3A,ELSET=COUPLING
*ELGEN,ELSET=COUPLING
**
** 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
*INTERFACE,ELSET=COUPLING
**
** BAYLISS, GUNZBERGER, TURKEL FIRST ORDER CONDITION. EXACT
** FOR MODE 0 OF A SPHERE.
** SEE SIAM JAM 1982, V42 N2 P430FF.
**
** 1/C1 = 1/SQRT(K * RHO) = 1/2/PI
** 1/K1 =-1/(R *RHO * (2 PI F)^2 )
**
*IMPEDANCEPROPERTY, NAME=BGT_AUTO,TYPE=SPHERE
*IMPEDANCE PROPERTY,NAME=BGT_TABLE
*NSET,NSET=QA_TEST
*TRANSFORM,NSET=INNER,TYPE=S
*STEP
*STEADY STATE DYNAMIC,DIRECT
*SIMPEDANCE
RADIATING, BGT_AUTO
**
** SURFACE RADIAL ACCELERATION OF 1/(RHO C) MAGNITUDE
**
*BOUNDARY,TYPE=VELOCITY
INNER,1,1,0.159154943
INNER,2,2,0.0
*NODE PRINT,NSET=QA_TEST,FREQUENCY=10000
POR,v
*NODE FILE,NSET=QA_TEST,FREQUENCY=10000
POR,v
*OUTPUT,FIELD,FREQUENCY=10000
*NODE OUTPUT
POR,v
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=QA_TEST
POR,v
*END STEP