Simulation of the deep-drawing process

Hi everyone,
I’m new to PrePoMax. I’ve been working through the example problems to get familiar with the program, and recently I started trying to set up my own simulation.

My goal is to use a deep-drawing process to test PrePoMax’s capabilities for simulating contact and elastoplastic material behavior. To validate the simulation, I based my setup on the paper:
J. Danckert, “Experimental investigation of a square-cup deep-drawing process,” Journal of Materials Processing Technology 50 (1995) 375–384. (Redirecting)

Unfortunately, the simulation failed to run.

Here are the details of my setup:

  • *Material:
    • Mild steel
    • Young’s modulus: 206 GPa
    • Poisson’s ratio: 0.3
    • Swift hardening law: σ = 565.32(0.007117 + ε)^0.2589 (MPa)
  • Boundary conditions:
    • Symmetry: 1/4 model with symmetry boundary conditions
    • Load: 19.6 kN applied to the blank holder
    • Punch displacement: 15 mm and 40 mm (two separate simulations)
  • Contact:
    • Friction coefficient (μ): 0.144
    • Default surface-to-surface contact settings (because i dont have more info)

I used shell elements for all bodies. The tools are modeled as rigid, and using shell elements allowed me to achieve a finer mesh at the contact surfaces without an excessive number of elements.

If anyone has advice on why the simulation might not be running or suggestions for improving the setup, I would really appreciate it!

I would load the files I used for the simulation below but i cant (im new user).

Any error messages ?

Try now.

There are the files as Zip:
DeepSquare.zip (4.9 MB)

Files:

  • ModeloShell.pmx – Simulation file (this file has the error)
  • squarecup (MeshGod) – CAD file modeled as a solid body
  • squarecup (ShellGod) – CAD file modeled as shell elements (used for ModeloShell.pmx)
  • epsilon_sigma.csv – Data points for the Swift hardening law
  • danckert1995.pdf – Scientific paper reference

The error i get is:

04/28/2025 16:52:19

########   Starting run step number: 1   Increment number: 1   ########

Running command: C:\Users\Jose Ocampo\Desktop\PrePoMax v2.3.0\Solver\ccx_dynamic.exe ModeloShell

************************************************************

CalculiX Version 2.22, Copyright(C) 1998-2024 Guido Dhondt
CalculiX comes with ABSOLUTELY NO WARRANTY. This is free
software, and you are welcome to redistribute it under
certain conditions, see gpl.htm

************************************************************

You are using an executable made on Sun Aug  4 19:44:24     2024

  The numbers below are estimated upper bounds

  number of:

   nodes:      1839388
   elements:        52568
   one-dimensional elements:            0
   two-dimensional elements:        24443
   integration points per element:           27
   degrees of freedom per node:            3
   layers per element:            1

   distributed facial loads:            0
   distributed volumetric loads:            0
   concentrated loads:            4
   single point constraints:      2666380
   multiple point constraints:      3156457
   terms in all multiple point constraints:     20112001
   tie constraints:            3
   dependent nodes tied by cyclic constraints:            0
   dependent nodes in pre-tension constraints:            0

   sets:           37
   terms in all sets:       415934

   materials:            2
   constants per material and temperature:            8
   temperature points per material:            1
   plastic data points per material:          201

   orientations:        24443
   amplitudes:            9
   data points in all amplitudes:            9
   print requests:            0
   transformations:            0
   property cards:            0

 *WARNING reading *FRICTION: stick slope
          must be strictly positive
          the following default will be used:   103000.00000000000     
          the user is advised to analyze the results
          carefully and, if possible, to come up with
          a experimentally based stick slope
 *WARNING reading *FRICTION. Card image:
          0.144


 STEP            1

 Static analysis was selected

 Nonlinear material laws are taken into account

 Newton-Raphson iterative procedure is active

 Decascading the MPC's

 *ERROR in cascade: the DOF corresponding to 
 node 1 in direction 1 is detected on the 
 dependent side of a MPC and a SPC


 Job failed - no results exist.

