*Heading
Test of Surface generated by Cutting plane method (ABAQUS/Explicit)
Test Integrated output
Elements Tested: B21
*Node, nset=nall
      1,          14.,          10.
      2,          14.,         -10.
      3,          34.,         -10.
      4,          34.,          10.
      5,          14.,   3.33333325
      6,          14.,  -3.33333325
      7,    20.666666,         -10.
      8,    27.333334,         -10.
      9,          34.,  -3.33333325
     10,          34.,   3.33333325
     11,    27.333334,          10.
     12,    20.666666,          10.
     13,          24.,  -4.22650003
     14,   26.5841846,    2.5841856
     15,   20.2642021,   2.63760972
     16,   19.3745747,  -3.40534234
     17,   28.6254253,  -3.40534234
     18,          36.,          10.
     19,          36.,   3.33333325
     20,          36.,  -3.33333325
     21,          36.,         -10.
     22,          38.,          10.
     23,          38.,   3.33333325
     24,          38.,  -3.33333325
     25,          38.,         -10.
     26,          38.,   6.66666651
     27,          38.,           0.
     28,          38.,  -6.66666651
     29,          42.,          10.
     30,          42.,   3.33333325
     31,          42.,  -3.33333325
     32,          42.,         -10.
     33,          44.,          10.
     34,          44.,         -10.
     35,          64.,         -10.
     36,          64.,          10.
     37,          44.,   3.33333325
     38,          44.,  -3.33333325
     39,   50.6666679,         -10.
     40,   57.3333321,         -10.
     41,          64.,  -3.33333325
     42,          64.,   3.33333325
     43,   57.3333321,          10.
     44,   50.6666679,          10.
     45,          54.,  -4.22650003
     46,   56.5841866,    2.5841856
     47,   50.2642021,   2.63760972
     48,   49.3745766,  -3.40534234
     49,   58.6254234,  -3.40534234
     50,   52.3333321,  -7.11324978
     51,          54.,         -10.
     52,   55.6666679,  -7.11324978
     53,    50.465435,   6.31880474
     54,   47.3333321,   6.66666651
     55,    47.132103,   2.98547149
     56,   50.0206223,  -6.70267105
     57,   51.6872864,  -3.81592131
     58,   46.6872864,  -6.70267105
     59,   47.3333321,         -10.
     60,   61.3127136, -0.0360045433
     61,   60.2920914,   2.95875931
     62,    57.604805,  -0.41057837
     63,   61.3127136,  -6.70267105
     64,          64.,  -6.66666651
     65,   61.3127136,   -3.3693378
     66,   55.2920914, -0.821157217
     67,   56.3127136,  -3.81592131
     68,   60.6666679,   6.66666651
     69,          64.,   6.66666651
     70,   60.6666679,          10.
     71,   53.4241943,   2.61089754
     72,   53.6254272,    6.2920928
     73,   47.3333321,          10.
     74,          44.,   6.66666651
     75,   49.8193893,  -0.38386631
     76,   46.6872864, -0.0360045433
     77,   57.9793777,  -6.70267105
     78,   60.6666679,         -10.
     79,   56.9587593,    6.2920928
     80,          54.,          10.
     81,          44.,           0.
     82,   46.6872864,   -3.3693378
     83,    52.132103, -0.794445157
     84,          44.,  -6.66666651
     85,          64.,           0.
     86,          40.,          10.
     87,          40.,   3.33333325
     88,          40.,  -3.33333325
     89,          40.,         -10.
     90,          40.,   6.66666651
     91,          40.,           0.
     92,          40.,  -6.66666651
*Element, type=B21,elset=b12
21, 18, 19
22, 19, 20
23, 20, 21
*elset,elset=allbeams
b12,
*Beam Section, elset=allbeams, material=STEEL, section=CIRC
1.
0.,0.,-1.
**
*SURFACE, TYPE=CUTTING SURFACE, NAME=PARALLEL_beam_POS
40,0,0,0,1,0
allbeams
*SURFACE, TYPE=CUTTING SURFACE, NAME=PARALLEL_beam_neg
40,0,0,0,-1,0
allbeams
*SURFACE, TYPE=CUTTING SURFACE, NAME=angle_beam_pos
40,0,0,.5,1,0
allbeams
*SURFACE, TYPE=CUTTING SURFACE, NAME=angle_beam_neg
40,0,0,-.5,-.8,0
allbeams
*preprint,model=yes,history=yes
** 
** MATERIALS
** 
*Material, name=STEEL
*Elastic
 2.1e+11, 0.3
*density
7800
*nset,nset=fix
21,25,32,18,22,2,7,8,3
34,59,39,51,40,78,35
*INITIAL CONDITION, TYPE=STRESS
allbeams, 30000.0
******************
*nset, nset=nref
9999,
*nset, nset=ioref
9998,
*nset, nset=ioref_move
9997,
*rigid body, tie nset=nall, ref=nref
*BOUNDARY
nref,1,2
*AMPLITUDE,NAME=SMOOTH, DEFINITION=SMOOTH STEP
 0.,0.,1.,0.025
*integrated output section, surf=PARALLEL_beam_POS,  name=PARALLEL_beam_POS_sec,
ref node=ioref
*integrated output section, surf=PARALLEL_beam_neg, name=PARALLEL_beam_neg_sec
*integrated output section, surf=angle_beam_POS, name=angle_beam_POS_sec
ref node=ioref_move,  position=center
*integrated output section, surf=angle_beam_neg, name=angle_beam_neg_sec
****************
*step
*DYNAMIC,EXPLICIT
,1.
*****************
*BOUNDARY,TYPE=DISPLACEMENT,AMPLITUDE=SMOOTH
nref,6,6,62.8318531
*output, history,freq=1
**
*node output, nset=nall
u,
*integrated output, surf=PARALLEL_beam_POS
SOF,SOM,SOAREA
*integrated output, surf=PARALLEL_beam_neg
SOF,SOM,SOAREA
*integrated output, surf=angle_beam_POS
SOF,SOM,SOAREA
*integrated output, surf=angle_beam_neg
SOF,SOM,SOAREA
**
*integrated output, sec=PARALLEL_beam_POS_sec
SOF,SOM,SOAREA
*integrated output, sec=PARALLEL_beam_neg_sec
SOF,SOM,SOAREA
*integrated output, sec=angle_beam_POS_sec
SOF,SOM,SOAREA
*integrated output, sec=angle_beam_neg_sec
SOF,SOM,SOAREA
*element output, elset=allbeams
sf,se
*****************
*end step