*HEADING
SIMPLE CONTACT TEST WITH C3D10MT ELEMENTS,
DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
*RESTART,WRITE,NUM=1
*NODE
1, 1., 5.
2, 1.5, 5.
3, 2., 5.
4, 1., 7.5
5, 2., 7.5
6, 1., 10.
7, 1.5, 10.
8, 2., 10.
9, 1., 5., 0.5
10, 2., 5., 0.5
11, 1., 10., 0.5
12, 2., 10., 0.5
13, 1., 5., 1.
14, 1.5, 5., 1.
15, 2., 5., 1.
16, 1., 7.5, 1.
17, 2., 7.5, 1.
18, 1., 10., 1.
19, 1.5, 10., 1.
20, 2., 10., 1.
21, 1.5, 5., 0.5
22, 1.5, 10., 0.5
23, 1.5, 7.5,
24, 1.5, 7.5, 1.
25, 2., 7.5, 0.5
26, 1., 7.5, 0.5
121, 1., 12.5
122, 1., 11.25
123, 1., 10.,
124, 1.25, 12.5,
125, 1.25, 10.,
126, 1.5, 12.5
127, 1.5, 11.25
128, 1.5, 10.
129, 1., 12.5, 0.5
130, 1., 10., 0.5
131, 1.5, 12.5, 0.5
132, 1.5, 10., 0.5
133, 1., 12.5, 1.
134, 1., 11.25, 1.
135, 1., 10., 1.
136, 1.25, 12.5, 1.
137, 1.25, 10., 1.
138, 1.5, 12.5, 1.
139, 1.5, 11.25, 1.
140, 1.5, 10., 1.,
144, 1.75, 12.5,
145, 1.75, 10.
146, 2., 12.5
147, 2., 11.25
148, 2., 10.
151, 2., 12.5, 0.5
152, 2., 10., 0.5
156, 1.75, 12.5, 1.
157, 1.75, 10., 1.
158, 2., 12.5, 1.
159, 2., 11.25, 1.
160, 2., 10., 1.
161, 1.25, 10., 0.5
162, 1.25, 12.5, 0.5
163, 1.25, 11.25,
164, 1.25, 11.25, 1.
165, 1.5, 11.25, 0.5
166, 1., 11.25, 0.5
167, 1.75, 10., 0.5
168, 1.75, 12.5, 0.5
169, 1.75, 11.25,
170, 1.75, 11.25, 1.
171, 2., 11.25, 0.5
*NSET,NSET=BOTTOM
1,2,3,9,10,13,14,15,21
*NSET,NSET=TOP
121,124,126,144,146,129,
162,131,168,151,133,136,
138,156,158
*NSET,NSET=BACK
1,2,3,
121,124,126,144,146
*ELEMENT, TYPE=C3D10MT, ELSET=SOLID
1, 15, 3, 1, 8, 10, 2, 21,
25, 5, 23
2, 1, 13, 15, 18, 9, 14, 21,
26, 16, 24
3, 8, 20, 18, 15, 12, 19, 22,
25, 17, 24
4, 18, 6, 8, 1, 11, 7, 22,
26, 4, 23
5, 15, 8, 1, 18, 25, 23, 21,
24, 22, 26
6, 140, 128, 123, 126, 132, 125, 161,
165, 127, 163
7, 123, 135, 140, 133, 130, 137, 161,
166, 134, 164
8, 126, 138, 133, 140, 131, 136, 162,
165, 139, 164
9, 133, 121, 126, 123, 129, 124 ,162,
166, 122, 163
10, 140, 126, 123, 133, 165, 163, 161,
164, 162, 166
11, 160, 148, 128, 146, 152, 145, 167,
171, 147, 169
12, 128, 140, 160, 138, 132, 157, 167,
165, 139, 170
13, 146, 158, 138, 160, 151, 156, 168,
171, 159 ,170
14, 138, 126, 146, 128, 131, 144, 168,
165, 127, 169
15, 160, 146, 128, 138, 171, 169, 167,
170, 168, 165
*SOLID SECTION, ELSET=SOLID, MATERIAL=SOLMAT
*MATERIAL, NAME=SOLMAT
*ELASTIC
30.E6,0.3
*EXPANSION
10.E-6,
*CONDUCTIVITY
43.0,
*DENSITY
7700.,
*SPECIFIC HEAT
600.,
**
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16,STEFANB=5.6697E-8
**
*** nset definitions for prescribing temperatures.
**
*NSET,NSET=NELEM15,GEN
1,26
*NSET, NSET=QA_TEST
NELEM15,
*NSET, NSET=NELEM620
122,123,125,127,128,130
132,134,135,137,139,140
141,142,143,145,147,148,149,150
152,153,154,155,157,159,160
161,163,164,165,166,167,169,170
** surfaces
**
*surface, NAME=ASURF
3, S1
4, S1
*surface, NAME=BSURF
6, S1
7, S1
11,S1
12,S1
**
*** Apply different temperatures to the nodes above and
*** below the contact surface
**
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT,EXPLICIT
, 500.
*fixed mass scaling, elset=solid,factor=100.0
*CONTACT PAIR, SMALL SLIDING, INTERACTION=GAP, WEIGHT=0.0
BSURF, ASURF
*SURFACE INTERACTION, NAME=GAP
1.0,
*FRICTION
0.01,
*GAP CONDUCTANCE
1.,0.
1.,1.
*GAP RADIATION
3.4664e-2, 3.4664e-2
1.,0.
1.,1.
*GAP HEAT GENERATION
0.5,0.5
*BOUNDARY,type=velocity
BACK,3,3
BOTTOM,2,2
TOP,2,2
1, 1, 1
9, 1, 1
13,1, 1
*AMPLITUDE,NAME=RAMP, DEFINITION=SMOOTH STEP
0.,0.,100.,1.,500.,1.
*BOUNDARY,AMPLITUDE=RAMP
NELEM15, 11, 11, 700.0
NELEM620, 11, 11, 500.0
************************
*FILE OUTPUT, NUM=1
*NODE FILE,NSET=NELEM620
NT,U,
************************
*OUTPUT,FIELD,NUM=1
*NODE OUTPUT,NSET=QA_TEST
NT,U,
*************
*NODE OUTPUT
U,NT
*ELEMENT OUTPUT
S
**********************************
*OUTPUT,HISTORY,TIME INTERVAL=500.0
*ENERGY OUTPUT, VAR=ALL
**********************************
*END STEP
**
**** STEP 2. Push top block to the right.
**
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT,EXPLICIT
,500.0
*fixed mass scaling, elset=solid,factor=100.0
**
*amplitude,name=ramp2,definition=smooth step
0.,0., 100.,1.0, 400.,1.0, 500.,0.
*BOUNDARY, OP=NEW, amplitude=ramp2,type=velocity
123, 1, 1, 0.000025
130, 1, 1, 0.000025
135, 1, 1, 0.000025
122, 1, 1, 0.000025
134, 1, 1, 0.000025
166, 1, 1, 0.000025
*BOUNDARY, OP=NEW
NELEM15, 11, 11, 700
NELEM620, 11, 11, 500
*BOUNDARY, OP=NEW,type=velocity
BACK,3,3
TOP, 2, 2
BOTTOM, 2, 2
1, 1, 1
9, 1, 1
13,1, 1
4,1,1
16,1,1
6,1,1
11,1,1
18,1,1
26,1,1
*END STEP
**
*** STEP 3. Push top block back into position.
**
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT,EXPLICIT
,500.0
*fixed mass scaling, elset=solid,factor=100.0
*BOUNDARY, OP=NEW, amplitude=ramp2,type=velocity
123, 1, 1, -0.000025
130, 1, 1, -0.000025
135, 1, 1, -0.000025
122, 1, 1, -0.000025
134, 1, 1, -0.000025
166, 1, 1, -0.000025
*BOUNDARY, OP=NEW
NELEM15, 11, 11, 700
NELEM620, 11, 11, 500
*BOUNDARY, OP=NEW,type=velocity
BACK,3,3
TOP, 2, 2
BOTTOM, 2, 2
1, 1, 1
9, 1, 1
13,1, 1
4,1,1
16,1,1
6,1,1
11,1,1
18,1,1
26,1,1
*END STEP