=========================================================================== ________ _________ _________ __ _________ | |\ | |\ | |\ | |\ / \\ | ____|| | ____|| | __ || |__|| | ___ || | |\___\| | |\___\| | |\_| || \_\| | // \ || | || | || | || | || __ | || | || | ||__ | ||__ | ||_| || | |\ | || | || | |\ | |\ | \\ | || | || | || | ___|| | ___|| | ___ || | || | || | || | |\__\| | |\__\| | |\__| || | || | || | || | || | || | || | || | || | || | || | || | ||___ | ||__| || | || | \\__/ || | || | |\ | || | || | || |___|| |________|| |__________|| |__|| \_________// F I N I T E E L E M E N T S F O R B I O M E C H A N I C S --- v e r s i o n - 2 . 4 . 3 . 8286 --- Musculoskeletal Research Laboratory University of Utah http://febio.org FEBio is a registered trademark. copyright (c) 2006-2015 - All rights reserved This is the NON-COMMERCIAL version of FEBio or the commercial license key file could not be found. If you have a key file make sure it is placed in the same directory as the executable. This version may only be used for non-commercial purposes as described in the license agreement. The functionality of this version may be limited compared to the commercial version. If you wish to obtain a valid commercial license file, please contact the developers. =========================================================================== tf_joint FILES USED =========================================================================== Input file : model_MOD.feb Plot file : model_MOD.xplt Log file : model_MOD.log MESH INFO =========================================================================== Number of materials ............................ : 4 Number of domains .............................. : 4 Number of nodes ................................ : 29877 Number of solid elements ....................... : 4096 Number of shell elements ....................... : 25220 Number of truss elements ....................... : 0 CONTROL DATA =========================================================================== Module type .................................... : solid Analysis type .................................. : quasi-static Plane strain mode .............................. : no Number of timesteps ............................ : 10 Time step size ................................. : 0.1 Auto time stepper activated .................... : yes Optimal nr of iterations ..................... : 10 Minimum allowable step size .................. : 0.001 Maximum allowable step size .................. : 0 Number of loadcurves ........................... : 7 Number of displacement boundary conditions ..... : 0 Number of concentrated nodal forces ............ : 0 Max nr of stiffness reformations ............... : 10 per time steps Max nr of Quasi-Newton iterations .............. : 0 between stiffness matrix reformations Displacement convergence tolerance ............. : 0.001 Energy convergence tolerance ................... : 0.01 Residual convergence tolerance ................. : 0 Minimal residual value ......................... : 0.001 Linesearch convergence tolerance ............... : 0.9 Minimum line search size ....................... : 0.01 Maximum number of line search iterations ....... : 5 Max condition number ........................... : 100000 OUTPUT DATA =========================================================================== plot level ................................ : must points only plotfile format ........................... : FEBIO MATERIAL DATA =========================================================================== 1 - femur (type: rigid body) density ........................................... : 1.132e-06 density ........................................... : 1.132e-06 E ................................................. : 1 v ................................................. : 0 parent_id ......................................... : -1 center_of_mass .................................... : 0,0,0 --------------------------------------------------------------------------- 2 - tibia (type: rigid body) density ........................................... : 1.132e-06 density ........................................... : 1.132e-06 E ................................................. : 1 v ................................................. : 0 parent_id ......................................... : -1 center_of_mass .................................... : 0,0,0 --------------------------------------------------------------------------- 3 - acl (type: trans iso Mooney-Rivlin) density ........................................... : 1.5e-09 k ................................................. : 73.2 laugon ............................................ : 0 atol .............................................. : 0.01 minaug ............................................ : 0 maxaug ............................................ : 0 c1 ................................................ : 1.95 c2 ................................................ : 0 c3 ................................................ : 0.0139 c4 ................................................ : 116.22 c5 ................................................ : 535.039 lam_max ........................................... : 1.046 --------------------------------------------------------------------------- 4 - aclfiber (type: trans iso Mooney-Rivlin) density ........................................... : 1.5e-09 k ................................................. : 73.2 laugon ............................................ : 0 atol .............................................. : 0.01 minaug ............................................ : 0 maxaug ............................................ : 0 c1 ................................................ : 1.95 c2 ................................................ : 0 c3 ................................................ : 0.0139 c4 ................................................ : 116.22 c5 ................................................ : 535.039 lam_max ........................................... : 1.046 RIGID BODY DATA =========================================================================== Rigid Body 1: material id : 1 center of mass : 0, 0, 0 --------------------------------------------------------------------------- Rigid Body 2: material id : 2 center of mass : 0, 0, 0 LOADCURVE DATA =========================================================================== 1 0 0 0.1 0.1 0.2 0.1 0.3 0.1 0.4 0.1 0.5 0.1 0.6 0.1 0.7 0.1 0.8 0.1 0.9 0.1 1 0.1 --------------------------------------------------------------------------- 2 0 0 0.1 0.1 0.9 0.9 1 1 --------------------------------------------------------------------------- 3 0 0 0.1 0.1 0.9 0.9 1 1 --------------------------------------------------------------------------- 4 0 0 0.1 0.1 0.9 0.9 1 1 --------------------------------------------------------------------------- 5 0 0 0.1 0.1 0.9 0.9 1 1 --------------------------------------------------------------------------- 6 0 0 0.1 0.1 0.9 0.9 1 1 --------------------------------------------------------------------------- 7 0 0 0.1 0.1 0.9 0.9 1 1 LINEAR SOLVER DATA =========================================================================== Solver type ............................... : Pardiso Matrix format ............................. : symmetric ===== reforming stiffness matrix: Nr of equations ........................... : 13272 Nr of nonzeroes in stiffness matrix ....... : 501420 ===== beginning time step 1 : 0.1 ===== 1 Nonlinear solution status: time= 0.1 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.119771e+07 3.550517e+03 0.000000e+00 energy 7.591579e+04 3.489976e+02 7.591579e+02 displacement 1.021312e+04 1.021312e+04 1.021312e-02 Reforming stiffness matrix: reformation #1 2 ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29053 at gauss point 3 * * jacobian = -0.078882 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29054 at gauss point 3 * * jacobian = -3.84286 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29055 at gauss point 1 * * jacobian = -45.8763 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29056 at gauss point 1 * * jacobian = -10.4284 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29057 at gauss point 2 * * jacobian = -2.98102 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29058 at gauss point 1 * * jacobian = -1.83668 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29059 at gauss point 1 * * jacobian = -2.49203 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29062 at gauss point 3 * * jacobian = -0.202932 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29063 at gauss point 1 * * jacobian = -1.70922 * * * ************************************************************************* ------- failed to converge at time : 0.1 Retrying time step. Retry attempt 1 of max 5 AUTO STEPPER: retry step, dt = 0.05 ===== beginning time step 1 : 0.05 ===== 1 Nonlinear solution status: time= 0.05 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.299428e+06 3.768513e+02 0.000000e+00 energy 1.897895e+04 1.185372e+02 1.897895e+02 displacement 2.553281e+03 2.553281e+03 2.553281e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.05 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.299428e+06 2.073480e+02 0.000000e+00 energy 1.897895e+04 2.539441e+00 1.897895e+02 displacement 2.553281e+03 1.139351e+01 2.450219e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.05 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.299428e+06 3.759802e+00 0.000000e+00 energy 1.897895e+04 6.360964e-02 1.897895e+02 displacement 2.553281e+03 1.230892e+00 2.459157e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.05 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.299428e+06 5.791063e-01 0.000000e+00 energy 1.897895e+04 6.510548e-04 1.897895e+02 displacement 2.553281e+03 6.377532e-02 2.458981e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.05 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.299428e+06 1.031658e-01 0.000000e+00 energy 1.897895e+04 5.187794e-03 1.897895e+02 displacement 2.553281e+03 5.382285e-03 2.459571e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.05 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.299428e+06 9.972869e-03 0.000000e+00 energy 1.897895e+04 2.139570e-03 1.897895e+02 displacement 2.553281e+03 3.820679e-03 2.460096e-03 Reforming stiffness matrix: reformation #6 7 Nonlinear solution status: time= 0.05 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 7 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.299428e+06 3.021032e-05 0.000000e+00 energy 1.897895e+04 8.341374e-05 1.897895e+02 displacement 2.553281e+03 1.177220e-03 2.460392e-03 convergence summary number of iterations : 7 number of reformations : 7 ------- converged at time : 0.05 Data Record #1 =========================================================================== Step = 1 Time = 0.05 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 1 Time = 0.05 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 1 Time = 0.05 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 1 Time = 0.05 Data = tibia motion File = tibia_kinematics.txt ===== beginning time step 2 : 0.1 ===== 1 Nonlinear solution status: time= 0.1 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.023477e+06 3.690449e+02 0.000000e+00 energy 1.795164e+04 1.149726e+02 1.795164e+02 displacement 2.392603e+03 2.392603e+03 2.392603e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.1 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.023477e+06 2.136778e+02 0.000000e+00 energy 1.795164e+04 2.604348e+00 1.795164e+02 displacement 2.392603e+03 7.617154e+00 2.315486e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.1 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.023477e+06 3.851055e+00 0.000000e+00 energy 1.795164e+04 6.440561e-02 1.795164e+02 displacement 2.392603e+03 1.420764e-01 2.327535e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.1 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.023477e+06 6.059343e-01 0.000000e+00 energy 1.795164e+04 1.211257e-03 1.795164e+02 displacement 2.392603e+03 6.561069e-02 2.327371e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.1 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.023477e+06 1.062106e-01 0.000000e+00 energy 1.795164e+04 5.322433e-03 1.795164e+02 displacement 2.392603e+03 5.529987e-03 2.327959e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.1 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.023477e+06 1.041255e-02 0.000000e+00 energy 1.795164e+04 2.229620e-03 1.795164e+02 displacement 2.392603e+03 3.968233e-03 2.328477e-03 Reforming stiffness matrix: reformation #6 7 Nonlinear solution status: time= 0.1 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 7 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.023477e+06 3.206105e-05 0.000000e+00 energy 1.795164e+04 8.956693e-05 1.795164e+02 displacement 2.392603e+03 1.245822e-03 2.328772e-03 convergence summary number of iterations : 7 number of reformations : 7 ------- converged at time : 0.1 Data Record #1 =========================================================================== Step = 2 Time = 0.1 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 2 Time = 0.1 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 2 Time = 0.1 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 2 Time = 0.1 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0597614 ===== beginning time step 3 : 0.159761 ===== 1 Nonlinear solution status: time= 0.159761 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.805453e+06 6.496639e+02 0.000000e+00 energy 2.428155e+04 1.688429e+02 2.428155e+02 displacement 3.247616e+03 3.247616e+03 3.247616e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.159761 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 2 step from line search = 0.467000 convergence norms : INITIAL CURRENT REQUIRED residual 6.805453e+06 1.071409e+02 0.000000e+00 energy 2.428155e+04 1.860504e+00 2.428155e+02 displacement 3.247616e+03 5.405519e+00 3.182727e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.159761 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.805453e+06 4.665696e+00 0.000000e+00 energy 2.428155e+04 6.741708e-02 2.428155e+02 displacement 3.247616e+03 6.086044e+00 3.145134e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.159761 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.805453e+06 9.909260e-02 0.000000e+00 energy 2.428155e+04 6.183422e-03 2.428155e+02 displacement 3.247616e+03 6.363738e-02 3.152966e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.159761 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.805453e+06 7.548589e-03 0.000000e+00 energy 2.428155e+04 1.845653e-03 2.428155e+02 displacement 3.247616e+03 1.375013e-02 3.153889e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.159761 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.805453e+06 3.559182e-05 0.000000e+00 energy 2.428155e+04 4.444179e-05 2.428155e+02 displacement 3.247616e+03 4.658949e-04 3.154042e-03 convergence summary number of iterations : 6 number of reformations : 6 ------- converged at time : 0.159761 Data Record #1 =========================================================================== Step = 3 Time = 0.15976143 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 3 Time = 0.15976143 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 3 Time = 0.15976143 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 3 Time = 0.15976143 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0678091 MUST POINT CONTROLLER: adjusting time step. dt = 0.0402386 ===== beginning time step 4 : 0.2 ===== 1 Nonlinear solution status: time= 0.2 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.901490e+06 1.747434e+02 0.000000e+00 energy 1.034139e+04 6.274594e+01 1.034139e+02 displacement 1.393822e+03 1.393822e+03 1.393822e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.2 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.901490e+06 5.547120e+00 0.000000e+00 energy 1.034139e+04 2.363870e-01 1.034139e+02 displacement 1.393822e+03 2.918227e+00 1.362557e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.2 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.901490e+06 1.071189e-01 0.000000e+00 energy 1.034139e+04 6.651018e-03 1.034139e+02 displacement 1.393822e+03 1.774648e-02 1.365452e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.2 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.901490e+06 2.550991e-03 0.000000e+00 energy 1.034139e+04 1.058003e-03 1.034139e+02 displacement 1.393822e+03 3.892584e-03 1.366076e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.2 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.901490e+06 1.542816e-06 0.000000e+00 energy 1.034139e+04 1.041768e-05 1.034139e+02 displacement 1.393822e+03 2.898158e-04 1.366214e-03 convergence summary number of iterations : 5 number of reformations : 5 ------- converged at time : 0.2 Data Record #1 =========================================================================== Step = 4 Time = 0.2 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 4 Time = 0.2 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 4 Time = 0.2 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 4 Time = 0.2 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0521909 ===== beginning time step 5 : 0.252191 ===== 1 Nonlinear solution status: time= 0.252191 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.691462e+06 4.169512e+02 0.000000e+00 energy 1.670163e+04 1.169247e+02 1.670163e+02 displacement 2.253209e+03 2.253209e+03 2.253209e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.252191 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.691462e+06 5.148950e+02 0.000000e+00 energy 1.670163e+04 5.086161e+00 1.670163e+02 displacement 2.253209e+03 7.703718e+00 2.176970e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.252191 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.691462e+06 2.129008e+01 0.000000e+00 energy 1.670163e+04 2.540562e-01 1.670163e+02 displacement 2.253209e+03 1.106891e-01 2.183286e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.252191 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.691462e+06 9.716640e-01 0.000000e+00 energy 1.670163e+04 2.801969e-02 1.670163e+02 displacement 2.253209e+03 4.498872e-02 2.185196e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.252191 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.691462e+06 2.030151e-01 0.000000e+00 energy 1.670163e+04 9.499418e-03 1.670163e+02 displacement 2.253209e+03 9.430929e-03 2.185877e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.252191 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.691462e+06 3.301416e-02 0.000000e+00 energy 1.670163e+04 3.800699e-03 1.670163e+02 displacement 2.253209e+03 5.201007e-03 2.186456e-03 Reforming stiffness matrix: reformation #6 7 Nonlinear solution status: time= 0.252191 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 7 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.691462e+06 7.763943e-04 0.000000e+00 energy 1.670163e+04 6.673799e-04 1.670163e+02 displacement 2.253209e+03 3.189908e-03 2.187011e-03 ************************************************************************* * WARNING * * * * No force acting on the system. * ************************************************************************* convergence summary number of iterations : 7 number of reformations : 7 ------- converged at time : 0.252191 Data Record #1 =========================================================================== Step = 5 Time = 0.252190856 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 5 Time = 0.252190856 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 5 Time = 0.252190856 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 5 Time = 0.252190856 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0615246 MUST POINT CONTROLLER: adjusting time step. dt = 0.0478091 ===== beginning time step 6 : 0.3 ===== 1 Nonlinear solution status: time= 0.3 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.752317e+06 3.121701e+02 0.000000e+00 energy 1.333294e+04 9.077781e+01 1.333294e+02 displacement 1.783329e+03 1.783329e+03 1.783329e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.3 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.752317e+06 1.013196e+02 0.000000e+00 energy 1.333294e+04 1.557876e+00 1.333294e+02 displacement 1.783329e+03 5.020086e+00 1.727201e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.3 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.752317e+06 2.138696e+00 0.000000e+00 energy 1.333294e+04 4.782958e-02 1.333294e+02 displacement 1.783329e+03 4.169828e-02 1.731217e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.3 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.752317e+06 1.945769e-01 0.000000e+00 energy 1.333294e+04 6.982627e-03 1.333294e+02 displacement 1.783329e+03 2.982879e-02 1.731702e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.3 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.752317e+06 4.855548e-02 0.000000e+00 energy 1.333294e+04 3.992061e-03 1.333294e+02 displacement 1.783329e+03 4.218284e-03 1.732183e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.3 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.752317e+06 1.316286e-03 0.000000e+00 energy 1.333294e+04 8.983833e-04 1.333294e+02 displacement 1.783329e+03 3.530081e-03 1.732737e-03 Reforming stiffness matrix: reformation #6 7 Nonlinear solution status: time= 0.3 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 7 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.752317e+06 1.238040e-06 0.000000e+00 energy 1.333294e+04 8.964665e-06 1.333294e+02 displacement 1.783329e+03 2.278609e-04 1.732891e-03 convergence summary number of iterations : 7 number of reformations : 7 ------- converged at time : 0.3 Data Record #1 =========================================================================== Step = 6 Time = 0.3 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 6 Time = 0.3 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 6 Time = 0.3 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 6 Time = 0.3 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0579983 ===== beginning time step 7 : 0.357998 ===== 1 Nonlinear solution status: time= 0.357998 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.288960e+06 5.789398e+02 0.000000e+00 energy 1.877909e+04 1.357945e+02 1.877909e+02 displacement 2.481532e+03 2.481532e+03 2.481532e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.357998 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 2 step from line search = 0.505675 convergence norms : INITIAL CURRENT REQUIRED residual 5.288960e+06 7.537524e+01 0.000000e+00 energy 1.877909e+04 1.467448e+00 1.877909e+02 displacement 2.481532e+03 3.006097e+00 2.428585e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.357998 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.288960e+06 3.231768e+00 0.000000e+00 energy 1.877909e+04 3.903268e-02 1.877909e+02 displacement 2.481532e+03 2.481893e+00 2.392144e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.357998 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.288960e+06 7.174120e-02 0.000000e+00 energy 1.877909e+04 5.367794e-03 1.877909e+02 displacement 2.481532e+03 5.043422e-03 2.393631e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.357998 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.288960e+06 4.216067e-03 0.000000e+00 energy 1.877909e+04 1.416719e-03 1.877909e+02 displacement 2.481532e+03 2.597461e-03 2.394196e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.357998 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 5.288960e+06 8.088394e-06 0.000000e+00 energy 1.877909e+04 2.107230e-05 1.877909e+02 displacement 2.481532e+03 3.543645e-04 2.394376e-03 convergence summary number of iterations : 6 number of reformations : 6 ------- converged at time : 0.357998 Data Record #1 =========================================================================== Step = 7 Time = 0.357998293 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 7 Time = 0.357998293 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 7 Time = 0.357998293 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 7 Time = 0.357998293 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0663986 MUST POINT CONTROLLER: adjusting time step. dt = 0.0420017 ===== beginning time step 8 : 0.4 ===== 1 Nonlinear solution status: time= 0.4 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.643939e+06 2.012753e+02 0.000000e+00 energy 9.377972e+03 5.935865e+01 9.377972e+01 displacement 1.214745e+03 1.214745e+03 1.214745e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.4 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.643939e+06 9.926716e+00 0.000000e+00 energy 9.377972e+03 3.249460e-01 9.377972e+01 displacement 1.214745e+03 3.360450e+00 1.181727e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.4 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.643939e+06 2.917654e-01 0.000000e+00 energy 9.377972e+03 1.050375e-02 9.377972e+01 displacement 1.214745e+03 1.001627e-02 1.183188e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.4 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.643939e+06 3.527442e-02 0.000000e+00 energy 9.377972e+03 3.381740e-03 9.377972e+01 displacement 1.214745e+03 4.475625e-03 1.183684e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.4 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.643939e+06 6.281272e-04 0.000000e+00 energy 9.377972e+03 5.318169e-04 9.377972e+01 displacement 1.214745e+03 2.696143e-03 1.184161e-03 ************************************************************************* * WARNING * * * * No force acting on the system. * ************************************************************************* convergence summary number of iterations : 5 number of reformations : 5 ------- converged at time : 0.4 Data Record #1 =========================================================================== Step = 8 Time = 0.4 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 8 Time = 0.4 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 8 Time = 0.4 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 8 Time = 0.4 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0536014 ===== beginning time step 9 : 0.453601 ===== 1 Nonlinear solution status: time= 0.453601 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.175353e+06 4.543686e+02 0.000000e+00 energy 1.479393e+04 1.027287e+02 1.479393e+02 displacement 1.884507e+03 1.884507e+03 1.884507e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.453601 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.175353e+06 6.656360e+02 0.000000e+00 energy 1.479393e+04 6.072591e+00 1.479393e+02 displacement 1.884507e+03 1.064416e+01 1.821236e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.453601 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.175353e+06 3.546862e+01 0.000000e+00 energy 1.479393e+04 3.468740e-01 1.479393e+02 displacement 1.884507e+03 9.415434e-02 1.825474e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.453601 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.175353e+06 1.250570e+00 0.000000e+00 energy 1.479393e+04 4.082930e-02 1.479393e+02 displacement 1.884507e+03 4.823918e-02 1.827296e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.453601 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.175353e+06 2.528630e-01 0.000000e+00 energy 1.479393e+04 1.115824e-02 1.479393e+02 displacement 1.884507e+03 1.152977e-02 1.827602e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.453601 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.175353e+06 4.282179e-02 0.000000e+00 energy 1.479393e+04 3.948311e-03 1.479393e+02 displacement 1.884507e+03 5.451960e-03 1.828203e-03 Reforming stiffness matrix: reformation #6 7 Nonlinear solution status: time= 0.453601 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 7 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.175353e+06 2.643772e-03 0.000000e+00 energy 1.479393e+04 1.108157e-03 1.479393e+02 displacement 1.884507e+03 3.050817e-03 1.828820e-03 Reforming stiffness matrix: reformation #7 8 Nonlinear solution status: time= 0.453601 stiffness updates = 0 right hand side evaluations = 9 stiffness matrix reformations = 8 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.175353e+06 1.116894e-06 0.000000e+00 energy 1.479393e+04 1.009997e-05 1.479393e+02 displacement 1.884507e+03 3.456382e-04 1.829032e-03 convergence summary number of iterations : 8 number of reformations : 8 ------- converged at time : 0.453601 Data Record #1 =========================================================================== Step = 9 Time = 0.453601366 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 9 Time = 0.453601366 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 9 Time = 0.453601366 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 9 Time = 0.453601366 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.059078 MUST POINT CONTROLLER: adjusting time step. dt = 0.0463986 ===== beginning time step 10 : 0.5 ===== 1 Nonlinear solution status: time= 0.5 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.014898e+06 2.840843e+02 0.000000e+00 energy 1.068297e+04 6.961720e+01 1.068297e+02 displacement 1.338239e+03 1.338239e+03 1.338239e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.5 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.014898e+06 5.299693e+01 0.000000e+00 energy 1.068297e+04 9.638893e-01 1.068297e+02 displacement 1.338239e+03 7.000468e+00 1.310765e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.5 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.014898e+06 1.029842e+00 0.000000e+00 energy 1.068297e+04 3.493649e-02 1.068297e+02 displacement 1.338239e+03 3.089738e-02 1.313375e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.5 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.014898e+06 1.094778e-01 0.000000e+00 energy 1.068297e+04 4.899482e-03 1.068297e+02 displacement 1.338239e+03 1.769348e-02 1.313884e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.5 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.014898e+06 1.604757e-02 0.000000e+00 energy 1.068297e+04 2.340909e-03 1.068297e+02 displacement 1.338239e+03 3.233347e-03 1.314452e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.5 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.014898e+06 9.493743e-05 0.000000e+00 energy 1.068297e+04 1.664867e-04 1.068297e+02 displacement 1.338239e+03 1.501841e-03 1.314899e-03 ************************************************************************* * WARNING * * * * No force acting on the system. * ************************************************************************* convergence summary number of iterations : 6 number of reformations : 6 ------- converged at time : 0.5 Data Record #1 =========================================================================== Step = 10 Time = 0.5 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 10 Time = 0.5 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 10 Time = 0.5 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 10 Time = 0.5 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0571189 ===== beginning time step 11 : 0.557119 ===== 1 Nonlinear solution status: time= 0.557119 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.448182e+06 5.603953e+02 0.000000e+00 energy 1.576201e+04 1.069718e+02 1.576201e+02 displacement 1.973428e+03 1.973428e+03 1.973428e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.557119 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 2 step from line search = 0.600764 convergence norms : INITIAL CURRENT REQUIRED residual 4.448182e+06 4.833017e+01 0.000000e+00 energy 1.576201e+04 1.064365e+00 1.576201e+02 displacement 1.973428e+03 6.696403e+00 1.957946e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.557119 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.448182e+06 2.631010e+00 0.000000e+00 energy 1.576201e+04 3.423056e-03 1.576201e+02 displacement 1.973428e+03 3.024656e+00 1.963839e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.557119 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.448182e+06 8.818881e-02 0.000000e+00 energy 1.576201e+04 4.869575e-03 1.576201e+02 displacement 1.973428e+03 3.326110e-02 1.966329e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.557119 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.448182e+06 8.570295e-03 0.000000e+00 energy 1.576201e+04 1.723996e-03 1.576201e+02 displacement 1.973428e+03 2.343952e-03 1.966996e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.557119 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.448182e+06 3.384153e-05 0.000000e+00 energy 1.576201e+04 7.190236e-05 1.576201e+02 displacement 1.973428e+03 8.063764e-04 1.967371e-03 convergence summary number of iterations : 6 number of reformations : 6 ------- converged at time : 0.557119 Data Record #1 =========================================================================== Step = 11 Time = 0.557118907 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 11 Time = 0.557118907 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 11 Time = 0.557118907 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 11 Time = 0.557118907 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0656951 MUST POINT CONTROLLER: adjusting time step. dt = 0.0428811 ===== beginning time step 12 : 0.6 ===== 1 Nonlinear solution status: time= 0.6 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.438109e+06 2.228650e+02 0.000000e+00 energy 8.648570e+03 5.324162e+01 8.648570e+01 displacement 1.122668e+03 1.122668e+03 1.122668e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.6 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.438109e+06 1.481887e+01 0.000000e+00 energy 8.648570e+03 4.021506e-01 8.648570e+01 displacement 1.122668e+03 5.794220e+00 1.136046e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.6 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.438109e+06 3.469646e-01 0.000000e+00 energy 8.648570e+03 1.500583e-02 8.648570e+01 displacement 1.122668e+03 8.207453e-02 1.139819e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.6 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.438109e+06 3.819748e-02 0.000000e+00 energy 8.648570e+03 3.137525e-03 8.648570e+01 displacement 1.122668e+03 5.788614e-03 1.140404e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.6 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.438109e+06 1.468911e-03 0.000000e+00 energy 8.648570e+03 6.897853e-04 8.648570e+01 displacement 1.122668e+03 2.224643e-03 1.140997e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.6 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.438109e+06 7.953227e-07 0.000000e+00 energy 8.648570e+03 5.897000e-06 8.648570e+01 displacement 1.122668e+03 2.022479e-04 1.141194e-03 convergence summary number of iterations : 6 number of reformations : 6 ------- converged at time : 0.6 Data Record #1 =========================================================================== Step = 12 Time = 0.6 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 12 Time = 0.6 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 12 Time = 0.6 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 12 Time = 0.6 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0543049 ===== beginning time step 13 : 0.654305 ===== 1 Nonlinear solution status: time= 0.654305 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.846918e+06 5.019918e+02 0.000000e+00 energy 1.366725e+04 9.004664e+01 1.366725e+02 displacement 1.874002e+03 1.874002e+03 1.874002e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.654305 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.846918e+06 4.842376e+02 0.000000e+00 energy 1.366725e+04 4.921591e+00 1.366725e+02 displacement 1.874002e+03 1.302175e+01 1.919993e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.654305 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.846918e+06 2.595396e+01 0.000000e+00 energy 1.366725e+04 2.937009e-01 1.366725e+02 displacement 1.874002e+03 4.070811e-01 1.933386e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.654305 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.846918e+06 1.120843e+00 0.000000e+00 energy 1.366725e+04 4.429017e-02 1.366725e+02 displacement 1.874002e+03 2.795522e-02 1.933880e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.654305 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.846918e+06 1.513553e-01 0.000000e+00 energy 1.366725e+04 9.065091e-03 1.366725e+02 displacement 1.874002e+03 1.316842e-02 1.934108e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.654305 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.846918e+06 3.071737e-02 0.000000e+00 energy 1.366725e+04 2.914475e-03 1.366725e+02 displacement 1.874002e+03 3.871500e-03 1.934817e-03 Reforming stiffness matrix: reformation #6 7 Nonlinear solution status: time= 0.654305 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 7 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.846918e+06 1.321323e-03 0.000000e+00 energy 1.366725e+04 6.346319e-04 1.366725e+02 displacement 1.874002e+03 2.007373e-03 1.935493e-03 Reforming stiffness matrix: reformation #7 8 Nonlinear solution status: time= 0.654305 stiffness updates = 0 right hand side evaluations = 9 stiffness matrix reformations = 8 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.846918e+06 4.579605e-07 0.000000e+00 energy 1.366725e+04 3.749853e-06 1.366725e+02 displacement 1.874002e+03 1.521255e-04 1.935679e-03 convergence summary number of iterations : 8 number of reformations : 8 ------- converged at time : 0.654305 Data Record #1 =========================================================================== Step = 13 Time = 0.654304874 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 13 Time = 0.654304874 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 13 Time = 0.654304874 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 13 Time = 0.654304874 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0596985 MUST POINT CONTROLLER: adjusting time step. dt = 0.0456951 ===== beginning time step 14 : 0.7 ===== 1 Nonlinear solution status: time= 0.7 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.684652e+06 2.973860e+02 0.000000e+00 energy 9.563600e+03 5.919346e+01 9.563600e+01 displacement 1.423169e+03 1.423169e+03 1.423169e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.7 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.684652e+06 3.952568e+01 0.000000e+00 energy 9.563600e+03 7.817107e-01 9.563600e+01 displacement 1.423169e+03 5.354417e+00 1.456907e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.7 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.684652e+06 7.486901e-01 0.000000e+00 energy 9.563600e+03 2.927725e-02 9.563600e+01 displacement 1.423169e+03 7.224221e-02 1.461453e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.7 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.684652e+06 8.664005e-02 0.000000e+00 energy 9.563600e+03 3.944251e-03 9.563600e+01 displacement 1.423169e+03 1.360824e-02 1.461792e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.7 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.684652e+06 1.015777e-02 0.000000e+00 energy 9.563600e+03 1.726278e-03 9.563600e+01 displacement 1.423169e+03 2.576080e-03 1.462451e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.7 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.684652e+06 4.870690e-05 0.000000e+00 energy 9.563600e+03 9.160301e-05 9.563600e+01 displacement 1.423169e+03 9.243083e-04 1.462924e-03 convergence summary number of iterations : 6 number of reformations : 6 ------- converged at time : 0.7 Data Record #1 =========================================================================== Step = 14 Time = 0.7 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 14 Time = 0.7 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 14 Time = 0.7 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 14 Time = 0.7 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0565561 ===== beginning time step 15 : 0.756556 ===== 1 Nonlinear solution status: time= 0.756556 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.084959e+06 6.207618e+02 0.000000e+00 energy 1.459323e+04 9.287355e+01 1.459323e+02 displacement 2.305371e+03 2.305371e+03 2.305371e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.756556 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.084959e+06 9.261170e+02 0.000000e+00 energy 1.459323e+04 7.895306e+00 1.459323e+02 displacement 2.305371e+03 1.154531e+01 2.368062e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.756556 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 3 step from line search = 0.296503 convergence norms : INITIAL CURRENT REQUIRED residual 4.084959e+06 3.269706e+02 0.000000e+00 energy 1.459323e+04 1.549111e-01 1.459323e+02 displacement 2.305371e+03 2.255762e-01 2.360298e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.756556 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 4 step from line search = 0.166715 convergence norms : INITIAL CURRENT REQUIRED residual 4.084959e+06 1.615391e+02 0.000000e+00 energy 1.459323e+04 1.314260e-01 1.459323e+02 displacement 2.305371e+03 1.004233e-01 2.360086e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.756556 stiffness updates = 0 right hand side evaluations = 9 stiffness matrix reformations = 5 step from line search = 0.247447 convergence norms : INITIAL CURRENT REQUIRED residual 4.084959e+06 3.335104e+01 0.000000e+00 energy 1.459323e+04 2.139162e-01 1.459323e+02 displacement 2.305371e+03 1.862433e-01 2.367679e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.756556 stiffness updates = 0 right hand side evaluations = 10 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.084959e+06 2.247177e+00 0.000000e+00 energy 1.459323e+04 7.494771e-02 1.459323e+02 displacement 2.305371e+03 9.653313e-02 2.370836e-03 Reforming stiffness matrix: reformation #6 7 Nonlinear solution status: time= 0.756556 stiffness updates = 0 right hand side evaluations = 11 stiffness matrix reformations = 7 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.084959e+06 2.830671e-01 0.000000e+00 energy 1.459323e+04 1.384950e-02 1.459323e+02 displacement 2.305371e+03 4.136034e-02 2.370609e-03 Reforming stiffness matrix: reformation #7 8 Nonlinear solution status: time= 0.756556 stiffness updates = 0 right hand side evaluations = 12 stiffness matrix reformations = 8 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.084959e+06 8.125053e-02 0.000000e+00 energy 1.459323e+04 4.986062e-03 1.459323e+02 displacement 2.305371e+03 5.519589e-03 2.371076e-03 Reforming stiffness matrix: reformation #8 9 Nonlinear solution status: time= 0.756556 stiffness updates = 0 right hand side evaluations = 13 stiffness matrix reformations = 9 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.084959e+06 1.290562e-02 0.000000e+00 energy 1.459323e+04 2.011687e-03 1.459323e+02 displacement 2.305371e+03 2.825995e-03 2.371832e-03 Reforming stiffness matrix: reformation #9 10 Nonlinear solution status: time= 0.756556 stiffness updates = 0 right hand side evaluations = 14 stiffness matrix reformations = 10 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.084959e+06 1.935130e-04 0.000000e+00 energy 1.459323e+04 1.961747e-04 1.459323e+02 displacement 2.305371e+03 1.187226e-03 2.372515e-03 convergence summary number of iterations : 10 number of reformations : 10 ------- converged at time : 0.756556 Data Record #1 =========================================================================== Step = 15 Time = 0.756556101 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 15 Time = 0.756556101 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 15 Time = 0.756556101 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 15 Time = 0.756556101 Data = tibia motion File = tibia_kinematics.txt MUST POINT CONTROLLER: adjusting time step. dt = 0.0434439 ===== beginning time step 16 : 0.8 ===== 1 Nonlinear solution status: time= 0.8 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.407163e+06 2.559297e+02 0.000000e+00 energy 8.636281e+03 4.958494e+01 8.636281e+01 displacement 1.443037e+03 1.443037e+03 1.443037e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.8 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.407163e+06 1.833501e+01 0.000000e+00 energy 8.636281e+03 4.605272e-01 8.636281e+01 displacement 1.443037e+03 3.681962e+00 1.472327e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.8 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.407163e+06 3.891145e-01 0.000000e+00 energy 8.636281e+03 1.753276e-02 8.636281e+01 displacement 1.443037e+03 2.513001e-02 1.471914e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.8 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.407163e+06 4.055544e-02 0.000000e+00 energy 8.636281e+03 3.143215e-03 8.636281e+01 displacement 1.443037e+03 7.027677e-03 1.472519e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.8 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.407163e+06 1.912663e-03 0.000000e+00 energy 8.636281e+03 7.502949e-04 8.636281e+01 displacement 1.443037e+03 2.109597e-03 1.473198e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.8 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.407163e+06 1.670369e-06 0.000000e+00 energy 8.636281e+03 9.030607e-06 8.636281e+01 displacement 1.443037e+03 2.395743e-04 1.473457e-03 convergence summary number of iterations : 6 number of reformations : 6 ------- converged at time : 0.8 Data Record #1 =========================================================================== Step = 16 Time = 0.8 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 16 Time = 0.8 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 16 Time = 0.8 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 16 Time = 0.8 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0547551 ===== beginning time step 17 : 0.854755 ===== 1 Nonlinear solution status: time= 0.854755 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.838636e+06 5.447577e+02 0.000000e+00 energy 1.382685e+04 8.115204e+01 1.382685e+02 displacement 2.394146e+03 2.394146e+03 2.394146e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.854755 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.838636e+06 5.311442e+02 0.000000e+00 energy 1.382685e+04 5.273446e+00 1.382685e+02 displacement 2.394146e+03 9.492820e+00 2.460539e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.854755 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.838636e+06 2.905148e+01 0.000000e+00 energy 1.382685e+04 3.177747e-01 1.382685e+02 displacement 2.394146e+03 1.198945e-01 2.459928e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.854755 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.838636e+06 1.232444e+00 0.000000e+00 energy 1.382685e+04 4.775826e-02 1.382685e+02 displacement 2.394146e+03 2.987792e-02 2.460634e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.854755 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.838636e+06 1.571212e-01 0.000000e+00 energy 1.382685e+04 9.691929e-03 1.382685e+02 displacement 2.394146e+03 1.430699e-02 2.460727e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.854755 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.838636e+06 3.212874e-02 0.000000e+00 energy 1.382685e+04 2.927905e-03 1.382685e+02 displacement 2.394146e+03 4.026158e-03 2.461478e-03 Reforming stiffness matrix: reformation #6 7 Nonlinear solution status: time= 0.854755 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 7 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.838636e+06 1.616996e-03 0.000000e+00 energy 1.382685e+04 6.855641e-04 1.382685e+02 displacement 2.394146e+03 1.961799e-03 2.462242e-03 convergence summary number of iterations : 7 number of reformations : 7 ------- converged at time : 0.854755 Data Record #1 =========================================================================== Step = 17 Time = 0.85475512 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 17 Time = 0.85475512 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 17 Time = 0.85475512 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 17 Time = 0.85475512 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0635882 MUST POINT CONTROLLER: adjusting time step. dt = 0.0452449 ===== beginning time step 18 : 0.9 ===== 1 Nonlinear solution status: time= 0.9 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.653769e+06 2.825496e+02 0.000000e+00 energy 9.606330e+03 5.158952e+01 9.606330e+01 displacement 1.734993e+03 1.734993e+03 1.734993e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.9 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.653769e+06 3.414148e+01 0.000000e+00 energy 9.606330e+03 7.052637e-01 9.606330e+01 displacement 1.734993e+03 4.292667e+00 1.779102e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.9 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.653769e+06 6.452721e-01 0.000000e+00 energy 9.606330e+03 2.669294e-02 9.606330e+01 displacement 1.734993e+03 2.760879e-02 1.779057e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.9 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.653769e+06 7.984094e-02 0.000000e+00 energy 9.606330e+03 3.601186e-03 9.606330e+01 displacement 1.734993e+03 1.180133e-02 1.779468e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.9 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.653769e+06 8.566548e-03 0.000000e+00 energy 9.606330e+03 1.565035e-03 9.606330e+01 displacement 1.734993e+03 2.404372e-03 1.780181e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.9 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.653769e+06 3.557538e-05 0.000000e+00 energy 9.606330e+03 7.182228e-05 9.606330e+01 displacement 1.734993e+03 7.832364e-04 1.780664e-03 convergence summary number of iterations : 6 number of reformations : 6 ------- converged at time : 0.9 Data Record #1 =========================================================================== Step = 18 Time = 0.9 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 18 Time = 0.9 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 18 Time = 0.9 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 18 Time = 0.9 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0561959 ===== beginning time step 19 : 0.956196 ===== 1 Nonlinear solution status: time= 0.956196 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.153036e+06 5.685667e+02 0.000000e+00 energy 1.510964e+04 7.945091e+01 1.510964e+02 displacement 2.824731e+03 2.824731e+03 2.824731e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.956196 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.153036e+06 8.036233e+02 0.000000e+00 energy 1.510964e+04 7.133843e+00 1.510964e+02 displacement 2.824731e+03 1.043864e+01 2.922516e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.956196 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.153036e+06 4.431147e+01 0.000000e+00 energy 1.510964e+04 4.042744e-01 1.510964e+02 displacement 2.824731e+03 1.388590e-01 2.919964e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.956196 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.153036e+06 1.644661e+00 0.000000e+00 energy 1.510964e+04 4.721368e-02 1.510964e+02 displacement 2.824731e+03 6.388664e-02 2.923017e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.956196 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.153036e+06 2.378687e-01 0.000000e+00 energy 1.510964e+04 1.384120e-02 1.510964e+02 displacement 2.824731e+03 1.498807e-02 2.922265e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.956196 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.153036e+06 4.814170e-02 0.000000e+00 energy 1.510964e+04 3.404602e-03 1.510964e+02 displacement 2.824731e+03 5.569485e-03 2.923063e-03 Reforming stiffness matrix: reformation #6 7 Nonlinear solution status: time= 0.956196 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 7 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.153036e+06 5.321018e-03 0.000000e+00 energy 1.510964e+04 1.236445e-03 1.510964e+02 displacement 2.824731e+03 2.201354e-03 2.923911e-03 convergence summary number of iterations : 7 number of reformations : 7 ------- converged at time : 0.956196 Data Record #1 =========================================================================== Step = 19 Time = 0.956195904 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 19 Time = 0.956195904 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 19 Time = 0.956195904 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 19 Time = 0.956195904 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0647477 MUST POINT CONTROLLER: adjusting time step. dt = 0.0438041 ===== beginning time step 20 : 1 ===== 1 Nonlinear solution status: time= 1 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.590641e+06 2.477242e+02 0.000000e+00 energy 9.477557e+03 4.558822e+01 9.477557e+01 displacement 1.844426e+03 1.844426e+03 1.844426e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 1 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.590641e+06 2.079315e+01 0.000000e+00 energy 9.477557e+03 5.007699e-01 9.477557e+01 displacement 1.844426e+03 3.696348e+00 1.896107e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 1 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.590641e+06 4.251899e-01 0.000000e+00 energy 9.477557e+03 1.923814e-02 9.477557e+01 displacement 1.844426e+03 1.697393e-02 1.896225e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 1 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.590641e+06 4.314569e-02 0.000000e+00 energy 9.477557e+03 3.229999e-03 9.477557e+01 displacement 1.844426e+03 7.888579e-03 1.896831e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 1 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.590641e+06 2.274406e-03 0.000000e+00 energy 9.477557e+03 8.087324e-04 9.477557e+01 displacement 1.844426e+03 2.118018e-03 1.897572e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 1 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 6 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.590641e+06 2.527422e-06 0.000000e+00 energy 9.477557e+03 1.171909e-05 9.477557e+01 displacement 1.844426e+03 2.722329e-04 1.897880e-03 convergence summary number of iterations : 6 number of reformations : 6 ------- converged at time : 1 Data Record #1 =========================================================================== Step = 20 Time = 1 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 20 Time = 1 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 20 Time = 1 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 20 Time = 1 Data = tibia motion File = tibia_kinematics.txt N O N L I N E A R I T E R A T I O N I N F O R M A T I O N Number of time steps completed .................... : 20 Total number of equilibrium iterations ............ : 133 Average number of equilibrium iterations .......... : 6.65 Total number of right hand evaluations ............ : 160 Total number of stiffness reformations ............ : 134 Time in solver: 0:01:17 Elapsed time : 0:01:38 N O R M A L T E R M I N A T I O N