Author Topic: using prtrea.f in umacr  (Read 3811 times)

PaulM

  • New Member
  • *
  • Posts: 4
using prtrea.f in umacr
« on: June 30, 2020, 09:25:46 AM »
Dear FEAP Community,

i want to save the sums of the reaction forces and moments similiar to the command REAC COOR d x tol.
Since I need a different output format than the one in the output file, I wanted to implement this with a user macro.

In this macro I call the routine "prtrea.f" to calculate the corresponding forces and moments for previously defined nodes.
First of all I want to make sure that the output of the umacr.f is correct and therefore I let the results be stored in the output file and compare them with the standard commands REAC ALL and REAC COOR.

Unfortunately I had to realize that the output of my macro depends on the previously used standard commands.
This is expressed by the fact that only by using REAC ALL or REAC COOR all relevant nodes are considered in the user macro, regardless of which I pass to the routine "prtrea.f".
If I do not use the standard commands, the last nodes of the mesh are not considered for the force calculation, which leads to wrong results.

This error also occurs when I save the output of my macro in a different file as recommended by Prof. Taylor in this post:
http://feap.berkeley.edu/forum/index.php?topic=1427.msg5169#msg5169

I have attached the relevant files and would be very grateful for any advice to solve the problem.

Best regards,
Paul

Prof. R.L. Taylor

  • Administrator
  • FEAP Guru
  • *****
  • Posts: 2649
Re: using prtrea.f in umacr
« Reply #1 on: June 30, 2020, 01:29:58 PM »
If you look in pmacr1.f you will observe that we call formfe  "6" with flags set to compute all the reactions not just those for
the active equations.  I suggest you study that part of the code in pmacr1.f to see what you will need  but at least the pload and formfe parts are required.