How can I get the global coordinate of Integrated Point

24 views
Skip to first unread message

Aquaman

unread,
Jul 17, 2019, 5:57:16 PM7/17/19
to deal.II User Group
Dear all

What I'm looking for is the global coordinate of the this Integrated Point which the stress, strain, residual, and tangent are calculated based on.

I guess a command like in Abaqus UMAT: COORD(3)

But How I can get the X,Y,Z components in DEALII?

Thanks in advance,

Jacob

Jean-Paul Pelteret

unread,
Jul 18, 2019, 1:18:23 AM7/18/19
to dea...@googlegroups.com
Dear Jacob,

You should have a look at step-8, for example. In assemble_system() some local material properties are computed by evaluating some functions at all cell quadrature points:

FEValues<dim> fe_values(fe, quadrature_formula, update_quadrature_points | …);

for (const auto &cell : dof_handler.active_cell_iterators())
{
  fe_values.reinit(cell);
  …
  lambda.value_list(fe_values.get_quadrature_points(), lambda_values);
}

So fe_values.get_quadrature_points() returns the real positions of the coordinates at the integration points.

Best,
Jean-Paul


--
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+un...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/f11ba808-7e70-4c90-882d-7978317665f8%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Reply all
Reply to author
Forward
0 new messages