feap  ** Pre-Stressing Using the LINK **
0 0 0 2 2 4

block          ! Block 1
 cart 5 5
  1  0 0
  2 .9 0
  3 .9 1
  4  0 1

block         ! Block 2
 cart  5 5
  1 1  0
  2 2  0
  3 2  1
  4 1  1

eboun         ! Fix left and right edges
1 0 1 0
1 2 1 0

cboun          ! Pin corner nodes
 node 0 0 1 1 
 node 2 0 1 1 

eboun          ! Drive right side of block 1
 1 0.9 1 1


edisp          ! Stretch block 1 to meet block 2 node for node
 1 0.9 0.1 0.0



mate
 solid
  finite
  elastic iso 10 .1


end


batch             ! Stretch block 1 and compute equilibrium state
  plot,defo
  plot,frame, 1   ! upper left frame
  plot,mesh
  plot,outl
  loop,,10
   tang,,1
  next
  plot,frame,2   ! upper right frame
  plot,stre,1
end


batch              ! Remove driven boundary codes on right side of block 1
  mesh
end
eboun,set          ! Use 'set' to actually set the BC codes 0 before linking default is add
 1 .9 0 0 0

end  ! end the mesh to reset bcs



tie    ! force the flag to reevaluate the equation numbering etc (should be fixed in
       ! future versions, and then this can be removed)

elink  ! link the nodes on the interface (glue them together)
 1 0.9 1 



batch       ! Solve for the new equilibrium
  loop,,10
   tang,,1
  next
  plot,frame,3      ! lower left frame
  plot,stre,1
end

inte
stop