FEAP User Forum
FEAP => Programming => Topic started by: ankitk on January 22, 2020, 04:42:56 PM
-
All,
I found a user macro (originally created for FEAP 8.2) that provides a TECPLOT output (http://tuprints.ulb.tu-darmstadt.de/1336/1/main.pdf). I have modified the it mainly for syntax and added a few more formats.
When I compile FEAP 8.5, there are no issues however when I call the user macro, the files created are blank. Any assistance would be greatly provided.
Thank you,
AK
-
It looks like one issue is an extra "close(99)" just before the formats. With that removed the data goes to the correct file.
-
there are quite a few typographical errors in the file. Also the pointer to the stress values is off by "1" so that v. Mises appears in J3 instead of J2.
Be cautious also in later feap versions the force and displacement values are separated in the np(27) array. Maybe using "paraview" is better?
-
Sir,
Thank you for the reply. I recompiled the program without the close(99) and am slowly working through the typographical errors. However the .dat files that are created are still empty. The batch command I use is:
BATCh
DT,,1
LOOP,,5
TECP, init ! open output unit
TIME
LOOP,,1
TANG,,1
NEXT
STRE,NODE ! compute nodal stresses (mandatory!)
TECP,write ! write postprocessing data
TECP,close ! close output unit
NEXT
I have attached the Input file as well. I am using paraview at this time however my professor would prefer if we can get tecplot outputs.
Once again thank you for your assistance.
-
You need TECP INIT before each write
-
What version of TECPLPOT are you using. New versions are supposed to be able to read VTU files (which is what we output for paraview). You should be able to use our paraview routines to create VTU files and just read them in TECPLOT