Polycrystal on deal ii

46 views
Skip to first unread message

A.Z Ihsan

unread,
Apr 7, 2020, 4:30:34 AM4/7/20
to deal.II User Group
Hi, 

Would anyone give me a hint on how to implement the polycrystal problem, i.e. more than one grain boundary?
I already have the tensor mechanics solver written in deal.ii, but now it needs to be extended solveing also polycrystal problem. 

THank you. 

Best, 
ihsan

Wolfgang Bangerth

unread,
Apr 7, 2020, 12:17:59 PM4/7/20
to dea...@googlegroups.com
On 4/7/20 2:30 AM, A.Z Ihsan wrote:
>
> Would anyone give me a hint on how to implement the polycrystal problem, i.e.
> more than one grain boundary?
> I already have the tensor mechanics solver written in deal.ii, but now it
> needs to be extended solveing also polycrystal problem.

Ihsan,
that's a rather unspecific question. What have you tried already, and what
step specifically is it that you have trouble with?

Best
WB


--
------------------------------------------------------------------------
Wolfgang Bangerth email: bang...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/

A.Z Ihsan

unread,
Apr 7, 2020, 1:36:08 PM4/7/20
to deal.II User Group
Hi Wolfgang, 

alright, suppose i have a box that divided into two regions, and each region has different elastic constant (anisotropy material). 
Perhaps the idea is when calculating the cell_matrix, inside the for-loop there is a function for choosing in which elastic constant to be used depending on the region. 

I already have the elastic solver written specifically for one elastic constant.  
Do you have any hint where should i start? 

BR, 
Ihsan

Andrew McBride

unread,
Apr 7, 2020, 1:41:42 PM4/7/20
to dea...@googlegroups.com
In the classical continuum approach to plasticity there is no grain boundary. The elastic anisotropy of the crystal is captured at the Gauss point level. So you would start by associating different materials properties (elastic / plastic) to Gauss points depending on their spatial location. 

Have a look at the CellDataStorage class and how it is used in step-44.

Best
Andrew

-- 
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+un...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/aa866c2c-6d5b-4798-a3df-abcb6a811628%40googlegroups.com.

A.Z Ihsan

unread,
Apr 7, 2020, 1:59:55 PM4/7/20
to deal.II User Group
Hi Andrew, 

Yes, thats what i want, just simply put the elastic anisotropy at the Gauss point level. 
Thank you for the answer. 

Best, 
Ihsan
Andrew

To unsubscribe from this group and stop receiving emails from it, send an email to dea...@googlegroups.com.

Wolfgang Bangerth

unread,
Apr 7, 2020, 2:16:10 PM4/7/20
to dea...@googlegroups.com
On 4/7/20 11:59 AM, A.Z Ihsan wrote:
>
> Yes, thats what i want, just simply put the elastic anisotropy at the Gauss
> point level.

step-6 already does that: It uses different material parameters depending on
the location of the quadrature point.

Best
W.
Reply all
Reply to author
Forward
0 new messages