=========================================================================== ________ _________ _________ __ _________ | |\ | |\ | |\ | |\ / \\ | ____|| | ____|| | __ || |__|| | ___ || | |\___\| | |\___\| | |\_| || \_\| | // \ || | || | || | || | || __ | || | || | ||__ | ||__ | ||_| || | |\ | || | || | |\ | |\ | \\ | || | || | || | ___|| | ___|| | ___ || | || | || | || | |\__\| | |\__\| | |\__| || | || | || | || | || | || | || | || | || | || | || | || | ||___ | ||__| || | || | \\__/ || | || | |\ | || | || | || |___|| |________|| |__________|| |__|| \_________// 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 . 5 . 0 . 8514 --- 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. This version may only be used for non-commercial purposes as described in the license agreement. 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 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 ............... : 1090021888 per time steps Max nr of Quasi-Newton iterations .............. : 0 between stiffness matrix reformations Linesearch convergence tolerance ............... : 0.9 Minimum line search size ....................... : 0.01 Maximum number of line search iterations ....... : 5 Max condition number ........................... : 0.01 Displacement convergence tolerance ............. : 0.01 Energy convergence tolerance ................... : 0.1 Residual convergence tolerance ................. : 0 Minimal residual value ......................... : 0.001 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 ===== beginning time step 1 : 0.1 ===== ===== reforming stiffness matrix: Nr of equations ........................... : 13272 Nr of nonzeroes in stiffness matrix ....... : 501420 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+03 displacement 1.021312e+04 1.021312e+04 1.021312e+00 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+03 displacement 2.553281e+03 2.553281e+03 2.553281e-01 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+03 displacement 2.553281e+03 1.139351e+01 2.450219e-01 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+03 displacement 2.553281e+03 1.230892e+00 2.459157e-01 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+03 displacement 2.553281e+03 6.377532e-02 2.458981e-01 convergence summary number of iterations : 4 number of reformations : 4 ------- 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.261481e+06 3.639929e+02 0.000000e+00 energy 1.841734e+04 1.129830e+02 1.841734e+03 displacement 2.393129e+03 2.393129e+03 2.393129e-01 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.261481e+06 3.240456e+01 0.000000e+00 energy 1.841734e+04 7.106155e-01 1.841734e+03 displacement 2.393129e+03 7.518111e+00 2.318102e-01 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.261481e+06 6.686045e-01 0.000000e+00 energy 1.841734e+04 2.822264e-02 1.841734e+03 displacement 2.393129e+03 8.942496e-02 2.328690e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- 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.06 ===== beginning time step 3 : 0.16 ===== 1 Nonlinear solution status: time= 0.16 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 7.108926e+06 6.526147e+02 0.000000e+00 energy 2.497489e+04 1.684988e+02 2.497489e+03 displacement 3.274611e+03 3.274611e+03 3.274611e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.16 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 7.108926e+06 4.410430e+02 0.000000e+00 energy 2.497489e+04 4.958060e+00 2.497489e+03 displacement 3.274611e+03 2.125944e+01 3.148924e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.16 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 7.108926e+06 2.225112e+01 0.000000e+00 energy 2.497489e+04 3.055150e-01 2.497489e+03 displacement 3.274611e+03 1.030563e+00 3.171729e-01 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.16 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 7.108926e+06 1.435031e+00 0.000000e+00 energy 2.497489e+04 3.304153e-02 2.497489e+03 displacement 3.274611e+03 1.749898e+00 3.177988e-01 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.16 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 7.108926e+06 7.650535e-02 0.000000e+00 energy 2.497489e+04 6.367108e-03 2.497489e+03 displacement 3.274611e+03 6.485824e-02 3.179256e-01 convergence summary number of iterations : 5 number of reformations : 5 ------- converged at time : 0.16 Data Record #1 =========================================================================== Step = 3 Time = 0.16 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 3 Time = 0.16 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 3 Time = 0.16 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 3 Time = 0.16 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.068 MUST POINT CONTROLLER: adjusting time step. dt = 0.04 ===== 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.906001e+06 1.726337e+02 0.000000e+00 energy 1.029281e+04 6.185789e+01 1.029281e+03 displacement 1.378032e+03 1.378032e+03 1.378032e-01 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.906001e+06 4.234413e+00 0.000000e+00 energy 1.029281e+04 1.952688e-01 1.029281e+03 displacement 1.378032e+03 2.848369e+00 1.347580e-01 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.906001e+06 7.056323e-02 0.000000e+00 energy 1.029281e+04 5.374527e-03 1.029281e+03 displacement 1.378032e+03 1.891704e-02 1.350449e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- 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.0744 ===== beginning time step 5 : 0.2744 ===== 1 Nonlinear solution status: time= 0.2744 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 9.591787e+06 1.200182e+03 0.000000e+00 energy 3.404352e+04 2.397484e+02 3.404352e+03 displacement 4.580263e+03 4.580263e+03 4.580263e-01 Reforming stiffness matrix: reformation #1 2 ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29056 at gauss point 1 * * jacobian = -0.648697 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29057 at gauss point 1 * * jacobian = -1.10737 * * * ************************************************************************* ------- failed to converge at time : 0.2744 Retrying time step. Retry attempt 1 of max 5 AUTO STEPPER: retry step, dt = 0.0372 ===== beginning time step 5 : 0.2372 ===== 1 Nonlinear solution status: time= 0.2372 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.398104e+06 1.331289e+02 0.000000e+00 energy 8.511500e+03 4.987783e+01 8.511500e+02 displacement 1.145354e+03 1.145354e+03 1.145354e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.2372 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.398104e+06 2.226929e+00 0.000000e+00 energy 8.511500e+03 1.185463e-01 8.511500e+02 displacement 1.145354e+03 1.876733e+00 1.119899e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.2372 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.398104e+06 1.459424e-02 0.000000e+00 energy 8.511500e+03 2.309951e-03 8.511500e+02 displacement 1.145354e+03 1.147117e-02 1.121583e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.2372 Data Record #1 =========================================================================== Step = 5 Time = 0.2372 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 5 Time = 0.2372 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 5 Time = 0.2372 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 5 Time = 0.2372 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.04976 ===== beginning time step 6 : 0.28696 ===== 1 Nonlinear solution status: time= 0.28696 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.125543e+06 3.562210e+02 0.000000e+00 energy 1.465761e+04 1.013963e+02 1.465761e+03 displacement 1.965827e+03 1.965827e+03 1.965827e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.28696 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.125543e+06 2.120703e+02 0.000000e+00 energy 1.465761e+04 2.645502e+00 1.465761e+03 displacement 1.965827e+03 6.034531e+00 1.901030e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.28696 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.125543e+06 3.908673e+00 0.000000e+00 energy 1.465761e+04 6.677981e-02 1.465761e+03 displacement 1.965827e+03 8.370962e-02 1.906514e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.28696 Data Record #1 =========================================================================== Step = 6 Time = 0.28696 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 6 Time = 0.28696 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 6 Time = 0.28696 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 6 Time = 0.28696 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.059808 MUST POINT CONTROLLER: adjusting time step. dt = 0.01304 ===== beginning time step 7 : 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 2.997094e+05 4.433607e+00 0.000000e+00 energy 1.016888e+03 2.395879e+00 1.016888e+02 displacement 1.275576e+02 1.275576e+02 1.275576e-02 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 2.997094e+05 1.146654e-01 0.000000e+00 energy 1.016888e+03 4.762937e-03 1.016888e+02 displacement 1.275576e+02 2.954349e-02 1.267181e-02 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 2.997094e+05 9.237181e-03 0.000000e+00 energy 1.016888e+03 2.325967e-03 1.016888e+02 displacement 1.275576e+02 4.294581e-03 1.269206e-02 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.3 Data Record #1 =========================================================================== Step = 7 Time = 0.3 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 7 Time = 0.3 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 7 Time = 0.3 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 7 Time = 0.3 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0678464 ===== beginning time step 8 : 0.367846 ===== 1 Nonlinear solution status: time= 0.367846 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 7.255761e+06 9.171381e+02 0.000000e+00 energy 2.572868e+04 1.851246e+02 2.572868e+03 displacement 3.396382e+03 3.396382e+03 3.396382e-01 Reforming stiffness matrix: reformation #1 2 ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29056 at gauss point 3 * * jacobian = -0.0319881 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29057 at gauss point 4 * * jacobian = -0.000703481 * * * ************************************************************************* ------- failed to converge at time : 0.367846 Retrying time step. Retry attempt 1 of max 5 AUTO STEPPER: retry step, dt = 0.0339232 ===== beginning time step 8 : 0.333923 ===== 1 Nonlinear solution status: time= 0.333923 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.813976e+06 9.570188e+01 0.000000e+00 energy 6.432353e+03 3.592970e+01 6.432353e+02 displacement 8.492271e+02 8.492271e+02 8.492271e-02 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.333923 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.813976e+06 9.342703e-01 0.000000e+00 energy 6.432353e+03 6.123152e-02 6.432353e+02 displacement 8.492271e+02 1.245479e+00 8.310576e-02 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.333923 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.813976e+06 1.922400e-03 0.000000e+00 energy 6.432353e+03 5.504632e-04 6.432353e+02 displacement 8.492271e+02 6.357142e-03 8.319277e-02 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.333923 Data Record #1 =========================================================================== Step = 8 Time = 0.3339232 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 8 Time = 0.3339232 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 8 Time = 0.3339232 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 8 Time = 0.3339232 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0471386 ===== beginning time step 9 : 0.381062 ===== 1 Nonlinear solution status: time= 0.381062 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.396321e+06 2.971680e+02 0.000000e+00 energy 1.204947e+04 8.107007e+01 1.204947e+03 displacement 1.574576e+03 1.574576e+03 1.574576e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.381062 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.396321e+06 7.447911e+01 0.000000e+00 energy 1.204947e+04 1.243948e+00 1.204947e+03 displacement 1.574576e+03 5.128707e+00 1.524968e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.381062 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.396321e+06 1.543960e+00 0.000000e+00 energy 1.204947e+04 4.101060e-02 1.204947e+03 displacement 1.574576e+03 2.807995e-02 1.528031e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.381062 Data Record #1 =========================================================================== Step = 9 Time = 0.38106176 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 9 Time = 0.38106176 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 9 Time = 0.38106176 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 9 Time = 0.38106176 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0577108 MUST POINT CONTROLLER: adjusting time step. dt = 0.0189382 ===== beginning time step 10 : 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 5.686684e+05 1.141018e+01 0.000000e+00 energy 1.953093e+03 6.400822e+00 1.953093e+02 displacement 2.406425e+02 2.406425e+02 2.406425e-02 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 5.686684e+05 6.081368e-02 0.000000e+00 energy 1.953093e+03 3.042543e-03 1.953093e+02 displacement 2.406425e+02 1.450969e-01 2.381219e-02 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 5.686684e+05 1.511516e-03 0.000000e+00 energy 1.953093e+03 9.541596e-04 1.953093e+02 displacement 2.406425e+02 3.467680e-03 2.384022e-02 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.4 Data Record #1 =========================================================================== Step = 10 Time = 0.4 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 10 Time = 0.4 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 10 Time = 0.4 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 10 Time = 0.4 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0661687 ===== beginning time step 11 : 0.466169 ===== 1 Nonlinear solution status: time= 0.466169 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.366207e+06 8.562153e+02 0.000000e+00 energy 2.254967e+04 1.565554e+02 2.254967e+03 displacement 2.871871e+03 2.871871e+03 2.871871e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.466169 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 2 step from line search = 0.367277 convergence norms : INITIAL CURRENT REQUIRED residual 6.366207e+06 2.191896e+02 0.000000e+00 energy 2.254967e+04 2.266601e+00 2.254967e+03 displacement 2.871871e+03 3.567759e+00 2.821457e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.466169 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.366207e+06 1.433473e+01 0.000000e+00 energy 2.254967e+04 2.282479e-01 2.254967e+03 displacement 2.871871e+03 9.117900e+00 2.762996e-01 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.466169 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 6.366207e+06 4.020234e-01 0.000000e+00 energy 2.254967e+04 1.212877e-02 2.254967e+03 displacement 2.871871e+03 4.116457e-02 2.767484e-01 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 0.466169 Data Record #1 =========================================================================== Step = 11 Time = 0.466168678 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 11 Time = 0.466168678 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 11 Time = 0.466168678 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 11 Time = 0.466168678 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0729349 MUST POINT CONTROLLER: adjusting time step. dt = 0.0338313 ===== beginning time step 12 : 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 1.614096e+06 9.410821e+01 0.000000e+00 energy 5.682197e+03 3.059462e+01 5.682197e+02 displacement 7.048905e+02 7.048905e+02 7.048905e-02 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 1.614096e+06 8.191089e-01 0.000000e+00 energy 5.682197e+03 4.767499e-02 5.682197e+02 displacement 7.048905e+02 2.032830e+00 6.964400e-02 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 1.614096e+06 4.209266e-03 0.000000e+00 energy 5.682197e+03 1.084838e-03 5.682197e+02 displacement 7.048905e+02 5.504761e-03 6.972743e-02 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.5 Data Record #1 =========================================================================== Step = 12 Time = 0.5 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 12 Time = 0.5 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 12 Time = 0.5 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 12 Time = 0.5 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.078348 ===== beginning time step 13 : 0.578348 ===== 1 Nonlinear solution status: time= 0.578348 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 8.377337e+06 1.445590e+03 0.000000e+00 energy 2.966943e+04 1.729998e+02 2.966943e+03 displacement 3.713182e+03 3.713182e+03 3.713182e-01 Reforming stiffness matrix: reformation #1 2 ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29055 at gauss point 2 * * jacobian = -0.110107 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29056 at gauss point 1 * * jacobian = -3.3718 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29057 at gauss point 1 * * jacobian = -5.00833 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29058 at gauss point 1 * * jacobian = -1.00803 * * * ************************************************************************* ------- failed to converge at time : 0.578348 Retrying time step. Retry attempt 1 of max 5 AUTO STEPPER: retry step, dt = 0.039174 ===== beginning time step 13 : 0.539174 ===== 1 Nonlinear solution status: time= 0.539174 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.094369e+06 1.604248e+02 0.000000e+00 energy 7.417499e+03 4.402103e+01 7.417499e+02 displacement 9.283791e+02 9.283791e+02 9.283791e-02 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.539174 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.094369e+06 3.839205e+00 0.000000e+00 energy 7.417499e+03 1.705852e-01 7.417499e+02 displacement 9.283791e+02 4.018155e+00 9.230658e-02 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.539174 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.094369e+06 7.043857e-02 0.000000e+00 energy 7.417499e+03 4.847768e-03 7.417499e+02 displacement 9.283791e+02 2.130917e-02 9.246698e-02 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.539174 Data Record #1 =========================================================================== Step = 13 Time = 0.539173977 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 13 Time = 0.539173977 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 13 Time = 0.539173977 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 13 Time = 0.539173977 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0513392 ===== beginning time step 14 : 0.590513 ===== 1 Nonlinear solution status: time= 0.590513 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.543041e+06 4.055618e+02 0.000000e+00 energy 1.252992e+04 8.276860e+01 1.252992e+03 displacement 1.595855e+03 1.595855e+03 1.595855e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.590513 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.543041e+06 1.903279e+02 0.000000e+00 energy 1.252992e+04 2.440705e+00 1.252992e+03 displacement 1.595855e+03 1.216477e+01 1.607436e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.590513 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.543041e+06 3.049130e+00 0.000000e+00 energy 1.252992e+04 5.067185e-02 1.252992e+03 displacement 1.595855e+03 2.559070e-01 1.615977e-01 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.590513 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.543041e+06 6.878980e-01 0.000000e+00 energy 1.252992e+04 4.310695e-03 1.252992e+03 displacement 1.595855e+03 7.294261e-02 1.615959e-01 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 0.590513 Data Record #1 =========================================================================== Step = 14 Time = 0.590513159 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 14 Time = 0.590513159 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 14 Time = 0.590513159 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 14 Time = 0.590513159 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0610713 MUST POINT CONTROLLER: adjusting time step. dt = 0.00948684 ===== beginning time step 15 : 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 1.230285e+05 6.168157e-01 0.000000e+00 energy 4.296571e+02 5.994670e-01 4.296571e+01 displacement 5.663300e+01 5.663300e+01 5.663300e-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 1.230285e+05 8.397062e-03 0.000000e+00 energy 4.296571e+02 1.912692e-03 4.296571e+01 displacement 5.663300e+01 1.893198e-02 5.706581e-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 1.230285e+05 4.450914e-05 0.000000e+00 energy 4.296571e+02 9.766839e-05 4.296571e+01 displacement 5.663300e+01 1.058499e-03 5.716581e-03 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.6 Data Record #1 =========================================================================== Step = 15 Time = 0.6 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 15 Time = 0.6 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 15 Time = 0.6 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 15 Time = 0.6 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0688571 ===== beginning time step 16 : 0.668857 ===== 1 Nonlinear solution status: time= 0.668857 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.185802e+06 1.036785e+03 0.000000e+00 energy 2.197492e+04 1.374424e+02 2.197492e+03 displacement 3.012979e+03 3.012979e+03 3.012979e-01 Reforming stiffness matrix: reformation #1 2 ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29056 at gauss point 7 * * jacobian = -0.0719052 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29057 at gauss point 8 * * jacobian = -0.035071 * * * ************************************************************************* ------- failed to converge at time : 0.668857 Retrying time step. Retry attempt 1 of max 5 AUTO STEPPER: retry step, dt = 0.0344285 ===== beginning time step 16 : 0.634429 ===== 1 Nonlinear solution status: time= 0.634429 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.546452e+06 1.055730e+02 0.000000e+00 energy 5.493742e+03 2.959747e+01 5.493742e+02 displacement 7.532586e+02 7.532586e+02 7.532586e-02 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.634429 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.546452e+06 1.012452e+00 0.000000e+00 energy 5.493742e+03 6.550601e-02 5.493742e+02 displacement 7.532586e+02 2.142980e+00 7.663974e-02 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.634429 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.546452e+06 4.084625e-03 0.000000e+00 energy 5.493742e+03 8.871366e-04 5.493742e+02 displacement 7.532586e+02 3.491253e-02 7.687857e-02 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.634429 Data Record #1 =========================================================================== Step = 16 Time = 0.634428538 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 16 Time = 0.634428538 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 16 Time = 0.634428538 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 16 Time = 0.634428538 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0475428 ===== beginning time step 17 : 0.681971 ===== 1 Nonlinear solution status: time= 0.681971 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.921379e+06 3.338630e+02 0.000000e+00 energy 1.039193e+04 6.569843e+01 1.039193e+03 displacement 1.500440e+03 1.500440e+03 1.500440e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.681971 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.921379e+06 7.265550e+01 0.000000e+00 energy 1.039193e+04 1.204782e+00 1.039193e+03 displacement 1.500440e+03 6.748813e+00 1.537396e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.681971 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.921379e+06 1.523443e+00 0.000000e+00 energy 1.039193e+04 4.107727e-02 1.039193e+03 displacement 1.500440e+03 1.546260e-01 1.544740e-01 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.681971 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.921379e+06 1.258777e-01 0.000000e+00 energy 1.039193e+04 5.810525e-03 1.039193e+03 displacement 1.500440e+03 2.372750e-02 1.544901e-01 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 0.681971 Data Record #1 =========================================================================== Step = 17 Time = 0.681971369 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 17 Time = 0.681971369 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 17 Time = 0.681971369 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 17 Time = 0.681971369 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0580343 MUST POINT CONTROLLER: adjusting time step. dt = 0.0180286 ===== beginning time step 18 : 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 4.239993e+05 9.741104e+00 0.000000e+00 energy 1.502292e+03 4.618020e+00 1.502292e+02 displacement 2.296199e+02 2.296199e+02 2.296199e-02 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 4.239993e+05 3.832845e-03 0.000000e+00 energy 1.502292e+03 7.366980e-05 1.502292e+02 displacement 2.296199e+02 1.252724e-01 2.321724e-02 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 4.239993e+05 5.597111e-06 0.000000e+00 energy 1.502292e+03 6.431097e-06 1.502292e+02 displacement 2.296199e+02 5.142259e-04 2.323495e-02 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.7 Data Record #1 =========================================================================== Step = 18 Time = 0.7 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 18 Time = 0.7 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 18 Time = 0.7 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 18 Time = 0.7 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0664274 ===== beginning time step 19 : 0.766427 ===== 1 Nonlinear solution status: time= 0.766427 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.634636e+06 1.024248e+03 0.000000e+00 energy 2.013077e+04 1.230274e+02 2.013077e+03 displacement 3.180278e+03 3.180278e+03 3.180278e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.766427 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 2 step from line search = 0.206477 convergence norms : INITIAL CURRENT REQUIRED residual 5.634636e+06 5.816755e+02 0.000000e+00 energy 2.013077e+04 2.950100e+00 2.013077e+03 displacement 3.180278e+03 1.000118e+00 3.198895e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.766427 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.634636e+06 4.183455e+01 0.000000e+00 energy 2.013077e+04 9.026593e-01 2.013077e+03 displacement 3.180278e+03 1.248739e+01 3.278159e-01 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.766427 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.634636e+06 6.198951e-01 0.000000e+00 energy 2.013077e+04 2.793196e-02 2.013077e+03 displacement 3.180278e+03 1.325067e-01 3.285832e-01 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 0.766427 Data Record #1 =========================================================================== Step = 19 Time = 0.766427412 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 19 Time = 0.766427412 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 19 Time = 0.766427412 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 19 Time = 0.766427412 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0731419 MUST POINT CONTROLLER: adjusting time step. dt = 0.0335726 ===== beginning time step 20 : 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 1.473495e+06 9.937677e+01 0.000000e+00 energy 5.246643e+03 2.510285e+01 5.246643e+02 displacement 8.708827e+02 8.708827e+02 8.708827e-02 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 1.473495e+06 7.993681e-01 0.000000e+00 energy 5.246643e+03 4.466596e-02 5.246643e+02 displacement 8.708827e+02 1.291587e+00 8.851832e-02 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 1.473495e+06 7.224873e-03 0.000000e+00 energy 5.246643e+03 1.259014e-03 5.246643e+02 displacement 8.708827e+02 5.872987e-03 8.853672e-02 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.8 Data Record #1 =========================================================================== Step = 20 Time = 0.8 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 20 Time = 0.8 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 20 Time = 0.8 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 20 Time = 0.8 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0785135 ===== beginning time step 21 : 0.878514 ===== 1 Nonlinear solution status: time= 0.878514 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 7.905111e+06 1.672318e+03 0.000000e+00 energy 2.845365e+04 1.280297e+02 2.845365e+03 displacement 4.923286e+03 4.923286e+03 4.923286e-01 Reforming stiffness matrix: reformation #1 2 ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29055 at gauss point 2 * * jacobian = -0.259046 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29056 at gauss point 1 * * jacobian = -4.07704 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29057 at gauss point 1 * * jacobian = -5.95805 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29058 at gauss point 1 * * jacobian = -1.2049 * * * ************************************************************************* ------- failed to converge at time : 0.878514 Retrying time step. Retry attempt 1 of max 5 AUTO STEPPER: retry step, dt = 0.0392568 ===== beginning time step 21 : 0.839257 ===== 1 Nonlinear solution status: time= 0.839257 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.976318e+06 1.765640e+02 0.000000e+00 energy 7.113601e+03 3.761512e+01 7.113601e+02 displacement 1.230993e+03 1.230993e+03 1.230993e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.839257 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.976318e+06 3.978819e+00 0.000000e+00 energy 7.113601e+03 1.749119e-01 7.113601e+02 displacement 1.230993e+03 2.419489e+00 1.255347e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.839257 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.976318e+06 8.382761e-02 0.000000e+00 energy 7.113601e+03 4.928070e-03 7.113601e+02 displacement 1.230993e+03 1.690508e-02 1.255226e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.839257 Data Record #1 =========================================================================== Step = 21 Time = 0.839256772 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 21 Time = 0.839256772 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 21 Time = 0.839256772 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 21 Time = 0.839256772 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0514054 ===== beginning time step 22 : 0.890662 ===== 1 Nonlinear solution status: time= 0.890662 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.434913e+06 4.359065e+02 0.000000e+00 energy 1.238266e+04 6.916461e+01 1.238266e+03 displacement 2.201576e+03 2.201576e+03 2.201576e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.890662 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.434913e+06 1.349373e+02 0.000000e+00 energy 1.238266e+04 1.930833e+00 1.238266e+03 displacement 2.201576e+03 7.239666e+00 2.263358e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.890662 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.434913e+06 2.476870e+00 0.000000e+00 energy 1.238266e+04 5.320943e-02 1.238266e+03 displacement 2.201576e+03 8.516003e-02 2.264264e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.890662 Data Record #1 =========================================================================== Step = 22 Time = 0.890662189 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 22 Time = 0.890662189 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 22 Time = 0.890662189 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 22 Time = 0.890662189 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0611243 MUST POINT CONTROLLER: adjusting time step. dt = 0.00933781 ===== beginning time step 23 : 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 1.251327e+05 1.398814e+00 0.000000e+00 energy 4.381698e+02 5.347256e-01 4.381698e+01 displacement 7.718095e+01 7.718095e+01 7.718095e-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 1.251327e+05 6.778547e-02 0.000000e+00 energy 4.381698e+02 3.737527e-03 4.381698e+01 displacement 7.718095e+01 1.076380e-02 7.767658e-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 1.251327e+05 5.304041e-03 0.000000e+00 energy 4.381698e+02 1.301626e-03 4.381698e+01 displacement 7.718095e+01 2.522744e-03 7.785498e-03 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.9 Data Record #1 =========================================================================== Step = 23 Time = 0.9 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 23 Time = 0.9 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 23 Time = 0.9 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 23 Time = 0.9 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0688995 ===== beginning time step 24 : 0.968899 ===== 1 Nonlinear solution status: time= 0.968899 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.251563e+06 1.059525e+03 0.000000e+00 energy 2.272914e+04 1.056259e+02 2.272914e+03 displacement 4.246732e+03 4.246732e+03 4.246732e-01 Reforming stiffness matrix: reformation #1 2 ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29056 at gauss point 7 * * jacobian = -0.0862162 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29057 at gauss point 8 * * jacobian = -0.0525339 * * * ************************************************************************* ------- failed to converge at time : 0.968899 Retrying time step. Retry attempt 1 of max 5 AUTO STEPPER: retry step, dt = 0.0344497 ===== beginning time step 24 : 0.93445 ===== 1 Nonlinear solution status: time= 0.93445 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.562915e+06 1.076845e+02 0.000000e+00 energy 5.682418e+03 2.564996e+01 5.682418e+02 displacement 1.061830e+03 1.061830e+03 1.061830e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.93445 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.562915e+06 1.026045e+00 0.000000e+00 energy 5.682418e+03 6.632742e-02 5.682418e+02 displacement 1.061830e+03 1.397881e+00 1.083990e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.93445 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.562915e+06 4.480201e-03 0.000000e+00 energy 5.682418e+03 9.423614e-04 5.682418e+02 displacement 1.061830e+03 7.769086e-03 1.084510e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.93445 Data Record #1 =========================================================================== Step = 24 Time = 0.934449734 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 24 Time = 0.934449734 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 24 Time = 0.934449734 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 24 Time = 0.934449734 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0475598 ===== beginning time step 25 : 0.98201 ===== 1 Nonlinear solution status: time= 0.98201 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.020239e+06 3.280234e+02 0.000000e+00 energy 1.102581e+04 5.551947e+01 1.102581e+03 displacement 2.113473e+03 2.113473e+03 2.113473e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.98201 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.020239e+06 7.147250e+01 0.000000e+00 energy 1.102581e+04 1.192561e+00 1.102581e+03 displacement 2.113473e+03 5.203931e+00 2.176525e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.98201 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.020239e+06 1.496161e+00 0.000000e+00 energy 1.102581e+04 4.077403e-02 1.102581e+03 displacement 2.113473e+03 3.818281e-02 2.177085e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.98201 Data Record #1 =========================================================================== Step = 25 Time = 0.98200952 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 25 Time = 0.98200952 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 25 Time = 0.98200952 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 25 Time = 0.98200952 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0580478 MUST POINT CONTROLLER: adjusting time step. dt = 0.0179905 ===== beginning time step 26 : 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 4.681066e+05 9.849158e+00 0.000000e+00 energy 1.681778e+03 4.044829e+00 1.681778e+02 displacement 3.220569e+02 3.220569e+02 3.220569e-02 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 4.681066e+05 4.084355e-02 0.000000e+00 energy 1.681778e+03 2.163562e-03 1.681778e+02 displacement 3.220569e+02 1.042937e-01 3.260664e-02 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 4.681066e+05 9.598761e-04 0.000000e+00 energy 1.681778e+03 5.721371e-04 1.681778e+02 displacement 3.220569e+02 2.145156e-03 3.264041e-02 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 1 Data Record #1 =========================================================================== Step = 26 Time = 1 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 26 Time = 1 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 26 Time = 1 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 26 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 .................... : 26 Total number of equilibrium iterations ............ : 92 Average number of equilibrium iterations .......... : 3.53846 Total number of right hand evaluations ............ : 129 Total number of stiffness reformations ............ : 99 Time in linear solver: 0:00:53 T I M I N G I N F O R M A T I O N Input time ...................... : 0:00:00 (0 sec) Initialization time ............. : 0:00:00 (0 sec) Solve time ...................... : 0:01:20 (80 sec) IO-time (plot, dmp, data) .... : 0:00:00 (0 sec) reforming stiffness .......... : 0:00:00 (0 sec) evaluating stiffness ......... : 0:00:11 (11 sec) evaluating residual .......... : 0:00:04 (4 sec) model update ................. : 0:00:03 (3 sec) time in linear solver ........ : 0:00:53 (53 sec) Total elapsed time .............. : 0:01:20 (80 sec) N O R M A L T E R M I N A T I O N