Unusual Stress Output for ShellDKGQ

[epsilon]

Hi all,
I am trying to analyze a I beam model with ShellDKGQ element, Here is images of my model generated in another software:
- Model overall look: https://ibb.co/kHK7gk
- Model with showing nodes (red) and supports (green) and applied forces (pink): https://ibb.co/edVzSQ

- input file model.tcl: https://pastebin.com/VQ7RXsr6
- output file elements.xml: https://pastebin.com/NfB0f8FD
- output file nodes.xml: https://pastebin.com/kibTVejP

I want to get Nodal Displacements and Elements Internal Stresses at Gaussian point into xml files, so I use these TCL commands:

recorder Node -xml "nodes.xml" -dof 1 2 3 4 5 6 disp
recorder Element -xml "elements.xml" -time -closeOnWrite stresses

When i analyze the file like this in OpenSEES command line:

OpenSees > source model.tcl

The nodal displacements and element forces are written successfully.
The problem is that shell elements internal force are very odd, only 0 and 1 in elements.xml file! (please have a look at above link for elements.xml)
Am i making a mistake in definition or other places or elements have this very odd type of internal force?

Any comment would be appreciated.
Thanks