Question about calling subroutine to get strain and stress
I am writing a subroutine DFLUX, the FLUX(1) is dependent on strain and stress. But I don't know how to get the strain and stress. I use SPRINC in the code but it seems to be wrong.
The equation for the problem is
q=(α^2XΩXε^2XE)/2
E=(σ/ε)*sin(d)
σ is stress and ε is strain