*heading
Test for viscoelasticity with 3D cohesive elements (COH3D8) with
traction-separation response, including damage and failure.
*preprint, model=yes
*node, nset=nall
 1, 0., 0., 0.
 2, 1., 0., 0.
 3, 1., 1., 0.
 4, 0., 1., 0.
 5, 0., 0., 0.1
 6, 1., 0., 0.1
 7, 1., 1., 0.1
 8, 0., 1., 0.1
**
11, 0., 0., 1.
12, 1., 0., 1.
13, 1., 1., 1.
14, 0., 1., 1.
15, 0., 0., 1.1
16, 1., 0., 1.1
17, 1., 1., 1.1
18, 0., 1., 1.1
**
101, 2., 0., 0.
102, 3., 0., 0.
103, 3., 1., 0.
104, 2., 1., 0.
105, 2., 0., 0.1
106, 3., 0., 0.1
107, 3., 1., 0.1
108, 2., 1., 0.1
**
111, 2., 0., 1.
112, 3., 0., 1.
113, 3., 1., 1.
114, 2., 1., 1.
115, 2., 0., 1.1
116, 3., 0., 1.1
117, 3., 1., 1.1
118, 2., 1., 1.1
*nset, nset=upper1
5,6,7,8
105,106,107,108
*nset, nset=upper2
15,16,17,18
115,116,117,118
*nset, nset=lower
1,2,3,4
11,12,13,14
101,102,103,104
111,112,113,114
*element, type=coh3d8, elset=COH1
1,1,2,3,4,5,6,7,8
*element, type=coh3d8, elset=COH2
10,11,12,13,14,15,16,17,18
*element, type=coh3d8, elset=COH3
101,101,102,103,104,105,106,107,108
*element, type=coh3d8, elset=COH4
110,111,112,113,114,115,116,117,118
*elset, elset=QA_TEST_EALL
 1,10,101,110
*boundary
lower, 1,3
*cohesive section, elset=COH1, material=mat1, response=traction separation, thickness=geometry
*cohesive section, elset=COH2, material=mat1, response=traction separation, thickness=geometry
*cohesive section, elset=COH3, material=mat2, response=traction separation, thickness=geometry
*cohesive section, elset=COH4, material=mat2, response=traction separation, thickness=geometry
**
*material, name=mat1
*density
1.
*elastic, type=traction, moduli=long term
2.e6, 1.e6, 1.e6
*damage initiation, criterion=maxe
 0.1, 0.1, 0.1
*damage evolution,type=displacement,mode mix ratio=energy, mixed mode behavior=tabular, softening=linear
 0.2, 0., 0.
 0.2, 1., 0.
 0.2, 0., 1.
 0.2, 1., 1.
*viscoelastic,time=prony
 0.6,0.6,0.1
*trs, definition=wlf
 5.25, 0.163277, 269.23
**
*material, name=mat2
*density
1.
*elastic, type=coupled, moduli=long term
2.e6, 1.e6, 1.e6, 5.e5, 5.e5, 2.e5
*damage initiation, criterion=maxe
 0.1, 0.1, 0.1
*damage evolution,type=displacement,mode mix ratio=energy, mixed mode behavior=tabular, softening=linear
 0.2, 0., 0.
 0.2, 1., 0.
 0.2, 0., 1.
 0.2, 1., 1.
*viscoelastic,time=prony
 0.6,0.6,0.1
*trs, definition=wlf
 5.25, 0.163277, 269.23
**
*initial conditions, type=temperature
nall, 615.0
**
*amplitude, name=disp
0., 0., 1., .4
*step, nlgeom=yes
*dynamic, explicit
,1.
*boundary, type=displacement, amplitude=disp
upper1,3,3,1.
upper1,1,2,0.
upper2,3,3,0.
upper2,1,2,1.
*temperature
nall, 615.0
*output, field, number interval=20
*element output
s, le, sdeg, dmicrt, status
*node output
u,
*output, history, time interval=0.001
*element output, elset=qa_test_eall
s, sdeg, le, ne, dmicrt, status
*energy output, variable=all
*end step