=========================================================================== ________ _________ _________ __ _________ | |\ | |\ | |\ | |\ / \\ | ____|| | ____|| | __ || |__|| | ___ || | |\___\| | |\___\| | |\_| || \_\| | // \ || | || | || | || | || __ | || | || | ||__ | ||__ | ||_| || | |\ | || | || | |\ | |\ | \\ | || | || | || | ___|| | ___|| | ___ || | || | || | || | |\__\| | |\__\| | |\__| || | || | || | || | || | || | || | || | || | || | || | || | ||___ | ||__| || | || | \\__/ || | || | |\ | || | || | || |___|| |________|| |__________|| |__|| \_________// 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.001 Energy convergence tolerance ................... : 0.01 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+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.0742473 ===== beginning time step 5 : 0.274247 ===== 1 Nonlinear solution status: time= 0.274247 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.494692e+06 1.191390e+03 0.000000e+00 energy 3.380114e+04 2.387623e+02 3.380114e+02 displacement 4.560094e+03 4.560094e+03 4.560094e-03 Reforming stiffness matrix: reformation #1 2 ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29056 at gauss point 1 * * jacobian = -0.634017 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29057 at gauss point 1 * * jacobian = -1.15373 * * * ************************************************************************* ------- failed to converge at time : 0.274247 Retrying time step. Retry attempt 1 of max 5 AUTO STEPPER: retry step, dt = 0.0371237 ===== beginning time step 5 : 0.237124 ===== 1 Nonlinear solution status: time= 0.237124 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.373673e+06 1.313664e+02 0.000000e+00 energy 8.450289e+03 4.942657e+01 8.450289e+01 displacement 1.140026e+03 1.140026e+03 1.140026e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.237124 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.373673e+06 2.261951e+00 0.000000e+00 energy 8.450289e+03 1.192875e-01 8.450289e+01 displacement 1.140026e+03 1.857488e+00 1.114722e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.237124 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.373673e+06 1.497132e-02 0.000000e+00 energy 8.450289e+03 2.298471e-03 8.450289e+01 displacement 1.140026e+03 1.099029e-02 1.116362e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.237124 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.373673e+06 1.308617e-05 0.000000e+00 energy 8.450289e+03 1.902615e-05 8.450289e+01 displacement 1.140026e+03 4.877656e-04 1.116559e-03 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 0.237124 Data Record #1 =========================================================================== Step = 5 Time = 0.237123658 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 5 Time = 0.237123658 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 5 Time = 0.237123658 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 5 Time = 0.237123658 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0496989 ===== beginning time step 6 : 0.286823 ===== 1 Nonlinear solution status: time= 0.286823 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.107790e+06 3.546060e+02 0.000000e+00 energy 1.460884e+04 1.010002e+02 1.460884e+02 displacement 1.960830e+03 1.960830e+03 1.960830e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.286823 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.107790e+06 2.080487e+02 0.000000e+00 energy 1.460884e+04 2.607873e+00 1.460884e+02 displacement 1.960830e+03 5.993855e+00 1.896368e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.286823 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.107790e+06 3.862625e+00 0.000000e+00 energy 1.460884e+04 6.653300e-02 1.460884e+02 displacement 1.960830e+03 8.104542e-02 1.901772e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.286823 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.107790e+06 4.989412e-01 0.000000e+00 energy 1.460884e+04 2.002391e-03 1.460884e+02 displacement 1.960830e+03 5.946422e-02 1.902010e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.286823 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.107790e+06 1.056809e-01 0.000000e+00 energy 1.460884e+04 4.992603e-03 1.460884e+02 displacement 1.960830e+03 5.394683e-03 1.902529e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.286823 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.107790e+06 9.588650e-03 0.000000e+00 energy 1.460884e+04 2.288689e-03 1.460884e+02 displacement 1.960830e+03 4.140897e-03 1.903118e-03 Reforming stiffness matrix: reformation #6 7 Nonlinear solution status: time= 0.286823 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.107790e+06 2.197874e-05 0.000000e+00 energy 1.460884e+04 7.912883e-05 1.460884e+02 displacement 1.960830e+03 1.219732e-03 1.903452e-03 convergence summary number of iterations : 7 number of reformations : 7 ------- converged at time : 0.286823 Data Record #1 =========================================================================== Step = 6 Time = 0.286822584 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 6 Time = 0.286822584 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 6 Time = 0.286822584 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 6 Time = 0.286822584 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0595191 MUST POINT CONTROLLER: adjusting time step. dt = 0.0131774 ===== 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.761573e+05 2.726958e+00 0.000000e+00 energy 9.808341e+02 2.498588e+00 9.808341e+00 displacement 1.301083e+02 1.301083e+02 1.301083e-04 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.761573e+05 3.267329e-04 0.000000e+00 energy 9.808341e+02 1.450528e-04 9.808341e+00 displacement 1.301083e+02 2.742044e-02 1.290899e-04 ************************************************************************* * WARNING * * * * No force acting on the system. * ************************************************************************* convergence summary number of iterations : 2 number of reformations : 2 ------- 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.0676153 ===== beginning time step 8 : 0.367615 ===== 1 Nonlinear solution status: time= 0.367615 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.189821e+06 9.057794e+02 0.000000e+00 energy 2.552551e+04 1.837615e+02 2.552551e+02 displacement 3.372749e+03 3.372749e+03 3.372749e-03 Reforming stiffness matrix: reformation #1 2 ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29056 at gauss point 7 * * jacobian = -0.22557 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29057 at gauss point 8 * * jacobian = -0.17132 * * * ************************************************************************* ------- failed to converge at time : 0.367615 Retrying time step. Retry attempt 1 of max 5 AUTO STEPPER: retry step, dt = 0.0338077 ===== beginning time step 8 : 0.333808 ===== 1 Nonlinear solution status: time= 0.333808 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.797473e+06 9.401500e+01 0.000000e+00 energy 6.381440e+03 3.549902e+01 6.381440e+01 displacement 8.431975e+02 8.431975e+02 8.431975e-04 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.333808 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.797473e+06 9.005280e-01 0.000000e+00 energy 6.381440e+03 5.945738e-02 6.381440e+01 displacement 8.431975e+02 1.227189e+00 8.252123e-04 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.333808 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.797473e+06 1.767238e-03 0.000000e+00 energy 6.381440e+03 5.059465e-04 6.381440e+01 displacement 8.431975e+02 5.991108e-03 8.260598e-04 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.333808 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.797473e+06 1.528778e-07 0.000000e+00 energy 6.381440e+03 2.071368e-07 6.381440e+01 displacement 8.431975e+02 3.350183e-05 8.261057e-04 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 0.333808 Data Record #1 =========================================================================== Step = 8 Time = 0.333807654 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 8 Time = 0.333807654 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 8 Time = 0.333807654 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 8 Time = 0.333807654 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0470461 ===== beginning time step 9 : 0.380854 ===== 1 Nonlinear solution status: time= 0.380854 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.382202e+06 2.953459e+02 0.000000e+00 energy 1.200149e+04 8.067835e+01 1.200149e+02 displacement 1.568544e+03 1.568544e+03 1.568544e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.380854 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.382202e+06 7.193150e+01 0.000000e+00 energy 1.200149e+04 1.213987e+00 1.200149e+02 displacement 1.568544e+03 5.085092e+00 1.519246e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.380854 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.382202e+06 1.490354e+00 0.000000e+00 energy 1.200149e+04 4.022975e-02 1.200149e+02 displacement 1.568544e+03 2.726608e-02 1.522263e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.380854 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.382202e+06 1.390172e-01 0.000000e+00 energy 1.200149e+04 6.176656e-03 1.200149e+02 displacement 1.568544e+03 2.258315e-02 1.522758e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.380854 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.382202e+06 2.746419e-02 0.000000e+00 energy 1.200149e+04 3.194312e-03 1.200149e+02 displacement 1.568544e+03 3.897528e-03 1.523297e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.380854 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.382202e+06 2.750169e-04 0.000000e+00 energy 1.200149e+04 3.645048e-04 1.200149e+02 displacement 1.568544e+03 2.496441e-03 1.523797e-03 ************************************************************************* * WARNING * * * * No force acting on the system. * ************************************************************************* convergence summary number of iterations : 6 number of reformations : 6 ------- converged at time : 0.380854 Data Record #1 =========================================================================== Step = 9 Time = 0.380853777 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 9 Time = 0.380853777 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 9 Time = 0.380853777 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 9 Time = 0.380853777 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0576369 MUST POINT CONTROLLER: adjusting time step. dt = 0.0191462 ===== 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.401376e+05 1.148121e+01 0.000000e+00 energy 1.914553e+03 6.790518e+00 1.914553e+01 displacement 2.457663e+02 2.457663e+02 2.457663e-04 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.401376e+05 6.597118e-03 0.000000e+00 energy 1.914553e+03 1.444548e-03 1.914553e+01 displacement 2.457663e+02 1.484440e-01 2.429849e-04 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.401376e+05 1.315528e-07 0.000000e+00 energy 1.914553e+03 6.798403e-07 1.914553e+01 displacement 2.457663e+02 8.998235e-05 2.430357e-04 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.0661095 ===== beginning time step 11 : 0.46611 ===== 1 Nonlinear solution status: time= 0.46611 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.344917e+06 8.527566e+02 0.000000e+00 energy 2.249303e+04 1.561753e+02 2.249303e+02 displacement 2.866427e+03 2.866427e+03 2.866427e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.46611 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 2 step from line search = 0.112769 convergence norms : INITIAL CURRENT REQUIRED residual 6.344917e+06 6.485833e+02 0.000000e+00 energy 2.249303e+04 1.585913e+00 2.249303e+02 displacement 2.866427e+03 3.341638e-01 2.850231e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.46611 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.344917e+06 1.103773e+02 0.000000e+00 energy 2.249303e+04 1.825493e+00 2.249303e+02 displacement 2.866427e+03 1.864092e+01 2.749379e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.46611 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.344917e+06 2.250688e+00 0.000000e+00 energy 2.249303e+04 6.091619e-02 2.249303e+02 displacement 2.866427e+03 1.461623e-01 2.760134e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.46611 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.344917e+06 1.323980e-01 0.000000e+00 energy 2.249303e+04 6.968794e-03 2.249303e+02 displacement 2.866427e+03 7.949716e-03 2.761645e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.46611 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.344917e+06 2.106547e-02 0.000000e+00 energy 2.249303e+04 2.643548e-03 2.249303e+02 displacement 2.866427e+03 2.530886e-03 2.762267e-03 convergence summary number of iterations : 6 number of reformations : 6 ------- converged at time : 0.46611 Data Record #1 =========================================================================== Step = 11 Time = 0.466109519 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 11 Time = 0.466109519 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 11 Time = 0.466109519 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 11 Time = 0.466109519 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0728876 MUST POINT CONTROLLER: adjusting time step. dt = 0.0338905 ===== 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.602691e+06 9.621465e+01 0.000000e+00 energy 5.668189e+03 3.096698e+01 5.668189e+01 displacement 7.074811e+02 7.074811e+02 7.074811e-04 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.602691e+06 8.522627e-01 0.000000e+00 energy 5.668189e+03 5.791791e-02 5.668189e+01 displacement 7.074811e+02 2.045644e+00 6.988640e-04 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.602691e+06 2.225580e-03 0.000000e+00 energy 5.668189e+03 6.327595e-04 5.668189e+01 displacement 7.074811e+02 7.721639e-03 6.997709e-04 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 1.602691e+06 3.403701e-07 0.000000e+00 energy 5.668189e+03 9.034282e-07 5.668189e+01 displacement 7.074811e+02 7.139943e-05 6.998472e-04 convergence summary number of iterations : 4 number of reformations : 4 ------- 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.0783101 ===== beginning time step 13 : 0.57831 ===== 1 Nonlinear solution status: time= 0.57831 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.358779e+06 1.442659e+03 0.000000e+00 energy 2.962311e+04 1.727849e+02 2.962311e+02 displacement 3.709210e+03 3.709210e+03 3.709210e-03 Reforming stiffness matrix: reformation #1 2 ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29055 at gauss point 2 * * jacobian = -0.0105199 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29056 at gauss point 1 * * jacobian = -2.92233 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29057 at gauss point 1 * * jacobian = -4.41544 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29058 at gauss point 1 * * jacobian = -0.885257 * * * ************************************************************************* ------- failed to converge at time : 0.57831 Retrying time step. Retry attempt 1 of max 5 AUTO STEPPER: retry step, dt = 0.039155 ===== beginning time step 13 : 0.539155 ===== 1 Nonlinear solution status: time= 0.539155 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.089695e+06 1.599520e+02 0.000000e+00 energy 7.405779e+03 4.391260e+01 7.405779e+01 displacement 9.273031e+02 9.273031e+02 9.273031e-04 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.539155 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.089695e+06 3.793616e+00 0.000000e+00 energy 7.405779e+03 1.691184e-01 7.405779e+01 displacement 9.273031e+02 4.009801e+00 9.220024e-04 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.539155 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.089695e+06 6.835546e-02 0.000000e+00 energy 7.405779e+03 4.782409e-03 7.405779e+01 displacement 9.273031e+02 2.124153e-02 9.236068e-04 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.539155 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.089695e+06 1.300937e-03 0.000000e+00 energy 7.405779e+03 6.150126e-04 7.405779e+01 displacement 9.273031e+02 2.500195e-03 9.240979e-04 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.539155 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.089695e+06 3.083472e-07 0.000000e+00 energy 7.405779e+03 3.057500e-06 7.405779e+01 displacement 9.273031e+02 1.491189e-04 9.242350e-04 convergence summary number of iterations : 5 number of reformations : 5 ------- converged at time : 0.539155 Data Record #1 =========================================================================== Step = 13 Time = 0.539155046 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 13 Time = 0.539155046 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 13 Time = 0.539155046 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 13 Time = 0.539155046 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.051324 ===== beginning time step 14 : 0.590479 ===== 1 Nonlinear solution status: time= 0.590479 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.520596e+06 4.040446e+02 0.000000e+00 energy 1.248326e+04 8.259639e+01 1.248326e+02 displacement 1.594150e+03 1.594150e+03 1.594150e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.590479 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.520596e+06 2.525167e+02 0.000000e+00 energy 1.248326e+04 2.976297e+00 1.248326e+02 displacement 1.594150e+03 1.215919e+01 1.605323e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.590479 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.520596e+06 2.690311e+00 0.000000e+00 energy 1.248326e+04 8.289345e-03 1.248326e+02 displacement 1.594150e+03 2.964187e-01 1.614276e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.590479 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.520596e+06 1.942558e+00 0.000000e+00 energy 1.248326e+04 1.869916e-02 1.248326e+02 displacement 1.594150e+03 1.097701e-01 1.614227e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.590479 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.520596e+06 7.765584e-02 0.000000e+00 energy 1.248326e+04 6.352284e-03 1.248326e+02 displacement 1.594150e+03 7.114320e-03 1.614265e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.590479 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.520596e+06 9.778293e-03 0.000000e+00 energy 1.248326e+04 2.309575e-03 1.248326e+02 displacement 1.594150e+03 4.074343e-03 1.614894e-03 Reforming stiffness matrix: reformation #6 7 Nonlinear solution status: time= 0.590479 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.520596e+06 1.325189e-04 0.000000e+00 energy 1.248326e+04 1.879084e-04 1.248326e+02 displacement 1.594150e+03 1.349942e-03 1.615347e-03 convergence summary number of iterations : 7 number of reformations : 7 ------- converged at time : 0.590479 Data Record #1 =========================================================================== Step = 14 Time = 0.590479083 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 14 Time = 0.590479083 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 14 Time = 0.590479083 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 14 Time = 0.590479083 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.060827 MUST POINT CONTROLLER: adjusting time step. dt = 0.00952092 ===== 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.186630e+05 7.958784e-01 0.000000e+00 energy 4.213670e+02 7.732073e-01 4.213670e+00 displacement 5.699785e+01 5.699785e+01 5.699785e-05 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.186630e+05 2.551761e-05 0.000000e+00 energy 4.213670e+02 2.069355e-05 4.213670e+00 displacement 5.699785e+01 1.415118e-02 5.729060e-05 ************************************************************************* * WARNING * * * * No force acting on the system. * ************************************************************************* convergence summary number of iterations : 2 number of reformations : 2 ------- 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.0686616 ===== beginning time step 16 : 0.668662 ===== 1 Nonlinear solution status: time= 0.668662 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.150129e+06 1.027791e+03 0.000000e+00 energy 2.184945e+04 1.368998e+02 2.184945e+02 displacement 2.995896e+03 2.995896e+03 2.995896e-03 Reforming stiffness matrix: reformation #1 2 ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29056 at gauss point 7 * * jacobian = -0.0455604 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29057 at gauss point 8 * * jacobian = -0.0071746 * * * ************************************************************************* ------- failed to converge at time : 0.668662 Retrying time step. Retry attempt 1 of max 5 AUTO STEPPER: retry step, dt = 0.0343308 ===== beginning time step 16 : 0.634331 ===== 1 Nonlinear solution status: time= 0.634331 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.537537e+06 1.044762e+02 0.000000e+00 energy 5.462379e+03 2.937040e+01 5.462379e+01 displacement 7.489872e+02 7.489872e+02 7.489872e-04 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.634331 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.537537e+06 9.846280e-01 0.000000e+00 energy 5.462379e+03 6.416721e-02 5.462379e+01 displacement 7.489872e+02 2.119087e+00 7.620217e-04 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.634331 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.537537e+06 3.822406e-03 0.000000e+00 energy 5.462379e+03 8.478275e-04 5.462379e+01 displacement 7.489872e+02 3.442677e-02 7.643865e-04 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.634331 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.537537e+06 1.207093e-06 0.000000e+00 energy 5.462379e+03 3.174425e-06 5.462379e+01 displacement 7.489872e+02 1.389749e-04 7.645023e-04 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 0.634331 Data Record #1 =========================================================================== Step = 16 Time = 0.634330791 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 16 Time = 0.634330791 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 16 Time = 0.634330791 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 16 Time = 0.634330791 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0474646 ===== beginning time step 17 : 0.681795 ===== 1 Nonlinear solution status: time= 0.681795 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.909721e+06 3.323756e+02 0.000000e+00 energy 1.035415e+04 6.541783e+01 1.035415e+02 displacement 1.495524e+03 1.495524e+03 1.495524e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.681795 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.909721e+06 7.078590e+01 0.000000e+00 energy 1.035415e+04 1.182409e+00 1.035415e+02 displacement 1.495524e+03 6.699333e+00 1.532197e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.681795 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.909721e+06 1.485444e+00 0.000000e+00 energy 1.035415e+04 4.051090e-02 1.035415e+02 displacement 1.495524e+03 1.380368e-01 1.539204e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.681795 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.909721e+06 1.231316e-01 0.000000e+00 energy 1.035415e+04 5.754400e-03 1.035415e+02 displacement 1.495524e+03 2.312153e-02 1.539382e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.681795 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.909721e+06 2.270209e-02 0.000000e+00 energy 1.035415e+04 2.460811e-03 1.035415e+02 displacement 1.495524e+03 3.118006e-03 1.539944e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.681795 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.909721e+06 3.716919e-04 0.000000e+00 energy 1.035415e+04 3.236795e-04 1.035415e+02 displacement 1.495524e+03 1.665617e-03 1.540567e-03 ************************************************************************* * WARNING * * * * No force acting on the system. * ************************************************************************* convergence summary number of iterations : 6 number of reformations : 6 ------- converged at time : 0.681795 Data Record #1 =========================================================================== Step = 17 Time = 0.681795424 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 17 Time = 0.681795424 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 17 Time = 0.681795424 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 17 Time = 0.681795424 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0579717 MUST POINT CONTROLLER: adjusting time step. dt = 0.0182046 ===== 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.242960e+05 1.038906e+01 0.000000e+00 energy 1.513643e+03 4.872956e+00 1.513643e+01 displacement 2.338886e+02 2.338886e+02 2.338886e-04 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.242960e+05 4.472606e-03 0.000000e+00 energy 1.513643e+03 1.083255e-03 1.513643e+01 displacement 2.338886e+02 1.245723e-01 2.362810e-04 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.242960e+05 3.796443e-08 0.000000e+00 energy 1.513643e+03 1.595701e-07 1.513643e+01 displacement 2.338886e+02 1.977584e-04 2.363779e-04 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.0663774 ===== beginning time step 19 : 0.766377 ===== 1 Nonlinear solution status: time= 0.766377 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.626101e+06 1.021815e+03 0.000000e+00 energy 2.010037e+04 1.228885e+02 2.010037e+02 displacement 3.175479e+03 3.175479e+03 3.175479e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.766377 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 2 step from line search = 0.211845 convergence norms : INITIAL CURRENT REQUIRED residual 5.626101e+06 5.692442e+02 0.000000e+00 energy 2.010037e+04 2.987760e+00 2.010037e+02 displacement 3.175479e+03 1.049062e+00 3.194560e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.766377 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.626101e+06 4.006579e+01 0.000000e+00 energy 2.010037e+04 8.700460e-01 2.010037e+02 displacement 3.175479e+03 1.226064e+01 3.273139e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.766377 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.626101e+06 5.944648e-01 0.000000e+00 energy 2.010037e+04 2.692641e-02 2.010037e+02 displacement 3.175479e+03 1.300365e-01 3.280832e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.766377 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.626101e+06 5.204947e-02 0.000000e+00 energy 2.010037e+04 4.393026e-03 2.010037e+02 displacement 3.175479e+03 6.414114e-03 3.282044e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.766377 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.626101e+06 2.490341e-03 0.000000e+00 energy 2.010037e+04 7.893427e-04 2.010037e+02 displacement 3.175479e+03 1.426284e-03 3.282595e-03 convergence summary number of iterations : 6 number of reformations : 6 ------- converged at time : 0.766377 Data Record #1 =========================================================================== Step = 19 Time = 0.766377365 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 19 Time = 0.766377365 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 19 Time = 0.766377365 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 19 Time = 0.766377365 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0731019 MUST POINT CONTROLLER: adjusting time step. dt = 0.0336226 ===== 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.443341e+06 1.027583e+02 0.000000e+00 energy 5.181649e+03 2.562405e+01 5.181649e+01 displacement 8.728379e+02 8.728379e+02 8.728379e-04 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.443341e+06 8.129848e-01 0.000000e+00 energy 5.181649e+03 5.609203e-02 5.181649e+01 displacement 8.728379e+02 1.300414e+00 8.869619e-04 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.443341e+06 2.588320e-03 0.000000e+00 energy 5.181649e+03 6.476187e-04 5.181649e+01 displacement 8.728379e+02 9.081880e-03 8.871510e-04 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 1.443341e+06 6.433094e-07 0.000000e+00 energy 5.181649e+03 1.930210e-06 5.181649e+01 displacement 8.728379e+02 9.644786e-05 8.872505e-04 convergence summary number of iterations : 4 number of reformations : 4 ------- 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.0784815 ===== beginning time step 21 : 0.878482 ===== 1 Nonlinear solution status: time= 0.878482 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.885996e+06 1.668789e+03 0.000000e+00 energy 2.840573e+04 1.278914e+02 2.840573e+02 displacement 4.918516e+03 4.918516e+03 4.918516e-03 Reforming stiffness matrix: reformation #1 2 ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29055 at gauss point 2 * * jacobian = -0.081039 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29056 at gauss point 1 * * jacobian = -3.2822 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29057 at gauss point 1 * * jacobian = -4.90987 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29058 at gauss point 1 * * jacobian = -0.997477 * * * ************************************************************************* ------- failed to converge at time : 0.878482 Retrying time step. Retry attempt 1 of max 5 AUTO STEPPER: retry step, dt = 0.0392408 ===== beginning time step 21 : 0.839241 ===== 1 Nonlinear solution status: time= 0.839241 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.971499e+06 1.759180e+02 0.000000e+00 energy 7.101434e+03 3.751271e+01 7.101434e+01 displacement 1.229631e+03 1.229631e+03 1.229631e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.839241 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.971499e+06 3.957465e+00 0.000000e+00 energy 7.101434e+03 1.742467e-01 7.101434e+01 displacement 1.229631e+03 2.414881e+00 1.253941e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.839241 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.971499e+06 8.235944e-02 0.000000e+00 energy 7.101434e+03 4.889742e-03 7.101434e+01 displacement 1.229631e+03 1.722931e-02 1.253816e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.839241 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.971499e+06 2.976726e-03 0.000000e+00 energy 7.101434e+03 8.979094e-04 7.101434e+01 displacement 1.229631e+03 2.387971e-03 1.254403e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.839241 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.971499e+06 3.363375e-06 0.000000e+00 energy 7.101434e+03 1.385198e-05 7.101434e+01 displacement 1.229631e+03 3.093192e-04 1.254671e-03 convergence summary number of iterations : 5 number of reformations : 5 ------- converged at time : 0.839241 Data Record #1 =========================================================================== Step = 21 Time = 0.839240757 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 21 Time = 0.839240757 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 21 Time = 0.839240757 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 21 Time = 0.839240757 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0513926 ===== beginning time step 22 : 0.890633 ===== 1 Nonlinear solution status: time= 0.890633 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.407674e+06 4.354101e+02 0.000000e+00 energy 1.232080e+04 6.932575e+01 1.232080e+02 displacement 2.199533e+03 2.199533e+03 2.199533e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.890633 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.407674e+06 2.359760e+02 0.000000e+00 energy 1.232080e+04 2.842582e+00 1.232080e+02 displacement 2.199533e+03 7.268049e+00 2.260581e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.890633 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.407674e+06 2.700510e+00 0.000000e+00 energy 1.232080e+04 1.806820e-02 1.232080e+02 displacement 2.199533e+03 1.471134e-01 2.261777e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.890633 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.407674e+06 1.584796e+00 0.000000e+00 energy 1.232080e+04 1.618474e-02 1.232080e+02 displacement 2.199533e+03 1.034347e-01 2.262064e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.890633 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.407674e+06 7.157531e-02 0.000000e+00 energy 1.232080e+04 5.744769e-03 1.232080e+02 displacement 2.199533e+03 6.590670e-03 2.261996e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.890633 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.407674e+06 9.096181e-03 0.000000e+00 energy 1.232080e+04 1.938267e-03 1.232080e+02 displacement 2.199533e+03 3.604017e-03 2.262786e-03 Reforming stiffness matrix: reformation #6 7 Nonlinear solution status: time= 0.890633 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.407674e+06 1.029560e-04 0.000000e+00 energy 1.232080e+04 1.428653e-04 1.232080e+02 displacement 2.199533e+03 1.071908e-03 2.263328e-03 convergence summary number of iterations : 7 number of reformations : 7 ------- converged at time : 0.890633 Data Record #1 =========================================================================== Step = 22 Time = 0.890633362 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 22 Time = 0.890633362 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 22 Time = 0.890633362 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 22 Time = 0.890633362 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0608822 MUST POINT CONTROLLER: adjusting time step. dt = 0.00936664 ===== 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.149855e+05 7.596638e-01 0.000000e+00 energy 4.179128e+02 6.481787e-01 4.179128e+00 displacement 7.758832e+01 7.758832e+01 7.758832e-05 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.149855e+05 2.238158e-05 0.000000e+00 energy 4.179128e+02 1.929196e-05 4.179128e+00 displacement 7.758832e+01 7.564855e-03 7.803156e-05 ************************************************************************* * WARNING * * * * No force acting on the system. * ************************************************************************* convergence summary number of iterations : 2 number of reformations : 2 ------- 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.0687057 ===== beginning time step 24 : 0.968706 ===== 1 Nonlinear solution status: time= 0.968706 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.207529e+06 1.048934e+03 0.000000e+00 energy 2.258489e+04 1.052575e+02 2.258489e+02 displacement 4.222274e+03 4.222274e+03 4.222274e-03 Reforming stiffness matrix: reformation #1 2 ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29056 at gauss point 7 * * jacobian = -0.0400968 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 29057 at gauss point 8 * * jacobian = -0.00285081 * * * ************************************************************************* ------- failed to converge at time : 0.968706 Retrying time step. Retry attempt 1 of max 5 AUTO STEPPER: retry step, dt = 0.0343529 ===== beginning time step 24 : 0.934353 ===== 1 Nonlinear solution status: time= 0.934353 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.551886e+06 1.062509e+02 0.000000e+00 energy 5.646238e+03 2.540527e+01 5.646238e+01 displacement 1.055571e+03 1.055571e+03 1.055571e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.934353 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.551886e+06 9.960171e-01 0.000000e+00 energy 5.646238e+03 6.482866e-02 5.646238e+01 displacement 1.055571e+03 1.381488e+00 1.077537e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.934353 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.551886e+06 4.136299e-03 0.000000e+00 energy 5.646238e+03 8.885523e-04 5.646238e+01 displacement 1.055571e+03 7.544984e-03 1.078049e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.934353 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 1.551886e+06 1.879213e-06 0.000000e+00 energy 5.646238e+03 4.787667e-06 5.646238e+01 displacement 1.055571e+03 1.635683e-04 1.078206e-03 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 0.934353 Data Record #1 =========================================================================== Step = 24 Time = 0.934352864 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 24 Time = 0.934352864 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 24 Time = 0.934352864 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 24 Time = 0.934352864 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0474823 ===== beginning time step 25 : 0.981835 ===== 1 Nonlinear solution status: time= 0.981835 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.007993e+06 3.261170e+02 0.000000e+00 energy 1.098503e+04 5.527513e+01 1.098503e+02 displacement 2.106073e+03 2.106073e+03 2.106073e-03 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.981835 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.007993e+06 6.984781e+01 0.000000e+00 energy 1.098503e+04 1.172748e+00 1.098503e+02 displacement 2.106073e+03 5.168766e+00 2.168813e-03 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.981835 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.007993e+06 1.463207e+00 0.000000e+00 energy 1.098503e+04 4.025582e-02 1.098503e+02 displacement 2.106073e+03 3.765651e-02 2.169366e-03 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.981835 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.007993e+06 1.198215e-01 0.000000e+00 energy 1.098503e+04 5.657123e-03 1.098503e+02 displacement 2.106073e+03 2.284987e-02 2.169625e-03 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.981835 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.007993e+06 2.166571e-02 0.000000e+00 energy 1.098503e+04 2.349375e-03 1.098503e+02 displacement 2.106073e+03 3.007409e-03 2.170237e-03 Reforming stiffness matrix: reformation #5 6 Nonlinear solution status: time= 0.981835 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.007993e+06 3.653824e-04 0.000000e+00 energy 1.098503e+04 3.069569e-04 1.098503e+02 displacement 2.106073e+03 1.539898e-03 2.170946e-03 convergence summary number of iterations : 6 number of reformations : 6 ------- converged at time : 0.981835 Data Record #1 =========================================================================== Step = 25 Time = 0.981835155 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 25 Time = 0.981835155 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 25 Time = 0.981835155 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 25 Time = 0.981835155 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0579858 MUST POINT CONTROLLER: adjusting time step. dt = 0.0181648 ===== 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.512099e+05 9.820966e+00 0.000000e+00 energy 1.656173e+03 4.321620e+00 1.656173e+01 displacement 3.281216e+02 3.281216e+02 3.281216e-04 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.512099e+05 4.468018e-03 0.000000e+00 energy 1.656173e+03 1.090304e-03 1.656173e+01 displacement 3.281216e+02 1.060559e-01 3.320891e-04 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.512099e+05 3.664532e-08 0.000000e+00 energy 1.656173e+03 1.719953e-07 1.656173e+01 displacement 3.281216e+02 9.458134e-05 3.321200e-04 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 ............ : 132 Average number of equilibrium iterations .......... : 5.07692 Total number of right hand evaluations ............ : 168 Total number of stiffness reformations ............ : 139 Time in linear solver: 0:01:12 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:48 (108 sec) IO-time (plot, dmp, data) .... : 0:00:00 (0 sec) reforming stiffness .......... : 0:00:00 (0 sec) evaluating stiffness ......... : 0:00:17 (17 sec) evaluating residual .......... : 0:00:04 (4 sec) model update ................. : 0:00:03 (3 sec) time in linear solver ........ : 0:01:12 (72 sec) Total elapsed time .............. : 0:01:48 (108 sec) N O R M A L T E R M I N A T I O N