Zhang Shuaihao

We still need to revise it for better this issue. This method is quite good, the result is stable even when time t=100. Is there any reference about this method?...

> We still need to revise it for better this issue. > Sure, I will firstly check the current version in detail.

This is the method we discussed just now. We only calculate the inverse of shear strain for `particle i` when the return mapping is called, i.e., when `particle i` is...

> Are you sure that you have 1/2 in front of (\phi_i + \phi_j)? Yes, I do put 1/2 in the code. This actually doesn't matter because we have this...

> So the parameter \xi has not effect to the solution? The parameter `𝜉` do have effects. But the simulation still crash even I used a large `𝜉` . I...

> > Also, in the first part of the equation should it be \tild{\sigma_i} + \tild{\sigma_j}? Yes, this should be the shear stress after return mapping. This is a typo...

> I have checked you code and found, besides quite intervened coding, a serious consistent issue. That is. If I comment line https://github.com/Shuaihao-Zhang/SPHinXsys/blob/06d36e86df5301e93e3381d6e9ed49e9a845ff88/tests/user_examples/extra_src/shared/general_continuum/continuum_dynamics/continuum_dynamics_inner.cpp#L533 to skip return mapping and change https://github.com/Shuaihao-Zhang/SPHinXsys/blob/06d36e86df5301e93e3381d6e9ed49e9a845ff88/tests/user_examples/extra_src/shared/general_continuum/continuum_dynamics/continuum_dynamics_inner.cpp#L555...

> Also. in the TL version, we use full strain to obtain scaling factor. Please have a look. https://github.com/Xiangyu-Hu/SPHinXsys/blob/bbe70dc9411ab91c98b5c08fa3d8848b65112bef/src/shared/particle_dynamics/solid_dynamics/inelastic_dynamics.cpp#L44-L47C49 I will check this in detail tomorrow.

> If I do as you said by: return shear_stress_rate_elastic; Then I do not comments line 533 but move it before, Mat3d deviatoric_stress = shear_stress_3D_[index_i]; The simulation now should not...

> @Shuaihao-Zhang, the returning map you wrote is to return the deviatoric part. The inconsistency between the elastic and plastic behaviors may arise due to the incompressible assumption. Pls, check...