Jean-Paul
Any examples programming-wise? Like a deal.II manual/tutorial...
Jie
On Thursday, August 24, 2017 at 11:27:40 AM UTC-4, Jean-Paul Pelteret wrote:
>
> Hi Jie,
>
> Sure, we used it in this paper
>
>> @Article{vogel2014a-preprint,
>> author = {Vogel, F. and Pelteret, J-P. V. and Kaessmair, S. and
>> Steinmann, P.},
>> title = {Magnetic force and torque on particles subject to a
>> magnetic field},
>> journal = {European Journal of Mechanics A/Solids},
>> year = {2014},
>> volume = {48},
>> pages = {23--31},
>> month = {November--December}
>> }
>
>
> but, to be honest, almost any finite element textbook should have a blurb
> on it. E.g. page 227 of
>
>> @Book{Hughes2000a,
>> Title = {The Finite Element Method: Linear Static
>> and Dynamic Finite Element Analysis},
>> Author = {Hughes, T. J.},
>> Publisher = {Dover Publications Inc.},
>> Year = {2000},
>> Address = {New York, USA},
>> Note = {ISBN: 978-0486411811},
>> }
>
>
> Regards,
> Jean-Paul
>
> On Thursday, August 24, 2017 at 4:57:23 PM UTC+2, Jie Cheng wrote:
>>
>> Hi Jean-Paul
>>
>> Thanks for your help! The specific add_data_vector() function you shared
>> does solve my problem! For your second suggestion on L2-smoothing, I'd like
>> to give it a try. Can you please point me to some detailed references?
>> Examples or documentations would be great!
>>
>> Jie
>>
>> On Thursday, August 24, 2017 at 1:53:36 AM UTC-4, Jean-Paul Pelteret
>> wrote:
>>>
>>> Dear Jie,
>>>
>>> I'm glad to hear that you managed to find a solution to your problem and
>>> for posting what it was. Two quick points:
>>> 1. There's actually another add_data_vector()
>>> <https://www.dealii.org/8.5.0/doxygen/deal.II/classDataOut__DoFData.html#a4bce7f46c9ad03ecdd71a86d47f8fb11>
>>>
>>> member function of DataOut that takes in another DoFHandler as an argument,
>>> so you could compute the stress fields using another DoFHandler, with
>>> another DataComponentInterpretation etc. This will help remedy your last
>>> concern, because you will have only as many components as you need.
>>> 2. You might want to consider solving for the smoothed stress fields via
>>> global L2-smoothing rather than manual averaging. To the best of my
>>> understanding, it should produce a more accurate representation of the
>>> stress field than your approach. This should also be a relatively
>>> inexpensive operation because solving the linear system only involves
>>> inverting a mass matrix.
>>>
>>> Regards,
>>> Jean-Paul
>>>
>>> On Thursday, August 24, 2017 at 5:13:29 AM UTC+2, Jie Cheng wrote:
>>>>
>>>> Update:
>>>>
>>>> I switched to my old trick of nodal averaging method: take the average
>>>> of the values at quadrature points, then add this average to the vertex
>>>> dofs. And if a vertex is shared by n cells, divide the values at the
>>>> related dofs by n.
>>>> I am happy with the results because they look quite smooth.
>>>>
>>>> The only concern is that: for a stress component, say Sxx, appears in
>>>> my resulting vtu files as Sxx_0, Sxx_1 (in 2D) because the variable Sxx is
>>>> a vector of size n_dofs as required by the add_data_vector function. Is
>>>> there a way to get rid of the unnecessary duplicates? Essentially I want
>>>> Sxx to be a vector of size n_vertices.
>>>>
>>>> Thank you
>>>> Jie
>>>>
>>>>
>>>>
--
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.
