belondo,
I just took a quick look in the NASTRAN Quick Reference Guide for PSHELL card entries. The 7th entry on the PSHELL card (MID3) determines whether or not transverse shear is included.
If you are using a preprocessor, look at the bulk data file it generates to see if the PSHELL card contains the material ID for tranverse shear.
For example, there should be a line in the bulk data file that looks like:
PSHELL PID MID1 T MID2 I21/T3 MID3 TS/T
Where,
PID = property ID
MID1 = Material ID 1 for membrane stiffness
T = thickness
MID2 = Material ID 2 for bending
12I/T3 = bending moment inertia ratio (default is 1.0 when this entry is blank)
MID3 = Material ID for transverse shear
TS/T = transverse shear thickness ratio (default is 0.8333 if this entry is blank).
Actual line might look like
PSHELL 5 1 .050 1 1
or
PSHELL,5,1,.050,1,,1,
Here the property ID is 5, the material ID for membrane, bending and transverse shear is 1, the thickness is 0.050 and defaults are used for 12I/T3 and TS/T. Note that the material properties are specified by the MAT1 card (for linear isotropic materials).
Bottom line, I think you're going to want to check the reference guide for the PSHELL and MAT entries required to get transverse shear, and also check the theoretical manual to make sure you are getting an element that behaves the way you want it to behave. Lastly, make sure the case control section of the BDF is requesting all available stress output.
Good luck