Dear delii members
While modifying the material model in the step-44.cc, implementation of the
function *get_Jc_vol(). *I found the implementation in the code is
different than the formula provided. The worse part is that the solution
stops converging after correcting the implementation. I will be really
grateful if some one can hint me as, I am stuck between convergence and
divergence of material models from quite some time.
The given implementation is:
SymmetricTensor<4, dim>
<https://www.dealii.org/current/doxygen/deal.II/classSymmetricTensor.html>
get_Jc_vol() const
{
return p_tilde * det_F *
(Physics::Elasticity::StandardTensors<dim>::IxI
<https://www.dealii.org/current/doxygen/deal.II/classPhysics_1_1Elasticity_1_1StandardTensors.html>
-
(2.0 * Physics::Elasticity::StandardTensors<dim>::S
<https://www.dealii.org/current/doxygen/deal.II/classPhysics_1_1Elasticity_1_1StandardTensors.html>
))
}
The implementation as per formula, should have been:
SymmetricTensor<4, dim> get_Jc_vol() const
{
return det_F * ( (get_dPsi_vol_dJ() +
det_F*get_d2Psi_vol_dJ2() )
*Physics::Elasticity::StandardTensors<dim>::IxI
- (2.0 * get_dPsi_vol_dJ() *
Physics::Elasticity::StandardTensors<dim>::S) );
}
Thank you,
Regards,
Navneet R
--
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 [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/CAK9McD0Rm7GotVhecV4%3D2ExxGcZPUtPLVQNVi652skkXYG0mOQ%40mail.gmail.com.
