Hamed,
assuming that both DoFHandlers are based on the same triangulation, you can
do something like the following:
FEValues fe_stress_values(...);
stress_cell = dof_handler_stress.begin_active();
temperature_cell = dof_handler_temperature.begin_active();
for (; stress_cell != dof_handler_stress.end(); ++stress_cell,
++temperature_cell)
{
fe_stress_values.reinit(stress_cell);
fe_stress_values.get_function_values(solution_stress, local_stress);
//do all the temperature related stuff and use local_stress
}
Instead of having two iterators separated, you can also use
SynchronousIterators [1]. They arec partly used in step-35 for that purpose.
Best,
Daniel
[1] https://dealii.org/8.4.1/doxygen/deal.II/structSynchronousIterators.html
Am Dienstag, 7. Juni 2016 20:30:27 UTC+2 schrieb Hamed Babaei:
>
> Hello Daniel,
>
> Since although we can solve two elastic and heat equations separately, we
> still need to enter stress computed from elastic equation into the heat
> equation in every time step and keep updating stress field.
> I was wondering when we have two different DoFHandlers for each of the
> equations, how to use, for instance, stress field stored on a quadrature
> point related to elastic problem DofHandler, in order to solve
> stress-induced heat equation problem. In other words, I don't know how to
> exchange information between two equation.
>
> Thanks,
> Hamed
>
> On Friday, June 3, 2016 at 3:40:28 AM UTC-5, Daniel Arndt wrote:
>>
>> You're right, Hamed. There is FEValuesViews::Vector< dim, spacedim
>> >::symmetric_gradient [1]
>> and FEValuesViews::Vector< dim, spacedim >::gradient [2] instead. The
>> first should be exactly what you are looking for.
>>
>> Code that I have written normally uses multiple DoFHandlers if I solve
>> for multiple equations separately.
>>
>> Best,
>> Daniel
>>
>> [1]
>> https://dealii.org/8.4.1/doxygen/deal.II/classFEValuesViews_1_1Vector.html#a4e5dfbb49d284a368acac6ef5dd4b2ef
>>
>> <https://dealii.org/8.4.1/doxygen/deal.II/classFEValuesViews_1_1Vector.html>
>> [2]
>> https://dealii.org/8.4.1/doxygen/deal.II/classFEValuesViews_1_1Vector.html#ad7b4df64147989229f566eff541ddebc
>>
>> <https://dealii.org/8.4.1/doxygen/deal.II/classFEValuesViews_1_1Vector.html>
>>
>>
>>
>>
>>
>>
>> Am Freitag, 3. Juni 2016 01:33:34 UTC+2 schrieb Hamed Babaei:
>> Daniel,
>>
>> Since in the FEValuesViews::Vector< dim, spacedim > Class there isn't any
>> shape_grad_component function, I am not sure that I can update the
>> get_strain as I mentioned.
>> I was wondering if you have already written a code for a thermoelastic
>> problem in dealii, solving heat equation and elastic equation separately.
>>
>> Best,
>> Hamed
>>
>> On Thursday, June 2, 2016 at 3:31:48 PM UTC-5, Daniel Arndt wrote:
>> Hamed,
>>
>>
>> There shouldn't be a problem to do what you following step-22. In
>> particular, you can modify step-18::get_strain as you proposed.
>> Did you face any difficulties so far?
>>
>>
>> Best,
>> Daniel
>>
>
--
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.
