*HEADING
**
*PREPRINT, ECHO=NO, MODEL=NO, HISTORY=NO
*NODE,NSET=BEAMS,SYSTEM=C
********************************
*NODE,NSET=STIFFNER
*****
*NGEN,NSET=BEAMS,LINE=C
*NSET,NSET=OUT
*NSET,NSET=P_OUT
****************************
*NODE,NSET=INNER,SYSTEM=C
*NGEN,NSET=INNER,LINE=C
*****************************
*NODE,NSET=OUTER,SYSTEM=C
*NGEN,NSET=OUTER,LINE=C
*****************************
** USE "*nmap" TO CREATE ELLIPTICAL
** OUTER BOUNDARY
*NMAP,NSET=OUTER,TYPE=RECTANGULAR
** points a, b, c: see documentation
** stretch factors for x, y, z:
*****************************
*NFILL,NSET=ACOUSTIC
INNER,OUTER,20,100
**********************************
*ELEMENT,TYPE=B21,ELSET=BEAMS
*ELGEN,ELSET=BEAMS
*ELEMENT,TYPE=B21,ELSET=STIFFNER
******************
*ELEMENT, TYPE=AC2D4,ELSET=ACOUSTIC
*ELGEN, ELSET=ACOUSTIC
****************************
*ELGEN,ELSET=SLAVE
*SURFACE, TYPE=ELEMENT, NAME=SLAVE
SLAVE,S1
*SURFACE,TYPE=ELEMENT,NAME=MASTER
BEAMS,SPOS
*TIE,NAME=COUPLING
SLAVE,MASTER
************************************
*MATERIAL,NAME=WATER
*DENSITY
1000.0,
*ACOUSTIC MEDIUM,BULK MODULUS
*SOLID SECTION,MATERIAL=WATER,ELSET=ACOUSTIC
0.1,
*******************
*BEAM SECTION,MATERIAL=MAT1,ELSET=BEAMS,SECTION=RECT
0.1,0.1
*BEAM SECTION,MATERIAL=MAT1,ELSET=STIFFNER,SECTION=RECT
0.1,0.1
*MATERIAL,NAME=MAT1
*ELASTIC
207E9,0.3
*DENSITY
7800
*************************
*ELSET,ELSET=OUTER,GENERATE
*SURFACE,NAME=OUTER
OUTER,S3
**********************
*STEP,INC=1000
*STEADY STATE DYNAMIC,DIRECT
*IMPEDANCE,NON=E
OUTER,I3,   9.6,0.5527708,0,0,0,0,1,0
*CLOAD
55,2,100.0
*BOUNDARY
25001,1,1,0.0
25071,1,1,0.0
25001,6,6,0.0
25071,6,6,0.0
55,1,1,0.0
55,6,6,0.0
*ENERGYPRINT
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=BEAMS
U,
*NODE OUTPUT,NSET=ACOUSTIC
POR,
*ELEMENTOUTPUT,ELSET=ACOUSTIC
ENER,ELEN,INTEN
*OUTPUT,HISTORY,FREQ=1
*NODE OUTPUT,NSET=BEAMS
U,V,A
*END STEP
*STEP,INC=1000
*STEADY STATE DYNAMIC,DIRECT
*IMPEDANCE,NONREF=E,OP=NEW
*SIMPEDANCE,OP=NEW,NONREF=E
OUTER,9.6,0.5527708,0,0,0,0,1,0
*CLOAD
55,2,100.0
*BOUNDARY
25001,1,1,0.0
25071,1,1,0.0
25001,6,6,0.0
25071,6,6,0.0
55,1,1,0.0
55,6,6,0.0
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=BEAMS
U,
*NODE OUTPUT,NSET=ACOUSTIC
POR,
*ELEMENT OUTPUT, ELSET=ACOUSTIC
ENER,ELEN,INTEN
*OUTPUT,HISTORY,FREQ=1
*NODE OUTPUT,NSET=BEAMS
U,V,A
*END STEP
****************************
*STEP,INC=1000
*Frequency, eigensolver=Lanczos, acoustic coupling=on, normalization=MASS
100, , 3000., -1.0, ,
*IMPEDANCE,NONREFLECTING=ELLI,OP=NEW
*SIMPEDANCE,NONREFLECTING=ELLI,OP=NEW
OUTER,9.6,0.5527708,0,0,0,0,1,0
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=BEAMS
U
*NODE OUTPUT,NSET=ACOUSTIC
POR,
*END STEP
************************