There is an error in subroutine ./elements/plate2t.f around line 347 the statements read:
tt(1) = hh(1,1)*gamm(1) + hh(1,2)*gamm(2)
tt(2) = hh(2,1)*gamm(1) + hh(2,2)*gamm(2)
they should read
tt(4) = hh(1,1)*gamm(1) + hh(1,2)*gamm(2)
tt(5) = hh(2,1)*gamm(1) + hh(2,2)*gamm(2)
This occurs in ver8.4, ver8.5, ver8.6. It does not occur for the quadrilateral element.
The interpolation for element shear values may not be very accurate