Thank you Jean-Paul. I did not hear about L2-projection method before, but it seems that the L2-projection method also evaluates the stress tensor at the quadrature points and then finds the cell-averaged stress. The only difference compared with the heuristic method is that it solves AX = Y where A is mass matrix, Y is integration of cell-averaged stress component and X is the nodal stress component that we want. Is that right?
On Thursday, August 24, 2017 at 12:37:33 PM UTC-4, Jean-Paul Pelteret wrote: > > There's also the MeshWorker > <https://www.dealii.org/8.5.0/doxygen/deal.II/namespaceLocalIntegrators_1_1L2.html> > > framework which has some integrators that can help. But I don't know how > that works so I can't help you further other than to point it out to you. > > On Thursday, August 24, 2017 at 6:35:18 PM UTC+2, Jean-Paul Pelteret wrote: >> >> I can't think of any off of the top of my head, but I think that >> MatrixCreator::create_mass_matrix() >> <https://www.dealii.org/8.5.0/doxygen/deal.II/namespaceMatrixCreator.html#aa95a194001f23f96327fb7a630c6b036> >> plus >> VectorTools::create_right_hand_side() >> <https://www.dealii.org/8.5.0/doxygen/deal.II/namespaceVectorTools.html#ad18037ddbd9cc65a77ffca01ddc347f5> >> >> collectively assemble what you're wanting. Its really not difficult to do >> by hand - you should just give it a go! >> > -- 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.