FEAP User Forum

FEAP => Programming => Topic started by: JStorm on June 11, 2020, 05:31:54 AM

Title: segmantation fault caused by plast_apl
Post by: JStorm on June 11, 2020, 05:31:54 AM

The return mapping algorithm in the implementation of Hill's plasticity model (plast_apl, FEAP 8.4) raises an error, when the maximum iteration number is reached (100). Than a message is printed out and plstop is called. But somehow FEAP is not terminating. I get the same error message again and again until the segmentation fault is raised. Forcing the termination by executing "stop" instead of calling plstop works as expected.

Title: Re: segmantation fault caused by plast_apl
Post by: JStorm on June 11, 2020, 05:42:33 AM

some more informations from the debugger:

• segmentation fault is caused by stack overflow due to extensive recursion
• recursion is caused by "last tplot" in plstop, where ptimpl is called and then formfe, etc.

Title: Re: segmantation fault caused by plast_apl
Post by: JStorm on June 11, 2020, 05:44:47 AM

found the solution, due to comments in the FEAP code  ;): prob_on = .true. musst be set before calling ptimpl in plstop.

Title: Re: segmantation fault caused by plast_apl
Post by: Prof. R.L. Taylor on June 11, 2020, 08:41:49 AM

Thank you