User element problrem with feap8.6

[Rahul]

Hello everyone,

I wrote some user elements. Those were working fine with feappv, but now I am using those element routines with feap 8.6, then it's not working.

Could someone please help me with this? What are the changes I need to make?

Thank you
Regards
Rahul Jangid

[Prof. R.L. Taylor]

Post one of the elements

[Prof. R.L. Taylor]

missing routines iom_delta , iom_det , iom_invert

[Rahul]

All these routines i copied in program folder of feap 8.6 .

[Prof. R.L. Taylor]

They are your routines, I cannot check your element without them. 

In general it is not good practice to add things like this to the main program but keep separate in your modules.

[Prof. R.L. Taylor]

I do not see any obvious reason element should not work in ver8.6.

Can you describe what happens?

Do all routines compile?

Are all inputs acknowledged in the output file?

I would suggest you run the problem in interactive mode.

Use commands like: EPRInt    to see the element stiffness produced.

Give one TANG to see what happens.

You may need to use a debugger to find what is the problem.

[Rahul]

Feappv solving the problem whereas feap 8.6 not solving the problem. (please check the attachment for the output files of the above-attached input file from both versions )

Yes, I compile all the programs.

please have look at the output files.(attached)

thank you

[Prof. R.L. Taylor]

The only error I see is the lack of setting "pstyp = 3" under isw.eq.1

You do not set the "inord" array in feap -- instead you tell the program that elements are 3-d solides using the pstyp parameter.  Since you have NOPRINT you do not see anything.

It runs on my computer and says it converged but I do not know if the answers are correct.

[Rahul]

Thank you very much !!
Now it's working fine on my pc also.