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

CalculiX Version 2.20, Copyright(C) 1998-2022 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 Oct 15 13:18:38 BST 2023

  The numbers below are estimated upper bounds

  number of:

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

   distributed facial loads:        10980
   distributed volumetric loads:            3
   concentrated loads:            0
   single point constraints:        35196
   multiple point constraints:            1
   terms in all multiple point constraints:            1
   tie constraints:           20
   dependent nodes tied by cyclic constraints:            0
   dependent nodes in pre-tension constraints:            0

   sets:           64
   terms in all sets:       431734

   materials:           30
   constants per material and temperature:            8
   temperature points per material:            1
   plastic data points per material:            0

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

 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *SURFACE BEHAVIOR: the tension
       value at large clearances
       (only relevant for node-to-surface contact)
       must be strictly positive
       the following default will be used:   2857142.8571428573     
       the user is advised to analyze the results
       carefully
 *INFO reading *STEP: nonlinear geometric
       effects are turned on


 STEP            1

 Dynamic analysis was selected

 Newton-Raphson iterative procedure is active

 Nonlinear geometric effects are taken into account

 Decascading the MPC's

 Determining the structure of the matrix:
 Using up to 12 cpu(s) for setting up the structure of the matrix.
 number of equations
 1951308
 number of nonzero lower triangular matrix elements
 132163945

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

 Using up to 12 cpu(s) for the energy calculation.

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

 Factoring the system of equations using the symmetric pardiso solver
 number of threads = 12

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

 Using up to 12 cpu(s) for the energy calculation.

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

 iteration 1

 Number of contact spring elements=685424

 Determining the structure of the matrix:
maximal possible contact elements = 
1761143

 Using up to 12 cpu(s) for setting up the structure of the matrix.
 number of equations
 1951308
 number of nonzero lower triangular matrix elements
 138462736

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

 Using up to 12 cpu(s) for the energy calculation.

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

 Factoring the system of equations using the symmetric pardiso solver
 number of threads = 12

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

 Using up to 12 cpu(s) for the energy calculation.

 average force= 39.448875
 time avg. forc= 39.448875
 largest residual force= 30178.583259 in node 173253 and dof 3
 largest increment of disp= 2.956984e-03
 largest correction to disp= 2.956984e-03 in node 122409 and dof 3

 no convergence

 iteration 2

 Number of contact spring elements=1231318

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

 Using up to 12 cpu(s) for the energy calculation.

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

 Factoring the system of equations using the symmetric pardiso solver
 number of threads = 12

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

 Using up to 12 cpu(s) for the energy calculation.

 average force= 17.930798
 time avg. forc= 17.930798
 largest residual force= 6345.714156 in node 112433 and dof 3
 largest increment of disp= 2.956984e-03
 largest correction to disp= 5.702041e-04 in node 155312 and dof 3

 no convergence

 iteration 3

 Number of contact spring elements=926769

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

 Using up to 12 cpu(s) for the energy calculation.

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

 Factoring the system of equations using the symmetric pardiso solver
 number of threads = 12

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

 Using up to 12 cpu(s) for the energy calculation.

 average force= 22.218393
 time avg. forc= 22.218393
 largest residual force= 362.518010 in node 119582 and dof 3
 largest increment of disp= 2.956984e-03
 largest correction to disp= 6.825959e-05 in node 171785 and dof 3