Your stress-strain response is nearly linear but you have defined it by specifying every point in *Plastic. The .msg file obviously indicates a lot of yielding. (review with https://msgfile.info)
This will never solve efficiently. I'm not sure what you're trying to simulate but you should...
Use Tools->Job Diagnostics or msgfile.info to find what node and issue is preventing convergence. After you understand the problem then you can fix it. If you have no idea what is wrong then you (and we) are only guessing.
mm-tonne-sec is a very common unit system but many others are also possible as long as you are consistent. This tool can be very helpful to get the correct units defined: https://msgfile.info/fea-units.html
You cannot expedite convergence until you know the cause of the slow convergence. The web tool https://msgfile.info can help summarize issues from the .msg file so you can have some ideas. As FEA way suggested, sometimes Explicit is a more efficient option for buckling but it brings its own...
I suggest defining a spherical partition or surface near the core of the ball and use that in the definition of your rigid body or coupling. Having the rigid elements small and far away from the interaction you're interested in will minimize their effect on the results.
This happens because your mesh is so fine. You can disable the contact thickness reduction if you exclude self contact and perimeter as described here:
https://www.cati.com/blog/2019/02/abaqus-contact-quirks-automatic-shell-thickness-reduction-and-you/
Attach the parts with *Tie or merge the nodes but this won't tell you the load transmitted.
You can use Integrated Output to get the total force transmitted by the weld in .odb file. *Section print is similar for .dat file.
If you're post-processing with Abaqus CAE then you could also use Free...
The part.faces is a FaceArray object. I'm sure you can find a FaceArray method that will help in the Scripting manual. In the worst-case you can iterate through each face yourself and compare centroids, normals, sizes, etc.
The Abaqus solver is using all upper case as a method to make the names case-insensitive. To preserve case you may put "double quotes" around your names in the .inp file.
Abaqus CAE does not permit "double quotes" in names but it will automatically add them if your name contains spaces so you...
TF2 will be the sum of concentrated force you specified (CF2) and reaction forces due to boundary conditions (RF2). RF2 should be all zero, as you are showing, since you don't have a boundary condition in that direction. Please take a look at TF2 or CF2 instead.
It's not completely clear what you are trying to do but it sounds like this would be much simpler if you just write restart data on the first analysis and read it on the following analysis. The plastic strain and everything else will continue as if the steps of the second analysis were a...
I think the main problem is that the "S" in subroutine on your first line should be in column 7, not column 1.
Google "fortran fixed format" to learn more about it.
What I meant was you can subtract the tool RF from model without cutting anything from the tool RF when it is cutting your workpiece.
And/or the RF constraining the workpiece should be equal and opposite to the cutting force, right? If you don’t have it at a ref point then just make a kinematic...