*Heading
CPE4 elements
elastic materials
static procedure
**
*Node
      1,          30.,           0.
      2,          30.,          20.
      3,          10.,          20.
      4,          10.,           0.
      5,          50.,           0.
      6,          50.,          20.
      7,          30.,           4.
      8,          30.,           8.
      9,          30.,          12.
     10,          30.,          16.
     11,      21.2158,          20.
     12,      15.6821,          20.
     13,      12.1961,          20.
     14,          10.,          16.
     15,          10.,          12.
     16,          10.,           8.
     17,          10.,           4.
     18,      12.1961,           0.
     19,      15.6821,           0.
     20,      21.2158,           0.
     21,      38.7842,           0.
     22,      44.3179,           0.
     23,      47.8039,           0.
     24,          50.,           4.
     25,          50.,           8.
     26,          50.,          12.
     27,          50.,          16.
     28,      47.8039,          20.
     29,      44.3179,          20.
     30,      38.7842,          20.
     31,      21.2158,           4.
     32,      21.2158,           8.
     33,      21.2158,          12.
     34,      21.2158,          16.
     35,      15.6821,           4.
     36,      15.6821,           8.
     37,      15.6821,          12.
     38,      15.6821,          16.
     39,      12.1961,           4.
     40,      12.1961,           8.
     41,      12.1961,          12.
     42,      12.1961,          16.
     43,      38.7842,          16.
     44,      38.7842,          12.
     45,      38.7842,           8.
     46,      38.7842,           4.
     47,      44.3179,          16.
     48,      44.3179,          12.
     49,      44.3179,           8.
     50,      44.3179,           4.
     51,      47.8039,          16.
     52,      47.8039,          12.
     53,      47.8039,           8.
     54,      47.8039,           4.
*Nset, nset=top
  2,  3,  6, 11, 12, 13, 28, 29, 30
*Nset, nset=bottom
  1,  4,  5, 18, 19, 20, 21, 22, 23
*Nset, nset=left
  3,  4, 14, 15, 16, 17
*Nset, nset=right
  5,  6, 24, 25, 26, 27
*Nset, nset=centerpoint
 1,
*Nset, nset=centerline
  1,  2,  7,  8,  9, 10
*element,type=cpe4
 1,  1,  7, 31, 20
 2,  7,  8, 32, 31
 3,  8,  9, 33, 32
 4,  9, 10, 34, 33
 5, 10,  2, 11, 34
 6, 20, 31, 35, 19
 7, 31, 32, 36, 35
 8, 32, 33, 37, 36
 9, 33, 34, 38, 37
10, 34, 11, 12, 38
11, 19, 35, 39, 18
12, 35, 36, 40, 39
13, 36, 37, 41, 40
14, 37, 38, 42, 41
15, 38, 12, 13, 42
16, 18, 39, 17,  4
17, 39, 40, 16, 17
18, 40, 41, 15, 16
19, 41, 42, 14, 15
20, 42, 13,  3, 14
*Elset, elset=leftsection, generate
  1,  20,   1
*element,type=cpe4
21,  2, 10, 43, 30
22, 10,  9, 44, 43
23,  9,  8, 45, 44
24,  8,  7, 46, 45
25,  7,  1, 21, 46
26, 30, 43, 47, 29
27, 43, 44, 48, 47
28, 44, 45, 49, 48
29, 45, 46, 50, 49
30, 46, 21, 22, 50
31, 29, 47, 51, 28
32, 47, 48, 52, 51
33, 48, 49, 53, 52
34, 49, 50, 54, 53
35, 50, 22, 23, 54
36, 28, 51, 27,  6
37, 51, 52, 26, 27
38, 52, 53, 25, 26
39, 53, 54, 24, 25
40, 54, 23,  5, 24
*Elset, elset=rightsection, generate
 21,  40,   1
** 
*solid section, elset=rightsection,material=softelpl
1.0, 
*solid section, elset=leftsection,material=hardelpl
1.0, 
*Material, name=softelpl
*Elastic
1.E4, 0.3
*Plastic
8.E3, 
*conductivity
1.0E6, 
*density
1.0, 
*specific heat
1.0, 
*Material, name=hardelpl
*Elastic
1.E5, 0.3
*Plastic
8.E4, 
*conductivity
1.0E6, 
*density
1.0, 
*specific heat
1.0, 
** 
*map solution,step=1,inc=10
*Step, nlgeom, inc=100
Hold
*Static
1.,1.
** 
*boundary,op=new
bottom,ysymm
centerline,xsymm
top,2,2,0.0
right,1,1,0.0
left,1,1,0.0
*Restart, write, frequency=1
*Output, field, op=NEW, frequency=1
*Node Output
U, 
*Element Output
S, E
*node file,freq=9999
U, 
*el file,freq=9999
S, E
*End Step