*HEADING
**
 FOAM PROBLEM -- 3D
 TEST OF COHESIVE BEHAVIOR WITH DAMAGE
*PREPRINT, MODEL=YES,HISTORY=YES
*NODE 
101,-.6,0.05,-.0225
109,-.4,0.05,-.0225
501,-.6,.1,-.0225
509,-.4,.1,-.0225
4501,-.6,0.05,.1775
4509,-.4,0.05,.1775
4901,-.6,.1,.1775
4909,-.4,.1,.1775
*NGEN,NSET=BACKBOT
101,109
*NGEN,NSET=BACKTOP
501,509
*NGEN,NSET=FRONTBOT
4501,4509
*NGEN,NSET=FRONTTOP
4901,4909
*NFILL,NSET=BACK
BACKBOT,BACKTOP,4,100
*NFILL,NSET=FRONT
FRONTBOT,FRONTTOP,4,100
*NFILL,NSET=ALL
BACK,FRONT,4,1100
*NSET,NSET=SYMMX, GEN
109,4509,1100
209,4609,1100
309,4709,1100
409,4809,1100
509,4909,1100
*NSET,NSET=SYMMZ
FRONT
*ELEMENT,TYPE=C3D8R,ELSET=SOLID
101, 101,102,202,201,1201,1202,1302,1301
*ELGEN,ELSET=SOLID
101, 8,1,1, 4,100,100,4,1100,1100
*ELSET,ELSET=TOP
401,402,403,404,405,406,407,408
**
**
*NODE
50101,-.59,-.15,-0.0125
50118,-.41,-.15,-0.0125
50501,-.59,0.05,-0.0125
50518,-.41,0.05,-0.0125
61101,-.59,-.15,.1675
61118,-.41,-.15,.1675
61501,-.59,0.05,.1675
61518,-.41,0.05,.1675
*NGEN,NSET=BACKBOT1
50101,50118
*NGEN,NSET=BACKTOP1
50501,50518
*NGEN,NSET=FRONTBOT1
61101,61118
*NGEN,NSET=FRONTTOP1
61501,61518
*NFILL,NSET=BACK1
BACKBOT1,BACKTOP1,4,100
*NFILL,NSET=FRONT1
FRONTBOT1,FRONTTOP1,4,100
*NFILL,NSET=ALL1
BACK1,FRONT1,10,1100
*NSET,NSET=SYMMX, GEN
50118,61118,1100
50218,61218,1100
50318,61318,1100
50418,61418,1100
50518,61518,1100
*NSET,NSET=SYMMZ
FRONT1
*NSET, NSET=GROUND,GEN
50101,50118,1
51201,51218,1
52301,52318,1
53401,53418,1
54501,54518,1
55601,55618,1
56701,56718,1
57801,57818,1
58901,58918,1
60001,60018,1
61101,61118,1
*ELEMENT,TYPE=C3D8R,ELSET=SOLID1
50101, 50101,50102,50202,50201,51201,51202,51302,51301
*ELGEN,ELSET=SOLID1
50101, 17,1,1, 4,100,100,10,1100,1100
**10,1000,1000
*ELSET,ELSET=TOP1,GEN
50401,50417,1
*SOLID SECTION,ELSET=SOLID,MATERIAL=MAT, CONTROLS=CONT
*SOLID SECTION,ELSET=SOLID1,MATERIAL=MAT, CONTROLS=CONT
*SECTION CONTROLS, NAME=CONT,HOURGLASS=ENHANCED, DISTORTION CONTROL=NO
*MATERIAL,NAME=MAT
*ELASTIC
 3000.E+6, 0.0
*CRUSHABLE FOAM,HARDENING=ISOTROPIC
 1.1, 0.2983
** plastic Poisson's ratio is chosen in such a way 
** that the plastic flow is associative.
*CRUSHABLE FOAM HARDENING
 220.E+6,  0.0
 246.5121E+6,  0.1
 272.9380E+6,  0.2
 299.0192E+6,  0.3
 324.5491E+6,  0.4
 349.3453E+6,  0.5
 373.2598E+6,  0.6
 396.1709E+6,  0.7
 418.0111E+6,  0.8
 438.7210E+6,  0.9
 458.2682E+6,  1.0
 493.8430E+6,  1.2
 524.8363E+6,  1.4
 551.5284E+6,  1.6
 574.3073E+6,  1.8
 593.5948E+6,  2.0
 629.3596E+6,  2.5
 651.9947E+6,  3.0
 683.3387E+6,  5.0
 688.3322E+6, 10.0
*DENSITY
 500.,
