User element stress projection - which array is used?

[FEAP_Admin]

What version of the code are you using?

[Jfrisch]

FEAP 8.5 and I tried "numnp" and "numnm" as the amount of nodes of mesh.

[Prof. R.L. Taylor]

Did you look carefully at the problem and example I posted a couple of days ago.  This small problem shows exactly what the arrays look like for the stress projection.  If you study it carefully you should be able to figure out the steps  you need to follow.

[FEAP_Admin]

If you have version 8.5, why are you not following the example given in the programmers manual, where you just set p( ) and s( ) that are passed to your element from FEAP?  FEAP then takes care of all the rest, namely dealing with the assembly into the proper locations and the solution of the equations.

Is there something that the standard set up will not accommodate that you are interested in?

[Jfrisch]

I wrote a subroutine following the manual (please see the initial post of this thread e.g.), but I need to get the numerical values of the projected stresses for external post-processing. So the main question of this thread is how to achieve this because during isw=8 the element-projection-array p(.) and s(.) are set (I only renamed them to residual and stiff).

To get the stresses at integration points e.g. I need to set the set of equations during isw=3 (and/or isw=6) and afterwards are able to access the solution by calling isw=4 and using the ul-array to calculate the stresses using the shape functions and then write them into an output-file e.g.

The first step of this analogy is done in isw=8 but I have no idea how to achieve the latter step of accessing the "solution" (the projected values) and writing them to an output-file by the element routine I coded.

[Prof. R.L. Taylor]

I looked back at your original posting  and the question you ask.  We already have a routine which writes to a file -- it is for the paraview use. 

Look at the routines for paraview and see how the stress is written to the file.

[Jfrisch]

Ah ok, I will do that! Many thanks for your explanations, answers and patience!

[FEAP_Admin]

see program/p_paraview.f and program/pparaview.f