On 05/31/2017 05:00 PM, 'Seyed Ali Mohseni' via deal.II User Group wrote:
I did another trick: Increasing G_c has the same effect to achieve
purely elastic behavior. So, I chose a high enough G_c value and it
works!
I obtain the same results and the duh values grow correctly.
Excellent. This is indeed also possible. Very good idea.
Still I cannot understand fully why when cracking is initiated, there
is no increase of the strains except the strain in the loaded direction.
Maybe you are right, I try to compare a bigger example, but I already
compared two phase-field codes.
I have my own version with staggered scheme and it was not identical
to your results after crack initiation.
As mentioned before, in the elastic regime everything is identical and
fine.
I made the experience that there can be large differences between
staggered, quasi-monolithic (Heister/Wick) and
fully monolithic (I have some recent studies on this).
Moreover, sometimes it makes a difference how you impose the
irreversibility condition: via strain history or
penalization.
The only thing I changed a bit is the integration from quadratic to
linear, but I am not sure, if seeting quadrature_formula(2) alone is
enough for linear shape functions. You also had
face_quadrature_formula or this Lobatto integration.
Would you be so kind and explain what lines I have to change in order
to obtain linear shape functions and not quadratic like set in the
original version?
You need to change the number in FE_Q( ....) for another degree. But in
our code,
it should be already linear if I remember well.
Best regards,
Thomas
How many changes and the code would help a lot.
Just to be sure I did it correct.
Thank you.
Best,
Seyed Ali
--
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
<mailto:dealii+unsubscr...@googlegroups.com>.
For more options, visit https://groups.google.com/d/optout.
--
++--------------------------------------------++
Dr. Thomas Wick
Maitre de conferences / Assistant Professor
Centre de Mathematiques Appliquees (CMAP)
Ecole Polytechnique
91128 Palaiseau cedex, France
Email: thomas.w...@cmap.polytechnique.fr
Phone: 0033 1 69 33 4579
www: http://www.cmap.polytechnique.fr/~wick/
++--------------------------------------------++
--
--
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.