Process elapsed time:       9.24 s

¡Thanks for reply!

This error indicates overconstraint. I see that you’ve already been checking the location of node1. The problem is that you have both rigid body constraint and symmetry BC applied to the same nodes. Remove rigid parts’ edges from symmetry BCs and you won’t get this error.

If I do that, can I simulate it as a symmetric body?

I deactivated the symmetry boundary condition and encountered the following error:

04/28/2025 19:15:35

########   Starting run step number: 1   Increment number: 1   ########

Running command: C:\Users\Jose Ocampo\Desktop\PrePoMax v2.3.0\Solver\ccx_dynamic.exe ModeloShell

************************************************************

CalculiX Version 2.22, Copyright(C) 1998-2024 Guido Dhondt
CalculiX comes with ABSOLUTELY NO WARRANTY. This is free
software, and you are welcome to redistribute it under
certain conditions, see gpl.htm

************************************************************

You are using an executable made on Sun Aug  4 19:44:24     2024

  The numbers below are estimated upper bounds

  number of:

   nodes:      1839388
   elements:        52568
   one-dimensional elements:            0
   two-dimensional elements:        24443
   integration points per element:           27
   degrees of freedom per node:            3
   layers per element:            1

   distributed facial loads:            0
   distributed volumetric loads:            0
   concentrated loads:            4
   single point constraints:      2653148
   multiple point constraints:      3156457
   terms in all multiple point constraints:     20112001
   tie constraints:            3
   dependent nodes tied by cyclic constraints:            0
   dependent nodes in pre-tension constraints:            0

   sets:           37
   terms in all sets:       415760

   materials:            2
   constants per material and temperature:            8
   temperature points per material:            1
   plastic data points per material:          201

   orientations:        24443
   amplitudes:            7
   data points in all amplitudes:            7
   print requests:            0
   transformations:            0
   property cards:            0

 *WARNING reading *FRICTION: stick slope
          must be strictly positive
          the following default will be used:   103000.00000000000     
          the user is advised to analyze the results
          carefully and, if possible, to come up with
          a experimentally based stick slope
 *WARNING reading *FRICTION. Card image:
          0.144


 STEP            1

 Static analysis was selected

 Nonlinear material laws are taken into account

 Newton-Raphson iterative procedure is active

 Decascading the MPC's

 *INFO in cascade: linear MPCs and
       nonlinear MPCs depend on each other
       common node: 33902 in direction 1

 increment 1 attempt 1 
 increment size= 1.000000e-05
 sum of previous increments=0.000000e+00
 actual step time=1.000000e-05
 actual total time=1.000000e-05

 iteration 1

 Number of contact spring elements=1725987

 Decascading the MPC's

 *INFO in cascade: linear MPCs and
       nonlinear MPCs depend on each other
       common node: 33902 in direction 1

 Determining the structure of the matrix:

 Using up to 1 cpu(s) for setting up the structure of the matrix.
 number of equations
 625237
 number of nonzero lower triangular matrix elements
 43408170

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000006
 time avg. forc= 0.000006
 largest residual force= 0.005703 in node 199236 and dof 1
 largest increment of disp= 1.502096e-04
 largest correction to disp= 1.502096e-04 in node 136960 and dof 3

 no convergence

 iteration 2

 Number of contact spring elements=1603821

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=0.289300

 Using up to 1 cpu(s) for the stress calculation.

 average force= 319866406287833628672.000000
 time avg. forc= 319866406287833628672.000000
 largest residual force= 2060825728400171141693440.000000 in node 5381 and dof 3
 largest increment of disp= 5.880000e+12
 largest correction to disp= 5.880000e+12 in node 200597 and dof 3

