Stress output issue(s) in version 8.6?

[Marin]

Hello,

I believe there is a small issue with the stress output routine in the 8.6 version. Please refer to the attached image that shows differences in the stress output (for the exact same input) for the 8.5 and the 8.6 version. I’ve noticed this happens with two different elements: 3-node shell and 3-node 2D solid. Attached image should explain everything.

Best regards,
Marin

[Prof. R.L. Taylor]

Thank you for noting the error.  It is necessary to change the write for 2003 in shl3ds.f (located in ./elements/shells) to read

            write(iow,2003) n_el,xx,norm(1:6),ma,yy,momt(1:6),
     &                      zz,eps,sigt,sigb
            if(ior.lt.0) then
              write(*,2003) n_el,xx,norm(1:6),ma,yy,momt(1:6),
     &                      zz,eps,sigt,sigb

Also change the Surv.. to Curv.. in the format 2002 in 3 places

We will correct in future releases.

[Marin]

Ok, great.
Also, angle format could be changed for consistency e.g.

from
2003  format(/i5,0p,1f8.3,1p,5e11.3,0p,1f8.2/i5,0p,1f8.3,1p,5e11.3,
     &       0p,1f8.2/5x,0p,1f8.3,1p,6e11.3/(13x,1p,6e11.3))

to
2003  format(/i5,0p,1f8.3,1p,5e11.3,0p,1f8.2/i5,0p,1f8.3,1p,5e11.3,
     &       0p,1f8.2/5x,0p,1f8.3,1p,6e11.3/(13x,1p,5e11.3,0p,1f8.2))

Additionally, 3-node 2D solid routine should also be modified for correct stress output (as shown on the right hand side of the previously attached image)