PARAVIEW Error

[Yaakov]

Dear Prof. Taylor & FEAP TEAM,

Thanks very much for your warm help!

I will give my feed back as soon as possible!

Kind Regards, Yaakov

[Yaakov]

Dear Prof. Taylor, Dear Prof. Govindjee & Dear FEAP TEAM,

I put the three new files (Reply #14) into feap user order ( make clean & install).
Now It works as well as I expect that PLOT PVIE can stay after PLOT CONT without bugs(e.g. for 8000 elements, or more elements) !

Kind Regards & Nice weekend, Yaakov


 

   

[resammc]

Can you try the attached files for paraview.  They are supposed to output all the stress values even if you have not projected them yet.  A projection is forced if an iteration has been performed and no stresses plotted or nodal stresses computed.

I have a problem with my non-mechanical user element using PVIE with the new FEAP8.5. Whenever I try to make a paraview output, the program exits.

My problem is a phase-field problem and essentially I just need to output the DOF fields and there is no need to project anything. I thought maybe it is because I did not define isw.eq.8 in my element, and added that and also set the ISTE value, then I got another error that the assigned size to NDER and NDNP are zero...

I believe the problem comes from the if condition for FL(11) in program/pmacr6.f (feap 8.5). It forces the program to project stresses, while there isn't any, thus the error.

I wrote a small user-macro where I temporarily switch the value of FL(11) to TRUE before calling PVIE and switch it back after the call. This way, the routine works just fine.

Is there any better work-around for this?

RM

[Prof. R.L. Taylor]

One other option is to change the coding in pmacr6.f  to


!     [pvie,<time,file> : Paraview interface

      elseif(j.eq.12) then

!       No stress projection case

        istvsv = istv
        if(ct(1).ne.0.0d0) then
          istv = 0

!       Compute and project stresses

        elseif(.not.fl(11)) then
          call pjstrs(.false.)
        endif

        call pparaview(lct)
        istv = istvsv

be sure to add dimension for 'istvsv'

Then if you use: PVIEw ,, 1

you get no stress output.

[resammc]

thank you very much, everything works fine now.