*HEADING
: BONDED RUBBER CYLINDRICAL BLOCK: CAXA42 ELEMENTS
  AXIAL COMPRESSION + RADIAL EXCITATION
*RESTART,WRITE,FREQ=100
*NODE,NSET=NALL
*NGEN,NSET=BOT
*NGEN,NSET=TOP
*NSET,NSET=OUTWALL,GEN
*NSET,NSET=NOUT
OUTWALL,200
*NSET,NSET=CENT,GEN
*NGEN,NSET=TOP1
*NFILL,NSET=NALL
BOT,TOP,8,10
*ELEMENT,TYPE=CAXA4H2, OFFSET=1000
*ELGEN,ELSET=EALL
*ELSET,ELSET=MID,GEN
*SOLID SECTION,ELSET=EALL,MATERIAL=MOONEY
*MATERIAL,NAME=MOONEY
*HYPERELASTIC,N=1
80.0, 20.0
*DENSITY
1.E-4,
*NSET, NSET=ATIE, GEN
*NSET, NSET=BTIE, GEN
*NSET, NSET=CTIE, GEN
*EQUATION
BTIE, 1, 1, ATIE, 1, -1
CTIE, 1, 1, ATIE, 1, -1
BTIE, 2, 1, ATIE, 2, -1
CTIE, 2, 1, ATIE, 2, -1
**
** INTERFACE ELEMENT DEFINITIONS
**
*ELEMENT, TYPE=IRS21A
*ELGEN, ELSET=T0
*ELGEN, ELSET=T90
*ELGEN, ELSET=T180
*INTERFACE, ELSET=T0, NAME=ST0
*SURFACE BEHAVIOR, NO SEPARATION
*FRICTION, ROUGH
*ASYMMETRIC-AXISYMMETRIC, ANGLE=0.0, MODE=2
*INTERFACE, ELSET=T90, NAME=ST90
*FRICTION, ROUGH
*SURFACE BEHAVIOR, NO SEPARATION
*ASYMMETRIC-AXISYMMETRIC, ANGLE=90.0, MODE=2
*INTERFACE, ELSET=T180, NAME=ST180
*SURFACE BEHAVIOR, NO SEPARATION
*FRICTION, ROUGH
*ASYMMETRIC-AXISYMMETRIC, ANGLE=180.0, MODE=2
*ELSET, ELSET=RSURF
T0, T90, T180
*RIGID SURFACE, ELSET=T0, TYPE=SEGMENTS
START, 1.5, 1.029
LINE, -0.1, 1.029
*RIGID SURFACE, ELSET=T90, TYPE=SEGMENTS
START, 1.5, 1.029
LINE, -0.1, 1.029
*RIGID SURFACE, ELSET=T180, TYPE=SEGMENTS
START, 1.5, 1.029
LINE, -0.1, 1.029
*MPC
TIE, 10000, 50000
TIE, 20000, 50000
TIE, 30000, 50000
*NSET, NSET=BASE, GEN
*AMPLITUDE,VALUE=R,NAME=ZERO
 0.0000000,  0.0000000,  6.2831853,  0.0000000
*AMPLITUDE,VALUE=R,NAME=SINEWAVE
 0.0000000,  0.0000000,  0.0628319,  0.0627905,  0.1256637,  0.1253332,  0.1884956,  0.1873813
 0.2513274,  0.2486899,  0.3141593,  0.3090170,  0.3769911,  0.3681246,  0.4398230,  0.4257793
 0.5026548,  0.4817537,  0.5654867,  0.5358268,  0.6283185,  0.5877853,  0.6911504,  0.6374240
 0.7539822,  0.6845471,  0.8168141,  0.7289686,  0.8796459,  0.7705132,  0.9424778,  0.8090170
 1.0053096,  0.8443279,  1.0681415,  0.8763067,  1.1309734,  0.9048271,  1.1938052,  0.9297765
 1.2566371,  0.9510565,  1.3194689,  0.9685832,  1.3823008,  0.9822873,  1.4451326,  0.9921147
 1.5079645,  0.9980267,  1.5707963,  1.0000000,  1.6336282,  0.9980267,  1.6964600,  0.9921147
 1.7592919,  0.9822873,  1.8221237,  0.9685832,  1.8849556,  0.9510565,  1.9477874,  0.9297765
*BASELINE CORRECTION
*STEP,NLGEOM,INC=20
APPLY AXIAL COMPRESSION BY RIGID PLATEN
*STATIC
0.1 , 1.0
*BOUNDARY
BASE, 2
1, 1, 1
1001, 1, 1
2001, 1, 1
50000,1,1
50000,6,6
50000,2,2,-0.3
*NODE PRINT,FREQ=999
U,RF
*EL PRINT,ELSET=MID,POSITION=AVERAGED AT NODES,SUMMARY=NO,FREQ=999
S,
E,
*NODE FILE,NSET=NOUT,FREQ=999
U,RF
*OUTPUT,FIELD, FREQ=999
*NODE OUTPUT
U,RF
*OUTPUT,HISTORY,FREQ=999
*NODE OUTPUT,NSET=NOUT
U,RF
*END STEP
*STEP,NLGEOM
 FREQUENCY EXTRACTION FOR COMPRESSED CYLINDER
*FREQUENCY, eigensolver=lanczos
5, , ,-10000.
*BOUNDARY,BASENAME=RSNODE
50000,2
*EL PRINT,FREQ=0
*NODE PRINT
U,
RF,
*NODE FILE, FREQ=0
U,
RF,
*MODAL FILE
*OUTPUT,FIELD
*NODE OUTPUT, NSET=NOUT
U,
RF,
*OUTPUT,HISTORY,FREQ=0
*NODE OUTPUT, NSET=NALL
U,
RF,
*OUTPUT,FIELD
*MODAL OUTPUT
*OUTPUT,HISTORY
*MODAL OUTPUT
*END STEP
*STEP
RESPONSE TO RADIAL EXCITATION AT OUTSIDE WALL
*MODAL DYNAMIC
0.01,6.28
*BASE MOTION,AMP=ZERO,DOF=2
*BASE MOTION,AMP=SINEWAVE,DOF=2,BASENAME=RSNODE
*MODAL DAMPING
*SELECT EIGENMODES, GENERATE
*MODAL PRINT,FREQ=100
GU,
BM,
*NODE PRINT,NSET=NOUT,FREQ=999
U,TU
V,TV
A,TA
RF,
*NODE FILE,NSET=NOUT,FREQ=100
U,TU
*NODE FILE,NSET=NOUT,FREQ=400
V,
TV,
A,
TA,
RF,
*MODAL FILE,FREQ=400
GU,
BM,
*OUTPUT,FIELD, FREQ=100
*NODE OUTPUT
U,TU
*OUTPUT,HISTORY,FREQ=100
*NODE OUTPUT,NSET=NOUT
U,TU
*OUTPUT,FIELD, FREQ=400
*NODE OUTPUT
V,
TV,
A,
TA,
RF,
*MODAL OUTPUT
GU,
BM,
*OUTPUT,HISTORY,FREQ=400
*NODE OUTPUT,NSET=NOUT
V,
TV,
A,
TA,
RF,
*MODAL OUTPUT
GU,
BM,
*END STEP