****master
***geometry
 **point point0
  *position 0.000000000000000e+00 -5.000000000000000e-01 0.000000000000000e+00
 **point point1
  *position 5.000000000000000e+00 -5.000000000000000e-01 0.000000000000000e+00
 **point point2
  *position 5.000000000000000e+00 5.000000000000000e-01 0.000000000000000e+00
 **point point3
  *position 0.000000000000000e+00 5.000000000000000e-01 0.000000000000000e+00
 **point point4
  *position 5.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
 **point point5
  *position 0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
 **line line0
  *p1 point0
  *p2 point1
  *cuts 40
 **line line1
  *p1 point2
  *p2 point3
  *cuts 40
 **line line2
  *p1 point2
  *p2 point4
  *cuts 5
  *progression 9.090909090909091e-01
 **line line3
  *p1 point4
  *p2 point1
  *cuts 5
  *progression 1.100000000000000e+00
 **line line4
  *p1 point3
  *p2 point5
  *cuts 5
  *progression 9.090909090909091e-01
 **line line5
  *p1 point5
  *p2 point0
  *cuts 5
  *progression 1.100000000000000e+00
 **line line6
  *p1 point5
  *p2 point4
  *cuts 40
 **domain domain1 Ruled1
  *name beam
  *element c2d8r
  *method 1
  *b1 line0 
  *b2 line5 
  *b3 line6 
  *b4 line3 
 **domain domain1 Ruled2
  *name beam
  *element c2d8r
  *method 1
  *b1 line1 
  *b2 line4 
  *b3 line6 
  *b4 line2 
 **domain bset right
  *name right
  *order 2
  *b line3 line2 
 **domain bset left
  *name left
  *order 2
  *b line4 line5 
***mesher
  **nset 
    *type cartesian
    *name load_pt 
    *axes 0 1 2 
    *limit 0.00100000 
    *func (abs(x-4.)<0.01)*(y>0.4999);
  **nset 
    *type cartesian
    *name sup_bottom 
    *axes 0 1 2 
    *limit 0.00100000 
    *func (abs(x-4.)<0.001)*(y<-0.4999);
  **extension 
    *elset ALL_ELEMENT 
    *new_elset ALL_ELEMENT 
    *prog 1.00000 
    *double_prog 0.00000 
    *fusion 0.00100000 
    *distance 2.00000 
    *num 10 
    *dir ( 0.00000 0.00000 1.00000 )
    *liset_names right left 
    *elset_names beam 
  **nset 
    *type cartesian
    *name sup_mid 
    *axes 0 1 2 
    *limit 0.00100000 
    *func (abs(x-5.)<0.001)*(abs(y)<0.001)*(abs(z-1.)<0.001);
  **nset 
    *type cartesian
    *name left_mid 
    *axes 0 1 2 
    *limit 0.00100000 
    *func (x<0.001)*(abs(y)<0.001)*(abs(z-1.)<0.001);
  **nset 
    *type cartesian
    *name left-midy 
    *axes 0 1 2 
    *limit 0.00100000 
    *func (x<0.001)*(abs(y)<0.001);
  **nset 
    *type cartesian
    *name left-midz 
    *axes 0 1 2 
    *limit 0.00100000 
    *func (x<0.001)*(abs(z-1.)<0.001);
***plotting
  **default_options 1 0 1
****return