*HEADING
ADAPTIVE MESHING EXAMPLE
WATER SLOSHING IN A BAFFLED TANK
Unit lb, in, sec
*RESTART, WRITE, NUMBER=30
*NODE
1,   10.00, 0.0, 0.0
21,  20.00, 0.0, 0.0
127, 10.00, 0.0, 1.2
147, 20.00, 0.0, 1.2
316, 10.00, 0.0, 4.0
336, 20.00, 0.0, 4.0
5041,10.00, 6.0, 0.0 
5061,20.00, 6.0, 0.0
5167,10.00, 6.0, 1.2
5187,20.00, 6.0, 1.2
5356,10.00, 6.0, 4.0 
5376,20.00, 6.0, 4.0 
75000,  25.00, 20.0
*NGEN, NSET=A1
1, 127,21
*NGEN, NSET=A2
5041, 5167,21
*NFILL, NSET=A
A1,A2,15,336
*NGEN, NSET=B1
21,147,21
*NGEN, NSET=B2
5061,5187,21
*NFILL, NSET=B
B1,B2,15,336
*NFILL,NSET=RIGHT1,BIAS=.9
A, B, 20, 1
*NGEN, NSET=C1
127, 316,21
*NGEN, NSET=C2
5167, 5356,21
*NFILL, NSET=C
C1,C2,15,336
*NGEN, NSET=D1
147,336,21
*NGEN, NSET=D2
5187,5376,21
*NFILL, NSET=D
D1,D2,15,336
*NFILL,NSET=RIGHT2,BIAS=.9
C, D, 20, 1
*NSET,NSET=RIGHT
RIGHT1,RIGHT2
*ELEMENT, TYPE=C3D8R
1,1,2,338,337,22,23,359,358
*ELGEN, ELSET=WATER1
1,20,1,1,15,21,20,15,336,300,
*****************************************
*NCOPY,CHANGE NUMBER=5376,OLD SET=RIGHT,
NEW SET=LEFT,REFLECT=MIRROR
10.,0., 20.,10.,3., 20.
10.,0., -5.,
*ELCOPY,ELEMENT SHIFT=4500,OLD SET=WATER1,
SHIFT NODES=5376,REFLECT,NEW SET=WATER2
*ELSET,ELSET=WATER
WATER1,WATER2
*NODE
10753, 0.0,0.0,0.0
10754,20.0,0.0,0.0
10755,20.0,0.0,6.0,
10756,0.00,0.0,6.0,
10757, 0.0,6.0,0.0,
10758,20.0,6.0,0.0,
10759,20.0,6.0,6.0,
10760,0.00,6.0,6.0,
*ELEMENT, TYPE=R3D4, ELSET=RIGID
9001,10756,10755,10754,10753
9002,10757,10758,10759,10760
9003,10755,10759,10758,10754
9004,10753,10757,10760,10756
9005,10753,10754,10758,10757
*NODE
10761,10.00,-0.1,1.2
10776,10.00,6.1,1.2
11001,10.00,-0.1,6.0
11016,10.00,6.1,6.0
*NGEN, NSET=S1
10761,10776,1
*NGEN, NSET=S2
11001,11016,1
*NFILL,NSET=S
S1,S2,15, 16
*NSET,NSET=EDGES, GEN
10777,11001,16
10792,11016,16
11001,11016,1
*NSET,NSET=OUT
316,5356,5692,10732
*ELEMENT, TYPE=S4R
9006,10761,10777,10778,10762
*ELGEN, ELSET=BAFFLE
9006,15,1,1,15,16,15
*INCLUDE, INPUT=ale_water_sloshingel.inp
*SHELL SECTION, ELSET=BAFFLE, MATERIAL=BAFFLE
.12,
*MATERIAL, NAME=BAFFLE
*HYPERELASTIC,POLYNOMIAL, N=1
100., 25.
*DENSITY
1.07E-3, 
*SOLID SECTION, ELSET=WATER, MATERIAL=WATER,
CONTROLS=SECT
*SECTION CONTROLS, NAME=SECT,
KINEMATICS=ORTHOGONAL, HOURGLASS=STIFFNESS
10000.,
***K=300,K/G=1.E+6, (Actual K of water = 300000)
*MATERIAL, NAME=WATER
*EOS, TYPE=USUP
1805.7878,0.,0.
*VISCOSITY
1.5E-8,
*DENSITY
.92E-4, 
*BOUNDARY
75000,3,6,0.0
EDGES,3,6,0.0
*AMPLITUDE,NAME=AMP,DEFINITION=SMOOTH STEP
0.,0.,.25,1.76777,.5,2.5,.75,1.76777,
1.0,0.0,10.0,0.0
*INITIAL CONDITIONS, TYPE=STRESS, GEOSTATIC
WATER,0.,4.,-.1420,0.,1.,1.,
*SURFACE,TYPE=ELEMENT, NAME=BAFFLESURF1,NOTHICK
BAFFLE,SNEG
*SURFACE,TYPE=ELEMENT, NAME=BAFFLESURF2,NOTHICK
BAFFLE,SPOS
*SURFACE,TYPE=ELEMENT, NAME=RSURF
RIGID,SPOS
*SURFACE,TYPE=ELEMENT, NAME=WATSURF1
WATER,
*SURFACE,TYPE=ELEMENT, NAME=WATSURF2
OUT1,S6
OUT11,S2
*SURFACE,TYPE=ELEMENT, NAME=WATSURF3
OUT2,S3
OUT21,S2
*RIGID BODY, REF NODE=75000,
POSITION=CENTER OF MASS,ELSET=RIGID
*STEP
*DYNAMIC, EXPLICIT
,2.0
*DLOAD
WATER, GRAV, 386., 0., 0.,-1.
*BOUNDARY,TYPE=VELOCITY,AMPLITUDE=AMP
75000,1,2,1.
EDGES,1,2,1.
*CONTACT PAIR,INTERACTION=WATRIG
RSURF,WATSURF1
*CONTACT PAIR,INTERACTION=WATBAFFLE,WEIGHT=1.0,
MECHANICAL CONSTRAINT=PENALTY
BAFFLESURF1,WATSURF2
*CONTACT PAIR,INTERACTION=WATBAFFLE,WEIGHT=1.0,
MECHANICAL CONSTRAINT=PENALTY
BAFFLESURF2,WATSURF3
*SURFACE INTERACTION, NAME=WATRIG
*SURFACE INTERACTION, NAME=WATBAFFLE
*FILE OUTPUT, NUMBER=1, TIMEMARKS=YES
*EL FILE,ELSET=OUT
 LE,
