Hi,
thanks in advance for your suggestions.
I have been working on a small strain version of the three-fields
formulation proposed in step 44. To this aim, I have been extended step 18,
including some non linear material models.
It is my understanding that there are two ways of assembling and solving
the system. One way simply extends the data structure of step 18, by using
a different FESystem and solving a larger linear system at the end.
The second way is using block data-structures and solving via Shur
complement. I aim at doing both, to get to learn.
The first strategy seems a little easier. I defined the FESystem following
step 44
fe(FE_Q<dim>(poly_degree), dim, // displacement
FE_DGPMonomial<dim>(poly_degree - 1), 1, // pressure
FE_DGPMonomial<dim>(poly_degree - 1), 1), // dilatation
and with some minor adjustments I got to solve the system. I am having a
little trouble in output data, since I want to separate displacements,
pressure, and dilatation from the solution vector and I am not sure how
this can be achieved.
This set of instructions:
DataOut<dim> data_out;
data_out.attach_dof_handler(dof_handler);
// displacement output
const std::vector<DataComponentInterpretation::DataComponentInterpretation
>
data_component_interpretation(dim, DataComponentInterpretation::
component_is_part_of_vector);
data_out.add_data_vector(system_solution,
std::vector<std::string> (dim, "displacement"),
DataOut<dim>::type_dof_data,
data_component_interpretation);
is inappropriate as it stands, since the system_solution contains u,p,J and
I want to extract u, p and J separately. Any hint?
Thanks
Alberto
--
Informativa sulla Privacy: http://www.unibs.it/node/8155
--
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.
