Hi On Thursday, December 22, 2016 at 5:47:24 AM UTC+1, cmha...@gmail.com wrote: > > I have gone through some of the tutorials previously and thought it was > about time I actually tried to implement something on my own in the > library. > > I am attempting to implement a finite strain code using the F-Bar method. > Essentially this is a single field formulation which requires the > calculation the of the deformation gradient at the element center. The > trick here is though that the element centroid isn't a vertex. > > My question is, is there a way to "interpolate" the displacement gradient > at the element center without explicitly making it a quadrature point or > will I have to implement my own quadrature formula for this method? >
Yes, employ a custom quadrature formula with a single point with coordinates {0.5, 0.5, 0.5} to feed it to FEValues. There is no reason to figure out how to interpolate things, just use Quadrautre<dim> and FEValues<dim> together to evaluate your solution at quadrature point(s). Regards, Denis. > > Thank you for any assistance that can be provided. > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.