divergence allowed: residual force too large
 divergence; the increment size is decreased to 1.000000e-05
 the increment is reattempted


 reducing the constant stiffnesses by a factor of 100 

 increment 1 attempt 2 
 increment size= 1.000000e-05
 sum of previous increments=0.000000e+00
 actual step time=1.000000e-05
 actual total time=1.000000e-05

 iteration 1

 Number of contact spring elements=1725987

 Decascading the MPC's

 *INFO in cascade: linear MPCs and
       nonlinear MPCs depend on each other
       common node: 33902 in direction 1

 Determining the structure of the matrix:

 Using up to 1 cpu(s) for setting up the structure of the matrix.
 number of equations
 625237
 number of nonzero lower triangular matrix elements
 43408170

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000005
 time avg. forc= 0.000005
 largest residual force= 0.000024 in node 48546 and dof 2
 largest increment of disp= 1.505522e-04
 largest correction to disp= 1.505522e-04 in node 153335 and dof 3

 no convergence

 iteration 2

 Number of contact spring elements=1584983

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=0.260629

 Using up to 1 cpu(s) for the stress calculation.

 average force= 291458381319277215744.000000
 time avg. forc= 291458381319277215744.000000
 largest residual force= 1856587901966470302662656.000000 in node 5381 and dof 3
 largest increment of disp= 5.880000e+12
 largest correction to disp= 5.880000e+12 in node 200597 and dof 3

divergence allowed: residual force too large
 no convergence

 iteration 3

 Number of contact spring elements=1558035

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=0.250000

 Using up to 1 cpu(s) for the stress calculation.

 average force= 405499279103936233472.000000
 time avg. forc= 405499279103936233472.000000
 largest residual force= 1939167864791244543623168.000000 in node 15987 and dof 3
 largest increment of disp= 9.649544e+18
 largest correction to disp= 9.649544e+18 in node 3 and dof 1

divergence allowed: residual force too large
 divergence; the increment size is decreased to 2.500000e-06
 the increment is reattempted


 *ERROR: increment size smaller than minimum
 best solution and residuals are in the frd file



Process elapsed time:       5932.641 s

You should leave the symmetry BCs only on the deformable part. Rigid parts don’t need them - all their DOFs are constrained by the rigid body constraint and you should apply BCs only directly to it.

This is typical error message indicating non-convergence. I was expecting it in this case. It’s very common in nonlinear analyses, especially those involving contact (even worse if there’s friction) and large deformations. There’s no single universal solution, you should examine the available results prior to non-convergence if available (check the deformed shape and output values) and try debugging it. Disable some advanced features, see if it works and so on.

When it comes to contact (as I said, that’s the most common cause of non-convergence), you should remember about the rules of master-slave assignment (master should be the rigid/stiffer and coarser one) and try playing with the contact stiffness setting - default is often too stiff.

In this case, I don’t have any previous results. The simulation fails before convergence, so I have to change the contact parameters first, and then continue from there.

One more thing — are there too many nodes? Maybe the mesh is not good enough?

In this case, there are three contact pairs: square matrix, square blank holder, and square punch.
In all of these contacts, the square part is the master region. Is that correct?

The main problem here is rather different. Since it doesn’t converge from the beginning, it’s likely due to initial rigid body motions before contact is established. Normally, such simulations are done (e.g. in Abaqus) with the dynamic explicit solver where RBMs are not a problem. But in CalculiX, dynamic explicit simulations are very limited and problematic so you usually have to resort to regular static analyses. To avoid such issues, you could try dividing the analysis into two steps and establishing contact in the first step by moving the blank holder with prescribed displacement instead of force (you can measure the reaction in the reference node to control the amount of prescribed displacement).

As I said, rigid/stiffer part should be master. Here only the blank (sheet metal part being formed) is not rigid so the choice is simple - it should be slave in contact pairs.

