Dear community,
I have encountered an issue when modeling a cylindrical shell with S4R elements in Abaqus.
The model is shown below (the platform does not allow me to attach the .inp).
Model description:
-
Cylindrical shell (L = 52.501E-03 mm, D = 9.648E-03 mm, th = 0.001 mm) with a perfectly structured mesh quad mesh with elm size of 180.E-06 mm.
-
Compressible isotropic hyperelastic Yeoh material.
*Hyperelastic, yeoh
0.000615, 0.002488, -0.001183, 162., 0., 0. -
Static General step with NLGEOM=ON.
-
Axial displacement-controlled pure tension.
The response is expected to be purely membrane, with no bending.
The prescribed displacement 0.0315 mm is applied using an amplitude function:
1. From 0 s to 200 s, the displacement increases rapidly.
2. From 200 s to 1000 s, the displacement is held constant.
*Amplitude, name=Amp_paper, time=TOTAL TIME
1., 0.5, 10., 0.658333333, 22., 0.725, 44., 0.791666667
74., 0.866666667, 106., 0.933333333, 148., 0.975, 200., 1.
1000., 1.
Unexpectedly, the analysis diverges at approximately 720 s.
There is no additional loading applied after 200 s, the displacement increment is 0.
The model remains stable with perfectly uniform stress distribution until around 586 s.
Then, a localized stress concentration begins to develop spontaneously; eventually the solution diverges at about 720 s.
Has anyone encountered similar behavior with S4R shell elements or hyperelastic materials?
Any suggestions on possible causes would be greatly appreciated.
Thank you.


