=========================================================================== ________ _________ _________ __ _________ | |\ | |\ | |\ | |\ / \\ | ____|| | ____|| | __ || |__|| | ___ || | |\___\| | |\___\| | |\_| || \_\| | // \ || | || | || | || | || __ | || | || | ||__ | ||__ | ||_| || | |\ | || | || | |\ | |\ | \\ | || | || | || | ___|| | ___|| | ___ || | || | || | || | |\__\| | |\__\| | |\__| || | || | || | || | || | || | || | || | || | || | || | || | ||___ | ||__| || | || | \\__/ || | || | |\ | || | || | || |___|| |________|| |__________|| |__|| \_________// F I N I T E E L E M E N T S F O R B I O M E C H A N I C S --- v e r s i o n - 2 . 4 . 3 . 8286 --- Musculoskeletal Research Laboratory University of Utah http://febio.org FEBio is a registered trademark. copyright (c) 2006-2015 - All rights reserved This is the NON-COMMERCIAL version of FEBio or the commercial license key file could not be found. If you have a key file make sure it is placed in the same directory as the executable. This version may only be used for non-commercial purposes as described in the license agreement. The functionality of this version may be limited compared to the commercial version. If you wish to obtain a valid commercial license file, please contact the developers. =========================================================================== tf_joint FILES USED =========================================================================== Input file : model_MOD.feb Plot file : model_MOD.xplt Log file : model_MOD.log MESH INFO =========================================================================== Number of materials ............................ : 4 Number of domains .............................. : 4 Number of nodes ................................ : 29877 Number of solid elements ....................... : 4096 Number of shell elements ....................... : 25220 Number of truss elements ....................... : 0 CONTROL DATA =========================================================================== Module type .................................... : solid Analysis type .................................. : quasi-static Plane strain mode .............................. : no Number of timesteps ............................ : 10 Time step size ................................. : 0.1 Auto time stepper activated .................... : yes Optimal nr of iterations ..................... : 10 Minimum allowable step size .................. : 0.001 Maximum allowable step size .................. : 0 Number of loadcurves ........................... : 7 Number of displacement boundary conditions ..... : 0 Number of concentrated nodal forces ............ : 0 Max nr of stiffness reformations ............... : 10 per time steps Max nr of Quasi-Newton iterations .............. : 0 between stiffness matrix reformations Displacement convergence tolerance ............. : 0.01 Energy convergence tolerance ................... : 0.1 Residual convergence tolerance ................. : 0 Minimal residual value ......................... : 0.001 Linesearch convergence tolerance ............... : 0.9 Minimum line search size ....................... : 0.01 Maximum number of line search iterations ....... : 5 Max condition number ........................... : 100000 OUTPUT DATA =========================================================================== plot level ................................ : must points only plotfile format ........................... : FEBIO MATERIAL DATA =========================================================================== 1 - femur (type: rigid body) density ........................................... : 1.132e-06 density ........................................... : 1.132e-06 E ................................................. : 1 v ................................................. : 0 parent_id ......................................... : -1 center_of_mass .................................... : 0,0,0 --------------------------------------------------------------------------- 2 - tibia (type: rigid body) density ........................................... : 1.132e-06 density ........................................... : 1.132e-06 E ................................................. : 1 v ................................................. : 0 parent_id ......................................... : -1 center_of_mass .................................... : 0,0,0 --------------------------------------------------------------------------- 3 - acl (type: trans iso Mooney-Rivlin) density ........................................... : 1.5e-09 k ................................................. : 73.2 laugon ............................................ : 0 atol .............................................. : 0.01 minaug ............................................ : 0 maxaug ............................................ : 0 c1 ................................................ : 1.95 c2 ................................................ : 0 c3 ................................................ : 0.0139 c4 ................................................ : 116.22 c5 ................................................ : 535.039 lam_max ........................................... : 1.046 --------------------------------------------------------------------------- 4 - aclfiber (type: trans iso Mooney-Rivlin) density ........................................... : 1.5e-09 k ................................................. : 73.2 laugon ............................................ : 0 atol .............................................. : 0.01 minaug ............................................ : 0 maxaug ............................................ : 0 c1 ................................................ : 1.95 c2 ................................................ : 0 c3 ................................................ : 0.0139 c4 ................................................ : 116.22 c5 ................................................ : 535.039 lam_max ........................................... : 1.046 RIGID BODY DATA =========================================================================== Rigid Body 1: material id : 1 center of mass : 0, 0, 0 --------------------------------------------------------------------------- Rigid Body 2: material id : 2 center of mass : 0, 0, 0 LOADCURVE DATA =========================================================================== 1 0 0 0.1 0.1 0.2 0.1 0.3 0.1 0.4 0.1 0.5 0.1 0.6 0.1 0.7 0.1 0.8 0.1 0.9 0.1 1 0.1 --------------------------------------------------------------------------- 2 0 0 0.1 0.1 0.9 0.9 1 1 --------------------------------------------------------------------------- 3 0 0 0.1 0.1 0.9 0.9 1 1 --------------------------------------------------------------------------- 4 0 0 0.1 0.1 0.9 0.9 1 1 --------------------------------------------------------------------------- 5 0 0 0.1 0.1 0.9 0.9 1 1 --------------------------------------------------------------------------- 6 0 0 0.1 0.1 0.9 0.9 1 1 --------------------------------------------------------------------------- 7 0 0 0.1 0.1 0.9 0.9 1 1 LINEAR SOLVER DATA =========================================================================== Solver type ............................... : Pardiso Matrix format ............................. : symmetric ===== reforming stiffness matrix: Nr of equations ........................... : 13272 Nr of nonzeroes in stiffness matrix ....... : 501420 ===== beginning time step 1 : 0.1 ===== 1 Nonlinear solution status: time= 0.1 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.119771e+07 3.550517e+03 0.000000e+00 energy 7.591579e+04 3.489976e+02 7.591579e+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.052 ===== beginning time step 5 : 0.252 ===== 1 Nonlinear solution status: time= 0.252 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.685680e+06 4.139116e+02 0.000000e+00 energy 1.663064e+04 1.162200e+02 1.663064e+03 displacement 2.237681e+03 2.237681e+03 2.237681e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.252 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.685680e+06 3.721095e+02 0.000000e+00 energy 1.663064e+04 4.044601e+00 1.663064e+03 displacement 2.237681e+03 7.574972e+00 2.162567e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.252 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 3 step from line search = 0.571809 convergence norms : INITIAL CURRENT REQUIRED residual 4.685680e+06 4.454727e+01 0.000000e+00 energy 1.663064e+04 1.428094e-02 1.663064e+03 displacement 2.237681e+03 2.397076e-01 2.162503e-01 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.252 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.685680e+06 4.627110e+00 0.000000e+00 energy 1.663064e+04 5.651471e-02 1.663064e+03 displacement 2.237681e+03 2.474993e-01 2.170235e-01 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.252 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.685680e+06 1.010770e-01 0.000000e+00 energy 1.663064e+04 1.176424e-02 1.663064e+03 displacement 2.237681e+03 2.644757e-02 2.171073e-01 convergence summary number of iterations : 5 number of reformations : 5 ------- converged at time : 0.252 Data Record #1 =========================================================================== Step = 5 Time = 0.252 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 5 Time = 0.252 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 5 Time = 0.252 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 5 Time = 0.252 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0616 MUST POINT CONTROLLER: adjusting time step. dt = 0.048 ===== beginning time step 6 : 0.3 ===== 1 Nonlinear solution status: time= 0.3 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.863224e+06 3.213228e+02 0.000000e+00 energy 1.360402e+04 9.144055e+01 1.360402e+03 displacement 1.798485e+03 1.798485e+03 1.798485e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.3 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.863224e+06 3.923393e+01 0.000000e+00 energy 1.360402e+04 8.121881e-01 1.360402e+03 displacement 1.798485e+03 5.069223e+00 1.742976e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.3 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.863224e+06 8.093544e-01 0.000000e+00 energy 1.360402e+04 3.127688e-02 1.360402e+03 displacement 1.798485e+03 2.757858e-02 1.746498e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.3 Data Record #1 =========================================================================== Step = 6 Time = 0.3 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 6 Time = 0.3 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 6 Time = 0.3 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 6 Time = 0.3 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0584 ===== beginning time step 7 : 0.3584 ===== 1 Nonlinear solution status: time= 0.3584 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.618402e+06 5.873079e+02 0.000000e+00 energy 1.955061e+04 1.358908e+02 1.955061e+03 displacement 2.517176e+03 2.517176e+03 2.517176e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.3584 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.618402e+06 2.104103e+02 0.000000e+00 energy 1.955061e+04 2.933200e+00 1.955061e+03 displacement 2.517176e+03 1.168306e+01 2.418766e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.3584 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.618402e+06 8.379445e+00 0.000000e+00 energy 1.955061e+04 1.436960e-01 1.955061e+03 displacement 2.517176e+03 1.267977e-01 2.427215e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.3584 Data Record #1 =========================================================================== Step = 7 Time = 0.3584 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 7 Time = 0.3584 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 7 Time = 0.3584 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 7 Time = 0.3584 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.06672 MUST POINT CONTROLLER: adjusting time step. dt = 0.0416 ===== beginning time step 8 : 0.4 ===== 1 Nonlinear solution status: time= 0.4 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.917733e+06 1.944340e+02 0.000000e+00 energy 9.798361e+03 5.681585e+01 9.798361e+02 displacement 1.191767e+03 1.191767e+03 1.191767e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.4 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.917733e+06 3.414187e+00 0.000000e+00 energy 9.798361e+03 1.296781e-01 9.798361e+02 displacement 1.191767e+03 3.174064e+00 1.160664e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.4 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.917733e+06 6.209949e-02 0.000000e+00 energy 9.798361e+03 4.161671e-03 9.798361e+02 displacement 1.191767e+03 8.189515e-03 1.162245e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.4 Data Record #1 =========================================================================== Step = 8 Time = 0.4 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 8 Time = 0.4 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 8 Time = 0.4 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 8 Time = 0.4 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.05328 ===== beginning time step 9 : 0.45328 ===== 1 Nonlinear solution status: time= 0.45328 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.158538e+06 4.465396e+02 0.000000e+00 energy 1.468150e+04 1.012220e+02 1.468150e+03 displacement 1.862668e+03 1.862668e+03 1.862668e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.45328 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.158538e+06 3.009722e+02 0.000000e+00 energy 1.468150e+04 3.484344e+00 1.468150e+03 displacement 1.862668e+03 1.032000e+01 1.801503e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.45328 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 3 step from line search = 0.248535 convergence norms : INITIAL CURRENT REQUIRED residual 4.158538e+06 1.112129e+02 0.000000e+00 energy 1.468150e+04 9.457704e-02 1.468150e+03 displacement 1.862668e+03 1.458713e-01 1.804789e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.45328 Data Record #1 =========================================================================== Step = 9 Time = 0.45328 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 9 Time = 0.45328 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 9 Time = 0.45328 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 9 Time = 0.45328 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.062624 MUST POINT CONTROLLER: adjusting time step. dt = 0.04672 ===== beginning time step 10 : 0.5 ===== 1 Nonlinear solution status: time= 0.5 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.945678e+06 6.862842e+02 0.000000e+00 energy 1.227788e+04 8.905515e+01 1.227788e+03 displacement 1.362317e+03 1.362317e+03 1.362317e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.5 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 2 step from line search = 0.241177 convergence norms : INITIAL CURRENT REQUIRED residual 3.945678e+06 3.224999e+02 0.000000e+00 energy 1.227788e+04 7.737420e-01 1.227788e+03 displacement 1.362317e+03 4.153535e-01 1.354059e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.5 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 3 step from line search = 0.479261 convergence norms : INITIAL CURRENT REQUIRED residual 3.945678e+06 4.562808e+01 0.000000e+00 energy 1.227788e+04 5.347444e-01 1.227788e+03 displacement 1.362317e+03 9.468475e-01 1.344139e-01 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.5 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.945678e+06 1.064518e+00 0.000000e+00 energy 1.227788e+04 2.253045e-02 1.227788e+03 displacement 1.362317e+03 1.023446e+00 1.336343e-01 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.5 stiffness updates = 0 right hand side evaluations = 8 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 3.945678e+06 1.476315e-01 0.000000e+00 energy 1.227788e+04 8.909225e-03 1.227788e+03 displacement 1.362317e+03 4.270354e-03 1.337115e-01 convergence summary number of iterations : 5 number of reformations : 5 ------- converged at time : 0.5 Data Record #1 =========================================================================== Step = 10 Time = 0.5 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 10 Time = 0.5 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 10 Time = 0.5 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 10 Time = 0.5 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.057376 ===== beginning time step 11 : 0.557376 ===== 1 Nonlinear solution status: time= 0.557376 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.531626e+06 5.529427e+02 0.000000e+00 energy 1.599270e+04 1.062949e+02 1.599270e+03 displacement 1.991651e+03 1.991651e+03 1.991651e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.557376 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.531626e+06 3.155128e+02 0.000000e+00 energy 1.599270e+04 3.581643e+00 1.599270e+03 displacement 1.991651e+03 1.865987e+01 1.974122e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.557376 stiffness updates = 0 right hand side evaluations = 5 stiffness matrix reformations = 3 step from line search = 0.468454 convergence norms : INITIAL CURRENT REQUIRED residual 4.531626e+06 6.031846e+01 0.000000e+00 energy 1.599270e+04 4.688352e-02 1.599270e+03 displacement 1.991651e+03 2.110440e-01 1.976847e-01 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.557376 stiffness updates = 0 right hand side evaluations = 6 stiffness matrix reformations = 4 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.531626e+06 4.614521e+00 0.000000e+00 energy 1.599270e+04 6.310905e-02 1.599270e+03 displacement 1.991651e+03 2.576300e-01 1.984470e-01 Reforming stiffness matrix: reformation #4 5 Nonlinear solution status: time= 0.557376 stiffness updates = 0 right hand side evaluations = 7 stiffness matrix reformations = 5 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 4.531626e+06 1.844697e-01 0.000000e+00 energy 1.599270e+04 1.170261e-02 1.599270e+03 displacement 1.991651e+03 1.590449e-02 1.984688e-01 convergence summary number of iterations : 5 number of reformations : 5 ------- converged at time : 0.557376 Data Record #1 =========================================================================== Step = 11 Time = 0.557376 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 11 Time = 0.557376 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 11 Time = 0.557376 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 11 Time = 0.557376 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0659008 MUST POINT CONTROLLER: adjusting time step. dt = 0.042624 ===== beginning time step 12 : 0.6 ===== 1 Nonlinear solution status: time= 0.6 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.471500e+06 2.210232e+02 0.000000e+00 energy 8.678397e+03 5.217662e+01 8.678397e+02 displacement 1.109796e+03 1.109796e+03 1.109796e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.6 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.471500e+06 7.856303e+00 0.000000e+00 energy 8.678397e+03 2.543781e-01 8.678397e+02 displacement 1.109796e+03 5.655221e+00 1.123770e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.6 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.471500e+06 2.426155e-01 0.000000e+00 energy 8.678397e+03 8.579017e-03 8.678397e+02 displacement 1.109796e+03 8.055279e-02 1.127472e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.6 Data Record #1 =========================================================================== Step = 12 Time = 0.6 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 12 Time = 0.6 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 12 Time = 0.6 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 12 Time = 0.6 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0540992 ===== beginning time step 13 : 0.654099 ===== 1 Nonlinear solution status: time= 0.654099 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.910936e+06 4.918059e+02 0.000000e+00 energy 1.376457e+04 8.805502e+01 1.376457e+03 displacement 1.859993e+03 1.859993e+03 1.859993e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.654099 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.910936e+06 1.263082e+02 0.000000e+00 energy 1.376457e+04 1.877071e+00 1.376457e+03 displacement 1.859993e+03 1.272149e+01 1.907092e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.654099 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.910936e+06 3.018883e+00 0.000000e+00 energy 1.376457e+04 7.194075e-02 1.376457e+03 displacement 1.859993e+03 4.382417e-01 1.921540e-01 Reforming stiffness matrix: reformation #3 4 Nonlinear solution status: time= 0.654099 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.910936e+06 2.021025e-01 0.000000e+00 energy 1.376457e+04 5.835816e-03 1.376457e+03 displacement 1.859993e+03 3.372034e-02 1.921261e-01 convergence summary number of iterations : 4 number of reformations : 4 ------- converged at time : 0.654099 Data Record #1 =========================================================================== Step = 13 Time = 0.6540992 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 13 Time = 0.6540992 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 13 Time = 0.6540992 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 13 Time = 0.6540992 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0632794 MUST POINT CONTROLLER: adjusting time step. dt = 0.0459008 ===== beginning time step 14 : 0.7 ===== 1 Nonlinear solution status: time= 0.7 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.809276e+06 2.994160e+02 0.000000e+00 energy 9.864326e+03 5.880828e+01 9.864326e+02 displacement 1.435957e+03 1.435957e+03 1.435957e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.7 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.809276e+06 1.508401e+01 0.000000e+00 energy 9.864326e+03 4.000222e-01 9.864326e+02 displacement 1.435957e+03 5.434105e+00 1.471242e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.7 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.809276e+06 3.167808e-01 0.000000e+00 energy 9.864326e+03 1.395780e-02 9.864326e+02 displacement 1.435957e+03 7.211429e-02 1.475871e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.7 Data Record #1 =========================================================================== Step = 14 Time = 0.7 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 14 Time = 0.7 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 14 Time = 0.7 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 14 Time = 0.7 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0567206 ===== beginning time step 15 : 0.756721 ===== 1 Nonlinear solution status: time= 0.756721 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.219899e+06 6.199743e+02 0.000000e+00 energy 1.492482e+04 9.184946e+01 1.492482e+03 displacement 2.318664e+03 2.318664e+03 2.318664e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.756721 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.219899e+06 2.179612e+02 0.000000e+00 energy 1.492482e+04 2.829771e+00 1.492482e+03 displacement 2.318664e+03 1.148293e+01 2.383736e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.756721 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.219899e+06 7.071198e+00 0.000000e+00 energy 1.492482e+04 1.065997e-01 1.492482e+03 displacement 2.318664e+03 1.596786e-01 2.387853e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.756721 Data Record #1 =========================================================================== Step = 15 Time = 0.75672064 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 15 Time = 0.75672064 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 15 Time = 0.75672064 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 15 Time = 0.75672064 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0653765 MUST POINT CONTROLLER: adjusting time step. dt = 0.0432794 ===== beginning time step 16 : 0.8 ===== 1 Nonlinear solution status: time= 0.8 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.673056e+06 2.686214e+02 0.000000e+00 energy 9.124119e+03 4.885738e+01 9.124119e+02 displacement 1.432180e+03 1.432180e+03 1.432180e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.8 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.673056e+06 4.911652e+00 0.000000e+00 energy 9.124119e+03 1.859729e-01 9.124119e+02 displacement 1.432180e+03 3.554115e+00 1.461909e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.8 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.673056e+06 9.815682e-02 0.000000e+00 energy 9.124119e+03 4.091572e-03 9.124119e+02 displacement 1.432180e+03 2.167472e-02 1.462004e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.8 Data Record #1 =========================================================================== Step = 16 Time = 0.8 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 16 Time = 0.8 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 16 Time = 0.8 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 16 Time = 0.8 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0546235 ===== beginning time step 17 : 0.854623 ===== 1 Nonlinear solution status: time= 0.854623 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.881004e+06 5.370861e+02 0.000000e+00 energy 1.389965e+04 7.962745e+01 1.389965e+03 displacement 2.383445e+03 2.383445e+03 2.383445e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.854623 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.881004e+06 2.084579e+02 0.000000e+00 energy 1.389965e+04 2.697968e+00 1.389965e+03 displacement 2.383445e+03 9.291165e+00 2.450447e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.854623 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.881004e+06 4.517996e+00 0.000000e+00 energy 1.389965e+04 7.458937e-02 1.389965e+03 displacement 2.383445e+03 1.391651e-01 2.450710e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.854623 Data Record #1 =========================================================================== Step = 17 Time = 0.854623488 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 17 Time = 0.854623488 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 17 Time = 0.854623488 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 17 Time = 0.854623488 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0636988 MUST POINT CONTROLLER: adjusting time step. dt = 0.0453765 ===== beginning time step 18 : 0.9 ===== 1 Nonlinear solution status: time= 0.9 stiffness updates = 0 right hand side evaluations = 2 stiffness matrix reformations = 1 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.943579e+06 3.169255e+02 0.000000e+00 energy 1.021344e+04 5.260959e+01 1.021344e+03 displacement 1.745433e+03 1.745433e+03 1.745433e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.9 stiffness updates = 0 right hand side evaluations = 3 stiffness matrix reformations = 2 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.943579e+06 9.243050e+00 0.000000e+00 energy 1.021344e+04 3.049535e-01 1.021344e+03 displacement 1.745433e+03 4.266442e+00 1.790663e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.9 stiffness updates = 0 right hand side evaluations = 4 stiffness matrix reformations = 3 step from line search = 1.000000 convergence norms : INITIAL CURRENT REQUIRED residual 2.943579e+06 1.739225e-01 0.000000e+00 energy 1.021344e+04 6.933333e-03 1.021344e+03 displacement 1.745433e+03 2.550691e-02 1.790901e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.9 Data Record #1 =========================================================================== Step = 18 Time = 0.9 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 18 Time = 0.9 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 18 Time = 0.9 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 18 Time = 0.9 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0563012 ===== beginning time step 19 : 0.956301 ===== 1 Nonlinear solution status: time= 0.956301 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.254427e+06 5.673916e+02 0.000000e+00 energy 1.536988e+04 7.816599e+01 1.536988e+03 displacement 2.836150e+03 2.836150e+03 2.836150e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.956301 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.254427e+06 2.734386e+02 0.000000e+00 energy 1.536988e+04 3.297840e+00 1.536988e+03 displacement 2.836150e+03 1.036403e+01 2.935863e-01 Reforming stiffness matrix: reformation #2 3 ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 28968 at gauss point 2 * * jacobian = -0.0589503 * * * ************************************************************************* ************************************************************************* * ERROR * * * * Negative jacobian was detected at element 28969 at gauss point 1 * * jacobian = -0.20685 * * * ************************************************************************* ------- failed to converge at time : 0.956301 Retrying time step. Retry attempt 1 of max 5 AUTO STEPPER: retry step, dt = 0.0281506 ===== beginning time step 19 : 0.928151 ===== 1 Nonlinear solution status: time= 0.928151 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.063766e+06 5.123909e+01 0.000000e+00 energy 3.843172e+03 1.480035e+01 3.843172e+02 displacement 7.092087e+02 7.092087e+02 7.092087e-02 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.928151 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.063766e+06 1.422136e-01 0.000000e+00 energy 3.843172e+03 1.312521e-02 3.843172e+02 displacement 7.092087e+02 6.262655e-01 7.217006e-02 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.928151 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.063766e+06 1.744146e-04 0.000000e+00 energy 3.843172e+03 1.484916e-04 3.843172e+02 displacement 7.092087e+02 2.844750e-03 7.221038e-02 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.928151 Data Record #1 =========================================================================== Step = 19 Time = 0.928150605 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 19 Time = 0.928150605 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 19 Time = 0.928150605 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 19 Time = 0.928150605 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0425205 ===== beginning time step 20 : 0.970671 ===== 1 Nonlinear solution status: time= 0.970671 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.405778e+06 2.242038e+02 0.000000e+00 energy 8.779191e+03 4.289738e+01 8.779191e+02 displacement 1.675616e+03 1.675616e+03 1.675616e-01 Reforming stiffness matrix: reformation #1 2 Nonlinear solution status: time= 0.970671 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.405778e+06 1.311451e+01 0.000000e+00 energy 8.779191e+03 3.687095e-01 8.779191e+02 displacement 1.675616e+03 3.271715e+00 1.719819e-01 Reforming stiffness matrix: reformation #2 3 Nonlinear solution status: time= 0.970671 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.405778e+06 3.146254e-01 0.000000e+00 energy 8.779191e+03 1.359912e-02 8.779191e+02 displacement 1.675616e+03 1.482920e-02 1.720047e-01 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 0.970671 Data Record #1 =========================================================================== Step = 20 Time = 0.970671089 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 20 Time = 0.970671089 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 20 Time = 0.970671089 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 20 Time = 0.970671089 Data = tibia motion File = tibia_kinematics.txt AUTO STEPPER: increasing time step, dt = 0.0540164 MUST POINT CONTROLLER: adjusting time step. dt = 0.0293289 ===== beginning time step 21 : 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 1.198407e+06 5.952869e+01 0.000000e+00 energy 4.355306e+03 1.620623e+01 4.355306e+02 displacement 8.431587e+02 8.431587e+02 8.431587e-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 1.198407e+06 2.069531e-01 0.000000e+00 energy 4.355306e+03 1.658710e-02 4.355306e+02 displacement 8.431587e+02 7.327976e-01 8.596767e-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 1.198407e+06 4.265507e-04 0.000000e+00 energy 4.355306e+03 2.828357e-04 4.355306e+02 displacement 8.431587e+02 3.728636e-03 8.601644e-02 convergence summary number of iterations : 3 number of reformations : 3 ------- converged at time : 1 Data Record #1 =========================================================================== Step = 21 Time = 1 Data = femur load File = femur_kinetics.txt Data Record #2 =========================================================================== Step = 21 Time = 1 Data = tibia load File = tibia_kinetics.txt Data Record #3 =========================================================================== Step = 21 Time = 1 Data = femur motion File = femur_kinematics.txt Data Record #4 =========================================================================== Step = 21 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 .................... : 21 Total number of equilibrium iterations ............ : 76 Average number of equilibrium iterations .......... : 3.61905 Total number of right hand evaluations ............ : 104 Total number of stiffness reformations ............ : 78 Time in solver: 0:00:46 Elapsed time : 0:01:05 N O R M A L T E R M I N A T I O N