Please confirm that you have NIP=3 in *SECTION_SHELL and MAXINT=3 in *DATABASE_EXTENT_ BINARY.  If both are true, then I suspect the problem is due to averaging of the in-plane integration points.  Either set ELFORM to 2 (which has only one in-plane integration point) or else try setting MAXINT= -3 so that all 4 in-plane integration points of ELFORM 16 are output.  Otherwise, you're getting an average of the 4 in-plane integration points at each through-thickness location.
If NIP=3 and ELFORM=16, there are actually 3 x 4 = 12 integration points.
You are navigating away from the AIS Discovery experience