I understand, but I’ve watched your video " PrePoMax (CalculiX FEA) - Tutorial 43 - Snap-fit" (https://youtu.be/IsSsUDkVlpM), and in it, the two bodies are not in contact either.

i would try!

In the CAD model, the blank holder and the square appear to be in contact. Are you referring to the simulation instead?

Yes, but there I use prescribed displacement instead of force. Displacement control is one (very good) way to get rid of initial RBMs and aid convergence in nonlinear analyses not only when contact is the issue.

Yes, I mean establishing contact in the simulation - making contact elements work. Only when contact is properly established, it may prevent RBMs.

In the previous video, ‘PrePoMax (CalculiX FEA) - Tutorial 43 - Snap-fit,’ the beam is flexible and is defined as the master. Shouldn’t it be the slave, as you mentioned in the quote?

I’m used to ANSYS, where these things are more straightforward, so I never really had to worry about them. That’s the reason for so many questions — sorry about that.

That was an extremely problematic case, I had to play a lot with the settings and use some non-standard approaches to make this contact problem work in CalculiX (there was some buggy behavior). You can read about it here: Snap-fit contact snagging problem - CalculiX (official versions are on www.calculix.de, the official GitHub repository is at https://github.com/Dhondtguido/CalculiX).

1 Like

’ll try to explain what I did, and if something is wrong, feel free to let me know:

  1. Symmetry conditions are removed.
  2. I simulate one step first, where the displacement of the blank holder is imposed as a boundary condition (not as a force this time).
  3. The punch displacement is set to (0, 0, 0).
  4. The blank holder moves to (0, 0, -2 mm).
  5. The die (matrix) is fixed.
  6. The square part has no boundary conditions.
  7. Regarding contact, the stiffer/rigid parts (all tools) are defined as the master surfaces, and the square is the slave surface.
  8. A rigid body constraint is applied to the tools (die, punch, and blank holder).
  9. I don’t know how to properly define the force history output.

This time, the simulation produced results; however, the monitor displayed the following message:

04/29/2025 17:47:17

########   Starting run step number: 1   Increment number: 1   ########

Running command: C:\Users\Jose Ocampo\Desktop\PrePoMax v2.3.0\Solver\ccx_dynamic.exe ModeloShell

************************************************************

CalculiX Version 2.22, Copyright(C) 1998-2024 Guido Dhondt
CalculiX comes with ABSOLUTELY NO WARRANTY. This is free
software, and you are welcome to redistribute it under
certain conditions, see gpl.htm

************************************************************

You are using an executable made on Sun Aug  4 19:44:24     2024

  The numbers below are estimated upper bounds

  number of:

   nodes:      1839388
   elements:        52568
   one-dimensional elements:            0
   two-dimensional elements:        24443
   integration points per element:           27
   degrees of freedom per node:            3
   layers per element:            1

   distributed facial loads:            0
   distributed volumetric loads:            0
   concentrated loads:            0
   single point constraints:      2653152
   multiple point constraints:      3156457
   terms in all multiple point constraints:     20112001
   tie constraints:            3
   dependent nodes tied by cyclic constraints:            0
   dependent nodes in pre-tension constraints:            0

   sets:           36
   terms in all sets:       429802

   materials:            2
   constants per material and temperature:            8
   temperature points per material:            1
   plastic data points per material:          201

   orientations:        24443
   amplitudes:            6
   data points in all amplitudes:            6
   print requests:            1
   transformations:            0
   property cards:            0


 STEP            1

 Static analysis was selected

 Nonlinear material laws are taken into account

 Newton-Raphson iterative procedure is active

 Decascading the MPC's

 *INFO in cascade: linear MPCs and
       nonlinear MPCs depend on each other
       common node: 16726 in direction 1

 increment 1 attempt 1 
 increment size= 1.000000e-05
 sum of previous increments=0.000000e+00
 actual step time=1.000000e-05
 actual total time=1.000000e-05

 iteration 1

 Number of contact spring elements=2346243

 Decascading the MPC's

 *INFO in cascade: linear MPCs and
       nonlinear MPCs depend on each other
       common node: 16726 in direction 1

 Determining the structure of the matrix:

 Using up to 1 cpu(s) for setting up the structure of the matrix.
 number of equations
 625236
 number of nonzero lower triangular matrix elements
 44742993

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000017
 time avg. forc= 0.000017
 largest residual force= 0.015673 in node 196906 and dof 3
 largest increment of disp= 2.000000e-05
 largest correction to disp= 2.000000e-05 in node 11664 and dof 3

 no convergence

 iteration 2

 Number of contact spring elements=2255427

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=0.979427

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000018
 time avg. forc= 0.000018
 largest residual force= 0.000048 in node 74427 and dof 2
 largest increment of disp= 2.000000e-05
 largest correction to disp= 1.008850e-06 in node 154683 and dof 3

 no convergence

 iteration 3

 Number of contact spring elements=2238489

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.010000

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000018
 time avg. forc= 0.000018
 largest residual force= 0.000044 in node 73331 and dof 1
 largest increment of disp= 2.000000e-05
 largest correction to disp= 8.262513e-07 in node 137561 and dof 3

 no convergence

 iteration 4

 Number of contact spring elements=2211751

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.010000

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000018
 time avg. forc= 0.000018
 largest residual force= 0.000039 in node 74315 and dof 2
 largest increment of disp= 2.000000e-05
 largest correction to disp= 7.546043e-07 in node 143428 and dof 3

 no convergence

 iteration 5

 Number of contact spring elements=2221938

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.010000

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000018
 time avg. forc= 0.000018
 largest residual force= 0.000032 in node 189358 and dof 2
 largest increment of disp= 2.000000e-05
 largest correction to disp= 7.003347e-07 in node 137564 and dof 3

 no convergence

 iteration 6

 Number of contact spring elements=2188230

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.010000

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000018
 time avg. forc= 0.000018
 largest residual force= 0.000032 in node 74173 and dof 2
 largest increment of disp= 2.000000e-05
 largest correction to disp= 6.303740e-07 in node 143437 and dof 3

 no convergence

 iteration 7

 Number of contact spring elements=2181554

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.010000

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000018
 time avg. forc= 0.000018
 largest residual force= 0.000029 in node 74049 and dof 1
 largest increment of disp= 2.000000e-05
 largest correction to disp= 5.594512e-07 in node 137567 and dof 3

 no convergence

 iteration 8

 Number of contact spring elements=2171292

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.010000

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000019
 time avg. forc= 0.000019
 largest residual force= 0.000034 in node 73958 and dof 1
 largest increment of disp= 2.000000e-05
 largest correction to disp= 4.895065e-07 in node 137568 and dof 3

 no convergence

 iteration 9

 Number of contact spring elements=2179910

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.010000

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000018
 time avg. forc= 0.000018
 largest residual force= 0.000038 in node 73916 and dof 2
 largest increment of disp= 2.000000e-05
 largest correction to disp= 4.209594e-07 in node 143449 and dof 3

 no convergence

 iteration 10

 Number of contact spring elements=2164321

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.010000

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000019
 time avg. forc= 0.000019
 largest residual force= 0.000032 in node 73826 and dof 2
 largest increment of disp= 2.000000e-05
 largest correction to disp= 4.278042e-07 in node 150815 and dof 3

 no convergence

 iteration 11

 Number of contact spring elements=2157893

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.010000

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000019
 time avg. forc= 0.000019
 largest residual force= 0.000033 in node 73734 and dof 1
 largest increment of disp= 2.000000e-05
 largest correction to disp= 4.364700e-07 in node 141031 and dof 3

 no convergence

 iteration 12

 Number of contact spring elements=2154976

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.010000

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000019
 time avg. forc= 0.000019
 largest residual force= 0.000031 in node 192903 and dof 2
 largest increment of disp= 2.000000e-05
 largest correction to disp= 4.399200e-07 in node 149687 and dof 3

 no convergence

 iteration 13

 Number of contact spring elements=2149506

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.010000

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000019
 time avg. forc= 0.000019
 largest residual force= 0.000031 in node 73601 and dof 2
 largest increment of disp= 2.000000e-05
 largest correction to disp= 4.519246e-07 in node 149279 and dof 3

 no convergence

 iteration 14

 Number of contact spring elements=2145638

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.010000

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000019
 time avg. forc= 0.000019
 largest residual force= 0.000026 in node 73509 and dof 2
 largest increment of disp= 2.000000e-05
 largest correction to disp= 4.724667e-07 in node 148413 and dof 3

 no convergence

 iteration 15

 Number of contact spring elements=2142223

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.010000

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000019
 time avg. forc= 0.000019
 largest residual force= 0.000022 in node 73516 and dof 1
 largest increment of disp= 2.000000e-05
 largest correction to disp= 4.913104e-07 in node 146991 and dof 3

 no convergence

 iteration 16

 Number of contact spring elements=2140141

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.010000

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000019
 time avg. forc= 0.000019
 largest residual force= 0.000023 in node 73370 and dof 2
 largest increment of disp= 2.000000e-05
 largest correction to disp= 5.900948e-07 in node 137578 and dof 3

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 17

 Number of contact spring elements=2139123

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.010000

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000019
 time avg. forc= 0.000019
 largest residual force= 0.000018 in node 73372 and dof 2
 largest increment of disp= 2.000000e-05
 largest correction to disp= 4.316528e-07 in node 143989 and dof 3

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 18

 Number of contact spring elements=2138709

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.009319

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000019
 time avg. forc= 0.000019
 largest residual force= 0.000013 in node 73374 and dof 2
 largest increment of disp= 2.000000e-05
 largest correction to disp= 2.505242e-07 in node 137747 and dof 3

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 19

 Number of contact spring elements=2138528

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.001193

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000019
 time avg. forc= 0.000019
 largest residual force= 0.000012 in node 53830 and dof 2
 largest increment of disp= 2.000000e-05
 largest correction to disp= 2.770535e-07 in node 154916 and dof 3

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 20

 Number of contact spring elements=2138485

 Using up to 1 cpu(s) for the stress calculation.

 Using up to 1 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 1 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 1

 Using up to 1 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 1 cpu(s) for the stress calculation.

 average force= 0.000019
 time avg. forc= 0.000019
 largest residual force= 0.000000 in node 703 and dof 3
 largest increment of disp= 2.000000e-05
 largest correction to disp= 7.267564e-06 in node 183418 and dof 3

 convergence

The simulation is still running, but the first step is taking a long time to converge.

Mi matrix tools is rigid but in the result show deformation (very small), it is normal?

Do you mean the blank (sheet metal part being formed) ? It should have symmetry BCs. They need to be removed only from the rigid tools (die, punch and holder).

You can select a reference point (the one where you prescribe displacement) for RF history output.

You should change Tools → Settings → Calculix → Number of processors to the maximum number of CPU threads you can use on your computer to speed up the calculations.

Rigid parts can have displacements from rigid body motion. They just don’t have strains and stresses.

I do it, thank .

What should I use in total or global parameter?

Thanks!

Totals is only relevant if you select multiple nodes and want to sum the reaction forces from them. Here you have only a single node so no need to change that.

Global is only relevant if you use local coordinate systems in BC/load definitions and you don’t do it here so no need to change this one either.

1 Like

What do you thing about the following?:

04/30/2025 13:32:45

########   Starting run step number: 1   Increment number: 1   ########

Running command: C:\Users\JoseOcampo\Desktop\PrePoMax v2.3.0\Solver\ccx_dynamic.exe Modelo_result

************************************************************

CalculiX Version 2.22, Copyright(C) 1998-2024 Guido Dhondt
CalculiX comes with ABSOLUTELY NO WARRANTY. This is free
software, and you are welcome to redistribute it under
certain conditions, see gpl.htm

************************************************************

You are using an executable made on Sun Aug  4 19:44:24     2024

  The numbers below are estimated upper bounds

  number of:

   nodes:      1839388
   elements:        52568
   one-dimensional elements:            0
   two-dimensional elements:        24443
   integration points per element:           27
   degrees of freedom per node:            3
   layers per element:            1

   distributed facial loads:            0
   distributed volumetric loads:            0
   concentrated loads:            0
   single point constraints:      2661624
   multiple point constraints:      3156457
   terms in all multiple point constraints:     20112001
   tie constraints:            3
   dependent nodes tied by cyclic constraints:            0
   dependent nodes in pre-tension constraints:            0

   sets:           37
   terms in all sets:       414744

   materials:            2
   constants per material and temperature:            8
   temperature points per material:            1
   plastic data points per material:          201

   orientations:        24443
   amplitudes:            8
   data points in all amplitudes:            8
   print requests:            2
   transformations:            0
   property cards:            0


 STEP            1

 *INFO reading *STEP: nonlinear geometric
       effects are turned on

 Static analysis was selected

 Nonlinear material laws are taken into account

 Newton-Raphson iterative procedure is active

 Nonlinear geometric effects are taken into account

 Decascading the MPC's

 *INFO in cascade: linear MPCs and
       nonlinear MPCs depend on each other
       common node: 16726 in direction 1

 increment 1 attempt 1 
 increment size= 1.000000e-05
 sum of previous increments=0.000000e+00
 actual step time=1.000000e-05
 actual total time=1.000000e-05

 iteration 1

 Number of contact spring elements=2346243

 Decascading the MPC's

 *INFO in cascade: linear MPCs and
       nonlinear MPCs depend on each other
       common node: 16726 in direction 1

 Determining the structure of the matrix:

 Using up to 8 cpu(s) for setting up the structure of the matrix.
 number of equations
 623118
 number of nonzero lower triangular matrix elements
 44515588

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

 average force= 1288.609287
 time avg. forc= 1288.609287
 largest residual force= 12110031.866343 in node 5381 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

 no convergence

 iteration 2

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

 no convergence

 iteration 3

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

 no convergence

 iteration 4

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 13

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 39

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 40

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 41

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 42

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence
 iteration 47

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 48

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 49

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 50

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 51

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 52

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 53

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 54

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 55

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 56

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 57

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 58

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 59

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 60

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
maximum number of iterations for face-to-face contact reached
 the increment size is decreased to 1.000000e-05
 the increment is reattempted


 reducing the constant stiffnesses by a factor of 100 

 increment 1 attempt 2 
 increment size= 1.000000e-05
 sum of previous increments=0.000000e+00
 actual step time=1.000000e-05
 actual total time=1.000000e-05

 iteration 1

 Number of contact spring elements=2346243

 Decascading the MPC's

 *INFO in cascade: linear MPCs and
       nonlinear MPCs depend on each other
       common node: 16726 in direction 1

 Determining the structure of the matrix:

 Using up to 8 cpu(s) for setting up the structure of the matrix.
 number of equations
 623118
 number of nonzero lower triangular matrix elements
 44515588

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

 average force= 1288.609287
 time avg. forc= 1288.609287
 largest residual force= 12110031.866343 in node 5381 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

 no convergence

 iteration 2

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

 no convergence

 iteration 3

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

 no convergence

 iteration 4

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 5

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 6

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 7

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 8

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 9

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 10

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 11

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 27

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized
 no convergence

 iteration 32

 Number of contact spring elements=2338056

 Using up to 8 cpu(s) for the stress calculation.

 Using up to 8 cpu(s) for the symmetric stiffness/mass contributions.

 Using up to 8 cpu(s) for the asymmetric stiffness/mass contributions.

 Factoring the system of equations using the unsymmetric pardiso solver
 number of threads = 8

 Using up to 8 cpu(s) for the stress calculation.

line search factor=1.000000

 Using up to 8 cpu(s) for the stress calculation.

 average force= 872.619783
 time avg. forc= 872.619783
 largest residual force= 8040784.822977 in node 15987 and dof 3
 largest increment of disp= 0.000000e+00
 largest correction to disp= 0.000000e+00

divergence allowed: number of contact elements stabilized

My output history isn’t showing.