*BOUNDARY
GROUND,1,3
*SURFACE INTERACTION, NAME=TIE
*COHESIVE BEHAVIOR, Eligibility = Current Contacts 
**1.E+11,1.E+10,1.E+10
*DAMAGE INITIATION, CRITERION=MAXS
2.5E+8,2.5E+8,2.5E+8
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT,SOFTENING=LINEAR
.05,
***contact damping, definition=critical damping fraction
**.1,
*SURFACE, NAME=TOPSURF
SOLID,
*SURFACE, NAME=BOTSURF
SOLID1,
*ELSET, ELSET=BOTSURFTIE
50401, 50402, 50403, 50404, 50405, 50406, 50407, 
50408, 50409, 50410, 50411, 50412, 50413, 50414, 50415, 
50416, 50417, 51501, 51502, 51503, 51504, 51505, 
51506, 51507, 51508, 51509, 51510, 51511, 51512, 51513, 
51514, 51515, 51516, 51517, 52601, 52602, 52603, 
52604, 52605, 52606, 52607, 52608, 52609, 52610, 52611, 
52612, 52613, 52614, 52615, 52616, 52617, 53701, 
53702, 53703, 53704, 53705, 53706, 53707, 53708, 53709, 
53710, 53711, 53712, 53713, 53714, 53715, 53716, 
53717, 54801, 54802, 54803, 54804, 54805, 54806, 54807, 
54808, 54809, 54810, 54811, 54812, 54813, 54814, 
54815, 54816, 54817, 55901, 55902, 55903, 55904, 55905, 
55906, 55907, 55908, 55909, 55910, 55911, 55912, 
55913, 55914, 55915, 55916, 55917, 57001, 57002, 57003, 
57004, 57005, 57006, 57007, 57008, 57009, 57010, 
57011, 57012, 57013, 57014, 57015, 57016, 57017, 58101, 
58102, 58103, 58104, 58105, 58106, 58107, 58108, 
58109, 58110, 58111, 58112, 58113, 58114, 58115, 58116, 
58117, 59201, 59202, 59203, 59204, 59205, 59206, 
59207, 59208, 59209, 59210, 59211, 59212, 59213, 59214, 
59215, 59216, 59217, 60301, 60302, 60303, 60304, 
60305, 60306, 60307, 60308, 60309, 60310, 60311, 60312, 
60313, 60314, 60315, 60316, 60317
*ELSET, ELSET=TOPSURFTIE
101, 102, 103, 104, 105, 106, 107, 108, 1201, 1202, 
1203, 1204, 1205, 1206, 1207, 1208, 2301, 2302, 
2303, 2304, 2305, 2306, 2307, 2308, 3401, 3402, 
3403, 3404, 3405, 3406, 3407, 3408
*ELSET, ELSET=TOPSURFPULL
 401, 402, 403, 404, 405, 406, 407, 408, 1501, 1502, 
1503, 1504, 1505, 1506, 1507, 1508, 2601, 2602, 
2603, 2604, 2605, 2606, 2607, 2608, 3701, 3702, 3703, 
3704, 3705, 3706, 3707, 3708
*SURFACE, NAME=TOPSURFTIE
TOPSURFTIE,S3
*SURFACE, NAME=BOTSURFTIE
BOTSURFTIE,S5
*SURFACE, NAME=TOPSURFPULL
TOPSURFPULL,S5
*NODE, NSET=REFN
100000,-.5,.3,.0775
*COUPLING, CONSTRAINT NAME=KC,REFNODE=100000,SURFACE=TOPSURFPULL
*KINEMATIC
1,6
*AMPLITUDE, NAME=AMP,DEFINITION=SMOOTH STEP
0.0,0.0,.6,1.0
*timepoints, name=timepoints, gen
.3,.35,.05E-3
*BOUNDARY
100000,1,1,0.0
100000,3,6,0.0
*STEP
*DYNAMIC,EXPLICIT, element by element
**,.1
**,2.E-4
**,4.E-5
,.14
*BOUNDARY, AMPLITUDE=AMP
100000,2,2,1.0
*CONTACT
*CONTACT INCLUSIONS, ALL EXTERIOR
*CONTACT PROPERTY ASSIGNMENT
BOTSURFTIE,TOPSURFTIE,TIE
*contact controls assignment, type=scale penalty
BOTSURFTIE,TOPSURFTIE,1.
*OUTPUT, FIELD, NUMBER INTERVAL=4
***output, field, timepoints=timepoints
*NODE OUTPUT
U,
***ELEMENT OUTPUT
**S,
***CONTACT OUTPUT
**CFORCE,
*OUTPUT, HISTORY,FREQ=100
*NODE OUTPUT, NSET=REFN
U,
**RF
*ENERGY OUTPUT
ALLIE,ALLAE,ALLWK,ALLKE,ALLVD,ETOTAL,ALLPW
*END STEP