Reply To: Discrepancy between LSPrepost and hand calculations -Shell element stress output

[mmia1]

Thank you for the reply. 

What I have understood yesterday while working is that the discrepancy occurred due to maxint in *Database_Extent_Binary and NIP(number of integration point) in *Section_Shell. I considered default maxint =3 while in *Section_Shell I considered nip as 2. 

 

For shell fully integrated element, LSPREPOST showed results at reference surface which is middle surface and this is corresponding to first integration point. In case of maxint=3, there are middle surface(first integration point), inner surface (second integration point) and outer surface (third integration point).

 

So I revised the k file and used nip=3 and maxint =3 and when I read the first integration point results only it matched with lsprepost values. 

However, if elements are not distorted or for small displacement case, then the integration point results do not vary too much and the above thing does not matter.

 

I will ask a last thing now. When I checked the output file, the inner surface/lower surface (second integration point) always have larger stress compared to middle surface (first integration point). 

In this regard, what to do? Should I read only first integration point (middle surface) stress value to be consistent with lsprepost values or I should read the second integration point i.e. maximum stress value out of three integration points?

 

Thanks again.