FEAP User Forum

FEAP => Programming => Topic started by: K.Li on June 15, 2015, 07:48:29 AM

Title: Ordering of Gauss points
Post by: K.Li on June 15, 2015, 07:48:29 AM

Dear FEAP users:

It seems the ordering (numbering) of Gauss points in a 8-node brick element doesn't follow element node numbering in FEAP. The first Gauss point seems to be located near the 5th node of the element. This is different from other FE software. Could anyone please tell me how this ordering is determined or point me to the routine where this was defined? Can this ordering be changed for a specific element? 


Thanks a lot!

Title: Re: Ordering of Gauss points
Post by: Prof. R.L. Taylor on June 15, 2015, 08:16:46 AM

The quadrature order is defined in int3d.f (/program directory).

For the 2x2x2 quadrature the points would follow node numbering for the 8-node brick by setting:
        s(3,i    ) = -g
        s(3,i+4) =  g

Title: Re: Ordering of Gauss points
Post by: K.Li on June 15, 2015, 09:03:15 AM

Thank you for your help, Prof Taylor!