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).
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
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.
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.
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.
’ll try to explain what I did, and if something is wrong, feel free to let me know:
Symmetry conditions are removed.
I simulate one step first, where the displacement of the blank holder is imposed as a boundary condition (not as a force this time).
The punch displacement is set to (0, 0, 0).
The blank holder moves to (0, 0, -2 mm).
The die (matrix) is fixed.
The square part has no boundary conditions.
Regarding contact, the stiffer/rigid parts (all tools) are defined as the master surfaces, and the square is the slave surface.
A rigid body constraint is applied to the tools (die, punch, and blank holder).
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.
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.
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.
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