Hi, everyone:
I am doing a solid mechanics problem (plane stress), and I have already got the displacement solution. Now I need to use them to compute stress distribution with the constitutive relation, namely: <https://lh3.googleusercontent.com/-VMQA7z7hj0U/V0EdaSD1iXI/AAAAAAAAACY/uC6BlstmQ5MA0dlu8xAEw3nNARDsRABuQCLcB/s1600/constitutive.png> So, I also need to compute shape gradient values at the nodal points, and then attach these stress solutions to the corresponding nodes when output results. But the problem here is that I need to know shape gradient values at triangulation nodal points. It seems that I can only use FEValues to achieve this, but FEValues can only give me shape gradient values at quadrature points. So I tried to create a second order QGaussLobatto object, so that the quadrature points are the same with nodal points, and the quadrature points of index 0, 1, 2 and 3 corresponds respectively to each node of one cell. However it's still not that easy to control the which node this quadrature object actually refers to when I want to compute the stress at a given triangullation node. Dose anyone have some suggestions or better ideas about this? Any help would be greatly appreciated! Best regards, David. -- 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 [email protected]. For more options, visit https://groups.google.com/d/optout.
