I am using the following failure function inside my own 2D FE code for analysis of laminated composites:
F=(1/Xt-1/Xc)*S11+(1/Yt-1/Yc)*S22 +S11^2/(Xt*Xc)+S22^2/(Yt*Yc)+S12^2/(XY^2) -S11*S22/sqrt(Xt*Xc*Yt*Yc)
I sometimes get negative F values. Numerically, it seems to be possible to have negative F values, however, since the theory is originated from energy based theories, negative F values do not seem to make any sense.
Any thoughts on this?