*heading
Test for viscoelasticity with 3D cohesive elements (COH3D8) with
traction-separation response.
*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
*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
*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
*node output
u,
*output, history, time interval=0.001
*element output, elset=qa_test_eall
s, le, ne
*energy output, variable=all
*end step
*step, nlgeom=yes
** Relaxation step
*dynamic, explicit
,1.
*boundary, type=velocity
upper1,3,3,0.
upper1,1,2,0.
upper2,3,3,0.
upper2,1,2,0.
*temperature
nall, 615.0
*end step