FEAP User Forum

FEAP => nurbFEAP => Topic started by: dz on December 18, 2020, 06:16:21 AM

Title: Projection of stresses and strains in v8.6
Post by: dz on December 18, 2020, 06:16:21 AM

Dear professors and FEAP users,

I am new to iga feap, and am currently running some simple tests(with both v8.5 and v8.6 on Ubuntu).

In my case, it seems that the projection of stresses and strains in v8.6 gives results different to that in v8.5. i.e. with a symmetric setup and loading, I would expect the stresses to be symmetric, but this is only the case in v8.5. I compared the element subroutine(sld3d1) and the projection subroutine(slcn3d), and did not find any difference that I think would lead to the different results.

Attached are my input file and the plot results.

Any suggestion is very much appreciated.

Best,
dz

Title: Re: Projection of stresses and strains in v8.6
Post by: Prof. R.L. Taylor on December 18, 2020, 10:45:24 AM

The but is in ./iga/plot/nslcn3d.f

Change "100" to "125" on the call to invert for 'matnp' -- this allows IgA 3d problems with quadrature up to 5x5x5.

Thank you for reporting the bug.

Title: Re: Projection of stresses and strains in v8.6
Post by: dz on December 21, 2020, 12:49:16 AM

I changed the code according to your suggestion, now it works fine. Thank you for the clarification, Prof. Taylor.