Linear and Powerlaw strain hardening question? #16449
-
Hi Moose Experts, I'm new to the MOOSE environment and trying to run some Tensor Mechanics examples. I need to verify analytical solution of this paper "Plastic Relaxation of the Transformation Strain Energy of a Misfitting Spherical Precipitate: Linear and Power-Law Strain Hardening" with moose numerical solution. Can anyone suggest me where should I start? Is this sufficient to start with this examples projects/moose/modules/tensor_mechanics/test/tests/recompute_radial_return/isotropic_plasticity_finite_strain.i and projects/moose/modules/combined/test/tests/power_law_hardening/PowrLawHardening.i ? Paper Attached in this discussion. Thanks for the assistance! |
Beta Was this translation helpful? Give feedback.
Replies: 4 comments 5 replies
-
If you are new to MOOSE, I suggest you start to read the Tensor Mechanics documentation, i.e. https://mooseframework.inl.gov/source/materials/IsotropicPowerLawHardeningStressUpdate.html |
Beta Was this translation helpful? Give feedback.
-
Hi Moose Experts, I could create the input file for the problem as mentioned in the last post (for isotropic plasticity ), however the solution is not converging. I have read this information here :https://mooseframework.inl.gov/modules/tensor_mechanics/Convergence.html , but could not understand the process of setting up the convergence criteria? Can someone please explain the convergence criteria in detail. Thanks in advanced, |
Beta Was this translation helpful? Give feedback.
-
Dear @naveeniisc I noticed you fixed the displacement only on one node, called boundary = 100 Can you try more standard boundary conditions for mechanics problems, [./z0_bot] [./y0_bot] [./x0_bot] Convergence problems are likely to appear if you don't constrain properly your system. |
Beta Was this translation helpful? Give feedback.
Dear @naveeniisc
I noticed you fixed the displacement only on one node, called boundary = 100
Can you try more standard boundary conditions for mechanics problems,
like constraining your cubic RVE in a corner,
something like:
[./z0_bot]
type = DirichletBC
variable = disp_z
boundary = back
value = 0.0
[../]
[./y0_bot]
type = DirichletBC
variable = disp_y
boundary = bottom
value = 0.0
[../]
[./x0_bot]
type = DirichletBC
variable = disp_x
boundary = left
value = 0.0
[../]
Convergence problems are likely to appear if you don't constrain properly your system.
Best Regards,
Nicolò Grilli
National University of Singapore