FEAP User Forum
FEAP => FEAPpv => Topic started by: tuyet on June 25, 2013, 02:04:26 AM
-
Dear professor Taylor,
Dear Feap Admin,
I tested a dynamics example in Feappv. This example ran well in Feappv version 2005. However it cannot run in the new version of Feappv posted in your website.
There are two problems:
1. After the first inte command, I tried to plot mesh. However Feappv quit the program immediately.
2. To solve the problem program the following command is used
batch ! MACRO n.3: solve problem
dt,,.4
tran,newm ! Transient analysis
loop,,200
loop,,4
time
loop,,10 ! Convergence loop
tang,,1
next
next
plot,wipe
plot,defo,1,1 ! Plot on deformed conf.
plot,fill ! Plot filled mesh
plot,unde,,1 ! plot on undefor. conf.
plot,outl ! Plot outlined mesh
next
end
But it has an error in output file:
**ERROR** Insufficient storage to allocate TANGS
Required size = 17910
Check data or choose other option.
This example also runs well in Feap version 8.2.
Could you please show me what is the problem of current version Feappv?
Thank you very much. I am looking forward to hearing from you.
Best regards,
Tuyet
-
Please post your input file so we can check what is happening.
-
1. I have tested your batch program with my input file and it works fine.
2. Have you solved other problems (static ones ) with no problems? If so, do a static solution
on your mesh to see if the allocation is O.K.
3. Did you use solution command "CHECk" to make sure mesh is o.k.?
-
Thank you very much for your response.
I attach here the input files: Imens, Imens.mac, Imens.msh. This example was introduced by L. Petrini & A. Reali as benchmark example for feappv.
Here are some problems
1. Check mesh command: it is fine.
2 Plot mesh after TIE command: it is fine
Plot mesh after
batch ! MACRO n.1: check mesh
noprint ! Reduce output printing
plot,defa,off
plot,pers
plot,hide
plot,fill
end
0
2000,300,2000
0,0,1
end
+ works in Feap v8.2
+ does not work under running my project in microsoft studio 8 + Intel
+ works under the built Feappv.exe
3. Static case:
+ works in Feap v8.2
+ does not work under running ''Start debugging'' icon under my project in microsoft studio 8 + Intel
+ works under the built Feappv.exe
4. Dynamics case:
+ works in Feap v8.2
+ does not work under running ''Start debugging'' icon under my project in microsoft studio 8 + Intel as well as the built Feappv.exe
Thank you very much.
My best regards,
Tuyet
-
There are two known errors in memory allocation that probably affect your solution. One is for proportional loads (which is probably the one affecting you now) and the other for compressed profile which will affect eigen-computations, etc.
A new version for the program should be posted soon and hopefully will fix your problem.
-
Thank you very much for your quick reply.
I am looking forward to the new version.
-
Tuyet, a new version has been uploaded to the FEAPpv web page.
-
Thank you very much for your updated code.
You program runs all right with the test problem posted in Feappv website. Unfortunately the program cannot run well for dynamics in my computer. I ran program in
1. Windows Xp, Microsoft Visual Studio 8 + Intel Fortrans in the first computer, which has the sam problem as I solve you before:
Computing solution at time 4.0000E-01: Total proportional load 4.0000E-01
Individual factors:
1 = 4.0000E-01 2 = 0.0000E+00 3 = 0.0000E+00
**ERROR** Insufficient storage to allocate TANGS
Required size = 17910
Check data or choose other option.
2. Windows 7, Microsoft Visual Studio 8 + Intel Fortrans Second computer, which happens another problem:
Computing solution at time 4.0000E-01: Total proportional load 4.0000E-01
Individual factors:
1 = 4.0000E-01 2 = 0.0000E+00 3 = 0.0000E+00
*ERROR* DCHECK: Attempt to input 1 value(s). Input is:
1
Inconsistency occurred from INPUTS in data file Imens.mac
at or near record number 65. Input record is:
end
If this record is correct error may result from
missing blank record before new command type.
I don't know what happens right now.
-
Which input file are you using now? Send again please.
-
Here are the attached files. Thank you very much.
-
I am not able to test your program with a windows machine at this time.
The only thing I can see in your files relates to use of END records. You have extra ones in the Imens.mac file (see attached).
I do not see how you can be getting the allocation problem unless you have a bad compile, that is you still have the old program being loaded. Otherwise there is no way that I can understand receiving this message still. I suggest you clean everything out, download the program again and recompile from the very beginning.
-
Thank you very much.
WIndows: I tried to recompile the new feappv program in windows and use new Imens.mac. It happens the same problem.
Linux: With the new feappv program, the program immediately interrupted after input the mesh. In the feappv screen, it appears
feappv: malloc.c:2451: sYSMALLOc: Assertion `(old_top == (((mbinptr) (((char *) &((av)->bins[((1) - 1) * 2])) - __builtin_offsetof (struct malloc_chunk, fd)))) && old_size == 0) || ((unsigned long) (old_size) >= (unsigned long)((((__builtin_offsetof (struct malloc_chunk, fd_nextsize))+((2 * (sizeof(size_t))) - 1)) & ~((2 * (sizeof(size_t))) - 1))) && ((old_top)->size & 0x1) && ((unsigned long)old_end & pagemask) == 0)' failed.
Aborted (core dumped)
-
tuyet,
I tried your problem on our Windows machine at it fails (though with a different error).
It seems that your problem is causing FEAPpv to generate a memory. We will look into
it further and get back to you.
-
Thank you very much, Feap Admin.
I am looking forward to hearing from you.
-
A new version has been posted. Your problem ran on our windows machine with the changes -- please give it a try now.
-
Thanks for new version.
Your program works now. However there are still some following errors:
1. In input file we use:
plot,wipe
plot,defo,1,1 ! Plot on deformed conf.
plot,fill ! Plot filled mesh
plot,unde,,1 ! Plot on undefor. conf
plot,outl ! Plot outlined mesh
--> feappv appears an error on screen:
**ERROR** No 3d-Outline
Feappv can show only the new deformation shape, and cannot show the undefomation shape.
2. After solving the problem, I use the command: plot mesh
Then Feappv is interrupted immediately.
-
Thanks for find another error. I have verified the error and will issue a patch shortly.
-
The "ERROR" message is informing you that FEAPpv does not have 3-d outline capability, only FEAP has this feature. We will make the message clearer.
-
An error existed in xcompp.f in the plot directory. Try replacing with the attached file and rerun problems to see if all errors are now corrected.
Also instead of plot outl you could use plot mesh and at least follow where the deformed position is.
-
Thank you very much for the file. With your file Feappv works now with plot undefo + plot mesh. However I cannot plot elem. The program is interrupted again at line 110 in pltord.f.
-
All the problems you are finding relate to the flag we set for element faces that are not visible. For this array we set the material number to be negative. You are locating places where the material number is negative but we do not ignore the element. The current one is in 'pltelm.f' which is located in both the /unix and the /windows directories. Attached is the corrected one for the /windows directory. The but is in the /unix one too and you can correct in the same way as in the attached file.
Note the bug is in FEAP too and needs to be corrected in the same two directories.
If you are in perspective view, the element numbers printed will be face numbers that we construct and have nothing to do with the real element numbers you described.
-
Thank you very much for the file and your explanation. Feappv works fine now.