GEOS-DEV / GEOS

GEOS Simulation Framework
GNU Lesser General Public License v2.1
210 stars 84 forks source link

bug about Implicit Small Strain Newmark method #3077

Open Jianming-Pu opened 5 months ago

Jianming-Pu commented 5 months ago

Describe the bug When I used the solid implicit dynamic method to calculate problems, the calculation did not converge. So I checked the relevant code. Here, you can see that detJ is not reassigned after its initial value of 0, so none of the following calculations seem to make sense.

https://github.com/GEOS-DEV/GEOS/blob/4696a656f7f586c824eb8106341e9af676981362/src/coreComponents/physicsSolvers/solidMechanics/kernels/ImplicitSmallStrainNewmark_impl.hpp#L109

I think there might be something wrong with the code.