It looks like you are passing "nes" as "1" when you call gaussValueProjection. If you only want to project the one value you need to ensure that all the numbers match correctly since the stress is dimensioned sig(1,*) with what you have. What is the real dimension for the stress and other arrays??
Dear Prof. Taylor,
Thanks for your reply. Yes indeed, nes is 1, as I only want to plot the cohesive damage. Therefore, sig should be (1,*). As I want to plot it in the 16th slot, I do s(1:nel,16). However, since the other 15 are left blank, I have modified the code to plot the sig(1,*) into s(1:nel,1), the first slot. Still, I cannot manage to plot the cohesive damage neither in the built-in FEAP visualization tool nor in Paraview.
A couple of comments.
(1) Under isw.eq.1 when you set istv you should do it as istve=max(16,istv)
in case there are other elements in your mesh that need more slots.
(2) In your gaussValueProjection routine, it looks like you are putting stress-1 into plotting stress-16. You probably want something more likexg = xsj(l)*shp3(i,l)
r(i) = r(i) + xg
s(i,1:16) = s(i,1:16) + sig(1:16,l)*xg
Dear FEAP_Admin,
Thanks for your reply. I have tried doing both things (as I mentioned Prof. Taylor, I have modified the code to have only one slot. As nes = 1, I try to plot it in s(1:nel,1)), but still I cannot manage to have the results plotted.
In the new code, I am writing p(1:nel) and s(1:nel,*) into the terminal and as you can see, the results are being properly added into the matrices. Therefore, my guess is that the problem has to be in the input file.
JReinoso.