Hello everybody, i am trying to output my stresses by using a similar approach like in Step 18:
*Vector<double> <https://www.dealii.org/developer/doxygen/deal.II/classVector.html> norm_of_stress(triangulation <https://www.dealii.org/developer/doxygen/deal.II/p4est__wrappers_8cc.html#ace00f2f80d9780ef9aa1007e1c22c6a4>.n_active_cells());* * {* *for (auto &cell : triangulation <https://www.dealii.org/developer/doxygen/deal.II/p4est__wrappers_8cc.html#ace00f2f80d9780ef9aa1007e1c22c6a4>.active_cell_iterators())* *if (cell->is_locally_owned())* * {* *SymmetricTensor<2, dim> <https://www.dealii.org/developer/doxygen/deal.II/classSymmetricTensor.html> accumulated_stress;* *for (unsigned int q = 0; q < quadrature_formula.size(); ++q)* * accumulated_stress +=* *reinterpret_cast<PointHistory<dim> *>(cell->user_pointer())[q]* * .old_stress;* * norm_of_stress(cell->active_cell_index()) =* * (accumulated_stress / quadrature_formula.size()).norm <https://www.dealii.org/developer/doxygen/deal.II/namespaceLocalIntegrators_1_1Divergence.html#a8bcfc37d2a2be8faa18628a601ecf112>();* * }* *else* * norm_of_stress(cell->active_cell_index()) = -1e <https://www.dealii.org/developer/doxygen/deal.II/namespacePhysics_1_1Elasticity_1_1Kinematics.html#a9587d5229555daa5b1fa1ba2f8a40adb>+20;* * }* * data_out.add_data_vector <https://www.dealii.org/developer/doxygen/deal.II/classDataOut__DoFData.html#ac43d3b1f6e67424f36e474627fe8e401>(norm_of_stress, "norm_of_stress");* I have implemented my problem like in Step 44 by using the Workstream object. How can i implement this part of code in my case ? How can i call my "old_stress" from the PointHistory for each cell ? -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/9acd761a-b83f-484b-be89-d2ac7e186b15n%40googlegroups.com.
