FEAP User Forum

FEAP => Input File Issues => Topic started by: Jacopo on March 25, 2021, 11:34:53 AM

Title: CSURface in FEAP86
Post by: Jacopo on March 25, 2021, 11:34:53 AM

Dear FEAP community,
I am experiencing an issue with the CSURface command in FEAP86.
I can not apply a uniform normal traction distribution over a straight edge in a 2D problem.
Issuing the command results in a wavy distribution, instead of constant. 
Please find below a MWE that reproduces the problem. I tested the same toy example in FEAP83 and FEAPpv51,
where the results are fine. I hope I am not missing something trivial,
Thanks for your help and consideration,
Jacopo

Title: Re: CSURface in FEAP86
Post by: Prof. S. Govindjee on March 25, 2021, 12:03:44 PM

There is a bug in the code; the distribution you see is an axis-symmetric distribution. 

To work around, remove your statements for global plane strain and you will get what you want (plus you problem will be plane strain  as that is the default in FEAP).

I'll post a true fix later.

Title: Re: CSURface in FEAP86
Post by: Jacopo on March 25, 2021, 02:01:32 PM

Dear Professor,
Thank you for your reply, now it works properly.
Jacopo

Title: Re: CSURface in FEAP86
Post by: Prof. R.L. Taylor on March 25, 2021, 06:12:41 PM

The bug is in ./program/pesurf.f at about line 238

Delete:   axifl = gtypfl