*NODE FILE,NSET=EDGES
U,
***OUTPUT,HISTORY,TIME INTERVAL=1.0E-2
***NODE OUTPUT,NSET=EDGES
**U,V
***NODE OUTPUT,NSET=OUT
**U,V
***ELEMENT OUTPUT,ELSET=OUT
**S,
**PRESS,
**SHRMOD,DILMOD
***ENERGY OUTPUT
**ALLKE,ALLIE,ALLAE,ALLVD,ALLWK,ETOTAL,
**DT,
*****
*NSET, NSET=QA_TEST_REFN
7500,
*ELSET, ELSET=QA_TEST
OUT,
*NSET, NSET=QA_TEST
OUT,
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTERVAL=1
*ELEMENT OUTPUT, ELSET=QA_TEST
LE,
*NODE OUTPUT, NSET=QA_TEST
U,
*NODE OUTPUT, NSET=QA_TEST_REFN
RF,
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTERVAL=5
*ELEMENT OUTPUT
LE,
*NODE OUTPUT
U,RF
*OUTPUT, HIST, FREQ=9999
*ENERGY OUTPUT, VAR=PRESELECT
******
*ADAPTIVE MESH,ELSET=WATER,CONTROLS=CONTROLS,
FREQUENCY=5,MESH SWEEPS=3
*ADAPTIVE MESH CONTROLS,NAME=CONTROLS,
SMOOTHING=GRADED,TRANSITION FEATURE ANGLE=0.0
*END STEP