FEAP User Forum

FEAP => General questions => Topic started by: Rahul on September 07, 2022, 01:02:12 AM

Title: User element problrem with feap8.6
Post by: Rahul on September 07, 2022, 01:02:12 AM
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
Title: Re: User element problrem with feap8.6
Post by: Prof. R.L. Taylor on September 07, 2022, 03:55:14 AM
Post one of the elements
Title: Re: User element problrem with feap8.6
Post by: Prof. R.L. Taylor on September 07, 2022, 05:26:15 AM
missing routines iom_delta , iom_det , iom_invert
Title: Re: User element problrem with feap8.6
Post by: Rahul on September 07, 2022, 05:27:59 AM
All these routines i copied in program folder of feap 8.6 .
Title: Re: User element problrem with feap8.6
Post by: Prof. R.L. Taylor on September 07, 2022, 05:41:55 AM
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.
Title: Re: User element problrem with feap8.6
Post by: Prof. R.L. Taylor on September 07, 2022, 05:52:03 AM
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.

Title: Re: User element problrem with feap8.6
Post by: Rahul on September 07, 2022, 06:41:04 AM
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
Title: Re: User element problrem with feap8.6
Post by: Prof. R.L. Taylor on September 07, 2022, 07:49:06 AM
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.
Title: Re: User element problrem with feap8.6
Post by: Rahul on September 07, 2022, 10:18:20 AM
Thank you very much !!
Now it's working fine on my pc also.